1
|
Leko K, Hanževački M, Brkljača Z, Pičuljan K, Ribić R, Požar J. Solvophobically Driven Complexation of Adamantyl Mannoside with β‐Cyclodextrin in Water and Structured Organic Solvents. Chemistry 2020; 26:5208-5219. [DOI: 10.1002/chem.202000282] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/17/2020] [Indexed: 01/20/2023]
Affiliation(s)
- Katarina Leko
- Department of ChemistryFaculty of ScienceUniversity of Zagreb Horvatovac 102a 10000 Zagreb Croatia
| | - Marko Hanževački
- Division of Physical ChemistryRuđer Bošković Institute Bijenička 54 10000 Zagreb Croatia
- Department of Chemical and Environmental EngineeringThe University of Nottingham University Park Nottingham NG7 2RD UK
| | - Zlatko Brkljača
- Division of Organic Chemistry and BiochemistryRuđer Bošković Institute Bijenička 54 10000 Zagreb Croatia
| | - Katarina Pičuljan
- Department of ChemistryFaculty of ScienceUniversity of Zagreb Horvatovac 102a 10000 Zagreb Croatia
| | - Rosana Ribić
- Department of ChemistryFaculty of ScienceUniversity of Zagreb Horvatovac 102a 10000 Zagreb Croatia
- University Center VaraždinUniversity North Jurja Križanića 31b 42000 Varaždin Croatia
| | - Josip Požar
- Department of ChemistryFaculty of ScienceUniversity of Zagreb Horvatovac 102a 10000 Zagreb Croatia
| |
Collapse
|
2
|
Hashikawa Y, Murata Y. A Single H2O Molecule inside Hydrophobic Carbon Nanocavities: Effect of Local Electrostatic Potential. CHEM LETT 2020. [DOI: 10.1246/cl.190874] [Citation(s) in RCA: 9] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/12/2022]
Affiliation(s)
- Yoshifumi Hashikawa
- Institute for Chemical Research, Kyoto University, Uji, Kyoto 611-0011, Japan
| | - Yasujiro Murata
- Institute for Chemical Research, Kyoto University, Uji, Kyoto 611-0011, Japan
| |
Collapse
|
3
|
Ahalawat N, Bandyopadhyay S, Mondal J. On the role of solvent in hydrophobic cavity–ligand recognition kinetics. J Chem Phys 2020; 152:074104. [DOI: 10.1063/1.5139584] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022] Open
Affiliation(s)
- Navjeet Ahalawat
- Tata Institute of Fundamental Research, Center for Interdisciplinary Sciences, Hyderabad 500107, India
- Department of Molecular Biology, Biotechnology and Bioinformatics, Chaudhary Charan Singh Haryana Agricultural University, Hisar 125004, India
| | - Satyabrata Bandyopadhyay
- Tata Institute of Fundamental Research, Center for Interdisciplinary Sciences, Hyderabad 500107, India
| | - Jagannath Mondal
- Tata Institute of Fundamental Research, Center for Interdisciplinary Sciences, Hyderabad 500107, India
| |
Collapse
|
4
|
Escalante DE, Aksan A. Role of Water Hydrogen Bonding on Transport of Small Molecules inside Hydrophobic Channels. J Phys Chem B 2019; 123:6673-6685. [PMID: 31310534 DOI: 10.1021/acs.jpcb.9b03060] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
We conducted a systematic analysis of water networking inside smooth hyperboloid hydrophobic structures (cylindrical, barrel, and hourglass shapes) to elucidate the role of water hydrogen bonding on the transport of small hydrophobic molecules (ligands). Through a series of molecular dynamics simulations, we established that a hydrogen-bonded network forming along the centerline resulted in a water exclusion zone adjacent to the walls. The size of the exclusion zone is a function of the geometry and the nonbonded interaction strength, defining the effective hydrophobicity of the structure. Exclusion of water molecules from this zone results in lower apparent viscosity, leading to acceleration of ligand transport up to 7 times faster than that measured in the bulk. Transport of ligands into and out of the hydrophobic structures was shown to be controlled by a single water molecule that capped the narrow regions in the structure. This mechanism provides physical insights into the behavior and role of water in the bottleneck regions of hydrophobic enzyme channels. These findings were then used in a sister publication [ Escalante , D. E. , Comput. Struct. Biotechnol. J. 2019 17 757 760 ] to develop a model that can accurately predict the transport of ligands along nanochannels of broad-substrate specificity enzymes.
Collapse
Affiliation(s)
- Diego E Escalante
- Department of Mechanical Engineering , University of Minnesota , Minneapolis , Minnesota 55455 , United States
| | - Alptekin Aksan
- Department of Mechanical Engineering , University of Minnesota , Minneapolis , Minnesota 55455 , United States.,BioTechnology Institute , University of Minnesota , St. Paul , Minnesota 55108 , United States
| |
Collapse
|
5
|
Variational implicit-solvent predictions of the dry-wet transition pathways for ligand-receptor binding and unbinding kinetics. Proc Natl Acad Sci U S A 2019; 116:14989-14994. [PMID: 31270236 DOI: 10.1073/pnas.1902719116] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022] Open
Abstract
Ligand-receptor binding and unbinding are fundamental biomolecular processes and particularly essential to drug efficacy. Environmental water fluctuations, however, impact the corresponding thermodynamics and kinetics and thereby challenge theoretical descriptions. Here, we devise a holistic, implicit-solvent, multimethod approach to predict the (un)binding kinetics for a generic ligand-pocket model. We use the variational implicit-solvent model (VISM) to calculate the solute-solvent interfacial structures and the corresponding free energies, and combine the VISM with the string method to obtain the minimum energy paths and transition states between the various metastable ("dry" and "wet") hydration states. The resulting dry-wet transition rates are then used in a spatially dependent multistate continuous-time Markov chain Brownian dynamics simulation and the related Fokker-Planck equation calculations of the ligand stochastic motion, providing the mean first-passage times for binding and unbinding. We find the hydration transitions to significantly slow down the binding process, in semiquantitative agreement with existing explicit-water simulations, but significantly accelerate the unbinding process. Moreover, our methods allow the characterization of nonequilibrium hydration states of pocket and ligand during the ligand movement, for which we find substantial memory and hysteresis effects for binding vs. unbinding. Our study thus provides a significant step forward toward efficient, physics-based interpretation and predictions of the complex kinetics in realistic ligand-receptor systems.
Collapse
|
6
|
Islam N, Flint M, Rick SW. Water hydrogen degrees of freedom and the hydrophobic effect. J Chem Phys 2019; 150:014502. [DOI: 10.1063/1.5053239] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Naeyma Islam
- Department of Chemistry, University of New Orleans, New Orleans, Louisiana 70148, USA
| | - Mahalia Flint
- Department of Chemistry, University of New Orleans, New Orleans, Louisiana 70148, USA
| | - Steven W. Rick
- Department of Chemistry, University of New Orleans, New Orleans, Louisiana 70148, USA
| |
Collapse
|
7
|
Weiß RG, Chudoba R, Setny P, Dzubiella J. Affinity, kinetics, and pathways of anisotropic ligands binding to hydrophobic model pockets. J Chem Phys 2018; 149:094902. [DOI: 10.1063/1.5025118] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/01/2023] Open
Affiliation(s)
- R. Gregor Weiß
- Institut für Physik, Humboldt-Universität zu Berlin, Newtonstraße 15, D-12489 Berlin, Germany
- Laboratory of Physical Chemistry, ETH Zürich, CH-8093 Zürich, Switzerland
| | - Richard Chudoba
- Institut für Physik, Humboldt-Universität zu Berlin, Newtonstraße 15, D-12489 Berlin, Germany
- Research Group Simulations of Energy Materials, Helmholtz-Zentrum Berlin, Hahn-Meitner-Platz 1, D-14109 Berlin, Germany
- Physikalisches Institut, Albert-Ludwigs-Universität Freiburg, Hermann-Herder Straße 3, D-79104 Freiburg, Germany
| | - Piotr Setny
- Centre of New Technologies, University of Warsaw, Stefana Banacha 2c, 00-927 Warsaw, Poland
| | - Joachim Dzubiella
- Institut für Physik, Humboldt-Universität zu Berlin, Newtonstraße 15, D-12489 Berlin, Germany
- Research Group Simulations of Energy Materials, Helmholtz-Zentrum Berlin, Hahn-Meitner-Platz 1, D-14109 Berlin, Germany
- Physikalisches Institut, Albert-Ludwigs-Universität Freiburg, Hermann-Herder Straße 3, D-79104 Freiburg, Germany
| |
Collapse
|
8
|
Ricci CG, Li B, Cheng LT, Dzubiella J, McCammon JA. Tailoring the Variational Implicit Solvent Method for New Challenges: Biomolecular Recognition and Assembly. Front Mol Biosci 2018; 5:13. [PMID: 29484300 PMCID: PMC5816062 DOI: 10.3389/fmolb.2018.00013] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/08/2017] [Accepted: 01/26/2018] [Indexed: 01/12/2023] Open
Abstract
Predicting solvation free energies and describing the complex water behavior that plays an important role in essentially all biological processes is a major challenge from the computational standpoint. While an atomistic, explicit description of the solvent can turn out to be too expensive in large biomolecular systems, most implicit solvent methods fail to capture “dewetting” effects and heterogeneous hydration by relying on a pre-established (i.e., guessed) solvation interface. Here we focus on the Variational Implicit Solvent Method, an implicit solvent method that adds water “plasticity” back to the picture by formulating the solvation free energy as a functional of all possible solvation interfaces. We survey VISM's applications to the problem of molecular recognition and report some of the most recent efforts to tailor VISM for more challenging scenarios, with the ultimate goal of including thermal fluctuations into the framework. The advances reported herein pave the way to make VISM a uniquely successful approach to characterize complex solvation properties in the recognition and binding of large-scale biomolecular complexes.
Collapse
Affiliation(s)
- Clarisse Gravina Ricci
- Department of Pharmacology and Department of Chemistry and Biochemistry, National Biomedical Computation Resource, University of California, San Diego, La Jolla, CA, United States
| | - Bo Li
- Department of Mathematics, University of California, San Diego, La Jolla, CA, United States.,Quantitative Biology Graduate Program, University of California, San Diego, La Jolla, CA, United States
| | - Li-Tien Cheng
- Department of Mathematics, University of California, San Diego, La Jolla, CA, United States
| | - Joachim Dzubiella
- Institut für Physik, Humboldt-Universität zu Berlin, Berlin, Germany.,Soft Matter and Functional Materials, Helmholtz-Zentrum Berlin, Berlin, Germany
| | - J Andrew McCammon
- Department of Pharmacology and Department of Chemistry and Biochemistry, National Biomedical Computation Resource, University of California, San Diego, La Jolla, CA, United States
| |
Collapse
|
9
|
Jalalitalab E, Abbaspour M, Akbarzadeh H. Thermodynamic, structural, and dynamical properties of nano-confined water using SPC/E and TIP4P models by molecular dynamics simulations. NEW J CHEM 2018. [DOI: 10.1039/c8nj01185k] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/31/2023]
Abstract
Different morphologies of water molecules are confined between two parallel graphene surfaces.
Collapse
|
10
|
Abstract
This review focuses on papers published since 2000 on the topic of the properties of solutes in water. More specifically, it evaluates the state of the art of our understanding of the complex relationship between the shape of a hydrophobe and the hydrophobic effect. To highlight this, we present a selection of references covering both empirical and molecular dynamics studies of small (molecular-scale) solutes. These include empirical studies of small molecules, synthetic hosts, crystalline monolayers, and proteins, as well as in silico investigations of entities such as idealized hard and soft spheres, small solutes, hydrophobic plates, artificial concavity, molecular hosts, carbon nanotubes and spheres, and proteins.
Collapse
Affiliation(s)
- Matthew B Hillyer
- Department of Chemistry, Tulane University, New Orleans, Louisiana 70118;
| | - Bruce C Gibb
- Department of Chemistry, Tulane University, New Orleans, Louisiana 70118;
| |
Collapse
|
11
|
Weiß RG, Setny P, Dzubiella J. Principles for Tuning Hydrophobic Ligand–Receptor Binding Kinetics. J Chem Theory Comput 2017; 13:3012-3019. [DOI: 10.1021/acs.jctc.7b00216] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/29/2023]
Affiliation(s)
- R. Gregor Weiß
- Institut
für Physik, Humboldt-Universität zu Berlin, Newtonstr. 15, D-12489 Berlin, Germany
- Institut
für Weiche Materie and Funktionale Materialen, Helmholtz-Zentrum Berlin, Hahn-Meitner Platz 1, D-14109 Berlin, Germany
| | - Piotr Setny
- Centre
of New Technologies, University of Warsaw, Banacha 2c, 02-097 Warsaw, Poland
| | - Joachim Dzubiella
- Institut
für Physik, Humboldt-Universität zu Berlin, Newtonstr. 15, D-12489 Berlin, Germany
- Institut
für Weiche Materie and Funktionale Materialen, Helmholtz-Zentrum Berlin, Hahn-Meitner Platz 1, D-14109 Berlin, Germany
| |
Collapse
|
12
|
Dominguez H, Pizio O, Pusztai L, Sokolowski S. The Structural Properties and Diffusion of a Three-Dimensional Isotropic Core-Softened Model Fluid in Disordered Porous Media. Molecular Dynamics Simulation. ADSORPT SCI TECHNOL 2016. [DOI: 10.1260/0263-6174.25.7.479] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/03/2022] Open
Abstract
The microscopic structure and dynamic properties of an isotropic three-dimensional core-softened model fluid in disordered matrices of Lennard-Jones particles have been studied. Molecular dynamics computer simulations in Grand Canonical ensemble were used as the methodological tools. It was shown that the microscopic structure of the fluid is characterized by anomalies similar to those found in a bulk model, but that it is affected by the fluid-matrix interactions. The dynamic properties also exhibit anomalous dependence on fluid density, but the magnitude of these anomalies is suppressed in comparison to the bulk fluid model. The anomalous behaviour of the diffusion coefficient is attributed to structural changes in the first coordination shell of a given fluid particle. It seems that the anomalies can only be suppressed at matrix densities which are higher than those studied in the present work.
Collapse
Affiliation(s)
- Hector Dominguez
- Instituto de Investigaciones en Materiales, UNAM, Coyoacan 04510, Mexico, D.F
| | - Orest Pizio
- Instituto de Quimica de la UNAM, Coyoacan 04510, Mexico, D.F
| | - Laszlo Pusztai
- Institute for Solid State Physics and Optics, Hungarian Academy of Sciences, Budapest, H-1121, Hungary
| | - Stefan Sokolowski
- Department for the Modelling of Physico-Chemical Processes, Maria Curie-Skłodowska University, Lublin 20031, Poland
| |
Collapse
|
13
|
Shevkunov SV. Computer simulation of the hydration of a chloride anion in a nanopore with hydrophilic walls. RUSSIAN JOURNAL OF PHYSICAL CHEMISTRY A 2016. [DOI: 10.1134/s0036024416050290] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/23/2022]
|
14
|
Sun H, Wen J, Zhao Y, Li B, McCammon JA. A self-consistent phase-field approach to implicit solvation of charged molecules with Poisson-Boltzmann electrostatics. J Chem Phys 2015; 143:243110. [PMID: 26723595 DOI: 10.1063/1.4932336] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/03/2023] Open
Abstract
Dielectric boundary based implicit-solvent models provide efficient descriptions of coarse-grained effects, particularly the electrostatic effect, of aqueous solvent. Recent years have seen the initial success of a new such model, variational implicit-solvent model (VISM) [Dzubiella, Swanson, and McCammon Phys. Rev. Lett. 96, 087802 (2006) and J. Chem. Phys. 124, 084905 (2006)], in capturing multiple dry and wet hydration states, describing the subtle electrostatic effect in hydrophobic interactions, and providing qualitatively good estimates of solvation free energies. Here, we develop a phase-field VISM to the solvation of charged molecules in aqueous solvent to include more flexibility. In this approach, a stable equilibrium molecular system is described by a phase field that takes one constant value in the solute region and a different constant value in the solvent region, and smoothly changes its value on a thin transition layer representing a smeared solute-solvent interface or dielectric boundary. Such a phase field minimizes an effective solvation free-energy functional that consists of the solute-solvent interfacial energy, solute-solvent van der Waals interaction energy, and electrostatic free energy described by the Poisson-Boltzmann theory. We apply our model and methods to the solvation of single ions, two parallel plates, and protein complexes BphC and p53/MDM2 to demonstrate the capability and efficiency of our approach at different levels. With a diffuse dielectric boundary, our new approach can describe the dielectric asymmetry in the solute-solvent interfacial region. Our theory is developed based on rigorous mathematical studies and is also connected to the Lum-Chandler-Weeks theory (1999). We discuss these connections and possible extensions of our theory and methods.
Collapse
Affiliation(s)
- Hui Sun
- Department of Mathematics, University of California, San Diego, 9500 Gilman Drive, Mail Code 0112. La Jolla, California 92093-0112, USA
| | - Jiayi Wen
- Department of Mathematics, University of California, San Diego, 9500 Gilman Drive, Mail Code 0112. La Jolla, California 92093-0112, USA
| | - Yanxiang Zhao
- Department of Mathematics, the George Washington University, Monroe Hall, 2115 G St. NW, Washington, DC 20052, USA
| | - Bo Li
- Department of Mathematics and Quantitative Biology Graduate Program, University of California, San Diego, 9500 Gilman Drive, Mail Code 0112. La Jolla, California 92093-0112, USA
| | - J Andrew McCammon
- Department of Chemistry and Biochemistry, Department of Pharmacology, and Howard Hughes Medical Institute, University of California, San Diego, La Jolla, California 92093-0365, USA
| |
Collapse
|
15
|
Abstract
In our previous work, we introduced a solvation model based on discrete solvent representation and demonstrated its ability to estimate hydration free energies for neutral solutes. Here, we present modifications extending the applicability of the model to charged solutes. They include improvements in the representation of the first hydration shell and systematic treatment of long-range interactions. While sharing computational efficiency of implicit solvent models, our approach avoids some of their important limitations, both in the context of electrostatic and nonpolar hydration effects: it naturally captures hydration asymmetry of opposite charges, it relies on solute description by standard all atom force fields instead of utilizing specialized sets of atomic parameters, it predicts solvent distribution in space without the need to geometrically define solvent accessible surface. By combining discrete solvent representation in vicinity of a solute with continuum description of long-range interactions, the model addresses two distinct aspects of biomolecular hydration: complex, short-range effects arising due to molecular nature of aqueous solvent, and bulk contributions. We demonstrate that the model provides good agreement with experimental results for an extensive set of roughly 700 diverse compounds, including neutral and charged solutes with hydration free energies ranging from +3.4 to -536 kcal/mol.
Collapse
Affiliation(s)
- Piotr Setny
- Centre of New Technologies, University of Warsaw, Banacha 2c, 02-097 Warsaw, Poland
| |
Collapse
|
16
|
Guo Z, Li B, Dzubiella J, Cheng LT, McCammon JA, Che J. Heterogeneous Hydration of p53/MDM2 Complex. J Chem Theory Comput 2014; 10:1302-1313. [PMID: 24803860 PMCID: PMC3958133 DOI: 10.1021/ct400967m] [Citation(s) in RCA: 18] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/05/2013] [Indexed: 12/23/2022]
Abstract
![]()
Water-mediated
interactions play critical roles in biomolecular
recognition processes. Explicit solvent molecular dynamics (MD) simulations
and the variational implicit-solvent model (VISM) are used to study
those hydration properties during binding for the biologically important
p53/MDM2 complex. Unlike simple model solutes, in such a realistic
and heterogeneous solute–solvent system with both geometrical
and chemical complexity, the local water distribution sensitively
depends on nearby amino acid properties and the geometric shape of
the protein. We show that the VISM can accurately describe the locations
of high and low density solvation shells identified by the MD simulations
and can explain them by a local coupling balance of solvent–solute
interaction potentials and curvature. In particular, capillary transitions
between local dry and wet hydration states in the binding pocket are
captured for interdomain distance between 4 to 6 Å, right at
the onset of binding. The underlying physical connection between geometry
and polarity is illustrated and quantified. Our study offers a microscopic
and physical insight into the heterogeneous hydration behavior of
the biologically highly relevant p53/MDM2 system and demonstrates
the fundamental importance of hydrophobic effects for biological binding
processes. We hope our study can help to establish new design rules
for drugs and medical substances.
Collapse
Affiliation(s)
- Zuojun Guo
- Genomics Institute of the Novartis Research Foundation , 10675 John Jay Hopkins Drive, San Diego, California 92121, United States
| | - Bo Li
- Department of Mathematics and Center for Theoretical Biological Physics, University of California, San Diego , 9500 Gilman Drive, La Jolla, California 92093-0112, United States
| | - Joachim Dzubiella
- Department of Physics, Humboldt University of Berlin , Newtonstr. 15, 12489 Berlin, Germany ; Soft Matter and Functional Materials, Helmholtz-Zentrum Berlin , Hahn-Meitner Platz 1, 14109 Berlin, Germany
| | - Li-Tien Cheng
- Department of Mathematics, University of California, San Diego , 9500 Gilman Drive, La Jolla, California 92093-0112, United States
| | - J Andrew McCammon
- Department of Chemistry and Biochemistry, Department of Pharmacology, Howard Hughes Medical Institute, and Center for Theoretical Biological Physics, University of California , San Diego, 9500 Gilman Drive, La Jolla, California 92093-0365, United States
| | - Jianwei Che
- Genomics Institute of the Novartis Research Foundation , 10675 John Jay Hopkins Drive, San Diego, California 92121, United States
| |
Collapse
|
17
|
Solvent fluctuations in hydrophobic cavity-ligand binding kinetics. Proc Natl Acad Sci U S A 2013; 110:1197-202. [PMID: 23297241 DOI: 10.1073/pnas.1221231110] [Citation(s) in RCA: 77] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022] Open
Abstract
Water plays a crucial part in virtually all protein-ligand binding processes in and out of equilibrium. Here, we investigate the role of water in the binding kinetics of a ligand to a prototypical hydrophobic pocket by explicit-water molecular dynamics (MD) simulations and implicit diffusional approaches. The concave pocket in the unbound state exhibits wet/dry hydration oscillations whose magnitude and time scale are significantly amplified by the approaching ligand. In turn, the ligand's stochastic motion intimately couples to the slow hydration fluctuations, leading to a sixfold-enhanced friction in the vicinity of the pocket entrance. The increased friction considerably decelerates association in the otherwise barrierless system, indicating the importance of molecular-scale hydrodynamic effects in cavity-ligand binding arising due to capillary fluctuations. We derive and analyze the diffusivity profile and show that the mean first passage time distribution from the MD simulation can be accurately reproduced by a standard Brownian dynamics simulation if the appropriate position-dependent friction profile is included. However, long-time decays in the water-ligand (random) force autocorrelation demonstrate violation of the Markovian assumption, challenging standard diffusive approaches for rate prediction. Remarkably, the static friction profile derived from the force correlations strongly resembles the profile derived on the Markovian assumption apart from a simple shift in space, which can be rationalized by a time-space retardation in the ligand's downhill dynamics toward the pocket. The observed spatiotemporal hydrodynamic coupling may be of biological importance providing the time needed for conformational receptor-ligand adjustments, typical of the induced-fit paradigm.
Collapse
|
18
|
Setny P, Baron R, McCammon JA. Comment on ‘Molecular driving forces of the pocket-ligand hydrophobic association’ by G. Graziano, Chem. Phys. Lett. 533 (2012) 95. Chem Phys Lett 2013. [DOI: 10.1016/j.cplett.2012.11.004] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/25/2022]
|
19
|
Baron R, Setny P, Paesani F. Water Structure, Dynamics, and Spectral Signatures: Changes Upon Model Cavity–Ligand Recognition. J Phys Chem B 2012; 116:13774-80. [DOI: 10.1021/jp309373q] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Riccardo Baron
- Department of Chemistry and
Biochemistry, University of California, San Diego, La Jolla, California 92093, United States
- Center for Theoretical Biological
Physics, University of California, San Diego, La Jolla, California 92093, United States
| | - Piotr Setny
- Physics Department T38, Technical University Munich, James Franck Str. 1, 85748
Garching, Germany
| | - Francesco Paesani
- Department of Chemistry and
Biochemistry, University of California, San Diego, La Jolla, California 92093, United States
| |
Collapse
|
20
|
Baron R, Molinero V. Water-Driven Cavity–Ligand Binding: Comparison of Thermodynamic Signatures from Coarse-Grained and Atomic-Level Simulations. J Chem Theory Comput 2012; 8:3696-704. [DOI: 10.1021/ct300121r] [Citation(s) in RCA: 42] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Riccardo Baron
- Department of Medicinal Chemistry,
College of Pharmacy, and The Henry Eyring Center for Theoretical Chemistry,
The University of Utah, Salt Lake City, Utah 84112-5820, United States
| | - Valeria Molinero
- Department of Chemistry and
The Henry Eyring Center for Theoretical Chemistry, The University
of Utah, 315 South 1400 East, Salt Lake City, Utah 84112-0850, United
States
| |
Collapse
|
21
|
Gibb CLD, Gibb BC. The Thermodynamics of Molecular Recognition. Supramol Chem 2012. [DOI: 10.1002/9780470661345.smc005] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Indexed: 11/10/2022]
|
22
|
Gibb BC. Supramolecular Assembly and Binding in Aqueous Solution: Useful Tips Regarding the Hofmeister and Hydrophobic Effects. Isr J Chem 2011. [DOI: 10.1002/ijch.201100058] [Citation(s) in RCA: 29] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/11/2022]
|
23
|
Lawrenz M, Baron R, Wang Y, McCammon JA. Effects of Biomolecular Flexibility on Alchemical Calculations of Absolute Binding Free Energies. J Chem Theory Comput 2011; 7:2224-2232. [PMID: 21811708 PMCID: PMC3146769 DOI: 10.1021/ct200230v] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
The independent trajectory thermodynamic integration (IT-TI) approach (Lawrenz et. al J. Chem. Theory. Comput. 2009, 5:1106-1116(1)) for free energy calculations with distributed computing is employed to study two distinct cases of protein-ligand binding: first, the influenza surface protein N1 neuraminidase bound to the inhibitor oseltamivir, and second, the M. tuberculosis enzyme RmlC complexed with the molecule CID 77074. For both systems, finite molecular dynamics (MD) sampling and varied molecular flexibility give rise to IT-TI free energy distributions that are remarkably centered on the target experimental values, with a spread directly related to protein, ligand, and solvent dynamics. Using over 2 μs of total MD simulation, alternative protocols for the practical, general implementation of IT-TI are investigated, including the optimal use of distributed computing, the total number of alchemical intermediates, and the procedure to perturb electrostatics and van der Waals interactions. A protocol that maximizes predictive power and computational efficiency is proposed. IT-TI outperforms traditional TI predictions and allows a straightforward evaluation of the reliability of free energy estimates. Our study has broad implications for the use of distributed computing in free energy calculations of macromolecular systems.
Collapse
Affiliation(s)
- Morgan Lawrenz
- Department of Chemistry and Biochemistry, Center for Theoretical Biological Physics, University of California, San Diego, La Jolla, California, U.S.A
| | - Riccardo Baron
- Department of Chemistry and Biochemistry, Center for Theoretical Biological Physics, University of California, San Diego, La Jolla, California, U.S.A
| | - Yi Wang
- Department of Chemistry and Biochemistry, Center for Theoretical Biological Physics, University of California, San Diego, La Jolla, California, U.S.A
| | - J. Andrew McCammon
- Department of Chemistry and Biochemistry, Center for Theoretical Biological Physics, University of California, San Diego, La Jolla, California, U.S.A
- Department of Pharmacology, Howard Hughes Medical Institute, University of California, San Diego, La Jolla, California, U.S.A
| |
Collapse
|
24
|
Laughrey Z, Gibb BC. Water-soluble, self-assembling container molecules: an update. Chem Soc Rev 2011; 40:363-86. [DOI: 10.1039/c0cs00030b] [Citation(s) in RCA: 145] [Impact Index Per Article: 11.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/14/2023]
|
25
|
Setny P, Baron R, McCammon JA. How Can Hydrophobic Association Be Enthalpy Driven? J Chem Theory Comput 2010; 6:2866-2871. [PMID: 20844599 PMCID: PMC2938831 DOI: 10.1021/ct1003077] [Citation(s) in RCA: 176] [Impact Index Per Article: 12.6] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/05/2010] [Indexed: 01/19/2023]
Abstract
Hydrophobic association is often recognized as being driven by favorable entropic contributions. Here, using explicit solvent molecular dynamics simulations we investigate binding in a model hydrophobic receptor-ligand system which appears, instead, to be driven by enthalpy and opposed by entropy. We use the temperature dependence of the potential of mean force to analyze the thermodynamic contributions along the association coordinate. Relating such contributions to the ongoing changes in system hydration allows us to demonstrate that the overall binding thermodynamics is determined by the expulsion of disorganized water from the receptor cavity. Our model study sheds light on the solvent-induced driving forces for receptor-ligand association of general, transferable relevance for biological systems with poorly hydrated binding sites.
Collapse
Affiliation(s)
- Piotr Setny
- Corresponding authors. E-mail: (P.S.), (R.B.)
| | | | | |
Collapse
|
26
|
Cheng LT, Wang Z, Setny P, Dzubiella J, Li B, McCammon JA. Interfaces and hydrophobic interactions in receptor-ligand systems: A level-set variational implicit solvent approach. J Chem Phys 2010; 131:144102. [PMID: 19831428 DOI: 10.1063/1.3242274] [Citation(s) in RCA: 36] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/19/2023] Open
Abstract
A model nanometer-sized hydrophobic receptor-ligand system in aqueous solution is studied by the recently developed level-set variational implicit solvent model (VISM). This approach is compared to all-atom computer simulations. The simulations reveal complex hydration effects within the (concave) receptor pocket, sensitive to the distance of the (convex) approaching ligand. The ligand induces and controls an intermittent switching between dry and wet states of the hosting pocket, which determines the range and magnitude of the pocket-ligand attraction. In the level-set VISM, a geometric free-energy functional of all possible solute-solvent interfaces coupled to the local dispersion potential is minimized numerically. This approach captures the distinct metastable states that correspond to topologically different solute-solvent interfaces, and thereby reproduces the bimodal hydration behavior observed in the all-atom simulation. Geometrical singularities formed during the interface relaxation are found to contribute significantly to the energy barrier between different metastable states. While the hydration phenomena can thus be explained by capillary effects, the explicit inclusion of dispersion and curvature corrections seems to be essential for a quantitative description of hydrophobically confined systems on nanoscales. This study may shed more light onto the tight connection between geometric and energetic aspects of biomolecular hydration and may represent a valuable step toward the proper interpretation of experimental receptor-ligand binding rates.
Collapse
Affiliation(s)
- Li-Tien Cheng
- Department of Mathematics, University of California San Diego, La Jolla, California 92093, USA.
| | | | | | | | | | | |
Collapse
|
27
|
Forces mediating protein–protein interactions: a computational study of p53 “approaching” MDM2. Theor Chem Acc 2009. [DOI: 10.1007/s00214-009-0682-1] [Citation(s) in RCA: 10] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/20/2022]
|
28
|
Abstract
Hydrophobicity manifests itself differently on large and small length scales. This review focuses on large-length-scale hydrophobicity, particularly on dewetting at single hydrophobic surfaces and drying in regions bounded on two or more sides by hydrophobic surfaces. We review applicable theories, simulations, and experiments pertaining to large-scale hydrophobicity in physical and biomolecular systems and clarify some of the critical issues pertaining to this subject. Given space constraints, we cannot review all the significant and interesting work in this active field.
Collapse
Affiliation(s)
- Bruce J Berne
- Department of Chemistry, Columbia University, New York, New York 10027, USA.
| | | | | |
Collapse
|
29
|
Abstract
The structure and dynamics of water inside a water-soluble, bowl-shaped cavitand molecule with a hydrophobic interior are studied using molecular dynamics computer simulations. The simulations find that the number of inside water molecules is about 4.5, but it fluctuates from being completely empty to full on a time scale of tens of nanoseconds. The transition from empty to full is energetically favorable and entropically unfavorable. The water molecules inside have fewer hydrogen bonds than the bulk and in general weaker interactions; the lower energy results from the nearest-neighbor interactions with the cavitand atoms and the water molecules at the entrance of the cavitand, interactions that are lost upon dewetting. An analysis of translational and rotational motion suggests that the lower entropy of the inside water molecules is due to decreased translational entropy, which outweighs an increased orientational entropy. The cavitand molecule acts as a host binding hydrophobic guests, and dewetting can be induced by the presence of a hydrophobic guest molecule about 3 A above the entrance. At this position, the guest displaces the water molecules which stabilize the inside water molecules and the empty cavitand becomes more stable than the full.
Collapse
Affiliation(s)
- Jeffrey Ewell
- Department of Chemistry, University of New Orleans, New Orleans, LA 70148, USA
| | | | | |
Collapse
|
30
|
Setny P. Hydrophobic interactions between methane and a nanoscopic pocket: Three dimensional distribution of potential of mean force revealed by computer simulations. J Chem Phys 2008; 128:125105. [DOI: 10.1063/1.2839885] [Citation(s) in RCA: 20] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
31
|
Abstract
We suggest that bubbles are the bistable hydrophobic gates responsible for the on-off transitions of single channel currents. In this view, many types of channels gate by the same physical mechanism-dewetting by capillary evaporation-but different types of channels use different sensors to modulate hydrophobic properties of the channel wall and thereby trigger and control bubbles and gating. Spontaneous emptying of channels has been seen in many simulations. Because of the physics involved, such phase transitions are inherently sensitive, unstable threshold phenomena that are difficult to simulate reproducibly and thus convincingly. We present a thermodynamic analysis of a bubble gate using morphometric density functional theory of classical (not quantum) mechanics. Thermodynamic analysis of phase transitions is generally more reproducible and less sensitive to details than simulations. Anesthetic actions of inert gases-and their interactions with hydrostatic pressure (e.g., nitrogen narcosis)-can be easily understood by actions on bubbles. A general theory of gas anesthesia may involve bubbles in channels. Only experiments can show whether, or when, or which channels actually use bubbles as hydrophobic gates: direct observation of bubbles in channels is needed. Existing experiments show thin gas layers on hydrophobic surfaces in water and suggest that bubbles nearly exist in bulk water.
Collapse
|
32
|
Setny P. Water properties and potential of mean force for hydrophobic interactions of methane and nanoscopic pockets studied by computer simulations. J Chem Phys 2007; 127:054505. [PMID: 17688347 DOI: 10.1063/1.2749250] [Citation(s) in RCA: 23] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We consider model systems consisting of a methane molecule and hemispherical pockets of subnanometer radii whose walls are made of hydrophobic material. The potential of mean force for process of translocation of the methane molecule from bulk water into the pockets' interior is obtained, based on an explicit solvent molecular dynamics simulations. Accompanying changes in water density around the interacting objects and spatial distribution of solvent's potential energy are analyzed, allowing for interpretation of details of hydrophobic interactions in relation to hydrophobic hydration properties. Applicability of surface area-based models of hydrophobic effect for systems of interest is also investigated. A total work for the translocation process is not dependent on pocket's size, indicating that pocket desolvation has little contribution to free energy changes, which is consistent with the observation that solvent density is significantly reduced inside "unperturbed" pockets. Substantial solvent effects are shown to have a longer range than in case of a well investigated methane pair. A desolvation barrier is present in a smaller pocket system but disappears in the larger one, suggesting that a form of a "hydrophobic collapse" is observed.
Collapse
Affiliation(s)
- Piotr Setny
- Department of Biophysics, Institute of Experimental Physics, University of Warsaw, 02-089 Warsaw, Poland
| |
Collapse
|