Search

We provide an unsupervised adaptive sampling strategy capable of producing μs-timescale molecular dynamics (MD) simulations of large biosystems using many-body polarizable force fields (PFFs). The global exploration problem is decomposed into a set of separate MD trajectories that can be restarted within a selective process to achieve sufficient phase-space sampling. Accurate statistical properties can be obtained through reweighting. Within this highly parallel setup, the Tinker-HP package can be powered by an arbitrary large number of GPUs on supercomputers, reducing exploration time from years to days. This approach is used to tackle the urgent modeling problem of the SARS-CoV-2 Main Protease (Mpro) producing more than 38 μs of all-atom simulations of its apo (ligand-free) dimer using the high-resolution AMOEBA PFF. The first 15.14 μs simulation (physiological pH) is compared to available non-PFF long-timescale simulation data. A detailed clustering analysis exhibits striking differences between FFs, with AMOEBA showing a richer conformational space. Focusing on key structural markers related to the oxyanion hole stability, we observe an asymmetry between protomers. One of them appears less structured resembling the experimentally inactive monomer for which a 6 μs simulation was performed as a basis for comparison. Results highlight the plasticity of the Mpro active site. The C-terminal end of its less structured protomer is shown to oscillate between several states, being able to interact with the other protomer, potentially modulating its activity. Active and distal site volumes are found to be larger in the most active protomer within our AMOEBA simulations compared to non-PFFs as additional cryptic pockets are uncovered. A second 17 μs AMOEBA simulation is performed with protonated His172 residues mimicking lower pH. Data show the protonation impact on the destructuring of the oxyanion loop. We finally analyze the solvation patterns around key histidine residues. The confined AMOEBA polarizable water molecules are able to explore a wide range of dipole moments, going beyond bulk values, leading to a water molecule count consistent with experimental data. Results suggest that the use of PFFs could be critical in drug discovery to accurately model the complexity of the molecular interactions structuring Mpro., We provide an unsupervised adaptive sampling strategy capable of producing μs-timescale molecular dynamics (MD) simulations of large biosystems using many-body polarizable force fields (PFFs).

The estrogen receptors α (ERα) are transcription factors involved in several physiological processes belonging to the nuclear receptors (NRs) protein family. Besides the endogenous ligands, several other chemicals are able to bind to those receptors. Among them are endocrine disrupting chemicals (EDCs) that can trigger toxicological pathways. Many studies have focused on predicting EDCs based on their ability to bind NRs, mainly, estrogen receptors (ER), thyroid hormones receptors (TR), androgen receptors (AR), glucocorticoid receptors (GR), and peroxisome proliferator-activated receptors gamma (PPARγ). In this work, we suggest a pipeline designed for the prediction of ERα binding activity. The flagged compounds can be further explored using experimental techniques to assess their potential to be EDCs. The pipeline is a combination of structure based (docking and pharmacophore models) and ligand based (pharmacophore models) methods. The models have been constructed using the Environmental Protection Agency (EPA) data encompassing a large number of structurally diverse compounds. A validation step was then achieved using two external databases: the NR-DBIND (Nuclear Receptors DataBase Including Negative Data) and the EADB (Estrogenic Activity DataBase). Different combination protocols were explored. Results showed that the combination of models performed better than each model taken individually. The consensus protocol that reached values of 0.81 and 0.54 for sensitivity and specificity, respectively, was the best suited for our toxicological study. Insights and recommendations were drawn to alleviate the screening quality of other projects focusing on ERα binding predictions.

The Drug Design Data Resource (D3R) Grand Challenges are blind contests organized to assess the state-of-the-art methods accuracy in predicting binding modes and relative binding free energies of experimentally validated ligands for a given target. The second stage of the D3R Grand Challenge 2 (GC2) was focused on ranking 102 compounds according to their predicted affinity for Farnesoid X Receptor. In this task, our workflow was ranked 5th out of the 77 submissions in the structure-based category. Our strategy consisted in (1) a combination of molecular docking using AutoDock 4.2 and manual edition of available structures for binding poses generation using SeeSAR, (2) the use of HYDE scoring for pose selection, and (3) a hierarchical ranking using HYDE and MM/GBSA. In this report, we detail our pose generation and ligands ranking protocols and provide guidelines to be used in a prospective computer aided drug design program.

International audience; The androgen receptor (AR) is a transcription factor that plays a key role in sexual phenotype and neuromuscular development. AR can be modulated by exogenous compounds such as pharmaceuticals or chemicals present in the environment, and particularly by AR agonist compounds that mimic the action of endogenous agonist ligands and whether restore or alter the AR endocrine system functions. The activation of AR must be correctly balanced and identifying potent AR agonist compounds is of high interest to both propose treatments for certain diseases, or to predict the risk related to agonist chemicals exposure. The development of in silico approaches and the publication of structural, affinity and activity data provide a good framework to develop rational AR hits prediction models. Herein, we present a docking and a pharmacophore modeling strategy to help identifying AR agonist compounds. All models were trained on the NR-DBIND that provides high quality binding data on AR and tested on AR-agonist activity assays from the Tox21 initiative. Both methods display high performance on the NR-DBIND set and could serve as starting point for biologists and toxicologists. Yet, the pharmacophore models still need data feeding to be used as large scope undesired effect prediction models.

Screening Explorer is a web-based application that allows for an intuitive evaluation of the results of screening experiments using complementary metrics in the field. The usual evaluation of screening results implies the separate generation and apprehension of the ROC, predictiveness, and enrichment curves and their global metrics. Similarly, partial metrics need to be calculated repeatedly for different fractions of a data set and there exists no handy tool that allows reading partial metrics simultaneously on different charts. For a deeper understanding of the results of screening experiments, we rendered their analysis straightforward by linking these metrics interactively in an interactive usable web-based application. We also implemented simple consensus scoring methods based on scores normalization, standardization (z-scores), and compounds ranking to evaluate the enrichments that can be expected through methods combination. Two demonstration data sets allow the users to easily apprehend the functions of this tool that can be applied to the analysis of virtual and experimental screening results. Screening Explorer is freely accessible at http://stats.drugdesign.fr .

International audience; In this study, a surface acoustic wave (SAW) biosensor has been investigated to quantify the affinity and to estimate the binding kinetic between tumor necrosis factor α (TNF-α) and its inhibitor: SPD-304. To the best of our knowledge this is first ever report on real time and label free monitoring of TNF-α/ SPD-304 affinity. Effects of the SPD-304 solvent and DMSO cosolvent on TNFα/SPD-304 interaction have been investigated. Gravimetric results indicate that the limit of detection (LOD) of the developed biosensor was of order of 10 nM and that DMSO cosolvent influences the kinetic of interaction, the saturation value, and thus the value of the dissociation constant (Kd) of the TNF-α / SPD-304 system.

Nuclear receptors (NRs) are transcription factors that regulate gene expression in various physiological processes through their interactions with small hydrophobic molecules. They constitute an important class of targets for drugs and endocrine disruptors and are widely studied for both health and environment concerns. Since the integration of negative data can be critical for accurate modeling of ligand activity profiles, we manually collected and annotated NRs interaction data (positive and negative) through a sharp review of the corresponding literature. 15 116 positive and negative interactions data are provided for 28 NRs together with 593 PDB structures in the freely available Nuclear Receptors Database Including Negative Data ( http://nr-dbind.drugdesign.fr ). The NR-DBIND contains the most extensive information about interaction data on NRs, which should bring valuable information to chemists, biologists, pharmacologists and toxicologists.

International audience; In this work, a serie of cyclocoumarol derivatives was designed, synthesized, characterized and studied for their potentialities as selective inhibitors of COX-2. All target compounds have been screened for their anti-inflammatory activity by the assay of PGE2 production. Among them, compound 5d exhibited the most potent inhibitory activity with a PGE2 inhibition compared to NS-398 (79% and 88% respectively) and showed non-inhibitory activity towards the COX-1 enzyme. Docking studies revealed the capacity of this compound to occupy the selective COX-2 cavity establishing additional hydrogen bonds between the oxygen of the methoxy group and the His90 and Arg513 of the binding site of the enzyme.

Nuclear receptors (NRs) constitute an important class of therapeutic targets. During the last 4 years, we tackled the pharmacological profile assessment of NR ligands for which we constructed the NRLiSt BDB. We evaluated and compared the performance of different virtual screening approaches: mean of molecular descriptor distribution values, molecular docking and 3D pharmacophore models. The simple comparison of the distribution profiles of 4885 molecular descriptors between the agonist and antagonist datasets didn't provide satisfying results. We obtained an overall good performance with the docking method we used, Surflex-Dock which was able to discriminate agonist from antagonist ligands. But the availability of PDB structures in the "pharmacological-profile-to-predict-bound-state" (agonist-bound or antagonist-bound) and the availability of enough ligands of both pharmacological profiles constituted limits to generalize this protocol for all NRs. Finally, the 3D pharmacophore modeling approach, allowed us to generate selective agonist pharmacophores and selective antagonist pharmacophores that covered more than 99 % of the whole NRLiSt BDB. This study allowed a better understanding of the pharmacological modulation of NRs with small molecules and could be extended to other therapeutic classes.

International audience; Inflammation is a complex phenomenon necessary in human defense mechanisms but also involved in the development of some human diseases. The discovery of cyclooxygenase-2 (COX-2) improved the pharmacology of nonsteroidal anti-inflammatory drugs (NSAID) giving a clear mechanism for prostaglandin regulation in vivo and providing a new target for the development of COX-2-selective drugs without gastrointestinal side-effects. Keeping in view the importance of this pharmacological class, several literature reports have underlined the impact of these antiinflammatory compounds in therapeutics. The present review considers the most recently published literature concerning COX-2 inhibitors until 2016. Through a wide chemical classification, the last developments concerning this therapeutic family by highlighting structure-activity relationships insights and mechanisms are presented. A summary of the principal adverse effects observed and an overview of the new potential therapeutic indications for COX-2 inhibitors are also reported.

In the present work, we describe and discuss two interactive visual­ ization techniques recently added to the Udock software. First, we propose to display a spherical projection of the protein's molecular surface properties. The resulting 2D map allows to get in one glance a global view of the protein. Second, we propose to let the user choose a specific level of detail when visualizing the protein's 3D surface. We describe a simple smoothing algorithm that achieves this goal in real time. These techniques are designed to enhance the usability of molecular visualization and interactive simulation software.

International audience; Proteasome inhibition is a promising strategy for treating cancers. Herein, we report the discovery of novel drug-like inhibitors of mammalian proteasome 20S using a multistep structure-based virtual ligand screening strategy. Sulfone- or piperazine-containing hits essentially belong to the under-represented class of noncovalent and nonpeptidic proteasome inhibitors. Several of our compounds act in the micromolar range and are cytotoxic on human tumoral cell lines. Optimization of these molecules could lead to better anticancer therapy.

In silico screening based on the structures of the ligands or of the receptors has become an essential tool to facilitate the drug discovery process but compound collections are needed to carry out such in silico experiments. It has been recognized that absorption, distribution, metabolism, excretion and toxicity (ADME/tox) are key properties that need to be considered early on, even during the database preparation stage. FAF-Drugs is an online service based on Frowns (a chemoinformatics toolkit) that allows users to process their own compound collections via simple ADME/Tox filtering rules such as molecular weight, polar surface area, logP or number of rotatable bonds. SMILES (Simplified Molecular Input Line Entry System), CANSMILES (canonical smiles) or SDF (structure data file) files are required as input and molecules that pass or do not pass the filters are sent back in CANSMILES format. This service should thus help scientists engaging in drug discovery campaigns. Other utilities and several compound collections suitable for in silico screening are available at our site. FAF-Drugs can be accessed at http://bioserv.rpbs.jussieu.fr/FAFDrugs.html.

CDC25 dual-specificity phosphatases are essential key regulators of eukaryotic cell cycle progression and the CDC25A and B isoforms are over-expressed in different tumors and related cancer cell lines. CDC25s are now considered to be interesting targets in the search for novel anticancer agents. We describe new compounds derived from vitamin K3 that inhibit CDC25B activity with IC50 values in the low micromolar range. These naphthoquinone derivatives also display antiproliferative activity on HeLa cells as expected for CDC25 inhibitors and inhibit cell growth in a clonogenic assay at submicromolar concentrations. They increase inhibitory tyrosine 15 phosphorylation of CDK and induce the cleavage of PARP, a hallmark of apoptosis.

Nuclear receptors (NRs) constitute an important class of therapeutic targets. We evaluated the performance of 3D structure-based and ligand-based pharmacophore models in predicting the pharmacological profile of NRs ligands using the NRLiSt BDB database. We could generate selective pharmacophores for agonist and antagonist ligands and we found that the best performances were obtained by combining the structure-based and the ligand-based approaches. The combination of pharmacophores that were generated allowed to cover most of the chemical space of the NRLiSt BDB datasets. By screening the whole NRLiSt BDB on our 3D pharmacophores, we demonstrated their selectivity towards their dedicated NRs ligands. The 3D pharmacophores herein presented can thus be used as a predictor of the pharmacological activity of NRs ligands.Graphical AbstractUsing a combination of structure-based and ligand-based pharmacophores, agonist and antagonist ligands of the Nuclear Receptors included in the NRLiSt BDB database could be separated Electronic supplementary material The online version of this article (doi:10.1186/s13321-016-0154-2) contains supplementary material, which is available to authorized users.