RESUMEN
We present a neural-network-based high-throughput molecular conformer-generation algorithm. A chemical graph-convolutional network is trained to predict low-energy conformers in internal coordinate representation (bond lengths, bond, and torsion angles), starting from two-dimensional (2D) chemical topology. Generative neural network (NN) architecture performs denoising from torsion space, producing conformer ensembles with populations that are well correlated with torsion energy profiles. Short force-field-based energy minimization is applied to refine final conformers. All computation-intensive stages of the algorithm are GPU-optimized. The procedure (termed GINGER) is benchmarked on a commonly used test set of bioactive three-dimensional (3D) conformers from the PDB. We demonstrate highly competitive results in conformer recovery and throughput rates suitable for giga-scale compound library processing. A web server that allows interactive conformer ensemble generation by GINGER and their viewing is made freely available at https://www.molsoft.com/gingerdemo.html.
RESUMEN
We present a graph-convolutional neural network (GCNN)-based method for learning and prediction of statistical torsional profiles (STP) in small organic molecules based on the experimental X-ray structure data. A specialized GCNN torsion profile model is trained using the structures in the Crystallography Open Database (COD). The GCNN-STP model captures torsional preferences over a wide range of torsion rotor chemotypes and correctly predicts a variety of effects from the vicinal atoms and moieties. GCNN-STP statistical profiles also show good agreement with quantum chemically (DFT) calculated torsion energy profiles. Furthermore, we demonstrate the application of the GCNN-STP statistical profiles for conformer generation. A web server that allows interactive profile prediction and viewing is made freely available at https://www.molsoft.com/tortool.html.
Asunto(s)
Redes Neurales de la Computación , Cristalografía , Bases de Datos FactualesRESUMEN
Communication of data and ideas within a medicinal chemistry project on a global as well as local level is a crucial aspect in the drug design cycle. Over a time frame of eight years, we built and optimized FOCUS, a platform to produce, visualize, and share information on various aspects of a drug discovery project such as cheminformatics, data analysis, structural information, and design. FOCUS is tightly integrated with internal services that involve-among others-data retrieval systems and in-silico models and provides easy access to automated modeling procedures such as pharmacophore searches, R-group analysis, and similarity searches. In addition, an interactive 3D editor was developed to assist users in the generation and docking of close analogues of a known lead. In this paper, we will specifically concentrate on issues we faced during development, deployment, and maintenance of the software and how we continually adapted the software in order to improve usability. We will provide usage examples to highlight the functionality as well as limitations of FOCUS at the various stages of the development process. We aim to make the discussion as independent of the software platform as possible, so that our experiences can be of more general value to the drug discovery community.
Asunto(s)
Química Farmacéutica/métodos , Comunicación , Simulación por Computador , Descubrimiento de Drogas/métodos , Biología Computacional , LigandosRESUMEN
The function of a protein is determined by its intrinsic activity in the context of its subcellular distribution. Membranes localize proteins within cellular compartments and govern their specific activities. Discovering such membrane-protein interactions is important for understanding biological mechanisms and could uncover novel sites for therapeutic intervention. We present a method for detecting membrane interactive proteins and their exposed residues that insert into lipid bilayers. Although the development process involved analysis of how C1b, C2, ENTH, FYVE, Gla, pleckstrin homology (PH), and PX domains bind membranes, the resulting membrane optimal docking area (MODA) method yields predictions for a given protein of known three-dimensional structures without referring to canonical membrane-targeting modules. This approach was tested on the Arf1 GTPase, ATF2 acetyltransferase, von Willebrand factor A3 domain, and Neisseria gonorrhoeae MsrB protein and further refined with membrane interactive and non-interactive FAPP1 and PKD1 pleckstrin homology domains, respectively. Furthermore we demonstrate how this tool can be used to discover unprecedented membrane binding functions as illustrated by the Bro1 domain of Alix, which was revealed to recognize lysobisphosphatidic acid (LBPA). Validation of novel membrane-protein interactions relies on other techniques such as nuclear magnetic resonance spectroscopy (NMR), which was used here to map the sites of micelle interaction. Together this indicates that genome-wide identification of known and novel membrane interactive proteins and sites is now feasible and provides a new tool for functional annotation of the proteome.
Asunto(s)
Membrana Celular/química , Proteínas de la Membrana/química , Anotación de Secuencia Molecular/métodos , Análisis de Secuencia de Proteína/métodos , Proteínas Bacterianas/química , Proteínas Bacterianas/genética , Proteínas Bacterianas/metabolismo , Membrana Celular/genética , Membrana Celular/metabolismo , Humanos , Proteínas de la Membrana/genética , Proteínas de la Membrana/metabolismo , Neisseria gonorrhoeae , Estructura Terciaria de Proteína , Proteoma/química , Proteoma/genética , Proteoma/metabolismoAsunto(s)
Biología Molecular , Publicaciones Periódicas como Asunto/tendencias , Conformación Proteica , Edición/tendencias , Interfaz Usuario-Computador , Algoritmos , Bases de Datos de Proteínas , Ensayos Analíticos de Alto Rendimiento , Humanos , Modelos Moleculares , Biología Molecular/métodos , Biología Molecular/tendencias , Programas InformáticosRESUMEN
SUMMARY: SimiCon is a web server designed for an automated identification of equivalent protein-ligand atomic contacts in different conformational models of a complex. The contacts are computed with internal coordinate mechanics (ICM) software with respect to molecular symmetry and the results are shown in the browser as text, tables and interactive 3D graphics. The web server can be executed remotely without a browser to allow users to automate multiple calculations. AVAILABILITY: SimiCon is freely available at http://abagyan.ucsd.edu/SimiCon
Asunto(s)
Proteínas/química , Programas Informáticos , Gráficos por Computador , Bases de Datos de Proteínas , Internet , Ligandos , Conformación ProteicaRESUMEN
Essential for viral replication and highly conserved among poxviridae, the vaccinia virus I7L ubiquitin-like proteinase (ULP) is an attractive target for development of smallpox antiviral drugs. At the same time, the I7L proteinase exemplifies several interesting challenges from the rational drug design perspective. In the absence of a published I7L X-ray structure, we have built a detailed 3D model of the I7L ligand binding site (S2-S2' pocket) based on exceptionally high structural conservation of this site in proteases of the ULP family. The accuracy and limitations of this model were assessed through comparative analysis of available X-ray structures of ULPs, as well as energy based conformational modeling. The 3D model of the I7L ligand binding site was used to perform covalent docking and VLS of a comprehensive library of about 230,000 available ketone and aldehyde compounds. Out of 456 predicted ligands, 97 inhibitors of I7L proteinase activity were confirmed in biochemical assays ( approximately 20% overall hit rate). These experimental results both validate our I7L ligand binding model and provide initial leads for rational optimization of poxvirus I7L proteinase inhibitors. Thus, fragments predicted to bind in the prime portion of the active site can be combined with fragments on non-prime side to yield compounds with improved activity and specificity.
Asunto(s)
Antivirales/química , Antivirales/farmacología , Cisteína Endopeptidasas/química , Inhibidores de Cisteína Proteinasa/química , Inhibidores de Cisteína Proteinasa/farmacología , Ubiquitinas/antagonistas & inhibidores , Virus Vaccinia/enzimología , Secuencia de Aminoácidos , Sitios de Unión , Simulación por Computador , Cisteína Endopeptidasas/genética , Diseño de Fármacos , Evaluación Preclínica de Medicamentos , Cetonas/química , Ligandos , Modelos Moleculares , Datos de Secuencia Molecular , Poxviridae/efectos de los fármacos , Poxviridae/enzimología , Poxviridae/genética , Homología de Secuencia de Aminoácido , Relación Estructura-Actividad , Especificidad por Sustrato , Ubiquitinas/química , Ubiquitinas/genética , Interfaz Usuario-Computador , Virus Vaccinia/efectos de los fármacos , Virus Vaccinia/genéticaRESUMEN
Recent advances in structural proteomics call for development of fast and reliable automatic methods for prediction of functional surfaces of proteins with known three-dimensional structure, including binding sites for known and unknown protein partners as well as oligomerization interfaces. Despite significant progress the problem is still far from being solved. Most existing methods rely, at least partially, on evolutionary information from multiple sequence alignments projected on protein surface. The common drawback of such methods is their limited applicability to the proteins with a sparse set of sequential homologs, as well as inability to detect interfaces in evolutionary variable regions. In this study, the authors developed an improved method for predicting interfaces from a single protein structure, which is based on local statistical properties of the protein surface derived at the level of atomic groups. The proposed Protein IntErface Recognition (PIER) method achieved the overall precision of 60% at the recall threshold of 50% at the residue level on a diverse benchmark of 490 homodimeric, 62 heterodimeric, and 196 transient interfaces (compared with 25% precision at 50% recall expected from random residue function assignment). For 70% of proteins in the benchmark, the binding patch residues were successfully detected with precision exceeding 50% at 50% recall. The calculation only took seconds for an average 300-residue protein. The authors demonstrated that adding the evolutionary conservation signal only marginally influenced the overall prediction performance on the benchmark; moreover, for certain classes of proteins, using this signal actually resulted in a deteriorated prediction. Thorough benchmarking using other datasets from literature showed that PIER yielded improved performance as compared with several alignment-free or alignment-dependent predictions. The accuracy, efficiency, and dependence on structure alone make PIER a suitable tool for automated high-throughput annotation of protein structures emerging from structural proteomics projects.