1
|
Kaiser S, Yue Z, Peng Y, Nguyen TD, Chen S, Teng D, Voth GA. Molecular Dynamics Simulation of Complex Reactivity with the Rapid Approach for Proton Transport and Other Reactions (RAPTOR) Software Package. J Phys Chem B 2024; 128:4959-4974. [PMID: 38742764 PMCID: PMC11129700 DOI: 10.1021/acs.jpcb.4c01987] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/26/2024] [Revised: 05/05/2024] [Accepted: 05/06/2024] [Indexed: 05/16/2024]
Abstract
Simulating chemically reactive phenomena such as proton transport on nanosecond to microsecond and beyond time scales is a challenging task. Ab initio methods are unable to currently access these time scales routinely, and traditional molecular dynamics methods feature fixed bonding arrangements that cannot account for changes in the system's bonding topology. The Multiscale Reactive Molecular Dynamics (MS-RMD) method, as implemented in the Rapid Approach for Proton Transport and Other Reactions (RAPTOR) software package for the LAMMPS molecular dynamics code, offers a method to routinely sample longer time scale reactive simulation data with statistical precision. RAPTOR may also be interfaced with enhanced sampling methods to drive simulations toward the analysis of reactive rare events, and a number of collective variables (CVs) have been developed to facilitate this. Key advances to this methodology, including GPU acceleration efforts and novel CVs to model water wire formation are reviewed, along with recent applications of the method which demonstrate its versatility and robustness.
Collapse
Affiliation(s)
- Scott Kaiser
- Department
of Chemistry, Chicago Center for Theoretical Chemistry, James Franck
Institute, and Institute for Biophysical Dynamics, The University of Chicago, Chicago, Illinois 60637, United States
| | - Zhi Yue
- Department
of Chemistry, Chicago Center for Theoretical Chemistry, James Franck
Institute, and Institute for Biophysical Dynamics, The University of Chicago, Chicago, Illinois 60637, United States
| | - Yuxing Peng
- NVIDIA
Corporation, Santa
Clara, California 95051, United States
| | - Trung Dac Nguyen
- Research
Computing Center, The University of Chicago, Chicago, Illinois 60637, United States
| | - Sijia Chen
- Department
of Chemistry, Chicago Center for Theoretical Chemistry, James Franck
Institute, and Institute for Biophysical Dynamics, The University of Chicago, Chicago, Illinois 60637, United States
| | - Da Teng
- Department
of Chemistry, Chicago Center for Theoretical Chemistry, James Franck
Institute, and Institute for Biophysical Dynamics, The University of Chicago, Chicago, Illinois 60637, United States
| | - Gregory A. Voth
- Department
of Chemistry, Chicago Center for Theoretical Chemistry, James Franck
Institute, and Institute for Biophysical Dynamics, The University of Chicago, Chicago, Illinois 60637, United States
| |
Collapse
|
2
|
Kratochvil HT, Watkins LC, Mravic M, Thomaston JL, Nicoludis JM, Somberg NH, Liu L, Hong M, Voth GA, DeGrado WF. Transient water wires mediate selective proton transport in designed channel proteins. Nat Chem 2023; 15:1012-1021. [PMID: 37308712 PMCID: PMC10475958 DOI: 10.1038/s41557-023-01210-4] [Citation(s) in RCA: 11] [Impact Index Per Article: 11.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/19/2022] [Accepted: 04/19/2023] [Indexed: 06/14/2023]
Abstract
Selective proton transport through proteins is essential for forming and using proton gradients in cells. Protons are conducted along hydrogen-bonded 'wires' of water molecules and polar side chains, which, somewhat surprisingly, are often interrupted by dry apolar stretches in the conduction pathways, inferred from static protein structures. Here we hypothesize that protons are conducted through such dry spots by forming transient water wires, often highly correlated with the presence of the excess protons in the water wire. To test this hypothesis, we performed molecular dynamics simulations to design transmembrane channels with stable water pockets interspersed by apolar segments capable of forming flickering water wires. The minimalist designed channels conduct protons at rates similar to viral proton channels, and they are at least 106-fold more selective for H+ over Na+. These studies inform the mechanisms of biological proton conduction and the principles for engineering proton-conductive materials.
Collapse
Affiliation(s)
- Huong T Kratochvil
- Department of Pharmaceutical Chemistry, University of California-San Francisco, San Francisco, CA, USA.
- Department of Chemistry, University of North Carolina at Chapel Hill, Chapel Hill, NC, USA.
| | - Laura C Watkins
- Department of Chemistry, Chicago Center for Theoretical Chemistry, Institute for Biophysical Dynamics and James Franck Institute, The University of Chicago, Chicago, IL, USA
- Kemper Insurance, Chicago, IL, USA
| | - Marco Mravic
- Department of Pharmaceutical Chemistry, University of California-San Francisco, San Francisco, CA, USA
- Department of Integrative Structural and Computational Biology Scripps Research Institute, La Jolla, CA, USA
| | - Jessica L Thomaston
- Department of Pharmaceutical Chemistry, University of California-San Francisco, San Francisco, CA, USA
| | - John M Nicoludis
- Department of Pharmaceutical Chemistry, University of California-San Francisco, San Francisco, CA, USA
- Genentech, San Francisco, CA, USA
| | - Noah H Somberg
- Department of Chemistry, Massachusetts Institute of Technology, Cambridge, MA, USA
| | | | - Mei Hong
- Department of Chemistry, Massachusetts Institute of Technology, Cambridge, MA, USA
| | - Gregory A Voth
- Department of Chemistry, Chicago Center for Theoretical Chemistry, Institute for Biophysical Dynamics and James Franck Institute, The University of Chicago, Chicago, IL, USA.
| | - William F DeGrado
- Department of Pharmaceutical Chemistry, University of California-San Francisco, San Francisco, CA, USA.
| |
Collapse
|
3
|
Chen S, Li Z, Voth GA. Acidic Conditions Impact Hydrophobe Transfer across the Oil-Water Interface in Unusual Ways. J Phys Chem B 2023; 127:3911-3918. [PMID: 37084419 PMCID: PMC10166083 DOI: 10.1021/acs.jpcb.3c00828] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/06/2023] [Revised: 04/09/2023] [Indexed: 04/23/2023]
Abstract
Molecular dynamics simulation and enhanced free energy sampling are used to study hydrophobic solute transfer across the water-oil interface with explicit consideration of the effect of different electrolytes: hydronium cation (hydrated excess proton) and sodium cation, both with chloride counterions (i.e., dissociated acid and salt, HCl and NaCl). With the Multistate Empirical Valence Bond (MS-EVB) methodology, we find that, surprisingly, hydronium can to a certain degree stabilize the hydrophobic solute, neopentane, in the aqueous phase and including at the oil-water interface. At the same time, the sodium cation tends to "salt out" the hydrophobic solute in the expected fashion. When it comes to the solvation structure of the hydrophobic solute in the acidic conditions, hydronium shows an affinity to the hydrophobic solute, as suggested by the radial distribution functions (RDFs). Upon consideration of this interfacial effect, we find that the solvation structure of the hydrophobic solute varies at different distances from the oil-liquid interface due to a competition between the bulk oil phase and the hydrophobic solute phase. Together with an observed orientational preference of the hydroniums and the lifetime of water molecules in the first solvation shell of neopentane, we conclude that hydronium stabilizes to a certain degree the dispersal of neopentane in the aqueous phase and eliminates any salting out effect in the acid solution; i.e., the hydronium acts like a surfactant. The present molecular dynamics study provides new insight into the hydrophobic solute transfer across the water-oil interface process, including for acid and salt solutions.
Collapse
Affiliation(s)
- Sijia Chen
- Department of Chemistry, Chicago Center
for Theoretical Chemistry, The James Franck Institute, and Institute
for Biophysical Dynamics, The University
of Chicago, Chicago, Illinois 60637, United States
| | - Zhefu Li
- Department of Chemistry, Chicago Center
for Theoretical Chemistry, The James Franck Institute, and Institute
for Biophysical Dynamics, The University
of Chicago, Chicago, Illinois 60637, United States
| | - Gregory A. Voth
- Department of Chemistry, Chicago Center
for Theoretical Chemistry, The James Franck Institute, and Institute
for Biophysical Dynamics, The University
of Chicago, Chicago, Illinois 60637, United States
| |
Collapse
|
4
|
Liu Y, Li C, Voth GA. Generalized Transition State Theory Treatment of Water-Assisted Proton Transport Processes in Proteins. J Phys Chem B 2022; 126:10452-10459. [PMID: 36459423 PMCID: PMC9762399 DOI: 10.1021/acs.jpcb.2c06703] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/20/2022] [Revised: 11/15/2022] [Indexed: 12/03/2022]
Abstract
Transition state theory (TST) is widely employed for estimating the transition rate of a reaction when combined with free energy sampling techniques. A derivation of the transition theory rate expression for a general n-dimensional case is presented in this work which specifically focuses on water-assisted proton transfer/transport reactions, especially for protein systems. Our work evaluates the TST prefactor calculated at the transition state dividing surface compared to one sampled, as an approximation, in the reactant state in four case studies of water-assisted proton transport inside membrane proteins and highlights the significant impact of the prefactor position dependence in proton transport processes.
Collapse
Affiliation(s)
- Yu Liu
- Department of Chemistry, Chicago Center
for Theoretical Chemistry, James Franck Institute, and Institute for
Biophysical Dynamics, The University of
Chicago, Chicago, Illinois60637, United States
| | - Chenghan Li
- Department of Chemistry, Chicago Center
for Theoretical Chemistry, James Franck Institute, and Institute for
Biophysical Dynamics, The University of
Chicago, Chicago, Illinois60637, United States
| | - Gregory A. Voth
- Department of Chemistry, Chicago Center
for Theoretical Chemistry, James Franck Institute, and Institute for
Biophysical Dynamics, The University of
Chicago, Chicago, Illinois60637, United States
| |
Collapse
|
5
|
Zuchniarz J, Liu Y, Li C, Voth GA. Accurate p Ka Calculations in Proteins with Reactive Molecular Dynamics Provide Physical Insight Into the Electrostatic Origins of Their Values. J Phys Chem B 2022; 126:7321-7330. [PMID: 36106487 PMCID: PMC9528908 DOI: 10.1021/acs.jpcb.2c04899] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/11/2022] [Revised: 08/28/2022] [Indexed: 11/29/2022]
Abstract
Classical molecular dynamics simulations are a versatile tool in the study of biomolecular systems, but they usually rely on a fixed bonding topology, precluding the explicit simulation of chemical reactivity. Certain modifications can permit the modeling of reactions. One such method, multiscale reactive molecular dynamics, makes use of a linear combination approach to describe condensed-phase free energy surfaces of reactive processes of biological interest. Before these simulations can be performed, models of the reactive moieties must first be parametrized using electronic structure data. A recent study demonstrated that gas-phase electronic structure data can be used to derive parameters for glutamate and lysine which reproduce experimental pKa values in both bulk water and the staphylococcal nuclease protein with remarkable accuracy and transferability between the water and protein environments. In this work, we first present a new model for aspartate derived in similar fashion and demonstrate that it too produces accurate pKa values in both bulk and protein contexts. We also describe a modification to the prior methodology, involving refitting some of the classical force field parameters to density functional theory calculations, which improves the transferability of the existing glutamate model. Finally and most importantly, this reactive molecular dynamics approach, based on rigorous statistical mechanics, allows one to specifically analyze the fundamental physical causes for the marked pKa shift of both aspartate and glutamate between bulk water and protein and also to demonstrate that local steric and electrostatic effects largely explain the observed differences.
Collapse
Affiliation(s)
- Joshua Zuchniarz
- Department of Chemistry, Chicago Center
for Theoretical Chemistry, James Franck Institute, and Institute for
Biophysical Dynamics, The University of
Chicago, Chicago, Illinois 60637, United States
| | - Yu Liu
- Department of Chemistry, Chicago Center
for Theoretical Chemistry, James Franck Institute, and Institute for
Biophysical Dynamics, The University of
Chicago, Chicago, Illinois 60637, United States
| | - Chenghan Li
- Department of Chemistry, Chicago Center
for Theoretical Chemistry, James Franck Institute, and Institute for
Biophysical Dynamics, The University of
Chicago, Chicago, Illinois 60637, United States
| | - Gregory A. Voth
- Department of Chemistry, Chicago Center
for Theoretical Chemistry, James Franck Institute, and Institute for
Biophysical Dynamics, The University of
Chicago, Chicago, Illinois 60637, United States
| |
Collapse
|
6
|
Proton Coupling and the Multiscale Kinetic Mechanism of a Peptide Transporter. Biophys J 2022; 121:2266-2278. [PMID: 35614850 DOI: 10.1016/j.bpj.2022.05.029] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/24/2022] [Revised: 05/13/2022] [Accepted: 05/19/2022] [Indexed: 11/02/2022] Open
Abstract
Proton coupled peptide transporters (POTs) are crucial for the uptake of di- and tri-peptides as well as drug and pro-drug molecules in prokaryotes and eukaryotic cells. We illustrate from multiscale modeling how transmembrane proton flux couples within a POT protein to drive essential steps of the full functional cycle: 1) protonation of a glutamate on transmembrane helix (TM) 7 opens the extracellular gate, allowing ligand entry; 2) inward proton flow induces the cytosolic release of ligand by varying the protonation state of a second conserved glutamate on TM10; 3) proton movement between TM7 and TM10 is thermodynamically driven and kinetically permissible via water proton shuttling without the participation of ligand. Our results, for the first time, give direct computational confirmation for the alternating access model of POTs, and point to a quantitative multiscale kinetic picture of the functioning protein mechanism.
Collapse
|
7
|
Watkins LC, DeGrado WF, Voth GA. Multiscale Simulation of an Influenza A M2 Channel Mutant Reveals Key Features of Its Markedly Different Proton Transport Behavior. J Am Chem Soc 2022; 144:769-776. [PMID: 34985907 PMCID: PMC8834648 DOI: 10.1021/jacs.1c09281] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/21/2023]
Abstract
The influenza A M2 channel, a prototype for viroporins, is an acid-activated viroporin that conducts protons across the viral membrane, a critical step in the viral life cycle. Four central His37 residues control channel activation by binding subsequent protons from the viral exterior, which opens the Trp41 gate and allows proton flux to the interior. Asp44 is essential for maintaining the Trp41 gate in a closed state at high pH, resulting in asymmetric conduction. The prevalent D44N mutant disrupts this gate and opens the C-terminal end of the channel, resulting in increased conduction and a loss of this asymmetric conduction. Here, we use extensive Multiscale Reactive Molecular Dynamics (MS-RMD) and quantum mechanics/molecular mechanics (QM/MM) molecular dynamics simulations with an explicit, reactive excess proton to calculate the free energy of proton transport in this M2 mutant and to study the dynamic molecular-level behavior of D44N M2. We find that this mutation significantly lowers the barrier of His37 deprotonation in the activated state and shifts the barrier for entry to the Val27 tetrad. These free energy changes are reflected in structural shifts. Additionally, we show that the increased hydration around the His37 tetrad diminishes the effect of the His37 charge on the channel's water structure, facilitating proton transport and enabling activation from the viral interior. Altogether, this work provides key insight into the fundamental characteristics of PT in WT M2 and how the D44N mutation alters this PT mechanism, and it expands understanding of the role of emergent mutations in viroporins.
Collapse
Affiliation(s)
- Laura C. Watkins
- Department of Chemistry, Chicago Center for Theoretical Chemistry, Institute for Biophysical Dynamics and James Franck Institute, The University of Chicago, Chicago, Illinois 60637, United States
| | - William F. DeGrado
- Department of Pharmaceutical Chemistry, University of California, San Francisco, San Francisco, California, 94158, United States
| | - Gregory A. Voth
- Department of Chemistry, Chicago Center for Theoretical Chemistry, Institute for Biophysical Dynamics and James Franck Institute, The University of Chicago, Chicago, Illinois 60637, United States,Corresponding Author
| |
Collapse
|
8
|
Hu Y, Wang S, He Y, An L. Evaluation of proton transport and solvation effect in hydrated Nafion membrane with degradation. Phys Chem Chem Phys 2022; 24:29024-29033. [DOI: 10.1039/d2cp02817d] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/26/2022]
Abstract
In proton exchange membrane fuel cells (PEMFCs), free radicals easily attack ionomers, resulting in membrane degradation.
Collapse
Affiliation(s)
- Yu Hu
- School of Energy Science and Engineering, Harbin Institute of Technology, Harbin, 150001, China
| | - Shuai Wang
- School of Energy Science and Engineering, Harbin Institute of Technology, Harbin, 150001, China
| | - Yurong He
- School of Energy Science and Engineering, Harbin Institute of Technology, Harbin, 150001, China
| | - Liang An
- Department of Mechanical Engineering, The Hong Kong Polytechnic University, Hung Hom, Kowloon, Hong Kong SAR, China
| |
Collapse
|
9
|
A quantitative paradigm for water-assisted proton transport through proteins and other confined spaces. Proc Natl Acad Sci U S A 2021; 118:2113141118. [PMID: 34857630 DOI: 10.1073/pnas.2113141118] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 10/15/2021] [Indexed: 11/18/2022] Open
Abstract
Water-assisted proton transport through confined spaces influences many phenomena in biomolecular and nanomaterial systems. In such cases, the water molecules that fluctuate in the confined pathways provide the environment and the medium for the hydrated excess proton migration via Grotthuss shuttling. However, a definitive collective variable (CV) that accurately couples the hydration and the connectivity of the proton wire with the proton translocation has remained elusive. To address this important challenge-and thus to define a quantitative paradigm for facile proton transport in confined spaces-a CV is derived in this work from graph theory, which is verified to accurately describe water wire formation and breakage coupled to the proton translocation in carbon nanotubes and the Cl-/H+ antiporter protein, ClC-ec1. Significant alterations in the conformations and thermodynamics of water wires are uncovered after introducing an excess proton into them. Large barriers in the proton translocation free-energy profiles are found when water wires are defined to be disconnected according to the new CV, even though the pertinent confined space is still reasonably well hydrated and-by the simple measure of the mere existence of a water structure-the proton transport would have been predicted to be facile via that oversimplified measure. In this paradigm, however, the simple presence of water is not sufficient for inferring proton translocation, since an excess proton itself is able to drive hydration, and additionally, the water molecules themselves must be adequately connected to facilitate any successful proton transport.
Collapse
|
10
|
Edwards T, Foloppe N, Harris SA, Wells G. The future of biomolecular simulation in the pharmaceutical industry: what we can learn from aerodynamics modelling and weather prediction. Part 1. understanding the physical and computational complexity of in silico drug design. Acta Crystallogr D Struct Biol 2021; 77:1348-1356. [PMID: 34726163 PMCID: PMC8561735 DOI: 10.1107/s2059798321009712] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/25/2020] [Accepted: 09/17/2021] [Indexed: 02/04/2023] Open
Abstract
The predictive power of simulation has become embedded in the infrastructure of modern economies. Computer-aided design is ubiquitous throughout industry. In aeronautical engineering, built infrastructure and materials manufacturing, simulations are routinely used to compute the performance of potential designs before construction. The ability to predict the behaviour of products is a driver of innovation by reducing the cost barrier to new designs, but also because radically novel ideas can be piloted with relatively little risk. Accurate weather forecasting is essential to guide domestic and military flight paths, and therefore the underpinning simulations are critical enough to have implications for national security. However, in the pharmaceutical and biotechnological industries, the application of computer simulations remains limited by the capabilities of the technology with respect to the complexity of molecular biology and human physiology. Over the last 30 years, molecular-modelling tools have gradually gained a degree of acceptance in the pharmaceutical industry. Drug discovery has begun to benefit from physics-based simulations. While such simulations have great potential for improved molecular design, much scepticism remains about their value. The motivations for such reservations in industry and areas where simulations show promise for efficiency gains in preclinical research are discussed. In this, the first of two complementary papers, the scientific and technical progress that needs to be made to improve the predictive power of biomolecular simulations, and how this might be achieved, is firstly discussed (Part 1). In Part 2, the status of computer simulations in pharma is contrasted with aerodynamics modelling and weather forecasting, and comments are made on the cultural changes needed for equivalent computational technologies to become integrated into life-science industries.
Collapse
Affiliation(s)
- Tom Edwards
- School of Molecular and Cellular Biology, University of Leeds, Leeds, United Kingdom
- Astbury Centre for Structural and Molecular Biology, University of Leeds, Leeds, United Kingdom
| | | | - Sarah Anne Harris
- Astbury Centre for Structural and Molecular Biology, University of Leeds, Leeds, United Kingdom
- School of Physics and Astronomy, University of Leeds, Leeds, United Kingdom
| | - Geoff Wells
- School of Pharmacy, University College London, London, United Kingdom
| |
Collapse
|
11
|
Li C, Voth GA. Accurate and Transferable Reactive Molecular Dynamics Models from Constrained Density Functional Theory. J Phys Chem B 2021; 125:10471-10480. [PMID: 34520198 PMCID: PMC8480781 DOI: 10.1021/acs.jpcb.1c05992] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
![]()
Chemical reactions
constitute the central feature of many liquid,
material, and biomolecular processes. Conventional molecular dynamics
(MD) is inadequate for simulating chemical reactions given the fixed
bonding topology of most force fields, while modeling chemical reactions
using ab initio molecular dynamics is limited to
shorter time and length scales given its high computational cost.
As such, the multiscale reactive molecular dynamics method provides
one promising alternative for simulating complex chemical systems
at atomistic detail on a reactive potential energy surface. However,
the parametrization of such models is a key barrier to their applicability
and success. In this work, we present reactive MD models derived from
constrained density functional theory that are both accurate and transferable.
We illustrate the features of these models for proton dissociation
reactions of amino acids in both aqueous and protein environments.
Specifically, we present models for ionizable glutamate and lysine
that predict accurate absolute pKa values
in water as well as their significantly shifted pKa in staphylococcal nuclease (SNase) without any modification
of the models. As one outcome of the new methodology, the simulations
show that the deprotonation of ionizable residues in SNase can be
closely coupled with side chain rotations, which is a concept likely
generalizable to many other proteins. Furthermore, the present approach
is not limited to only pKa prediction
but can enable the fully atomistic simulation of many other reactive
systems along with a determination of the key aspects of the reaction
mechanisms.
Collapse
Affiliation(s)
- Chenghan Li
- Department of Chemistry, Chicago Center for Theoretical Chemistry, James Franck Institute, and Institute for Biophysical Dynamics, University of Chicago, Chicago, Illinois 60637, United States
| | - Gregory A Voth
- Department of Chemistry, Chicago Center for Theoretical Chemistry, James Franck Institute, and Institute for Biophysical Dynamics, University of Chicago, Chicago, Illinois 60637, United States
| |
Collapse
|
12
|
Barry E, Burns R, Chen W, De Hoe GX, De Oca JMM, de Pablo JJ, Dombrowski J, Elam JW, Felts AM, Galli G, Hack J, He Q, He X, Hoenig E, Iscen A, Kash B, Kung HH, Lewis NHC, Liu C, Ma X, Mane A, Martinson ABF, Mulfort KL, Murphy J, Mølhave K, Nealey P, Qiao Y, Rozyyev V, Schatz GC, Sibener SJ, Talapin D, Tiede DM, Tirrell MV, Tokmakoff A, Voth GA, Wang Z, Ye Z, Yesibolati M, Zaluzec NJ, Darling SB. Advanced Materials for Energy-Water Systems: The Central Role of Water/Solid Interfaces in Adsorption, Reactivity, and Transport. Chem Rev 2021; 121:9450-9501. [PMID: 34213328 DOI: 10.1021/acs.chemrev.1c00069] [Citation(s) in RCA: 19] [Impact Index Per Article: 6.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
The structure, chemistry, and charge of interfaces between materials and aqueous fluids play a central role in determining properties and performance of numerous water systems. Sensors, membranes, sorbents, and heterogeneous catalysts almost uniformly rely on specific interactions between their surfaces and components dissolved or suspended in the water-and often the water molecules themselves-to detect and mitigate contaminants. Deleterious processes in these systems such as fouling, scaling (inorganic deposits), and corrosion are also governed by interfacial phenomena. Despite the importance of these interfaces, much remains to be learned about their multiscale interactions. Developing a deeper understanding of the molecular- and mesoscale phenomena at water/solid interfaces will be essential to driving innovation to address grand challenges in supplying sufficient fit-for-purpose water in the future. In this Review, we examine the current state of knowledge surrounding adsorption, reactivity, and transport in several key classes of water/solid interfaces, drawing on a synergistic combination of theory, simulation, and experiments, and provide an outlook for prioritizing strategic research directions.
Collapse
Affiliation(s)
- Edward Barry
- Advanced Materials for Energy-Water Systems (AMEWS) Energy Frontier Research Center (EFRC), Argonne National Laboratory, 9700 South Cass Avenue, Lemont, Illinois 60439 United States.,Applied Materials Division, Argonne National Laboratory, 9700 South Cass Avenue, Lemont, Illinois 60439 United States.,Center for Molecular Engineering, Argonne National Laboratory, 9700 South Cass Avenue, Lemont, Illinois 60439 United States
| | - Raelyn Burns
- Advanced Materials for Energy-Water Systems (AMEWS) Energy Frontier Research Center (EFRC), Argonne National Laboratory, 9700 South Cass Avenue, Lemont, Illinois 60439 United States.,Applied Materials Division, Argonne National Laboratory, 9700 South Cass Avenue, Lemont, Illinois 60439 United States
| | - Wei Chen
- Advanced Materials for Energy-Water Systems (AMEWS) Energy Frontier Research Center (EFRC), Argonne National Laboratory, 9700 South Cass Avenue, Lemont, Illinois 60439 United States.,Materials Science Division, Argonne National Laboratory, 9700 South Cass Avenue, Lemont, Illinois 60439 United States.,Center for Molecular Engineering, Argonne National Laboratory, 9700 South Cass Avenue, Lemont, Illinois 60439 United States.,Pritzker School of Molecular Engineering, University of Chicago, 5640 South Ellis Avenue, Chicago, Illinois 60637 United States
| | - Guilhem X De Hoe
- Advanced Materials for Energy-Water Systems (AMEWS) Energy Frontier Research Center (EFRC), Argonne National Laboratory, 9700 South Cass Avenue, Lemont, Illinois 60439 United States.,Materials Science Division, Argonne National Laboratory, 9700 South Cass Avenue, Lemont, Illinois 60439 United States.,Center for Molecular Engineering, Argonne National Laboratory, 9700 South Cass Avenue, Lemont, Illinois 60439 United States.,Pritzker School of Molecular Engineering, University of Chicago, 5640 South Ellis Avenue, Chicago, Illinois 60637 United States
| | - Joan Manuel Montes De Oca
- Advanced Materials for Energy-Water Systems (AMEWS) Energy Frontier Research Center (EFRC), Argonne National Laboratory, 9700 South Cass Avenue, Lemont, Illinois 60439 United States.,Pritzker School of Molecular Engineering, University of Chicago, 5640 South Ellis Avenue, Chicago, Illinois 60637 United States
| | - Juan J de Pablo
- Advanced Materials for Energy-Water Systems (AMEWS) Energy Frontier Research Center (EFRC), Argonne National Laboratory, 9700 South Cass Avenue, Lemont, Illinois 60439 United States.,Materials Science Division, Argonne National Laboratory, 9700 South Cass Avenue, Lemont, Illinois 60439 United States.,Pritzker School of Molecular Engineering, University of Chicago, 5640 South Ellis Avenue, Chicago, Illinois 60637 United States
| | - James Dombrowski
- Advanced Materials for Energy-Water Systems (AMEWS) Energy Frontier Research Center (EFRC), Argonne National Laboratory, 9700 South Cass Avenue, Lemont, Illinois 60439 United States.,Department of Chemical and Biological Engineering, Northwestern University, 2145 Sheridan Road, Evanston, Illinois 60208 United States
| | - Jeffrey W Elam
- Advanced Materials for Energy-Water Systems (AMEWS) Energy Frontier Research Center (EFRC), Argonne National Laboratory, 9700 South Cass Avenue, Lemont, Illinois 60439 United States.,Applied Materials Division, Argonne National Laboratory, 9700 South Cass Avenue, Lemont, Illinois 60439 United States.,Center for Molecular Engineering, Argonne National Laboratory, 9700 South Cass Avenue, Lemont, Illinois 60439 United States
| | - Alanna M Felts
- Advanced Materials for Energy-Water Systems (AMEWS) Energy Frontier Research Center (EFRC), Argonne National Laboratory, 9700 South Cass Avenue, Lemont, Illinois 60439 United States.,Department of Chemistry, Northwestern University, 2145 Sheridan Road, Evanston, Illinois 60208 United States
| | - Giulia Galli
- Advanced Materials for Energy-Water Systems (AMEWS) Energy Frontier Research Center (EFRC), Argonne National Laboratory, 9700 South Cass Avenue, Lemont, Illinois 60439 United States.,Materials Science Division, Argonne National Laboratory, 9700 South Cass Avenue, Lemont, Illinois 60439 United States.,Pritzker School of Molecular Engineering, University of Chicago, 5640 South Ellis Avenue, Chicago, Illinois 60637 United States
| | - John Hack
- Advanced Materials for Energy-Water Systems (AMEWS) Energy Frontier Research Center (EFRC), Argonne National Laboratory, 9700 South Cass Avenue, Lemont, Illinois 60439 United States.,Department of Chemistry, University of Chicago, 5735 South Ellis Avenue, Chicago, Illinois 60637 United States
| | - Qiming He
- Advanced Materials for Energy-Water Systems (AMEWS) Energy Frontier Research Center (EFRC), Argonne National Laboratory, 9700 South Cass Avenue, Lemont, Illinois 60439 United States.,Materials Science Division, Argonne National Laboratory, 9700 South Cass Avenue, Lemont, Illinois 60439 United States.,Center for Molecular Engineering, Argonne National Laboratory, 9700 South Cass Avenue, Lemont, Illinois 60439 United States.,Pritzker School of Molecular Engineering, University of Chicago, 5640 South Ellis Avenue, Chicago, Illinois 60637 United States
| | - Xiang He
- Advanced Materials for Energy-Water Systems (AMEWS) Energy Frontier Research Center (EFRC), Argonne National Laboratory, 9700 South Cass Avenue, Lemont, Illinois 60439 United States.,Chemical Sciences and Engineering Division, Argonne National Laboratory, 9700 South Cass Avenue, Lemont, Illinois 60439 United States
| | - Eli Hoenig
- Advanced Materials for Energy-Water Systems (AMEWS) Energy Frontier Research Center (EFRC), Argonne National Laboratory, 9700 South Cass Avenue, Lemont, Illinois 60439 United States.,Pritzker School of Molecular Engineering, University of Chicago, 5640 South Ellis Avenue, Chicago, Illinois 60637 United States
| | - Aysenur Iscen
- Advanced Materials for Energy-Water Systems (AMEWS) Energy Frontier Research Center (EFRC), Argonne National Laboratory, 9700 South Cass Avenue, Lemont, Illinois 60439 United States.,Department of Chemical and Biological Engineering, Northwestern University, 2145 Sheridan Road, Evanston, Illinois 60208 United States
| | - Benjamin Kash
- Pritzker School of Molecular Engineering, University of Chicago, 5640 South Ellis Avenue, Chicago, Illinois 60637 United States
| | - Harold H Kung
- Advanced Materials for Energy-Water Systems (AMEWS) Energy Frontier Research Center (EFRC), Argonne National Laboratory, 9700 South Cass Avenue, Lemont, Illinois 60439 United States.,Department of Chemical and Biological Engineering, Northwestern University, 2145 Sheridan Road, Evanston, Illinois 60208 United States
| | - Nicholas H C Lewis
- Advanced Materials for Energy-Water Systems (AMEWS) Energy Frontier Research Center (EFRC), Argonne National Laboratory, 9700 South Cass Avenue, Lemont, Illinois 60439 United States.,Department of Chemistry, University of Chicago, 5735 South Ellis Avenue, Chicago, Illinois 60637 United States
| | - Chong Liu
- Advanced Materials for Energy-Water Systems (AMEWS) Energy Frontier Research Center (EFRC), Argonne National Laboratory, 9700 South Cass Avenue, Lemont, Illinois 60439 United States.,Pritzker School of Molecular Engineering, University of Chicago, 5640 South Ellis Avenue, Chicago, Illinois 60637 United States
| | - Xinyou Ma
- Advanced Materials for Energy-Water Systems (AMEWS) Energy Frontier Research Center (EFRC), Argonne National Laboratory, 9700 South Cass Avenue, Lemont, Illinois 60439 United States.,Department of Chemistry, University of Chicago, 5735 South Ellis Avenue, Chicago, Illinois 60637 United States
| | - Anil Mane
- Advanced Materials for Energy-Water Systems (AMEWS) Energy Frontier Research Center (EFRC), Argonne National Laboratory, 9700 South Cass Avenue, Lemont, Illinois 60439 United States.,Applied Materials Division, Argonne National Laboratory, 9700 South Cass Avenue, Lemont, Illinois 60439 United States
| | - Alex B F Martinson
- Advanced Materials for Energy-Water Systems (AMEWS) Energy Frontier Research Center (EFRC), Argonne National Laboratory, 9700 South Cass Avenue, Lemont, Illinois 60439 United States.,Materials Science Division, Argonne National Laboratory, 9700 South Cass Avenue, Lemont, Illinois 60439 United States.,Center for Molecular Engineering, Argonne National Laboratory, 9700 South Cass Avenue, Lemont, Illinois 60439 United States
| | - Karen L Mulfort
- Advanced Materials for Energy-Water Systems (AMEWS) Energy Frontier Research Center (EFRC), Argonne National Laboratory, 9700 South Cass Avenue, Lemont, Illinois 60439 United States.,Chemical Sciences and Engineering Division, Argonne National Laboratory, 9700 South Cass Avenue, Lemont, Illinois 60439 United States
| | - Julia Murphy
- Advanced Materials for Energy-Water Systems (AMEWS) Energy Frontier Research Center (EFRC), Argonne National Laboratory, 9700 South Cass Avenue, Lemont, Illinois 60439 United States.,Department of Chemistry, University of Chicago, 5735 South Ellis Avenue, Chicago, Illinois 60637 United States
| | - Kristian Mølhave
- Advanced Materials for Energy-Water Systems (AMEWS) Energy Frontier Research Center (EFRC), Argonne National Laboratory, 9700 South Cass Avenue, Lemont, Illinois 60439 United States.,Technical University of Denmark, Anker Engelunds Vej 1 Bygning 101A, Kgs. Lyngby, Lyngby, Hovedstaden 2800, DK Denmark
| | - Paul Nealey
- Advanced Materials for Energy-Water Systems (AMEWS) Energy Frontier Research Center (EFRC), Argonne National Laboratory, 9700 South Cass Avenue, Lemont, Illinois 60439 United States.,Materials Science Division, Argonne National Laboratory, 9700 South Cass Avenue, Lemont, Illinois 60439 United States.,Pritzker School of Molecular Engineering, University of Chicago, 5640 South Ellis Avenue, Chicago, Illinois 60637 United States
| | - Yijun Qiao
- Advanced Materials for Energy-Water Systems (AMEWS) Energy Frontier Research Center (EFRC), Argonne National Laboratory, 9700 South Cass Avenue, Lemont, Illinois 60439 United States.,Materials Science Division, Argonne National Laboratory, 9700 South Cass Avenue, Lemont, Illinois 60439 United States.,Center for Molecular Engineering, Argonne National Laboratory, 9700 South Cass Avenue, Lemont, Illinois 60439 United States
| | - Vepa Rozyyev
- Advanced Materials for Energy-Water Systems (AMEWS) Energy Frontier Research Center (EFRC), Argonne National Laboratory, 9700 South Cass Avenue, Lemont, Illinois 60439 United States.,Applied Materials Division, Argonne National Laboratory, 9700 South Cass Avenue, Lemont, Illinois 60439 United States
| | - George C Schatz
- Advanced Materials for Energy-Water Systems (AMEWS) Energy Frontier Research Center (EFRC), Argonne National Laboratory, 9700 South Cass Avenue, Lemont, Illinois 60439 United States.,Department of Chemistry, Northwestern University, 2145 Sheridan Road, Evanston, Illinois 60208 United States
| | - Steven J Sibener
- Advanced Materials for Energy-Water Systems (AMEWS) Energy Frontier Research Center (EFRC), Argonne National Laboratory, 9700 South Cass Avenue, Lemont, Illinois 60439 United States.,Department of Chemistry, University of Chicago, 5735 South Ellis Avenue, Chicago, Illinois 60637 United States
| | - Dmitri Talapin
- Advanced Materials for Energy-Water Systems (AMEWS) Energy Frontier Research Center (EFRC), Argonne National Laboratory, 9700 South Cass Avenue, Lemont, Illinois 60439 United States.,Department of Chemistry, University of Chicago, 5735 South Ellis Avenue, Chicago, Illinois 60637 United States
| | - David M Tiede
- Advanced Materials for Energy-Water Systems (AMEWS) Energy Frontier Research Center (EFRC), Argonne National Laboratory, 9700 South Cass Avenue, Lemont, Illinois 60439 United States.,Chemical Sciences and Engineering Division, Argonne National Laboratory, 9700 South Cass Avenue, Lemont, Illinois 60439 United States
| | - Matthew V Tirrell
- Advanced Materials for Energy-Water Systems (AMEWS) Energy Frontier Research Center (EFRC), Argonne National Laboratory, 9700 South Cass Avenue, Lemont, Illinois 60439 United States.,Materials Science Division, Argonne National Laboratory, 9700 South Cass Avenue, Lemont, Illinois 60439 United States.,Center for Molecular Engineering, Argonne National Laboratory, 9700 South Cass Avenue, Lemont, Illinois 60439 United States.,Pritzker School of Molecular Engineering, University of Chicago, 5640 South Ellis Avenue, Chicago, Illinois 60637 United States
| | - Andrei Tokmakoff
- Advanced Materials for Energy-Water Systems (AMEWS) Energy Frontier Research Center (EFRC), Argonne National Laboratory, 9700 South Cass Avenue, Lemont, Illinois 60439 United States.,Department of Chemistry, University of Chicago, 5735 South Ellis Avenue, Chicago, Illinois 60637 United States
| | - Gregory A Voth
- Advanced Materials for Energy-Water Systems (AMEWS) Energy Frontier Research Center (EFRC), Argonne National Laboratory, 9700 South Cass Avenue, Lemont, Illinois 60439 United States.,Department of Chemistry, University of Chicago, 5735 South Ellis Avenue, Chicago, Illinois 60637 United States
| | - Zhongyang Wang
- Advanced Materials for Energy-Water Systems (AMEWS) Energy Frontier Research Center (EFRC), Argonne National Laboratory, 9700 South Cass Avenue, Lemont, Illinois 60439 United States.,Pritzker School of Molecular Engineering, University of Chicago, 5640 South Ellis Avenue, Chicago, Illinois 60637 United States
| | - Zifan Ye
- Advanced Materials for Energy-Water Systems (AMEWS) Energy Frontier Research Center (EFRC), Argonne National Laboratory, 9700 South Cass Avenue, Lemont, Illinois 60439 United States.,Pritzker School of Molecular Engineering, University of Chicago, 5640 South Ellis Avenue, Chicago, Illinois 60637 United States
| | - Murat Yesibolati
- Advanced Materials for Energy-Water Systems (AMEWS) Energy Frontier Research Center (EFRC), Argonne National Laboratory, 9700 South Cass Avenue, Lemont, Illinois 60439 United States.,Technical University of Denmark, Anker Engelunds Vej 1 Bygning 101A, Kgs. Lyngby, Lyngby, Hovedstaden 2800, DK Denmark
| | - Nestor J Zaluzec
- Advanced Materials for Energy-Water Systems (AMEWS) Energy Frontier Research Center (EFRC), Argonne National Laboratory, 9700 South Cass Avenue, Lemont, Illinois 60439 United States.,Photon Sciences Directorate, Argonne National Laboratory, 9700 South Cass Avenue, Lemont, Illinois 60439 United States
| | - Seth B Darling
- Advanced Materials for Energy-Water Systems (AMEWS) Energy Frontier Research Center (EFRC), Argonne National Laboratory, 9700 South Cass Avenue, Lemont, Illinois 60439 United States.,Center for Molecular Engineering, Argonne National Laboratory, 9700 South Cass Avenue, Lemont, Illinois 60439 United States.,Chemical Sciences and Engineering Division, Argonne National Laboratory, 9700 South Cass Avenue, Lemont, Illinois 60439 United States.,Pritzker School of Molecular Engineering, University of Chicago, 5640 South Ellis Avenue, Chicago, Illinois 60637 United States
| |
Collapse
|
13
|
Watkins LC, DeGrado WF, Voth GA. Influenza A M2 Inhibitor Binding Understood through Mechanisms of Excess Proton Stabilization and Channel Dynamics. J Am Chem Soc 2020; 142:17425-17433. [PMID: 32933245 PMCID: PMC7564090 DOI: 10.1021/jacs.0c06419] [Citation(s) in RCA: 16] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022]
Abstract
![]()
Prevalent resistance to inhibitors
that target the influenza A
M2 proton channel has necessitated a continued drug design effort,
supported by a sustained study of the mechanism of channel function
and inhibition. Recent high-resolution X-ray crystal structures present
the first opportunity to see how the adamantyl amine class of inhibitors
bind to M2 and disrupt and interact with the channel’s water
network, providing insight into the critical properties that enable
their effective inhibition in wild-type M2. In this work, we examine
the hypothesis that these drugs act primarily as mechanism-based inhibitors
by comparing hydrated excess proton stabilization during proton transport
in M2 with the interactions revealed in the crystal structures, using
the Multiscale Reactive Molecular Dynamics (MS-RMD) methodology. MS-RMD,
unlike classical molecular dynamics, models the hydrated proton (hydronium-like
cation) as a dynamic excess charge defect and allows bonds to break
and form, capturing the intricate interactions between the hydrated
excess proton, protein atoms, and water. Through this, we show that
the ammonium group of the inhibitors is effectively positioned to
take advantage of the channel’s natural ability to stabilize
an excess protonic charge and act as a hydronium mimic. Additionally,
we show that the channel is especially stable in the drug binding
region, highlighting the importance of this property for binding the
adamantane group. Finally, we characterize an additional hinge point
near Val27, which dynamically responds to charge and inhibitor binding.
Altogether, this work further illuminates a dynamic understanding
of the mechanism of drug inhibition in M2, grounded in the fundamental
properties that enable the channel to transport and stabilize excess
protons, with critical implications for future drug design efforts.
Collapse
Affiliation(s)
- Laura C Watkins
- Department of Chemistry, Chicago Center for Theoretical Chemistry, Institute for Biophysical Dynamics and James Franck Institute, The University of Chicago, Chicago, Illinois 60637, United States
| | - William F DeGrado
- Department of Pharmaceutical Chemistry, University of California, San Francisco, San Francisco, California 94158, United States
| | - Gregory A Voth
- Department of Chemistry, Chicago Center for Theoretical Chemistry, Institute for Biophysical Dynamics and James Franck Institute, The University of Chicago, Chicago, Illinois 60637, United States
| |
Collapse
|
14
|
Scivetti I, Sen K, Elena AM, Todorov I. Reactive Molecular Dynamics at Constant Pressure via Nonreactive Force Fields: Extending the Empirical Valence Bond Method to the Isothermal-Isobaric Ensemble. J Phys Chem A 2020; 124:7585-7597. [PMID: 32820921 DOI: 10.1021/acs.jpca.0c05461] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
The Empirical Valence Bond (EVB) method offers a suitable framework to obtain reactive potentials through the coupling of nonreactive force fields. In this formalism, most of the implemented coupling terms are built using functional forms that depend on spatial coordinates, while parameters are fitted against reference data to model the change of chemistry between the participating nonreactive states. In this work, we demonstrate that the use of such coupling terms precludes the computation of the stress tensor for condensed phase systems and prevents the possibility to carry out EVB molecular dynamics in the isothermal-isobaric (NPT) ensemble. Alternatively, we make use of coupling terms that depend on the energy gaps, defined as the energy differences between the participating nonreactive force fields, and derive a general expression for the EVB stress tensor suitable for computation. Implementation of this new methodology is tested for a model of a single reactive malonaldehyde solvated in nonreactive water. Mass densities and probability distributions for the values of the energy gaps computed in the NPT ensemble reveal a negligible role of the reactive potential in the limit of low concentrated solutions, thus corroborating for the first time the validity of approximations based on the canonical NVT ensemble, customarily adopted for EVB simulations. The presented formalism also aims to contribute to future implementations and extensions of the EVB method to research the limit of highly concentrated solutions.
Collapse
Affiliation(s)
- Ivan Scivetti
- Daresbury Laboratory, Sc. Tech., Keckwick Lane, Daresbury, Warrington WA4 4AD, U.K.,Department of Chemistry, University of Liverpool, Liverpool L69 3BX, U.K
| | - Kakali Sen
- Daresbury Laboratory, Sc. Tech., Keckwick Lane, Daresbury, Warrington WA4 4AD, U.K
| | - Alin M Elena
- Daresbury Laboratory, Sc. Tech., Keckwick Lane, Daresbury, Warrington WA4 4AD, U.K
| | - Ilian Todorov
- Daresbury Laboratory, Sc. Tech., Keckwick Lane, Daresbury, Warrington WA4 4AD, U.K
| |
Collapse
|
15
|
Khalid S, Newstead S. A Computational Swiss Army Knife Approach to Unraveling the Secrets of Proton Movement through SERCA. Biophys J 2020; 119:890-891. [DOI: 10.1016/j.bpj.2020.07.028] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/22/2020] [Accepted: 07/30/2020] [Indexed: 10/23/2022] Open
|
16
|
Li C, Yue Z, Espinoza-Fonseca LM, Voth GA. Multiscale Simulation Reveals Passive Proton Transport Through SERCA on the Microsecond Timescale. Biophys J 2020; 119:1033-1040. [PMID: 32814059 DOI: 10.1016/j.bpj.2020.07.027] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/17/2020] [Revised: 07/20/2020] [Accepted: 07/24/2020] [Indexed: 12/11/2022] Open
Abstract
The sarcoplasmic reticulum Ca2+-ATPase (SERCA) transports two Ca2+ ions from the cytoplasm to the reticulum lumen at the expense of ATP hydrolysis. In addition to transporting Ca2+, SERCA facilitates bidirectional proton transport across the sarcoplasmic reticulum to maintain the charge balance of the transport sites and to balance the charge deficit generated by the exchange of Ca2+. Previous studies have shown the existence of a transient water-filled pore in SERCA that connects the Ca2+ binding sites with the lumen, but the capacity of this pathway to sustain passive proton transport has remained unknown. In this study, we used the multiscale reactive molecular dynamics method and free energy sampling to quantify the free energy profile and timescale of the proton transport across this pathway while also explicitly accounting for the dynamically coupled hydration changes of the pore. We find that proton transport from the central binding site to the lumen has a microsecond timescale, revealing a novel passive cytoplasm-to-lumen proton flow beside the well-known inverse proton countertransport occurring in active Ca2+ transport. We propose that this proton transport mechanism is operational and serves as a functional conduit for passive proton transport across the sarcoplasmic reticulum.
Collapse
Affiliation(s)
- Chenghan Li
- Department of Chemistry, Chicago Center for Theoretical Chemistry, James Franck Institute, and Institute for Biophysical Dynamics, The University of Chicago, Chicago, Illinois
| | - Zhi Yue
- Department of Chemistry, Chicago Center for Theoretical Chemistry, James Franck Institute, and Institute for Biophysical Dynamics, The University of Chicago, Chicago, Illinois
| | - L Michel Espinoza-Fonseca
- Center for Arrhythmia Research, Department of Internal Medicine, Division of Cardiovascular Medicine, University of Michigan, Ann Arbor, Michigan
| | - Gregory A Voth
- Department of Chemistry, Chicago Center for Theoretical Chemistry, James Franck Institute, and Institute for Biophysical Dynamics, The University of Chicago, Chicago, Illinois.
| |
Collapse
|
17
|
Li C, Swanson JMJ. Understanding and Tracking the Excess Proton in Ab Initio Simulations; Insights from IR Spectra. J Phys Chem B 2020; 124:5696-5708. [PMID: 32515957 PMCID: PMC7448536 DOI: 10.1021/acs.jpcb.0c03615] [Citation(s) in RCA: 13] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/07/2023]
Abstract
Proton transport in aqueous media is ubiquitously important in chemical and biological processes. Although ab initio molecular dynamics (AIMD) simulations have made great progress in characterizing proton transport, there has been a long-standing challenge in defining and tracking the excess proton, or more properly, the center of excess charge (CEC) created when a hydrogen nucleus distorts the electron distributions of water molecules in a delocalized and highly dynamic nature. Yet, defining (and biasing) such a CEC is essential when combining AIMD with enhanced sampling methods to calculate the relevant macroscopic properties via free-energy landscapes, which is the standard practice for most processes of interest. Several CEC formulas have been proposed and used, but none have yet been systematically tested or rigorously derived. In this paper, we show that the CEC can be used as a computational tool to disentangle IR features of the solvated excess proton from its surrounding solvent, and in turn, how correlating the features in the excess charge spectrum with the behavior of CEC in simulations enables a systematic evaluation of various CEC definitions. We present a new definition of CEC and show how it overcomes the limitations of those currently available both from a spectroscopic point of view and from a practical perspective of performance in enhanced sampling simulations.
Collapse
Affiliation(s)
- Chenghan Li
- Department of Chemistry, Chicago Center for Theoretical Chemistry, James Franck Institute, and Institute for Biophysical Dynamics, The University of Chicago, Chicago, Illinois 60637, United States
| | - Jessica M. J. Swanson
- Department of Chemistry, Biological Chemistry Program, and Center for Cell and Genome Science, The University of Utah, Salt Lake City, Utah 84112, United States
| |
Collapse
|
18
|
Li Z, Li C, Wang Z, Voth GA. What Coordinate Best Describes the Affinity of the Hydrated Excess Proton for the Air-Water Interface? J Phys Chem B 2020; 124:5039-5046. [PMID: 32426982 DOI: 10.1021/acs.jpcb.0c03288] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Molecular dynamics simulations and free energy sampling are employed in this work to investigate the surface affinity of the hydrated excess proton with two definitions of the interface: The Gibbs dividing interface (GDI) and the Willard-Chandler interface (WCI). Both the multistate empirical valence bond (MS-EVB) reactive molecular dynamics method and the density functional theory-based ab initio molecular dynamics (AIMD) were used to describe the hydrated excess proton species, including "vehicular" (standard diffusion) transport and (Grotthuss) proton hopping transport and associated structures of the hydrated excess proton net positive charge defect. The excess proton is found to exhibit a similar trend and quantitative free energy behavior in terms of its surface affinity as a function of the GDI or WCI. Importantly, the definitions of the two interfaces in terms of the excess proton charge defect are highly correlated and far from independent of one another, thus undermining the argument that one interface is superior to the other when describing the proton interface affinity. Moreover, the hydrated excess proton and its solvation shell significantly influence the location and local curvature of the WCI, making it difficult to disentangle the interfacial thermodynamics of the excess proton from the influence of that species on the instantaneous surface curvature.
Collapse
Affiliation(s)
- Zhefu Li
- Department of Chemistry, Chicago Center for Theoretical Chemistry, James Franck Institute, and Institute of Biophysics Dynamics, The University of Chicago, 5735 South Ellis Avenue, Chicago, Illinois 60637, United States
| | - Chenghan Li
- Department of Chemistry, Chicago Center for Theoretical Chemistry, James Franck Institute, and Institute of Biophysics Dynamics, The University of Chicago, 5735 South Ellis Avenue, Chicago, Illinois 60637, United States
| | - Zhi Wang
- Department of Chemistry, Chicago Center for Theoretical Chemistry, James Franck Institute, and Institute of Biophysics Dynamics, The University of Chicago, 5735 South Ellis Avenue, Chicago, Illinois 60637, United States
| | - Gregory A Voth
- Department of Chemistry, Chicago Center for Theoretical Chemistry, James Franck Institute, and Institute of Biophysics Dynamics, The University of Chicago, 5735 South Ellis Avenue, Chicago, Illinois 60637, United States
| |
Collapse
|
19
|
Wang Z, Swanson JMJ, Voth GA. Local conformational dynamics regulating transport properties of a Cl - /H + antiporter. J Comput Chem 2020; 41:513-519. [PMID: 31633205 PMCID: PMC7184886 DOI: 10.1002/jcc.26093] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/08/2019] [Revised: 09/17/2019] [Accepted: 10/03/2019] [Indexed: 11/08/2022]
Abstract
ClC-ec1 is a Cl- /H+ antiporter that exchanges Cl- and H+ ions across the membrane. Experiments have demonstrated that several mutations, including I109F, decrease the Cl- and H+ transport rates by an order of magnitude. Using reactive molecular dynamics simulations of explicit proton transport across the central region in the I109F mutant, a two-dimensional free energy profile has been constructed that is consistent with the experimental transport rates. The importance of a phenylalanine gate formed by F109 and F357 and its influence on hydration connectivity through the central proton transport pathway is revealed. This work demonstrates how seemingly subtle changes in local conformational dynamics can dictate hydration changes and thus transport properties. © 2019 Wiley Periodicals, Inc.
Collapse
Affiliation(s)
- Zhi Wang
- Department of Chemistry, James Franck Institute, and Institute for Biophysical Dynamics, The University of Chicago, Chicago, Illinois 60637, United States
| | | | - Gregory A. Voth
- Department of Chemistry, James Franck Institute, and Institute for Biophysical Dynamics, The University of Chicago, Chicago, Illinois 60637, United States
| |
Collapse
|
20
|
Dannenhoffer-Lafage T, Voth GA. Reactive Coarse-Grained Molecular Dynamics. J Chem Theory Comput 2020; 16:2541-2549. [DOI: 10.1021/acs.jctc.9b01140] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/20/2023]
Affiliation(s)
- Thomas Dannenhoffer-Lafage
- Department of Chemistry, Chicago Center for Theoretical Chemistry, James Franck Institute, and Institute for Biophysical Dynamics, The University of Chicago, 5735 South Ellis Avenue, Chicago, Illinois 60637, United States
| | - Gregory A. Voth
- Department of Chemistry, Chicago Center for Theoretical Chemistry, James Franck Institute, and Institute for Biophysical Dynamics, The University of Chicago, 5735 South Ellis Avenue, Chicago, Illinois 60637, United States
| |
Collapse
|
21
|
Yue Z, Li C, Voth GA, Swanson JMJ. Dynamic Protonation Dramatically Affects the Membrane Permeability of Drug-like Molecules. J Am Chem Soc 2019; 141:13421-13433. [PMID: 31382734 DOI: 10.1021/jacs.9b04387] [Citation(s) in RCA: 44] [Impact Index Per Article: 8.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/06/2023]
Abstract
Permeability (Pm) across biological membranes is of fundamental importance and a key factor in drug absorption, distribution, and development. Although the majority of drugs will be charged at some point during oral delivery, our understanding of membrane permeation by charged species is limited. The canonical model assumes that only neutral molecules partition into and passively permeate across membranes, but there is mounting evidence that these processes are also facile for certain charged species. However, it is unknown whether such ionizable permeants dynamically neutralize at the membrane surface or permeate in their charged form. To probe protonation-coupled permeation in atomic detail, we herein apply continuous constant-pH molecular dynamics along with free energy sampling to study the permeation of a weak base propranolol (PPL), and evaluate the impact of including dynamic protonation on Pm. The simulations reveal that PPL dynamically neutralizes at the lipid-tail interface, which dramatically influences the permeation free energy landscape and explains why the conventional model overestimates the assigned intrinsic permeability. We demonstrate how fixed-charge-state simulations can account for this effect, and propose a revised model that better describes pH-coupled partitioning and permeation. Our results demonstrate how dynamic changes in protonation state may play a critical role in the permeation of ionizable molecules, including pharmaceuticals and drug-like molecules, thus requiring a revision of the standard picture.
Collapse
Affiliation(s)
- Zhi Yue
- Department of Chemistry, James Frank Institute, and Institute for Biophysical Dynamics , The University of Chicago , Chicago , Illinois 60637 , United States
| | - Chenghan Li
- Department of Chemistry, James Frank Institute, and Institute for Biophysical Dynamics , The University of Chicago , Chicago , Illinois 60637 , United States
| | - Gregory A Voth
- Department of Chemistry, James Frank Institute, and Institute for Biophysical Dynamics , The University of Chicago , Chicago , Illinois 60637 , United States
| | - Jessica M J Swanson
- Department of Chemistry, James Frank Institute, and Institute for Biophysical Dynamics , The University of Chicago , Chicago , Illinois 60637 , United States
| |
Collapse
|
22
|
Watkins LC, Liang R, Swanson JMJ, DeGrado WF, Voth GA. Proton-Induced Conformational and Hydration Dynamics in the Influenza A M2 Channel. J Am Chem Soc 2019; 141:11667-11676. [PMID: 31264413 DOI: 10.1021/jacs.9b05136] [Citation(s) in RCA: 21] [Impact Index Per Article: 4.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022]
Abstract
The influenza A M2 protein is an acid-activated proton channel responsible for acidification of the inside of the virus, a critical step in the viral life cycle. This channel has four central histidine residues that form an acid-activated gate, binding protons from the outside until an activated state allows proton transport to the inside. While previous work has focused on proton transport through the channel, the structural and dynamic changes that accompany proton flux and enable activation have yet to be resolved. In this study, extensive Multiscale Reactive Molecular Dynamics simulations with explicit Grotthuss-shuttling hydrated excess protons are used to explore detailed molecular-level interactions that accompany proton transport in the +0, + 1, and +2 histidine charge states. The results demonstrate how the hydrated excess proton strongly influences both the protein and water hydrogen-bonding network throughout the channel, providing further insight into the channel's acid-activation mechanism and rectification behavior. We find that the excess proton dynamically, as a function of location, shifts the protein structure away from its equilibrium distributions uniquely for different pH conditions consistent with acid-activation. The proton distribution in the xy-plane is also shown to be asymmetric about the channel's main axis, which has potentially important implications for the mechanism of proton conduction and future drug design efforts.
Collapse
Affiliation(s)
- Laura C Watkins
- Department of Chemistry, Institute for Biophysical Dynamics and James Franck Institute , The University of Chicago , Chicago , Illinois 60637 , United States
| | - Ruibin Liang
- Department of Chemistry, Institute for Biophysical Dynamics and James Franck Institute , The University of Chicago , Chicago , Illinois 60637 , United States
| | - Jessica M J Swanson
- Department of Chemistry, Institute for Biophysical Dynamics and James Franck Institute , The University of Chicago , Chicago , Illinois 60637 , United States
| | - William F DeGrado
- Department of Pharmaceutical Chemistry , University of California , San Francisco , California 94158 , United States
| | - Gregory A Voth
- Department of Chemistry, Institute for Biophysical Dynamics and James Franck Institute , The University of Chicago , Chicago , Illinois 60637 , United States
| |
Collapse
|
23
|
Wang Z, Swanson JMJ, Voth GA. Modulating the Chemical Transport Properties of a Transmembrane Antiporter via Alternative Anion Flux. J Am Chem Soc 2018; 140:16535-16543. [PMID: 30421606 PMCID: PMC6379079 DOI: 10.1021/jacs.8b07614] [Citation(s) in RCA: 21] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/16/2023]
Abstract
![]()
ClC-ec1 is a prototype of the ClC
antiporters, proteins that stoichiometrically
exchange Cl– and H+ ions in opposite
directions across a membrane. It has been shown that other polyatomic
anions, such as NO3– and SCN–, can also be transported by ClC-ec1, but with partially or completely
uncoupled proton flux. Herein, with the help of multiscale computer
simulations in which the Grotthuss mechanism of proton transport (PT)
is treated explicitly, we demonstrate how the chemical nature of these
anions alters the coupling mechanism and qualitatively explain the
shifts in the ion stoichiometry. Multidimensional free energy profiles
for PT and the coupled changes in hydration are presented for NO3– and SCN–. The calculated
proton conductances agree with experiment, showing reduced or abolished
proton flux. Surprisingly, the proton affinity of the anion is less
influential on the PT, while its size and interactions with the protein
significantly alter hydration and shift its influence on PT from facilitating
to inhibiting. We find that the hydration of the cavity below the
anion is relatively fast, but connecting the water network past the
steric hindrance of these polyatomic anions is quite slow. Hence,
the most relevant coordinate to the PT free energy barrier is the
water connectivity along the PT pathway, but importantly only in the
presence of the excess proton, and this coordinate is significantly
affected by the nature of the bound anion. This work again demonstrates
how PT is intrinsically coupled with protein cavity hydration changes
as well as influenced by the protein environment. It additionally
suggests ways in which ion exchange can be modulated and exchange
stoichiometries altered.
Collapse
Affiliation(s)
- Zhi Wang
- Department of Chemistry, James Franck Institute, and Institute for Biophysical Dynamics , The University of Chicago , Chicago , Illinois 60637 , United States
| | - Jessica M J Swanson
- Department of Chemistry, James Franck Institute, and Institute for Biophysical Dynamics , The University of Chicago , Chicago , Illinois 60637 , United States
| | - Gregory A Voth
- Department of Chemistry, James Franck Institute, and Institute for Biophysical Dynamics , The University of Chicago , Chicago , Illinois 60637 , United States
| |
Collapse
|
24
|
Abstract
An important limitation of standard classical molecular dynamics simulations is the inability to make or break chemical bonds. This restricts severely our ability to study processes that involve even the simplest of chemical reactions, the transfer of a proton. Existing approaches for allowing proton transfer in the context of classical mechanics are rather cumbersome and have not achieved widespread use and routine status. Here we reconsider the combination of molecular dynamics with periodic stochastic proton hops. To ensure computational efficiency, we propose a non-Boltzmann acceptance criterion that is heuristically adjusted to maintain the correct or desirable thermodynamic equilibria between different protonation states and proton transfer rates. Parameters are proposed for hydronium, Asp, Glu, and His. The algorithm is implemented in the program CHARMM and tested on proton diffusion in bulk water and carbon nanotubes and on proton conductance in the gramicidin A channel. Using hopping parameters determined from proton diffusion in bulk water, the model reproduces the enhanced proton diffusivity in carbon nanotubes and gives a reasonable estimate of the proton conductance in gramicidin A.
Collapse
Affiliation(s)
- Themis Lazaridis
- Department of Chemistry, City College of New York/CUNY , 160 Convent Avenue, New York, New York 10031, United States.,Graduate Programs in Chemistry, Biochemistry & Physics, Graduate Center, City University of New York , 365 Fifth Ave, New York, New York 10016, United States
| | - Gerhard Hummer
- Department of Theoretical Biophysics, Max Planck Institute of Biophysics , Max-von-Laue Strasse 3, 60438 Frankfurt am Main, Germany
| |
Collapse
|
25
|
Affiliation(s)
- Jesse G. McDaniel
- Department of Chemistry, University of Wisconsin, 1101 University Avenue, Madison, Wisconsin 53706, United States
| | - Arun Yethiraj
- Department of Chemistry, University of Wisconsin, 1101 University Avenue, Madison, Wisconsin 53706, United States
| |
Collapse
|
26
|
Liang R, Swanson JMJ, Wikström M, Voth GA. Understanding the essential proton-pumping kinetic gates and decoupling mutations in cytochrome c oxidase. Proc Natl Acad Sci U S A 2017; 114:5924-5929. [PMID: 28536198 PMCID: PMC5468613 DOI: 10.1073/pnas.1703654114] [Citation(s) in RCA: 32] [Impact Index Per Article: 4.6] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022] Open
Abstract
Cytochrome c oxidase (CcO) catalyzes the reduction of oxygen to water and uses the released free energy to pump protons against the transmembrane proton gradient. To better understand the proton-pumping mechanism of the wild-type (WT) CcO, much attention has been given to the mutation of amino acid residues along the proton translocating D-channel that impair, and sometimes decouple, proton pumping from the chemical catalysis. Although their influence has been clearly demonstrated experimentally, the underlying molecular mechanisms of these mutants remain unknown. In this work, we report multiscale reactive molecular dynamics simulations that characterize the free-energy profiles of explicit proton transport through several important D-channel mutants. Our results elucidate the mechanisms by which proton pumping is impaired, thus revealing key kinetic gating features in CcO. In the N139T and N139C mutants, proton back leakage through the D-channel is kinetically favored over proton pumping due to the loss of a kinetic gate in the N139 region. In the N139L mutant, the bulky L139 side chain inhibits timely reprotonation of E286 through the D-channel, which impairs both proton pumping and the chemical reaction. In the S200V/S201V double mutant, the proton affinity of E286 is increased, which slows down both proton pumping and the chemical catalysis. This work thus not only provides insight into the decoupling mechanisms of CcO mutants, but also explains how kinetic gating in the D-channel is imperative to achieving high proton-pumping efficiency in the WT CcO.
Collapse
Affiliation(s)
- Ruibin Liang
- Department of Chemistry, Institute for Biophysical Dynamics and James Franck Institute, The University of Chicago, Chicago, IL 60637
| | - Jessica M J Swanson
- Department of Chemistry, Institute for Biophysical Dynamics and James Franck Institute, The University of Chicago, Chicago, IL 60637;
| | - Mårten Wikström
- Helsinki Bioenergetics Group, Programme for Structural Biology and Biophysics, Institute of Biotechnology, University of Helsinki, FI-00014, Helsinki, Finland
| | - Gregory A Voth
- Department of Chemistry, Institute for Biophysical Dynamics and James Franck Institute, The University of Chicago, Chicago, IL 60637;
| |
Collapse
|
27
|
Sand AM, Truhlar DG, Gagliardi L. Efficient algorithm for multiconfiguration pair-density functional theory with application to the heterolytic dissociation energy of ferrocene. J Chem Phys 2017; 146:034101. [DOI: 10.1063/1.4973709] [Citation(s) in RCA: 26] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/2022] Open
Affiliation(s)
- Andrew M. Sand
- Department of Chemistry, Chemical Theory Center, and Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, USA
| | - Donald G. Truhlar
- Department of Chemistry, Chemical Theory Center, and Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, USA
| | - Laura Gagliardi
- Department of Chemistry, Chemical Theory Center, and Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, USA
| |
Collapse
|
28
|
Multiscale Simulations Reveal Key Aspects of the Proton Transport Mechanism in the ClC-ec1 Antiporter. Biophys J 2016; 110:1334-45. [PMID: 27028643 DOI: 10.1016/j.bpj.2016.02.014] [Citation(s) in RCA: 44] [Impact Index Per Article: 5.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/06/2015] [Revised: 02/05/2016] [Accepted: 02/08/2016] [Indexed: 11/22/2022] Open
Abstract
Multiscale reactive molecular dynamics simulations are used to study proton transport through the central region of ClC-ec1, a widely studied ClC transporter that enables the stoichiometric exchange of 2 Cl(-) ions for 1 proton (H(+)). It has long been known that both Cl(-) and proton transport occur through partially congruent pathways, and that their exchange is strictly coupled. However, the nature of this coupling and the mechanism of antiporting remain topics of debate. Here multiscale simulations have been used to characterize proton transport between E203 (Glu(in)) and E148 (Glu(ex)), the internal and external intermediate proton binding sites, respectively. Free energy profiles are presented, explicitly accounting for the binding of Cl(-) along the central pathway, the dynamically coupled hydration changes of the central region, and conformational changes of Glu(in) and Glu(ex). We find that proton transport between Glu(in) and Glu(ex) is possible in both the presence and absence of Cl(-) in the central binding site, although it is facilitated by the anion presence. These results support the notion that the requisite coupling between Cl(-) and proton transport occurs elsewhere (e.g., during proton uptake or release). In addition, proton transport is explored in the E203K mutant, which maintains proton permeation despite the substitution of a basic residue for Glu(in). This collection of calculations provides for the first time, to our knowledge, a detailed picture of the proton transport mechanism in the central region of ClC-ec1 at a molecular level.
Collapse
|
29
|
Liang R, Swanson JMJ, Madsen JJ, Hong M, DeGrado WF, Voth GA. Acid activation mechanism of the influenza A M2 proton channel. Proc Natl Acad Sci U S A 2016; 113:E6955-E6964. [PMID: 27791184 PMCID: PMC5111692 DOI: 10.1073/pnas.1615471113] [Citation(s) in RCA: 70] [Impact Index Per Article: 8.8] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022] Open
Abstract
The homotetrameric influenza A M2 channel (AM2) is an acid-activated proton channel responsible for the acidification of the influenza virus interior, an important step in the viral lifecycle. Four histidine residues (His37) in the center of the channel act as a pH sensor and proton selectivity filter. Despite intense study, the pH-dependent activation mechanism of the AM2 channel has to date not been completely understood at a molecular level. Herein we have used multiscale computer simulations to characterize (with explicit proton transport free energy profiles and their associated calculated conductances) the activation mechanism of AM2. All proton transfer steps involved in proton diffusion through the channel, including the protonation/deprotonation of His37, are explicitly considered using classical, quantum, and reactive molecular dynamics methods. The asymmetry of the proton transport free energy profile under high-pH conditions qualitatively explains the rectification behavior of AM2 (i.e., why the inward proton flux is allowed when the pH is low in viral exterior and high in viral interior, but outward proton flux is prohibited when the pH gradient is reversed). Also, in agreement with electrophysiological results, our simulations indicate that the C-terminal amphipathic helix does not significantly change the proton conduction mechanism in the AM2 transmembrane domain; the four transmembrane helices flanking the channel lumen alone seem to determine the proton conduction mechanism.
Collapse
Affiliation(s)
- Ruibin Liang
- Department of Chemistry, The University of Chicago, Chicago, IL 60637
- Institute for Biophysical Dynamics, The University of Chicago, Chicago, IL 60637
- James Franck Institute, The University of Chicago, Chicago, IL 60637
| | - Jessica M J Swanson
- Department of Chemistry, The University of Chicago, Chicago, IL 60637
- Institute for Biophysical Dynamics, The University of Chicago, Chicago, IL 60637
- James Franck Institute, The University of Chicago, Chicago, IL 60637
| | - Jesper J Madsen
- Department of Chemistry, The University of Chicago, Chicago, IL 60637
- Institute for Biophysical Dynamics, The University of Chicago, Chicago, IL 60637
- James Franck Institute, The University of Chicago, Chicago, IL 60637
| | - Mei Hong
- Department of Chemistry, Massachusetts Institute of Technology, Cambridge, MA 02139
| | - William F DeGrado
- Department of Pharmaceutical Chemistry, University of San Francisco, San Francisco, CA 94158
| | - Gregory A Voth
- Department of Chemistry, The University of Chicago, Chicago, IL 60637;
- Institute for Biophysical Dynamics, The University of Chicago, Chicago, IL 60637
- James Franck Institute, The University of Chicago, Chicago, IL 60637
| |
Collapse
|
30
|
Lee S, Mayes HB, Swanson JMJ, Voth GA. The Origin of Coupled Chloride and Proton Transport in a Cl -/H + Antiporter. J Am Chem Soc 2016; 138:14923-14930. [PMID: 27783900 PMCID: PMC5114699 DOI: 10.1021/jacs.6b06683] [Citation(s) in RCA: 31] [Impact Index Per Article: 3.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
![]()
The ClC family of transmembrane proteins
functions throughout nature
to control the transport of Cl– ions across biological
membranes. ClC-ec1 from Escherichia coli is an antiporter,
coupling the transport of Cl– and H+ ions
in opposite directions and driven by the concentration gradients of
the ions. Despite keen interest in this protein, the molecular mechanism
of the Cl–/H+ coupling has not been fully
elucidated. Here, we have used multiscale simulation to help identify
the essential mechanism of the Cl–/H+ coupling. We find that the highest barrier for proton transport
(PT) from the intra- to extracellular solution is attributable to
a chemical reaction, the deprotonation of glutamic acid 148 (E148).
This barrier is significantly reduced by the binding of Cl– in the “central” site (Cl–cen), which displaces E148 and thereby facilitates its deprotonation.
Conversely, in the absence of Cl–cen E148
favors the “down” conformation, which results in a much
higher cumulative rotation and deprotonation barrier that effectively
blocks PT to the extracellular solution. Thus, the rotation of E148
plays a critical role in defining the Cl–/H+ coupling. As a control, we have also simulated PT in the
ClC-ec1 E148A mutant to further understand the role of this residue.
Replacement with a non-protonatable residue greatly increases the
free energy barrier for PT from E203 to the extracellular solution,
explaining the experimental result that PT in E148A is blocked whether
or not Cl–cen is present. The results
presented here suggest both how a chemical reaction can control the
rate of PT and also how it can provide a mechanism for a coupling
of the two ion transport processes.
Collapse
Affiliation(s)
- Sangyun Lee
- Department of Chemistry, James Franck Institute, and Institute for Biophysical Dynamics, The University of Chicago , Chicago, Illinois 60637, United States
| | - Heather B Mayes
- Department of Chemistry, James Franck Institute, and Institute for Biophysical Dynamics, The University of Chicago , Chicago, Illinois 60637, United States
| | - Jessica M J Swanson
- Department of Chemistry, James Franck Institute, and Institute for Biophysical Dynamics, The University of Chicago , Chicago, Illinois 60637, United States
| | - Gregory A Voth
- Department of Chemistry, James Franck Institute, and Institute for Biophysical Dynamics, The University of Chicago , Chicago, Illinois 60637, United States
| |
Collapse
|
31
|
Multiscale simulations reveal key features of the proton-pumping mechanism in cytochrome c oxidase. Proc Natl Acad Sci U S A 2016; 113:7420-5. [PMID: 27339133 DOI: 10.1073/pnas.1601982113] [Citation(s) in RCA: 49] [Impact Index Per Article: 6.1] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022] Open
Abstract
Cytochrome c oxidase (CcO) reduces oxygen to water and uses the released free energy to pump protons across the membrane. We have used multiscale reactive molecular dynamics simulations to explicitly characterize (with free-energy profiles and calculated rates) the internal proton transport events that enable proton pumping during first steps of oxidation of the fully reduced enzyme. Our results show that proton transport from amino acid residue E286 to both the pump loading site (PLS) and to the binuclear center (BNC) are thermodynamically driven by electron transfer from heme a to the BNC, but that the former (i.e., pumping) is kinetically favored whereas the latter (i.e., transfer of the chemical proton) is rate-limiting. The calculated rates agree with experimental measurements. The backflow of the pumped proton from the PLS to E286 and from E286 to the inside of the membrane is prevented by large free-energy barriers for the backflow reactions. Proton transport from E286 to the PLS through the hydrophobic cavity and from D132 to E286 through the D-channel are found to be strongly coupled to dynamical hydration changes in the corresponding pathways and, importantly, vice versa.
Collapse
|
32
|
McDaniel JG, Yethiraj A. Importance of hydrophobic traps for proton diffusion in lyotropic liquid crystals. J Chem Phys 2016; 144:094705. [DOI: 10.1063/1.4943131] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Jesse G. McDaniel
- Department of Chemistry, University of Wisconsin, Madison, Wisconsin 53706, USA
| | - Arun Yethiraj
- Department of Chemistry, University of Wisconsin, Madison, Wisconsin 53706, USA
| |
Collapse
|
33
|
Lee S, Liang R, Voth GA, Swanson JMJ. Computationally Efficient Multiscale Reactive Molecular Dynamics to Describe Amino Acid Deprotonation in Proteins. J Chem Theory Comput 2016; 12:879-91. [PMID: 26734942 PMCID: PMC4750100 DOI: 10.1021/acs.jctc.5b01109] [Citation(s) in RCA: 42] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/21/2015] [Indexed: 11/30/2022]
Abstract
An important challenge in the simulation of biomolecular systems is a quantitative description of the protonation and deprotonation process of amino acid residues. Despite the seeming simplicity of adding or removing a positively charged hydrogen nucleus, simulating the actual protonation/deprotonation process is inherently difficult. It requires both the explicit treatment of the excess proton, including its charge defect delocalization and Grotthuss shuttling through inhomogeneous moieties (water and amino residues), and extensive sampling of coupled condensed phase motions. In a recent paper (J. Chem. Theory Comput. 2014, 10, 2729-2737), a multiscale approach was developed to map high-level quantum mechanics/molecular mechanics (QM/MM) data into a multiscale reactive molecular dynamics (MS-RMD) model in order to describe amino acid deprotonation in bulk water. In this article, we extend the fitting approach (called FitRMD) to create MS-RMD models for ionizable amino acids within proteins. The resulting models are shown to faithfully reproduce the free energy profiles of the reference QM/MM Hamiltonian for PT inside an example protein, the ClC-ec1 H(+)/Cl(-) antiporter. Moreover, we show that the resulting MS-RMD models are computationally efficient enough to then characterize more complex 2-dimensional free energy surfaces due to slow degrees of freedom such as water hydration of internal protein cavities that can be inherently coupled to the excess proton charge translocation. The FitRMD method is thus shown to be an effective way to map ab initio level accuracy into a much more computationally efficient reactive MD method in order to explicitly simulate and quantitatively describe amino acid protonation/deprotonation in proteins.
Collapse
Affiliation(s)
| | | | - Gregory A. Voth
- Department of Chemistry,
Institute for Biophysical Dynamics, James Franck Institute, and Computation
Institute, University of Chicago, 5735 S. Ellis Avenue, Chicago, Illinois 60637, United States
| | - Jessica M. J. Swanson
- Department of Chemistry,
Institute for Biophysical Dynamics, James Franck Institute, and Computation
Institute, University of Chicago, 5735 S. Ellis Avenue, Chicago, Illinois 60637, United States
| |
Collapse
|
34
|
Jose SP, Tiwary CS, Kosolwattana S, Raghavan P, Machado LD, Gautam C, Prasankumar T, Joyner J, Ozden S, Galvao DS, Ajayan PM. Enhanced supercapacitor performance of a 3D architecture tailored using atomically thin rGO–MoS2 2D sheets. RSC Adv 2016. [DOI: 10.1039/c6ra20960b] [Citation(s) in RCA: 30] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/21/2022] Open
Abstract
A stable, conductive, additive-free and scalable 3D architecture supercapacitor electrode fabricated by atomically thin 2D sheets of GO and MoS2 shows superior electrochemical properties which are further substantiated using MD simulations.
Collapse
Affiliation(s)
- Sujin P. Jose
- Department of Materials Science and Nano Engineering
- Rice University
- Houston
- USA-77005
- School of Physics
| | | | | | - Prasanth Raghavan
- Department of Materials Science and Nano Engineering
- Rice University
- Houston
- USA-77005
| | - Leonardo D. Machado
- Department of Materials Science and Nano Engineering
- Rice University
- Houston
- USA-77005
- Department of Applied Physics
| | - Chandkiram Gautam
- Department of Materials Science and Nano Engineering
- Rice University
- Houston
- USA-77005
- Department of Physics
| | - T. Prasankumar
- School of Physics
- Madurai Kamaraj University
- Madurai-625021
- India
| | - Jarin Joyner
- Department of Materials Science and Nano Engineering
- Rice University
- Houston
- USA-77005
| | - Sehmus Ozden
- Department of Materials Science and Nano Engineering
- Rice University
- Houston
- USA-77005
| | - Douglas S. Galvao
- Department of Applied Physics
- State University of Campinas
- Campinas
- Brazil
| | - P. M. Ajayan
- Department of Materials Science and Nano Engineering
- Rice University
- Houston
- USA-77005
| |
Collapse
|
35
|
Soniat M, Kumar R, Rick SW. Hydrated proton and hydroxide charge transfer at the liquid/vapor interface of water. J Chem Phys 2015; 143:044702. [DOI: 10.1063/1.4926831] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Marielle Soniat
- Department of Chemistry, University of New Orleans, New Orleans, Louisiana 70148, USA
| | - Revati Kumar
- Department of Chemistry, Louisiana State University, Baton Rouge, Louisiana 70808, USA
| | - Steven W. Rick
- Department of Chemistry, University of New Orleans, New Orleans, Louisiana 70148, USA
| |
Collapse
|
36
|
Peng Y, Swanson JMJ, Kang SG, Zhou R, Voth GA. Hydrated Excess Protons Can Create Their Own Water Wires. J Phys Chem B 2015; 119:9212-8. [PMID: 25369445 PMCID: PMC4515783 DOI: 10.1021/jp5095118] [Citation(s) in RCA: 70] [Impact Index Per Article: 7.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/19/2014] [Revised: 11/01/2014] [Indexed: 11/30/2022]
Abstract
Grotthuss shuttling of an excess proton charge defect through hydrogen bonded water networks has long been the focus of theoretical and experimental studies. In this work we show that there is a related process in which water molecules move ("shuttle") through a hydrated excess proton charge defect in order to wet the path ahead for subsequent proton charge migration. This process is illustrated through reactive molecular dynamics simulations of proton transport through a hydrophobic nanotube, which penetrates through a hydrophobic region. Surprisingly, before the proton enters the nanotube, it starts "shooting" water molecules into the otherwise dry space via Grotthuss shuttling, effectively creating its own water wire where none existed before. As the proton enters the nanotube (by 2-3 Å), it completes the solvation process, transitioning the nanotube to the fully wet state. By contrast, other monatomic cations (e.g., K(+)) have just the opposite effect, by blocking the wetting process and making the nanotube even drier. As the dry nanotube gradually becomes wet when the proton charge defect enters it, the free energy barrier of proton permeation through the tube via Grotthuss shuttling drops significantly. This finding suggests that an important wetting mechanism may influence proton translocation in biological systems, i.e., one in which protons "create" their own water structures (water "wires") in hydrophobic spaces (e.g., protein pores) before migrating through them. An existing water wire, e.g., one seen in an X-ray crystal structure or MD simulations without an explicit excess proton, is therefore not a requirement for protons to transport through hydrophobic spaces.
Collapse
Affiliation(s)
- Yuxing Peng
- †Department of Chemistry, James Franck Institute, Computation Institute, The University of Chicago, 5735 South Ellis Avenue, Chicago, Illinois 60637, United States
| | - Jessica M J Swanson
- †Department of Chemistry, James Franck Institute, Computation Institute, The University of Chicago, 5735 South Ellis Avenue, Chicago, Illinois 60637, United States
| | - Seung-gu Kang
- ‡Computational Biology Center, IBM Thomas J. Watson Research Center, Yorktown Heights, New York 10598, United States
| | - Ruhong Zhou
- ‡Computational Biology Center, IBM Thomas J. Watson Research Center, Yorktown Heights, New York 10598, United States
| | - Gregory A Voth
- †Department of Chemistry, James Franck Institute, Computation Institute, The University of Chicago, 5735 South Ellis Avenue, Chicago, Illinois 60637, United States
| |
Collapse
|
37
|
Moore SG, Crozier PS. Extension and evaluation of the multilevel summation method for fast long-range electrostatics calculations. J Chem Phys 2015; 140:234112. [PMID: 24952528 DOI: 10.1063/1.4883695] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/24/2022] Open
Abstract
Several extensions and improvements have been made to the multilevel summation method (MSM) of computing long-range electrostatic interactions. These include pressure calculation, an improved error estimator, faster direct part calculation, extension to non-orthogonal (triclinic) systems, and parallelization using the domain decomposition method. MSM also allows fully non-periodic long-range electrostatics calculations which are not possible using traditional Ewald-based methods. In spite of these significant improvements to the MSM algorithm, the particle-particle particle-mesh (PPPM) method was still found to be faster for the periodic systems we tested on a single processor. However, the fast Fourier transforms (FFTs) that PPPM relies on represent a major scaling bottleneck for the method when running on many cores (because the many-to-many communication pattern of the FFT becomes expensive) and MSM scales better than PPPM when using a large core count for two test problems on Sandia's Redsky machine. This FFT bottleneck can be reduced by running PPPM on only a subset of the total processors. MSM is most competitive for relatively low accuracy calculations. On Sandia's Chama machine, however, PPPM is found to scale better than MSM for all core counts that we tested. These results suggest that PPPM is usually more efficient than MSM for typical problems running on current high performance computers. However, further improvements to MSM algorithm could increase its competitiveness for calculation of long-range electrostatic interactions.
Collapse
Affiliation(s)
- Stan G Moore
- Sandia National Laboratories, P.O. Box 5800, MS 1322, Albuquerque, New Mexico 87185-1322, USA
| | - Paul S Crozier
- Sandia National Laboratories, P.O. Box 5800, MS 1322, Albuquerque, New Mexico 87185-1322, USA
| |
Collapse
|
38
|
Akimov AV, Prezhdo OV. Large-Scale Computations in Chemistry: A Bird’s Eye View of a Vibrant Field. Chem Rev 2015; 115:5797-890. [DOI: 10.1021/cr500524c] [Citation(s) in RCA: 159] [Impact Index Per Article: 17.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/08/2023]
Affiliation(s)
- Alexey V. Akimov
- Department
of Chemistry, University of South California, Los Angeles, California 90089, United States
| | - Oleg V. Prezhdo
- Department
of Chemistry, University of South California, Los Angeles, California 90089, United States
| |
Collapse
|
39
|
Sode O, Voth GA. Electron transfer activation of a second water channel for proton transport in [FeFe]-hydrogenase. J Chem Phys 2014; 141:22D527. [DOI: 10.1063/1.4902236] [Citation(s) in RCA: 17] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/31/2022] Open
Affiliation(s)
- Olaseni Sode
- Department of Chemistry, James Franck Institute, Institute for Biophysical Dynamics, Computation Institute, The University of Chicago, Chicago, Illinois 60637, USA and Computing, Environment and Life Sciences, Argonne National Laboratory, Argonne, Illinois 60439, USA
| | - Gregory A. Voth
- Department of Chemistry, James Franck Institute, Institute for Biophysical Dynamics, Computation Institute, The University of Chicago, Chicago, Illinois 60637, USA and Computing, Environment and Life Sciences, Argonne National Laboratory, Argonne, Illinois 60439, USA
| |
Collapse
|
40
|
Multiscale simulation reveals a multifaceted mechanism of proton permeation through the influenza A M2 proton channel. Proc Natl Acad Sci U S A 2014; 111:9396-401. [PMID: 24979779 DOI: 10.1073/pnas.1401997111] [Citation(s) in RCA: 66] [Impact Index Per Article: 6.6] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/26/2022] Open
Abstract
The influenza A virus M2 channel (AM2) is crucial in the viral life cycle. Despite many previous experimental and computational studies, the mechanism of the activating process in which proton permeation acidifies the virion to release the viral RNA and core proteins is not well understood. Herein the AM2 proton permeation process has been systematically characterized using multiscale computer simulations, including quantum, classical, and reactive molecular dynamics methods. We report, to our knowledge, the first complete free-energy profiles for proton transport through the entire AM2 transmembrane domain at various pH values, including explicit treatment of excess proton charge delocalization and shuttling through the His37 tetrad. The free-energy profiles reveal that the excess proton must overcome a large free-energy barrier to diffuse to the His37 tetrad, where it is stabilized in a deep minimum reflecting the delocalization of the excess charge among the histidines and the cost of shuttling the proton past them. At lower pH values the His37 tetrad has a larger total charge that increases the channel width, hydration, and solvent dynamics, in agreement with recent 2D-IR spectroscopic studies. The proton transport barrier becomes smaller, despite the increased charge repulsion, due to backbone expansion and the more dynamic pore water molecules. The calculated conductances are in quantitative agreement with recent experimental measurements. In addition, the free-energy profiles and conductances for proton transport in several mutants provide insights for explaining our findings and those of previous experimental mutagenesis studies.
Collapse
|
41
|
Nelson JG, Peng Y, Silverstein DW, Swanson JMJ. Multiscale Reactive Molecular Dynamics for Absolute p Ka Predictions and Amino Acid Deprotonation. J Chem Theory Comput 2014; 10:2729-2737. [PMID: 25061442 PMCID: PMC4095931 DOI: 10.1021/ct500250f] [Citation(s) in RCA: 37] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/23/2014] [Indexed: 01/16/2023]
Abstract
Accurately calculating a weak acid's pKa from simulations remains a challenging task. We report a multiscale theoretical approach to calculate the free energy profile for acid ionization, resulting in accurate absolute pKa values in addition to insights into the underlying mechanism. Importantly, our approach minimizes empiricism by mapping electronic structure data (QM/MM forces) into a reactive molecular dynamics model capable of extensive sampling. Consequently, the bulk property of interest (the absolute pKa) is the natural consequence of the model, not a parameter used to fit it. This approach is applied to create reactive models of aspartic and glutamic acids. We show that these models predict the correct pKa values and provide ample statistics to probe the molecular mechanism of dissociation. This analysis shows changes in the solvation structure and Zundel-dominated transitions between the protonated acid, contact ion pair, and bulk solvated excess proton.
Collapse
Affiliation(s)
- J Gard Nelson
- Department of Chemistry, Institute for Biophysical Dynamics, and Computation Institute, University of Chicago , 5735 S. Ellis Ave., Chicago, Illinois 60637, United States
| | - Yuxing Peng
- Department of Chemistry, Institute for Biophysical Dynamics, and Computation Institute, University of Chicago , 5735 S. Ellis Ave., Chicago, Illinois 60637, United States
| | - Daniel W Silverstein
- Department of Chemistry, Institute for Biophysical Dynamics, and Computation Institute, University of Chicago , 5735 S. Ellis Ave., Chicago, Illinois 60637, United States
| | - Jessica M J Swanson
- Department of Chemistry, Institute for Biophysical Dynamics, and Computation Institute, University of Chicago , 5735 S. Ellis Ave., Chicago, Illinois 60637, United States
| |
Collapse
|
42
|
Kumar R, Knight C, Voth GA. Exploring the behaviour of the hydrated excess proton at hydrophobic interfaces. Faraday Discuss 2013; 167:263-78. [DOI: 10.1039/c3fd00087g] [Citation(s) in RCA: 31] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/18/2022]
|