1
|
Majumder A, Kwon S, Straub JE. On Computing Equilibrium Binding Constants for Protein-Protein Association in Membranes. J Chem Theory Comput 2022; 18:3961-3971. [PMID: 35580264 PMCID: PMC11260346 DOI: 10.1021/acs.jctc.2c00106] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Protein association in lipid membranes is fundamental to membrane protein function and of great biomedical relevance. All-atom and coarse-grained models have been extensively used to understand the protein-protein interactions in the membrane and to compute equilibrium association constants. However, slow translational and rotational diffusion of protein in membrane presents challenges to the effective sampling of conformations defining the ensembles of free and bound states contributing to the association equilibrium and the free energy of dimerization. We revisit the homodimerization equilibrium of the TM region of glycophorin A. Conformational sampling is performed using umbrella sampling along previously proposed one-dimensional collective variables and compared with sampling over a two-dimensional collective variable space using the MARTINI v2.2 force field. We demonstrate that the one-dimensional collective variables suffer from restricted sampling of the native homodimer conformations leading to a biased free energy landscape. Conversely, simulations along the two-dimensional collective variable effectively characterize the thermodynamically relevant native and non-native interactions contributing to the association equilibrium. These results demonstrate the challenges associated with accurately characterizing binding equilibria when multiple poses contribute to the bound state ensemble.
Collapse
Affiliation(s)
- Ayan Majumder
- Department of Chemistry, Boston University, 590 Commonwealth Avenue, Boston, Massachusetts 02215, United States
| | - Seulki Kwon
- Department of Chemistry, Boston University, 590 Commonwealth Avenue, Boston, Massachusetts 02215, United States
| | - John E Straub
- Department of Chemistry, Boston University, 590 Commonwealth Avenue, Boston, Massachusetts 02215, United States
| |
Collapse
|
2
|
Gallicchio E. Free Energy-Based Computational Methods for the Study of Protein-Peptide Binding Equilibria. Methods Mol Biol 2022; 2405:303-334. [PMID: 35298820 DOI: 10.1007/978-1-0716-1855-4_15] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/14/2023]
Abstract
This chapter discusses the theory and application of physics-based free energy methods to estimate protein-peptide binding free energies. It presents a statistical mechanics formulation of molecular binding, which is then specialized in three methodologies: (1) alchemical absolute binding free energy estimation with implicit solvation, (2) alchemical relative binding free energy estimation with explicit solvation, and (3) potential of mean force binding free energy estimation. Case studies of protein-peptide binding application taken from the recent literature are discussed for each method.
Collapse
Affiliation(s)
- Emilio Gallicchio
- Department of Chemistry, Ph.D. Program in Biochemistry and Ph.D. Program in Chemistry at The Graduate Center of the City University of New York, Brooklyn College of the City University of New York, New York, NY, USA.
| |
Collapse
|
3
|
Duboué-Dijon E, Hénin J. Building intuition for binding free energy calculations: Bound state definition, restraints, and symmetry. J Chem Phys 2021; 154:204101. [PMID: 34241173 DOI: 10.1063/5.0046853] [Citation(s) in RCA: 20] [Impact Index Per Article: 6.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/15/2022] Open
Abstract
The theory behind computation of absolute binding free energies using explicit-solvent molecular simulations is well-established, yet somewhat complex, with counter-intuitive aspects. This leads to frequent frustration, common misconceptions, and sometimes erroneous numerical treatment. To improve this, we present the main practically relevant segments of the theory with constant reference to physical intuition. We pinpoint the role of the implicit or explicit definition of the bound state (or the binding site) to make a robust link between an experimental measurement and a computational result. We clarify the role of symmetry and discuss cases where symmetry number corrections have been misinterpreted. In particular, we argue that symmetry corrections as classically presented are a source of confusion and could be advantageously replaced by restraint free energy contributions. We establish that contrary to a common intuition, partial or missing sampling of some modes of symmetric bound states does not affect the calculated decoupling free energies. Finally, we review these questions and pitfalls in the context of a few common practical situations: binding to a symmetric receptor (equivalent binding sites), binding of a symmetric ligand (equivalent poses), and formation of a symmetric complex, in the case of homodimerization.
Collapse
Affiliation(s)
- E Duboué-Dijon
- CNRS, Université de Paris, UPR 9080, Laboratoire de Biochimie Théorique, 13 rue Pierre et Marie Curie, 75005 Paris, France
| | - J Hénin
- CNRS, Université de Paris, UPR 9080, Laboratoire de Biochimie Théorique, 13 rue Pierre et Marie Curie, 75005 Paris, France
| |
Collapse
|
4
|
Tahir U, Hussam A, Roy P, Hashmi I. Noncovalent Association and Partitioning of Some Perfume Components at Infinite Dilution with Myelin Basic Protein Pseudophase in Normal Saline. LANGMUIR : THE ACS JOURNAL OF SURFACES AND COLLOIDS 2021; 37:4793-4801. [PMID: 33851853 DOI: 10.1021/acs.langmuir.0c03475] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/12/2023]
Abstract
Myelin basic protein (MBP), one of the major protein constituents of the myelin sheath, possesses unique ligand-binding features. We present a novel equilibrium headspace gas chromatographic technique to examine the thermodynamics of noncovalent interactions between common perfume components: Lilial, Hedione, Hexylcinnamic aldehyde, and Versalide with MBP monomers and its hexameric MBP-pseudophase. A general theoretical model is used to calculate the critical aggregation concentration (cac) of MBP, perfume component binding constants with monomeric MBP, K11, and MBP as pseudophase, Kn1, and free energies for perfume component binding with monomeric MBP, ΔGb,11, and MBP as pseudophase, ΔGb,n1. In addition, the pseudophase-water partition coefficients, Kx, the free energies of transfer of perfume from bulk water to the MBP-pseudophase, ΔGt, and the intra-aggregate activity coefficients, γm∞, at infinite dilution were also determined. The cac value measured by the method of fractional distribution is a unique and precise approach in understanding the aggregation phenomenon. Within the experimental error, the 1:1 binding free energies did not differ by more than 1 kJ/mol among the perfume components but favored the MBP pseudophase binding by 6 kJ/mol. Therefore, that protein aggregation can enhance the binding of small molecules is probably a general conclusion. While the magnitudes of K11, Kn1, ΔGb, Kx, and ΔGt show weak trends, the γm∞ values show a strong and distinct trend in interaction, spanning 4 orders of magnitude among the perfume components.
Collapse
Affiliation(s)
| | | | | | - Irina Hashmi
- Department of Information Science and Technology, George Mason University, Fairfax, Virginia 22030, United States
| |
Collapse
|
5
|
Judy E, Kishore N. Discrepancies in Thermodynamic Information Obtained from Calorimetry and Spectroscopy in Ligand Binding Reactions: Implications on Correct Analysis in Systems of Biological Importance. BULLETIN OF THE CHEMICAL SOCIETY OF JAPAN 2021. [DOI: 10.1246/bcsj.20200248] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/09/2023]
Affiliation(s)
- Eva Judy
- Department of Chemistry, Indian Institute of Technology Bombay, Powai, Mumbai – 400 076, India
| | - Nand Kishore
- Department of Chemistry, Indian Institute of Technology Bombay, Powai, Mumbai – 400 076, India
| |
Collapse
|
6
|
Sen Chowdhury M, Sarkar A, Rai SR, Dasgupta S, Majumder I, Bhattacharya A, Das D, Bose D, Mukhopadhyay J, Mukhopadhyay M. Probing the binding interaction of zinc (II) Schiff bases with bovine serum albumin: A spectroscopic and molecular docking study. Appl Organomet Chem 2021. [DOI: 10.1002/aoc.6164] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/09/2022]
Affiliation(s)
- Megha Sen Chowdhury
- Department of Chemistry, Amity Institute of Applied Sciences (AIAS) Amity University Kolkata 700156 India
| | - Anwita Sarkar
- Amity Institute of Biotechnology Amity University Kolkata 700156 India
| | - Sristi Raj Rai
- Amity Institute of Biotechnology Amity University Kolkata 700156 India
| | - Sanchari Dasgupta
- Department of Chemistry University of Calcutta 92, A. P. C. Road Kolkata West Bengal 700 009 India
| | - Ishani Majumder
- Department of Chemistry University of Calcutta 92, A. P. C. Road Kolkata West Bengal 700 009 India
| | - Abir Bhattacharya
- Department of Physics, The Bhawanipur Education Society College University of Calcutta Kolkata 700020 India
| | - Debasis Das
- Department of Chemistry University of Calcutta 92, A. P. C. Road Kolkata West Bengal 700 009 India
| | - Debosreeta Bose
- Department of Chemistry, Amity Institute of Applied Sciences (AIAS) Amity University Kolkata 700156 India
| | - Jayanta Mukhopadhyay
- Energy Materials and Devices Division CSIR‐Central Glass and Ceramic Research Institute Kolkata 700 032 India
| | - Madhumita Mukhopadhyay
- Department of Chemistry, Amity Institute of Applied Sciences (AIAS) Amity University Kolkata 700156 India
| |
Collapse
|
7
|
Procacci P. Methodological uncertainties in drug-receptor binding free energy predictions based on classical molecular dynamics. Curr Opin Struct Biol 2020; 67:127-134. [PMID: 33220532 DOI: 10.1016/j.sbi.2020.08.001] [Citation(s) in RCA: 17] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/24/2020] [Revised: 08/02/2020] [Accepted: 08/02/2020] [Indexed: 12/13/2022]
Abstract
Computational approaches are becoming an essential tool in modern drug design and discovery, with fast compound triaging using a combination of machine learning and docking techniques followed by molecular dynamics binding free energies assessment using alchemical techniques. The traditional MD-based alchemical free energy perturbation (FEP) method faces severe sampling issues that may limits its reliability in automated workflows. Here we review the major sources of uncertainty in FEP protocols for drug discovery, showing how the sampling problem can be effectively tackled by switching to nonequilibrium alchemical techniques.
Collapse
Affiliation(s)
- Piero Procacci
- Dipartimento di Chimica "Ugo Schiff", Università degli Studi di Firenze, dVia della Lastruccia 3, 50019 Sesto Fiorentino, Italy.
| |
Collapse
|
8
|
Tanida Y, Matsuura A. Alchemical free energy calculations via metadynamics: Application to the theophylline-RNA aptamer complex. J Comput Chem 2020; 41:1804-1819. [PMID: 32449538 DOI: 10.1002/jcc.26221] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/01/2019] [Revised: 04/03/2020] [Accepted: 04/26/2020] [Indexed: 11/08/2022]
Abstract
We propose a computational workflow for robust and accurate prediction of both binding poses and their affinities at early stage in designing drug candidates. Small, rigid ligands with few intramolecular degrees of freedom, for example, fragment-like molecules, have multiple binding poses, even at a single binding site, and their affinities are often close to each other. We explore various structures of ligand binding to a target through metadynamics using a small number of collective variables, followed by reweighting to obtain the atomic coordinates. After identifying each binding pose by cluster analysis, we perform alchemical free energy calculations on each structure to obtain the overall value. We applied this protocol in computing free energy of binding for the theophylline-RNA aptamer complex. Of the six (meta)stable structures found, the most favorable binding structure is consistent with the structure obtained by NMR. The overall free energy of binding reproduces the experimental values very well.
Collapse
|
9
|
Limongelli V. Ligand binding free energy and kinetics calculation in 2020. WILEY INTERDISCIPLINARY REVIEWS-COMPUTATIONAL MOLECULAR SCIENCE 2020. [DOI: 10.1002/wcms.1455] [Citation(s) in RCA: 50] [Impact Index Per Article: 12.5] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Affiliation(s)
- Vittorio Limongelli
- Faculty of Biomedical Sciences, Institute of Computational Science – Center for Computational Medicine in Cardiology Università della Svizzera italiana (USI) Lugano Switzerland
- Department of Pharmacy University of Naples “Federico II” Naples Italy
| |
Collapse
|
10
|
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
|
11
|
Modern drug design: the implication of using artificial neuronal networks and multiple molecular dynamic simulations. J Comput Aided Mol Des 2017; 32:299-311. [PMID: 29134430 PMCID: PMC5767208 DOI: 10.1007/s10822-017-0085-7] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/23/2017] [Accepted: 11/08/2017] [Indexed: 01/31/2023]
Abstract
We report the implementation of molecular modeling approaches developed as a part of the 2016 Grand Challenge 2, the blinded competition of computer aided drug design technologies held by the D3R Drug Design Data Resource ( https://drugdesigndata.org/ ). The challenge was focused on the ligands of the farnesoid X receptor (FXR), a highly flexible nuclear receptor of the cholesterol derivative chenodeoxycholic acid. FXR is considered an important therapeutic target for metabolic, inflammatory, bowel and obesity related diseases (Expert Opin Drug Metab Toxicol 4:523-532, 2015), but in the context of this competition it is also interesting due to the significant ligand-induced conformational changes displayed by the protein. To deal with these conformational changes we employed multiple simulations of molecular dynamics (MD). Our MD-based protocols were top-ranked in estimating the free energy of binding of the ligands and FXR protein. Our approach was ranked second in the prediction of the binding poses where we also combined MD with molecular docking and artificial neural networks. Our approach showed mediocre results for high-throughput scoring of interactions.
Collapse
|
12
|
Lapelosa M. Free Energy of Binding and Mechanism of Interaction for the MEEVD-TPR2A Peptide-Protein Complex. J Chem Theory Comput 2017; 13:4514-4523. [PMID: 28723223 DOI: 10.1021/acs.jctc.7b00105] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/04/2023]
Abstract
The association between the MEEVD C-terminal peptide from the heat shock protein 90 (Hsp90) and tetratricopeptide repeat A (TPR2A) domain of the heat shock organizing protein (Hop) is a useful prototype to study the fundamental molecular details about the Hop-Hsp90 interaction. We study here the mechanism of binding/unbinding and compute the standard binding free energy and potential of mean force for the association of the MEEVD peptide to the TPR2A domain using the Adaptive Biasing Force (ABF) methodology. We observe conformational changes of the peptide and the protein receptor induced by binding. We measure the binding free energy of -8.4 kcal/mol, which is consistent with experimental estimates. The simulations achieve multiple unbinding and rebinding events along a consistent pathway connecting the binding site to solvent. The MEEVD peptide slowly dissociates disrupting the hydrogen bonds first, then tilting on the side while preserving the interaction with the side chain of residue Asp 5 of the peptide. After this initial displacement, the peptide completely dissociates and moves into the solvent. Rebinding of the MEEVD peptide from the solvent to the receptor binding site occurs slowly through the portal of entry. Unbinding and rebinding go through intermediate states characterized by the peptide interacting with a lateral helix, helix A1, of the receptor with mainly Asp 5, Val 4, and Glu 3 of the peptide. This newly discovered intermediate structure is characterized by numerous contacts with the receptor which lead to complete formation of the bound complex. The structure of the bound complex obtained after rebinding is structurally very similar to the crystal structure of the complex (0.48 Å root-mean square deviation). The residues Asp 5, Val 4, and Glu 3 adopt conformations and intermolecular contacts with excellent structural similarity with the native ones. Finally, the dissociation and reassociation of MEEVD induce hydration/dehydration transitions, which provide insights on the role of desolvation and solvation processes in protein-peptide binding.
Collapse
Affiliation(s)
- Mauro Lapelosa
- Department of Drug Discovery and Development, Italian Institute of Technology , Via Morego 30, Genova 16163, Italy
| |
Collapse
|
13
|
Alchemical determination of drug-receptor binding free energy: Where we stand and where we could move to. J Mol Graph Model 2017; 71:233-241. [DOI: 10.1016/j.jmgm.2016.11.018] [Citation(s) in RCA: 17] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/21/2016] [Revised: 11/24/2016] [Accepted: 11/29/2016] [Indexed: 01/05/2023]
|
14
|
Procacci P. Reformulating the entropic contribution in molecular docking scoring functions. J Comput Chem 2016; 37:1819-27. [DOI: 10.1002/jcc.24397] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/20/2016] [Revised: 03/30/2016] [Accepted: 04/10/2016] [Indexed: 11/08/2022]
|
15
|
Procacci P. I. Dissociation free energies of drug–receptor systems via non-equilibrium alchemical simulations: a theoretical framework. Phys Chem Chem Phys 2016; 18:14991-5004. [DOI: 10.1039/c5cp05519a] [Citation(s) in RCA: 25] [Impact Index Per Article: 3.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
In this contribution I critically discuss the alchemical approach for evaluating binding free energies in drug–receptor systems, placing this methodology into the broader context of non-equilibrium thermodynamics.
Collapse
|
16
|
Chang CE, Chen W, Gilson MK. Evaluating the Accuracy of the Quasiharmonic Approximation. J Chem Theory Comput 2015; 1:1017-28. [PMID: 26641917 DOI: 10.1021/ct0500904] [Citation(s) in RCA: 144] [Impact Index Per Article: 16.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
The quasiharmonic approximation (QH) allows the configurational entropy of a molecule to be estimated on the basis of a molecular dynamics simulation, through construction of a Gaussian probability distribution of conformations with variances equal to those provided by the simulation. At its introduction in 1981, the QH method was successfully applied to simple molecular systems with only one highly occupied energy well, and fluctuations were analyzed in a system of internal bond-angle-torsion coordinates. However, more recent studies have applied the QH method to complex biomolecular systems and have relied upon Cartesian coordinates. The present study evaluates the accuracy of the QH method through comparisons with more detailed methods. The chief findings are that the QH method can markedly overestimate the configurational entropy for systems with multiple occupied energy wells and that such errors tend to be magnified by the use of Cartesian coordinates instead of bond-angle-torsion coordinates.
Collapse
Affiliation(s)
- Chia-En Chang
- Department of Chemistry and Biochemistry, University of California at San Diego, La Jolla, California, and Center for Advanced Research in Biotechnology, 9600 Gudelsky Drive, Rockville, Maryland 20850
| | - Wei Chen
- Department of Chemistry and Biochemistry, University of California at San Diego, La Jolla, California, and Center for Advanced Research in Biotechnology, 9600 Gudelsky Drive, Rockville, Maryland 20850
| | - Michael K Gilson
- Department of Chemistry and Biochemistry, University of California at San Diego, La Jolla, California, and Center for Advanced Research in Biotechnology, 9600 Gudelsky Drive, Rockville, Maryland 20850
| |
Collapse
|
17
|
Deng Y, Roux B. Calculation of Standard Binding Free Energies: Aromatic Molecules in the T4 Lysozyme L99A Mutant. J Chem Theory Comput 2015; 2:1255-73. [PMID: 26626834 DOI: 10.1021/ct060037v] [Citation(s) in RCA: 224] [Impact Index Per Article: 24.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Calculations of the binding free energy of various nonpolar aromatic ligands with the L99A mutant of T4 lysozyme using molecular dynamics (MD) simulation are presented. To ensure better convergence, biasing potentials are used to restrain the ligand orientation and center-of-mass movement relative to the binding site when the ligand is decoupled from its environment in the binding pocket. The bias introduced by the restraint potentials is removed once the ligand fully interacts with the rest of the system and the calculated binding free energy is independent of the applied restraints. To decrease the computational cost, the simulations are generated with a reduced system in which protein and water atoms within a 15 Å-radius sphere around the ligand are included explicitly, while the rest of the system is treated with the generalized solvent boundary potential (GSBP). For all the ligands, the precision of the calculated free energy is less than 0.5 kcal/mol. For small nonpolar ligands such as benzene, toluene, and ethylbenzene, the calculated binding free energies are within 1.1 kcal/mol of the experimental values. For larger ligands, the computed binding free energies are slightly more favorable than the experimental values. The nonbinding polar molecule, phenol, has a calculated binding free energy of -0.88 kcal/mol. The simulation protocol presented here provides a way to calculate the binding free energy of small molecules to receptors at moderate computational cost.
Collapse
Affiliation(s)
- Yuqing Deng
- Biosciences Division, Argonne National Laboratory, 9700 South Cass Avenue, Building 221, Argonne, Illinois 60439, and Institute of Molecular Pediatric Sciences, Gordon Center for Integrative Sciences, University of Chicago, 929 East 57th Street, Room W323B, Chicago, Illinois 60637
| | - Benoît Roux
- Biosciences Division, Argonne National Laboratory, 9700 South Cass Avenue, Building 221, Argonne, Illinois 60439, and Institute of Molecular Pediatric Sciences, Gordon Center for Integrative Sciences, University of Chicago, 929 East 57th Street, Room W323B, Chicago, Illinois 60637
| |
Collapse
|
18
|
Hansen N, van Gunsteren WF. Practical Aspects of Free-Energy Calculations: A Review. J Chem Theory Comput 2014; 10:2632-47. [PMID: 26586503 DOI: 10.1021/ct500161f] [Citation(s) in RCA: 289] [Impact Index Per Article: 28.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/26/2022]
Abstract
Free-energy calculations in the framework of classical molecular dynamics simulations are nowadays used in a wide range of research areas including solvation thermodynamics, molecular recognition, and protein folding. The basic components of a free-energy calculation, that is, a suitable model Hamiltonian, a sampling protocol, and an estimator for the free energy, are independent of the specific application. However, the attention that one has to pay to these components depends considerably on the specific application. Here, we review six different areas of application and discuss the relative importance of the three main components to provide the reader with an organigram and to make nonexperts aware of the many pitfalls present in free energy calculations.
Collapse
Affiliation(s)
- Niels Hansen
- Institute of Thermodynamics and Thermal Process Engineering, University of Stuttgart , D-70569 Stuttgart, Germany.,Laboratory of Physical Chemistry, Swiss Federal Institute of Technology, ETH , CH-8093 Zürich, Switzerland
| | - Wilfred F van Gunsteren
- Laboratory of Physical Chemistry, Swiss Federal Institute of Technology, ETH , CH-8093 Zürich, Switzerland
| |
Collapse
|
19
|
Purisima EO, Hogues H. Protein-ligand binding free energies from exhaustive docking. J Phys Chem B 2012; 116:6872-9. [PMID: 22432509 DOI: 10.1021/jp212646s] [Citation(s) in RCA: 10] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
We explore the use of exhaustive docking as an alternative to Monte Carlo and molecular dynamics sampling for the direct integration of the partition function for protein-ligand binding. We enumerate feasible poses for the ligand and calculate the Boltzmann factor contribution of each pose to the partition function. From the partition function, the free energy, enthalpy, and entropy can be derived. All our calculations are done with a continuum solvation model that includes solving the Poisson equation. In contrast to Monte Carlo and molecular dynamics simulations, exhaustive docking avoids (within the limitations of a discrete sampling) the question of "Have we run long enough?" due to its deterministic complete enumeration of states. We tested the method on the T4 lysozyme L99A mutant, which has a nonpolar cavity that can accommodate a number of small molecules. We tested two electrostatic models. Model 1 used a solute dielectric of 2.25 for the complex apoprotein and free ligand and 78.5 for the solvent. Model 2 used a solute dielectric of 2.25 for the complex and apoprotein but 1.0 for the free ligand. For our test set of eight molecules, we obtain a reasonable correlation with a Pearson r(2) = 0.66 using model 1. The trend in binding affinity ranking is generally preserved with a Kendall τ = 0.64 and Spearman ρ = 0.83. With model 2, the correlation is improved with a Pearson r(2) = 0.83, Kendall τ = 0.93, and Spearman ρ = 0.98. This suggests that the energy function and sampling method adequately captured most of the thermodynamics of binding of the nonpolar ligands to T4 lysozyme L99A.
Collapse
Affiliation(s)
- Enrico O Purisima
- Biotechnology Research Institute, National Research Council of Canada, 6100 Royalmount Avenue, Montreal, Quebec H4P 2R2, Canada.
| | | |
Collapse
|
20
|
Blind prediction of host-guest binding affinities: a new SAMPL3 challenge. J Comput Aided Mol Des 2012; 26:475-87. [PMID: 22366955 DOI: 10.1007/s10822-012-9554-1] [Citation(s) in RCA: 100] [Impact Index Per Article: 8.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/03/2012] [Accepted: 02/13/2012] [Indexed: 10/28/2022]
Abstract
The computational prediction of protein-ligand binding affinities is of central interest in early-stage drug-discovery, and there is a widely recognized need for improved methods. Low molecular weight receptors and their ligands--i.e., host-guest systems--represent valuable test-beds for such affinity prediction methods, because their small size makes for fast calculations and relatively facile numerical convergence. The SAMPL3 community exercise included the first ever blind prediction challenge for host-guest binding affinities, through the incorporation of 11 new host-guest complexes. Ten participating research groups addressed this challenge with a variety of approaches. Statistical assessment indicates that, although most methods performed well at predicting some general trends in binding affinity, overall accuracy was not high, as all the methods suffered from either poor correlation or high RMS errors or both. There was no clear advantage in using explicit versus implicit solvent models, any particular force field, or any particular approach to conformational sampling. In a few cases, predictions using very similar energy models but different sampling and/or free-energy methods resulted in significantly different results. The protonation states of one host and some guest molecules emerged as key uncertainties beyond the choice of computational approach. The present results have implications for methods development and future blind prediction exercises.
Collapse
|
21
|
Gallicchio E. Role of Ligand Reorganization and Conformational Restraints on the Binding Free Energies of DAPY Non-Nucleoside Inhibitors to HIV Reverse Transcriptase. ACTA ACUST UNITED AC 2012; 2:7-22. [PMID: 22708073 DOI: 10.4236/cmb.2012.21002] [Citation(s) in RCA: 16] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/01/2023]
Abstract
The results of computer simulations of the binding of etravirine (TMC125) and rilpivirine (TMC278) to HIV reverse transcriptase are reported. It is confirmed that consistent binding free energy estimates are obtained with or without the application of torsional restraints when the free energies of imposing the restraints are taken into account. The restraints have a smaller influence on the thermodynamics and apparent kinetics of binding of TMC125 compared to the more flexible TMC278 inhibitor. The concept of the reorganization free energy of binding is useful to understand and categorize these effects. Contrary to expectations, the use of conformational restraints did not consistently enhance convergence of binding free energy estimates due to suppression of binding/unbinding pathways and due to the influence of rotational degrees of freedom not directly controlled by the restraints. Physical insights concerning the thermodynamic driving forces for binding and the role of "jiggling" and "wiggling" motion of the ligands are discussed. Based on these insights we conclude that an ideal inhibitor, if chemically realizable, would possess the electrostatic charge distribution of TMC125, so as to form strong interactions with the receptor, and the larger and more flexible substituents of TMC278, so as to minimize reorganization free energy penalties and the effects of resistance mutations, suitably modified, as in TMC125, so as to disfavor the formation of non-binding competent extended conformations when free in solution.
Collapse
Affiliation(s)
- Emilio Gallicchio
- BioMaPS Institute for Quantitative Biology and Department of Chemistry and Chemical Biology, Rutgers the State University of New Jersey, Piseatawav, NJ 08854
| |
Collapse
|
22
|
Lapelosa M, Gallicchio E, Levy RM. Conformational Transitions and Convergence of Absolute Binding Free Energy Calculations. J Chem Theory Comput 2011; 8:47-60. [PMID: 22368530 DOI: 10.1021/ct200684b] [Citation(s) in RCA: 55] [Impact Index Per Article: 4.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/03/2023]
Abstract
The Binding Energy Distribution Analysis Method (BEDAM) is employed to compute the standard binding free energies of a series of ligands to a FK506 binding protein (FKBP12) with implicit solvation. Binding free energy estimates are in reasonably good agreement with experimental affinities. The conformations of the complexes identified by the simulations are in good agreement with crystallographic data, which was not used to restrain ligand orientations. The BEDAM method is based on λ -hopping Hamiltonian parallel Replica Exchange (HREM) molecular dynamics conformational sampling, the OPLS-AA/AGBNP2 effective potential, and multi-state free energy estimators (MBAR). Achieving converged and accurate results depends on all of these elements of the calculation. Convergence of the binding free energy is tied to the level of convergence of binding energy distributions at critical intermediate states where bound and unbound states are at equilibrium, and where the rate of binding/unbinding conformational transitions is maximal. This finding mirrors similar observations in the context of order/disorder transitions as for example in protein folding. Insights concerning the physical mechanism of ligand binding and unbinding are obtained. Convergence for the largest FK506 ligand is achieved only after imposing strict conformational restraints, which however require accurate prior structural knowledge of the structure of the complex. The analytical AGBNP2 model is found to underestimate the magnitude of the hydrophobic driving force towards binding in these systems characterized by loosely packed protein-ligand binding interfaces. Rescoring of the binding energies using a numerical surface area model corrects this deficiency. This study illustrates the complex interplay between energy models, exploration of conformational space, and free energy estimators needed to obtain robust estimates from binding free energy calculations.
Collapse
Affiliation(s)
- Mauro Lapelosa
- BioMaPS Institute for Quantitative Biology and Department of Chemistry and Chemical Biology, Rutgers the State University of New Jersey, Piscataway, NJ 08854
| | | | | |
Collapse
|
23
|
Odinokov AV, Basilevsky MV, Petrov NK. Decay kinetics of the excited S1 state of the cyanine dye Cy+I− (thiacarbocyanine iodide): The computation of quantum yields for different pathways. J Chem Phys 2011; 135:144504. [DOI: 10.1063/1.3647956] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/17/2023] Open
|
24
|
Gallicchio E, Levy RM. Advances in all atom sampling methods for modeling protein-ligand binding affinities. Curr Opin Struct Biol 2011; 21:161-6. [PMID: 21339062 DOI: 10.1016/j.sbi.2011.01.010] [Citation(s) in RCA: 88] [Impact Index Per Article: 6.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/15/2010] [Revised: 01/25/2011] [Accepted: 01/27/2011] [Indexed: 01/23/2023]
Abstract
Conformational dynamics plays a fundamental role in the regulation of molecular recognition processes. Conformational heterogeneity and entropy variations upon binding, although not always evident from the analysis of structural data, can substantially affect affinity and specificity. Computer modeling is able to provide some of the most direct insights into these aspects of molecular recognition. We review recent physics-based computational studies that employ advanced conformational sampling algorithms and effective potentials to model the three main classes of degrees of freedom relevant to the binding process: ligand positioning relative to the receptor, ligand and receptor internal reorganization, and hydration. Collectively these studies show that all of these elements are important for proper modeling of protein-ligand interactions.
Collapse
Affiliation(s)
- Emilio Gallicchio
- BioMaPS Institute for Quantitative Biology and Department of Chemistry and Chemical Biology, Rutgers University, Piscataway, NJ 08854, USA
| | | |
Collapse
|
25
|
Gallicchio E, Levy RM. Recent theoretical and computational advances for modeling protein-ligand binding affinities. ADVANCES IN PROTEIN CHEMISTRY AND STRUCTURAL BIOLOGY 2011; 85:27-80. [PMID: 21920321 DOI: 10.1016/b978-0-12-386485-7.00002-8] [Citation(s) in RCA: 93] [Impact Index Per Article: 7.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Abstract
We review recent theoretical and algorithmic advances for the modeling of protein ligand binding free energies. We first describe a statistical mechanics theory of noncovalent association, with particular focus on deriving the fundamental formulas on which computational methods are based. The second part reviews the main computational models and algorithms in current use or development, pointing out the relations with each other and with the theory developed in the first part. Particular emphasis is given to the modeling of conformational reorganization and entropic effect. The methods reviewed are free energy perturbation, double decoupling, the Binding Energy Distribution Analysis Method, the potential of mean force method, mining minima and MM/PBSA. These models have different features and limitations, and their ranges of applicability vary correspondingly. Yet their origins can all be traced back to a single fundamental theory.
Collapse
Affiliation(s)
- Emilio Gallicchio
- Department of Chemistry and Chemical Biology, BioMaPS Institute for Quantitative Biology, Rutgers University, Piscataway, New Jersey, USA
| | | |
Collapse
|
26
|
Janosi L, Prakash A, Doxastakis M. Lipid-modulated sequence-specific association of glycophorin A in membranes. Biophys J 2010; 99:284-92. [PMID: 20655857 DOI: 10.1016/j.bpj.2010.04.005] [Citation(s) in RCA: 49] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/08/2010] [Revised: 04/01/2010] [Accepted: 04/05/2010] [Indexed: 11/26/2022] Open
Abstract
Protein association in lipid membranes is a complex process with thermodynamics directed by a multitude of different factors. Amino-acid sequence is a molecular parameter that affects dimerization as shown by limited directed mutations along the transmembrane domains. Membrane-mediated interactions are also important although details of such contributions remain largely unclear. In this study, we probe directly the free energy of association of Glycophorin A by means of extensive parallel Monte Carlo simulations with recently developed methods and a model that accounts for sequence-specificity while representing lipid membranes faithfully. We find that lipid-induced interactions are significant both at short and intermediate separations. The ability of molecules to tilt in a specific hydrophobic environment extends their accessible interfaces, leading to intermittent contacts during protein recognition. The dimer with the lowest free energy is largely determined by the favorable lipid-induced attractive interactions at the closest distance. Finally, the coarse-grained model employed herein, together with the extensive sampling performed, provides estimates of the free energy of association that are in excellent agreement with existing data.
Collapse
Affiliation(s)
- Lorant Janosi
- Department of Chemical and Biomolecular Engineering, University of Houston, Houston, Texas, USA
| | | | | |
Collapse
|
27
|
Gallicchio E, Lapelosa M, Levy RM. The Binding Energy Distribution Analysis Method (BEDAM) for the Estimation of Protein-Ligand Binding Affinities. J Chem Theory Comput 2010; 6:2961-2977. [PMID: 21116484 DOI: 10.1021/ct1002913] [Citation(s) in RCA: 125] [Impact Index Per Article: 8.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/26/2023]
Abstract
The Binding Energy Distribution Analysis Method (BEDAM) for the computation of receptor-ligand standard binding free energies with implicit solvation is presented. The method is based on a well established statistical mechanics theory of molecular association. It is shown that, in the context of implicit solvation, the theory is homologous to the test particle method of solvation thermodynamics with the solute-solvent potential represented by the effective binding energy of the protein-ligand complex. Accordingly, in BEDAM the binding constant is computed by means of a weighted integral of the probability distribution of the binding energy obtained in the canonical ensemble in which the ligand is positioned in the binding site but the receptor and the ligand interact only with the solvent continuum. It is shown that the binding energy distribution encodes all of the physical effects of binding. The balance between binding enthalpy and entropy is seen in our formalism as a balance between favorable and unfavorable binding modes which are coupled through the normalization of the binding energy distribution function. An efficient computational protocol for the binding energy distribution based on the AGBNP2 implicit solvent model, parallel Hamiltonian replica exchange sampling and histogram reweighting is developed. Applications of the method to a set of known binders and non-binders of the L99A and L99A/M102Q mutants of T4 lysozyme receptor are illustrated. The method is able to discriminate without error binders from non-binders, and the computed standard binding free energies of the binders are found to be in good agreement with experimental measurements. Analysis of the results reveals that the binding affinities of these systems reflect the contributions from multiple conformations spanning a wide range of binding energies.
Collapse
Affiliation(s)
- Emilio Gallicchio
- BioMaPS Institute for Quantitative Biology and Department of Chemistry and Chemical Biology, Rutgers the State University of New Jersey, Piscataway, NJ 08854
| | | | | |
Collapse
|
28
|
Thalmann F. A schematic model for molecular affinity and binding with Ising variables. THE EUROPEAN PHYSICAL JOURNAL. E, SOFT MATTER 2010; 31:441-454. [PMID: 20480964 DOI: 10.1140/epje/i2010-10600-9] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/29/2023]
Abstract
After discussing the relevance of statistical physics in molecular recognition processes, we present a schematic model for ligand-receptor association based on an Ising chain. We discuss the possible behaviors of the affinity when the stiffness of the ligand increases. We also consider the case of flexible receptors. A variety of interesting behaviors is obtained, including some affinity modulation upon bond hardening or softening. The affinity of a ligand for its receptor is shown to depend on the details of its rigidity profile, and we question the possibility of encoding information in the rigidities as well as in the shape. An exhaustive study of the selectivity of patterns with length n < 8 is carried out. Connection with other spin models, in particular spin glasses is mentioned in the conclusion.
Collapse
Affiliation(s)
- F Thalmann
- Institut Charles Sadron, Université de Strasbourg, CNRS UPR 22, 23 rue du Loess, BP 84047, F-67034 Strasbourg Cedex, France.
| |
Collapse
|
29
|
Affiliation(s)
- Huan-Xiang Zhou
- Department of Physics and Institute of Molecular Biophysics, Florida State University, Tallahassee, Florida 32306, USA.
| | | |
Collapse
|
30
|
Irudayam SJ, Henchman RH. Entropic cost of protein-ligand binding and its dependence on the entropy in solution. J Phys Chem B 2009; 113:5871-84. [PMID: 19351118 DOI: 10.1021/jp809968p] [Citation(s) in RCA: 53] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Two theoretical formulations are proposed and compared for the loss of translational and rotational entropy upon protein-ligand binding in water. The two theories share the same approach to evaluate the translational and rotational entropy of the ligand when bound. The potential of the bound ligand is modeled by six harmonic oscillators that are parametrized from the force and torque magnitudes measured in a molecular dynamics simulation, yielding vibrational and librational entropies. In the aqueous phase, the theories differ because there is no unique way to assign the total entropy to molecules in solution. In one approach, the ligand is allowed unrestricted access to the full solution volume at the standard concentration and is assigned the same translational and rotational entropy as if it were an ideal gas. We term this a "molecule-frame" (MF) theory because it considers configurational space in the reference frame of the molecule of interest. The entropy of the solvent is penalized because it is excluded from the molecule's volume. In the second theory, all molecules including the solvent are confined by their neighbors in mean-field configurational volumes. This we term a "system-frame" (SF) theory because the configurational space available to all molecules is considered in the reference frame of the whole system. Molecules have vibrational and librational entropy in the same way as they do when bound. In addition, the discrete size of the solvent molecules quantizes the configurational space into an effective number of minima according to the solute molecule's standard concentration and the mean volume of a solvent molecule. This leads to the cratic entropy expressed in terms of the solute molecule's mole fraction. The equivalent number of minima in rotational space depends on both the solute molecule's volume and the solvent molecule's volume. This leads to an equation for the orientational entropy based on the proposed concept of "angle fraction". The MF and SF theories are applied to calculate the translational and rotational entropy losses involved in the formation of six different protein-ligand complexes, in two of which the ligand is water. The MF entropy losses range from -80 to -142 J K(-1) mol(-1) for ligands at the 1 M standard-state concentration and from -52 to -63 J K(-1) mol(-1) for water at the 55.6 M standard-state concentration. They depend logarithmically on both the number and strength of interactions between the ligand and protein through the forces and torques. This is observed to lead to moderate dependencies on the ligands' moments of inertia and masses. The SF entropy losses are smaller and range from -50 to -75 J K(-1) mol(-1) for ligands at the 1 M standard-state concentration and from 0 to -12 J K(-1) mol(-1) for water. They depend logarithmically on the ligand solvent's molecular volume and weakly on the relative strengths of the ligand's interactions with the protein and water. The cratic entropy loss in water at the standard concentration is constant and is also demonstrated to be implicit in MF theories. Entropy losses from the two approaches are also compared with those from other computational approaches and with experiment. The use of the force and torque magnitudes leads to smaller bound volumes than are obtained from ligand-displacement approaches. The general agreement of the SF entropy losses with those from experiment suggests that the SF theory is more consistent with the assumptions made in experimental measurements than the MF solvation theories, which would require a compensating entropy gain in the solvent in order to agree.
Collapse
Affiliation(s)
- Sheeba Jem Irudayam
- Manchester Interdisciplinary Biocentre, The University of Manchester, 131 Princess Street, Manchester M1 7DN, United Kingdom
| | | |
Collapse
|
31
|
Doudou S, Burton NA, Henchman RH. Standard Free Energy of Binding from a One-Dimensional Potential of Mean Force. J Chem Theory Comput 2009; 5:909-18. [PMID: 26609600 DOI: 10.1021/ct8002354] [Citation(s) in RCA: 152] [Impact Index Per Article: 10.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
A practical approach that enables one to calculate the standard free energy of binding from a one-dimensional potential of mean force (PMF) is proposed. Umbrella sampling and the weighted histogram analysis method are used to generate a PMF along the reaction coordinate of binding. At each point, a restraint is applied orthogonal to the reaction coordinate to make possible the determination of the volume sampled by the ligand. The free energy of binding from an arbitrary unbound volume to the restrained bound form is calculated from the ratio of the PMF integrated over the bound region to that of the unbound. Adding the free energy changes from the standard-state volume to the unbound volume and from the restrained to the unrestrained bound state gives the standard free energy of binding. Exploration of the best choice of binding paths is also made. This approach is first demonstrated on a model binding system and then tested on the benzamidine-trypsin system for which reasonable agreement with experiment is found. A comparison is made with other methods to obtain the standard free energy of binding from the PMF.
Collapse
Affiliation(s)
- Slimane Doudou
- School of Chemistry, The University of Manchester, Oxford Road, Manchester, M13 9PL, United Kingdom, and Manchester Interdisciplinary Biocentre, The University of Manchester, 131 Princess Street, Manchester M1 7DN, United Kingdom
| | - Neil A Burton
- School of Chemistry, The University of Manchester, Oxford Road, Manchester, M13 9PL, United Kingdom, and Manchester Interdisciplinary Biocentre, The University of Manchester, 131 Princess Street, Manchester M1 7DN, United Kingdom
| | - Richard H Henchman
- School of Chemistry, The University of Manchester, Oxford Road, Manchester, M13 9PL, United Kingdom, and Manchester Interdisciplinary Biocentre, The University of Manchester, 131 Princess Street, Manchester M1 7DN, United Kingdom
| |
Collapse
|
32
|
Deng Y, Roux B. Computations of standard binding free energies with molecular dynamics simulations. J Phys Chem B 2009; 113:2234-46. [PMID: 19146384 PMCID: PMC3837708 DOI: 10.1021/jp807701h] [Citation(s) in RCA: 404] [Impact Index Per Article: 26.9] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
An increasing number of studies have reported computations of the standard (absolute) binding free energy of small ligands to proteins using molecular dynamics (MD) simulations and explicit solvent molecules that are in good agreement with experiments. This encouraging progress suggests that physics-based approaches hold the promise of making important contributions to the process of drug discovery and optimization in the near future. Two types of approaches are principally used to compute binding free energies with MD simulations. The most widely known is the alchemical double decoupling method, in which the interaction of the ligand with its surroundings are progressively switched off. It is also possible to use a potential of mean force (PMF) method, in which the ligand is physically separated from the protein receptor. For both of these computational approaches, restraining potentials may be activated and released during the simulation for sampling efficiently the changes in translational, rotational, and conformational freedom of the ligand and protein upon binding. Because such restraining potentials add bias to the simulations, it is important that their effects be rigorously removed to yield a binding free energy that is properly unbiased with respect to the standard state. A review of recent results is presented, and differences in computational methods are discussed. Examples of computations with T4-lysozyme mutants, FKBP12, SH2 domain, and cytochrome P450 are discussed and compared. Remaining difficulties and challenges are highlighted.
Collapse
Affiliation(s)
- Yuqing Deng
- Biosciences Division, Argonne National Laboratory, 9700 South Cass Avenue, Argonne, Illinois, USA
| | | |
Collapse
|
33
|
Lee MS, Olson MA. Calculation of absolute ligand binding free energy to a ribosome-targeting protein as a function of solvent model. J Phys Chem B 2008; 112:13411-7. [PMID: 18821791 DOI: 10.1021/jp802460p] [Citation(s) in RCA: 11] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
A comparative analysis is provided of the effect of different solvent models on the calculation of a potential of mean force (PMF) for determining the absolute binding affinity of the small molecule inhibitor pteroic acid bound to ricin toxin A-chain (RTA). Solvent models include the distance-dependent dielectric constant, several different generalized Born (GB) approximations, and a hybrid explicit/GB-based implicit solvent model. We found that the simpler approximation of dielectric screening and a GB model, with Born radii fitted to a switching-window dielectric-boundary surface Poisson solvent model, severely overpredicted the binding affinity as compared to the experimental value, estimated to range from -4.4 to -6.0 kcal/mol. In contrast, GB models that are parametrized to fit the Lee-Richards molecular surface performed much better, predicting binding free energy within 1-3 kcal/mol of experimental estimates. However, the predicted free-energy profiles of these GB models displayed alternative binding modes not observed in the crystal structure. Finally, the most rigorous and computationally costly approach in this work, which used a hybrid explicit/implicit solvent model, correctly determined a binding funnel in the PMF near the crystallographic bound state and predicted an absolute binding affinity that was 2 kcal/mol more favorable than the estimated experimental binding affinity.
Collapse
Affiliation(s)
- Michael S Lee
- Computational Sciences and Engineering Branch, U.S. Army Research Laboratory, Aberdeen Proving Ground, Maryland 21005, USA.
| | | |
Collapse
|
34
|
Noskov SY. Molecular mechanism of substrate specificity in the bacterial neutral amino acid transporter LeuT. Proteins 2008; 73:851-63. [DOI: 10.1002/prot.22108] [Citation(s) in RCA: 32] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/05/2022]
|
35
|
Gitlin I, Carbeck JD, Whitesides GM. Why are proteins charged? Networks of charge-charge interactions in proteins measured by charge ladders and capillary electrophoresis. Angew Chem Int Ed Engl 2007; 45:3022-60. [PMID: 16619322 DOI: 10.1002/anie.200502530] [Citation(s) in RCA: 196] [Impact Index Per Article: 11.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/07/2022]
Abstract
Almost all proteins contain charged amino acids. While the function in catalysis or binding of individual charges in the active site can often be identified, it is less clear how to assign function to charges beyond this region. Are they necessary for solubility? For reasons other than solubility? Can manipulating these charges change the properties of proteins? A combination of capillary electrophoresis (CE) and protein charge ladders makes it possible to study the roles of charged residues on the surface of proteins outside the active site. This method involves chemical modification of those residues to generate a large number of derivatives of the protein that differ in charge. CE separates those derivatives into groups with the same number of modified charged groups. By studying the influence of charge on the properties of proteins using charge ladders, it is possible to estimate the net charge and hydrodynamic radius and to infer the role of charged residues in ligand binding and protein folding.
Collapse
Affiliation(s)
- Irina Gitlin
- Department of Chemistry and Chemical Biology, Harvard University, 12 Oxford St., Cambridge, MA 02138, USA
| | | | | |
Collapse
|
36
|
Abstract
Accurate methods of computing the affinity of a small molecule with a protein are needed to speed the discovery of new medications and biological probes. This paper reviews physics-based models of binding, beginning with a summary of the changes in potential energy, solvation energy, and configurational entropy that influence affinity, and a theoretical overview to frame the discussion of specific computational approaches. Important advances are reported in modeling protein-ligand energetics, such as the incorporation of electronic polarization and the use of quantum mechanical methods. Recent calculations suggest that changes in configurational entropy strongly oppose binding and must be included if accurate affinities are to be obtained. The linear interaction energy (LIE) and molecular mechanics Poisson-Boltzmann surface area (MM-PBSA) methods are analyzed, as are free energy pathway methods, which show promise and may be ready for more extensive testing. Ultimately, major improvements in modeling accuracy will likely require advances on multiple fronts, as well as continued validation against experiment.
Collapse
Affiliation(s)
- Michael K Gilson
- Center for Advanced Research in Biotechnology, University of Maryland Biotechnology Institute, Rockville, Maryland 20850, USA.
| | | |
Collapse
|
37
|
Chang CEA, Trylska J, Tozzini V, McCammon JA. Binding pathways of ligands to HIV-1 protease: coarse-grained and atomistic simulations. Chem Biol Drug Des 2007; 69:5-13. [PMID: 17313452 DOI: 10.1111/j.1747-0285.2007.00464.x] [Citation(s) in RCA: 59] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Multiscale simulations (coarse-grained Brownian dynamics simulations and all-atom molecular dynamics simulations in implicit solvent) were applied to reveal the binding processes of ligands as they enter the binding site of the HIV-1 protease. The initial structures used for the molecular dynamics simulations were generated based on the Brownian dynamics trajectories, and this is the first molecular dynamics simulation of modeling the association of a ligand with the protease. We found that a protease substrate successfully binds to the protein when the flaps are fully open. Surprisingly, a smaller cyclic urea inhibitor (XK263) can reach the binding site when the flaps are not fully open. However, if the flaps are nearly closed, the inhibitor must rearrange or binding can fail because the inhibitor cannot attain proper conformations to enter the binding site. Both the peptide substrate and XK263 can also affect the protein's internal motion, which may help the flaps to open. Simulations allow us to efficiently study the ligand binding processes and may help those who study drug discovery to find optimal association pathways and to design those ligands with the best binding kinetics.
Collapse
Affiliation(s)
- Chia-En A Chang
- Department of Chemistry and Biochemistry, University of California at San Diego, La Jolla, California 92093, USA.
| | | | | | | |
Collapse
|
38
|
Abstract
Formation of a stereospecific protein complex is favored by specific interactions between two proteins but disfavored by the loss of translational and rotational freedom. Echoing the protein folding process, we have previously proposed a transition state for protein-protein association. Here we clarify the specification of the transition state by working with two types of toy models for protein association. A "hemisphere" model consists of two matching hemispheres as associating proteins, and a "crater" model consists of a spherical protein with a crater to which another spherical protein fits snugly. Short-range pairwise interactions between sites across the interface hold together the bound complex. Small relative translation and rotation between the subunits quickly destroy the interactions, leading to a sharp transition between the bound state with numerous short-range interactions but restricted translation and rotational freedom and the unbound state with, at most, a small number of interactions but expanded configurational freedom. This transition sets the outer boundary of the bound state as well as the transition state for association. The energy landscape is funnel-like, with the deep well of the bound state surrounded by a broad shallow basin. Calculations with the toy models suggest that mutational change in the interaction energy in the x-ray structure of a protein-protein complex, commonly used to approximate the effect on the association constant, overestimates the effect of mutation by 10-20%. With an eye toward specifying the transition states of actual protein complexes, we find that the total number of contacts between the subunits serves as a good surrogate of the interaction energy. This formalism of protein-protein association is applied to the barnase-barstar complex, reproducing the experimental results for the association rate over a wide range of ionic strength.
Collapse
Affiliation(s)
- Ramzi Alsallaq
- Department of Physics and Institute of Molecular Biophysics and School of Computational Science, Florida State University, Tallahassee, Florida, USA
| | | |
Collapse
|
39
|
Patriciu A, Chirikjian GS, Pappu RV. Analysis of the conformational dependence of mass-metric tensor determinants in serial polymers with constraints. J Chem Phys 2006; 121:12708-20. [PMID: 15606297 DOI: 10.1063/1.1821492] [Citation(s) in RCA: 22] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
It is well known that mass-metric tensor determinants det(G(s)) influence the equilibrium statistics and the rates of conformational transitions for polymers with constrained bond lengths and bond angles. It is now standard practice to include a Fixman-style compensating potential of the form U(c)(q(s)) proportional, variant(-k(B)T/2)ln[det(G(s))] as part of algorithms for torsional space molecular dynamics. This elegant strategy helps eliminate unwarranted biases that arise due to the imposition of holonomic constraints. However, the precise nature and extent of variation of det(G(s)) and hence ln[det(G(s))] with chain conformation and chain length has never been quantified. This type of analysis is crucial for understanding the nature of the conformational bias that the introduction of a Fixman potential aims to eliminate. Additionally, a detailed analysis of the conformational dependence of det(G(s)) will help resolve ambiguities regarding suggestions for incorporating terms related to det(G(s)) in the design of move sets in torsional space Monte Carlo simulations. In this work, we present results from a systematic study of the variation of det(G(s)) for a serial polymer with fixed bond lengths and bond angles as a function of chain conformation and chain length. This analysis requires an algorithm designed for rapid computation of det(G(s)) which simultaneously allows for a physical/geometric interpretation of the conformational dependence of det(G(s)). Consequently, we provide a detailed discussion of our adaptation of an O(n) algorithm from the robotics literature, which leads to simple recursion relations for direct evaluation of det(G(s)). Our analysis of the conformational dependence of det(G(s)) yields the following insights. (1) det(G(s)) is maximized for spatial conformers and minimized for planar conformations. (2) Previous work suggests that it is logical to expect that the conformational dependence of det(G(s)) becomes more pronounced with increase in chain length. Confirming this expectation, we provide systematic quantification of the nature of this dependency and show that the difference in det(G(s)) between spatial and planar conformers, i.e., between the maxima and minima of det(G(s)) grows systematically with chain length. Finally, we provide a brief discussion of implications of our analysis for the design of move sets in Monte Carlo simulations.
Collapse
Affiliation(s)
- Alexandru Patriciu
- Department of Mechanical Engineering, Johns Hopkins University, Baltimore, Maryland 21218, USA
| | | | | |
Collapse
|
40
|
Gitlin I, Carbeck JD, Whitesides GM. Warum sind Proteine geladen? Netzwerke aus Ladungs-Ladungs-Wechselwirkungen in Proteinen, analysiert über Ladungsleitern und Kapillarelektrophorese. Angew Chem Int Ed Engl 2006. [DOI: 10.1002/ange.200502530] [Citation(s) in RCA: 19] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022]
|
41
|
Lee MS, Olson MA. Calculation of absolute protein-ligand binding affinity using path and endpoint approaches. Biophys J 2005; 90:864-77. [PMID: 16284269 PMCID: PMC1367111 DOI: 10.1529/biophysj.105.071589] [Citation(s) in RCA: 142] [Impact Index Per Article: 7.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022] Open
Abstract
A comparative analysis is provided of rigorous and approximate methods for calculating absolute binding affinities of two protein-ligand complexes: the FKBP protein bound with small molecules 4-hydroxy-2-butanone and FK506. Our rigorous approach is an umbrella sampling technique where a potential of mean force is determined by pulling the ligand out of the protein active site over several simulation windows. The results of this approach agree well with experimentally observed binding affinities. Also assessed is a commonly used approximate endpoint approach, which separately estimates enthalpy, solvation free energy, and entropy. We show that this endpoint approach has numerous variations, all of which are prone to critical shortcomings. For example, conventional harmonic and quasiharmonic entropy estimation procedures produce disparate results for the relatively simple protein-ligand systems studied in this work.
Collapse
Affiliation(s)
- Michael S Lee
- Department of Cell Biology and Biochemistry, U.S. Army Medical Research Institute of Infectious Diseases, Frederick, Maryland 21702, USA.
| | | |
Collapse
|
42
|
Chang CE, Gilson MK. Free energy, entropy, and induced fit in host-guest recognition: calculations with the second-generation mining minima algorithm. J Am Chem Soc 2005; 126:13156-64. [PMID: 15469315 DOI: 10.1021/ja047115d] [Citation(s) in RCA: 176] [Impact Index Per Article: 9.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
This study applies a novel computational method to study molecular recognition for three sets of synthetic hosts: molecular clips, molecular tweezers, and a synthetic barbiturate receptor. The computed standard free energies of binding for the 12 binding reactions agree closely with experiment and provide insight into the roles of configurational entropy, preorganization, and induced fit in the systems studied. The computed changes in configurational entropy are comparable in magnitude to the changes in mean potential plus solvation energy, and they result primarily from changes in the average width of the energy wells upon binding. A strong correlation is observed between the changes in configurational energy and configurational entropy upon binding, resulting in near-linear compensation analogous to classical entropy-enthalpy compensation.
Collapse
Affiliation(s)
- Chia-En Chang
- Contribution from the Center for Advanced Research in Biotechnology, 9600 Gudelsky Drive, Rockville, MD 20850, USA
| | | |
Collapse
|