DG-AMMOS: A New tool to generate 3D conformation of small molecules using Distance Geometry and Automated Molecular Mechanics Optimization for in silico Screening

<p>Abstract</p> <p>Background</p> <p>Discovery of new bioactive molecules that could enter drug discovery programs or that could serve as chemical probes is a very complex and costly endeavor. Structure-based and ligand-based <it>in silico </it>screening approaches are nowadays extensively used to complement experimental screening approaches in order to increase the effectiveness of the process and facilitating the screening of thousands or millions of small molecules against a biomolecular target. Both <it>in silico </it>screening methods require as input a suitable chemical compound collection and most often the 3D structure of the small molecules has to be generated since compounds are usually delivered in 1D SMILES, CANSMILES or in 2D SDF formats.</p> <p>Results</p> <p>Here, we describe the new open source program DG-AMMOS which allows the generation of the 3D conformation of small molecules using Distance Geometry and their energy minimization via Automated Molecular Mechanics Optimization. The program is validated on the Astex dataset, the ChemBridge Diversity database and on a number of small molecules with known crystal structures extracted from the Cambridge Structural Database. A comparison with the free program Balloon and the well-known commercial program Omega generating the 3D of small molecules is carried out. The results show that the new free program DG-AMMOS is a very efficient 3D structure generator engine.</p> <p>Conclusion</p> <p>DG-AMMOS provides fast, automated and reliable access to the generation of 3D conformation of small molecules and facilitates the preparation of a compound collection prior to high-throughput virtual screening computations. The validation of DG-AMMOS on several different datasets proves that generated structures are generally of equal quality or sometimes better than structures obtained by other tested methods.</p

MS-DOCK: Accurate multiple conformation generator and rigid docking protocol for multi-step virtual ligand screening

<p>Abstract</p> <p>Background</p> <p>The number of protein targets with a known or predicted tri-dimensional structure and of drug-like chemical compounds is growing rapidly and so is the need for new therapeutic compounds or chemical probes. Performing flexible structure-based virtual screening computations on thousands of targets with millions of molecules is intractable to most laboratories nor indeed desirable. Since shape complementarity is of primary importance for most protein-ligand interactions, we have developed a tool/protocol based on rigid-body docking to select compounds that fit well into binding sites.</p> <p>Results</p> <p>Here we present an efficient multiple conformation rigid-body docking approach, MS-DOCK, which is based on the program DOCK. This approach can be used as the first step of a multi-stage docking/scoring protocol. First, we developed and validated the Multiconf-DOCK tool that generates several conformers per input ligand. Then, each generated conformer (bioactives and 37970 decoys) was docked rigidly using DOCK6 with our optimized protocol into seven different receptor-binding sites. MS-DOCK was able to significantly reduce the size of the initial input library for all seven targets, thereby facilitating subsequent more CPU demanding flexible docking procedures.</p> <p>Conclusion</p> <p>MS-DOCK can be easily used for the generation of multi-conformer libraries and for shape-based filtering within a multi-step structure-based screening protocol in order to shorten computation times.</p

FAF-Drugs2: Free ADME/tox filtering tool to assist drug discovery and chemical biology projects

<p>Abstract</p> <p>Background</p> <p>Drug discovery and chemical biology are exceedingly complex and demanding enterprises. In recent years there are been increasing awareness about the importance of predicting/optimizing the absorption, distribution, metabolism, excretion and toxicity (ADMET) properties of small chemical compounds along the search process rather than at the final stages. Fast methods for evaluating ADMET properties of small molecules often involve applying a set of simple empirical rules (educated guesses) and as such, compound collections' property profiling can be performed <it>in silico</it>. Clearly, these rules cannot assess the full complexity of the human body but can provide valuable information and assist decision-making.</p> <p>Results</p> <p>This paper presents FAF-Drugs2, a free adaptable tool for ADMET filtering of electronic compound collections. FAF-Drugs2 is a command line utility program (e.g., written in Python) based on the open source chemistry toolkit OpenBabel, which performs various physicochemical calculations, identifies key functional groups, some toxic and unstable molecules/functional groups. In addition to filtered collections, FAF-Drugs2 can provide, via Gnuplot, several distribution diagrams of major physicochemical properties of the screened compound libraries.</p> <p>Conclusion</p> <p>We have developed FAF-Drugs2 to facilitate compound collection preparation, prior to (or after) experimental screening or virtual screening computations. Users can select to apply various filtering thresholds and add rules as needed for a given project. As it stands, FAF-Drugs2 implements numerous filtering rules (23 physicochemical rules and 204 substructure searching rules) that can be easily tuned.</p

AMMOS: Automated Molecular Mechanics Optimization tool for in silico Screening

<p>Abstract</p> <p>Background</p> <p>Virtual or <it>in silico </it>ligand screening combined with other computational methods is one of the most promising methods to search for new lead compounds, thereby greatly assisting the drug discovery process. Despite considerable progresses made in virtual screening methodologies, available computer programs do not easily address problems such as: structural optimization of compounds in a screening library, receptor flexibility/induced-fit, and accurate prediction of protein-ligand interactions. It has been shown that structural optimization of chemical compounds and that post-docking optimization in multi-step structure-based virtual screening approaches help to further improve the overall efficiency of the methods. To address some of these points, we developed the program AMMOS for refining both, the 3D structures of the small molecules present in chemical libraries and the predicted receptor-ligand complexes through allowing partial to full atom flexibility through molecular mechanics optimization.</p> <p>Results</p> <p>The program AMMOS carries out an automatic procedure that allows for the structural refinement of compound collections and energy minimization of protein-ligand complexes using the open source program AMMP. The performance of our package was evaluated by comparing the structures of small chemical entities minimized by AMMOS with those minimized with the Tripos and MMFF94s force fields. Next, AMMOS was used for full flexible minimization of protein-ligands complexes obtained from a mutli-step virtual screening. Enrichment studies of the selected pre-docked complexes containing 60% of the initially added inhibitors were carried out with or without final AMMOS minimization on two protein targets having different binding pocket properties. AMMOS was able to improve the enrichment after the pre-docking stage with 40 to 60% of the initially added active compounds found in the top 3% to 5% of the entire compound collection.</p> <p>Conclusion</p> <p>The open source AMMOS program can be helpful in a broad range of <it>in silico </it>drug design studies such as optimization of small molecules or energy minimization of pre-docked protein-ligand complexes. Our enrichment study suggests that AMMOS, designed to minimize a large number of ligands pre-docked in a protein target, can successfully be applied in a final post-processing step and that it can take into account some receptor flexibility within the binding site area.</p

In Silico Modeling the Effects of Missense Mutations Causing Snyder-Robinson Syndrome and Rescuing the Effects by Small Molecules Binding

Snyder-Robinson Syndrome (SRS) is an X-linked mental retardation disorder. Three missense mutations (G56S, V132G and I150T) on human spermine synthase (SMS) were reported to cause SRS. SMS is an important enzyme which converts spermindine into spermine, both of which are two polyamines controlling the normal cell growth and development. In vitro experiments showed that the dimer conformation played a crucial role on the SMS function. Our in silico studies including energy calculation, pKa calculation and molecular dynamics (MD) simulation based on the available 3D structure of SMS revealed that these mutations affected SMS function by affecting the dimer affinity, monomer stability or hydrogen bond network. One of the above sites, G56S, is accessible from the water phase, thus it provides the opportunity to rescue the disease-causing effect by binding an appropriate small molecule to the vicinity of the mutation site. Currently we run MD simulation to generate multiple receptor conformations and identified two potent binding pockets. Then two programs, Surflex and Autodock Vina, were applied for structure-based virtual screening (SBVS) and a consensus list of about 200 common compounds selected by both of the programs was created, and these compounds were tested experimentally by our collaborators

Serca1 Truncated Proteins Unable to Pump Calcium Reduce the Endoplasmic Reticulum Calcium Concentration and Induce Apoptosis

By pumping calcium from the cytosol to the ER, sarco/endoplasmic reticulum calcium ATPases (SERCAs) play a major role in the control of calcium signaling. We describe two SERCA1 splice variants (S1Ts) characterized by exon 4 and/or exon 11 splicing, encoding COOH terminally truncated proteins, having only one of the seven calcium-binding residues, and thus unable to pump calcium. As shown by semiquantitative RT-PCR, S1T transcripts are differentially expressed in several adult and fetal human tissues, but not in skeletal muscle and heart. S1T proteins expression was detected by Western blot in nontransfected cell lines. In transiently transfected cells, S1T homodimers were revealed by Western blot using mildly denaturing conditions. S1T proteins were shown, by confocal scanning microscopy, to colocalize with endogenous SERCA2b into the ER membrane. Using ER-targeted aequorin (erAEQ), we have found that S1T proteins reduce ER calcium and reverse elevation of ER calcium loading induced by SERCA1 and SERCA2b. Our results also show that SERCA1 variants increase ER calcium leakage and are consistent with the hypothesis of a cation channel formed by S1T homodimers. Finally, when overexpressed in liver-derived cells, S1T proteins significantly induce apoptosis. These data reveal a further mechanism modulating Ca2+ accumulation into the ER of nonmuscle cells and highlight the relevance of S1T proteins to the control of apoptosis

Tyrosine Kinase Syk Non-Enzymatic Inhibitors and Potential Anti-Allergic Drug-Like Compounds Discovered by Virtual and In Vitro Screening

In the past decade, the spleen tyrosine kinase (Syk) has shown a high potential for the discovery of new treatments for inflammatory and autoimmune disorders. Pharmacological inhibitors of Syk catalytic site bearing therapeutic potential have been developed, with however limited specificity towards Syk. To address this topic, we opted for the design of drug-like compounds that could impede the interaction of Syk with its cellular partners while maintaining an active kinase protein. To achieve this challenging task, we used the powerful potential of intracellular antibodies for the modulation of cellular functions in vivo, combined to structure-based in silico screening. In our previous studies, we reported the anti-allergic properties of the intracellular antibody G4G11. With the aim of finding functional mimics of G4G11, we developed an Antibody Displacement Assay and we isolated the drug-like compound C-13, with promising in vivo anti-allergic activity. The likely binding cavity of this compound is located at the close vicinity of G4G11 epitope, far away from the catalytic site of Syk. Here we report the virtual screen of a collection of 500,000 molecules against this new cavity, which led to the isolation of 1000 compounds subsequently evaluated for their in vitro inhibitory effects using the Antibody Displacement Assay. Eighty five compounds were selected and evaluated for their ability to inhibit the liberation of allergic mediators from mast cells. Among them, 10 compounds inhibited degranulation with IC50 values ≤10 µM. The most bioactive compounds combine biological activity, significant inhibition of antibody binding and strong affinity for Syk. Moreover, these molecules show a good potential for oral bioavailability and are not kinase catalytic site inhibitors. These bioactive compounds could be used as starting points for the development of new classes of non-enzymatic inhibitors of Syk and for drug discovery endeavour in the field of inflammation related disorders

The Human Phenotype Ontology in 2024: phenotypes around the world.

The Human Phenotype Ontology (HPO) is a widely used resource that comprehensively organizes and defines the phenotypic features of human disease, enabling computational inference and supporting genomic and phenotypic analyses through semantic similarity and machine learning algorithms. The HPO has widespread applications in clinical diagnostics and translational research, including genomic diagnostics, gene-disease discovery, and cohort analytics. In recent years, groups around the world have developed translations of the HPO from English to other languages, and the HPO browser has been internationalized, allowing users to view HPO term labels and in many cases synonyms and definitions in ten languages in addition to English. Since our last report, a total of 2239 new HPO terms and 49235 new HPO annotations were developed, many in collaboration with external groups in the fields of psychiatry, arthrogryposis, immunology and cardiology. The Medical Action Ontology (MAxO) is a new effort to model treatments and other measures taken for clinical management. Finally, the HPO consortium is contributing to efforts to integrate the HPO and the GA4GH Phenopacket Schema into electronic health records (EHRs) with the goal of more standardized and computable integration of rare disease data in EHRs

Solve-RD: systematic pan-European data sharing and collaborative analysis to solve rare diseases.

For the first time in Europe hundreds of rare disease (RD) experts team up to actively share and jointly analyse existing patient's data. Solve-RD is a Horizon 2020-supported EU flagship project bringing together >300 clinicians, scientists, and patient representatives of 51 sites from 15 countries. Solve-RD is built upon a core group of four European Reference Networks (ERNs; ERN-ITHACA, ERN-RND, ERN-Euro NMD, ERN-GENTURIS) which annually see more than 270,000 RD patients with respective pathologies. The main ambition is to solve unsolved rare diseases for which a molecular cause is not yet known. This is achieved through an innovative clinical research environment that introduces novel ways to organise expertise and data. Two major approaches are being pursued (i) massive data re-analysis of >19,000 unsolved rare disease patients and (ii) novel combined -omics approaches. The minimum requirement to be eligible for the analysis activities is an inconclusive exome that can be shared with controlled access. The first preliminary data re-analysis has already diagnosed 255 cases form 8393 exomes/genome datasets. This unprecedented degree of collaboration focused on sharing of data and expertise shall identify many new disease genes and enable diagnosis of many so far undiagnosed patients from all over Europe

Twist exome capture allows for lower average sequence coverage in clinical exome sequencing

Background Exome and genome sequencing are the predominant techniques in the diagnosis and research of genetic disorders. Sufficient, uniform and reproducible/consistent sequence coverage is a main determinant for the sensitivity to detect single-nucleotide (SNVs) and copy number variants (CNVs). Here we compared the ability to obtain comprehensive exome coverage for recent exome capture kits and genome sequencing techniques. Results We compared three different widely used enrichment kits (Agilent SureSelect Human All Exon V5, Agilent SureSelect Human All Exon V7 and Twist Bioscience) as well as short-read and long-read WGS. We show that the Twist exome capture significantly improves complete coverage and coverage uniformity across coding regions compared to other exome capture kits. Twist performance is comparable to that of both short- and long-read whole genome sequencing. Additionally, we show that even at a reduced average coverage of 70× there is only minimal loss in sensitivity for SNV and CNV detection. Conclusion We conclude that exome sequencing with Twist represents a significant improvement and could be performed at lower sequence coverage compared to other exome capture techniques