1
|
Basit A, Mishra RK, Bandyopadhyay P. Calcium ion binding to calmodulin: binding free energy calculation using the molecular mechanics Poisson-Boltzmann surface area (MM-PBSA) method by incorporating implicit polarization. J Biomol Struct Dyn 2020; 39:7213-7222. [PMID: 32835611 DOI: 10.1080/07391102.2020.1810125] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/04/2023]
Abstract
Binding of calcium ion to calcium-binding proteins (CBP) triggers a large number of biological processes in a cell. CBP are known to play important roles in various diseases, such as cancer, alzheimer, and neuronal problems. However, the calculation of the binding affinity of calcium ion to CBP still possesses a significant challenge to the computational investigators. One of the main reasons for this difficulty is the polarization of CBP due to the binding of calcium. In the current work, we have used the implicit polarization method of Leontyev et al. (PCCP, 13.7 (2011): 2613-2626) to calculate the binding free energy of calcium ion binding to calmodulin, an important CBP. We have used the widely used MM-PBSA method to find a good protocol of calculation with implicit polarization. We have also optimized the best value of the calcium radius to match the experimental results. Our results show incorporation of polarization improves the agreement between the calculated and experimental results, although still, some discrepancy remains. On the whole, this work shows implicit polarization when combined with the MM-PBSA method can give results better than calculation without any polarization, and further improvement is necessary to get a quantitative match with experiments.Communicated by Ramaswamy H. Sarma.
Collapse
Affiliation(s)
- Abdul Basit
- School of computational and integrative sciences, Jawaharlal Nehru University, New Delhi, India
| | - Rakesh Kumar Mishra
- School of computational and integrative sciences, Jawaharlal Nehru University, New Delhi, India
| | - Pradipta Bandyopadhyay
- School of computational and integrative sciences, Jawaharlal Nehru University, New Delhi, India
| |
Collapse
|
2
|
Miteva MA, Villoutreix BO. Computational Biology and Chemistry in MTi: Emphasis on the Prediction of Some ADMET Properties. Mol Inform 2017; 36. [DOI: 10.1002/minf.201700008] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/17/2017] [Accepted: 02/03/2017] [Indexed: 12/21/2022]
Affiliation(s)
- Maria A. Miteva
- Université Paris Diderot, Sorbonne Paris Cité, Molécules Thérapeutiques In Silico , Inserm UMR−S 973; 35 rue Helene Brion 75013 Paris France
- INSERM, U973; F-75205 Paris France
| | - Bruno O. Villoutreix
- Université Paris Diderot, Sorbonne Paris Cité, Molécules Thérapeutiques In Silico , Inserm UMR−S 973; 35 rue Helene Brion 75013 Paris France
- INSERM, U973; F-75205 Paris France
| |
Collapse
|
3
|
Offutt TL, Swift RV, Amaro RE. Enhancing Virtual Screening Performance of Protein Kinases with Molecular Dynamics Simulations. J Chem Inf Model 2016; 56:1923-1935. [PMID: 27662181 DOI: 10.1021/acs.jcim.6b00261] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022]
Abstract
In silico virtual screening (VS) is a powerful hit identification technique used in drug discovery projects that aims to effectively distinguish true actives from inactive or decoy molecules. To better capture the dynamic behavior of protein drug targets, compound databases may be screened against an ensemble of protein conformations, which may be experimentally determined or generated computationally, i.e. via molecular dynamics (MD) simulations. Several studies have shown that conformations generated by MD are useful in identifying novel hit compounds, in part because structural rearrangements sampled during MD can provide novel targetable areas. However, it remains difficult to predict a priori when an MD conformation will outperform a VS against the crystal structure alone. Here, we assess whether MD conformations result in improved VS performance for six protein kinases. MD conformations are selected using three different methods, and their VS performances are compared to the corresponding crystal structures. Additionally, these conformations are used to train ensembles, and their VS performance is compared to the individual MD conformations and the corresponding crystal structures using receiver operating characteristic curve (ROC) metrics. We show that performing MD results in at least one conformation that offers better VS performance than the crystal structure, and that, while it is possible to train ensembles to outperform the crystal structure alone, the extent of this enhancement is target dependent. Lastly, we show that the optimal structural selection method is also target dependent and recommend optimizing virtual screens on a kinase-by-kinase basis to improve the likelihood of success.
Collapse
Affiliation(s)
- Tavina L Offutt
- Department of Chemistry and Biochemistry, University of California, San Diego , 9500 Gilman Drive, La Jolla, California 92092-0340, United States
| | - Robert V Swift
- Department of Chemistry and Biochemistry, University of California, San Diego , 9500 Gilman Drive, La Jolla, California 92092-0340, United States
| | - Rommie E Amaro
- Department of Chemistry and Biochemistry, University of California, San Diego , 9500 Gilman Drive, La Jolla, California 92092-0340, United States
| |
Collapse
|
4
|
Sampling of conformational ensemble for virtual screening using molecular dynamics simulations and normal mode analysis. Future Med Chem 2015; 7:2317-31. [DOI: 10.4155/fmc.15.150] [Citation(s) in RCA: 20] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/03/2023] Open
Abstract
Aim: Molecular dynamics simulations and normal mode analysis are well-established approaches to generate receptor conformational ensembles (RCEs) for ligand docking and virtual screening. Here, we report new fast molecular dynamics-based and normal mode analysis-based protocols combined with conformational pocket classifications to efficiently generate RCEs. Materials & Methods: We assessed our protocols on two well-characterized protein targets showing local active site flexibility, dihydrofolate reductase and large collective movements, CDK2. The performance of the RCEs was validated by distinguishing known ligands of dihydrofolate reductase and CDK2 among a dataset of diverse chemical decoys. Results & discussion: Our results show that different simulation protocols can be efficient for generation of RCEs depending on different kind of protein flexibility.[Formula: see text]
Collapse
|
5
|
Surfing the Protein-Protein Interaction Surface Using Docking Methods: Application to the Design of PPI Inhibitors. Molecules 2015; 20:11569-603. [PMID: 26111183 PMCID: PMC6272567 DOI: 10.3390/molecules200611569] [Citation(s) in RCA: 44] [Impact Index Per Article: 4.9] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/27/2015] [Revised: 06/02/2015] [Accepted: 06/15/2015] [Indexed: 02/06/2023] Open
Abstract
Blocking protein-protein interactions (PPI) using small molecules or peptides modulates biochemical pathways and has therapeutic significance. PPI inhibition for designing drug-like molecules is a new area that has been explored extensively during the last decade. Considering the number of available PPI inhibitor databases and the limited number of 3D structures available for proteins, docking and scoring methods play a major role in designing PPI inhibitors as well as stabilizers. Docking methods are used in the design of PPI inhibitors at several stages of finding a lead compound, including modeling the protein complex, screening for hot spots on the protein-protein interaction interface and screening small molecules or peptides that bind to the PPI interface. There are three major challenges to the use of docking on the relatively flat surfaces of PPI. In this review we will provide some examples of the use of docking in PPI inhibitor design as well as its limitations. The combination of experimental and docking methods with improved scoring function has thus far resulted in few success stories of PPI inhibitors for therapeutic purposes. Docking algorithms used for PPI are in the early stages, however, and as more data are available docking will become a highly promising area in the design of PPI inhibitors or stabilizers.
Collapse
|
6
|
Xu D, Wang B, Meroueh SO. Structure-based computational approaches for small-molecule modulation of protein-protein interactions. Methods Mol Biol 2015; 1278:77-92. [PMID: 25859944 DOI: 10.1007/978-1-4939-2425-7_5] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/04/2023]
Abstract
Three-dimensional structures of proteins offer an opportunity for the rational design of small molecules to modulate protein-protein interactions. The presence of a well-defined binding pocket on the surface of protein complexes, particularly at their interface, can be used for docking-based virtual screening of chemical libraries. Several approaches have been developed to identify binding pockets that are implemented in programs such as SiteMap, fpocket, and FTSite. These programs enable the scoring of these pockets to determine whether they are suitable to accommodate high-affinity small molecules. Virtual screening of commercial or combinatorial libraries can be carried out to enrich these libraries and select compounds for further experimental validation. In virtual screening, a compound library is docked to the target protein. The resulting structures are scored and ranked for the selection and experimental validation of top candidates. Molecular docking has been implemented in a number of computer programs such as AutoDock Vina. We select a set of protein-protein interactions that have been successfully inhibited with small molecules in the past. Several computer programs are applied to identify pockets on the surface, and molecular docking is conducted in an attempt to reproduce the binding pose of the inhibitors. The results highlight the strengths and limitations of computational methods for the design of PPI inhibitors.
Collapse
Affiliation(s)
- David Xu
- Center for Computational Biology and Bioinformatics, Indiana University School of Medicine, 410 W. 10th Street, Indianapolis, IN, 46202, USA
| | | | | |
Collapse
|
7
|
Wang B, Buchman CD, Li L, Hurley TD, Meroueh SO. Enrichment of chemical libraries docked to protein conformational ensembles and application to aldehyde dehydrogenase 2. J Chem Inf Model 2014; 54:2105-16. [PMID: 24856086 PMCID: PMC4114474 DOI: 10.1021/ci5002026] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/30/2022]
Abstract
Molecular recognition is a complex process that involves a large ensemble of structures of the receptor and ligand. Yet, most structure-based virtual screening is carried out on a single structure typically from X-ray crystallography. Explicit-solvent molecular dynamics (MD) simulations offer an opportunity to sample multiple conformational states of a protein. Here we evaluate our recently developed scoring method SVMSP in its ability to enrich chemical libraries docked to MD structures of seven proteins from the Directory of Useful Decoys (DUD). SVMSP is a target-specific rescoring method that combines machine learning with statistical potentials. We find that enrichment power as measured by the area under the ROC curve (ROC-AUC) is not affected by increasing the number of MD structures. Among individual MD snapshots, many exhibited enrichment that was significantly better than the crystal structure, but no correlation between enrichment and structural deviation from crystal structure was found. We followed an innovative approach by training SVMSP scoring models using MD structures (SVMSPMD). The resulting models were applied to two difficult cases (p38 and CDK2) for which enrichment was not better than random. We found remarkable increase in enrichment power, particularly for p38, where the ROC-AUC increased by 0.30 to 0.85. Finally, we explored approaches for a priori identification of MD snapshots with high enrichment power from an MD simulation in the absence of active compounds. We found that the use of randomly selected compounds docked to the target of interest using SVMSP led to notable enrichment for EGFR and Src MD snapshots. SVMSP rescoring of protein-compound MD structures was applied for the search of small-molecule inhibitors of the mitochondrial enzyme aldehyde dehydrogenase 2 (ALDH2). Rank-ordering of a commercial library of 50 000 compounds docked to MD structures of ALDH2 led to five small-molecule inhibitors. Four compounds had IC50s below 5 μM. These compounds serve as leads for the design and synthesis of more potent and selective ALDH2 inhibitors.
Collapse
Affiliation(s)
- Bo Wang
- Department of Biochemistry and Molecular Biology, ‡Melvin and Bren Simon Cancer Center, §Center for Computational Biology and Bioinformatics, and ⊥Stark Neurosciences Institute, Indiana University School of Medicine , 535 Barnhill Drive, Indianapolis, Indiana 46202, United States
| | | | | | | | | |
Collapse
|
8
|
Jamroz M, Kolinski A, Kmiecik S. CABS-flex predictions of protein flexibility compared with NMR ensembles. ACTA ACUST UNITED AC 2014; 30:2150-4. [PMID: 24735558 PMCID: PMC4103595 DOI: 10.1093/bioinformatics/btu184] [Citation(s) in RCA: 58] [Impact Index Per Article: 5.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/23/2022]
Abstract
Motivation: Identification of flexible regions of protein structures is important for understanding of their biological functions. Recently, we have developed a fast approach for predicting protein structure fluctuations from a single protein model: the CABS-flex. CABS-flex was shown to be an efficient alternative to conventional all-atom molecular dynamics (MD). In this work, we evaluate CABS-flex and MD predictions by comparison with protein structural variations within NMR ensembles. Results: Based on a benchmark set of 140 proteins, we show that the relative fluctuations of protein residues obtained from CABS-flex are well correlated to those of NMR ensembles. On average, this correlation is stronger than that between MD and NMR ensembles. In conclusion, CABS-flex is useful and complementary to MD in predicting protein regions that undergo conformational changes as well as the extent of such changes. Availability and implementation: The CABS-flex is freely available to all users at http://biocomp.chem.uw.edu.pl/CABSflex. Contact: sekmi@chem.uw.edu.pl Supplementary information:Supplementary data are available at Bioinformatics online.
Collapse
Affiliation(s)
- Michal Jamroz
- Laboratory of Theory of Biopolymers, Faculty of Chemistry, University of Warsaw, Pasteura 1, Warsaw 02-093, Poland
| | - Andrzej Kolinski
- Laboratory of Theory of Biopolymers, Faculty of Chemistry, University of Warsaw, Pasteura 1, Warsaw 02-093, Poland
| | - Sebastian Kmiecik
- Laboratory of Theory of Biopolymers, Faculty of Chemistry, University of Warsaw, Pasteura 1, Warsaw 02-093, Poland
| |
Collapse
|
9
|
In silico mechanistic profiling to probe small molecule binding to sulfotransferases. PLoS One 2013; 8:e73587. [PMID: 24039991 PMCID: PMC3765257 DOI: 10.1371/journal.pone.0073587] [Citation(s) in RCA: 20] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/21/2013] [Accepted: 07/28/2013] [Indexed: 01/01/2023] Open
Abstract
Drug metabolizing enzymes play a key role in the metabolism, elimination and detoxification of xenobiotics, drugs and endogenous molecules. While their principal role is to detoxify organisms by modifying compounds, such as pollutants or drugs, for a rapid excretion, in some cases they render their substrates more toxic thereby inducing severe side effects and adverse drug reactions, or their inhibition can lead to drug–drug interactions. We focus on sulfotransferases (SULTs), a family of phase II metabolizing enzymes, acting on a large number of drugs and hormones and showing important structural flexibility. Here we report a novel in silico structure-based approach to probe ligand binding to SULTs. We explored the flexibility of SULTs by molecular dynamics (MD) simulations in order to identify the most suitable multiple receptor conformations for ligand binding prediction. Then, we employed structure-based docking-scoring approach to predict ligand binding and finally we combined the predicted interaction energies by using a QSAR methodology. The results showed that our protocol successfully prioritizes potent binders for the studied here SULT1 isoforms, and give new insights on specific molecular mechanisms for diverse ligands’ binding related to their binding sites plasticity. Our best QSAR models, introducing predicted protein-ligand interaction energy by using docking, showed accuracy of 67.28%, 78.00% and 75.46%, for the isoforms SULT1A1, SULT1A3 and SULT1E1, respectively. To the best of our knowledge our protocol is the first in silico structure-based approach consisting of a protein-ligand interaction analysis at atomic level that considers both ligand and enzyme flexibility, along with a QSAR approach, to identify small molecules that can interact with II phase dug metabolizing enzymes.
Collapse
|
10
|
Isvoran A, Craciun D, Martiny V, Sperandio O, Miteva MA. Computational analysis of protein-protein interfaces involving an alpha helix: insights for terphenyl-like molecules binding. BMC Pharmacol Toxicol 2013; 14:31. [PMID: 23768251 PMCID: PMC3689098 DOI: 10.1186/2050-6511-14-31] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/23/2013] [Accepted: 06/11/2013] [Indexed: 12/23/2022] Open
Abstract
Background Protein-Protein Interactions (PPIs) are key for many cellular processes. The characterization of PPI interfaces and the prediction of putative ligand binding sites and hot spot residues are essential to design efficient small-molecule modulators of PPI. Terphenyl and its derivatives are small organic molecules known to mimic one face of protein-binding alpha-helical peptides. In this work we focus on several PPIs mediated by alpha-helical peptides. Method We performed computational sequence- and structure-based analyses in order to evaluate several key physicochemical and surface properties of proteins known to interact with alpha-helical peptides and/or terphenyl and its derivatives. Results Sequence-based analysis revealed low sequence identity between some of the analyzed proteins binding alpha-helical peptides. Structure-based analysis was performed to calculate the volume, the fractal dimension roughness and the hydrophobicity of the binding regions. Besides the overall hydrophobic character of the binding pockets, some specificities were detected. We showed that the hydrophobicity is not uniformly distributed in different alpha-helix binding pockets that can help to identify key hydrophobic hot spots. Conclusions The presence of hydrophobic cavities at the protein surface with a more complex shape than the entire protein surface seems to be an important property related to the ability of proteins to bind alpha-helical peptides and low molecular weight mimetics. Characterization of similarities and specificities of PPI binding sites can be helpful for further development of small molecules targeting alpha-helix binding proteins.
Collapse
Affiliation(s)
- Adriana Isvoran
- Department of Biology and Chemistry, West University of Timisoara, 16 Pestalozzi, Timisoara 300115, Romania
| | | | | | | | | |
Collapse
|
11
|
Connors CR, Rosenman DJ, Lopes DHJ, Mittal S, Bitan G, Sorci M, Belfort G, Garcia A, Wang C. Tranilast binds to aβ monomers and promotes aβ fibrillation. Biochemistry 2013; 52:3995-4002. [PMID: 23679559 PMCID: PMC4082028 DOI: 10.1021/bi400426t] [Citation(s) in RCA: 10] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/13/2023]
Abstract
The antiallergy and potential anticancer drug tranilast has been patented for treating Alzheimer's disease (AD), in which amyloid β-protein (Aβ) plays a key pathogenic role. We used solution NMR to determine that tranilast binds to Aβ40 monomers with ∼300 μM affinity. Remarkably, tranilast increases Aβ40 fibrillation more than 20-fold in the thioflavin T assay at a 1:1 molar ratio, as well as significantly reducing the lag time. Tranilast likely promotes fibrillation by shifting Aβ monomer conformations to those capable of seed formation and fibril elongation. Molecular docking results qualitatively agree with NMR chemical shift perturbation, which together indicate that hydrophobic interactions are the major driving force of the Aβ-tranilast interaction. These data suggest that AD may be a potential complication for tranilast usage in elderly patients.
Collapse
Affiliation(s)
- Christopher R. Connors
- Center for Biotechnology and Interdisciplinary Studies, Rensselaer Polytechnic Institute (RPI), Troy, New York 12180, United States
- Graduate Program in Biochemistry and Biophysics, Rensselaer Polytechnic Institute (RPI), Troy, New York 12180, United States
| | - David J. Rosenman
- Center for Biotechnology and Interdisciplinary Studies, Rensselaer Polytechnic Institute (RPI), Troy, New York 12180, United States
- Department of Biology, Rensselaer Polytechnic Institute (RPI), Troy, New York 12180, United States
| | - Dahabada H. J. Lopes
- Department of Neurology, David Geffen School of Medicine, University of California at Los Angeles, Los Angeles, California 90095, United States
| | - Shivina Mittal
- Department of Neurology, David Geffen School of Medicine, University of California at Los Angeles, Los Angeles, California 90095, United States
| | - Gal Bitan
- Department of Neurology, David Geffen School of Medicine, University of California at Los Angeles, Los Angeles, California 90095, United States
- Brain Research Institute, University of California at Los Angeles, Los Angeles, California 90095, United States
- Molecular Biology Institute, University of California at Los Angeles, Los Angeles, California 90095, United States
| | - Mirco Sorci
- Center for Biotechnology and Interdisciplinary Studies, Rensselaer Polytechnic Institute (RPI), Troy, New York 12180, United States
- Department of Chemical and Biological Engineering, Rensselaer Polytechnic Institute (RPI), Troy, New York 12180, United States
| | - Georges Belfort
- Center for Biotechnology and Interdisciplinary Studies, Rensselaer Polytechnic Institute (RPI), Troy, New York 12180, United States
- Department of Chemical and Biological Engineering, Rensselaer Polytechnic Institute (RPI), Troy, New York 12180, United States
| | - Angel Garcia
- Center for Biotechnology and Interdisciplinary Studies, Rensselaer Polytechnic Institute (RPI), Troy, New York 12180, United States
- Department of Physics, Rensselaer Polytechnic Institute (RPI), Troy, New York 12180, United States
| | - Chunyu Wang
- Center for Biotechnology and Interdisciplinary Studies, Rensselaer Polytechnic Institute (RPI), Troy, New York 12180, United States
- Graduate Program in Biochemistry and Biophysics, Rensselaer Polytechnic Institute (RPI), Troy, New York 12180, United States
- Department of Biology, Rensselaer Polytechnic Institute (RPI), Troy, New York 12180, United States
| |
Collapse
|
12
|
Moroy G, Martiny VY, Vayer P, Villoutreix BO, Miteva MA. Toward in silico structure-based ADMET prediction in drug discovery. Drug Discov Today 2011; 17:44-55. [PMID: 22056716 DOI: 10.1016/j.drudis.2011.10.023] [Citation(s) in RCA: 167] [Impact Index Per Article: 12.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/27/2011] [Revised: 10/07/2011] [Accepted: 10/21/2011] [Indexed: 12/12/2022]
Abstract
Quantitative structure-activity relationship (QSAR) methods and related approaches have been used to investigate the molecular features that influence the absorption, distribution, metabolism, excretion and toxicity (ADMET) of drugs. As the three-dimensional structures of several major ADMET proteins become available, structure-based (docking-scoring) computations can be carried out to complement or to go beyond QSAR studies. Applying docking-scoring methods to ADMET proteins is a challenging process because they usually have a large and flexible binding cavity; however, promising results relating to metabolizing enzymes have been reported. After reviewing current trends in the field we applied structure-based methods in the context of receptor flexibility in a case study involving the phase II metabolizing sulfotransferases. Overall, the explored concepts and results suggested that structure-based ADMET profiling will probably join the mainstream during the coming years.
Collapse
Affiliation(s)
- Gautier Moroy
- Inserm UMR-S 973, Molécules Thérapeutiques In Silico, Université Paris Diderot, Sorbonne Paris Cité, 35 Rue Helene Brion, 75013 Paris, France
| | | | | | | | | |
Collapse
|