1
|
Mirgaux M, Leherte L, Wouters J. Human indoleamine-2,3-dioxygenase 2 cofactor lability and low substrate affinity explained by homology modeling, molecular dynamics and molecular docking. J Biomol Struct Dyn 2024; 42:4475-4488. [PMID: 37301605 DOI: 10.1080/07391102.2023.2220830] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/07/2023] [Accepted: 05/28/2023] [Indexed: 06/12/2023]
Abstract
The human indoleamine-2,3-dioxygenase 2 (hIDO2) protein is growing of interest as it is increasingly implicated in multiple diseases (cancer, autoimmune diseases, COVID-19). However, it is only poorly reported in the literature. Its mode of action remains unknown because it does not seem to catalyze the reaction for which it is attributed: the degradation of the L-Tryptophan into N-formyl-kynurenine. This contrasts with its paralog, the human indoleamine-2,3-dioxygenase 1 (hIDO1), which has been extensively studied in the literature and for which several inhibitors are already in clinical trials. Yet, the recent failure of one of the most advanced hIDO1 inhibitors, the Epacadostat, could be caused by a still unknown interaction between hIDO1 and hIDO2. In order to better understand the mechanism of hIDO2, and in the absence of experimental structural data, a computational study mixing homology modeling, Molecular Dynamics, and molecular docking was conducted. The present article highlights an exacerbated lability of the cofactor as well as an inadequate positioning of the substrate in the active site of hIDO2, which might bring part of an answer to its lack of activity.Communicated by Ramaswamy H. Sarma.
Collapse
Affiliation(s)
- Manon Mirgaux
- Department of Chemistry, Laboratoire de Chimie Biologique Structurale (CBS), Namur Institute of Structured Matter (NISM), Namur Research Institute for Life Science (NARILIS), University of Namur (UNamur), Namur, Belgium
| | - Laurence Leherte
- Department of Chemistry, Laboratoire de Chimie Biologique Structurale (CBS), Namur Institute of Structured Matter (NISM), Namur Research Institute for Life Science (NARILIS), University of Namur (UNamur), Namur, Belgium
| | - Johan Wouters
- Department of Chemistry, Laboratoire de Chimie Biologique Structurale (CBS), Namur Institute of Structured Matter (NISM), Namur Research Institute for Life Science (NARILIS), University of Namur (UNamur), Namur, Belgium
| |
Collapse
|
2
|
Wesołowski P, Wales DJ, Pracht P. Multilevel Framework for Analysis of Protein Folding Involving Disulfide Bond Formation. J Phys Chem B 2024; 128:3145-3156. [PMID: 38512062 PMCID: PMC11000224 DOI: 10.1021/acs.jpcb.4c00104] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/05/2024] [Revised: 03/06/2024] [Accepted: 03/06/2024] [Indexed: 03/22/2024]
Abstract
In this study, a three-layered multicenter ONIOM approach is implemented to characterize the naive folding pathway of bovine pancreatic trypsin inhibitor (BPTI). Each layer represents a distinct level of theory, where the initial layer, encompassing the entire protein, is modeled by a general all-atom force-field GFN-FF. An intermediate electronic structure layer consisting of three multicenter fragments is introduced with the state-of-the-art semiempirical tight-binding method GFN2-xTB. Higher accuracy, specifically addressing the breaking and formation of the three disulfide bonds, is achieved at the innermost layer using the composite DFT method r2SCAN-3c. Our analysis sheds light on the structural stability of BPTI, particularly the significance of interlinking disulfide bonds. The accuracy and efficiency of the multicenter QM/SQM/MM approach are benchmarked using the oxidative formation of cystine. For the folding pathway of BPTI, relative stabilities are investigated through the calculation of free energy contributions for selected intermediates, focusing on the impact of the disulfide bond. Our results highlight the intricate trade-off between accuracy and computational cost, demonstrating that the multicenter ONIOM approach provides a well-balanced and comprehensive solution to describe electronic structure effects in biomolecular systems. We conclude that multiscale energy landscape exploration provides a robust methodology for the study of intriguing biological targets.
Collapse
Affiliation(s)
- Patryk
A. Wesołowski
- Yusuf Hamied Department of Chemistry, University of Cambridge, Lensfield Road, Cambridge CB2 1EW, U.K.
| | - David J. Wales
- Yusuf Hamied Department of Chemistry, University of Cambridge, Lensfield Road, Cambridge CB2 1EW, U.K.
| | - Philipp Pracht
- Yusuf Hamied Department of Chemistry, University of Cambridge, Lensfield Road, Cambridge CB2 1EW, U.K.
| |
Collapse
|
3
|
Csóka J, Hégely B, Nagy PR, Kállay M. Development of analytic gradients for the Huzinaga quantum embedding method and its applications to large-scale hybrid and double hybrid DFT forces. J Chem Phys 2024; 160:124113. [PMID: 38530010 DOI: 10.1063/5.0194463] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/27/2023] [Accepted: 03/06/2024] [Indexed: 03/27/2024] Open
Abstract
The theory of analytic gradients is presented for the projector-based density functional theory (DFT) embedding approach utilizing the Huzinaga-equation. The advantages of the Huzinaga-equation-based formulation are demonstrated. In particular, it is shown that the projector employed does not appear in the Lagrangian, and the potential risk of numerical problems is avoided at the evaluation of the gradients. The efficient implementation of the analytic gradient theory is presented for approaches where hybrid DFT, second-order Møller-Plesset perturbation theory, or double hybrid DFT are embedded in lower-level DFT environments. To demonstrate the applicability of the method and to gain insight into its accuracy, it is applied to equilibrium geometry optimizations, transition state searches, and potential energy surface scans. Our results show that bond lengths and angles converge rapidly with the size of the embedded system. While providing structural parameters close to high-level quality for the embedded atoms, the embedding approach has the potential to relax the coordinates of the environment as well. Our demonstrations on a 171-atom zeolite and a 570-atom protein system show that the Huzinaga-equation-based embedding can accelerate (double) hybrid gradient computations by an order of magnitude with sufficient active regions and enables affordable force evaluations or geometry optimizations for molecules of hundreds of atoms.
Collapse
Affiliation(s)
- József Csóka
- Department of Physical Chemistry and Materials Science, Faculty of Chemical Technology and Biotechnology, Budapest University of Technology and Economics, Műegyetem rkp. 3., H-1111 Budapest, Hungary
- HUN-REN-BME Quantum Chemistry Research Group, Műegyetem rkp. 3., H-1111 Budapest, Hungary
- MTA-BME Lendület Quantum Chemistry Research Group, Műegyetem rkp. 3., H-1111 Budapest, Hungary
| | - Bence Hégely
- Department of Physical Chemistry and Materials Science, Faculty of Chemical Technology and Biotechnology, Budapest University of Technology and Economics, Műegyetem rkp. 3., H-1111 Budapest, Hungary
- HUN-REN-BME Quantum Chemistry Research Group, Műegyetem rkp. 3., H-1111 Budapest, Hungary
- MTA-BME Lendület Quantum Chemistry Research Group, Műegyetem rkp. 3., H-1111 Budapest, Hungary
| | - Péter R Nagy
- Department of Physical Chemistry and Materials Science, Faculty of Chemical Technology and Biotechnology, Budapest University of Technology and Economics, Műegyetem rkp. 3., H-1111 Budapest, Hungary
- HUN-REN-BME Quantum Chemistry Research Group, Műegyetem rkp. 3., H-1111 Budapest, Hungary
- MTA-BME Lendület Quantum Chemistry Research Group, Műegyetem rkp. 3., H-1111 Budapest, Hungary
| | - Mihály Kállay
- Department of Physical Chemistry and Materials Science, Faculty of Chemical Technology and Biotechnology, Budapest University of Technology and Economics, Műegyetem rkp. 3., H-1111 Budapest, Hungary
- HUN-REN-BME Quantum Chemistry Research Group, Műegyetem rkp. 3., H-1111 Budapest, Hungary
- MTA-BME Lendület Quantum Chemistry Research Group, Műegyetem rkp. 3., H-1111 Budapest, Hungary
| |
Collapse
|
4
|
Maghsoud Y, Vázquez-Montelongo EA, Yang X, Liu C, Jing Z, Lee J, Harger M, Smith AK, Espinoza M, Guo HF, Kurie JM, Dalby KN, Ren P, Cisneros GA. Computational Investigation of a Series of Small Molecules as Potential Compounds for Lysyl Hydroxylase-2 (LH2) Inhibition. J Chem Inf Model 2023; 63:986-1001. [PMID: 36779232 PMCID: PMC10233724 DOI: 10.1021/acs.jcim.2c01448] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/01/2023]
Abstract
The catalytic function of lysyl hydroxylase-2 (LH2), a member of the Fe(II)/αKG-dependent oxygenase superfamily, is to catalyze the hydroxylation of lysine to hydroxylysine in collagen, resulting in stable hydroxylysine aldehyde-derived collagen cross-links (HLCCs). Reports show that high amounts of LH2 lead to the accumulation of HLCCs, causing fibrosis and specific types of cancer metastasis. Some members of the Fe(II)/αKG-dependent family have also been reported to have intramolecular O2 tunnels, which aid in transporting one of the required cosubstrates into the active site. While LH2 can be a promising target to combat these diseases, efficacious inhibitors are still lacking. We have used computational simulations to investigate a series of 44 small molecules as lead compounds for LH2 inhibition. Tunneling analyses indicate the existence of several intramolecular tunnels. The lengths of the calculated O2-transporting tunnels in holoenzymes are relatively longer than those in the apoenzyme, suggesting that the ligands may affect the enzyme's structure and possibly block (at least partially) the tunnels. The sequence alignment analysis between LH enzymes from different organisms shows that all of the amino acid residues with the highest occurrence rate in the oxygen tunnels are conserved. Our results suggest that the enolate form of diketone compounds establishes stronger interactions with the Fe(II) in the active site. Branching the enolate compounds with functional groups such as phenyl and pyridinyl enhances the interaction with various residues around the active site. Our results provide information about possible leads for further LH2 inhibition design and development.
Collapse
Affiliation(s)
- Yazdan Maghsoud
- Department of Chemistry and Biochemistry, The University of Texas at Dallas, Richardson, Texas 75080, United States
| | - Erik Antonio Vázquez-Montelongo
- Department of Physical Medicine and Rehabilitation, University of Texas Southwestern Medical Center, Dallas, Texas 75390, United States
| | - Xudong Yang
- Department of Biomedical Engineering, The University of Texas at Austin, Austin, Texas 78712, United States
| | - Chengwen Liu
- Department of Biomedical Engineering, The University of Texas at Austin, Austin, Texas 78712, United States
| | - Zhifeng Jing
- Department of Biomedical Engineering, The University of Texas at Austin, Austin, Texas 78712, United States
| | - Juhoon Lee
- Division of Chemical Biology and Medicinal Chemistry, College of Pharmacy, The University of Texas at Austin, Austin, Texas 78712, United States
| | - Matthew Harger
- Department of Biomedical Engineering, The University of Texas at Austin, Austin, Texas 78712, United States
| | - Ally K Smith
- Department of Chemistry, University of North Texas, Denton, Texas 76201, United States
| | - Miguel Espinoza
- Department of Chemistry, University of North Texas, Denton, Texas 76201, United States
| | - Hou-Fu Guo
- Department of Molecular and Cellular Biochemistry, College of Medicine, The University of Kentucky, Lexington, Kentucky 40536, United States
| | - Jonathan M Kurie
- Department of Thoracic/Head and Neck Medical Oncology, The University of Texas MD Anderson Cancer Center, Houston, Texas 77005, United States
| | - Kevin N Dalby
- Division of Chemical Biology and Medicinal Chemistry, College of Pharmacy, The University of Texas at Austin, Austin, Texas 78712, United States
| | - Pengyu Ren
- Department of Biomedical Engineering, The University of Texas at Austin, Austin, Texas 78712, United States
| | - G Andrés Cisneros
- Department of Chemistry and Biochemistry, The University of Texas at Dallas, Richardson, Texas 75080, United States
- Department of Physics, The University of Texas at Dallas, Richardson, Texas 75080, United States
| |
Collapse
|
5
|
Bím D, Navrátil M, Gutten O, Konvalinka J, Kutil Z, Culka M, Navrátil V, Alexandrova AN, Bařinka C, Rulíšek L. Predicting Effects of Site-Directed Mutagenesis on Enzyme Kinetics by QM/MM and QM Calculations: A Case of Glutamate Carboxypeptidase II. J Phys Chem B 2022; 126:132-143. [PMID: 34978450 DOI: 10.1021/acs.jpcb.1c09240] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/06/2023]
Abstract
Quantum and molecular mechanics (QM/MM) and QM-only (cluster model) modeling techniques represent the two workhorses in mechanistic understanding of enzyme catalysis. One of the stringent tests for QM/MM and/or QM approaches is to provide quantitative answers to real-world biochemical questions, such as the effect of single-point mutations on enzyme kinetics. This translates into predicting the relative activation energies to 1-2 kcal·mol-1 accuracy; such predictions can be used for the rational design of novel enzyme variants with desired/improved characteristics. Herein, we employ glutamate carboxypeptidase II (GCPII), a dizinc metallopeptidase, also known as the prostate specific membrane antigen, as a model system. The structure and activity of this major cancer antigen have been thoroughly studied, both experimentally and computationally, which makes it an ideal model system for method development. Its reaction mechanism is quite well understood: the reaction coordinate comprises a "tetrahedral intermediate" and two transition states and experimental activation Gibbs free energy of ∼17.5 kcal·mol-1 can be inferred for the known kcat ≈ 1 s-1. We correlate experimental kinetic data (including the E424H variant, newly characterized in this work) for various GCPII mutants (kcat = 8.6 × 10-5 s-1 to 2.7 s-1) with the energy profiles calculated by QM/MM and QM-only (cluster model) approaches. We show that the near-quantitative agreement between the experimental values and the calculated activation energies (ΔH⧧) can be obtained and recommend the combination of the two protocols: QM/MM optimized structures and cluster model (QM) energetics. The trend in relative activation energies is mostly independent of the QM method (DFT functional) used. Last but not least, a satisfactory correlation between experimental and theoretical data allows us to provide qualitative and fairly simple explanations of the observed kinetic effects which are thus based on a rigorous footing.
Collapse
Affiliation(s)
- Daniel Bím
- Institute of Organic Chemistry and Biochemistry of the Czech Academy of Sciences, Flemingovo náměstí 2, 166 10 Praha 6, Czech Republic.,Department of Chemistry and Biochemistry, University of California, Los Angeles, California 90095-1569, United States
| | - Michal Navrátil
- Institute of Organic Chemistry and Biochemistry of the Czech Academy of Sciences, Flemingovo náměstí 2, 166 10 Praha 6, Czech Republic
| | - Ondrej Gutten
- Institute of Organic Chemistry and Biochemistry of the Czech Academy of Sciences, Flemingovo náměstí 2, 166 10 Praha 6, Czech Republic
| | - Jan Konvalinka
- Institute of Organic Chemistry and Biochemistry of the Czech Academy of Sciences, Flemingovo náměstí 2, 166 10 Praha 6, Czech Republic.,Department of Biochemistry, Faculty of Science, Charles University, Hlavova 2030, 2120 00 Prague, Czech Republic
| | - Zsófia Kutil
- Institute of Biotechnology of the Czech Academy of Sciences, Průmyslová 595, 252 50 Vestec, Czech Republic
| | - Martin Culka
- Institute of Organic Chemistry and Biochemistry of the Czech Academy of Sciences, Flemingovo náměstí 2, 166 10 Praha 6, Czech Republic
| | - Václav Navrátil
- Institute of Organic Chemistry and Biochemistry of the Czech Academy of Sciences, Flemingovo náměstí 2, 166 10 Praha 6, Czech Republic
| | - Anastassia N Alexandrova
- Department of Chemistry and Biochemistry, University of California, Los Angeles, California 90095-1569, United States
| | - Cyril Bařinka
- Institute of Biotechnology of the Czech Academy of Sciences, Průmyslová 595, 252 50 Vestec, Czech Republic
| | - Lubomír Rulíšek
- Institute of Organic Chemistry and Biochemistry of the Czech Academy of Sciences, Flemingovo náměstí 2, 166 10 Praha 6, Czech Republic
| |
Collapse
|
6
|
Xue Y, Wang JN, Hu W, Zheng J, Li Y, Pan X, Mo Y, Shao Y, Wang L, Mei Y. Affordable Ab Initio Path Integral for Thermodynamic Properties via Molecular Dynamics Simulations Using Semiempirical Reference Potential. J Phys Chem A 2021; 125:10677-10685. [PMID: 34894680 PMCID: PMC9108008 DOI: 10.1021/acs.jpca.1c07727] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/25/2022]
Abstract
Path integral molecular dynamics (PIMD) is becoming a routinely applied method for incorporating the nuclear quantum effect in computer simulations. However, direct PIMD simulations at an ab initio level of theory are formidably expensive. Using the protonated 1,8-bis(dimethylamino)naphthalene molecule as an example, we show in this work that the computational expense for the intramolecular proton transfer between the two nitrogen atoms can be remarkably reduced by implementing the idea of reference-potential methods. The simulation time can be easily extended to a scale of nanoseconds while maintaining the accuracy on an ab initio level of theory for thermodynamic properties. In addition, postprocessing can be carried out in parallel on massive computer nodes. A 545-fold reduction in the total CPU time can be achieved in this way as compared to a direct PIMD simulation at the same ab initio level of theory.
Collapse
Affiliation(s)
- Yuanfei Xue
- State Key Laboratory of Precision Spectroscopy, School of Physics and Electronic Science, East China Normal University, Shanghai 200062, China
| | - Jia-Ning Wang
- State Key Laboratory of Precision Spectroscopy, School of Physics and Electronic Science, East China Normal University, Shanghai 200062, China
| | - Wenxin Hu
- The Computer Center, School of Data Science & Engineering, East China Normal University, Shanghai 200062, China
| | - Jun Zheng
- The Computer Center, School of Data Science & Engineering, East China Normal University, Shanghai 200062, China
| | - Yongle Li
- Department of Physics, International Center of Quantum and Molecular Structure, and Shanghai Key Laboratory of High Temperature Superconductors, Shanghai University, Shanghai 200444, China
| | - Xiaoliang Pan
- Department of Chemistry and Biochemistry, University of Oklahoma, Norman, Oklahoma 73019, United States
| | - Yan Mo
- State Key Laboratory of Precision Spectroscopy, School of Physics and Electronic Science, East China Normal University, Shanghai 200062, China,NYU–ECNU Center for Computational Chemistry at NYU Shanghai, Shanghai 200062, China,Collaborative Innovation Center of Extreme Optics, Shanxi University, Taiyuan, Shanxi 030006, China
| | - Yihan Shao
- Department of Chemistry and Biochemistry, University of Oklahoma, Norman, Oklahoma 73019, United States
| | - Lu Wang
- Department of Chemistry and Chemical Biology, Institute for Quantitative Biomedicine, Rutgers University, Piscataway, New Jersey 08854, United States
| | - Ye Mei
- State Key Laboratory of Precision Spectroscopy, School of Physics and Electronic Science, East China Normal University, Shanghai 200062, China,NYU–ECNU Center for Computational Chemistry at NYU Shanghai, Shanghai 200062, China,Collaborative Innovation Center of Extreme Optics, Shanxi University, Taiyuan, Shanxi 030006, China
| |
Collapse
|
7
|
Kulkarni PU, Shah H, Vyas VK. Hybrid Quantum Mechanics/Molecular Mechanics (QM/MM) Simulation: A Tool for Structure-based Drug Design and Discovery. Mini Rev Med Chem 2021; 22:1096-1107. [PMID: 34620049 DOI: 10.2174/1389557521666211007115250] [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: 12/29/2020] [Revised: 04/22/2021] [Accepted: 08/10/2021] [Indexed: 11/22/2022]
Abstract
Quantum mechanics (QM) is physics based theory which explains the physical properties of nature at the level of atoms and sub-atoms. Molecular mechanics (MM) construct molecular systems through the use of classical mechanics. So, hybrid quantum mechanics and molecular mechanics (QM/MM) when combined together can act as computer-based methods which can be used to calculate structure and property data of molecular structures. Hybrid QM/MM combines the strengths of QM with accuracy and MM with speed. QM/MM simulation can also be applied for the study of chemical process in solutions as well as in the proteins, and has a great scope in structure-based drug design (CADD) and discovery. Hybrid QM/MM also applied to HTS, to derive QSAR models and due to availability of many protein crystal structures; it has a great role in computational chemistry, especially in structure- and fragment-based drug design. Fused QM/MM simulations have been developed as a widespread method to explore chemical reactions in condensed phases. In QM/MM simulations, the quantum chemistry theory is used to treat the space in which the chemical reactions occur; however the rest is defined through molecular mechanics force field (MMFF). In this review, we have extensively reviewed recent literature pertaining to the use and applications of hybrid QM/MM simulations for ligand and structure-based computational methods for the design and discovery of therapeutic agents.
Collapse
Affiliation(s)
- Prajakta U Kulkarni
- School of Pharmacy, ITM (SLS) Baroda University, Vadodara 391510, Gujarat. India
| | - Harshil Shah
- Department of Pharmaceutical Chemistry, Sardar Patel College of Pharmacy, Bakrol, Anand 388315, Gujarat. India
| | - Vivek K Vyas
- Department of Pharmaceutical Chemistry, Institute of Pharmacy, Nirma University, Ahmedabad 382481, Gujarat. India
| |
Collapse
|
8
|
Nochebuena J, Naseem-Khan S, Cisneros GA. Development and application of quantum mechanics/molecular mechanics methods with advanced polarizable potentials. WILEY INTERDISCIPLINARY REVIEWS. COMPUTATIONAL MOLECULAR SCIENCE 2021; 11:e1515. [PMID: 34367343 PMCID: PMC8341087 DOI: 10.1002/wcms.1515] [Citation(s) in RCA: 12] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/28/2020] [Accepted: 12/19/2020] [Indexed: 01/02/2023]
Abstract
Quantum mechanics/molecular mechanics (QM/MM) simulations are a popular approach to study various features of large systems. A common application of QM/MM calculations is in the investigation of reaction mechanisms in condensed-phase and biological systems. The combination of QM and MM methods to represent a system gives rise to several challenges that need to be addressed. The increase in computational speed has allowed the expanded use of more complicated and accurate methods for both QM and MM simulations. Here, we review some approaches that address several common challenges encountered in QM/MM simulations with advanced polarizable potentials, from methods to account for boundary across covalent bonds and long-range effects, to polarization and advanced embedding potentials.
Collapse
Affiliation(s)
- Jorge Nochebuena
- Department of Chemistry, University of North Texas, Denton, Texas, USA
| | - Sehr Naseem-Khan
- Department of Chemistry, University of North Texas, Denton, Texas, USA
| | - G Andrés Cisneros
- Department of Chemistry, University of North Texas, Denton, Texas, USA
| |
Collapse
|
9
|
Elmi A, Cockroft SL. Quantifying Interactions and Solvent Effects Using Molecular Balances and Model Complexes. Acc Chem Res 2021; 54:92-103. [PMID: 33315374 DOI: 10.1021/acs.accounts.0c00545] [Citation(s) in RCA: 17] [Impact Index Per Article: 5.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022]
Abstract
Where the basic units of molecular chemistry are the bonds within molecules, supramolecular chemistry is based on the interactions that occur between molecules. Understanding the "how" and "why" of the processes that govern molecular self-assembly remains an open challenge to the supramolecular community. While many interactions are readily examined in silico through electronic structure calculations, such insights may not be directly applicable to experimentalists. The practical limitations of computationally accounting for solvation is perhaps the largest bottleneck in this regard, with implicit solvation models failing to comprehensively account for the specific nature of solvent effects and explicit models incurring a prohibitively high computational cost. Since molecular recognition processes usually occur in solution, insight into the nature and effect of solvation is imperative not only for understanding these phenomena but also for the rational design of systems that exploit them.Molecular balances and supramolecular complexes have emerged as useful tools for the experimental dissection of the physicochemical basis of various noncovalent interactions, but they have historically been underexploited as a platform for the evaluation of solvent effects. Contrasting with large biological complexes, smaller synthetic model systems enable combined experimental and computational analyses, often facilitating theoretical analyses that can work in concert with experiment.Our research has focused on the development of supramolecular systems to evaluate the role of solvents in molecular recognition, and further characterize the underlying mechanisms by which molecules associate. In particular, the use of molecular balances has provided a framework to measure the magnitude of solvent effects and to examine the accuracy of solvent models. Such approaches have revealed how solvation can modulate the electronic landscape of a molecule and how competitive solvation and solvent cohesion can provide thermodynamic driving forces for association. Moreover, the use of simple model systems facilitates the interrogation and further dissection of the physicochemical origins of molecular recognition. This tandem experimental/computational approach has married less common computational techniques, like symmetry adapted perturbation theory (SAPT) and natural bonding orbital (NBO) analysis, with experimental observations to elucidate the influence of effects that are difficult to resolve experimentally (e.g., London dispersion and electron delocalization).
Collapse
Affiliation(s)
- Alex Elmi
- EaStCHEM School of Chemistry, University of Edinburgh, David Brewster Road, Edinburgh EH9 3FJ, United Kingdom
| | - Scott L. Cockroft
- EaStCHEM School of Chemistry, University of Edinburgh, David Brewster Road, Edinburgh EH9 3FJ, United Kingdom
| |
Collapse
|
10
|
Lee TS, Allen BK, Giese TJ, Guo Z, Li P, Lin C, McGee TD, Pearlman DA, Radak BK, Tao Y, Tsai HC, Xu H, Sherman W, York DM. Alchemical Binding Free Energy Calculations in AMBER20: Advances and Best Practices for Drug Discovery. J Chem Inf Model 2020; 60:5595-5623. [PMID: 32936637 PMCID: PMC7686026 DOI: 10.1021/acs.jcim.0c00613] [Citation(s) in RCA: 163] [Impact Index Per Article: 40.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/03/2023]
Abstract
Predicting protein-ligand binding affinities and the associated thermodynamics of biomolecular recognition is a primary objective of structure-based drug design. Alchemical free energy simulations offer a highly accurate and computationally efficient route to achieving this goal. While the AMBER molecular dynamics package has successfully been used for alchemical free energy simulations in academic research groups for decades, widespread impact in industrial drug discovery settings has been minimal because of the previous limitations within the AMBER alchemical code, coupled with challenges in system setup and postprocessing workflows. Through a close academia-industry collaboration we have addressed many of the previous limitations with an aim to improve accuracy, efficiency, and robustness of alchemical binding free energy simulations in industrial drug discovery applications. Here, we highlight some of the recent advances in AMBER20 with a focus on alchemical binding free energy (BFE) calculations, which are less computationally intensive than alternative binding free energy methods where full binding/unbinding paths are explored. In addition to scientific and technical advances in AMBER20, we also describe the essential practical aspects associated with running relative alchemical BFE calculations, along with recommendations for best practices, highlighting the importance not only of the alchemical simulation code but also the auxiliary functionalities and expertise required to obtain accurate and reliable results. This work is intended to provide a contemporary overview of the scientific, technical, and practical issues associated with running relative BFE simulations in AMBER20, with a focus on real-world drug discovery applications.
Collapse
Affiliation(s)
- Tai-Sung Lee
- Rutgers, the State University of New Jersey, Laboratory for Biomolecular Simulation Research, and Department of Chemistry and Chemical Biology, United States
| | - Bryce K. Allen
- Silicon Therapeutics, Boston, Massachusetts 02210, United States
| | - Timothy J. Giese
- Rutgers, the State University of New Jersey, Laboratory for Biomolecular Simulation Research, and Department of Chemistry and Chemical Biology, United States
| | - Zhenyu Guo
- Silicon Therapeutics, Boston, Massachusetts 02210, United States
| | - Pengfei Li
- Silicon Therapeutics, Boston, Massachusetts 02210, United States
| | - Charles Lin
- Silicon Therapeutics, Boston, Massachusetts 02210, United States
| | - T. Dwight McGee
- Silicon Therapeutics, Boston, Massachusetts 02210, United States
| | - David A. Pearlman
- QSimulate Incorporated, Cambridge, Massachusetts 02139, United States
| | - Brian K. Radak
- Silicon Therapeutics, Boston, Massachusetts 02210, United States
| | - Yujun Tao
- Rutgers, the State University of New Jersey, Laboratory for Biomolecular Simulation Research, and Department of Chemistry and Chemical Biology, United States
| | - Hsu-Chun Tsai
- Rutgers, the State University of New Jersey, Laboratory for Biomolecular Simulation Research, and Department of Chemistry and Chemical Biology, United States
| | - Huafeng Xu
- Silicon Therapeutics, Boston, Massachusetts 02210, United States
| | - Woody Sherman
- Silicon Therapeutics, Boston, Massachusetts 02210, United States
| | - Darrin M. York
- Rutgers, the State University of New Jersey, Laboratory for Biomolecular Simulation Research, and Department of Chemistry and Chemical Biology, United States
| |
Collapse
|
11
|
Dong C, Montes M, Al-Sawai WM. Xanthine oxidoreductase inhibition – A review of computational aspect. JOURNAL OF THEORETICAL & COMPUTATIONAL CHEMISTRY 2020. [DOI: 10.1142/s0219633620400088] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/03/2023]
Abstract
Xanthine Oxidoreductase (XOR) exists in a variety of organisms from bacteria to humans and catalyzes the oxidation of hypoxanthine to xanthine and from xanthine to uric acid. Excessive uric acid could lead to gout and hyperuricemia. In this paper, we have reviewed the recent computational studies on xanthine oxidase inhibition. Computational methods, such as molecular dynamics (molecular mechanics), quantum mechanics, and quantum mechanics/molecular mechanics (QM/MM), have been employed to investigate the binding affinity of xanthine oxidase with synthesized and isolated nature inhibitors. The limitations of different computational methods for xanthine oxidase inhibition studies were also discussed. Implications of the computational approach could be used to help to understand the existing arguments on substrate/product orientation in xanthine oxidase inhibition, which allows designing new inhibitors with higher efficacy.
Collapse
Affiliation(s)
- Chao Dong
- Department of Chemistry, The University of Texas of the Permian Basin, Odessa, Texas 79762, USA
| | - Milka Montes
- Department of Chemistry, The University of Texas of the Permian Basin, Odessa, Texas 79762, USA
| | - Wael M. Al-Sawai
- Department of Mathematics & Physics, The University of Texas of the Permian Basin, Odessa, Texas 79762, USA
| |
Collapse
|
12
|
Xin X, Niu X, Liu W, Wang D. Hybrid Solvation Model with First Solvation Shell for Calculation of Solvation Free Energy. Chemphyschem 2020; 21:762-769. [PMID: 32154979 DOI: 10.1002/cphc.202000039] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/14/2020] [Revised: 02/15/2020] [Indexed: 02/03/2023]
Abstract
We present a hybrid solvation model with first solvation shell to calculate solvation free energies. This hybrid model combines the quantum mechanics and molecular mechanics methods with the analytical expression based on the Born solvation model to calculate solvation free energies. Based on calculated free energies of solvation and reaction profiles in gas phase, we set up a unified scheme to predict reaction profiles in solution. The predicted solvation free energies and reaction barriers are compared with experimental results for twenty bimolecular nucleophilic substitution reactions. These comparisons show that our hybrid solvation model can predict reliable solvation free energies and reaction barriers for chemical reactions of small molecules in aqueous solution.
Collapse
Affiliation(s)
- Xin Xin
- College of Physics and Electronics, Shandong Normal University, Jinan, Shandong, 250014, China
| | - Xiao Niu
- College of Physics and Electronics, Shandong Normal University, Jinan, Shandong, 250014, China
| | - Wanqi Liu
- College of Physics and Electronics, Shandong Normal University, Jinan, Shandong, 250014, China
| | - Dunyou Wang
- College of Physics and Electronics, Shandong Normal University, Jinan, Shandong, 250014, China
| |
Collapse
|
13
|
Saura P, Röpke M, Gamiz-Hernandez AP, Kaila VRI. Quantum Chemical and QM/MM Models in Biochemistry. Methods Mol Biol 2020; 2022:75-104. [PMID: 31396900 DOI: 10.1007/978-1-4939-9608-7_4] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/25/2023]
Abstract
Quantum chemical (QC) calculations provide a basis for deriving a microscopic understanding of enzymes and photobiological systems. Here we describe how QC models can be used to explore the electronic structure, dynamics, and energetics of biomolecules. We introduce the hybrid quantum mechanics/classical mechanics (QM/MM) approach, where a quantum mechanically described system of interest is embedded in a classically described force field representation of the biochemical surroundings. We also discuss the QM cluster model approach, as well as embedding theories, that provide complementary methodologies to model quantum mechanical effects in biomolecules. The chapter also provides some practical guides for building quantum biochemical models using the quinone reduction catalysis in respiratory complex I and a model reaction in solution as examples.
Collapse
Affiliation(s)
- Patricia Saura
- Department Chemie, Technische Universität München, Garching, Germany
| | - Michael Röpke
- Department Chemie, Technische Universität München, Garching, Germany
| | | | - Ville R I Kaila
- Department Chemie, Technische Universität München, Garching, Germany.
| |
Collapse
|
14
|
Chen J, Chan B, Shao Y, Ho J. How accurate are approximate quantum chemical methods at modelling solute-solvent interactions in solvated clusters? Phys Chem Chem Phys 2020; 22:3855-3866. [PMID: 32022044 PMCID: PMC7394230 DOI: 10.1039/c9cp06792b] [Citation(s) in RCA: 17] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/28/2022]
Abstract
In this paper, the performance of a wide range of DFT methods is assessed for the calculation of interaction energies of thermal clusters of a solute in water. Three different charge states (neutral, proton transfer transition state and zwitterion) of glycine were solvated by 1 to 40 water molecules as sampled from molecular dynamics simulations. While some ab initio composite methods that employ insufficiently large basis sets incurred significant errors even for a cluster containing only 5 water molecules relative to the W1X-2 benchmark, the DLPNO-CCSD(T)/CBS and DSD-PBEP86 (triple zeta basis set) levels of theory predicted very accurate interaction energies. These levels of theory were used to benchmark the performance of 16 density functionals from different rungs of Jacob's Ladder. Of the Rung 4 functionals examined, the ωB97M-V and ωB97X-V functionals stood out for predicting absolute interaction energies in 40-water clusters with mean absolute deviations (MAD) ∼4 kJ mol-1. The B3LYP-D3(BJ) functional performed exceptionally well with a MAD ∼1.7 kJ mol-1 and is the overall best performing method. Calculations of relative interaction energies allow for cancellation of systematic errors, including basis set truncation and superposition errors, and the ωB97M-V and B3LYP-D3(BJ) double zeta basis set calculations yielded relative interaction energies that are within ∼3 kJ mol-1 of the benchmark. The ONIOM approximation provides another strategy for accelerating the calculation of accurate absolute interaction energies provided that the calculations have converged with respect to the size of the "high-level-layer".
Collapse
Affiliation(s)
- Junbo Chen
- School of Chemistry, University of New South Wales, Sydney, NSW 2052, Australia.
| | - Bun Chan
- Graduate School of Engineering, Nagasaki University, Bunkyo-Machi 1-14, Nagasaki 852-8521, Japan.
| | - Yihan Shao
- Department of Chemistry and Biochemistry, University of Oklahoma, Norman, Oklahoma 73019, USA
| | - Junming Ho
- School of Chemistry, University of New South Wales, Sydney, NSW 2052, Australia.
| |
Collapse
|
15
|
Kelly BD, Smith WR. Alchemical Hydration Free-Energy Calculations Using Molecular Dynamics with Explicit Polarization and Induced Polarity Decoupling: An On–the–Fly Polarization Approach. J Chem Theory Comput 2020; 16:1146-1161. [DOI: 10.1021/acs.jctc.9b01139] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022]
Affiliation(s)
- Braden D. Kelly
- Department of Mathematics and Statistics, University of Guelph, Guelph, ON N1G 2W1, Canada
| | - William R. Smith
- Department of Mathematics and Statistics, University of Guelph, Guelph, ON N1G 2W1, Canada
- Department of Chemical Engineering, University of Waterloo, Waterloo, ON N2L 3G1, Canada
- Faculty of Science, Ontario Tech University, Oshawa, ON L1H 7K4, Canada
| |
Collapse
|
16
|
Collins MA, Ho J. Accelerating the Calculation of Solute-Solvent Interaction Energies through Systematic Molecular Fragmentation. J Phys Chem A 2019; 123:8476-8484. [PMID: 31509417 DOI: 10.1021/acs.jpca.9b06041] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
The method of systematic molecular fragmentation by annihilation (SMFA) is modified to apply to the interaction energy between a solute and solvent, where the solute is a pair of reacting molecules. For NH3 + CH3Cl as the solute, it is shown that SMFA can estimate (to chemical accuracy) the average binding energy of the solute in large water clusters containing up to 160 water molecules, at an appropriate level of electronic structure theory. The SMFA calculation can be carried out in a computation time that makes it feasible to estimate the solvation contribution to free energies of activation and reaction by ensemble averaging.
Collapse
Affiliation(s)
- Michael A Collins
- Research School of Chemistry , Australian National University , Canberra , ACT 2601 , Australia
| | - Junming Ho
- School of Chemistry , University of New South Wales , Sydney , NSW 2052 , Australia
| |
Collapse
|
17
|
Giese TJ, York DM. Development of a Robust Indirect Approach for MM → QM Free Energy Calculations That Combines Force-Matched Reference Potential and Bennett's Acceptance Ratio Methods. J Chem Theory Comput 2019; 15:5543-5562. [PMID: 31507179 DOI: 10.1021/acs.jctc.9b00401] [Citation(s) in RCA: 33] [Impact Index Per Article: 6.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/20/2022]
Abstract
We use the PBE0/6-31G* density functional method to perform ab initio quantum mechanical/molecular mechanical (QM/MM) molecular dynamics (MD) simulations under periodic boundary conditions with rigorous electrostatics using the ambient potential composite Ewald method in order to test the convergence of MM → QM/MM free energy corrections for the prediction of 17 small-molecule solvation free energies and eight ligand binding free energies to T4 lysozyme. The "indirect" thermodynamic cycle for calculating free energies is used to explore whether a series of reference potentials improve the statistical quality of the predictions. Specifically, we construct a series of reference potentials that optimize a molecular mechanical (MM) force field's parameters to reproduce the ab initio QM/MM forces from a QM/MM simulation. The optimizations form a systematic progression of successively expanded parameters that include bond, angle, dihedral, and charge parameters. For each reference potential, we calculate benchmark quality reference values for the MM → QM/MM correction by performing the mixed MM and QM/MM Hamiltonians at 11 intermediate states, each for 200 ps. We then compare forward and reverse application of Zwanzig's relation, thermodynamic integration (TI), and Bennett's acceptance ratio (BAR) methods as a function of reference potential, simulation time, and the number of simulated intermediate states. We find that Zwanzig's equation is inadequate unless a large number of intermediate states are explicitly simulated. The TI and BAR mean signed errors are very small even when only the end-state simulations are considered, and the standard deviations of the TI and BAR errors are decreased by choosing a reference potential that optimizes the bond and angle parameters. We find a robust approach for the data sets of fairly rigid molecules considered here is to use bond + angle reference potential together with the end-state-only BAR analysis. This requires QM/MM simulations to be performed in order to generate reference data to parametrize the bond + angle reference potential, and then this same simulation serves a dual purpose as the full QM/MM end state. The convergence of the results with respect to time suggests that computational resources may be used more efficiently by running multiple simulations for no more than 50 ps, rather than running one long simulation.
Collapse
Affiliation(s)
- Timothy J Giese
- Laboratory for Biomolecular Simulation Research, Center for Integrative Proteomics Research and Department of Chemistry and Chemical Biology , Rutgers University , Piscataway , New Jersey 08854-8087 , United States
| | - Darrin M York
- Laboratory for Biomolecular Simulation Research, Center for Integrative Proteomics Research and Department of Chemistry and Chemical Biology , Rutgers University , Piscataway , New Jersey 08854-8087 , United States
| |
Collapse
|
18
|
Kaliakin DS, Fedorov DG, Alexeev Y, Varganov SA. Locating Minimum Energy Crossings of Different Spin States Using the Fragment Molecular Orbital Method. J Chem Theory Comput 2019; 15:6074-6084. [DOI: 10.1021/acs.jctc.9b00641] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/13/2023]
Affiliation(s)
- Danil S. Kaliakin
- Department of Chemistry, University of Nevada, Reno, 1664 N. Virginia Street, Reno, Nevada 89557-0216, United States
| | - Dmitri G. Fedorov
- Research Center for Computational Design of Advanced Functional Materials (CD-FMat), National Institute of Advanced Industrial Science and Technology (AIST), Central 2, Umezono 1-1-1, Tsukuba 305-8568, Japan
| | - Yuri Alexeev
- Computational Science Division and Argonne Leadership Computing Facility, Argonne National Laboratory, Argonne, Illinois 60439, United States
| | - Sergey A. Varganov
- Department of Chemistry, University of Nevada, Reno, 1664 N. Virginia Street, Reno, Nevada 89557-0216, United States
| |
Collapse
|
19
|
Hudson PS, Woodcock HL, Boresch S. Use of Interaction Energies in QM/MM Free Energy Simulations. J Chem Theory Comput 2019; 15:4632-4645. [PMID: 31142113 DOI: 10.1021/acs.jctc.9b00084] [Citation(s) in RCA: 21] [Impact Index Per Article: 4.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/2022]
Abstract
The use of the most accurate (i.e., QM or QM/MM) levels of theory for free energy simulations (FES) is typically not possible. Primarily, this is because the computational cost associated with the extensive configurational sampling needed for converging FES is prohibitive. To ensure the feasibility of QM-based FES, the "indirect" approach is generally taken, necessitating a free energy calculation between the MM and QM/MM potential energy surfaces. Ideally, this step is performed with standard free energy perturbation (Zwanzig's equation) as it only requires simulations be carried out at the low level of theory; however, work from several groups over the past few years has conclusively shown that Zwanzig's equation is ill-suited to this task. As such, many approximations have arisen to mitigate difficulties with Zwanzig's equation. One particularly popular notion is that the convergence of Zwanzig's equation can be improved by using interaction energy differences instead of total energy differences. Although problematic numerical fluctuations (a major problem when using Zwanzig's equation) are indeed reduced, our results and analysis demonstrate that this "interaction energy approximation" (IEA) is theoretically incorrect, and the implicit approximation invoked is spurious at best. Herein, we demonstrate this via solvation free energy calculations using IEA from two different low levels of theory to the same target high level. Results from this proof-of-concept consistently yield the wrong results, deviating by ∼1.5 kcal/mol from the rigorously obtained value.
Collapse
Affiliation(s)
- Phillip S Hudson
- Department of Chemistry , University of South Florida , 4202 East Fowler Avenue, CHE205 , Tampa , Florida 33620-5250 , United States.,Laboratory of Computational Biology , National Institutes of Health, National Heart, Lung and Blood Institute , 12 South Drive, Rm 3053 , Bethesda , Maryland 20892-5690 , United States
| | - H Lee Woodcock
- Department of Chemistry , University of South Florida , 4202 East Fowler Avenue, CHE205 , Tampa , Florida 33620-5250 , United States
| | - Stefan Boresch
- Faculty of Chemistry, Department of Computational Biological Chemistry , University of Vienna , Währingerstraße 17 , Vienna A-1090 , Austria
| |
Collapse
|
20
|
Affiliation(s)
- Ida-Marie Høyvik
- Department of Chemistry, The Norwegian University of Science and Technology, Trondheim, Norway
| |
Collapse
|
21
|
Chen J, Shao Y, Ho J. Are Explicit Solvent Models More Accurate than Implicit Solvent Models? A Case Study on the Menschutkin Reaction. J Phys Chem A 2019; 123:5580-5589. [DOI: 10.1021/acs.jpca.9b03995] [Citation(s) in RCA: 20] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Junbo Chen
- School of Chemistry, University of New South Wales, Sydney, NSW 2052, Australia
| | - Yihan Shao
- Department of Chemistry and Biochemistry, University of Oklahoma, Norman, Oklahoma 73019, United States
| | - Junming Ho
- School of Chemistry, University of New South Wales, Sydney, NSW 2052, Australia
| |
Collapse
|
22
|
Duster AW, Garza CM, Aydintug BO, Negussie MB, Lin H. Adaptive Partitioning QM/MM for Molecular Dynamics Simulations: 6. Proton Transport through a Biological Channel. J Chem Theory Comput 2019; 15:892-905. [DOI: 10.1021/acs.jctc.8b01128] [Citation(s) in RCA: 20] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Adam W. Duster
- Chemistry Department, CB 194, University of Colorado, Denver, Colorado 80217, United States
| | - Christina M. Garza
- Chemistry Department, CB 194, University of Colorado, Denver, Colorado 80217, United States
| | - Baris O. Aydintug
- Chemistry Department, CB 194, University of Colorado, Denver, Colorado 80217, United States
| | - Mikias B. Negussie
- Chemistry Department, CB 194, University of Colorado, Denver, Colorado 80217, United States
| | - Hai Lin
- Chemistry Department, CB 194, University of Colorado, Denver, Colorado 80217, United States
| |
Collapse
|
23
|
Huggins DJ, Biggin PC, Dämgen MA, Essex JW, Harris SA, Henchman RH, Khalid S, Kuzmanic A, Laughton CA, Michel J, Mulholland AJ, Rosta E, Sansom MSP, van der Kamp MW. Biomolecular simulations: From dynamics and mechanisms to computational assays of biological activity. WILEY INTERDISCIPLINARY REVIEWS-COMPUTATIONAL MOLECULAR SCIENCE 2018. [DOI: 10.1002/wcms.1393] [Citation(s) in RCA: 71] [Impact Index Per Article: 11.8] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/01/2023]
Affiliation(s)
- David J. Huggins
- TCM Group, Cavendish Laboratory University of Cambridge Cambridge UK
- Unilever Centre, Department of Chemistry University of Cambridge Cambridge UK
- Department of Physiology and Biophysics Weill Cornell Medical College New York NY
| | | | - Marc A. Dämgen
- Department of Biochemistry University of Oxford Oxford UK
| | - Jonathan W. Essex
- School of Chemistry University of Southampton Southampton UK
- Institute for Life Sciences University of Southampton Southampton UK
| | - Sarah A. Harris
- School of Physics and Astronomy University of Leeds Leeds UK
- Astbury Centre for Structural and Molecular Biology University of Leeds Leeds UK
| | - Richard H. Henchman
- Manchester Institute of Biotechnology The University of Manchester Manchester UK
- School of Chemistry The University of Manchester Oxford UK
| | - Syma Khalid
- School of Chemistry University of Southampton Southampton UK
- Institute for Life Sciences University of Southampton Southampton UK
| | | | - Charles A. Laughton
- School of Pharmacy University of Nottingham Nottingham UK
- Centre for Biomolecular Sciences University of Nottingham Nottingham UK
| | - Julien Michel
- EaStCHEM school of Chemistry University of Edinburgh Edinburgh UK
| | - Adrian J. Mulholland
- Centre of Computational Chemistry, School of Chemistry University of Bristol Bristol UK
| | - Edina Rosta
- Department of Chemistry King's College London London UK
| | | | - Marc W. van der Kamp
- Centre of Computational Chemistry, School of Chemistry University of Bristol Bristol UK
- School of Biochemistry, Biomedical Sciences Building University of Bristol Bristol UK
| |
Collapse
|
24
|
Duster AW, Wang CH, Lin H. Adaptive QM/MM for Molecular Dynamics Simulations: 5. On the Energy-Conserved Permuted Adaptive-Partitioning Schemes. Molecules 2018; 23:E2170. [PMID: 30154373 PMCID: PMC6225285 DOI: 10.3390/molecules23092170] [Citation(s) in RCA: 16] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/24/2018] [Revised: 08/23/2018] [Accepted: 08/24/2018] [Indexed: 11/16/2022] Open
Abstract
In combined quantum-mechanical/molecular-mechanical (QM/MM) dynamics simulations, the adaptive-partitioning (AP) schemes reclassify atoms on-the-fly as QM or MM in a smooth manner. This yields a mobile QM subsystem with contents that are continuously updated as needed. Here, we tailor the Hamiltonian adaptive many-body correction (HAMBC) proposed by Boreboom et al. [J. Chem. Theory Comput.2016, 12, 3441] to the permuted AP (PAP) scheme. The treatments lead to the HAMBC-PAP method (HPAP), which both conserves energy and produces accurate solvation structures in the test of "water-in-water" model system.
Collapse
Affiliation(s)
- Adam W Duster
- Department of Chemistry, University of Colorado Denver, Denver, CO 80217, USA.
| | - Chun-Hung Wang
- Department of Chemistry, University of Colorado Denver, Denver, CO 80217, USA.
| | - Hai Lin
- Department of Chemistry, University of Colorado Denver, Denver, CO 80217, USA.
| |
Collapse
|
25
|
Saielli G. Computational Spectroscopy of Ionic Liquids for Bulk Structure Elucidation. ADVANCED THEORY AND SIMULATIONS 2018. [DOI: 10.1002/adts.201800084] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/21/2023]
Affiliation(s)
- Giacomo Saielli
- CNR Institute on Membrane Technology; Unit of Padova; Via Marzolo 1-35131 Padova Italy
- Department of Chemical Sciences; University of Padova; Via Marzolo 1-35131 Padova Italy
| |
Collapse
|
26
|
Zhang YJ, Khorshidi A, Kastlunger G, Peterson AA. The potential for machine learning in hybrid QM/MM calculations. J Chem Phys 2018; 148:241740. [DOI: 10.1063/1.5029879] [Citation(s) in RCA: 30] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/11/2022] Open
Affiliation(s)
- Yin-Jia Zhang
- Department of Chemistry, Brown University, Providence, Rhode Island 02912, USA
| | - Alireza Khorshidi
- School of Engineering, Brown University, Providence, Rhode Island 02912, USA
| | - Georg Kastlunger
- School of Engineering, Brown University, Providence, Rhode Island 02912, USA
| | - Andrew A. Peterson
- School of Engineering, Brown University, Providence, Rhode Island 02912, USA
| |
Collapse
|
27
|
Hofer TS, Hünenberger PH. Absolute proton hydration free energy, surface potential of water, and redox potential of the hydrogen electrode from first principles: QM/MM MD free-energy simulations of sodium and potassium hydration. J Chem Phys 2018; 148:222814. [DOI: 10.1063/1.5000799] [Citation(s) in RCA: 60] [Impact Index Per Article: 10.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/01/2023] Open
Affiliation(s)
- Thomas S. Hofer
- Theoretical Chemistry Division, Institute of General, Inorganic and Theoretical Chemistry, Centre for Chemistry and Biomedicine, University of Innsbruck, Innrain 80-82, A-6020 Innsbruck, Austria
| | | |
Collapse
|
28
|
Steinmann C, Olsson MA, Ryde U. Relative Ligand-Binding Free Energies Calculated from Multiple Short QM/MM MD Simulations. J Chem Theory Comput 2018; 14:3228-3237. [PMID: 29768915 DOI: 10.1021/acs.jctc.8b00081] [Citation(s) in RCA: 19] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022]
Abstract
We have devised a new efficient approach to compute combined quantum mechanical (QM) and molecular mechanical (MM, i.e. QM/MM) ligand-binding relative free energies. Our method employs the reference-potential approach with free-energy perturbation both at the MM level (between the two ligands) and from MM to QM/MM (for each ligand). To ensure that converged results are obtained for the MM → QM/MM perturbations, explicit QM/MM molecular dynamics (MD) simulations are performed with two intermediate mixed states. To speed up the calculations, we utilize the fact that the phase space can be extensively sampled at the MM level. Therefore, we run many short QM/MM MD simulations started from snapshots of the MM simulations, instead of a single long simulation. As a test case, we study the binding of nine cyclic carboxylate ligands to the octa-acid deep cavitand. Only the ligand is in the QM system, treated with the semiempirical PM6-DH+ method. We show that for eight of the ligands, we obtain well converged results with short MD simulations (1-15 ps). However, in one case, the convergence is slower (∼50 ps) owing to a mismatch between the conformational preferences of the MM and QM/MM potentials. We test the effect of initial minimization, the need of equilibration, and how many independent simulations are needed to reach a certain precision. The results show that the present approach is about four times faster than using standard MM → QM/MM free-energy perturbations with the same accuracy and precision.
Collapse
Affiliation(s)
- Casper Steinmann
- Department of Chemistry and Bioscience , Aalborg University , Frederik Bajers Vej 7H , DK-9220 Aalborg , Denmark.,Department of Theoretical Chemistry , Lund University , Chemical Centre , P.O. Box 124, SE-221 00 Lund , Sweden
| | - Martin A Olsson
- Department of Theoretical Chemistry , Lund University , Chemical Centre , P.O. Box 124, SE-221 00 Lund , Sweden
| | - Ulf Ryde
- Department of Theoretical Chemistry , Lund University , Chemical Centre , P.O. Box 124, SE-221 00 Lund , Sweden
| |
Collapse
|
29
|
Boulanger E, Harvey JN. QM/MM methods for free energies and photochemistry. Curr Opin Struct Biol 2018; 49:72-76. [PMID: 29414514 DOI: 10.1016/j.sbi.2018.01.003] [Citation(s) in RCA: 32] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/18/2017] [Revised: 12/17/2017] [Accepted: 01/02/2018] [Indexed: 11/27/2022]
Abstract
Hybrid computational methods describing a small region of a biomolecular system with quantum mechanics and the bulk with molecular mechanics, referred to as QM/MM methods, are now a central part of computational biochemistry. This review considers developments in the QM/MM approach that make it easier to calculate free energies using accurate QM-based potential energy expressions. We also describe techniques to treat electronic coupling between the core region and the MM environment. Polarizability of the protein matrix is important but so is electronic coupling. Applications of these new methods, especially to photochemistry, are discussed.
Collapse
Affiliation(s)
- Eliot Boulanger
- Department of Chemistry, KU Leuven, Celestijnenlaan 200F, B-3001 Leuven, Belgium.
| | - Jeremy N Harvey
- Department of Chemistry, KU Leuven, Celestijnenlaan 200F, B-3001 Leuven, Belgium.
| |
Collapse
|
30
|
Ainsley J, Lodola A, Mulholland AJ, Christov CZ, Karabencheva-Christova TG. Combined Quantum Mechanics and Molecular Mechanics Studies of Enzymatic Reaction Mechanisms. ADVANCES IN PROTEIN CHEMISTRY AND STRUCTURAL BIOLOGY 2018; 113:1-32. [PMID: 30149903 DOI: 10.1016/bs.apcsb.2018.07.001] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 02/02/2023]
Abstract
The combined quantum mechanics/molecular mechanics (QM/MM) methods have become a valuable tool in computational biochemistry and received versatile applications for studying the reaction mechanisms of enzymes. The approach combines the calculations of the electronic structure of the active site by QM, with modeling of the protein environment using MM force field, which allows the long-range electrostatics and steric effects on the enzyme reactivity to be accounted for. In this review, we review some key theoretical and computational aspects of the method and we also present some applications to particular enzymatic reactions such as tryptophan-7-halogenase, cyclooxygenase-1, and the epidermal growth factor receptor.
Collapse
Affiliation(s)
- Jon Ainsley
- Department of Applied Sciences, Faculty of Health and Life Sciences, Northumbria University, Newcastle upon Tyne, United Kingdom
| | | | - Adrian J Mulholland
- Centre for Computational Chemistry, School of Chemistry, University of Bristol, Bristol, United Kingdom
| | - Christo Z Christov
- Department of Chemistry, Michigan Technological University, Houghton, MI, United States.
| | | |
Collapse
|
31
|
DasGupta D, Mandalaparthy V, Jayaram B. A component analysis of the free energies of folding of 35 proteins: A consensus view on the thermodynamics of folding at the molecular level. J Comput Chem 2017; 38:2791-2801. [PMID: 28940242 DOI: 10.1002/jcc.25072] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/02/2017] [Revised: 07/27/2017] [Accepted: 09/01/2017] [Indexed: 02/05/2023]
Abstract
What factors favor protein folding? This is a textbook question. Parsing the experimental free energies of folding/unfolding into diverse enthalpic and entropic components of solute and solvent favoring or disfavoring folding is not an easy task. In this study, we present a computational protocol for estimating the free energy contributors to protein folding semi-quantitatively using ensembles of unfolded and native states generated via molecular dynamics simulations. We tested the methodology on 35 proteins with diverse structural motifs and sizes and found that the calculated free energies correlate well with experiment (correlation coefficient ∼ 0.85), enabling us to develop a consensus view of the energetics of folding. As a more sensitive test of the methodology, we also investigated the free energies of folding of an additional 33 single point mutants and obtained a correlation coefficient of 0.8. A notable observation is that the folding free energy components appear to carry signatures of the fold (SCOP classification) of the protein. © 2017 Wiley Periodicals, Inc.
Collapse
Affiliation(s)
- Debarati DasGupta
- Department of Chemistry, Indian Institute of Technology, New Delhi, 110016, India.,Supercomputing Facility for Bioinformatics and Computational Biology, Indian Institute of Technology, New Delhi, 110016, India
| | - Varun Mandalaparthy
- Department of Chemistry, Indian Institute of Technology, New Delhi, 110016, India
| | - Bhyravabhotla Jayaram
- Department of Chemistry, Indian Institute of Technology, New Delhi, 110016, India.,Supercomputing Facility for Bioinformatics and Computational Biology, Indian Institute of Technology, New Delhi, 110016, India.,Kusuma School of Biological Sciences, Indian Institute of Technology, New Delhi, 110016, India
| |
Collapse
|
32
|
Duboué-Dijon E, Pluhařová E, Domin D, Sen K, Fogarty AC, Chéron N, Laage D. Coupled Valence-Bond State Molecular Dynamics Description of an Enzyme-Catalyzed Reaction in a Non-Aqueous Organic Solvent. J Phys Chem B 2017; 121:7027-7041. [PMID: 28675789 DOI: 10.1021/acs.jpcb.7b03102] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Enzymes are widely used in nonaqueous solvents to catalyze non-natural reactions. While experimental measurements showed that the solvent nature has a strong effect on the reaction kinetics, the molecular details of the catalytic mechanism in nonaqueous solvents have remained largely elusive. Here we study the transesterification reaction catalyzed by the paradigm subtilisin Carlsberg serine protease in an organic apolar solvent. The rate-limiting acylation step involves a proton transfer between active-site residues and the nucleophilic attack of the substrate to form a tetrahedral intermediate. We design the first coupled valence-bond state model that simultaneously describes both reactions in the enzymatic active site. We develop a new systematic procedure to parametrize this model on high-level ab initio QM/MM free energy calculations that account for the molecular details of the active site and for both substrate and protein conformational fluctuations. Our calculations show that the reaction energy barrier changes dramatically with the solvent and protein conformational fluctuations. We find that the mechanism of the tetrahedral intermediate formation during the acylation step is similar to that determined under aqueous conditions, and that the proton transfer and nucleophilic attack reactions occur concertedly. We identify the reaction coordinate to be mostly due to the rearrangement of some residual water molecules close to the active site.
Collapse
Affiliation(s)
- Elise Duboué-Dijon
- École Normale Supérieure - PSL Research University , Département de Chimie, Sorbonne Universités - UPMC Univ Paris 06, CNRS UMR 8640 PASTEUR, 24 rue Lhomond, 75005 Paris, France
| | - Eva Pluhařová
- École Normale Supérieure - PSL Research University , Département de Chimie, Sorbonne Universités - UPMC Univ Paris 06, CNRS UMR 8640 PASTEUR, 24 rue Lhomond, 75005 Paris, France
| | - Dominik Domin
- École Normale Supérieure - PSL Research University , Département de Chimie, Sorbonne Universités - UPMC Univ Paris 06, CNRS UMR 8640 PASTEUR, 24 rue Lhomond, 75005 Paris, France
| | - Kakali Sen
- École Normale Supérieure - PSL Research University , Département de Chimie, Sorbonne Universités - UPMC Univ Paris 06, CNRS UMR 8640 PASTEUR, 24 rue Lhomond, 75005 Paris, France
| | - Aoife C Fogarty
- École Normale Supérieure - PSL Research University , Département de Chimie, Sorbonne Universités - UPMC Univ Paris 06, CNRS UMR 8640 PASTEUR, 24 rue Lhomond, 75005 Paris, France
| | - Nicolas Chéron
- École Normale Supérieure - PSL Research University , Département de Chimie, Sorbonne Universités - UPMC Univ Paris 06, CNRS UMR 8640 PASTEUR, 24 rue Lhomond, 75005 Paris, France
| | - Damien Laage
- École Normale Supérieure - PSL Research University , Département de Chimie, Sorbonne Universités - UPMC Univ Paris 06, CNRS UMR 8640 PASTEUR, 24 rue Lhomond, 75005 Paris, France
| |
Collapse
|
33
|
Ujiie Y, Tanaka W, Hanaoka K, Harada R, Kayanuma M, Shoji M, Murakawa T, Ishida T, Shigeta Y, Hayashi H. Molecular Mechanism of the Reaction Specificity in Threonine Synthase: Importance of the Substrate Conformations. J Phys Chem B 2017; 121:5536-5543. [PMID: 28489381 DOI: 10.1021/acs.jpcb.7b02932] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Threonine synthase (ThrS) catalyzes the final chemical reaction of l-threonine biosynthesis from its precursor, O-phospho-l-homoserine. As the phosphate ion generated in its former half reaction assists its latter reaction, ThrS is recognized as one of the best examples of product-assisted catalysis. In our previous QM/MM study, the chemical reactions for the latter half reactions, which are critical for the product-assisted catalysis, were revealed. However, accurate free energy changes caused by the conformational ensembles and entrance of water molecules into the active site are unknown. In the present study, by performing long-time scale MD simulations, the free energy changes by the divalent anions (phosphate or sulfate ions) and conformational states of the intermediate states were theoretically investigated. We found that the calculated free energy double differences are in good agreement with the experimental results. We also revealed that the phosphate ion contributes to forming hydrogen bonds that are suitable for the main reaction progress. This means that the conformation of the active site amino acid residues and the substrate, and hence, the tunable catalysis, are controlled by the product phosphate ion, and this clearly demonstrates a molecular mechanism of the product-assisted catalysis in ThrS.
Collapse
Affiliation(s)
- Yuzuru Ujiie
- Graduate School of Pure and Applied Sciences, University of Tsukuba , Tennodai 1-1-1, Tsukuba, Ibaraki 305-8571, Japan
| | - Wataru Tanaka
- Graduate School of Pure and Applied Sciences, University of Tsukuba , Tennodai 1-1-1, Tsukuba, Ibaraki 305-8571, Japan
| | - Kyohei Hanaoka
- Graduate School of Pure and Applied Sciences, University of Tsukuba , Tennodai 1-1-1, Tsukuba, Ibaraki 305-8571, Japan
| | - Ryuhei Harada
- Center for Computational Sciences, University of Tsukuba , Tennodai 1-1-1, Tsukuba, Ibaraki 305-8577, Japan
| | - Megumi Kayanuma
- Center for Computational Sciences, University of Tsukuba , Tennodai 1-1-1, Tsukuba, Ibaraki 305-8577, Japan
| | - Mitsuo Shoji
- Graduate School of Pure and Applied Sciences, University of Tsukuba , Tennodai 1-1-1, Tsukuba, Ibaraki 305-8571, Japan.,Center for Computational Sciences, University of Tsukuba , Tennodai 1-1-1, Tsukuba, Ibaraki 305-8577, Japan
| | - Takeshi Murakawa
- Department of Biochemistry, Osaka Medical College , Takatsuki 569-8686, Japan
| | - Toyokazu Ishida
- Research Center for Computational Design of Advanced Functional Materials (CD-FMat), National Institute of Advanced Industrial Science and Technology (AIST) , Tsukuba Central 2, 1-1-1 Umezono, Tsukuba, Ibaraki 305-8568, Japan
| | - Yasuteru Shigeta
- Graduate School of Pure and Applied Sciences, University of Tsukuba , Tennodai 1-1-1, Tsukuba, Ibaraki 305-8571, Japan.,Center for Computational Sciences, University of Tsukuba , Tennodai 1-1-1, Tsukuba, Ibaraki 305-8577, Japan
| | - Hideyuki Hayashi
- Department of Chemistry, Osaka Medical College , Takatsuki 569-8686, Japan
| |
Collapse
|
34
|
Hayashi S, Uchida Y, Hasegawa T, Higashi M, Kosugi T, Kamiya M. QM/MM Geometry Optimization on Extensive Free-Energy Surfaces for Examination of Enzymatic Reactions and Design of Novel Functional Properties of Proteins. Annu Rev Phys Chem 2017; 68:135-154. [DOI: 10.1146/annurev-physchem-052516-050827] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/09/2022]
Affiliation(s)
- Shigehiko Hayashi
- Department of Chemistry, Graduate School of Science, Kyoto University, Kyoto 606-8502, Japan;, , ,
| | - Yoshihiro Uchida
- Department of Chemistry, Graduate School of Science, Kyoto University, Kyoto 606-8502, Japan;, , ,
| | - Taisuke Hasegawa
- Department of Chemistry, Graduate School of Science, Kyoto University, Kyoto 606-8502, Japan;, , ,
| | - Masahiro Higashi
- Department of Chemistry, Biology and Marine Science, University of the Ryukyus, Nishihara, Okinawa 903-0213, Japan
| | - Takahiro Kosugi
- Research Center of Integrative Molecular Systems, Institute for Molecular Science, and Department of Structural Molecular Science, School of Physical Science, The Graduate University for Advanced Studies (SOKENDAI), Okazaki 444-8585, Japan
| | - Motoshi Kamiya
- Department of Chemistry, Graduate School of Science, Kyoto University, Kyoto 606-8502, Japan;, , ,
| |
Collapse
|
35
|
Nakano H, Sato H. Introducing the mean field approximation to CDFT/MMpol method: Statistically converged equilibrium and nonequilibrium free energy calculation for electron transfer reactions in condensed phases. J Chem Phys 2017; 146:154101. [DOI: 10.1063/1.4979895] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/04/2023] Open
Affiliation(s)
- Hiroshi Nakano
- Department of Molecular Engineering, Graduate School of Engineering, Kyoto University, Kyoto 615-8510, Japan
- Elements Strategy Initiative for Catalysts and Batteries, Kyoto University, Japan
| | - Hirofumi Sato
- Department of Molecular Engineering, Graduate School of Engineering, Kyoto University, Kyoto 615-8510, Japan
- Elements Strategy Initiative for Catalysts and Batteries, Kyoto University, Japan
| |
Collapse
|
36
|
Duster AW, Wang C, Garza CM, Miller DE, Lin H. Adaptive quantum/molecular mechanics: what have we learned, where are we, and where do we go from here? WILEY INTERDISCIPLINARY REVIEWS-COMPUTATIONAL MOLECULAR SCIENCE 2017. [DOI: 10.1002/wcms.1310] [Citation(s) in RCA: 43] [Impact Index Per Article: 6.1] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/11/2023]
Affiliation(s)
- Adam W. Duster
- Chemistry Department University of Colorado Denver Denver CO USA
| | - Chun‐Hung Wang
- Chemistry Department University of Colorado Denver Denver CO USA
| | | | | | - Hai Lin
- Chemistry Department University of Colorado Denver Denver CO USA
| |
Collapse
|
37
|
Olsson MA, Ryde U. Comparison of QM/MM Methods To Obtain Ligand-Binding Free Energies. J Chem Theory Comput 2017; 13:2245-2253. [DOI: 10.1021/acs.jctc.6b01217] [Citation(s) in RCA: 36] [Impact Index Per Article: 5.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/18/2022]
Affiliation(s)
- Martin A. Olsson
- Department of Theoretical
Chemistry, Lund University, Chemical Centre, P.O. Box 124, SE-221
00 Lund, Sweden
| | - Ulf Ryde
- Department of Theoretical
Chemistry, Lund University, Chemical Centre, P.O. Box 124, SE-221
00 Lund, Sweden
| |
Collapse
|
38
|
Dong L, Liu Y. Comparative studies of the catalytic mechanisms of two chorismatases: CH-fkbo and CH-Hyg5. Proteins 2017; 85:1146-1158. [DOI: 10.1002/prot.25279] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/11/2016] [Revised: 02/13/2017] [Accepted: 02/21/2017] [Indexed: 01/11/2023]
Affiliation(s)
- Lihua Dong
- School of Chemistry and Chemical Engineering; Shandong University; Jinan Shandong 250100 China
- School of Chemistry and Chemical Engineering; Qilu Normal University; Jinan Shandong 250013 China
| | - Yongjun Liu
- School of Chemistry and Chemical Engineering; Shandong University; Jinan Shandong 250100 China
| |
Collapse
|
39
|
Vaissier V, Van Voorhis T. Adiabatic Approximation in Explicit Solvent Models of RedOx Chemistry. J Chem Theory Comput 2016; 12:5111-5116. [DOI: 10.1021/acs.jctc.6b00746] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Valérie Vaissier
- Department of Chemistry, Massachusetts Institute of Technology, 77 Massachusetts Avenue, Cambridge, Massachusetts 02139, United States
| | - Troy Van Voorhis
- Department of Chemistry, Massachusetts Institute of Technology, 77 Massachusetts Avenue, Cambridge, Massachusetts 02139, United States
| |
Collapse
|
40
|
Olsson MA, Söderhjelm P, Ryde U. Converging ligand-binding free energies obtained with free-energy perturbations at the quantum mechanical level. J Comput Chem 2016; 37:1589-600. [PMID: 27117350 PMCID: PMC5074236 DOI: 10.1002/jcc.24375] [Citation(s) in RCA: 40] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/02/2016] [Revised: 03/08/2016] [Accepted: 03/09/2016] [Indexed: 12/15/2022]
Abstract
In this article, the convergence of quantum mechanical (QM) free-energy simulations based on molecular dynamics simulations at the molecular mechanics (MM) level has been investigated. We have estimated relative free energies for the binding of nine cyclic carboxylate ligands to the octa-acid deep-cavity host, including the host, the ligand, and all water molecules within 4.5 Å of the ligand in the QM calculations (158-224 atoms). We use single-step exponential averaging (ssEA) and the non-Boltzmann Bennett acceptance ratio (NBB) methods to estimate QM/MM free energy with the semi-empirical PM6-DH2X method, both based on interaction energies. We show that ssEA with cumulant expansion gives a better convergence and uses half as many QM calculations as NBB, although the two methods give consistent results. With 720,000 QM calculations per transformation, QM/MM free-energy estimates with a precision of 1 kJ/mol can be obtained for all eight relative energies with ssEA, showing that this approach can be used to calculate converged QM/MM binding free energies for realistic systems and large QM partitions. © 2016 The Authors. Journal of Computational Chemistry Published by Wiley Periodicals, Inc.
Collapse
Affiliation(s)
- Martin A Olsson
- Department of Theoretical Chemistry, Lund University, Chemical Centre, P. O. Box 124, Lund, SE-221 00, Sweden
| | - Pär Söderhjelm
- Department of Biophysical Chemistry, Lund University, Chemical Centre, P. O. Box 124, Lund, SE-221 00, Sweden
| | - Ulf Ryde
- Department of Theoretical Chemistry, Lund University, Chemical Centre, P. O. Box 124, Lund, SE-221 00, Sweden
| |
Collapse
|
41
|
Kmiecik S, Gront D, Kolinski M, Wieteska L, Dawid AE, Kolinski A. Coarse-Grained Protein Models and Their Applications. Chem Rev 2016; 116:7898-936. [DOI: 10.1021/acs.chemrev.6b00163] [Citation(s) in RCA: 555] [Impact Index Per Article: 69.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/07/2023]
Affiliation(s)
- Sebastian Kmiecik
- Faculty
of Chemistry, University of Warsaw, Pasteura 1, 02-093 Warsaw, Poland
| | - Dominik Gront
- Faculty
of Chemistry, University of Warsaw, Pasteura 1, 02-093 Warsaw, Poland
| | - Michal Kolinski
- Bioinformatics
Laboratory, Mossakowski Medical Research Center of the Polish Academy of Sciences, Pawinskiego 5, 02-106 Warsaw, Poland
| | - Lukasz Wieteska
- Faculty
of Chemistry, University of Warsaw, Pasteura 1, 02-093 Warsaw, Poland
- Department
of Medical Biochemistry, Medical University of Lodz, Mazowiecka 6/8, 92-215 Lodz, Poland
| | | | - Andrzej Kolinski
- Faculty
of Chemistry, University of Warsaw, Pasteura 1, 02-093 Warsaw, Poland
| |
Collapse
|
42
|
Ryde U, Söderhjelm P. Ligand-Binding Affinity Estimates Supported by Quantum-Mechanical Methods. Chem Rev 2016; 116:5520-66. [DOI: 10.1021/acs.chemrev.5b00630] [Citation(s) in RCA: 175] [Impact Index Per Article: 21.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Ulf Ryde
- Department of Theoretical
Chemistry and ‡Department of Biophysical Chemistry, Lund University, Chemical Centre, P.O. Box 124, SE-221 00 Lund, Sweden
| | - Pär Söderhjelm
- Department of Theoretical
Chemistry and ‡Department of Biophysical Chemistry, Lund University, Chemical Centre, P.O. Box 124, SE-221 00 Lund, Sweden
| |
Collapse
|
43
|
Barnard AS. Challenges in modelling nanoparticles for drug delivery. JOURNAL OF PHYSICS. CONDENSED MATTER : AN INSTITUTE OF PHYSICS JOURNAL 2016; 28:023002. [PMID: 26682622 DOI: 10.1088/0953-8984/28/2/023002] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/05/2023]
Abstract
Although there have been significant advances in the fields of theoretical condensed matter and computational physics, when confronted with the complexity and diversity of nanoparticles available in conventional laboratories a number of modeling challenges remain. These challenges are generally shared among application domains, but the impacts of the limitations and approximations we make to overcome them (or circumvent them) can be more significant one area than another. In the case of nanoparticles for drug delivery applications some immediate challenges include the incompatibility of length-scales, our ability to model weak interactions and solvation, the complexity of the thermochemical environment surrounding the nanoparticles, and the role of polydispersivity in determining properties and performance. Some of these challenges can be met with existing technologies, others with emerging technologies including the data-driven sciences; some others require new methods to be developed. In this article we will briefly review some simple methods and techniques that can be applied to these (and other) challenges, and demonstrate some results using nanodiamond-based drug delivery platforms as an exemplar.
Collapse
Affiliation(s)
- Amanda S Barnard
- CSIRO Virtual Nanoscience Laboratory, 343 Royal Parade, Parkville, Victoria 3052, Australia
| |
Collapse
|
44
|
Jia X, Wang M, Shao Y, König G, Brooks BR, Zhang JZH, Mei Y. Calculations of Solvation Free Energy through Energy Reweighting from Molecular Mechanics to Quantum Mechanics. J Chem Theory Comput 2016; 12:499-511. [DOI: 10.1021/acs.jctc.5b00920] [Citation(s) in RCA: 74] [Impact Index Per Article: 9.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/07/2023]
Affiliation(s)
- Xiangyu Jia
- State
Key Laboratory of Precision Spectroscopy and Department of Physics
and Institute of Theoretical and Computational Science, East China Normal University, Shanghai 200062, China
| | - Meiting Wang
- State
Key Laboratory of Precision Spectroscopy and Department of Physics
and Institute of Theoretical and Computational Science, East China Normal University, Shanghai 200062, China
| | - Yihan Shao
- Q-Chem Inc., 6601 Owens Drive, Suite
105, Pleasanton, California 94588, United States
| | - Gerhard König
- Laboratory
of Computational Biology, National Institutes of Health, National Heart, Lung and Blood Institute, 5635 Fishers Lane, T-900 Suite, Rockville, Maryland 20852, United States
| | - Bernard R. Brooks
- Laboratory
of Computational Biology, National Institutes of Health, National Heart, Lung and Blood Institute, 5635 Fishers Lane, T-900 Suite, Rockville, Maryland 20852, United States
| | - John Z. H. Zhang
- State
Key Laboratory of Precision Spectroscopy and Department of Physics
and Institute of Theoretical and Computational Science, East China Normal University, Shanghai 200062, China
- NYU-ECNU Center
for Computational Chemistry at NYU Shanghai, Shanghai 200062, China
| | - Ye Mei
- State
Key Laboratory of Precision Spectroscopy and Department of Physics
and Institute of Theoretical and Computational Science, East China Normal University, Shanghai 200062, China
- NYU-ECNU Center
for Computational Chemistry at NYU Shanghai, Shanghai 200062, China
| |
Collapse
|
45
|
Uritsky N, Shokhen M, Albeck A. Stepwise Versus Concerted Mechanisms in General-Base Catalysis by Serine Proteases. Angew Chem Int Ed Engl 2015; 55:1680-4. [DOI: 10.1002/anie.201507772] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/19/2015] [Indexed: 01/08/2023]
Affiliation(s)
- Neta Uritsky
- The Julius Spokojny Bioorganic Chemistry Laboratory; Department of Chemistry; Bar Ilan University; Ramat Gan 5290002 Israel
| | - Michael Shokhen
- The Julius Spokojny Bioorganic Chemistry Laboratory; Department of Chemistry; Bar Ilan University; Ramat Gan 5290002 Israel
| | - Amnon Albeck
- The Julius Spokojny Bioorganic Chemistry Laboratory; Department of Chemistry; Bar Ilan University; Ramat Gan 5290002 Israel
| |
Collapse
|
46
|
Uritsky N, Shokhen M, Albeck A. Stepwise Versus Concerted Mechanisms in General-Base Catalysis by Serine Proteases. Angew Chem Int Ed Engl 2015. [DOI: 10.1002/ange.201507772] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/08/2022]
Affiliation(s)
- Neta Uritsky
- The Julius Spokojny Bioorganic Chemistry Laboratory; Department of Chemistry; Bar Ilan University; Ramat Gan 5290002 Israel
| | - Michael Shokhen
- The Julius Spokojny Bioorganic Chemistry Laboratory; Department of Chemistry; Bar Ilan University; Ramat Gan 5290002 Israel
| | - Amnon Albeck
- The Julius Spokojny Bioorganic Chemistry Laboratory; Department of Chemistry; Bar Ilan University; Ramat Gan 5290002 Israel
| |
Collapse
|
47
|
|