1
|
Elsässer B, Goettig P. Mechanisms of Proteolytic Enzymes and Their Inhibition in QM/MM Studies. Int J Mol Sci 2021; 22:3232. [PMID: 33810118 PMCID: PMC8004986 DOI: 10.3390/ijms22063232] [Citation(s) in RCA: 17] [Impact Index Per Article: 5.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/18/2021] [Revised: 02/24/2021] [Accepted: 02/26/2021] [Indexed: 12/11/2022] Open
Abstract
Experimental evidence for enzymatic mechanisms is often scarce, and in many cases inadvertently biased by the employed methods. Thus, apparently contradictory model mechanisms can result in decade long discussions about the correct interpretation of data and the true theory behind it. However, often such opposing views turn out to be special cases of a more comprehensive and superior concept. Molecular dynamics (MD) and the more advanced molecular mechanical and quantum mechanical approach (QM/MM) provide a relatively consistent framework to treat enzymatic mechanisms, in particular, the activity of proteolytic enzymes. In line with this, computational chemistry based on experimental structures came up with studies on all major protease classes in recent years; examples of aspartic, metallo-, cysteine, serine, and threonine protease mechanisms are well founded on corresponding standards. In addition, experimental evidence from enzyme kinetics, structural research, and various other methods supports the described calculated mechanisms. One step beyond is the application of this information to the design of new and powerful inhibitors of disease-related enzymes, such as the HIV protease. In this overview, a few examples demonstrate the high potential of the QM/MM approach for sophisticated pharmaceutical compound design and supporting functions in the analysis of biomolecular structures.
Collapse
Affiliation(s)
| | - Peter Goettig
- Structural Biology Group, Department of Biosciences, University of Salzburg, Billrothstrasse 11, 5020 Salzburg, Austria;
| |
Collapse
|
2
|
Xia J, Flynn W, Gallicchio E, Uplinger K, Armstrong JD, Forli S, Olson AJ, Levy RM. Massive-Scale Binding Free Energy Simulations of HIV Integrase Complexes Using Asynchronous Replica Exchange Framework Implemented on the IBM WCG Distributed Network. J Chem Inf Model 2019; 59:1382-1397. [PMID: 30758197 PMCID: PMC6496938 DOI: 10.1021/acs.jcim.8b00817] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
To perform massive-scale replica exchange molecular dynamics (REMD) simulations for calculating binding free energies of protein-ligand complexes, we implemented the asynchronous replica exchange (AsyncRE) framework of the binding energy distribution analysis method (BEDAM) in implicit solvent on the IBM World Community Grid (WCG) and optimized the simulation parameters to reduce the overhead and improve the prediction power of the WCG AsyncRE simulations. We also performed the first massive-scale binding free energy calculations using the WCG distributed computing grid and 301 ligands from the SAMPL4 challenge for large-scale binding free energy predictions of HIV-1 integrase complexes. In total there are ∼10000 simulated complexes, ∼1 million replicas, and ∼2000 μs of aggregated MD simulations. Running AsyncRE MD simulations on the WCG requires accepting a trade-off between the number of replicas that can be run (breadth) and the number of full RE cycles that can be completed per replica (depth). As compared with synchronous Replica Exchange (SyncRE) running on tightly coupled clusters like XSEDE, on the WCG many more replicas can be launched simultaneously on heterogeneous distributed hardware, but each full RE cycle requires more overhead. We compared the WCG results with that from AutoDock and more advanced RE simulations including the use of flattening potentials to accelerate sampling of selected degrees of freedom of ligands and/or receptors related to slow dynamics due to high energy barriers. We propose a suitable strategy of RE simulations to refine high throughput docking results which can be matched to corresponding computing resources: from HPC clusters, to small or medium-size distributed campus grids, and finally to massive-scale computing networks including millions of CPUs like the resources available on the WCG.
Collapse
Affiliation(s)
- Junchao Xia
- Center for Biophysics and Computational Biology and Department of Physics , Temple University , Philadelphia , Pennsylvania 19122 , United States
| | - William Flynn
- Center for Biophysics and Computational Biology and Department of Chemistry , Temple University , Philadelphia , Pennsylvania 19122 , United States
| | - Emilio Gallicchio
- Department of Chemistry , CUNY Brooklyn College , Brooklyn , New York 11210 , United States
| | - Keith Uplinger
- IBM WCG Team, 1177 South Belt Line Road , Coppell , Texas 75019 , United States
| | - Jonathan D Armstrong
- IBM WCG Team, 11400 Burnet Road , 0453B129, Austin , Texas 78758 , United States
| | - Stefano Forli
- Department of Integrative Structural and Computational Biology , The Scripps Research Institute , La Jolla , California 92037-1000 , United States
| | - Arthur J Olson
- Department of Integrative Structural and Computational Biology , The Scripps Research Institute , La Jolla , California 92037-1000 , United States
| | - Ronald M Levy
- Center for Biophysics and Computational Biology and Department of Chemistry , Temple University , Philadelphia , Pennsylvania 19122 , United States
| |
Collapse
|
3
|
König G, Pickard FC, Huang J, Thiel W, MacKerell AD, Brooks BR, York DM. A Comparison of QM/MM Simulations with and without the Drude Oscillator Model Based on Hydration Free Energies of Simple Solutes. Molecules 2018; 23:E2695. [PMID: 30347691 PMCID: PMC6222909 DOI: 10.3390/molecules23102695] [Citation(s) in RCA: 25] [Impact Index Per Article: 4.2] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/08/2018] [Revised: 10/15/2018] [Accepted: 10/16/2018] [Indexed: 12/01/2022] Open
Abstract
Maintaining a proper balance between specific intermolecular interactions and non-specific solvent interactions is of critical importance in molecular simulations, especially when predicting binding affinities or reaction rates in the condensed phase. The most rigorous metric for characterizing solvent affinity are solvation free energies, which correspond to a transfer from the gas phase into solution. Due to the drastic change of the electrostatic environment during this process, it is also a stringent test of polarization response in the model. Here, we employ both the CHARMM fixed charge and polarizable force fields to predict hydration free energies of twelve simple solutes. The resulting classical ensembles are then reweighted to obtain QM/MM hydration free energies using a variety of QM methods, including MP2, Hartree⁻Fock, density functional methods (BLYP, B3LYP, M06-2X) and semi-empirical methods (OM2 and AM1 ). Our simulations test the compatibility of quantum-mechanical methods with molecular-mechanical water models and solute Lennard⁻Jones parameters. In all cases, the resulting QM/MM hydration free energies were inferior to purely classical results, with the QM/MM Drude force field predictions being only marginally better than the QM/MM fixed charge results. In addition, the QM/MM results for different quantum methods are highly divergent, with almost inverted trends for polarizable and fixed charge water models. While this does not necessarily imply deficiencies in the QM models themselves, it underscores the need to develop consistent and balanced QM/MM interactions. Both the QM and the MM component of a QM/MM simulation have to match, in order to avoid artifacts due to biased solute⁻solvent interactions. Finally, we discuss strategies to improve the convergence and efficiency of multi-scale free energy simulations by automatically adapting the molecular-mechanics force field to the target quantum method.
Collapse
Affiliation(s)
- Gerhard König
- Laboratory of Computational Biology, National Heart, Lung, and Blood Institute, National Institutes of Health, Bethesda, MD 20892, USA.
- Laboratory for Biomolecular Simulation Research, Institute for Quantitative Biomedicine, Department of Chemistry and Chemical Biology, Rutgers University, Piscataway, NJ 08854, USA.
- Max-Planck-Institut für Kohlenforschung, 45470 Mülheim an der Ruhr, Germany.
| | - Frank C Pickard
- Laboratory of Computational Biology, National Heart, Lung, and Blood Institute, National Institutes of Health, Bethesda, MD 20892, USA.
| | - Jing Huang
- Laboratory of Computational Biology, National Heart, Lung, and Blood Institute, National Institutes of Health, Bethesda, MD 20892, USA.
- Department of Pharmaceutical Science, School of Pharmacy, University of Maryland, 20 Penn Street, Baltimore, MD 21201, USA.
- School of Life Sciences, Westlake University, 18 Shilongshan Street, Hangzhou 310024, China.
| | - Walter Thiel
- Max-Planck-Institut für Kohlenforschung, 45470 Mülheim an der Ruhr, Germany.
| | - Alexander D MacKerell
- Department of Pharmaceutical Science, School of Pharmacy, University of Maryland, 20 Penn Street, Baltimore, MD 21201, USA.
| | - Bernard R Brooks
- Laboratory of Computational Biology, National Heart, Lung, and Blood Institute, National Institutes of Health, Bethesda, MD 20892, USA.
| | - Darrin M York
- Laboratory for Biomolecular Simulation Research, Institute for Quantitative Biomedicine, Department of Chemistry and Chemical Biology, Rutgers University, Piscataway, NJ 08854, USA.
| |
Collapse
|
4
|
König G, Brooks BR, Thiel W, York DM. On the convergence of multi-scale free energy simulations. MOLECULAR SIMULATION 2018; 44:1062-1081. [PMID: 30581251 DOI: 10.1080/08927022.2018.1475741] [Citation(s) in RCA: 40] [Impact Index Per Article: 6.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 10/16/2022]
Abstract
In this work we employ simple model systems to evaluate the relative performance of two of the most important free energy methods: The Zwanzig equation (also known as "Free energy perturbation") and Bennett's acceptance ratio method (BAR). Although our examples should be transferable to other kinds of free energy simulations, we focus on applications of multi-scale free energy simulations. Such calculations are especially complex, since they connect two different levels of theory with very different requirements in terms of speed, accuracy, sampling and parallelizability. We try to reconcile all those different factors by developing some simple criteria to guide the early stages of the development of a free energy protocol. This is accomplished by quantifying how many λ intermediate steps and how many potential energy evaluations are necessary in order to reach a certain level of convergence.
Collapse
Affiliation(s)
- Gerhard König
- Max-Planck-Institut für Kohlenforschung, 45470 Mülheim an der Ruhr, Germany, EU.,Laboratory for Biomolecular Simulation Research, Center for Integrative Proteomics Research, and Department of Chemistry and Chemical Biology, Rutgers University, Piscataway, New Jersey 08854, USA.,Laboratory of Computational Biology, National Heart Lung and Blood Institute, National Institutes of Health, Bethesda, Maryland 20892, USA
| | - Bernard R Brooks
- Laboratory of Computational Biology, National Heart Lung and Blood Institute, National Institutes of Health, Bethesda, Maryland 20892, USA
| | - Walter Thiel
- Max-Planck-Institut für Kohlenforschung, 45470 Mülheim an der Ruhr, Germany, EU
| | - 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, USA
| |
Collapse
|
5
|
König G, Pickard FC, Huang J, Simmonett AC, Tofoleanu F, Lee J, Dral PO, Prasad S, Jones M, Shao Y, Thiel W, Brooks BR. Calculating distribution coefficients based on multi-scale free energy simulations: an evaluation of MM and QM/MM explicit solvent simulations of water-cyclohexane transfer in the SAMPL5 challenge. J Comput Aided Mol Des 2016; 30:989-1006. [PMID: 27577746 DOI: 10.1007/s10822-016-9936-x] [Citation(s) in RCA: 22] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/18/2016] [Accepted: 08/09/2016] [Indexed: 11/29/2022]
Abstract
One of the central aspects of biomolecular recognition is the hydrophobic effect, which is experimentally evaluated by measuring the distribution coefficients of compounds between polar and apolar phases. We use our predictions of the distribution coefficients between water and cyclohexane from the SAMPL5 challenge to estimate the hydrophobicity of different explicit solvent simulation techniques. Based on molecular dynamics trajectories with the CHARMM General Force Field, we compare pure molecular mechanics (MM) with quantum-mechanical (QM) calculations based on QM/MM schemes that treat the solvent at the MM level. We perform QM/MM with both density functional theory (BLYP) and semi-empirical methods (OM1, OM2, OM3, PM3). The calculations also serve to test the sensitivity of partition coefficients to solute polarizability as well as the interplay of the quantum-mechanical region with the fixed-charge molecular mechanics environment. Our results indicate that QM/MM with both BLYP and OM2 outperforms pure MM. However, this observation is limited to a subset of cases where convergence of the free energy can be achieved.
Collapse
Affiliation(s)
- Gerhard König
- Laboratory of Computational Biology, National Heart, Lung and Blood Institute, National Institutes of Health, Bethesda, MD, 20892, USA. .,Max-Planck-Institut für Kohlenforschung, 45470, Mülheim an der Ruhr, Germany.
| | - Frank C Pickard
- Laboratory of Computational Biology, National Heart, Lung and Blood Institute, National Institutes of Health, Bethesda, MD, 20892, USA
| | - Jing Huang
- Laboratory of Computational Biology, National Heart, Lung and Blood Institute, National Institutes of Health, Bethesda, MD, 20892, USA
| | - Andrew C Simmonett
- Laboratory of Computational Biology, National Heart, Lung and Blood Institute, National Institutes of Health, Bethesda, MD, 20892, USA
| | - Florentina Tofoleanu
- Laboratory of Computational Biology, National Heart, Lung and Blood Institute, National Institutes of Health, Bethesda, MD, 20892, USA
| | - Juyong Lee
- Laboratory of Computational Biology, National Heart, Lung and Blood Institute, National Institutes of Health, Bethesda, MD, 20892, USA
| | - Pavlo O Dral
- Max-Planck-Institut für Kohlenforschung, 45470, Mülheim an der Ruhr, Germany
| | - Samarjeet Prasad
- Laboratory of Computational Biology, National Heart, Lung and Blood Institute, National Institutes of Health, Bethesda, MD, 20892, USA
| | - Michael Jones
- Laboratory of Computational Biology, National Heart, Lung and Blood Institute, National Institutes of Health, Bethesda, MD, 20892, USA
| | - Yihan Shao
- Laboratory of Computational Biology, National Heart, Lung and Blood Institute, National Institutes of Health, Bethesda, MD, 20892, USA
| | - Walter Thiel
- Max-Planck-Institut für Kohlenforschung, 45470, Mülheim an der Ruhr, Germany
| | - Bernard R Brooks
- Laboratory of Computational Biology, National Heart, Lung and Blood Institute, National Institutes of Health, Bethesda, MD, 20892, USA
| |
Collapse
|
6
|
Lv C, Aitchison EW, Wu D, Zheng L, Cheng X, Yang W. Comparative exploration of hydrogen sulfide and water transmembrane free energy surfaces via orthogonal space tempering free energy sampling. J Comput Chem 2016; 37:567-74. [PMID: 26119423 PMCID: PMC4959446 DOI: 10.1002/jcc.23982] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/28/2015] [Accepted: 05/24/2015] [Indexed: 01/22/2023]
Abstract
Hydrogen sulfide (H2 S), a commonly known toxic gas compound, possesses unique chemical features that allow this small solute molecule to quickly diffuse through cell membranes. Taking advantage of the recent orthogonal space tempering (OST) method, we comparatively mapped the transmembrane free energy landscapes of H2 S and its structural analogue, water (H2 O), seeking to decipher the molecular determinants that govern their drastically different permeabilities. As revealed by our OST sampling results, in contrast to the highly polar water solute, hydrogen sulfide is evidently amphipathic, and thus inside membrane is favorably localized at the interfacial region, that is, the interface between the polar head-group and nonpolar acyl chain regions. Because the membrane binding affinity of H2 S is mainly governed by its small hydrophobic moiety and the barrier height inbetween the interfacial region and the membrane center is largely determined by its moderate polarity, the transmembrane free energy barriers to encounter by this toxic molecule are very small. Moreover when H2 S diffuses from the bulk solution to the membrane center, the above two effects nearly cancel each other, so as to lead to a negligible free energy difference. This study not only explains why H2 S can quickly pass through cell membranes but also provides a practical illustration on how to use the OST free energy sampling method to conveniently analyze complex molecular processes. © 2015 Wiley Periodicals, Inc.
Collapse
Affiliation(s)
- Chao Lv
- Department of Chemistry & Biochemistry, Florida State University, Tallahassee, Florida 32306
| | - Erick W. Aitchison
- Institute of Molecular Biophysics, Florida State University, Tallahassee, Florida 32306
| | - Dongsheng Wu
- Institute of Molecular Biophysics, Florida State University, Tallahassee, Florida 32306
| | - Lianqing Zheng
- Institute of Molecular Biophysics, Florida State University, Tallahassee, Florida 32306
| | - Xiaolin Cheng
- UT-ORNL Center for Molecular Biophysics, Oak Ridge National Laboratory, Oak Ridge, Tennessee, 37830
- Department of Biochemistry, Cellular and Molecular Biology, the University of Tennessee, Knoxville, Tennessee, 37996
| | - Wei Yang
- Department of Chemistry & Biochemistry, Florida State University, Tallahassee, Florida 32306
- Institute of Molecular Biophysics, Florida State University, Tallahassee, Florida 32306
| |
Collapse
|
7
|
Lv C, Li X, Wu D, Zheng L, Yang W. Predictive Sampling of Rare Conformational Events in Aqueous Solution: Designing a Generalized Orthogonal Space Tempering Method. J Chem Theory Comput 2016; 12:41-52. [PMID: 26636477 PMCID: PMC4968881 DOI: 10.1021/acs.jctc.5b00953] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
In aqueous solution, solute conformational transitions are governed by intimate interplays of the fluctuations of solute-solute, solute-water, and water-water interactions. To promote molecular fluctuations to enhance sampling of essential conformational changes, a common strategy is to construct an expanded Hamiltonian through a series of Hamiltonian perturbations and thereby broaden the distribution of certain interactions of focus. Due to a lack of active sampling of configuration response to Hamiltonian transitions, it is challenging for common expanded Hamiltonian methods to robustly explore solvent mediated rare conformational events. The orthogonal space sampling (OSS) scheme, as exemplified by the orthogonal space random walk and orthogonal space tempering methods, provides a general framework for synchronous acceleration of slow configuration responses. To more effectively sample conformational transitions in aqueous solution, in this work, we devised a generalized orthogonal space tempering (gOST) algorithm. Specifically, in the Hamiltonian perturbation part, a solvent-accessible-surface-area-dependent term is introduced to implicitly perturb near-solute water-water fluctuations; more importantly in the orthogonal space response part, the generalized force order parameter is generalized as a two-dimension order parameter set, in which essential solute-solvent and solute-solute components are separately treated. The gOST algorithm is evaluated through a molecular dynamics simulation study on the explicitly solvated deca-alanine (Ala10) peptide. On the basis of a fully automated sampling protocol, the gOST simulation enabled repetitive folding and unfolding of the solvated peptide within a single continuous trajectory and allowed for detailed constructions of Ala10 folding/unfolding free energy surfaces. The gOST result reveals that solvent cooperative fluctuations play a pivotal role in Ala10 folding/unfolding transitions. In addition, our assessment analysis suggests that because essential conformational events are mainly driven by the compensating fluctuations of essential solute-solvent and solute-solute interactions, commonly employed "predictive" sampling methods are unlikely to be effective on this seemingly "simple" system. The gOST development presented in this paper illustrates how to employ the OSS scheme for physics-based sampling method designs.
Collapse
Affiliation(s)
- Chao Lv
- Department of Chemistry and Biochemistry, Florida State University, Tallahassee, FL, 32306
| | - Xubin Li
- Department of Chemistry and Biochemistry, Florida State University, Tallahassee, FL, 32306
| | - Dongsheng Wu
- Institute of Molecular Biophysics, Florida State University, Tallahassee, FL, 32306
| | - Lianqing Zheng
- Institute of Molecular Biophysics, Florida State University, Tallahassee, FL, 32306
| | - Wei Yang
- Department of Chemistry and Biochemistry, Florida State University, Tallahassee, FL, 32306
- Institute of Molecular Biophysics, Florida State University, Tallahassee, FL, 32306
| |
Collapse
|
8
|
Xia J, Flynn WF, Gallicchio E, Zhang BW, He P, Tan Z, Levy RM. Large-scale asynchronous and distributed multidimensional replica exchange molecular simulations and efficiency analysis. J Comput Chem 2015; 36:1772-85. [PMID: 26149645 PMCID: PMC4512903 DOI: 10.1002/jcc.23996] [Citation(s) in RCA: 18] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/20/2015] [Revised: 06/10/2015] [Accepted: 06/11/2015] [Indexed: 01/25/2023]
Abstract
We describe methods to perform replica exchange molecular dynamics (REMD) simulations asynchronously (ASyncRE). The methods are designed to facilitate large scale REMD simulations on grid computing networks consisting of heterogeneous and distributed computing environments as well as on homogeneous high-performance clusters. We have implemented these methods on NSF (National Science Foundation) XSEDE (Extreme Science and Engineering Discovery Environment) clusters and BOINC (Berkeley Open Infrastructure for Network Computing) distributed computing networks at Temple University and Brooklyn College at CUNY (the City University of New York). They are also being implemented on the IBM World Community Grid. To illustrate the methods, we have performed extensive (more than 60 ms in aggregate) simulations for the beta-cyclodextrin-heptanoate host-guest system in the context of one- and two-dimensional ASyncRE, and we used the results to estimate absolute binding free energies using the binding energy distribution analysis method. We propose ways to improve the efficiency of REMD simulations: these include increasing the number of exchanges attempted after a specified molecular dynamics (MD) period up to the fast exchange limit and/or adjusting the MD period to allow sufficient internal relaxation within each thermodynamic state. Although ASyncRE simulations generally require long MD periods (>picoseconds) per replica exchange cycle to minimize the overhead imposed by heterogeneous computing networks, we found that it is possible to reach an efficiency similar to conventional synchronous REMD, by optimizing the combination of the MD period and the number of exchanges attempted per cycle.
Collapse
Affiliation(s)
- Junchao Xia
- Center for Biophysics and Computational Biology, Department of Chemistry and Institute for Computational Molecular Science, Temple University, Philadelphia, PA 19122
| | - William F. Flynn
- Center for Biophysics and Computational Biology, Department of Chemistry and Institute for Computational Molecular Science, Temple University, Philadelphia, PA 19122
- Department of Physics & Astronomy, Rutgers University, Piscataway, NJ 08854
| | | | - Bin W. Zhang
- Center for Biophysics and Computational Biology, Department of Chemistry and Institute for Computational Molecular Science, Temple University, Philadelphia, PA 19122
| | - Peng He
- Center for Biophysics and Computational Biology, Department of Chemistry and Institute for Computational Molecular Science, Temple University, Philadelphia, PA 19122
| | - Zhiqiang Tan
- Department of Statistics, Rutgers University, Piscataway, NJ 08854
| | - Ronald M. Levy
- Center for Biophysics and Computational Biology, Department of Chemistry and Institute for Computational Molecular Science, Temple University, Philadelphia, PA 19122
| |
Collapse
|
9
|
Li X, Lv C, Corbett KM, Zheng L, Wu D, Yang W. Free energy landscape of a minimalist salt bridge model. Protein Sci 2015; 25:270-6. [PMID: 26300526 DOI: 10.1002/pro.2789] [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: 07/18/2015] [Revised: 08/19/2015] [Accepted: 08/20/2015] [Indexed: 11/06/2022]
Abstract
Salt bridges are essential to protein stability and dynamics. Despite the importance, there has been scarce of detailed discussion on how salt bridge partners interact with each other in distinct solvent exposed environments. In this study, employing a recent generalized orthogonal space tempering (gOST) method, we enabled efficient molecular dynamics simulation of repetitive breaking and reforming of salt bridge structures within a minimalist salt-bridge model, the Asp-Arg dipeptide and thereby were able to map its detailed free energy landscape in aqueous solution. Free energy surface analysis shows that although individually-solvated states are more favorable, salt-bridge states still occupy a noticeable portion of the overall population. Notably, the competing forces, e.g. intercharge attractions that drive the formation of salt bridges and solvation forces that pull the charged groups away from each other, are energetically comparable. As the result, the salt bridge stability is highly tunable by local environments; for instance when local water molecules are perturbed to interact more strongly with each other, the population of the salt-bridge states is likely to increase. Our results reveal the critical role of local solvent structures in modulating salt-bridge partner interactions and imply the importance of water fluctuations on conformational dynamics that involves solvent accessible salt bridge formations.
Collapse
Affiliation(s)
- Xubin Li
- Department of Chemistry and Biochemistry, Florida State University, Tallahassee, Florida, 32306
| | - Chao Lv
- Department of Chemistry and Biochemistry, Florida State University, Tallahassee, Florida, 32306
| | - Karen M Corbett
- Institute of Molecular Biophysics, Florida State University, Tallahassee, Florida, 32306
| | - Lianqing Zheng
- Institute of Molecular Biophysics, Florida State University, Tallahassee, Florida, 32306
| | - Dongsheng Wu
- Institute of Molecular Biophysics, Florida State University, Tallahassee, Florida, 32306
| | - Wei Yang
- Department of Chemistry and Biochemistry, Florida State University, Tallahassee, Florida, 32306.,Institute of Molecular Biophysics, Florida State University, Tallahassee, Florida, 32306
| |
Collapse
|
10
|
Khrenova M, Savitsky AP, Topol IA, Nemukhin AV. Exploration of the zinc finger motif in controlling activity of matrix metalloproteinases. J Phys Chem B 2014; 118:13505-12. [PMID: 25375834 PMCID: PMC4254000 DOI: 10.1021/jp5088702] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/02/2014] [Revised: 11/05/2014] [Indexed: 12/19/2022]
Abstract
Discovering ways to control the activity of matrix metalloproteinases (MMPs), zinc-dependent enzymes capable of degrading extracellular matrix proteins, is an important field of cancer research. We report here a novel strategy for assembling MMP inhibitors on the basis of oligopeptide ligands by exploring the pattern known as the zinc finger motif. Advanced molecular modeling tools were used to characterize the structural binding motifs of experimentally tested MMP inhibitors, as well as those of newly proposed peptidomimetics, in their zinc-containing active sites. The results of simulations based on the quantum mechanics/molecular mechanics (QM/MM) approach and Car-Parrinello molecular dynamics with QM/MM potentials demonstrate that, upon binding of Regasepin1, a known MMP-9 inhibitor, the Zn(2+)(His3) structural element is rearranged to the Zn(2+)(Cys2His2) zinc finger motif, in which two Cys residues are borrowed from the ligand. Following consideration of the crystal structure of MMP-2 with its inhibitor, the oligopeptide APP-IP, we proposed a new peptidomimetic with two replacements in the substrate, Tyr3Cys and Asp6Cys. Simulations show that this peptide variant blocks an enzyme active site by the Zn(2+)(Cys2His2) zinc finger construct. Similarly, a natural substrate of MMP-2, Ace-Gln-Gly ∼ Ile-Ala-Gly-Nme, can be converted to an inhibiting compound by two replacements, Ile by Cys and Gly by the d isomer of Cys, favoring formation of the zinc finger motif.
Collapse
Affiliation(s)
- Maria
G. Khrenova
- A.N. Bach Institute
of Biochemistry of the Russian Academy of Science, Leninsky Prospect, 33, Moscow 119071, Russian Federation
- Chemistry
Department, M.V. Lomonosov Moscow State
University, Leninskie
Gory 1/3, Moscow, 119991, Russian Federation
| | - Alexander P. Savitsky
- A.N. Bach Institute
of Biochemistry of the Russian Academy of Science, Leninsky Prospect, 33, Moscow 119071, Russian Federation
| | - Igor A. Topol
- Advanced
Biomedical Computing Center, Information Systems Program, Leidos Biomedical
Research Inc., Frederick National Laboratory
for Cancer Research, Frederick, Maryland 21702, United States
| | - Alexander V. Nemukhin
- Chemistry
Department, M.V. Lomonosov Moscow State
University, Leninskie
Gory 1/3, Moscow, 119991, Russian Federation
- N.M.
Emanuel Institute of Biochemical Physics, Russian Academy of Sciences, Kosygina 4, Moscow, 119334, Russian Federation
| |
Collapse
|
11
|
Hudson PS, White JK, Kearns FL, Hodoscek M, Boresch S, Lee Woodcock H. Efficiently computing pathway free energies: New approaches based on chain-of-replica and Non-Boltzmann Bennett reweighting schemes. Biochim Biophys Acta Gen Subj 2014; 1850:944-953. [PMID: 25239198 DOI: 10.1016/j.bbagen.2014.09.016] [Citation(s) in RCA: 29] [Impact Index Per Article: 2.9] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/14/2014] [Revised: 09/09/2014] [Accepted: 09/10/2014] [Indexed: 11/25/2022]
Abstract
BACKGROUND Accurately modeling condensed phase processes is one of computation's most difficult challenges. Include the possibility that conformational dynamics may be coupled to chemical reactions, where multiscale (i.e., QM/MM) methods are needed, and this task becomes even more daunting. METHODS Free energy simulations (i.e., molecular dynamics), multiscale modeling, and reweighting schemes. RESULTS Herein, we present two new approaches for mitigating the aforementioned challenges. The first is a new chain-of-replica method (off-path simulations, OPS) for computing potentials of mean force (PMFs) along an easily defined reaction coordinate. This development is coupled with a new distributed, highly-parallel replica framework (REPDstr) within the CHARMM package. Validation of these new schemes is carried out on two processes that undergo conformational changes. First is the simple torsional rotation of butane, while a much more challenging glycosidic rotation (in vacuo and solvated) is the second. Additionally, a new approach that greatly improves (i.e., possibly an order of magnitude) the efficiency of computing QM/MM PMFs is introduced and compared to standard schemes. Our efforts are grounded in the recently developed method for efficiently computing QM-based free energies (i.e., QM-Non-Boltzmann Bennett, QM-NBB). Again, we validate this new technique by computing the QM/MM PMF of butane's torsional rotation. CONCLUSIONS The OPS-REPDstr method is a promising new approach that overcomes many limitations of standard pathway simulations in CHARMM. The combination of QM-NBB with pathway techniques is very promising as it offers significant advantages over current procedures. GENERAL SIGNIFICANCE Efficiently computing potentials of mean force is a major, unresolved, area of interest. This article is part of a Special Issue entitled Recent developments of molecular dynamics.
Collapse
Affiliation(s)
- Phillip S Hudson
- Department of Chemistry, University of South Florida, 4202 E. Fowler Ave., CHE205, Tampa, FL 33620-5250, USA
| | - Justin K White
- Department of Chemistry, University of South Florida, 4202 E. Fowler Ave., CHE205, Tampa, FL 33620-5250, USA
| | - Fiona L Kearns
- Department of Chemistry, University of South Florida, 4202 E. Fowler Ave., CHE205, Tampa, FL 33620-5250, USA
| | - Milan Hodoscek
- Center for Molecular Modeling, National Institute of Chemistry, Hajdrihova 19, SI-1000 Ljubljana, Slovenia
| | - Stefan Boresch
- Department of Computational Biological Chemistry, Faculty of Chemistry, University of Vienna, Währingerstraße 17, A-1090 Vienna, Austria
| | - H Lee Woodcock
- Department of Chemistry, University of South Florida, 4202 E. Fowler Ave., CHE205, Tampa, FL 33620-5250, USA.
| |
Collapse
|
12
|
König G, Hudson PS, Boresch S, Woodcock HL. Multiscale Free Energy Simulations: An Efficient Method for Connecting Classical MD Simulations to QM or QM/MM Free Energies Using Non-Boltzmann Bennett Reweighting Schemes. J Chem Theory Comput 2014; 10:1406-1419. [PMID: 24803863 PMCID: PMC3985817 DOI: 10.1021/ct401118k] [Citation(s) in RCA: 107] [Impact Index Per Article: 10.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/31/2013] [Indexed: 11/28/2022]
Abstract
![]()
The reliability of free energy simulations
(FES) is limited by
two factors: (a) the need for correct sampling and (b) the accuracy
of the computational method employed. Classical methods (e.g., force
fields) are typically used for FES and present a myriad of challenges,
with parametrization being a principle one. On the other hand, parameter-free
quantum mechanical (QM) methods tend to be too computationally expensive
for adequate sampling. One widely used approach is a combination of
methods, where the free energy difference between the two end states
is computed by, e.g., molecular mechanics (MM), and the end states
are corrected by more accurate methods, such as QM or hybrid QM/MM
techniques. Here we report two new approaches that significantly improve
the aforementioned scheme; with a focus on how to compute corrections
between, e.g., the MM and the more accurate QM calculations. First,
a molecular dynamics trajectory that properly samples relevant conformational
degrees of freedom is generated. Next, potential energies of each
trajectory frame are generated with a QM or QM/MM Hamiltonian. Free
energy differences are then calculated based on the QM or QM/MM energies
using either a non-Boltzmann Bennett approach (QM-NBB) or non-Boltzmann
free energy perturbation (NB-FEP). Both approaches are applied to
calculate relative and absolute solvation free energies in explicit
and implicit solvent environments. Solvation free energy differences
(relative and absolute) between ethane and methanol in explicit solvent
are used as the initial test case for QM-NBB. Next, implicit solvent
methods are employed in conjunction with both QM-NBB and NB-FEP to
compute absolute solvation free energies for 21 compounds. These compounds
range from small molecules such as ethane and methanol to fairly large,
flexible solutes, such as triacetyl glycerol. Several technical aspects
were investigated. Ultimately some best practices are suggested for
improving methods that seek to connect MM to QM (or QM/MM) levels
of theory in FES.
Collapse
Affiliation(s)
- Gerhard König
- Laboratory of Computational Biology, National Heart Lung and Blood Institute, National Institutes of Health , Bethesda, Maryland 20892, United States
| | - Phillip S Hudson
- Department of Chemistry, University of South Florida , 4202 E. Fowler Avenue, CHE205, Tampa, Florida 33620-5250, United States
| | - Stefan Boresch
- Department of Computational Biological Chemistry, Faculty of Chemistry, University of Vienna , Währingerstraße 17, A-1090 Vienna, Austria
| | - H Lee Woodcock
- Department of Chemistry, University of South Florida , 4202 E. Fowler Avenue, CHE205, Tampa, Florida 33620-5250, United States
| |
Collapse
|
13
|
Lee TS, Radak BK, Huang M, Wong KY, York DM. Roadmaps through free energy landscapes calculated using the multi-dimensional vFEP approach. J Chem Theory Comput 2013; 10:24-34. [PMID: 24505217 DOI: 10.1021/ct400691f] [Citation(s) in RCA: 55] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/27/2022]
Abstract
The variational free energy profile (vFEP) method is extended to two dimensions and tested with molecular simulation applications. The proposed 2D-vFEP approach effectively addresses the two major obstacles to constructing free energy profiles from simulation data using traditional methods: the need for overlap in the re-weighting procedure and the problem of data representation. This is especially evident as these problems are shown to be more severe in two dimensions. The vFEP method is demonstrated to be highly robust and able to provide stable, analytic free energy profiles with only a paucity of sampled data. The analytic profiles can be analyzed with conventional search methods to easily identify stationary points (e.g. minima and first-order saddle points) as well as the pathways that connect these points. These "roadmaps" through the free energy surface are useful not only as a post-processing tool to characterize mechanisms, but can also serve as a basis from which to direct more focused "on-the-fly" sampling or adaptive force biasing. Test cases demonstrate that 2D-vFEP outperforms other methods in terms of the amount and sparsity of the data needed to construct stable, converged analytic free energy profiles. In a classic test case, the two dimensional free energy profile of the backbone torsion angles of alanine dipeptide, 2D-vFEP needs less than 1% of the original data set to reach a sampling accuracy of 0.5 kcal/mol in free energy shifts between windows. A new software tool for performing one and two dimensional vFEP calculations is herein described and made publicly available.
Collapse
Affiliation(s)
- Tai-Sung Lee
- Center for Integrative Proteomics Research and Department of Chemistry and Chemical Biology, Rutgers University, Piscataway, NJ 08854, USA, Department of Chemistry, University of Minnesota, Minneapolis, MN 55455, USA, Scientific Computation Program, University of Minnesota, Minneapolis, MN 55455, USA, and Department of Physics, High Performance Cluster Computing Centre, and Institute of Computational and Theoretical Studies, Hong Kong Baptist University, Kowloon Tong, Hong Kong
| | - Brian K Radak
- Center for Integrative Proteomics Research and Department of Chemistry and Chemical Biology, Rutgers University, Piscataway, NJ 08854, USA, Department of Chemistry, University of Minnesota, Minneapolis, MN 55455, USA, Scientific Computation Program, University of Minnesota, Minneapolis, MN 55455, USA, and Department of Physics, High Performance Cluster Computing Centre, and Institute of Computational and Theoretical Studies, Hong Kong Baptist University, Kowloon Tong, Hong Kong
| | - Ming Huang
- Center for Integrative Proteomics Research and Department of Chemistry and Chemical Biology, Rutgers University, Piscataway, NJ 08854, USA, Department of Chemistry, University of Minnesota, Minneapolis, MN 55455, USA, Scientific Computation Program, University of Minnesota, Minneapolis, MN 55455, USA, and Department of Physics, High Performance Cluster Computing Centre, and Institute of Computational and Theoretical Studies, Hong Kong Baptist University, Kowloon Tong, Hong Kong
| | - Kin-Yiu Wong
- Center for Integrative Proteomics Research and Department of Chemistry and Chemical Biology, Rutgers University, Piscataway, NJ 08854, USA, Department of Chemistry, University of Minnesota, Minneapolis, MN 55455, USA, Scientific Computation Program, University of Minnesota, Minneapolis, MN 55455, USA, and Department of Physics, High Performance Cluster Computing Centre, and Institute of Computational and Theoretical Studies, Hong Kong Baptist University, Kowloon Tong, Hong Kong
| | - Darrin M York
- Center for Integrative Proteomics Research and Department of Chemistry and Chemical Biology, Rutgers University, Piscataway, NJ 08854, USA, Department of Chemistry, University of Minnesota, Minneapolis, MN 55455, USA, Scientific Computation Program, University of Minnesota, Minneapolis, MN 55455, USA, and Department of Physics, High Performance Cluster Computing Centre, and Institute of Computational and Theoretical Studies, Hong Kong Baptist University, Kowloon Tong, Hong Kong
| |
Collapse
|
14
|
Lin Z, van Gunsteren WF. Enhanced conformational sampling using enveloping distribution sampling. J Chem Phys 2013; 139:144105. [DOI: 10.1063/1.4824391] [Citation(s) in RCA: 10] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
15
|
Cao L, Lv C, Yang W. Hidden Conformation Events in DNA Base Extrusions: A Generalized Ensemble Path Optimization and Equilibrium Simulation Study. J Chem Theory Comput 2013; 9:10.1021/ct400198q. [PMID: 24250279 PMCID: PMC3829643 DOI: 10.1021/ct400198q] [Citation(s) in RCA: 16] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/01/2023]
Abstract
DNA base extrusion is a crucial component of many biomolecular processes. Elucidating how bases are selectively extruded from the interiors of double-strand DNAs is pivotal to accurately understanding and efficiently sampling this general type of conformational transitions. In this work, the on-the-path random walk (OTPRW) method, which is the first generalized ensemble sampling scheme designed for finite-temperature-string path optimizations, was improved and applied to obtain the minimum free energy path (MFEP) and the free energy profile of a classical B-DNA major-groove base extrusion pathway. Along the MFEP, an intermediate state and the corresponding transition state were located and characterized. The MFEP result suggests that a base-plane-elongation event rather than the commonly focused base-flipping event is dominant in the transition state formation portion of the pathway; and the energetic penalty at the transition state is mainly introduced by the stretching of the Watson-Crick base pair. Moreover to facilitate the essential base-plane-elongation dynamics, the surrounding environment of the flipped base needs to be intimately involved. Further taking the advantage of the extended-dynamics nature of the OTPRW Hamiltonian, an equilibrium generalized ensemble simulation was performed along the optimized path; and based on the collected samples, several base-flipping (opening) angle collective variables were evaluated. In consistence with the MFEP result, the collective variable analysis result reveals that none of these commonly employed flipping (opening) angles alone can adequately represent the base extrusion pathway, especially in the pre-transition-state portion. As further revealed by the collective variable analysis, the base-pairing partner of the extrusion target undergoes a series of in-plane rotations to facilitate the base-plane-elongation dynamics. A base-plane rotation angle is identified to be a possible reaction coordinate to represent these in-plane rotations. Notably, these in-plane rotation motions may play a pivotal role in determining the base extrusion selectivity.
Collapse
Affiliation(s)
- Liaoran Cao
- Department of Chemistry and Biochemistry, Florida State University, Tallahassee, FL, 32306
| | - Chao Lv
- Department of Chemistry and Biochemistry, Florida State University, Tallahassee, FL, 32306
| | - Wei Yang
- Department of Chemistry and Biochemistry, Florida State University, Tallahassee, FL, 32306
- Institute of Molecular Biophysics, Florida State University, Tallahassee, FL, 32306
| |
Collapse
|
16
|
Lv C, Zheng L, Yang W. Generalized essential energy space random walks to more effectively accelerate solute sampling in aqueous environment. J Chem Phys 2012; 136:044103. [PMID: 22299857 DOI: 10.1063/1.3678220] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Molecular dynamics sampling can be enhanced via the promoting of potential energy fluctuations, for instance, based on a Hamiltonian modified with the addition of a potential-energy-dependent biasing term. To overcome the diffusion sampling issue, which reveals the fact that enlargement of event-irrelevant energy fluctuations may abolish sampling efficiency, the essential energy space random walk (EESRW) approach was proposed earlier. To more effectively accelerate the sampling of solute conformations in aqueous environment, in the current work, we generalized the EESRW method to a two-dimension-EESRW (2D-EESRW) strategy. Specifically, the essential internal energy component of a focused region and the essential interaction energy component between the focused region and the environmental region are employed to define the two-dimensional essential energy space. This proposal is motivated by the general observation that in different conformational events, the two essential energy components have distinctive interplays. Model studies on the alanine dipeptide and the aspartate-arginine peptide demonstrate sampling improvement over the original one-dimension-EESRW strategy; with the same biasing level, the present generalization allows more effective acceleration of the sampling of conformational transitions in aqueous solution. The 2D-EESRW generalization is readily extended to higher dimension schemes and employed in more advanced enhanced-sampling schemes, such as the recent orthogonal space random walk method.
Collapse
Affiliation(s)
- Chao Lv
- Department of Chemistry and Biochemistry, Florida State University, Tallahassee, Florida 32306, USA
| | | | | |
Collapse
|
17
|
Wu X, Damjanovic A, Brooks BR. Efficient and Unbiased Sampling of Biomolecular Systems in the Canonical Ensemble: A Review of Self-Guided Langevin Dynamics. ADVANCES IN CHEMICAL PHYSICS 2012; 150:255-326. [PMID: 23913991 PMCID: PMC3731171 DOI: 10.1002/9781118197714.ch6] [Citation(s) in RCA: 28] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 04/08/2023]
Abstract
This review provides a comprehensive description of the self-guided Langevin dynamics (SGLD) and the self-guided molecular dynamics (SGMD) methods and their applications. Example systems are included to provide guidance on optimal application of these methods in simulation studies. SGMD/SGLD has enhanced ability to overcome energy barriers and accelerate rare events to affordable time scales. It has been demonstrated that with moderate parameters, SGLD can routinely cross energy barriers of 20 kT at a rate that molecular dynamics (MD) or Langevin dynamics (LD) crosses 10 kT barriers. The core of these methods is the use of local averages of forces and momenta in a direct manner that can preserve the canonical ensemble. The use of such local averages results in methods where low frequency motion "borrows" energy from high frequency degrees of freedom when a barrier is approached and then returns that excess energy after a barrier is crossed. This self-guiding effect also results in an accelerated diffusion to enhance conformational sampling efficiency. The resulting ensemble with SGLD deviates in a small way from the canonical ensemble, and that deviation can be corrected with either an on-the-fly or a post processing reweighting procedure that provides an excellent canonical ensemble for systems with a limited number of accelerated degrees of freedom. Since reweighting procedures are generally not size extensive, a newer method, SGLDfp, uses local averages of both momenta and forces to preserve the ensemble without reweighting. The SGLDfp approach is size extensive and can be used to accelerate low frequency motion in large systems, or in systems with explicit solvent where solvent diffusion is also to be enhanced. Since these methods are direct and straightforward, they can be used in conjunction with many other sampling methods or free energy methods by simply replacing the integration of degrees of freedom that are normally sampled by MD or LD.
Collapse
Affiliation(s)
- Xiongwu Wu
- Laboratory of Computational Biology, National Heart, Lung, and Blood Institute (NHLBI), National Institutes of Health(NIH), 5635 Fishers Lane, Room T900, Bethesda, MD 20892-9314
| | - Ana Damjanovic
- Laboratory of Computational Biology, National Heart, Lung, and Blood Institute (NHLBI), National Institutes of Health(NIH), 5635 Fishers Lane, Room T900, Bethesda, MD 20892-9314
- Johns Hopkins University, Department of Biophysics, 3400 N. Charles Street, Baltimore, MD 21218
| | - Bernard R. Brooks
- Laboratory of Computational Biology, National Heart, Lung, and Blood Institute (NHLBI), National Institutes of Health(NIH), 5635 Fishers Lane, Room T900, Bethesda, MD 20892-9314
| |
Collapse
|
18
|
Gallicchio E. Role of Ligand Reorganization and Conformational Restraints on the Binding Free Energies of DAPY Non-Nucleoside Inhibitors to HIV Reverse Transcriptase. ACTA ACUST UNITED AC 2012; 2:7-22. [PMID: 22708073 DOI: 10.4236/cmb.2012.21002] [Citation(s) in RCA: 16] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/01/2023]
Abstract
The results of computer simulations of the binding of etravirine (TMC125) and rilpivirine (TMC278) to HIV reverse transcriptase are reported. It is confirmed that consistent binding free energy estimates are obtained with or without the application of torsional restraints when the free energies of imposing the restraints are taken into account. The restraints have a smaller influence on the thermodynamics and apparent kinetics of binding of TMC125 compared to the more flexible TMC278 inhibitor. The concept of the reorganization free energy of binding is useful to understand and categorize these effects. Contrary to expectations, the use of conformational restraints did not consistently enhance convergence of binding free energy estimates due to suppression of binding/unbinding pathways and due to the influence of rotational degrees of freedom not directly controlled by the restraints. Physical insights concerning the thermodynamic driving forces for binding and the role of "jiggling" and "wiggling" motion of the ligands are discussed. Based on these insights we conclude that an ideal inhibitor, if chemically realizable, would possess the electrostatic charge distribution of TMC125, so as to form strong interactions with the receptor, and the larger and more flexible substituents of TMC278, so as to minimize reorganization free energy penalties and the effects of resistance mutations, suitably modified, as in TMC125, so as to disfavor the formation of non-binding competent extended conformations when free in solution.
Collapse
Affiliation(s)
- Emilio Gallicchio
- BioMaPS Institute for Quantitative Biology and Department of Chemistry and Chemical Biology, Rutgers the State University of New Jersey, Piseatawav, NJ 08854
| |
Collapse
|
19
|
Lapelosa M, Gallicchio E, Levy RM. Conformational Transitions and Convergence of Absolute Binding Free Energy Calculations. J Chem Theory Comput 2011; 8:47-60. [PMID: 22368530 DOI: 10.1021/ct200684b] [Citation(s) in RCA: 55] [Impact Index Per Article: 4.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/03/2023]
Abstract
The Binding Energy Distribution Analysis Method (BEDAM) is employed to compute the standard binding free energies of a series of ligands to a FK506 binding protein (FKBP12) with implicit solvation. Binding free energy estimates are in reasonably good agreement with experimental affinities. The conformations of the complexes identified by the simulations are in good agreement with crystallographic data, which was not used to restrain ligand orientations. The BEDAM method is based on λ -hopping Hamiltonian parallel Replica Exchange (HREM) molecular dynamics conformational sampling, the OPLS-AA/AGBNP2 effective potential, and multi-state free energy estimators (MBAR). Achieving converged and accurate results depends on all of these elements of the calculation. Convergence of the binding free energy is tied to the level of convergence of binding energy distributions at critical intermediate states where bound and unbound states are at equilibrium, and where the rate of binding/unbinding conformational transitions is maximal. This finding mirrors similar observations in the context of order/disorder transitions as for example in protein folding. Insights concerning the physical mechanism of ligand binding and unbinding are obtained. Convergence for the largest FK506 ligand is achieved only after imposing strict conformational restraints, which however require accurate prior structural knowledge of the structure of the complex. The analytical AGBNP2 model is found to underestimate the magnitude of the hydrophobic driving force towards binding in these systems characterized by loosely packed protein-ligand binding interfaces. Rescoring of the binding energies using a numerical surface area model corrects this deficiency. This study illustrates the complex interplay between energy models, exploration of conformational space, and free energy estimators needed to obtain robust estimates from binding free energy calculations.
Collapse
Affiliation(s)
- Mauro Lapelosa
- BioMaPS Institute for Quantitative Biology and Department of Chemistry and Chemical Biology, Rutgers the State University of New Jersey, Piscataway, NJ 08854
| | | | | |
Collapse
|
20
|
Wu P, Hu X, Yang W. λ-Meta Dynamics Approach To Compute Absolute Solvation Free Energy. J Phys Chem Lett 2011; 2:2099-2103. [PMID: 23678385 PMCID: PMC3652470 DOI: 10.1021/jz200808x] [Citation(s) in RCA: 18] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/02/2023]
Abstract
We present a new approach to combine λ dynamics with meta-dynamics (named λ-meta dynamics) to compute free energy surface with respect to λ. Particularly, the λ-meta dynamics method extends meta-dynamics to a single virtual variable λ, i.e., the coupling parameter between solute and solvent, to compute absolute solvation free energy as an exemplary application. We demonstrate that λ-meta dynamics simulations can recover the accurate potential of mean force surface with respect to λ compared to the benchmark results from traditional λ-dynamics with umbrella sampling. The solvation free energy results for five small organic molecules from λ-meta dynamics simulations using the same filling scheme show that the statistical errors are within ±0.5 kcal/mol. The new λ-meta dynamics method is general and other variables such as order parameters to describe conformational changes can be easily combined with λ-meta dynamics. This should allow for efficient samplings on high-dimension free energy landscapes.
Collapse
|