1
|
Gilson MK, Stewart LE, Potter MJ, Webb SP. Rapid, Accurate, Ranking of Protein-Ligand Binding Affinities with VM2, the Second-Generation Mining Minima Method. J Chem Theory Comput 2024; 20:6328-6340. [PMID: 38989926 DOI: 10.1021/acs.jctc.4c00407] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/12/2024]
Abstract
The structure-based technologies most widely used to rank the affinities of candidate small molecule drugs for proteins range from faster but less reliable docking methods to slower but more accurate explicit solvent free energy methods. In recent years, we have advanced another technology, which is called mining minima because it "mines" out the main contributions to the chemical potentials of the free and bound molecular species by identifying and characterizing their main local energy minima. The present study provides systematic benchmarks of the accuracy and computational speed of mining minima, as implemented in the VeraChem Mining Minima Generation 2 (VM2) code, across two well-regarded protein-ligand benchmark data sets, for which there are already benchmark data for docking, free energy, and other computational methods. A core result is that VM2's accuracy approaches that of explicit solvent free energy methods at a far lower computational cost. In finer-grained analyses, we also examine the influence of various run settings, such as the treatment of crystallographic water molecules, on the accuracy, and define the costs in time and dollars of representative runs on Amazon Web Services (AWS) compute instances with various CPU and GPU combinations. We also use the benchmark data to determine the importance of VM2's correction from generalized Born to finite-difference Poisson-Boltzmann results for each energy well and find that this correction affords a remarkably consistent improvement in accuracy at a modest computational cost. The present results establish VM2 as a distinctive technology for early-stage drug discovery, which provides a strong combination of efficiency and predictivity.
Collapse
Affiliation(s)
- Michael K Gilson
- VeraChem LLC, 12850 Middlebrook Rd, Ste 205, Germantown, Maryland 20874, United States
- Skaggs School of Pharmacy and Pharmaceutical Sciences, University of California San Diego, 9255 Pharmacy Lane, La Jolla, California 92093, United States
| | - Lawrence E Stewart
- VeraChem LLC, 12850 Middlebrook Rd, Ste 205, Germantown, Maryland 20874, United States
| | - Michael J Potter
- VeraChem LLC, 12850 Middlebrook Rd, Ste 205, Germantown, Maryland 20874, United States
| | - Simon P Webb
- VeraChem LLC, 12850 Middlebrook Rd, Ste 205, Germantown, Maryland 20874, United States
| |
Collapse
|
2
|
Liu T, Simine L. DeltaGzip: Computing Biopolymer-Ligand Binding Affinity via Kolmogorov Complexity and Lossless Compression. J Chem Inf Model 2024; 64:5617-5623. [PMID: 38980667 DOI: 10.1021/acs.jcim.4c00461] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/10/2024]
Abstract
The design of biosequences for biosensing and therapeutics is a challenging multistep search and optimization task. In principle, computational modeling may speed up the design process by virtual screening of sequences based on their binding affinities to target molecules. However, in practice, existing machine-learned models trained to predict binding affinities lack the flexibility with respect to reaction conditions, and molecular dynamics simulations that can incorporate reaction conditions suffer from high computational costs. Here, we describe a computational approach called DeltaGzip that evaluates the free energy of binding in biopolymer-ligand complexes from ultrashort equilibrium molecular dynamics simulations. The entropy of binding is evaluated using the Kolmogorov complexity definition of entropy and approximated using a lossless compression algorithm, Gzip. We benchmark the method on a well-studied data set of protein-ligand complexes comparing the predictions of DeltaGzip to the free energies of binding obtained using Jarzynski equality and experimental measurements.
Collapse
Affiliation(s)
- Tao Liu
- Department of Chemistry, McGill University, Montreal, Quebec H3A 0B8, Canada
| | - Lena Simine
- Department of Chemistry, McGill University, Montreal, Quebec H3A 0B8, Canada
| |
Collapse
|
3
|
Jung J, Yagi K, Tan C, Oshima H, Mori T, Yu I, Matsunaga Y, Kobayashi C, Ito S, Ugarte La Torre D, Sugita Y. GENESIS 2.1: High-Performance Molecular Dynamics Software for Enhanced Sampling and Free-Energy Calculations for Atomistic, Coarse-Grained, and Quantum Mechanics/Molecular Mechanics Models. J Phys Chem B 2024; 128:6028-6048. [PMID: 38876465 PMCID: PMC11215777 DOI: 10.1021/acs.jpcb.4c02096] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/30/2024] [Revised: 05/15/2024] [Accepted: 05/21/2024] [Indexed: 06/16/2024]
Abstract
GENeralized-Ensemble SImulation System (GENESIS) is a molecular dynamics (MD) software developed to simulate the conformational dynamics of a single biomolecule, as well as molecular interactions in large biomolecular assemblies and between multiple biomolecules in cellular environments. To achieve the latter purpose, the earlier versions of GENESIS emphasized high performance in atomistic MD simulations on massively parallel supercomputers, with or without graphics processing units (GPUs). Here, we implemented multiscale MD simulations that include atomistic, coarse-grained, and hybrid quantum mechanics/molecular mechanics (QM/MM) calculations. They demonstrate high performance and are integrated with enhanced conformational sampling algorithms and free-energy calculations without using external programs except for the QM programs. In this article, we review new functions, molecular models, and other essential features in GENESIS version 2.1 and discuss ongoing developments for future releases.
Collapse
Affiliation(s)
- Jaewoon Jung
- Computational
Biophysics Research Team, RIKEN Center for
Computational Science, Kobe, Hyogo 650-0047, Japan
- Theoretical
Molecular Science Laboratory, RIKEN Cluster
for Pioneering Research, Wako, Saitama 351-0198, Japan
| | - Kiyoshi Yagi
- Theoretical
Molecular Science Laboratory, RIKEN Cluster
for Pioneering Research, Wako, Saitama 351-0198, Japan
| | - Cheng Tan
- Computational
Biophysics Research Team, RIKEN Center for
Computational Science, Kobe, Hyogo 650-0047, Japan
| | - Hiraku Oshima
- Laboratory
for Biomolecular Function Simulation, RIKEN
Center for Biosystems Dynamics Research, Kobe, Hyogo 650-0047, Japan
- Graduate
School of Life Science, University of Hyogo, Harima Science Park City, Hyogo 678-1297, Japan
| | - Takaharu Mori
- Theoretical
Molecular Science Laboratory, RIKEN Cluster
for Pioneering Research, Wako, Saitama 351-0198, Japan
- Department
of Chemistry, Tokyo University of Science, Shinjuku-ku, Tokyo 162-8601, Japan
| | - Isseki Yu
- Theoretical
Molecular Science Laboratory, RIKEN Cluster
for Pioneering Research, Wako, Saitama 351-0198, Japan
- Department
of Bioinformatics, Maebashi Institute of
Technology, Maebashi, Gunma 371-0816, Japan
| | - Yasuhiro Matsunaga
- Computational
Biophysics Research Team, RIKEN Center for
Computational Science, Kobe, Hyogo 650-0047, Japan
- Graduate
School of Science and Engineering, Saitama
University, Saitama 338-8570, Japan
| | - Chigusa Kobayashi
- Computational
Biophysics Research Team, RIKEN Center for
Computational Science, Kobe, Hyogo 650-0047, Japan
| | - Shingo Ito
- Theoretical
Molecular Science Laboratory, RIKEN Cluster
for Pioneering Research, Wako, Saitama 351-0198, Japan
| | - Diego Ugarte La Torre
- Computational
Biophysics Research Team, RIKEN Center for
Computational Science, Kobe, Hyogo 650-0047, Japan
| | - Yuji Sugita
- Computational
Biophysics Research Team, RIKEN Center for
Computational Science, Kobe, Hyogo 650-0047, Japan
- Theoretical
Molecular Science Laboratory, RIKEN Cluster
for Pioneering Research, Wako, Saitama 351-0198, Japan
- Laboratory
for Biomolecular Function Simulation, RIKEN
Center for Biosystems Dynamics Research, Kobe, Hyogo 650-0047, Japan
| |
Collapse
|
4
|
Verteramo ML, Ignjatović MM, Kumar R, Wernersson S, Ekberg V, Wallerstein J, Carlström G, Chadimová V, Leffler H, Zetterberg F, Logan DT, Ryde U, Akke M, Nilsson UJ. Interplay of halogen bonding and solvation in protein-ligand binding. iScience 2024; 27:109636. [PMID: 38633000 PMCID: PMC11021960 DOI: 10.1016/j.isci.2024.109636] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/21/2023] [Revised: 02/13/2024] [Accepted: 03/26/2024] [Indexed: 04/19/2024] Open
Abstract
Halogen bonding is increasingly utilized in efforts to achieve high affinity and selectivity of molecules designed to bind proteins, making it paramount to understand the relationship between structure, dynamics, and thermodynamic driving forces. We present a detailed analysis addressing this problem using a series of protein-ligand complexes involving single halogen substitutions - F, Cl, Br, and I - and nearly identical structures. Isothermal titration calorimetry reveals an increasingly favorable binding enthalpy from F to I that correlates with the halogen size and σ-hole electropositive character, but is partially counteracted by unfavorable entropy, which is constant from F to Cl and Br, but worse for I. Consequently, the binding free energy is roughly equal for Cl, Br, and I. QM and solvation-free-energy calculations reflect an intricate balance between halogen bonding, hydrogen bonds, and solvation. These advances have the potential to aid future drug design initiatives involving halogenated compounds.
Collapse
Affiliation(s)
| | | | - Rohit Kumar
- Department of Chemistry, Lund University, Lund, Sweden
| | | | | | | | | | | | - Hakon Leffler
- Microbiology, Immunology, and Glycobiology, Department of Experimental Medicine, Lund University, Lund, Sweden
| | | | | | - Ulf Ryde
- Department of Chemistry, Lund University, Lund, Sweden
| | - Mikael Akke
- Department of Chemistry, Lund University, Lund, Sweden
| | | |
Collapse
|
5
|
Brooks CL, MacKerell AD, Post CB, Nilsson L. Biomolecular dynamics in the 21st century. Biochim Biophys Acta Gen Subj 2024; 1868:130534. [PMID: 38065235 PMCID: PMC10842176 DOI: 10.1016/j.bbagen.2023.130534] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/26/2023] [Revised: 11/28/2023] [Accepted: 11/29/2023] [Indexed: 01/03/2024]
Abstract
The relevance of motions in biological macromolecules has been clear since the early structural analyses of proteins by X-ray crystallography. Computer simulations have been applied to provide a deeper understanding of the dynamics of biological macromolecules since 1976, and are now a standard tool in many labs working on the structure and function of biomolecules. In this mini-review we highlight some areas of current interest and active development for simulations, in particular all-atom molecular dynamics simulations.
Collapse
Affiliation(s)
- Charles L Brooks
- University of Michigan, Department of Chemistry, Ann Arbor, MI 48109, USA.
| | | | - Carol B Post
- Purdue University, Department of Medicinal Chemistry and Molecular Pharmacology, West Lafayette, IN 47907-2091, USA.
| | - Lennart Nilsson
- Karolinska Institutet, Department of Biosciences and Nutrition, SE-14183 Huddinge, Sweden.
| |
Collapse
|
6
|
Donyapour N, Fathi Niazi F, Roussey NM, Bose S, Dickson A. Flexible Topology: A Dynamic Model of a Continuous Chemical Space. J Chem Theory Comput 2023; 19:5088-5098. [PMID: 37487141 PMCID: PMC11060842 DOI: 10.1021/acs.jctc.3c00409] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/26/2023]
Abstract
Ligand design problems involve searching chemical space for a molecule with a set of desired properties. As chemical space is discrete, this search must be conducted in a pointwise manner, separately investigating one molecule at a time, which can be inefficient. We propose a method called "Flexible Topology", where a ligand is composed of a set of shapeshifting "ghost" atoms, whose atomic identities and connectivity can dynamically change over the course of a simulation. Ghost atoms are guided toward their target positions using a translation-, rotation-, and index-invariant restraint potential. This is the first step toward a continuous model of chemical space, where a dynamic simulation can move from one molecule to another by following gradients of a potential energy function. This builds on a substantial history of alchemy in the field of molecular dynamics simulation, including the Lambda dynamics method developed by Brooks and co-workers [X. Kong and C.L. Brooks III, J. Chem. Phys. 105, 2414 (1996)], but takes it to an extreme by associating a set of four dynamical attributes with each shapeshifting ghost atom that control not only its presence but also its atomic identity. Here, we outline the theoretical details of this method, its implementation using the OpenMM simulation package, and some preliminary studies of ghost particle assembly simulations in vacuum. We examine a set of 10 small molecules, ranging in size from 6 to 50 atoms, and show that Flexible Topology is able to consistently assemble all of these molecules to high accuracy, beginning from randomly initialized positions and attributes.
Collapse
Affiliation(s)
- Nazanin Donyapour
- Department of Computational Mathematics, Science & Engineering, Michigan State University, East Lansing, Michigan 48824, United States
| | - Fatemeh Fathi Niazi
- Department of Computational Mathematics, Science & Engineering, Michigan State University, East Lansing, Michigan 48824, United States
| | - Nicole M Roussey
- Department of Biochemistry & Molecular Biology, Michigan State University, East Lansing, Michigan 48824, United States
| | - Samik Bose
- Department of Biochemistry & Molecular Biology, Michigan State University, East Lansing, Michigan 48824, United States
| | - Alex Dickson
- Department of Computational Mathematics, Science & Engineering, Michigan State University, East Lansing, Michigan 48824, United States
- Department of Biochemistry & Molecular Biology, Michigan State University, East Lansing, Michigan 48824, United States
| |
Collapse
|
7
|
Li XS, Van Koten B, Dinner AR, Thiede EH. Understanding the sources of error in MBAR through asymptotic analysis. J Chem Phys 2023; 158:214107. [PMID: 37259996 PMCID: PMC10845960 DOI: 10.1063/5.0147243] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/20/2023] [Accepted: 05/08/2023] [Indexed: 06/02/2023] Open
Abstract
Many sampling strategies commonly used in molecular dynamics, such as umbrella sampling and alchemical free energy methods, involve sampling from multiple states. The Multistate Bennett Acceptance Ratio (MBAR) formalism is a widely used way of recombining the resulting data. However, the error of the MBAR estimator is not well-understood: previous error analyses of MBAR assumed independent samples. In this work, we derive a central limit theorem for MBAR estimates in the presence of correlated data, further justifying the use of MBAR in practical applications. Moreover, our central limit theorem yields an estimate of the error that can be decomposed into contributions from the individual Markov chains used to sample the states. This gives additional insight into how sampling in each state affects the overall error. We demonstrate our error estimator on an umbrella sampling calculation of the free energy of isomerization of the alanine dipeptide and an alchemical calculation of the hydration free energy of methane. Our numerical results demonstrate that the time required for the Markov chain to decorrelate in individual states can contribute considerably to the total MBAR error, highlighting the importance of accurately addressing the effect of sample correlation.
Collapse
Affiliation(s)
- Xiang Sherry Li
- Department of Chemistry and James Franck Institute, University of Chicago, Chicago, Illinois 60637, USA
| | - Brian Van Koten
- Department of Mathematics and Statistics, University of Massachusetts, Amherst, Massachusetts 01003, USA
| | - Aaron R. Dinner
- Department of Chemistry and James Franck Institute, University of Chicago, Chicago, Illinois 60637, USA
| | - Erik H. Thiede
- Center for Computational Mathematics, Flatiron Institute, New York, New York 10010, USA
| |
Collapse
|
8
|
Xu H. The slow but steady rise of binding free energy calculations in drug discovery. J Comput Aided Mol Des 2023; 37:67-74. [PMID: 36469232 DOI: 10.1007/s10822-022-00494-x] [Citation(s) in RCA: 6] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/01/2022] [Accepted: 11/28/2022] [Indexed: 12/12/2022]
Abstract
Binding free energy calculations are increasingly used in drug discovery research to predict protein-ligand binding affinities and to prioritize candidate drug molecules accordingly. It has taken decades of collective effort to transform this academic concept into a technology adopted by the pharmaceutical and biotech industry. Having personally witnessed and taken part in this transformation, here I recount the (incomplete) list of problems that had to be solved to make this computational tool practical and suggest areas of future development.
Collapse
Affiliation(s)
- Huafeng Xu
- Roivant Discovery, 151 West 42nd Street, New York, NY, 10036, USA.
| |
Collapse
|
9
|
Lee KH, Kuczera K. Effect of alanine versus serine at position 88 of human transthyretin mutants on the protein stability. Protein Eng Des Sel 2023; 36:6972274. [PMID: 36611015 DOI: 10.1093/protein/gzad001] [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: 09/21/2022] [Revised: 10/22/2022] [Accepted: 12/05/2022] [Indexed: 01/09/2023] Open
Abstract
Human transthyretin (TTR) is a homo-tetrameric plasma protein associated with a high percentage of β-sheet forming amyloid fibrils. It accumulates in tissues or extracellular matrices to cause amyloid diseases. Free energy simulations with thermodynamic integration based on all-atom molecular dynamics simulations have been carried out to analyze the effects of the His88 → Ala and Ser mutations on the stability of human TTR. The calculated free energy change differences (ΔΔG) caused by the His88 → Ala and His88 → Ser mutations are -1.84 ± 0.86 and 7.56 ± 0.55 kcal/mol, respectively, which are in excellent agreement with prior reported experimental values. The simulation results show that the H88A mutant is more stable than the wild type, whereas the H88S mutant is less stable than the wild type. The free energy component analysis shows that the contribution to the free energy change difference (ΔΔG) for the His88 → Ala and His88 → Ser mutations mainly arise from electrostatic and van der Waals interactions, respectively. The electrostatic term stabilizes the H88A mutant more than the wild type, but the van der Waals interaction destabilizes the H88S mutant relative to the wild type. Individual residue contributions to the free energy change show neighboring residues exert stabilizing and destabilizing influence on the mutants. The implications of the simulation results for understanding the stabilizing and destabilizing effect and its contribution to protein stability are discussed.
Collapse
Affiliation(s)
- Kyung-Hoon Lee
- Department of Biology, Chowan University, One University Place, Murfreesboro, NC 27855, USA
| | - Krzysztof Kuczera
- Department of Chemistry and Department of Molecular Biosciences, University of Kansas, 1567 Irving Hill Road, Lawrence, KS 66045, USA
| |
Collapse
|
10
|
Tsai HC, Lee TS, Ganguly A, Giese TJ, Ebert MCCJC, Labute P, Merz KM, York DM. AMBER Free Energy Tools: A New Framework for the Design of Optimized Alchemical Transformation Pathways. J Chem Theory Comput 2023; 19:10.1021/acs.jctc.2c00725. [PMID: 36622640 PMCID: PMC10329732 DOI: 10.1021/acs.jctc.2c00725] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/10/2023]
Abstract
We develop a framework for the design of optimized alchemical transformation pathways in free energy simulations using nonlinear mixing and a new functional form for so-called "softcore" potentials. We describe the implementation and testing of this framework in the GPU-accelerated AMBER software suite. The new optimized alchemical transformation pathways integrate a number of important features, including (1) the use of smoothstep functions to stabilize behavior near the transformation end points, (2) consistent power scaling of Coulomb and Lennard-Jones (LJ) interactions with unitless control parameters to maintain balance of electrostatic attractions and exchange repulsions, (3) pairwise form based on the LJ contact radius for the effective interaction distance with separation-shifted scaling, and (4) rigorous smoothing of the potential at the nonbonded cutoff boundary. The new softcore potential form is combined with smoothly transforming nonlinear λ weights for mixing specific potential energy terms, along with flexible λ-scheduling features, to enable robust and stable alchemical transformation pathways. The resulting pathways are demonstrated and tested, and shown to be superior to the traditional methods in terms of numerical stability and minimal variance of the free energy estimates for all cases considered. The framework presented here can be used to design new alchemical enhanced sampling methods, and leveraged in robust free energy workflows for large ligand data sets.
Collapse
Affiliation(s)
- Hsu-Chun Tsai
- Laboratory for Biomolecular Simulation Research, Institute for Quantitative Biomedicine and Department of Chemistry and Chemical Biology, Rutgers University, Piscataway, NJ 08854, USA
| | - Tai-Sung Lee
- Laboratory for Biomolecular Simulation Research, Institute for Quantitative Biomedicine and Department of Chemistry and Chemical Biology, Rutgers University, Piscataway, NJ 08854, USA
| | - Abir Ganguly
- Laboratory for Biomolecular Simulation Research, Institute for Quantitative Biomedicine and Department of Chemistry and Chemical Biology, Rutgers University, Piscataway, NJ 08854, USA
| | - Timothy J. Giese
- Laboratory for Biomolecular Simulation Research, Institute for Quantitative Biomedicine and Department of Chemistry and Chemical Biology, Rutgers University, Piscataway, NJ 08854, USA
| | - Maximilian CCJC Ebert
- Congruence Therapeutics, 7171 Rue Frederick Banting #117, Saint-Laurent, Quebec, Canada H4S 1Z9
| | - Paul Labute
- Chemical Computing Group ULC, 910-1010 Sherbrooke West, Montreal, Quebec, Canada H3A 2R7
| | - Kenneth M. Merz
- Department of Chemistry and Department of Biochemistry and Molecular Biology, Michigan State University, East Lansing, MI 48824, USA
| | - Darrin M. York
- Laboratory for Biomolecular Simulation Research, Institute for Quantitative Biomedicine and Department of Chemistry and Chemical Biology, Rutgers University, Piscataway, NJ 08854, USA
| |
Collapse
|
11
|
Jeong KJ, Jeong S, Lee S, Son CY. Predictive Molecular Models for Charged Materials Systems: From Energy Materials to Biomacromolecules. ADVANCED MATERIALS (DEERFIELD BEACH, FLA.) 2023; 35:e2204272. [PMID: 36373701 DOI: 10.1002/adma.202204272] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 05/11/2022] [Revised: 07/05/2022] [Indexed: 06/16/2023]
Abstract
Electrostatic interactions play a dominant role in charged materials systems. Understanding the complex correlation between macroscopic properties with microscopic structures is of critical importance to develop rational design strategies for advanced materials. But the complexity of this challenging task is augmented by interfaces present in the charged materials systems, such as electrode-electrolyte interfaces or biological membranes. Over the last decades, predictive molecular simulations that are founded in fundamental physics and optimized for charged interfacial systems have proven their value in providing molecular understanding of physicochemical properties and functional mechanisms for diverse materials. Novel design strategies utilizing predictive models have been suggested as promising route for the rational design of materials with tailored properties. Here, an overview of recent advances in the understanding of charged interfacial systems aided by predictive molecular simulations is presented. Focusing on three types of charged interfaces found in energy materials and biomacromolecules, how the molecular models characterize ion structure, charge transport, morphology relation to the environment, and the thermodynamics/kinetics of molecular binding at the interfaces is discussed. The critical analysis brings two prominent field of energy materials and biological science under common perspective, to stimulate crossover in both research field that have been largely separated.
Collapse
Affiliation(s)
- Kyeong-Jun Jeong
- Department of Chemistry, Pohang University of Science and Technology (POSTECH), Pohang, 790-784, South Korea
| | - Seungwon Jeong
- Department of Chemistry, Pohang University of Science and Technology (POSTECH), Pohang, 790-784, South Korea
| | - Sangmin Lee
- Department of Chemistry, Pohang University of Science and Technology (POSTECH), Pohang, 790-784, South Korea
| | - Chang Yun Son
- Department of Chemistry, Pohang University of Science and Technology (POSTECH), Pohang, 790-784, South Korea
| |
Collapse
|
12
|
Matsunaga Y, Kamiya M, Oshima H, Jung J, Ito S, Sugita Y. Use of multistate Bennett acceptance ratio method for free-energy calculations from enhanced sampling and free-energy perturbation. Biophys Rev 2022; 14:1503-1512. [PMID: 36659993 PMCID: PMC9842838 DOI: 10.1007/s12551-022-01030-9] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/24/2022] [Accepted: 12/01/2022] [Indexed: 12/15/2022] Open
Abstract
Multistate Bennett acceptance ratio (MBAR) works as a method to analyze molecular dynamics (MD) simulation data after the simulations have been finished. It is widely used to estimate free-energy changes between different states and averaged properties at the states of interest. MBAR allows us to treat a wide range of states from those at different temperature/pressure to those with different model parameters. Due to the broad applicability, the MBAR equations are rather difficult to apply for free-energy calculations using different types of MD simulations including enhanced conformational sampling methods and free-energy perturbation. In this review, we first summarize the basic theory of the MBAR equations and categorize the representative usages into the following four: (i) perturbation, (ii) scaling, (iii) accumulation, and (iv) full potential energy. For each, we explain how to prepare input data using MD simulation trajectories for solving the MBAR equations. MBAR is also useful to estimate reliable free-energy differences using MD trajectories based on a semi-empirical quantum mechanics/molecular mechanics (QM/MM) model and ab initio QM/MM energy calculations on the MD snapshots. We also explain how to use the MBAR software in the GENESIS package, which we call mbar_analysis, for the four representative cases. The proposed estimations of free-energy changes and thermodynamic averages are effective and useful for various biomolecular systems.
Collapse
Affiliation(s)
- Yasuhiro Matsunaga
- grid.263023.60000 0001 0703 3735Graduate School of Science and Engineering, Saitama University, Saitama, Saitama 338-8570 Japan
| | - Motoshi Kamiya
- grid.467196.b0000 0001 2285 6123Institute for Molecular Science, Myodaiji, Okazaki, Aichi 444-8585 Japan
| | - Hiraku Oshima
- grid.508743.dLaboratory for Biomolecular Function Simulation, RIKEN Center for Biosystems Dynamics Research, Kobe, Hyogo 650-0047 Japan
| | - Jaewoon Jung
- grid.474693.bComputational Biophysics Research Team, RIKEN Center for Computational Science, Kobe, Hyogo 650-0047 Japan
| | - Shingo Ito
- grid.7597.c0000000094465255Theoretical Molecular Science Laboratory, RIKEN Cluster for Pioneering Research, Wako, Saitama 351-0198 Japan
| | - Yuji Sugita
- grid.508743.dLaboratory for Biomolecular Function Simulation, RIKEN Center for Biosystems Dynamics Research, Kobe, Hyogo 650-0047 Japan ,grid.474693.bComputational Biophysics Research Team, RIKEN Center for Computational Science, Kobe, Hyogo 650-0047 Japan ,grid.7597.c0000000094465255Theoretical Molecular Science Laboratory, RIKEN Cluster for Pioneering Research, Wako, Saitama 351-0198 Japan
| |
Collapse
|
13
|
Karwounopoulos J, Wieder M, Boresch S. Relative binding free energy calculations with transformato: A molecular dynamics engine-independent tool. Front Mol Biosci 2022; 9:954638. [PMID: 36148009 PMCID: PMC9485484 DOI: 10.3389/fmolb.2022.954638] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/27/2022] [Accepted: 07/18/2022] [Indexed: 11/21/2022] Open
Abstract
We present the software package transformato for the setup of large-scale relative binding free energy calculations. Transformato is written in Python as an open source project (https://github.com/wiederm/transformato); in contrast to comparable tools, it is not closely tied to a particular molecular dynamics engine to carry out the underlying simulations. Instead of alchemically transforming a ligand L1 directly into another L2, the two ligands are mutated to a common core. Thus, while dummy atoms are required at intermediate states, in particular at the common core state, none are present at the physical endstates. To validate the method, we calculated 76 relative binding free energy differences ΔΔGL1→L2bind for five protein–ligand systems. The overall root mean squared error to experimental binding free energies is 1.17 kcal/mol with a Pearson correlation coefficient of 0.73. For selected cases, we checked that the relative binding free energy differences between pairs of ligands do not depend on the choice of the intermediate common core structure. Additionally, we report results with and without hydrogen mass reweighting. The code currently supports OpenMM, CHARMM, and CHARMM/OpenMM directly. Since the program logic to choose and construct alchemical transformation paths is separated from the generation of input and topology/parameter files, extending transformato to support additional molecular dynamics engines is straightforward.
Collapse
Affiliation(s)
- Johannes Karwounopoulos
- Faculty of Chemistry, Institute of Computational Biological Chemistry, University of Vienna, Vienna, Austria
- Vienna Doctoral School of Chemistry (DoSChem), University of Vienna, Vienna, Austria
| | - Marcus Wieder
- Department of Pharmaceutical Sciences, Faculty of Life Sciences, University of Vienna, Vienna, Austria
- *Correspondence: Marcus Wieder, ; Stefan Boresch,
| | - Stefan Boresch
- Faculty of Chemistry, Institute of Computational Biological Chemistry, University of Vienna, Vienna, Austria
- *Correspondence: Marcus Wieder, ; Stefan Boresch,
| |
Collapse
|
14
|
Lee KH, Kuczera K. Free energy simulations to study mutational effect of a conserved residue, Trp24, on stability of human serum retinol-binding protein. J Biomol Struct Dyn 2022:1-11. [PMID: 35899456 DOI: 10.1080/07391102.2022.2100829] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/16/2022]
Abstract
Human serum retinol-binding protein (RBP) is a plasma transport protein for vitamin A. RBP is a prime subclass of lipocalins, which bind nonpolar ligands within a β-barrel. To understand the role of Trp 24, one of the highly conserved residues in RBP, free energy simulations have been carried out to understand the effects of the mutations from Trp at position 24 to Leu, Phe, and Tyr in the apo-RBP on its thermal stability. We examine various unfolded systems to study the dependence of the free energy differences on the denatured structure. Our calculated free energy difference values for the three mutations are in excellent agreement with the experimental values when the initial coordinates of the seven-residue peptide segments truncated from the crystal structure are used for the denatured systems. Our free energy change differences for the Trp→Leu, Trp→Phe, and Trp→Tyr mutations are 2.50 ± 0.69, 2.58 ± 0.50, and 2.49 ± 0.48 kcal/mol, respectively, when the native-like seven-residue peptides are used as models for the denatured systems. The main contributions to the free energy change differences for the Trp24→Leu and Trp24→Phe mutations are mainly from van der Waals and covalent interactions, respectively. Electrostatic, van der Waals and covalent terms equally contribute to the free energy change difference for the Trp24→Tyr mutation. The free energy simulation helps understand the detailed microscopic mechanism of the stability of the RBP mutants relative to the wild type and the role of the highly conserved residue, Trp24, of the human RBP.Communicated by Ramaswamy H. Sarma.
Collapse
Affiliation(s)
- Kyung-Hoon Lee
- Department of Biology, Chowan University, Murfreesboro, NC, USA
| | - Krzysztof Kuczera
- Department of Chemistry and Department of Molecular Biosciences, University of Kansas, Lawrence, KS, USA
| |
Collapse
|
15
|
Lee K, Kuczera K. Modulation of human transthyretin stability by the mutations at histidine 88 studied by free energy simulation. Proteins 2022; 90:1825-1836. [DOI: 10.1002/prot.26353] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/16/2022] [Revised: 04/03/2022] [Accepted: 04/20/2022] [Indexed: 11/10/2022]
Affiliation(s)
- Kyung‐Hoon Lee
- Department of Biology Chowan University Murfreesboro North Carolina USA
| | - Krzysztof Kuczera
- Department of Chemistry and Department of Molecular Biosciences University of Kansas Lawrence Kansas USA
| |
Collapse
|
16
|
Lee KH, Kuczera K. Simulation analysis of selective alanine mutation effect on stability of human prion protein. J Biomol Struct Dyn 2022; 41:2619-2629. [PMID: 35176965 DOI: 10.1080/07391102.2022.2036237] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/19/2022]
Abstract
Prion diseases are neurodegenerative disorders caused by spongiform degeneration of the brain. Understanding the fundamental mechanism of prion protein aggregation caused by mutations is very crucial to resolve the pathology of prion diseases. To help understand the roles of individual residues on the stability of the human prion protein, the computational method of free energy simulations based on atomistic molecular dynamics trajectories is applied to Phe175 → Ala, Val180 → Ala, and Val209 → Ala mutations of the human prion protein. The simulations show that all three alanine mutations destabilize the human prion protein. The calculated free energy change differences, ΔΔG, for the Phe175 → Ala, Val180 → Ala, and Val209 → Ala mutations are in good agreement with the experimental values. The significant destabilizing effects on the mutants relative to the wild-type protein arise from van der Waals terms. Furthermore, our free energy decomposition analysis shows that the major contribution to destabilizing the V180A and V209A mutants relative to the wild-type protein is originated from van der Waals interactions from residues near the mutation sites. In contrast, the contribution to destabilizing the F175A mutant is mainly caused by van der Waals interactions from residues near and far away from the mutation site. Our results show that the free energy simulation with a thermodynamic integration approach for selected alanine scanning mutations is beneficial for understanding the detailed mechanism of human prion protein destabilization, specific residues' role, and the hydrophobic effect on protein stability.Communicated by Ramaswamy H. Sarma.
Collapse
Affiliation(s)
- Kyung-Hoon Lee
- Department of Biology, Chowan University, Murfreesboro, NC, USA
| | - Krzysztof Kuczera
- Department of Chemistry and Department of Molecular Biosciences, University of Kansas, Lawrence, KS, USA
| |
Collapse
|
17
|
Dickson CJ, Hornak V, Duca JS. Relative Binding Free-Energy Calculations at Lipid-Exposed Sites: Deciphering Hot Spots. J Chem Inf Model 2021; 61:5923-5930. [PMID: 34843243 DOI: 10.1021/acs.jcim.1c01147] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Relative binding free-energy (RBFE) calculations are experiencing resurgence in the computer-aided drug design of novel small molecules due to performance gains allowed by cutting-edge molecular mechanic force fields and computer hardware. Application of RBFE to soluble proteins is becoming a routine, while recent studies outline necessary steps to successfully apply RBFE at the orthosteric site of membrane-embedded G-protein-coupled receptors (GPCRs). In this work, we apply RBFE to a congeneric series of antagonists that bind to a lipid-exposed, extra-helical site of the P2Y1 receptor. We find promising performance of RBFE, such that it may be applied in a predictive manner on drug discovery programs targeting lipid-exposed sites. Further, by the application of the microkinetic model, binding at a lipid-exposed site can be split into (1) membrane partitioning of the drug molecule followed by (2) binding at the extra-helical site. We find that RBFE can be applied to calculate the free energy of each step, allowing the uncoupling of observed binding free energy from the influence of membrane affinity. This protocol may be used to identify binding hot spots at extra-helical sites and guide drug discovery programs toward optimizing intrinsic activity at the target.
Collapse
Affiliation(s)
- Callum J Dickson
- Computer-Aided Drug Discovery, Global Discovery Chemistry, Novartis Institutes for BioMedical Research, 181 Massachusetts Avenue, Cambridge, Massachusetts 02139, United States
| | - Viktor Hornak
- Computer-Aided Drug Discovery, Global Discovery Chemistry, Novartis Institutes for BioMedical Research, 181 Massachusetts Avenue, Cambridge, Massachusetts 02139, United States
| | - Jose S Duca
- Computer-Aided Drug Discovery, Global Discovery Chemistry, Novartis Institutes for BioMedical Research, 181 Massachusetts Avenue, Cambridge, Massachusetts 02139, United States
| |
Collapse
|
18
|
Ben-Shalom IY, Lin C, Radak BK, Sherman W, Gilson MK. Fast Equilibration of Water between Buried Sites and the Bulk by Molecular Dynamics with Parallel Monte Carlo Water Moves on Graphical Processing Units. J Chem Theory Comput 2021; 17:7366-7372. [PMID: 34762421 PMCID: PMC8716912 DOI: 10.1021/acs.jctc.1c00867] [Citation(s) in RCA: 9] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022]
Abstract
Molecular dynamics (MD) simulations of proteins are commonly used to sample from the Boltzmann distribution of conformational states, with wide-ranging applications spanning chemistry, biophysics, and drug discovery. However, MD can be inefficient at equilibrating water occupancy for buried cavities in proteins that are inaccessible to the surrounding solvent. Indeed, the time needed for water molecules to equilibrate between the bulk solvent and the binding site can be well beyond what is practical with standard MD, which typically ranges from hundreds of nanoseconds to a few microseconds. We recently introduced a hybrid Monte Carlo/MD (MC/MD) method, which speeds up the equilibration of water between buried cavities and the surrounding solvent, while sampling from the thermodynamically correct distribution of states. While the initial implementation of the MC functionality led to considerable slowing of the overall simulations, here we address this problem with a parallel MC algorithm implemented on graphical processing units. This results in speed-ups of 10-fold to 1000-fold over the original MC/MD algorithm, depending on the system and simulation parameters. The present method is available for use in the AMBER simulation software.
Collapse
Affiliation(s)
- Ido Y. Ben-Shalom
- Skaggs School of Pharmacy and Pharmaceutical Sciences, University of California, San Diego, 92093 La Jolla, California, USA
| | - Charles Lin
- Roivant Discovery, Boston, Massachusetts, 02110, USA
| | | | - Woody Sherman
- Roivant Discovery, Boston, Massachusetts, 02110, USA
| | - Michael K. Gilson
- Skaggs School of Pharmacy and Pharmaceutical Sciences, University of California, San Diego, 92093 La Jolla, California, USA,To whom correspondence should be addressed,
| |
Collapse
|
19
|
Thomaston JL, Samways ML, Konstantinidi A, Ma C, Hu Y, Bruce Macdonald HE, Wang J, Essex JW, DeGrado WF, Kolocouris A. Rimantadine Binds to and Inhibits the Influenza A M2 Proton Channel without Enantiomeric Specificity. Biochemistry 2021; 60:10.1021/acs.biochem.1c00437. [PMID: 34342217 PMCID: PMC8810914 DOI: 10.1021/acs.biochem.1c00437] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/08/2023]
Abstract
The influenza A M2 wild-type (WT) proton channel is the target of the anti-influenza drug rimantadine. Rimantadine has two enantiomers, though most investigations into drug binding and inhibition have used a racemic mixture. Solid-state NMR experiments using the full length-M2 WT have shown significant spectral differences that were interpreted to indicate tighter binding for (R)- vs (S)-rimantadine. However, it was unclear if this correlates with a functional difference in drug binding and inhibition. Using X-ray crystallography, we have determined that both (R)- and (S)-rimantadine bind to the M2 WT pore with slight differences in the hydration of each enantiomer. However, this does not result in a difference in potency or binding kinetics, as shown by similar values for kon, koff, and Kd in electrophysiological assays and for EC50 values in cellular assays. We concluded that the slight differences in hydration for the (R)- and (S)-rimantadine enantiomers are not relevant to drug binding or channel inhibition. To further explore the effect of the hydration of the M2 pore on binding affinity, the water structure was evaluated by grand canonical ensemble molecular dynamics simulations as a function of the chemical potential of the water. Initially, the two layers of ordered water molecules between the bound drug and the channel's gating His37 residues mask the drug's chirality. As the chemical potential becomes more unfavorable, the drug translocates down to the lower water layer, and the interaction becomes more sensitive to chirality. These studies suggest the feasibility of displacing the upper water layer and specifically recognizing the lower water layers in novel drugs.
Collapse
Affiliation(s)
- Jessica L Thomaston
- Laboratory of Medicinal Chemistry, Section of Pharmaceutical Chemistry, Department of Pharmacy, University of California, San Francisco, California 94158, United States
| | - Marley L Samways
- School of Chemistry, University of Southampton, Southampton SO17 1BJ, United Kingdom
| | - Athina Konstantinidi
- Department of Pharmaceutical Chemistry, School of Pharmacy, National and Kapodistrian University of Athens, 15771 Athens, Greece
| | - Chunlong Ma
- Department of Pharmacology and Toxicology, College of Pharmacy, University of Arizona, Tucson, Arizona 85721, United States
| | - Yanmei Hu
- Department of Pharmacology and Toxicology, College of Pharmacy, University of Arizona, Tucson, Arizona 85721, United States
| | - Hannah E Bruce Macdonald
- Computational and Systems Biology Program, Memorial Sloan Kettering Cancer Center, New York, New York 10065, United States
| | - Jun Wang
- Department of Pharmacology and Toxicology, College of Pharmacy, University of Arizona, Tucson, Arizona 85721, United States
| | - Jonathan W Essex
- School of Chemistry, University of Southampton, Southampton SO17 1BJ, United Kingdom
| | - William F DeGrado
- Laboratory of Medicinal Chemistry, Section of Pharmaceutical Chemistry, Department of Pharmacy, University of California, San Francisco, California 94158, United States
| | - Antonios Kolocouris
- Department of Pharmaceutical Chemistry, School of Pharmacy, National and Kapodistrian University of Athens, 15771 Athens, Greece
| |
Collapse
|
20
|
Ali HS, Chakravorty A, Kalayan J, de Visser SP, Henchman RH. Energy-entropy method using multiscale cell correlation to calculate binding free energies in the SAMPL8 host-guest challenge. J Comput Aided Mol Des 2021; 35:911-921. [PMID: 34264476 PMCID: PMC8367938 DOI: 10.1007/s10822-021-00406-5] [Citation(s) in RCA: 10] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/27/2021] [Accepted: 06/22/2021] [Indexed: 11/29/2022]
Abstract
Free energy drives a wide range of molecular processes such as solvation, binding, chemical reactions and conformational change. Given the central importance of binding, a wide range of methods exist to calculate it, whether based on scoring functions, machine-learning, classical or electronic structure methods, alchemy, or explicit evaluation of energy and entropy. Here we present a new energy-entropy (EE) method to calculate the host-guest binding free energy directly from molecular dynamics (MD) simulation. Entropy is evaluated using Multiscale Cell Correlation (MCC) which uses force and torque covariance and contacts at two different length scales. The method is tested on a series of seven host-guest complexes in the SAMPL8 (Statistical Assessment of the Modeling of Proteins and Ligands) "Drugs of Abuse" Blind Challenge. The EE-MCC binding free energies are found to agree with experiment with an average error of 0.9 kcal mol-1. MCC makes clear the origin of the entropy changes, showing that the large loss of positional, orientational, and to a lesser extent conformational entropy of each binding guest is compensated for by a gain in orientational entropy of water released to bulk, combined with smaller decreases in vibrational entropy of the host, guest and contacting water.
Collapse
Affiliation(s)
- Hafiz Saqib Ali
- Manchester Institute of Biotechnology, The University of Manchester, 131 Princess Street, Manchester, M1 7DN, UK
- Department of Chemistry, The University of Manchester, Oxford Road, Manchester, M13 9PL, UK
| | - Arghya Chakravorty
- Department of Chemistry, University of Michigan, Ann Arbor, MI, 48109, USA
| | - Jas Kalayan
- Manchester Institute of Biotechnology, The University of Manchester, 131 Princess Street, Manchester, M1 7DN, UK
- Department of Chemistry, The University of Manchester, Oxford Road, Manchester, M13 9PL, UK
| | - Samuel P de Visser
- Manchester Institute of Biotechnology, The University of Manchester, 131 Princess Street, Manchester, M1 7DN, UK
- Department of Chemical Engineering and Analytical Science, The University of Manchester, Oxford Road, Manchester, M13 9PL, UK
| | - Richard H Henchman
- Manchester Institute of Biotechnology, The University of Manchester, 131 Princess Street, Manchester, M1 7DN, UK.
- Department of Chemistry, The University of Manchester, Oxford Road, Manchester, M13 9PL, UK.
- Sydney Medical School, The University of Sydney, Sydney, NSW, 2006, Australia.
| |
Collapse
|
21
|
Fleck M, Wieder M, Boresch S. Dummy Atoms in Alchemical Free Energy Calculations. J Chem Theory Comput 2021; 17:4403-4419. [PMID: 34125525 PMCID: PMC8280730 DOI: 10.1021/acs.jctc.0c01328] [Citation(s) in RCA: 20] [Impact Index Per Article: 6.7] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/23/2020] [Indexed: 02/07/2023]
Abstract
In calculations of relative free energy differences, the number of atoms of the initial and final states is rarely the same. This necessitates the introduction of dummy atoms. These placeholders interact with the physical system only by bonded energy terms. We investigate the conditions necessary so that the presence of dummy atoms does not influence the result of a relative free energy calculation. On the one hand, one has to ensure that dummy atoms only give a multiplicative contribution to the partition function so that their contribution cancels from double-free energy differences. On the other hand, the bonded terms used to attach a dummy atom (or group of dummy atoms) to the physical system have to maintain it in a well-defined position and orientation relative to the physical system. A detailed theoretical analysis of both aspects is provided, illustrated by 24 calculations of relative solvation free energy differences, for which all four legs of the underlying thermodynamic cycle were computed. Cycle closure (or lack thereof) was used as a sensitive indicator to probing the effects of dummy atom treatment on the resulting free energy differences. We find that a naive (but often practiced) treatment of dummy atoms results in errors of up to kBT when calculating the relative solvation free energy difference between two small solutes, such as methane and ammonia. While our analysis focuses on the so-called single topology approach to set up alchemical transformations, similar considerations apply to dual topology, at least many widely used variants thereof.
Collapse
Affiliation(s)
- Markus Fleck
- Faculty
of Chemistry, Department of Computational Biological Chemistry, University of Vienna, Währingerstraße 17, A-1090 Vienna, Austria
| | - Marcus Wieder
- Department
of Pharmaceutical Sciences, Faculty of Life Sciences, University of Vienna, Althanstraße 14, 1090 Vienna, Austria
| | - Stefan Boresch
- Faculty
of Chemistry, Department of Computational Biological Chemistry, University of Vienna, Währingerstraße 17, A-1090 Vienna, Austria
| |
Collapse
|
22
|
Lee KH, Kuczera K. Free energy simulations to understand the effect of Met → Ala mutations at positions 205, 206 and 213 on stability of human prion protein. Biophys Chem 2021; 275:106620. [PMID: 34058726 DOI: 10.1016/j.bpc.2021.106620] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/10/2021] [Revised: 05/05/2021] [Accepted: 05/11/2021] [Indexed: 01/23/2023]
Abstract
Prion diseases are a family of infectious amyloid diseases affecting human and animals. Prion propagation in transmissible spongiform encephalopathies is associated with the unfolding and conversion of normal cellular prion protein into its pathogenic scrapie form. Understanding the fundamentals of prion protein aggregation caused by mutations is crucial to unravel the pathology of prion diseases. To help understand the contributions of individual residues to the stability of the human prion protein, we have carried out free energy simulations based on atomistic molecular dynamics trajectories. We focus on Met → Ala mutations at positions 205, 206 and 213, which are mostly buried residues located on helix 3 of the protein. The simulations predicted that all three mutations destabilize the prion protein. Changes in unfolding free energy upon mutation, ∆∆G, are 3.10 ± 0.79, 2.00 ± 0.26 and 3.06 ± 0.66 kcal/mol for M205A, M206A and M213A, respectively, in excellent agreement with the corresponding experimental values of 3.09 ± 0.28, 1.50 ± 0.34 and 3.12 ± 0.27 kcal/mol [T. Hart et al. (2009) PNAS 106, 5651-5656]. Component analysis indicates that the major contributions to the loss of protein stability arise from van der Waals interactions for the M205A and M206A mutations, and from van der Waals and covalent energy terms for M213A. Interestingly, while free energy contributions from a majority of residues neighboring the mutation sites tend to stabilize the wild type, there are a few residues stabilizing the mutant side chains. Our results show that this approach to free energy calculation can be very useful for understanding the detailed mechanism of human prion protein stability.
Collapse
Affiliation(s)
- Kyung-Hoon Lee
- Department of Biology, Chowan University, One University Drive, Murfreesboro, NC 27855, United States of America.
| | - Krzysztof Kuczera
- Department of Chemistry and Department of Molecular Biosciences, University of Kansas, 1567 Irving Hill Road, Lawrence, KS 66045, United States of America
| |
Collapse
|
23
|
Yagi K, Ito S, Sugita Y. Exploring the Minimum-Energy Pathways and Free-Energy Profiles of Enzymatic Reactions with QM/MM Calculations. J Phys Chem B 2021; 125:4701-4713. [PMID: 33914537 PMCID: PMC10986901 DOI: 10.1021/acs.jpcb.1c01862] [Citation(s) in RCA: 21] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Understanding molecular mechanisms of enzymatic reactions is of vital importance in biochemistry and biophysics. Here, we introduce new functions of hybrid quantum mechanical/molecular mechanical (QM/MM) calculations in the GENESIS program to compute the minimum-energy pathways (MEPs) and free-energy profiles of enzymatic reactions. For this purpose, an interface in GENESIS is developed to utilize a highly parallel electronic structure program, QSimulate-QM (https://qsimulate.com), calling it as a shared library from GENESIS. Second, algorithms to search the MEP are implemented, combining the string method (E et al. J. Chem. Phys. 2007, 126, 164103) with the energy minimization of the buffer MM region. The method implemented in GENESIS is applied to an enzyme, triosephosphate isomerase, which converts dihyroxyacetone phosphate to glyceraldehyde 3-phosphate in four proton-transfer processes. QM/MM-molecular dynamics simulations show performances of greater than 1 ns/day with the density functional tight binding (DFTB), and 10-30 ps/day with the hybrid density functional theory, B3LYP-D3. These performances allow us to compute not only MEP but also the potential of mean force (PMF) of the enzymatic reactions using the QM/MM calculations. The barrier height obtained as 13 kcal mol-1 with B3LYP-D3 in the QM/MM calculation is in agreement with the experimental results. The impact of conformational sampling in PMF calculations and the level of electronic structure calculations (DFTB vs B3LYP-D3) suggests reliable computational protocols for enzymatic reactions without high computational costs.
Collapse
Affiliation(s)
- Kiyoshi Yagi
- Theoretical
Molecular Science Laboratory, RIKEN Cluster
for Pioneering Research, 2-1 Hirosawa, Wako, Saitama 351-0198, Japan
| | - Shingo Ito
- Theoretical
Molecular Science Laboratory, RIKEN Cluster
for Pioneering Research, 2-1 Hirosawa, Wako, Saitama 351-0198, Japan
| | - Yuji Sugita
- Theoretical
Molecular Science Laboratory, RIKEN Cluster
for Pioneering Research, 2-1 Hirosawa, Wako, Saitama 351-0198, Japan
- Computational
Biophysics Research Team, RIKEN Center for
Computational Science, 7-1-26 minatojima-Minamimachi, Chuo-ku, Kobe, Hyogo 650-0047, Japan
- Laboratory
for Biomolecular Function Simulation, RIKEN
Center for Biosystems Dynamics Research, 1-6-5 minatojima-Minamimachi, Chuo-ku, Kobe, Hyogo 650-0047, Japan
| |
Collapse
|
24
|
Feng M, Song Y, Chen SH, Zhang Y, Zhou R. Molecular mechanism of secreted amyloid-β precursor protein in binding and modulating GABA BR1a. Chem Sci 2021; 12:6107-6116. [PMID: 33996007 PMCID: PMC8098695 DOI: 10.1039/d0sc06946a] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/31/2022] Open
Abstract
A recent phenomenal study discovered that the extension domain of secreted amyloid-β precursor protein (sAPP) can bind to the intrinsically disordered sushi 1 domain of the γ-aminobutyric acid type B receptor subunit 1a (GABABR1a) and modulate its synaptic transmission. The work provided an important structural foundation for the modulation of GABABR1a; however, the detailed molecular interaction mechanism, crucial for future drug design, remains elusive. Here, we further investigated the dynamical interactions between sAPP peptides and the natively unstructured sushi 1 domain using all-atom molecular dynamics simulations, for both the 17-residue sAPP peptide (APP 17-mer) and its minimally active 9 residue segment (APP 9-mer). We then explored mutations of the APP 9-mer with rigorous free energy perturbation (FEP) calculations. Our in silico mutagenesis studies revealed key residues (D4, W6, and W7) responsible for the binding with the sushi 1 domain. More importantly, one double mutation based on different vertebrate APP sequences from evolution exhibited a stronger binding (ΔΔG = −1.91 ± 0.66 kcal mol−1), indicating a potentially enhanced GABABR1a modulator. These large-scale simulations may provide new insights into the binding mechanism between sAPP and the sushi 1 domain, which could open new avenues in the development of future GABABR1a-specific therapeutics. A recent phenomenal study discovered that the extension domain of secreted amyloid-β precursor protein (sAPP) can bind to the intrinsically disordered sushi 1 domain of the γ-aminobutyric acid type B receptor subunit 1a (GABABR1a) and modulate its synaptic transmission.![]()
Collapse
Affiliation(s)
- Mei Feng
- Institute of Quantitative Biology, Shanghai Institute for Advanced Study, College of Life Sciences, Department of Physics, Zhejiang University Hangzhou 310027 China .,Lanzhou Center for Theoretical Physics, Key Laboratory of Theoretical Physics of Gansu Province, Lanzhou University Lanzhou Gansu 730000 China
| | - Yi Song
- Institute of Quantitative Biology, Shanghai Institute for Advanced Study, College of Life Sciences, Department of Physics, Zhejiang University Hangzhou 310027 China
| | - Serena H Chen
- Computational Sciences and Engineering Division, Oak Ridge National Laboratory Oak Ridge TN 37830 USA
| | - Yuanzhao Zhang
- Center for Applied Mathematics, Cornell University Ithaca NY 14583 USA
| | - Ruhong Zhou
- Institute of Quantitative Biology, Shanghai Institute for Advanced Study, College of Life Sciences, Department of Physics, Zhejiang University Hangzhou 310027 China .,Department of Chemistry, Columbia University New York NY 10027 USA
| |
Collapse
|
25
|
Heinzelmann G, Gilson MK. Automation of absolute protein-ligand binding free energy calculations for docking refinement and compound evaluation. Sci Rep 2021; 11:1116. [PMID: 33441879 PMCID: PMC7806944 DOI: 10.1038/s41598-020-80769-1] [Citation(s) in RCA: 40] [Impact Index Per Article: 13.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/20/2020] [Accepted: 12/24/2020] [Indexed: 02/06/2023] Open
Abstract
Absolute binding free energy calculations with explicit solvent molecular simulations can provide estimates of protein-ligand affinities, and thus reduce the time and costs needed to find new drug candidates. However, these calculations can be complex to implement and perform. Here, we introduce the software BAT.py, a Python tool that invokes the AMBER simulation package to automate the calculation of binding free energies for a protein with a series of ligands. The software supports the attach-pull-release (APR) and double decoupling (DD) binding free energy methods, as well as the simultaneous decoupling-recoupling (SDR) method, a variant of double decoupling that avoids numerical artifacts associated with charged ligands. We report encouraging initial test applications of this software both to re-rank docked poses and to estimate overall binding free energies. We also show that it is practical to carry out these calculations cheaply by using graphical processing units in common machines that can be built for this purpose. The combination of automation and low cost positions this procedure to be applied in a relatively high-throughput mode and thus stands to enable new applications in early-stage drug discovery.
Collapse
Affiliation(s)
- Germano Heinzelmann
- Departamento de Física, Universidade Federal de Santa Catarina, Florianópolis, Santa Catarina, Brazil.
| | - Michael K Gilson
- Skaggs School of Pharmacy and Pharmaceutical Sciences, University of California, San Diego, USA
| |
Collapse
|
26
|
Jespers W, Åqvist J, Gutiérrez-de-Terán H. Free Energy Calculations for Protein-Ligand Binding Prediction. Methods Mol Biol 2021; 2266:203-226. [PMID: 33759129 DOI: 10.1007/978-1-0716-1209-5_12] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/12/2023]
Abstract
Computational prediction of protein-ligand binding involves initial determination of the binding mode and subsequent evaluation of the strength of the protein-ligand interactions, which directly correlates with ligand binding affinities. As a consequence of increasing computer power, rigorous approaches to calculate protein-ligand binding affinities, such as free energy perturbation (FEP) methods, are becoming an essential part of the toolbox of computer-aided drug design. In this chapter, we provide a general overview of these methods and introduce the QFEP modules, which are open-source API workflows based on our molecular dynamics (MD) package Q. The module QligFEP allows estimation of relative binding affinities along ligand series, while QresFEP is a module to estimate binding affinity shifts caused by single-point mutations of the protein. We herein provide guidelines for the use of each of these modules based on data extracted from ligand-design projects. While these modules are stand-alone, the combined use of the two workflows in a drug-design project yields complementary perspectives of the ligand binding problem, providing two sides of the same coin. The selected case studies illustrate how to use QFEP to approach the two key questions associated with ligand binding prediction: identifying the most favorable binding mode from different alternatives and establishing structure-affinity relationships that allow the rational optimization of hit compounds.
Collapse
Affiliation(s)
- Willem Jespers
- Department of Cell and Molecular Biology, Biomedical Center, Uppsala University, Uppsala, Sweden
| | - Johan Åqvist
- Department of Cell and Molecular Biology, Biomedical Center, Uppsala University, Uppsala, Sweden
| | - Hugo Gutiérrez-de-Terán
- Department of Cell and Molecular Biology, Biomedical Center, Uppsala University, Uppsala, Sweden.
| |
Collapse
|
27
|
Ben-Shalom IY, Lin Z, Radak BK, Lin C, Sherman W, Gilson MK. Accounting for the Central Role of Interfacial Water in Protein-Ligand Binding Free Energy Calculations. J Chem Theory Comput 2020; 16:7883-7894. [PMID: 33206520 PMCID: PMC7725968 DOI: 10.1021/acs.jctc.0c00785] [Citation(s) in RCA: 13] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/07/2023]
Abstract
Rigorous binding free energy methods in drug discovery are growing in popularity because of a combination of methodological advances, improvements in computer hardware, and workflow automation. These calculations typically use molecular dynamics (MD) to sample from the Boltzmann distribution of conformational states. However, when part or all of the binding sites is inaccessible to the bulk solvent, the time needed for water molecules to equilibrate between bulk solvent and the binding site can be well beyond what is practical with standard MD. This sampling limitation is problematic in relative binding free energy calculations, which compute the reversible work of converting ligand 1 to ligand 2 within the binding site. Thus, if ligand 1 is smaller and/or more polar than ligand 2, the perturbation may allow additional water molecules to occupy a region of the binding site. However, this change in hydration may not be captured by standard MD simulations and may therefore lead to errors in the computed free energy. We recently developed a hybrid Monte Carlo/MD (MC/MD) method, which speeds up the equilibration of water between bulk solvent and buried cavities, while sampling from the intended distribution of states. Here, we report on the use of this approach in the context of alchemical binding free energy calculations. We find that using MC/MD markedly improves the accuracy of the calculations and also reduces hysteresis between the forward and reverse perturbations, relative to matched calculations using only MD with or without the crystallographic water molecules. The present method is available for use in AMBER simulation software.
Collapse
Affiliation(s)
- Ido Y Ben-Shalom
- Skaggs School of Pharmacy and Pharmaceutical Sciences, University of California, San Diego, 92093 La Jolla, California, United States
| | - Zhixiong Lin
- Silicon Therapeutics LLC, Boston, Massachusetts 02110, United States
| | - Brian K Radak
- Silicon Therapeutics LLC, Boston, Massachusetts 02110, United States
| | - Charles Lin
- Silicon Therapeutics LLC, Boston, Massachusetts 02110, United States
| | - Woody Sherman
- Silicon Therapeutics LLC, Boston, Massachusetts 02110, United States
| | - Michael K Gilson
- Skaggs School of Pharmacy and Pharmaceutical Sciences, University of California, San Diego, 92093 La Jolla, California, United States
| |
Collapse
|
28
|
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
|
29
|
Song LF, Merz KM. Evolution of Alchemical Free Energy Methods in Drug Discovery. J Chem Inf Model 2020; 60:5308-5318. [DOI: 10.1021/acs.jcim.0c00547] [Citation(s) in RCA: 30] [Impact Index Per Article: 7.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/20/2022]
Affiliation(s)
- Lin Frank Song
- Department of Chemistry and Department of Biochemistry and Molecular Biology, Michigan State University, 578 S. Shaw Lane, East Lansing, Michigan 48824, United States
| | - Kenneth M. Merz
- Department of Chemistry and Department of Biochemistry and Molecular Biology, Michigan State University, 578 S. Shaw Lane, East Lansing, Michigan 48824, United States
| |
Collapse
|
30
|
Perthold JW, Petrov D, Oostenbrink C. Toward Automated Free Energy Calculation with Accelerated Enveloping Distribution Sampling (A-EDS). J Chem Inf Model 2020; 60:5395-5406. [PMID: 32492343 PMCID: PMC7686955 DOI: 10.1021/acs.jcim.0c00456] [Citation(s) in RCA: 13] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/15/2022]
Abstract
![]()
Free-energy
perturbation (FEP) methods are commonly used in drug
design to calculate relative binding free energies of different ligands
to a common host protein. Alchemical ligand transformations are usually
performed in multiple steps which need to be chosen carefully to ensure
sufficient phase-space overlap between neighboring states. With one-step
or single-step FEP techniques, a single reference state is designed
that samples phase-space not only representative of a full transformation
but also ideally resembles multiple ligand end states and hence allows
for efficient multistate perturbations. Enveloping distribution sampling
(EDS) is one example for such a method in which the reference state
is created by a mathematical combination of the different ligand end
states based on solid statistical mechanics. We have recently proposed
a novel approach to EDS which enables efficient barrier crossing between
the different end states, termed accelerated EDS (A-EDS). In this
work, we further simplify the parametrization of the A-EDS reference
state and demonstrate the automated calculation of multiple free-energy
differences between different ligands from a single simulation in
three different well-described drug design model systems.
Collapse
Affiliation(s)
- Jan Walther Perthold
- Institute for Molecular Modeling and Simulation, Department for Material Sciences and Process Engineering, University of Natural Resources and Life Sciences (BOKU), Vienna, Muthgasse 18, 1190 Vienna, Austria
| | - Dražen Petrov
- Institute for Molecular Modeling and Simulation, Department for Material Sciences and Process Engineering, University of Natural Resources and Life Sciences (BOKU), Vienna, Muthgasse 18, 1190 Vienna, Austria
| | - Chris Oostenbrink
- Institute for Molecular Modeling and Simulation, Department for Material Sciences and Process Engineering, University of Natural Resources and Life Sciences (BOKU), Vienna, Muthgasse 18, 1190 Vienna, Austria
| |
Collapse
|
31
|
Hahn DF, Zarotiadis RA, Hünenberger PH. The Conveyor Belt Umbrella Sampling (CBUS) Scheme: Principle and Application to the Calculation of the Absolute Binding Free Energies of Alkali Cations to Crown Ethers. J Chem Theory Comput 2020; 16:2474-2493. [DOI: 10.1021/acs.jctc.9b00998] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- David F. Hahn
- Laboratory of Physical Chemistry, Department of Chemistry and Applied Biosciences, ETH Zürich, Vladimir-Prelog-Weg 2, 8093 Zürich, Switzerland
| | - Rhiannon A. Zarotiadis
- Laboratory of Physical Chemistry, Department of Chemistry and Applied Biosciences, ETH Zürich, Vladimir-Prelog-Weg 2, 8093 Zürich, Switzerland
| | - Philippe H. Hünenberger
- Laboratory of Physical Chemistry, Department of Chemistry and Applied Biosciences, ETH Zürich, Vladimir-Prelog-Weg 2, 8093 Zürich, Switzerland
| |
Collapse
|
32
|
Inoue M, Hayashi T, Hikiri S, Ikeguchi M, Kinoshita M. Hydration properties of a protein at low and high pressures: Physics of pressure denaturation. J Chem Phys 2020; 152:065103. [PMID: 32061219 DOI: 10.1063/1.5140499] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/30/2022] Open
Abstract
Using experimentally determined structures of ubiquitin at 1 and 3000 bar, we generate sufficiently large ensembles of model structures in the native and pressure-induced (denatured) states by means of molecular dynamics simulations with explicit water. We calculate the values of a free-energy function (FEF), which comprises the hydration free energy (HFE) and the intramolecular (conformational) energy and entropy, for the two states at 1 and 3000 bar. The HFE and the conformational entropy, respectively, are calculated using our statistical-mechanical method, which has recently been shown to be accurate, and the Boltzmann-quasi-harmonic method. The HFE is decomposed into a variety of physically insightful components. We show that the FEF of the native state is lower than that of the denatured state at 1 bar, whereas the opposite is true at 3000 bar, thus being successful in reproducing the pressure denaturation. We argue that the following two quantities of hydration play essential roles in the denaturation: the WASA-dependent term in the water-entropy loss upon cavity creation for accommodating the protein (WASA is the water-accessible surface area of the cavity) and the protein-water Lennard-Jones interaction energy. At a high pressure, the mitigation of the serious water crowding in the system is the most important, and the WASA needs to be sufficiently enlarged with the increase in the excluded-volume being kept as small as possible. The denatured structure thus induced is characterized by the water penetration into the protein interior. The pressure denaturation is accompanied by a significantly large gain of water entropy.
Collapse
Affiliation(s)
- Masao Inoue
- Institute of Advanced Energy, Kyoto University, Uji, Kyoto 611-0011, Japan
| | - Tomohiko Hayashi
- Institute of Advanced Energy, Kyoto University, Uji, Kyoto 611-0011, Japan
| | - Simon Hikiri
- Institute of Advanced Energy, Kyoto University, Uji, Kyoto 611-0011, Japan
| | - Mitsunori Ikeguchi
- Graduate School of Medical Life Science, Yokohama City University, 1-7-29, Suehiro-cho, Tsurumi-ku, Yokohama 230-0045, Japan
| | - Masahiro Kinoshita
- Institute of Advanced Energy, Kyoto University, Uji, Kyoto 611-0011, Japan
| |
Collapse
|
33
|
Hahn DF, König G, Hünenberger PH. Overcoming Orthogonal Barriers in Alchemical Free Energy Calculations: On the Relative Merits of λ-Variations, λ-Extrapolations, and Biasing. J Chem Theory Comput 2020; 16:1630-1645. [DOI: 10.1021/acs.jctc.9b00853] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/05/2023]
Affiliation(s)
- David F. Hahn
- Laboratory of Physical Chemistry, Department of Chemistry and Applied Biosciences, ETH Zürich, Vladimir-Prelog-Weg 2, 8093 Zürich, Switzerland
| | - Gerhard König
- Laboratory of Physical Chemistry, Department of Chemistry and Applied Biosciences, ETH Zürich, Vladimir-Prelog-Weg 2, 8093 Zürich, Switzerland
| | - Philippe H. Hünenberger
- Laboratory of Physical Chemistry, Department of Chemistry and Applied Biosciences, ETH Zürich, Vladimir-Prelog-Weg 2, 8093 Zürich, Switzerland
| |
Collapse
|
34
|
Öhlknecht C, Lier B, Petrov D, Fuchs J, Oostenbrink C. Correcting electrostatic artifacts due to net-charge changes in the calculation of ligand binding free energies. J Comput Chem 2020; 41:986-999. [PMID: 31930547 DOI: 10.1002/jcc.26143] [Citation(s) in RCA: 22] [Impact Index Per Article: 5.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/10/2019] [Revised: 12/16/2019] [Accepted: 12/22/2019] [Indexed: 01/06/2023]
Abstract
Alchemically derived free energies are artifacted when the perturbed moiety has a nonzero net charge. The source of the artifacts lies in the effective treatment of the electrostatic interactions within and between the perturbed atoms and remaining (partial) charges in the simulated system. To treat the electrostatic interactions effectively, lattice-summation (LS) methods or cutoff schemes in combination with a reaction-field contribution are usually employed. Both methods render the charging component of the calculated free energies sensitive to essential parameters of the system like the cutoff radius or the box side lengths. Here, we discuss the results of three previously published studies of ligand binding. These studies presented estimates of binding free energies that were artifacted due to the charged nature of the ligands. We show that the size of the artifacts can be efficiently calculated and raw simulation data can be corrected. We compare the corrected results with experimental estimates and nonartifacted estimates from path-sampling methods. Although the employed correction scheme involves computationally demanding continuum-electrostatics calculations, we show that the correction estimate can be deduced from a small sample of configurations rather than from the entire ensemble. This observation makes the calculations of correction terms feasible for complex biological systems. To show the general applicability of the proposed procedure, we also present results where the correction scheme was used to correct independent free energies obtained from simulations employing a cutoff scheme or LS electrostatics. In this work, we give practical guidelines on how to apply the appropriate corrections easily.
Collapse
Affiliation(s)
- Christoph Öhlknecht
- Institute of Molecular Modeling and Simulation, University of Natural Resources and Life Sciences, Vienna, Austria.,Austrian Centre of Industrial Biotechnology, Graz, Austria
| | - Bettina Lier
- Institute of Molecular Modeling and Simulation, University of Natural Resources and Life Sciences, Vienna, Austria
| | - Drazen Petrov
- Institute of Molecular Modeling and Simulation, University of Natural Resources and Life Sciences, Vienna, Austria
| | - Julian Fuchs
- Institute of Molecular Modeling and Simulation, University of Natural Resources and Life Sciences, Vienna, Austria.,Institute of General, Inorganic and Theoretical Chemistry, University of Innsbruck, Innsbruck, Austria
| | - Chris Oostenbrink
- Institute of Molecular Modeling and Simulation, University of Natural Resources and Life Sciences, Vienna, Austria
| |
Collapse
|
35
|
Cui G, Graves AP, Manas ES. GRAM: A True Null Model for Relative Binding Affinity Predictions. J Chem Inf Model 2020; 60:11-16. [PMID: 31874032 DOI: 10.1021/acs.jcim.9b00939] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Relative binding affinity prediction is a critical component in computer aided drug design. A significant amount of effort has been dedicated to developing rapid and reliable in silico methods. However, robust assessment of their performance is still a complicated issue, as it requires a performance measure applicable in the prospective setting and more importantly a true null model that defines the expected performance of being random in an objective manner. Although many performance metrics, such as the Pearson correlation coefficient (r), mean unsigned error (MUE), and root-mean-square error (RMSE), are frequently used in the literature, a true and nontrivial null model has yet been identified. To address this problem, here we introduce an interval estimate as an additional measure, namely, the prediction interval (PI), which can be estimated from the error distribution of the predictions. The benefits of using the interval estimate are (1) it provides the uncertainty range in the predicted activities, which is important in prospective applications, and (2) a true null model with well-defined PI can be established. We provide one such example termed the Gaussian Random Affinity Model (GRAM), which is based on the empirical observation that the affinity change in a typical lead optimization effort has the tendency to distribute normally N (0, σ). Having an analytically defined PI that only depends on the variation in the activities, GRAM should, in principle, allow us to compare the performance of relative binding affinity prediction methods in a standard way, ultimately critical to measuring the progress made in algorithm development.
Collapse
Affiliation(s)
- Guanglei Cui
- Computational and Modeling Science U.S., Platform Technology and Sciences , GlaxoSmithKline Pharmaceuticals , 1250 South Collegeville Road , Collegeville , Pennsylvania 19426 , United States
| | - Alan P Graves
- Computational and Modeling Science U.S., Platform Technology and Sciences , GlaxoSmithKline Pharmaceuticals , 1250 South Collegeville Road , Collegeville , Pennsylvania 19426 , United States
| | - Eric S Manas
- Computational and Modeling Science U.S., Platform Technology and Sciences , GlaxoSmithKline Pharmaceuticals , 1250 South Collegeville Road , Collegeville , Pennsylvania 19426 , United States
| |
Collapse
|
36
|
Abstract
When both the difference between two quantities and their individual values can be measured or computationally predicted, multiple quantities can be determined from the measurements or predictions of select individual quantities and select pairwise differences. These measurements and predictions form a network connecting the quantities through their differences. Here, I analyze the optimization of such networks, where the trace (A-optimal), the largest eigenvalue (E-optimal), or the determinant (D-optimal) of the covariance matrix associated with the estimated quantities are minimized with respect to the allocation of the measurement (or computational) cost to different measurements (or predictions). My statistical analysis of the performance of such optimal measurement networks-based on large sets of simulated data-suggests that they substantially accelerate the determination of the quantities and that they may be useful in applications such as the computational prediction of binding free energies of candidate drug molecules.
Collapse
Affiliation(s)
- Huafeng Xu
- Silicon Therapeutics , Boston , Massachusetts 02210 , United States
| |
Collapse
|
37
|
Lind C, Esguerra M, Jespers W, Satpati P, Gutierrez-de-Terán H, Åqvist J. Free energy calculations of RNA interactions. Methods 2019; 162-163:85-95. [DOI: 10.1016/j.ymeth.2019.02.014] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/03/2018] [Revised: 02/01/2019] [Accepted: 02/15/2019] [Indexed: 01/19/2023] Open
|
38
|
Erbaş A, Olvera de la Cruz M, Marko JF. Receptor-Ligand Rebinding Kinetics in Confinement. Biophys J 2019; 116:1609-1624. [PMID: 31029377 PMCID: PMC6506716 DOI: 10.1016/j.bpj.2019.02.033] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/16/2018] [Revised: 02/05/2019] [Accepted: 02/19/2019] [Indexed: 10/27/2022] Open
Abstract
Rebinding kinetics of molecular ligands plays a key role in the operation of biomachinery, from regulatory networks to protein transcription, and is also a key factor in design of drugs and high-precision biosensors. In this study, we investigate initial release and rebinding of ligands to their binding sites grafted on a planar surface, a situation commonly observed in single-molecule experiments and that occurs in vivo, e.g., during exocytosis. Via scaling arguments and molecular dynamic simulations, we analyze the dependence of nonequilibrium rebinding kinetics on two intrinsic length scales: the average separation distance between the binding sites and the total diffusible volume (i.e., height of the experimental reservoir in which diffusion takes place or average distance between receptor-bearing surfaces). We obtain time-dependent scaling laws for on rates and for the cumulative number of rebinding events. For diffusion-limited binding, the (rebinding) on rate decreases with time via multiple power-law regimes before the terminal steady-state (constant on-rate) regime. At intermediate times, when particle density has not yet become uniform throughout the diffusible volume, the cumulative number of rebindings exhibits a novel, to our knowledge, plateau behavior because of the three-dimensional escape process of ligands from binding sites. The duration of the plateau regime depends on the average separation distance between binding sites. After the three-dimensional diffusive escape process, a one-dimensional diffusive regime describes on rates. In the reaction-limited scenario, ligands with higher affinity to their binding sites (e.g., longer residence times) delay entry to the power-law regimes. Our results will be useful for extracting hidden timescales in experiments such as kinetic rate measurements for ligand-receptor interactions in microchannels, as well as for cell signaling via diffusing molecules.
Collapse
Affiliation(s)
- Aykut Erbaş
- UNAM-National Nanotechnology Research Center and Institute of Materials Science & Nanotechnology, Bilkent University, Ankara, Turkey.
| | - Monica Olvera de la Cruz
- Department of Materials Science and Engineering, Northwestern University, Evanston, Illinois; Department of Physics and Astronomy, Northwestern University, Evanston, Illinois; Department of Chemistry, Northwestern University, Evanston, Illinois
| | - John F Marko
- Department of Physics and Astronomy, Northwestern University, Evanston, Illinois; Department of Molecular Biosciences, Northwestern University, Evanston, Illinois.
| |
Collapse
|
39
|
Hikiri S, Hayashi T, Inoue M, Ekimoto T, Ikeguchi M, Kinoshita M. An accurate and rapid method for calculating hydration free energies of a variety of solutes including proteins. J Chem Phys 2019; 150:175101. [DOI: 10.1063/1.5093110] [Citation(s) in RCA: 16] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/06/2023] Open
Affiliation(s)
- Simon Hikiri
- Institute of Advanced Energy, Kyoto University, Uji, Kyoto 611-0011, Japan
- Graduate School of Science, Chiba University, 1-33 Yayoi-cho, Inage, Chiba 263-8522, Japan
| | - Tomohiko Hayashi
- Institute of Advanced Energy, Kyoto University, Uji, Kyoto 611-0011, Japan
| | - Masao Inoue
- Institute of Advanced Energy, Kyoto University, Uji, Kyoto 611-0011, Japan
| | - Toru Ekimoto
- Graduate School of Medical Life Science, Yokohama City University, 1-7-29, Suehiro-cho, Tsurumi-ku, Yokohama 230-0045, Japan
| | - Mitsunori Ikeguchi
- Graduate School of Medical Life Science, Yokohama City University, 1-7-29, Suehiro-cho, Tsurumi-ku, Yokohama 230-0045, Japan
- RIKEN Medical Sciences Innovation Hub Program, 1-7-22 Suehiro-cho, Tsurumi-ku, Yokohama 230-0045, Japan
| | - Masahiro Kinoshita
- Institute of Advanced Energy, Kyoto University, Uji, Kyoto 611-0011, Japan
| |
Collapse
|
40
|
Joshi DC, Lin J. Delineating Protein–Protein Curvilinear Dissociation Pathways and Energetics with Naïve Multiple‐Walker Umbrella Sampling Simulations. J Comput Chem 2019; 40:1652-1663. [DOI: 10.1002/jcc.25821] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/18/2018] [Revised: 01/03/2019] [Accepted: 02/10/2019] [Indexed: 01/01/2023]
Affiliation(s)
- Dhananjay C. Joshi
- Taiwan International Graduate Program (TIGP‐CBMB)Academia Sinica Taipei 11529 Taiwan
- Institute of Biochemical SciencesNational Taiwan University Taipei 10617 Taiwan
- Research Center for Applied SciencesAcademia Sinica Taipei 11529 Taiwan
| | - Jung‐Hsin Lin
- Research Center for Applied SciencesAcademia Sinica Taipei 11529 Taiwan
- Institute of Biomedical SciencesAcademia Sinica Taipei 11529 Taiwan
- College of Engineering SciencesChang Gung University Taoyuan 33302 Taiwan
- School of PharmacyNational Taiwan University Taipei 10050 Taiwan
| |
Collapse
|
41
|
Jespers W, Esguerra M, Åqvist J, Gutiérrez-de-Terán H. QligFEP: an automated workflow for small molecule free energy calculations in Q. J Cheminform 2019; 11:26. [PMID: 30941533 PMCID: PMC6444553 DOI: 10.1186/s13321-019-0348-5] [Citation(s) in RCA: 47] [Impact Index Per Article: 9.4] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/20/2018] [Accepted: 03/23/2019] [Indexed: 11/24/2022] Open
Abstract
The process of ligand binding to a biological target can be represented as the equilibrium between the relevant solvated and bound states of the ligand. This which is the basis of structure-based, rigorous methods such as the estimation of relative binding affinities by free energy perturbation (FEP). Despite the growing capacity of computing power and the development of more accurate force fields, a high throughput application of FEP is currently hampered due to the need, in the current schemes, of an expert user definition of the "alchemical" transformations between molecules in the series explored. Here, we present QligFEP, a solution to this problem using an automated workflow for FEP calculations based on a dual topology approach. In this scheme, the starting poses of each of the two ligands, for which the relative affinity is to be calculated, are explicitly present in the MD simulations associated with the (dual topology) FEP transformation, making the perturbation pathway between the two ligands univocal. We show that this generalized method can be applied to accurately estimate solvation free energies for amino acid sidechain mimics, as well as the binding affinity shifts due to the chemical changes typical of lead optimization processes. This is illustrated in a number of protein systems extracted from other FEP studies in the literature: inhibitors of CDK2 kinase and a series of A2A adenosine G protein-coupled receptor antagonists, where the results obtained with QligFEP are in excellent agreement with experimental data. In addition, our protocol allows for scaffold hopping perturbations to identify the binding affinities between different core scaffolds, which we illustrate with a series of Chk1 kinase inhibitors. QligFEP is implemented in the open-source MD package Q, and works with the most common family of force fields: OPLS, CHARMM and AMBER.
Collapse
Affiliation(s)
- Willem Jespers
- Department of Cell and Molecular Biology, Uppsala University, Uppsala, 75124 Sweden
| | - Mauricio Esguerra
- Department of Cell and Molecular Biology, Uppsala University, Uppsala, 75124 Sweden
| | - Johan Åqvist
- Department of Cell and Molecular Biology, Uppsala University, Uppsala, 75124 Sweden
| | - Hugo Gutiérrez-de-Terán
- Department of Cell and Molecular Biology, Uppsala University, Uppsala, 75124 Sweden
- Science for Life Laboratory, Uppsala University, Uppsala, Sweden
| |
Collapse
|
42
|
Wang M, Mei Y, Ryde U. Host-Guest Relative Binding Affinities at Density-Functional Theory Level from Semiempirical Molecular Dynamics Simulations. J Chem Theory Comput 2019; 15:2659-2671. [PMID: 30811192 DOI: 10.1021/acs.jctc.8b01280] [Citation(s) in RCA: 16] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/02/2023]
Abstract
Relative free energies for the binding of nine cyclic carboxylate ligands to the octa-acid deep-cavity host were calculated at the combined density-functional theory and molecular mechanics (DFT/MM) level of theory. The DFT calculations employed the BLYP functional and the 6-31G* basis set for the ligand. We employed free-energy perturbations (FEP) with the reference-potential approach and used molecular dynamics (MD) simulations with the semiempirical quantum mechanical (SQM) PM6-DH+ method for the ligand as an intermediate level between MM and DFT/MM to improve the convergence. Thus, the relative binding free energy of two ligands was first calculated at the MM level by an alchemical transformation from one ligand to another in both the bound and unbound states. Then, for each ligand the free-energy correction for going from the MM to the SQM/MM potentials was calculated using explicit SQM/MM MD simulations. Finally, the free-energy correction for going from the SQM/MM to the DFT/MM potentials was estimated with FEP without running any DFT/MM simulations. Instead, the free energy was calculated by single-step exponential averaging (ssEA) or employing the cumulant approximation to the second order (CA). The results show that CA converges much better than ssEA, and with 500-4500 DFT/MM single-point energy calculations, converged free energies with a precision of 0.3 kJ/mol can be obtained. These free energies reproduce the experimental binding free energy differences with a mean absolute deviation of 3.4 kJ/mol, a correlation ( R2) of 0.97, and correct signs for all of the eight free-energy differences. This is appreciably better than the results obtained at the SQM/MM level of theory and also slightly better than those obtained with MM. We show that the convergence of the SQM/MM → DFT/MM perturbations can be monitored by the use of Wu and Kofke's bias metric Π and by the standard deviation of the difference between the SQM/MM and DFT/MM energies. Finally, we show that the use of the intermediate SQM/MM MD simulations improves the convergence of the free energies by a factor of at least two, compared to doing direct MM → DFT/MM perturbations.
Collapse
Affiliation(s)
- Meiting Wang
- State Key Laboratory of Precision Spectroscopy, School of Physics and Materials Science , East China Normal University , Shanghai 200062 , China.,Department of Theoretical Chemistry , Lund University, Chemical Centre , P.O. Box 124, SE-221 00 Lund , Sweden
| | - Ye Mei
- State Key Laboratory of Precision Spectroscopy, School of Physics and Materials 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
| | - Ulf Ryde
- Department of Theoretical Chemistry , Lund University, Chemical Centre , P.O. Box 124, SE-221 00 Lund , Sweden
| |
Collapse
|
43
|
Ben-Shalom IY, Lin C, Kurtzman T, Walker RC, Gilson MK. Simulating Water Exchange to Buried Binding Sites. J Chem Theory Comput 2019; 15:2684-2691. [PMID: 30835999 DOI: 10.1021/acs.jctc.8b01284] [Citation(s) in RCA: 25] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/21/2022]
Abstract
Traditional molecular dynamics (MD) simulations of proteins, which relies on integration of Newton's equations of motion, cannot efficiently equilibrate water occupancy for buried cavities in proteins. This leads to slow convergence of thermodynamic averages for such systems. We have addressed this challenge by efficiently integrating standard Metropolis Monte Carlo (MC) translational water moves with MD in the AMBER simulation package. The translational moves allow water to easily enter or exit buried sites in a thermodynamically correct way during a simulation. To maximize efficiency, the algorithm avoids moves that only interchange waters within the bulk around the protein instead focusing on moves that can transfer water between bulk and the protein interior. In addition, a steric grid allows avoidance of moves that would lead to obvious steric clashes, and a fast grid-based energy evaluation is used to reduce the number of expensive full energy calculations. The potential energy distribution produced using MC/MD was found to be statistically indistinguishable from that of control simulations using only MD, and the algorithm effectively equilibrated water across steric barriers and into binding pockets that are not accessible with pure MD. The MC/MD method introduced here should be of increasing utility for applications spanning protein folding, the elucidation of protein mechanisms, and free energy calculations for computer-aided drug design. It is available in version 18 release of the widely disseminated AMBER simulation package.
Collapse
Affiliation(s)
- Ido Y Ben-Shalom
- Skaggs School of Pharmacy and Pharmaceutical Sciences , University of California, San Diego , La Jolla , California 92093 , United States
| | - Charles Lin
- Department of Chemistry and Biochemistry , University of California, San Diego , La Jolla , California 92093 , United States.,GlaxoSmithKline PLC, 1250 South Collegeville Road , Collegeville , Pennsylvania 19426 , United States
| | - Tom Kurtzman
- Department of Chemistry , Lehman College, The City University of New York , 250 Bedford Park Boulevard West , Bronx , New York 10468 , United States.,Ph.D. Programs in Biochemistry and Chemistry , The Graduate Center of The City University of New York , New York , New York 10016 , United States
| | - Ross C Walker
- Department of Chemistry and Biochemistry , University of California, San Diego , La Jolla , California 92093 , United States.,GlaxoSmithKline PLC, 1250 South Collegeville Road , Collegeville , Pennsylvania 19426 , United States
| | - Michael K Gilson
- Skaggs School of Pharmacy and Pharmaceutical Sciences , University of California, San Diego , La Jolla , California 92093 , United States
| |
Collapse
|
44
|
Chong S, Ham S. A New Computational Method for Protein–Ligand Binding Thermodynamics. B KOREAN CHEM SOC 2019. [DOI: 10.1002/bkcs.11681] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/12/2022]
Affiliation(s)
- Song‐Ho Chong
- Department of ChemistryThe Research Institute of Natural Sciences, Sookmyung Women's University, Cheongpa‐ro‐47‐gil 100, Yongsan‐ku Seoul 04310 South Korea
| | - Sihyun Ham
- Department of ChemistryThe Research Institute of Natural Sciences, Sookmyung Women's University, Cheongpa‐ro‐47‐gil 100, Yongsan‐ku Seoul 04310 South Korea
| |
Collapse
|
45
|
Feng M, Bell DR, Kang H, Shao Q, Zhou R. Exploration of HIV-1 fusion peptide–antibody VRC34.01 binding reveals fundamental neutralization sites. Phys Chem Chem Phys 2019; 21:18569-18576. [DOI: 10.1039/c9cp02909e] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/11/2022]
Abstract
VRC34.01 antibody binding to a vulnerable site of HIV envelope glycoprotein (Env), the gp41 fusion peptide, renders robust HIV neutralization, but several critical mutations decrease binding affinity and result in unbinding.
Collapse
Affiliation(s)
- Mei Feng
- Department of Physics
- Institute of Quantitative Biology
- Zhejiang University
- Hangzhou
- China
| | - David R. Bell
- Computational Biological Center
- IBM Thomas J. Watson Research Center
- Yorktown Heights
- USA
| | - Hongsuk Kang
- Computational Biological Center
- IBM Thomas J. Watson Research Center
- Yorktown Heights
- USA
| | - Qiwen Shao
- College of Nano Science and Technology
- Soochow University
- Suzhou 215123
- China
| | - Ruhong Zhou
- Department of Physics
- Institute of Quantitative Biology
- Zhejiang University
- Hangzhou
- China
| |
Collapse
|
46
|
Kilburg D, Gallicchio E. Analytical Model of the Free Energy of Alchemical Molecular Binding. J Chem Theory Comput 2018; 14:6183-6196. [DOI: 10.1021/acs.jctc.8b00967] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/16/2023]
Affiliation(s)
- Denise Kilburg
- Department of Chemistry, Brooklyn College of the City University of New York, Brooklyn, New York 11210, United States
| | - Emilio Gallicchio
- Department of Chemistry, Brooklyn College of the City University of New York, Brooklyn, New York 11210, United States
| |
Collapse
|
47
|
Loeffler HH, Bosisio S, Duarte Ramos Matos G, Suh D, Roux B, Mobley DL, Michel J. Reproducibility of Free Energy Calculations across Different Molecular Simulation Software Packages. J Chem Theory Comput 2018; 14:5567-5582. [PMID: 30289712 DOI: 10.1021/acs.jctc.8b00544] [Citation(s) in RCA: 57] [Impact Index Per Article: 9.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Alchemical free energy calculations are an increasingly important modern simulation technique to calculate free energy changes on binding or solvation. Contemporary molecular simulation software such as AMBER, CHARMM, GROMACS, and SOMD include support for the method. Implementation details vary among those codes, but users expect reliability and reproducibility, i.e., for a given molecular model and set of force field parameters, comparable free energy differences should be obtained within statistical bounds regardless of the code used. Relative alchemical free energy (RAFE) simulation is increasingly used to support molecule discovery projects, yet the reproducibility of the methodology has been less well tested than its absolute counterpart. Here we present RAFE calculations of hydration free energies for a set of small organic molecules and demonstrate that free energies can be reproduced to within about 0.2 kcal/mol with the aforementioned codes. Absolute alchemical free energy simulations have been carried out as a reference. Achieving this level of reproducibility requires considerable attention to detail and package-specific simulation protocols, and no universally applicable protocol emerges. The benchmarks and protocols reported here should be useful for the community to validate new and future versions of software for free energy calculations.
Collapse
Affiliation(s)
- Hannes H Loeffler
- Science & Technology Facilities Council , Daresbury, Warrington WA4 4AD , United Kingdom
| | - Stefano Bosisio
- EaStCHEM School of Chemistry , University of Edinburgh , David Brewster Road , Edinburgh EH9 3FJ , United Kingdom
| | | | - Donghyuk Suh
- University of Chicago , Chicago , Illinois 60637 , United States
| | - Benoit Roux
- University of Chicago , Chicago , Illinois 60637 , United States
| | - David L Mobley
- Departments of Pharmaceutical Sciences and Chemistry , University of California , Irvine , California 92697 , United States
| | - Julien Michel
- EaStCHEM School of Chemistry , University of Edinburgh , David Brewster Road , Edinburgh EH9 3FJ , United Kingdom
| |
Collapse
|
48
|
Baz J, Gebhardt J, Kraus H, Markthaler D, Hansen N. Insights into Noncovalent Binding Obtained from Molecular Dynamics Simulations. CHEM-ING-TECH 2018. [DOI: 10.1002/cite.201800050] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/11/2022]
Affiliation(s)
- Jörg Baz
- University of Stuttgart; Institute of Thermodynamics and Thermal Process Engineering; Pfaffenwaldring 9 70569 Stuttgart Germany
| | - Julia Gebhardt
- University of Stuttgart; Institute of Thermodynamics and Thermal Process Engineering; Pfaffenwaldring 9 70569 Stuttgart Germany
| | - Hamzeh Kraus
- University of Stuttgart; Institute of Thermodynamics and Thermal Process Engineering; Pfaffenwaldring 9 70569 Stuttgart Germany
| | - Daniel Markthaler
- University of Stuttgart; Institute of Thermodynamics and Thermal Process Engineering; Pfaffenwaldring 9 70569 Stuttgart Germany
| | - Niels Hansen
- University of Stuttgart; Institute of Thermodynamics and Thermal Process Engineering; Pfaffenwaldring 9 70569 Stuttgart Germany
| |
Collapse
|
49
|
Jackson NE, Webb MA, de Pablo JJ. Layered nested Markov chain Monte Carlo. J Chem Phys 2018; 149:072326. [PMID: 30134725 DOI: 10.1063/1.5030531] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
A configurational sampling algorithm based on nested layerings of Markov chains (Layered Nested Markov Chain Monte Carlo or L-NMCMC) is presented for simulations of systems characterized by rugged free energy landscapes. The layerings are generated using a set of auxiliary potential energy surfaces. The implementation of the method is demonstrated in the context of a rugged, two-dimensional potential energy surface. The versatility of the algorithm is next demonstrated on a simple, many-body system, namely, a canonical Lennard-Jones fluid in the liquid state. In that example, different layering schemes and auxiliary potentials are used, including variable cutoff distances and excluded-volume tempering. In addition to calculating a variety of properties of the system, it is also shown that L-NMCMC, when combined with a free-energy perturbation formalism, provides a straightforward means to construct approximate free-energy surfaces at no additional computational cost using the sampling distributions of each auxiliary Markov chain. The proposed L-NMCMC scheme is general in that it could be complementary to any number of methods that rely on sampling from a target distribution or methods that exploit a hierarchy of time scales and/or length scales through decomposition of the potential energy.
Collapse
Affiliation(s)
- Nicholas E Jackson
- Institute for Molecular Engineering, Argonne National Laboratory, Lemont, Illinois 06349, USA
| | - Michael A Webb
- Institute for Molecular Engineering, Argonne National Laboratory, Lemont, Illinois 06349, USA
| | - Juan J de Pablo
- Institute for Molecular Engineering, Argonne National Laboratory, Lemont, Illinois 06349, USA
| |
Collapse
|
50
|
Siebert DCB, Wieder M, Schlener L, Scholze P, Boresch S, Langer T, Schnürch M, Mihovilovic MD, Richter L, Ernst M, Ecker GF. SAR-Guided Scoring Function and Mutational Validation Reveal the Binding Mode of CGS-8216 at the α1+/γ2- Benzodiazepine Site. J Chem Inf Model 2018; 58:1682-1696. [PMID: 30028134 DOI: 10.1021/acs.jcim.8b00199] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/31/2022]
Abstract
The structural resolution of a bound ligand-receptor complex is a key asset to efficiently drive lead optimization in drug design. However, structural resolution of many drug targets still remains a challenging endeavor. In the absence of structural knowledge, scientists resort to structure-activity relationships (SARs) to promote compound development. In this study, we incorporated ligand-based knowledge to formulate a docking scoring function that evaluates binding poses for their agreement with a known SAR. We showcased this protocol by identifying the binding mode of the pyrazoloquinolinone (PQ) CGS-8216 at the benzodiazepine binding site of the GABAA receptor. Further evaluation of the final pose by molecular dynamics and free energy simulations revealed a close proximity between the pendent phenyl ring of the PQ and γ2D56, congruent with the low potency of carboxyphenyl analogues. Ultimately, we introduced the γ2D56A mutation and in fact observed a 10-fold potency increase in the carboxyphenyl analogue, providing experimental evidence in favor of our binding hypothesis.
Collapse
Affiliation(s)
- David C B Siebert
- Institute of Applied Synthetic Chemistry , TU Wien , Getreidemarkt 9/163 , 1060 Vienna , Austria
| | - Marcus Wieder
- Department of Pharmaceutical Chemistry , University of Vienna , Althanstrasse 14 , 1090 Vienna , Austria.,Faculty of Chemistry, Department of Computational Biological Chemistry , University of Vienna , Währingerstrasse 17 , 1090 Vienna , Austria
| | - Lydia Schlener
- Department of Pharmaceutical Chemistry , University of Vienna , Althanstrasse 14 , 1090 Vienna , Austria
| | - Petra Scholze
- Department of Pathobiology of the Nervous System, Center for Brain Research , Medical University of Vienna , Spitalgasse 4 , 1090 Vienna , Austria
| | - Stefan Boresch
- Faculty of Chemistry, Department of Computational Biological Chemistry , University of Vienna , Währingerstrasse 17 , 1090 Vienna , Austria
| | - Thierry Langer
- Department of Pharmaceutical Chemistry , University of Vienna , Althanstrasse 14 , 1090 Vienna , Austria
| | - Michael Schnürch
- Institute of Applied Synthetic Chemistry , TU Wien , Getreidemarkt 9/163 , 1060 Vienna , Austria
| | - Marko D Mihovilovic
- Institute of Applied Synthetic Chemistry , TU Wien , Getreidemarkt 9/163 , 1060 Vienna , Austria
| | - Lars Richter
- Department of Pharmaceutical Chemistry , University of Vienna , Althanstrasse 14 , 1090 Vienna , Austria
| | - Margot Ernst
- Department of Molecular Neurosciences, Center for Brain Research , Medical University of Vienna , Spitalgasse 4 , 1090 Vienna , Austria
| | - Gerhard F Ecker
- Department of Pharmaceutical Chemistry , University of Vienna , Althanstrasse 14 , 1090 Vienna , Austria
| |
Collapse
|