1
|
Zsidó BZ, Bayarsaikhan B, Börzsei R, Szél V, Mohos V, Hetényi C. The Advances and Limitations of the Determination and Applications of Water Structure in Molecular Engineering. Int J Mol Sci 2023; 24:11784. [PMID: 37511543 PMCID: PMC10381018 DOI: 10.3390/ijms241411784] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/20/2023] [Revised: 07/18/2023] [Accepted: 07/20/2023] [Indexed: 07/30/2023] Open
Abstract
Water is a key actor of various processes of nature and, therefore, molecular engineering has to take the structural and energetic consequences of hydration into account. While the present review focuses on the target-ligand interactions in drug design, with a focus on biomolecules, these methods and applications can be easily adapted to other fields of the molecular engineering of molecular complexes, including solid hydrates. The review starts with the problems and solutions of the determination of water structures. The experimental approaches and theoretical calculations are summarized, including conceptual classifications. The implementations and applications of water models are featured for the calculation of the binding thermodynamics and computational ligand docking. It is concluded that theoretical approaches not only reproduce or complete experimental water structures, but also provide key information on the contribution of individual water molecules and are indispensable tools in molecular engineering.
Collapse
Affiliation(s)
- Balázs Zoltán Zsidó
- Department of Pharmacology and Pharmacotherapy, Medical School, University of Pécs, Szigeti út 12, 7624 Pécs, Hungary
| | - Bayartsetseg Bayarsaikhan
- Department of Pharmacology and Pharmacotherapy, Medical School, University of Pécs, Szigeti út 12, 7624 Pécs, Hungary
| | - Rita Börzsei
- Department of Pharmacology and Pharmacotherapy, Medical School, University of Pécs, Szigeti út 12, 7624 Pécs, Hungary
| | - Viktor Szél
- Department of Pharmacology and Pharmacotherapy, Medical School, University of Pécs, Szigeti út 12, 7624 Pécs, Hungary
| | - Violetta Mohos
- Department of Pharmacology and Pharmacotherapy, Medical School, University of Pécs, Szigeti út 12, 7624 Pécs, Hungary
| | - Csaba Hetényi
- Department of Pharmacology and Pharmacotherapy, Medical School, University of Pécs, Szigeti út 12, 7624 Pécs, Hungary
| |
Collapse
|
2
|
Eberhardt J, Forli S. WaterKit: Thermodynamic Profiling of Protein Hydration Sites. J Chem Theory Comput 2023; 19:2535-2556. [PMID: 37094087 PMCID: PMC10732097 DOI: 10.1021/acs.jctc.2c01087] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/26/2023]
Abstract
Water desolvation is one of the key components of the free energy of binding of small molecules to their receptors. Thus, understanding the energetic balance of solvation and desolvation resulting from individual water molecules can be crucial when estimating ligand binding, especially when evaluating different molecules and poses as done in High-Throughput Virtual Screening (HTVS). Over the most recent decades, several methods were developed to tackle this problem, ranging from fast approximate methods (usually empirical functions using either discrete atom-atom pairwise interactions or continuum solvent models) to more computationally expensive and accurate ones, mostly based on Molecular Dynamics (MD) simulations, such as Grid Inhomogeneous Solvation Theory (GIST) or Double Decoupling. On one hand, MD-based methods are prohibitive to use in HTVS to estimate the role of waters on the fly for each ligand. On the other hand, fast and approximate methods show an unsatisfactory level of accuracy, with low agreement with results obtained with the more expensive methods. Here we introduce WaterKit, a new grid-based sampling method with explicit water molecules to calculate thermodynamic properties using the GIST method. Our results show that the discrete placement of water molecules is successful in reproducing the position of crystallographic waters with very high accuracy, as well as providing thermodynamic estimates with accuracy comparable to more expensive MD simulations. Unlike these methods, WaterKit can be used to analyze specific regions on the protein surface, (such as the binding site of a receptor), without having to hydrate and simulate the whole receptor structure. The results show the feasibility of a general and fast method to compute thermodynamic properties of water molecules, making it well-suited to be integrated in high-throughput pipelines such as molecular docking.
Collapse
Affiliation(s)
- Jerome Eberhardt
- Department of Integrative Structural and Computational Biology, Scripps Research, La Jolla, California 92037, United States
| | - Stefano Forli
- Department of Integrative Structural and Computational Biology, Scripps Research, La Jolla, California 92037, United States
| |
Collapse
|
3
|
Marie V, Gordon M. Understanding the co-evolutionary molecular mechanisms of resistance in the HIV-1 Gag and protease. J Biomol Struct Dyn 2022; 40:10852-10861. [PMID: 34253143 DOI: 10.1080/07391102.2021.1950569] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/29/2023]
Abstract
Due to high human immunodeficiency virus type 1 (HIV-1) subtype C infections coupled with increasing antiretroviral treatment failure, the elucidation of complex drug resistance mutational patterns arising through protein co-evolution is required. Despite the inclusion of potent protease inhibitors Lopinavir (LPV) and Darunavir (DRV) in second- and third-line therapies, many patients still fail treatment due to the accumulation of mutations in protease (PR) and recently, Gag. To understand the co-evolutionary molecular mechanisms of resistance in the HIV-1 PR and Gag, we performed 100 ns molecular dynamic simulations on multidrug resistant PR's when bound to LPV, DRV or a mutated A431V NC|p1 Gag cleavage site (CS). Here we showed that distinct changes in PR's active site, flap and elbow regions due to several PR resistance mutations (L10F, M46I, I54V, L76V, V82A) were found to alter LPV and DRV drug binding. However, binding was significantly exacerbated when the mutant PRs were bound to the NC|p1 Gag CS. Although A431V was shown to coordinate several residues in PR, the L76V PR mutation was found to have a significant role in substrate recognition. Consequently, a greater binding affinity was observed when the mutated substrate was bound to an L76V-inclusive PR mutant (Gbind: -62.46 ± 5.75 kcal/mol) than without (Gbind: -50.34 ± 6.28 kcal/mol). These data showed that the co-selection of resistance mutations in the enzyme and substrate can simultaneously constrict regions in PR's active site whilst flexing the flaps to allow flexible movement of the substrate and multiple, complex mechanisms of resistance to occur. Communicated by Ramaswamy H. Sarma.
Collapse
Affiliation(s)
- Veronna Marie
- KwaZulu-Natal Research Innovation & Sequencing Platform (KRISP), Nelson R Mandela School of Medicine, University of KwaZulu-Natal, DurbanSouth Africa
| | - Michelle Gordon
- KwaZulu-Natal Research Innovation & Sequencing Platform (KRISP), Nelson R Mandela School of Medicine, University of KwaZulu-Natal, DurbanSouth Africa
| |
Collapse
|
4
|
Sohrab SS, Kamal MA. Screening, Docking, and Molecular Dynamics Study of Natural Compounds as an Anti-HER2 for the Management of Breast Cancer. Life (Basel) 2022; 12:1729. [PMID: 36362883 PMCID: PMC9693058 DOI: 10.3390/life12111729] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/23/2022] [Revised: 10/19/2022] [Accepted: 10/25/2022] [Indexed: 08/27/2023] Open
Abstract
Breast cancer (BC) is one of the most frequent types of cancer that affect women. Human epidermal growth factor receptor-2 (HER2) is responsible for 20% of all BC cases. The use of anti-HER2 natural compounds in the cure of BC that is HER2-positive patients has resulted in significant increases in survival in both early and advanced stages. The findings of in-silico research support the use of ligands as possible HER2 inhibitors, and molecules with high free energy of binding may have considerable anti-BC action, making them candidates for future drug development. The inhibitory activity of selected ligands like ZINC43069427 and ZINC95918662 against HER2 was found to be -11.0 and -8.50 kcal/mol, respectively. The amino acid residues Leu726, Val734, Ala751, Lys753, Thr798, Gly804, Arg849, Leu852, Thr862, and Asp863 were found in common interaction as compared to the control compound Lapatinib. Molecular dynamics study calculations of these selected potent inhibitors were conducted and found to be stable over the 50 ns simulation time in terms of root mean square deviation (RMSD), root-mean square fluctuation (RMSF), radius of gyration (Rg), and solvent accessible surface area (SASA). In addition, there are several parameters such as absorption, distribution, metabolism, and excretion toxicity (ADMET), physicochemical, and drug-likeness that were checked and found in good range to be potential lead-like molecules. Several drug-likeness rules like Lipinski, Ghose, Veber, Egan, and Muegge were checked and found to be positive for these rules. Based on these calculations and different parameters, these top two selected natural compounds can be used as potential candidates for anti-HER2 for the management of BC.
Collapse
Affiliation(s)
- Sayed Sartaj Sohrab
- Special Infectious Agents Unit, King Fahd Medical Research Center, King Abdulaziz University, Jeddah 22254, Saudi Arabia
- Department of Medical Laboratory Sciences, Faculty of Applied Medical Sciences, King Abdulaziz University, Jeddah 22254, Saudi Arabia
| | - Mohammad Amjad Kamal
- Special Infectious Agents Unit, King Fahd Medical Research Center, King Abdulaziz University, Jeddah 22254, Saudi Arabia
- Institutes for Systems Genetics, Frontiers Science Center for Disease-Related Molecular Network, West China Hospital, Sichuan University, Chengdu 610041, China
- King Fahd Medical Research Center, King Abdulaziz University, Jeddah 22254, Saudi Arabia
- Department of Pharmacy, Faculty of Allied Health Sciences, Daffodil International University, Dhaka 1207, Bangladesh
- Enzymoics, 7 Peterlee Place, Novel Global Community Educational Foundation, Hebersham, NSW 2770, Australia
| |
Collapse
|
5
|
A Report on Multi-Target Anti-Inflammatory Properties of Phytoconstituents from Monochoria hastata (Family: Pontederiaceae). Molecules 2021; 26:molecules26237397. [PMID: 34885978 PMCID: PMC8658818 DOI: 10.3390/molecules26237397] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/14/2021] [Revised: 11/29/2021] [Accepted: 12/02/2021] [Indexed: 12/22/2022] Open
Abstract
This study aims to investigate the potential analgesic properties of the crude extract of Monochoria hastata (MH) leaves using in vivo experiments and in silico analysis. The extract, in a dose-dependent manner, exhibited a moderate analgesic property (~54% pain inhibition in acetic acid-induced writhing test), which is significant (** p < 0.001) as compared to the control group. The complex inflammatory mechanism involves diverse pathways and they are inter-connected. Therefore, multiple inflammatory modulator proteins were selected as the target for in silico analysis. Computational analysis suggests that all the selected targets had different degrees of interaction with the phytochemicals from the extract. Rutin (RU), protocatechuic acid (PA), vanillic acid (VA), and ferulic acid (FA) could regulate multiple targets with a robust efficiency. None of the compounds showed selectivity to Cyclooxygenase-2 (COX-2). However, regulation of COX and lipoxygenase (LOX) cascade by PA can reduce non-steroidal analgesic drugs (NSAIDs)-related side effects, including asthma. RU showed robust regulation of cytokine-mediated pathways like RAS/MAPK and PI3K/NF-kB by inhibition of EGFR and IKBα (IKK), which may prevent multi-organ failure due to cytokine storm in several microbial infections, for example, SARS-CoV-2. Further investigation, using in vivo and in vitro experiments, can be conducted to develop multi-target anti-inflammatory drugs using the isolated compounds from the extract.
Collapse
|
6
|
Tong AB, Burch JD, McKay D, Bustamante C, Crackower MA, Wu H. Could AlphaFold revolutionize chemical therapeutics? Nat Struct Mol Biol 2021; 28:771-772. [PMID: 34561631 DOI: 10.1038/s41594-021-00670-x] [Citation(s) in RCA: 18] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/09/2022]
Affiliation(s)
- Alexander B Tong
- Jason L. Choy Laboratory of Single-Molecule Biophysics, Institute for Quantitative Biosciences-QB3, and Chemistry Graduate Group, University of California Berkeley, Berkeley, CA, USA
| | - Jason D Burch
- Ventus Therapeutics, Waltham, MA, USA.,Ventus Therapeutics, Montreal, Canada
| | - Daniel McKay
- Ventus Therapeutics, Waltham, MA, USA.,Ventus Therapeutics, Montreal, Canada
| | - Carlos Bustamante
- Department of Molecular and Cell Biology, Department of Physics, and Department of Chemistry, Howard Hughes Medical Institute, University of California Berkeley, Berkeley, CA, USA
| | - Michael A Crackower
- Ventus Therapeutics, Waltham, MA, USA.,Ventus Therapeutics, Montreal, Canada
| | - Hao Wu
- Department of Biological Chemistry and Molecular Pharmacology, Harvard Medical School, Boston, MA, USA. .,Program in Cellular and Molecular Medicine, Boston Children's Hospital, Boston, MA, USA.
| |
Collapse
|
7
|
Abstract
Molecular docking is one of the most widely used computational tools in structure-based drug design and is critically dependent on accuracy and robustness of the scoring function. In this work, we introduce a new scoring function Lin_F9, which is a linear combination of nine empirical terms, including a unified metal bond term to specifically describe metal-ligand interactions. Parameters in Lin_F9 are obtained with a multistage fitting protocol using explicit water-included structures. For the CASF-2016 benchmark test set, Lin_F9 achieves the top scoring power among all 34 classical scoring functions for both original crystal poses and locally optimized poses with Pearson correlation coefficients (R) of 0.680 and 0.687, respectively. Meanwhile, in comparison with Vina, Lin_F9 achieves consistently better scoring power and ranking power with various types of protein-ligand complex structures that mimic real docking applications, including end-to-end flexible docking for the CASF-2016 benchmark test set using a single or an ensemble of protein receptor structures, as well as for D3R Grand Challenge (GC4) test sets. Lin_F9 has been implemented in a fork of Smina as an optional built-in scoring function that can be used for docking applications as well as for further improvement of scoring functions and docking protocols. Lin_F9 is accessible through https://yzhang.hpc.nyu.edu/Lin_F9/.
Collapse
Affiliation(s)
- Chao Yang
- Department of Chemistry, New York University, New York, New York 10003, United States
| | - Yingkai Zhang
- Department of Chemistry, New York University, New York, New York 10003, United States
- NYU-ECNU Center for Computational Chemistry at NYU Shanghai, Shanghai 200062, China
| |
Collapse
|
8
|
Kothandan R, Rajan CASG, Arjun J, Raj RRM, Syed S. Virtual screening of phytochemical compounds as potential inhibitors against SARS-CoV-2 infection. BENI-SUEF UNIVERSITY JOURNAL OF BASIC AND APPLIED SCIENCES 2021; 10:9. [PMID: 33527080 PMCID: PMC7841026 DOI: 10.1186/s43088-021-00095-x] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/02/2020] [Accepted: 01/04/2021] [Indexed: 12/28/2022] Open
Abstract
Background The present pandemic situation due to coronavirus has led to the search for newer prevention, diagnostic, and treatment methods. The onset of the corona infection in a human results in acute respiratory illness followed by death if not diagnosed and treated with suitable antiretroviral drugs. With the unavailability of the targeted drug treatment, several repurposed drugs are being used for treatment. However, the side-effects of the drugs urges us to move to a search for newer synthetic- or phytochemical-based drugs. The present study investigates the use of various phytochemicals virtually screened from various plant sources in Western Ghats, India, and subsequently molecular docking studies were performed to identify the efficacy of the drug in retroviral infection particularly coronavirus infection. Results Out of 57 phytochemicals screened initially based on the structural and physicochemical properties, 39 were effectively used for the docking analysis. Finally, 5 lead compounds with highest hydrophobic interaction and number of H-bonds were screened. Results from the interaction analysis suggest Piperolactam A to be pocketed well with good hydrophobic interaction with the residues in the binding region R1. ADME and toxicity profiling also reveals Piperolactam A with higher LogS values indicating higher permeation and hydrophilicity. Toxicity profiling suggests that the 5 screened compounds to be relatively safe. Conclusion The in silico methods used in this study suggests that the compound Piperolactam A to be the most effective inhibitor of S-protein from binding to the GRP78 receptor. By blocking the binding of the S-protein to the CS-GRP78 cell surface receptor, they can inhibit the binding of the virus to the host. Supplementary Information The online version contains supplementary material available at 10.1186/s43088-021-00095-x.
Collapse
Affiliation(s)
- Ram Kothandan
- Bioinformatics Laboratory, Department of Biotechnology, Kumaraguru College of Technology, Coimbatore, India
| | | | - Janamitra Arjun
- Bioinformatics Laboratory, Department of Biotechnology, Kumaraguru College of Technology, Coimbatore, India
| | - Rejoe Raymond Michael Raj
- Bioinformatics Laboratory, Department of Biotechnology, Kumaraguru College of Technology, Coimbatore, India
| | - Sowfia Syed
- Bioinformatics Laboratory, Department of Biotechnology, Kumaraguru College of Technology, Coimbatore, India
| |
Collapse
|
9
|
Inhibition of α-amylase by polyphenolic compounds: Substrate digestion, binding interactions and nutritional intervention. Trends Food Sci Technol 2020. [DOI: 10.1016/j.tifs.2020.08.003] [Citation(s) in RCA: 48] [Impact Index Per Article: 12.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/07/2023]
|
10
|
Pavlovicz RE, Park H, DiMaio F. Efficient consideration of coordinated water molecules improves computational protein-protein and protein-ligand docking discrimination. PLoS Comput Biol 2020; 16:e1008103. [PMID: 32956350 PMCID: PMC7529342 DOI: 10.1371/journal.pcbi.1008103] [Citation(s) in RCA: 24] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/07/2019] [Revised: 10/01/2020] [Accepted: 06/29/2020] [Indexed: 12/25/2022] Open
Abstract
Highly coordinated water molecules are frequently an integral part of protein-protein and protein-ligand interfaces. We introduce an updated energy model that efficiently captures the energetic effects of these ordered water molecules on the surfaces of proteins. A two-stage method is developed in which polar groups arranged in geometries suitable for water placement are first identified, then a modified Monte Carlo simulation allows highly coordinated waters to be placed on the surface of a protein while simultaneously sampling amino acid side chain orientations. This “semi-explicit” water model is implemented in Rosetta and is suitable for both structure prediction and protein design. We show that our new approach and energy model yield significant improvements in native structure recovery of protein-protein and protein-ligand docking discrimination tests. Well-coordinated water molecules—those forming multiple hydrogen bonds with nearby polar groups—play an important role in the structure of biomolecular systems, yet the effect of these waters is often not considered in molecular energy computations. In this paper, we describe a method to efficiently consider these water molecules both implicitly and explicitly at the interfaces formed by two polar molecules. In computations related to determining how a protein interacts with binding partners, we show that the use of this new method significantly improves results. Future application of this approach may improve the design of new protein and small molecule drugs.
Collapse
Affiliation(s)
- Ryan E. Pavlovicz
- Department of Biochemistry, University of Washington, Seattle, Washington, United States of America
- Institute for Protein Design, University of Washington, Seattle, Washington, United States of America
| | - Hahnbeom Park
- Department of Biochemistry, University of Washington, Seattle, Washington, United States of America
- Institute for Protein Design, University of Washington, Seattle, Washington, United States of America
| | - Frank DiMaio
- Department of Biochemistry, University of Washington, Seattle, Washington, United States of America
- Institute for Protein Design, University of Washington, Seattle, Washington, United States of America
- * E-mail:
| |
Collapse
|
11
|
Novel Sequence Feature of SecA Translocase Protein Unique to the Thermophilic Bacteria: Bioinformatics Analyses to Investigate Their Potential Roles. Microorganisms 2019; 8:microorganisms8010059. [PMID: 31905784 PMCID: PMC7023208 DOI: 10.3390/microorganisms8010059] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/19/2019] [Revised: 12/24/2019] [Accepted: 12/25/2019] [Indexed: 11/25/2022] Open
Abstract
SecA is an evolutionarily conserved protein that plays an indispensable role in the secretion of proteins across the bacterial cell membrane. Comparative analyses of SecA homologs have identified two large conserved signature inserts (CSIs) that are unique characteristics of thermophilic bacteria. A 50 aa conserved insert in SecA is exclusively present in the SecA homologs from the orders Thermotogales and Aquificales, while a 76 aa insert in SecA is specific for the order Thermales and Hydrogenibacillus schlegelii. Phylogenetic analyses on SecA sequences show that the shared presence of these CSIs in unrelated groups of thermophiles is not due to lateral gene transfers, but instead these large CSIs have likely originated independently in these lineages due to their advantageous function. Both of these CSIs are located in SecA protein in a surface exposed region within the ATPase domain. To gain insights into the functional significance of the 50 aa CSI in SecA, molecular dynamics (MD) simulations were performed at two different temperatures using ADP-bound SecA from Thermotoga maritima. These analyses have identified a conserved network of water molecules near the 50 aa insert in which the Glu185 residue from the CSI is found to play a key role towards stabilizing these interactions. The results provide evidence for the possible role of the 50 aa CSI in stabilizing the binding interaction of ADP/ATP, which is required for SecA function. Additionally, the surface-exposed CSIs in SecA, due to their potential to make novel protein-protein interactions, could also contribute to the thermostability of SecA from thermophilic bacteria.
Collapse
|
12
|
Orr AA, Yang J, Sule N, Chawla R, Hull KG, Zhu M, Romo D, Lele PP, Jayaraman A, Manson MD, Tamamis P. Molecular Mechanism for Attractant Signaling to DHMA by E. coli Tsr. Biophys J 2019; 118:492-504. [PMID: 31839263 DOI: 10.1016/j.bpj.2019.11.3382] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/19/2019] [Revised: 11/05/2019] [Accepted: 11/19/2019] [Indexed: 12/20/2022] Open
Abstract
The attractant chemotaxis response of Escherichia coli to norepinephrine requires that it be converted to 3,4-dihydroxymandelic acid (DHMA) by the monoamine oxidase TynA and the aromatic aldehyde dehydrogenase FeaB. DHMA is sensed by the serine chemoreceptor Tsr, and the attractant response requires that at least one subunit of the periplasmic domain of the Tsr homodimer (pTsr) has an intact serine-binding site. DHMA that is generated in vivo by E. coli is expected to be a racemic mixture of the (R) and (S) enantiomers, so it has been unclear whether one or both chiral forms are active. Here, we used a combination of state-of-the-art tools in molecular docking and simulations, including an in-house simulation-based docking protocol, to investigate the binding properties of (R)-DHMA and (S)-DHMA to E. coli pTsr. Our studies computationally predicted that (R)-DHMA should promote a stronger attractant response than (S)-DHMA because of a consistently greater-magnitude piston-like pushdown of the pTsr α-helix 4 toward the membrane upon binding of (R)-DHMA than upon binding of (S)-DHMA. This displacement is caused primarily by interaction of DHMA with Tsr residue Thr156, which has been shown by genetic studies to be critical for the attractant response to L-serine and DHMA. These findings led us to separate the two chiral species and test their effectiveness as chemoattractants. Both the tethered cell and motility migration coefficient assays validated the prediction that (R)-DHMA is a stronger attractant than (S)-DHMA. Our study demonstrates that refined computational docking and simulation studies combined with experiments can be used to investigate situations in which subtle differences between ligands may lead to diverse chemotactic responses.
Collapse
Affiliation(s)
- Asuka A Orr
- Artie McFerrin Department of Chemical Engineering, Texas A&M University, College Station, Texas
| | - Jingyun Yang
- Artie McFerrin Department of Chemical Engineering, Texas A&M University, College Station, Texas
| | - Nitesh Sule
- Artie McFerrin Department of Chemical Engineering, Texas A&M University, College Station, Texas
| | - Ravi Chawla
- Artie McFerrin Department of Chemical Engineering, Texas A&M University, College Station, Texas
| | - Kenneth G Hull
- Department of Chemistry & Biochemistry and CPRIT Synthesis and Drug-Lead Discovery Laboratory, Baylor University, Waco, Texas
| | - Mingzhao Zhu
- Department of Chemistry & Biochemistry and CPRIT Synthesis and Drug-Lead Discovery Laboratory, Baylor University, Waco, Texas
| | - Daniel Romo
- Department of Chemistry & Biochemistry and CPRIT Synthesis and Drug-Lead Discovery Laboratory, Baylor University, Waco, Texas
| | - Pushkar P Lele
- Artie McFerrin Department of Chemical Engineering, Texas A&M University, College Station, Texas
| | - Arul Jayaraman
- Artie McFerrin Department of Chemical Engineering, Texas A&M University, College Station, Texas
| | - Michael D Manson
- Department of Biology, Texas A&M University, College Station, Texas.
| | - Phanourios Tamamis
- Artie McFerrin Department of Chemical Engineering, Texas A&M University, College Station, Texas.
| |
Collapse
|
13
|
Lu J, Hou X, Wang C, Zhang Y. Incorporating Explicit Water Molecules and Ligand Conformation Stability in Machine-Learning Scoring Functions. J Chem Inf Model 2019; 59:4540-4549. [PMID: 31638801 PMCID: PMC6878146 DOI: 10.1021/acs.jcim.9b00645] [Citation(s) in RCA: 60] [Impact Index Per Article: 12.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/15/2022]
Abstract
Structure-based drug design is critically dependent on accuracy of molecular docking scoring functions, and there is of significant interest to advance scoring functions with machine learning approaches. In this work, by judiciously expanding the training set, exploring new features related to explicit mediating water molecules as well as ligand conformation stability, and applying extreme gradient boosting (XGBoost) with Δ-Vina parametrization, we have improved robustness and applicability of machine-learning scoring functions. The new scoring function ΔvinaXGB can not only perform consistently among the top compared to classical scoring functions for the CASF-2016 benchmark but also achieves significantly better prediction accuracy in different types of structures that mimic real docking applications.
Collapse
Affiliation(s)
- Jianing Lu
- Department of Chemistry, New York University, New York, New York 10003, United States
| | - Xuben Hou
- Department of Chemistry, New York University, New York, New York 10003, United States
- Department of Medicinal Chemistry, Key Laboratory of Chemical Biology (Ministry of Education), School of Pharmaceutical Science, Shandong University, Jinan, Shandong 250012, China
| | - Cheng Wang
- Department of Chemistry, New York University, New York, New York 10003, United States
| | - Yingkai Zhang
- Department of Chemistry, New York University, New York, New York 10003, United States
- NYU-ECNU Center for Computational Chemistry at NYU Shanghai, Shanghai 200062, China
| |
Collapse
|
14
|
Horváth I, Jeszenői N, Bálint M, Paragi G, Hetényi C. A Fragmenting Protocol with Explicit Hydration for Calculation of Binding Enthalpies of Target-Ligand Complexes at a Quantum Mechanical Level. Int J Mol Sci 2019; 20:ijms20184384. [PMID: 31489952 PMCID: PMC6770515 DOI: 10.3390/ijms20184384] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/05/2019] [Revised: 09/03/2019] [Accepted: 09/04/2019] [Indexed: 12/22/2022] Open
Abstract
Optimization of the enthalpy component of binding thermodynamics of drug candidates is a successful pathway of rational molecular design. However, the large size and missing hydration structure of target-ligand complexes often hinder such optimizations with quantum mechanical (QM) methods. At the same time, QM calculations are often necessitated for proper handling of electronic effects. To overcome the above problems, and help the QM design of new drugs, a protocol is introduced for atomic level determination of hydration structure and extraction of structures of target-ligand complex interfaces. The protocol is a combination of a previously published program MobyWat, an engine for assigning explicit water positions, and Fragmenter, a new tool for optimal fragmentation of protein targets. The protocol fostered a series of fast calculations of ligand binding enthalpies at the semi-empirical QM level. Ligands of diverse chemistry ranging from small aromatic compounds up to a large peptide helix of a molecular weight of 3000 targeting a leukemia protein were selected for systematic investigations. Comparison of various combinations of implicit and explicit water models demonstrated that the presence of accurately predicted explicit water molecules in the complex interface considerably improved the agreement with experimental results. A single scaling factor was derived for conversion of QM reaction heats into binding enthalpy values. The factor links molecular structure with binding thermodynamics via QM calculations. The new protocol and scaling factor will help automated optimization of binding enthalpy in future molecular design projects.
Collapse
Affiliation(s)
- István Horváth
- Chemistry Doctoral School, University of Szeged, Dugonics tér 13, 6720 Szeged, Hungary.
| | - Norbert Jeszenői
- Institute of Physiology, Medical School, University of Pécs, Szigeti út 12, 7624 Pécs, Hungary.
| | - Mónika Bálint
- Department of Pharmacology and Pharmacotherapy, Medical School, University of Pécs, Szigeti út 12, 7624 Pécs, Hungary.
| | - Gábor Paragi
- MTA-SZTE Biomimetic Systems Research Group, Dóm tér 8, 6720 Szeged, Hungary.
- Institute of Physics, University of Pécs, Ifjúság útja 6, 7624 Pécs, Hungary.
| | - Csaba Hetényi
- Department of Pharmacology and Pharmacotherapy, Medical School, University of Pécs, Szigeti út 12, 7624 Pécs, Hungary.
| |
Collapse
|
15
|
Ben-Shalom IY, Lin C, Kurtzman T, Walker RC, Gilson MK. Simulating Water Exchange to Buried Binding Sites. J Chem Theory Comput 2019; 15:2684-2691. [PMID: 30835999 DOI: 10.1021/acs.jctc.8b01284] [Citation(s) in RCA: 25] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/21/2022]
Abstract
Traditional molecular dynamics (MD) simulations of proteins, which relies on integration of Newton's equations of motion, cannot efficiently equilibrate water occupancy for buried cavities in proteins. This leads to slow convergence of thermodynamic averages for such systems. We have addressed this challenge by efficiently integrating standard Metropolis Monte Carlo (MC) translational water moves with MD in the AMBER simulation package. The translational moves allow water to easily enter or exit buried sites in a thermodynamically correct way during a simulation. To maximize efficiency, the algorithm avoids moves that only interchange waters within the bulk around the protein instead focusing on moves that can transfer water between bulk and the protein interior. In addition, a steric grid allows avoidance of moves that would lead to obvious steric clashes, and a fast grid-based energy evaluation is used to reduce the number of expensive full energy calculations. The potential energy distribution produced using MC/MD was found to be statistically indistinguishable from that of control simulations using only MD, and the algorithm effectively equilibrated water across steric barriers and into binding pockets that are not accessible with pure MD. The MC/MD method introduced here should be of increasing utility for applications spanning protein folding, the elucidation of protein mechanisms, and free energy calculations for computer-aided drug design. It is available in version 18 release of the widely disseminated AMBER simulation package.
Collapse
Affiliation(s)
- Ido Y Ben-Shalom
- Skaggs School of Pharmacy and Pharmaceutical Sciences , University of California, San Diego , La Jolla , California 92093 , United States
| | - Charles Lin
- Department of Chemistry and Biochemistry , University of California, San Diego , La Jolla , California 92093 , United States.,GlaxoSmithKline PLC, 1250 South Collegeville Road , Collegeville , Pennsylvania 19426 , United States
| | - Tom Kurtzman
- Department of Chemistry , Lehman College, The City University of New York , 250 Bedford Park Boulevard West , Bronx , New York 10468 , United States.,Ph.D. Programs in Biochemistry and Chemistry , The Graduate Center of The City University of New York , New York , New York 10016 , United States
| | - Ross C Walker
- Department of Chemistry and Biochemistry , University of California, San Diego , La Jolla , California 92093 , United States.,GlaxoSmithKline PLC, 1250 South Collegeville Road , Collegeville , Pennsylvania 19426 , United States
| | - Michael K Gilson
- Skaggs School of Pharmacy and Pharmaceutical Sciences , University of California, San Diego , La Jolla , California 92093 , United States
| |
Collapse
|
16
|
Brysbaert G, Blossey R, Lensink MF. The Inclusion of Water Molecules in Residue Interaction Networks Identifies Additional Central Residues. Front Mol Biosci 2018; 5:88. [PMID: 30364190 PMCID: PMC6193073 DOI: 10.3389/fmolb.2018.00088] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/29/2018] [Accepted: 09/18/2018] [Indexed: 12/13/2022] Open
Abstract
The relevance of water molecules for the recognition and the interaction of biomolecules is widely appreciated. In this paper we address the role that water molecules associated to protein complexes play for the functional relevance of residues by considering their residue interaction networks (RINs). These are commonly defined on the basis of the amino acid composition of the proteins themselves, disregarding the solvation state of the protein. We determine properties of the RINs of two protein complexes, colicin E2/Im2 and barnase/barstar, with and without associated water molecules, using a previously developed methodology and its associated application RINspector. We find that the inclusion of water molecules in RINs leads to an increase in the number of central residues which adds a novel mechanism to the relevance of water molecules for protein function.
Collapse
Affiliation(s)
- Guillaume Brysbaert
- CNRS UMR8576 UGSF, Institute for Structural and Functional Glycobiology, University of Lille, Lille, France
| | - Ralf Blossey
- CNRS UMR8576 UGSF, Institute for Structural and Functional Glycobiology, University of Lille, Lille, France
| | - Marc F Lensink
- CNRS UMR8576 UGSF, Institute for Structural and Functional Glycobiology, University of Lille, Lille, France
| |
Collapse
|
17
|
Aronica PGA, Fox SJ, Verma CS. Comparison of Charge Derivation Methods Applied to Amino Acid Parameterization. ACS OMEGA 2018; 3:4664-4673. [PMID: 31458687 PMCID: PMC6641686 DOI: 10.1021/acsomega.8b00438] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 03/08/2018] [Accepted: 04/13/2018] [Indexed: 06/10/2023]
Abstract
When using non-natural amino acids in computational simulations of proteins, it is necessary to ensure appropriate parameterization of the new amino acids toward the creation of appropriate input files. In particular, the charges on the atoms may have to be derived de novo and ad hoc for the new species. As there are many variables in the charge derivation process, an investigation was devised to compare different approaches and determine their effect on simulations. This was done with the purpose to identify the methods which produced results compatible with the existing parameters. It was found in this study that all analyzed charge derivation methods reproduce with sufficient accuracy the literature values and can be used with confidence when parameterizing novel species.
Collapse
Affiliation(s)
- Pietro G. A. Aronica
- Bioinformatics
Institute (A*STAR), 30
Biopolis Street, #07-01 Matrix, 138671, Singapore
| | - Stephen J. Fox
- Bioinformatics
Institute (A*STAR), 30
Biopolis Street, #07-01 Matrix, 138671, Singapore
| | - Chandra S. Verma
- Bioinformatics
Institute (A*STAR), 30
Biopolis Street, #07-01 Matrix, 138671, Singapore
- Department
of Biological Sciences, National University
of Singapore, 16 Science
Drive, 117558, Singapore
- School
of Biological Sciences, Nanyang Technological
University, 60 Nanyang
Drive, 637551, Singapore
| |
Collapse
|
18
|
Padhorny D, Hall DR, Mirzaei H, Mamonov AB, Moghadasi M, Alekseenko A, Beglov D, Kozakov D. Protein-ligand docking using FFT based sampling: D3R case study. J Comput Aided Mol Des 2018; 32:225-230. [PMID: 29101520 PMCID: PMC5767528 DOI: 10.1007/s10822-017-0069-7] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/01/2017] [Accepted: 09/16/2017] [Indexed: 12/15/2022]
Abstract
Fast Fourier transform (FFT) based approaches have been successful in application to modeling of relatively rigid protein-protein complexes. Recently, we have been able to adapt the FFT methodology to treatment of flexible protein-peptide interactions. Here, we report our latest attempt to expand the capabilities of the FFT approach to treatment of flexible protein-ligand interactions in application to the D3R PL-2016-1 challenge. Based on the D3R assessment, our FFT approach in conjunction with Monte Carlo minimization off-grid refinement was among the top performing methods in the challenge. The potential advantage of our method is its ability to globally sample the protein-ligand interaction landscape, which will be explored in further applications.
Collapse
Affiliation(s)
- Dzmitry Padhorny
- Department of Applied Mathematics and Statistics, Stony Brook University, Stony Brook, NY, 11794, USA
- Laufer Center for Physical and Quantitative Biology, Stony Brook University, Stony Brook, NY, 11794, USA
| | | | - Hanieh Mirzaei
- Department of Biomedical Engineering, Boston University, Boston, MA, 02215, USA
| | - Artem B Mamonov
- Department of Biomedical Engineering, Boston University, Boston, MA, 02215, USA
| | - Mohammad Moghadasi
- Department of Biomedical Engineering, Boston University, Boston, MA, 02215, USA
| | - Andrey Alekseenko
- Moscow Institute of Physics and Technology (State University), Institutskii per. 9, Dolgoprudny, Moscow Oblast, Russia, 141700
- Institute of Computer Aided Design of the Russian Academy of Sciences, 19/18, 2-nd Brestskaya St, Moscow, Russia, 123056
| | - Dmitri Beglov
- Department of Biomedical Engineering, Boston University, Boston, MA, 02215, USA.
| | - Dima Kozakov
- Department of Applied Mathematics and Statistics, Stony Brook University, Stony Brook, NY, 11794, USA.
- Laufer Center for Physical and Quantitative Biology, Stony Brook University, Stony Brook, NY, 11794, USA.
| |
Collapse
|
19
|
Gomes MN, Muratov EN, Pereira M, Peixoto JC, Rosseto LP, Cravo PVL, Andrade CH, Neves BJ. Chalcone Derivatives: Promising Starting Points for Drug Design. Molecules 2017; 22:E1210. [PMID: 28757583 PMCID: PMC6152227 DOI: 10.3390/molecules22081210] [Citation(s) in RCA: 199] [Impact Index Per Article: 28.4] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/16/2017] [Revised: 07/11/2017] [Accepted: 07/14/2017] [Indexed: 12/20/2022] Open
Abstract
Medicinal chemists continue to be fascinated by chalcone derivatives because of their simple chemistry, ease of hydrogen atom manipulation, straightforward synthesis, and a variety of promising biological activities. However, chalcones have still not garnered deserved attention, especially considering their high potential as chemical sources for designing and developing new effective drugs. In this review, we summarize current methodological developments towards the design and synthesis of new chalcone derivatives and state-of-the-art medicinal chemistry strategies (bioisosterism, molecular hybridization, and pro-drug design). We also highlight the applicability of computer-assisted drug design approaches to chalcones and address how this may contribute to optimizing research outputs and lead to more successful and cost-effective drug discovery endeavors. Lastly, we present successful examples of the use of chalcones and suggest possible solutions to existing limitations.
Collapse
Affiliation(s)
- Marcelo N Gomes
- Laboratory for Molecular Modeling and Drug Design, Faculty of Pharmacy, Universidade Federal de Goiás, Setor Leste Universitário, Goiânia 74605-510, Brazil.
| | - Eugene N Muratov
- Laboratory for Molecular Modeling, Eshelman School of Pharmacy, University of North Carolina, Chapel Hill, NC 27955-7568, USA.
| | - Maristela Pereira
- Laboratório de Biologia Molecular, Instituto de Ciências Biológicas, Universidade Federal de Goiás, Goiânia 74001-970, Brazil.
| | - Josana C Peixoto
- Programa de Pós-Graduação em Sociedade, Tecnologia e Meio Ambiente, Centro Universitário de Anápolis-UniEVANGÉLICA, Anápolis 75083-515, Brazil.
| | - Lucimar P Rosseto
- Programa de Pós-Graduação em Sociedade, Tecnologia e Meio Ambiente, Centro Universitário de Anápolis-UniEVANGÉLICA, Anápolis 75083-515, Brazil.
| | - Pedro V L Cravo
- Programa de Pós-Graduação em Sociedade, Tecnologia e Meio Ambiente, Centro Universitário de Anápolis-UniEVANGÉLICA, Anápolis 75083-515, Brazil.
- GHTM/Instituto de Higiene e Medicina Tropical, Universidade Nova de Lisboa, 1349-008 Lisboa, Portugal.
| | - Carolina H Andrade
- Laboratory for Molecular Modeling and Drug Design, Faculty of Pharmacy, Universidade Federal de Goiás, Setor Leste Universitário, Goiânia 74605-510, Brazil.
| | - Bruno J Neves
- Laboratory for Molecular Modeling and Drug Design, Faculty of Pharmacy, Universidade Federal de Goiás, Setor Leste Universitário, Goiânia 74605-510, Brazil.
- Laboratório de Biologia Molecular, Instituto de Ciências Biológicas, Universidade Federal de Goiás, Goiânia 74001-970, Brazil.
- Programa de Pós-Graduação em Sociedade, Tecnologia e Meio Ambiente, Centro Universitário de Anápolis-UniEVANGÉLICA, Anápolis 75083-515, Brazil.
| |
Collapse
|
20
|
Spyrakis F, Ahmed MH, Bayden AS, Cozzini P, Mozzarelli A, Kellogg GE. The Roles of Water in the Protein Matrix: A Largely Untapped Resource for Drug Discovery. J Med Chem 2017; 60:6781-6827. [PMID: 28475332 DOI: 10.1021/acs.jmedchem.7b00057] [Citation(s) in RCA: 98] [Impact Index Per Article: 14.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/01/2023]
Abstract
The value of thoroughly understanding the thermodynamics specific to a drug discovery/design study is well known. Over the past decade, the crucial roles of water molecules in protein structure, function, and dynamics have also become increasingly appreciated. This Perspective explores water in the biological environment by adopting its point of view in such phenomena. The prevailing thermodynamic models of the past, where water was seen largely in terms of an entropic gain after its displacement by a ligand, are now known to be much too simplistic. We adopt a set of terminology that describes water molecules as being "hot" and "cold", which we have defined as being easy and difficult to displace, respectively. The basis of these designations, which involve both enthalpic and entropic water contributions, are explored in several classes of biomolecules and structural motifs. The hallmarks for characterizing water molecules are examined, and computational tools for evaluating water-centric thermodynamics are reviewed. This Perspective's summary features guidelines for exploiting water molecules in drug discovery.
Collapse
Affiliation(s)
- Francesca Spyrakis
- Dipartimento di Scienza e Tecnologia del Farmaco, Università degli Studi di Torino , Via Pietro Giuria 9, 10125 Torino, Italy
| | - Mostafa H Ahmed
- Department of Medicinal Chemistry & Institute for Structural Biology, Drug Discovery and Development, Virginia Commonwealth University , Richmond, Virginia 23298-0540, United States
| | - Alexander S Bayden
- CMD Bioscience , 5 Science Park, New Haven, Connecticut 06511, United States
| | - Pietro Cozzini
- Dipartimento di Scienze degli Alimenti e del Farmaco, Laboratorio di Modellistica Molecolare, Università degli Studi di Parma , Parco Area delle Scienze 59/A, 43121 Parma, Italy
| | - Andrea Mozzarelli
- Dipartimento di Scienze degli Alimenti e del Farmaco, Laboratorio di Biochimica, Università degli Studi di Parma , Parco Area delle Scienze 23/A, 43121 Parma, Italy.,Istituto di Biofisica, Consiglio Nazionale delle Ricerche , Via Moruzzi 1, 56124 Pisa, Italy
| | - Glen E Kellogg
- Department of Medicinal Chemistry & Institute for Structural Biology, Drug Discovery and Development, Virginia Commonwealth University , Richmond, Virginia 23298-0540, United States
| |
Collapse
|
21
|
Allen WJ, Balius TE, Mukherjee S, Brozell SR, Moustakas DT, Lang PT, Case DA, Kuntz ID, Rizzo RC. DOCK 6: Impact of new features and current docking performance. J Comput Chem 2015; 36:1132-56. [PMID: 25914306 DOI: 10.1002/jcc.23905] [Citation(s) in RCA: 449] [Impact Index Per Article: 49.9] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/04/2014] [Revised: 03/01/2015] [Accepted: 03/07/2015] [Indexed: 12/11/2022]
Abstract
This manuscript presents the latest algorithmic and methodological developments to the structure-based design program DOCK 6.7 focused on an updated internal energy function, new anchor selection control, enhanced minimization options, a footprint similarity scoring function, a symmetry-corrected root-mean-square deviation algorithm, a database filter, and docking forensic tools. An important strategy during development involved use of three orthogonal metrics for assessment and validation: pose reproduction over a large database of 1043 protein-ligand complexes (SB2012 test set), cross-docking to 24 drug-target protein families, and database enrichment using large active and decoy datasets (Directory of Useful Decoys [DUD]-E test set) for five important proteins including HIV protease and IGF-1R. Relative to earlier versions, a key outcome of the work is a significant increase in pose reproduction success in going from DOCK 4.0.2 (51.4%) → 5.4 (65.2%) → 6.7 (73.3%) as a result of significant decreases in failure arising from both sampling 24.1% → 13.6% → 9.1% and scoring 24.4% → 21.1% → 17.5%. Companion cross-docking and enrichment studies with the new version highlight other strengths and remaining areas for improvement, especially for systems containing metal ions. The source code for DOCK 6.7 is available for download and free for academic users at http://dock.compbio.ucsf.edu/.
Collapse
Affiliation(s)
- William J Allen
- Department of Applied Mathematics & Statistics, Stony Brook University, Stony Brook, New York, 11794
| | | | | | | | | | | | | | | | | |
Collapse
|
22
|
Spyrakis F, Cavasotto CN. Open challenges in structure-based virtual screening: Receptor modeling, target flexibility consideration and active site water molecules description. Arch Biochem Biophys 2015; 583:105-19. [DOI: 10.1016/j.abb.2015.08.002] [Citation(s) in RCA: 63] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/30/2015] [Revised: 08/03/2015] [Accepted: 08/03/2015] [Indexed: 01/05/2023]
|
23
|
Danishuddin M, Khan AU. Structure based virtual screening to discover putative drug candidates: Necessary considerations and successful case studies. Methods 2015; 71:135-45. [DOI: 10.1016/j.ymeth.2014.10.019] [Citation(s) in RCA: 30] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/29/2014] [Revised: 09/25/2014] [Accepted: 10/17/2014] [Indexed: 12/19/2022] Open
|
24
|
Sabbadin D, Ciancetta A, Moro S. Perturbation of fluid dynamics properties of water molecules during G protein-coupled receptor-ligand recognition: the human A2A adenosine receptor as a key study. J Chem Inf Model 2014; 54:2846-55. [PMID: 25245783 DOI: 10.1021/ci500397y] [Citation(s) in RCA: 21] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/13/2023]
Abstract
Recent advances in structural biology revealed that water molecules play a crucial structural role in the protein architecture and ligand binding of G protein-coupled receptors. In this work, we present an alternative approach to monitor the time-dependent organization of water molecules during the final stage of the ligand-receptor recognition process by means of membrane molecular dynamics simulations. We inspect the variation of fluid dynamics properties of water molecules upon ligand binding with the aim to correlate the results with the binding affinities. The outcomes of this analysis are transferred into a bidimensional graph called water fluid dynamics maps, that allow a fast graphical identification of protein "hot-spots" characterized by peculiar shape and electrostatic properties that can play a critical role in ligand binding. We hopefully believe that the proposed approach might represent a valuable tool for structure-based drug discovery that can be extended to cases where crystal structures are not yet available, or have not been solved at high resolution.
Collapse
Affiliation(s)
- Davide Sabbadin
- Molecular Modeling Section (MMS), Dipartimento di Scienze del Farmaco, Università di Padova , via Marzolo 5, 35131 Padova, Italy
| | | | | |
Collapse
|
25
|
Grinter SZ, Zou X. Challenges, applications, and recent advances of protein-ligand docking in structure-based drug design. Molecules 2014; 19:10150-76. [PMID: 25019558 PMCID: PMC6270832 DOI: 10.3390/molecules190710150] [Citation(s) in RCA: 123] [Impact Index Per Article: 12.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/23/2014] [Revised: 06/13/2014] [Accepted: 07/02/2014] [Indexed: 11/16/2022] Open
Abstract
The docking methods used in structure-based virtual database screening offer the ability to quickly and cheaply estimate the affinity and binding mode of a ligand for the protein receptor of interest, such as a drug target. These methods can be used to enrich a database of compounds, so that more compounds that are subsequently experimentally tested are found to be pharmaceutically interesting. In addition, like all virtual screening methods used for drug design, structure-based virtual screening can focus on curated libraries of synthesizable compounds, helping to reduce the expense of subsequent experimental verification. In this review, we introduce the protein-ligand docking methods used for structure-based drug design and other biological applications. We discuss the fundamental challenges facing these methods and some of the current methodological topics of interest. We also discuss the main approaches for applying protein-ligand docking methods. We end with a discussion of the challenging aspects of evaluating or benchmarking the accuracy of docking methods for their improvement, and discuss future directions.
Collapse
Affiliation(s)
- Sam Z Grinter
- Informatics Institute, University of Missouri, Columbia, MO 65211, USA.
| | - Xiaoqin Zou
- Informatics Institute, University of Missouri, Columbia, MO 65211, USA.
| |
Collapse
|
26
|
Bodnarchuk MS, Viner R, Michel J, Essex JW. Strategies to calculate water binding free energies in protein-ligand complexes. J Chem Inf Model 2014; 54:1623-33. [PMID: 24684745 DOI: 10.1021/ci400674k] [Citation(s) in RCA: 46] [Impact Index Per Article: 4.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/07/2023]
Abstract
Water molecules are commonplace in protein binding pockets, where they can typically form a complex between the protein and a ligand or become displaced upon ligand binding. As a result, it is often of great interest to establish both the binding free energy and location of such molecules. Several approaches to predicting the location and affinity of water molecules to proteins have been proposed and utilized in the literature, although it is often unclear which method should be used under what circumstances. We report here a comparison between three such methodologies, Just Add Water Molecules (JAWS), Grand Canonical Monte Carlo (GCMC), and double-decoupling, in the hope of understanding the advantages and limitations of each method when applied to enclosed binding sites. As a result, we have adapted the JAWS scoring procedure, allowing the binding free energies of strongly bound water molecules to be calculated to a high degree of accuracy, requiring significantly less computational effort than more rigorous approaches. The combination of JAWS and GCMC offers a route to a rapid scheme capable of both locating and scoring water molecules for rational drug design.
Collapse
Affiliation(s)
- Michael S Bodnarchuk
- School of Chemistry, University of Southampton , Highfield, Southampton, SO17 1BJ, U.K
| | | | | | | |
Collapse
|
27
|
Lensink MF, Moal IH, Bates PA, Kastritis PL, Melquiond ASJ, Karaca E, Schmitz C, van Dijk M, Bonvin AMJJ, Eisenstein M, Jiménez-García B, Grosdidier S, Solernou A, Pérez-Cano L, Pallara C, Fernández-Recio J, Xu J, Muthu P, Praneeth Kilambi K, Gray JJ, Grudinin S, Derevyanko G, Mitchell JC, Wieting J, Kanamori E, Tsuchiya Y, Murakami Y, Sarmiento J, Standley DM, Shirota M, Kinoshita K, Nakamura H, Chavent M, Ritchie DW, Park H, Ko J, Lee H, Seok C, Shen Y, Kozakov D, Vajda S, Kundrotas PJ, Vakser IA, Pierce BG, Hwang H, Vreven T, Weng Z, Buch I, Farkash E, Wolfson HJ, Zacharias M, Qin S, Zhou HX, Huang SY, Zou X, Wojdyla JA, Kleanthous C, Wodak SJ. Blind prediction of interfacial water positions in CAPRI. Proteins 2014; 82:620-32. [PMID: 24155158 PMCID: PMC4582081 DOI: 10.1002/prot.24439] [Citation(s) in RCA: 44] [Impact Index Per Article: 4.4] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/24/2013] [Revised: 09/16/2013] [Accepted: 09/26/2013] [Indexed: 12/30/2022]
Abstract
We report the first assessment of blind predictions of water positions at protein-protein interfaces, performed as part of the critical assessment of predicted interactions (CAPRI) community-wide experiment. Groups submitting docking predictions for the complex of the DNase domain of colicin E2 and Im2 immunity protein (CAPRI Target 47), were invited to predict the positions of interfacial water molecules using the method of their choice. The predictions-20 groups submitted a total of 195 models-were assessed by measuring the recall fraction of water-mediated protein contacts. Of the 176 high- or medium-quality docking models-a very good docking performance per se-only 44% had a recall fraction above 0.3, and a mere 6% above 0.5. The actual water positions were in general predicted to an accuracy level no better than 1.5 Å, and even in good models about half of the contacts represented false positives. This notwithstanding, three hotspot interface water positions were quite well predicted, and so was one of the water positions that is believed to stabilize the loop that confers specificity in these complexes. Overall the best interface water predictions was achieved by groups that also produced high-quality docking models, indicating that accurate modelling of the protein portion is a determinant factor. The use of established molecular mechanics force fields, coupled to sampling and optimization procedures also seemed to confer an advantage. Insights gained from this analysis should help improve the prediction of protein-water interactions and their role in stabilizing protein complexes.
Collapse
Affiliation(s)
- Marc F Lensink
- Interdisciplinary Research Institute USR3078 CNRS, University Lille North of France, Villeneuve d'Ascq, France
| | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | |
Collapse
|
28
|
Michel J. Current and emerging opportunities for molecular simulations in structure-based drug design. Phys Chem Chem Phys 2014; 16:4465-77. [PMID: 24469595 PMCID: PMC4256725 DOI: 10.1039/c3cp54164a] [Citation(s) in RCA: 36] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/02/2013] [Accepted: 01/10/2014] [Indexed: 01/29/2023]
Abstract
An overview of the current capabilities and limitations of molecular simulation of biomolecular complexes in the context of computer-aided drug design is provided. Steady improvements in computer hardware coupled with more refined representations of energetics are leading to a new appreciation of the driving forces of molecular recognition. Molecular simulations are poised to more frequently guide the interpretation of biophysical measurements of biomolecular complexes. Ligand design strategies emerge from detailed analyses of computed structural ensembles. The feasibility of routine applications to ligand optimization problems hinges upon successful extensive large scale validation studies and the development of protocols to intelligently automate computations.
Collapse
Affiliation(s)
- Julien Michel
- EaStCHEM School of Chemistry, Joseph Black Building, The King's Buildings, Edinburgh, EH9 3JJ, UK.
| |
Collapse
|
29
|
Accounting for Target Flexibility and Water Molecules by Docking to Ensembles of Target Structures: The HCV NS5B Palm Site I Inhibitors Case Study. J Chem Inf Model 2013; 54:481-97. [DOI: 10.1021/ci400367m] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/05/2023]
|
30
|
Baker CM, Best RB. Matching of additive and polarizable force fields for multiscale condensed phase simulations. J Chem Theory Comput 2013; 9:2826-2837. [PMID: 23997691 PMCID: PMC3752912 DOI: 10.1021/ct400116g] [Citation(s) in RCA: 23] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/01/2023]
Abstract
Inclusion of electronic polarization effects is one of the key aspects in which the accuracy of current biomolecular force fields may be improved. The principal drawback of such approaches is the computational cost, which typically ranges from 3 - 10 times that of the equivalent additive model, and may be greater for more sophisticated treatments of polarization or other many-body effects. Here, we present a multiscale approach which may be used to enhance the sampling in simulations with polarizable models, by using the additive model as a tool to explore configuration space. We use a method based on information theory to determine the charges for an additive model that has optimal overlap with the polarizable one, and we demonstrate the feasibility of enhancing sampling via a hybrid replica exchange scheme for several model systems. An additional advantage is that, in the process, we obtain a systematic method for deriving charges for an additive model that will be the natural complement to its polarizable parent. The additive charges are found by an effective coarse-graining of the polarizable force field, rather than by ad hoc procedures.
Collapse
Affiliation(s)
- Christopher M. Baker
- University of Cambridge, Department of Chemistry, Lensfield Road, Cambridge, CB2 1EW, UK
| | - Robert B. Best
- University of Cambridge, Department of Chemistry, Lensfield Road, Cambridge, CB2 1EW, UK
- Laboratory of Chemical Physics, National Institute of Diabetes and Digestive and Kidney Diseases, National Institutes of Health, Bethesda, Maryland 20892-0520, U.S.A
| |
Collapse
|
31
|
García-Sosa AT. Hydration Properties of Ligands and Drugs in Protein Binding Sites: Tightly-Bound, Bridging Water Molecules and Their Effects and Consequences on Molecular Design Strategies. J Chem Inf Model 2013; 53:1388-405. [DOI: 10.1021/ci3005786] [Citation(s) in RCA: 54] [Impact Index Per Article: 4.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022]
|
32
|
Huang X, Yang J, Zhu Y. A solvated ligand rotamer approach and its application in computational protein design. J Mol Model 2012. [PMID: 23192355 DOI: 10.1007/s00894-012-1695-6] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/27/2022]
Abstract
The structure-based design of protein-ligand interfaces with respect to different small molecules is of great significance in the discovery of functional proteins. By statistical analysis of a set of protein-ligand complex structures, it was determined that water-mediated hydrogen bonding at the protein-ligand interface plays a crucial role in governing the binding between the protein and the ligand. Based on the novel statistic results, a solvated ligand rotamer approach was developed to explicitly describe the key water molecules at the protein-ligand interface and a water-mediated hydrogen bonding model was applied in the computational protein design context to complement the continuum solvent model. The solvated ligand rotamer approach produces only one additional solvated rotamer for each rotamer in the ligand rotamer library and does not change the number of side-chain rotamers at each protein design site. This has greatly reduced the total combinatorial number in sequence selection for protein design, and the accuracy of the model was confirmed by two tests. For the water placement test, 61% of the crystal water molecules were predicted correctly in five protein-ligand complex structures. For the sequence recapitulation test, 44.7% of the amino acid identities were recovered using the solvated ligand rotamer approach and the water-mediated hydrogen bonding model, while only 30.4% were recovered when the explicitly bound waters were removed. These results indicated that the developed solvated ligand rotamer approach is promising for functional protein design targeting novel protein-ligand interactions.
Collapse
Affiliation(s)
- Xiaoqiang Huang
- Department of Chemical Engineering, Tsinghua University, Beijing 100084, China
| | | | | |
Collapse
|
33
|
Matijssen C, Silva-Santisteban MC, Westwood IM, Siddique S, Choi V, Sheldrake P, van Montfort RL, Blagg J. Benzimidazole inhibitors of the protein kinase CHK2: clarification of the binding mode by flexible side chain docking and protein-ligand crystallography. Bioorg Med Chem 2012; 20:6630-9. [PMID: 23058106 PMCID: PMC3778940 DOI: 10.1016/j.bmc.2012.09.024] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/22/2012] [Revised: 09/07/2012] [Accepted: 09/13/2012] [Indexed: 11/23/2022]
Abstract
Two closely related binding modes have previously been proposed for the ATP-competitive benzimidazole class of checkpoint kinase 2 (CHK2) inhibitors; however, neither binding mode is entirely consistent with the reported SAR. Unconstrained rigid docking of benzimidazole ligands into representative CHK2 protein crystal structures reveals an alternative binding mode involving a water-mediated interaction with the hinge region; docking which incorporates protein side chain flexibility for selected residues in the ATP binding site resulted in a refinement of the water-mediated hinge binding mode that is consistent with observed SAR. The flexible docking results are in good agreement with the crystal structures of four exemplar benzimidazole ligands bound to CHK2 which unambiguously confirmed the binding mode of these inhibitors, including the water-mediated interaction with the hinge region, and which is significantly different from binding modes previously postulated in the literature.
Collapse
Key Words
- adp, adenosine diphosphate
- atm, ataxia telangiectasia mutated
- atp, adenosine triphosphate
- chk2, checkpoint kinase 2
- gold, genetic optimisation for ligand docking
- gst, glutathione s-transferase
- kd, kinase domain
- moe, molecular operating environment
- parp, poly adp-ribose polymerase
- pdb, protein data bank
- plif, protein ligand interaction fingerprints
- sar, structure activity relationship
- sift, structural interaction fingerprints
- kinase
- chk2
- flexible docking
- crystallography
- inhibitor
Collapse
Affiliation(s)
- Cornelis Matijssen
- Cancer Research UK Cancer Therapeutics Unit, Division of Cancer Therapeutics, Institute of Cancer Research, Haddow Laboratories, Sutton, Surrey SM2 5NG, UK
| | - M. Cris Silva-Santisteban
- Cancer Research UK Cancer Therapeutics Unit, Division of Cancer Therapeutics, Institute of Cancer Research, Haddow Laboratories, Sutton, Surrey SM2 5NG, UK
- Division of Structural Biology, Institute of Cancer Research, Chester Beatty Laboratories, Chelsea, London SW3 6JB, UK
| | - Isaac M. Westwood
- Cancer Research UK Cancer Therapeutics Unit, Division of Cancer Therapeutics, Institute of Cancer Research, Haddow Laboratories, Sutton, Surrey SM2 5NG, UK
- Division of Structural Biology, Institute of Cancer Research, Chester Beatty Laboratories, Chelsea, London SW3 6JB, UK
| | - Samerene Siddique
- Cancer Research UK Cancer Therapeutics Unit, Division of Cancer Therapeutics, Institute of Cancer Research, Haddow Laboratories, Sutton, Surrey SM2 5NG, UK
| | - Vanessa Choi
- Cancer Research UK Cancer Therapeutics Unit, Division of Cancer Therapeutics, Institute of Cancer Research, Haddow Laboratories, Sutton, Surrey SM2 5NG, UK
- Division of Structural Biology, Institute of Cancer Research, Chester Beatty Laboratories, Chelsea, London SW3 6JB, UK
| | - Peter Sheldrake
- Cancer Research UK Cancer Therapeutics Unit, Division of Cancer Therapeutics, Institute of Cancer Research, Haddow Laboratories, Sutton, Surrey SM2 5NG, UK
| | - Rob L.M. van Montfort
- Cancer Research UK Cancer Therapeutics Unit, Division of Cancer Therapeutics, Institute of Cancer Research, Haddow Laboratories, Sutton, Surrey SM2 5NG, UK
- Division of Structural Biology, Institute of Cancer Research, Chester Beatty Laboratories, Chelsea, London SW3 6JB, UK
| | - Julian Blagg
- Cancer Research UK Cancer Therapeutics Unit, Division of Cancer Therapeutics, Institute of Cancer Research, Haddow Laboratories, Sutton, Surrey SM2 5NG, UK
| |
Collapse
|
34
|
Designing electrostatic interactions in biological systems via charge optimization or combinatorial approaches: insights and challenges with a continuum electrostatic framework. Theor Chem Acc 2012. [DOI: 10.1007/s00214-012-1252-5] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
|
35
|
Taboureau O, Baell JB, Fernández-Recio J, Villoutreix BO. Established and emerging trends in computational drug discovery in the structural genomics era. ACTA ACUST UNITED AC 2012; 19:29-41. [PMID: 22284352 DOI: 10.1016/j.chembiol.2011.12.007] [Citation(s) in RCA: 52] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/25/2011] [Revised: 12/05/2011] [Accepted: 12/08/2011] [Indexed: 12/01/2022]
Abstract
Bioinformatics and chemoinformatics approaches contribute to hit discovery, hit-to-lead optimization, safety profiling, and target identification and enhance our overall understanding of the health and disease states. A vast repertoire of computational methods has been reported and increasingly combined in order to address more and more challenging targets or complex molecular mechanisms in the context of large-scale integration of structure and bioactivity data produced by private and public drug research. This review explores some key computational methods directly linked to drug discovery and chemical biology with a special emphasis on compound collection preparation, virtual screening, protein docking, and systems pharmacology. A list of generally freely available software packages and online resources is provided, and examples of successful applications are briefly commented upon.
Collapse
Affiliation(s)
- Olivier Taboureau
- Center for Biological Sequences Analysis, Department of Systems Biology, Technical University of Denmark, 2800 Kongens Lyngby, Denmark
| | | | | | | |
Collapse
|
36
|
Huggins DJ, Altman MD, Tidor B. Evaluation of an inverse molecular design algorithm in a model binding site. Proteins 2009; 75:168-86. [PMID: 18831031 DOI: 10.1002/prot.22226] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/07/2022]
Abstract
Computational molecular design is a useful tool in modern drug discovery. Virtual screening is an approach that docks and then scores individual members of compound libraries. In contrast to this forward approach, inverse approaches construct compounds from fragments, such that the computed affinity, or a combination of relevant properties, is optimized. We have recently developed a new inverse approach to drug design based on the dead-end elimination and A* algorithms employing a physical potential function. This approach has been applied to combinatorially constructed libraries of small-molecule ligands to design high-affinity HIV-1 protease inhibitors (Altman et al., J Am Chem Soc 2008;130:6099-6013). Here we have evaluated the new method using the well-studied W191G mutant of cytochrome c peroxidase. This mutant possesses a charged binding pocket and has been used to evaluate other design approaches. The results show that overall the new inverse approach does an excellent job of separating binders from nonbinders. For a few individual cases, scoring inaccuracies led to false positives. The majority of these involve erroneous solvation energy estimation for charged amines, anilinium ions, and phenols, which has been observed previously for a variety of scoring algorithms. Interestingly, although inverse approaches are generally expected to identify some but not all binders in a library, due to limited conformational searching, these results show excellent coverage of the known binders while still showing strong discrimination of the nonbinders.
Collapse
Affiliation(s)
- David J Huggins
- Computer Science and Artificial Intelligence Laboratory, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, USA
| | | | | |
Collapse
|