1
|
Wu Y, Cao S, Qiu Y, Huang X. Tutorial on how to build non-Markovian dynamic models from molecular dynamics simulations for studying protein conformational changes. J Chem Phys 2024; 160:121501. [PMID: 38516972 DOI: 10.1063/5.0189429] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/28/2023] [Accepted: 02/20/2024] [Indexed: 03/23/2024] Open
Abstract
Protein conformational changes play crucial roles in their biological functions. In recent years, the Markov State Model (MSM) constructed from extensive Molecular Dynamics (MD) simulations has emerged as a powerful tool for modeling complex protein conformational changes. In MSMs, dynamics are modeled as a sequence of Markovian transitions among metastable conformational states at discrete time intervals (called lag time). A major challenge for MSMs is that the lag time must be long enough to allow transitions among states to become memoryless (or Markovian). However, this lag time is constrained by the length of individual MD simulations available to track these transitions. To address this challenge, we have recently developed Generalized Master Equation (GME)-based approaches, encoding non-Markovian dynamics using a time-dependent memory kernel. In this Tutorial, we introduce the theory behind two recently developed GME-based non-Markovian dynamic models: the quasi-Markov State Model (qMSM) and the Integrative Generalized Master Equation (IGME). We subsequently outline the procedures for constructing these models and provide a step-by-step tutorial on applying qMSM and IGME to study two peptide systems: alanine dipeptide and villin headpiece. This Tutorial is available at https://github.com/xuhuihuang/GME_tutorials. The protocols detailed in this Tutorial aim to be accessible for non-experts interested in studying the biomolecular dynamics using these non-Markovian dynamic models.
Collapse
Affiliation(s)
- Yue Wu
- Department of Chemistry, Theoretical Chemistry Institute, University of Wisconsin-Madison, Madison, Wisconsin 53706, USA
| | - Siqin Cao
- Department of Chemistry, Theoretical Chemistry Institute, University of Wisconsin-Madison, Madison, Wisconsin 53706, USA
| | - Yunrui Qiu
- Department of Chemistry, Theoretical Chemistry Institute, University of Wisconsin-Madison, Madison, Wisconsin 53706, USA
| | - Xuhui Huang
- Department of Chemistry, Theoretical Chemistry Institute, University of Wisconsin-Madison, Madison, Wisconsin 53706, USA
- Data Science Institute, University of Wisconsin-Madison, Madison, Wisconsin 53706, USA
| |
Collapse
|
2
|
Deng J, Cui Q. Efficient Sampling of Cavity Hydration in Proteins with Nonequilibrium Grand Canonical Monte Carlo and Polarizable Force Fields. J Chem Theory Comput 2024; 20:1897-1911. [PMID: 38417108 DOI: 10.1021/acs.jctc.4c00013] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/01/2024]
Abstract
Prediction of the hydration levels of protein cavities and active sites is important to both mechanistic analysis and ligand design. Due to the unique microscopic environment of these buried water molecules, a polarizable model is expected to be crucial for an accurate treatment of protein internal hydration in simulations. Here we adapt a nonequilibrium candidate Monte Carlo approach for conducting grand canonical Monte Carlo simulations with the Drude polarizable force field. The GPU implementation enables the efficient sampling of internal cavity hydration levels in biomolecular systems. We also develop an enhanced sampling approach referred to as B-walking, which satisfies detailed balance and readily combines with grand canonical integration to efficiently calculate quantitative binding free energies of water to protein cavities. Applications of these developments are illustrated in a solvent box and the polar ligand binding site in trypsin. Our simulation results show that including electronic polarization leads to a modest but clear improvement in the description of water position and occupancy compared to the crystal structure. The B-walking approach enhances the range of water sampling in different chemical potential windows and thus improves the accuracy of water binding free energy calculations.
Collapse
Affiliation(s)
- Jiahua Deng
- Department of Chemistry, Boston University, 590 Commonwealth Avenue, Boston, Massachusetts 02215, United States
| | - Qiang Cui
- Department of Chemistry, Boston University, 590 Commonwealth Avenue, Boston, Massachusetts 02215, United States
- Department of Physics, Boston University, 590 Commonwealth Avenue, Boston, Massachusetts 02215, United States
- Department of Biomedical Engineering, Boston University, 44 Cummington Mall, Boston, Massachusetts 02215, United States
| |
Collapse
|
3
|
Baserga F, Storm J, Schlesinger R, Heberle J, Stripp ST. The catalytic reaction of cytochrome c oxidase probed by in situ gas titrations and FTIR difference spectroscopy. BIOCHIMICA ET BIOPHYSICA ACTA. BIOENERGETICS 2023; 1864:149000. [PMID: 37516233 DOI: 10.1016/j.bbabio.2023.149000] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/13/2022] [Revised: 06/22/2023] [Accepted: 07/24/2023] [Indexed: 07/31/2023]
Abstract
Cytochrome c oxidase (CcO) is a transmembrane heme‑copper metalloenzyme that catalyzes the reduction of O2 to H2O at the reducing end of the respiratory electron transport chain. To understand this reaction, we followed the conversion of CcO from Rhodobacter sphaeroides between several active-ready and carbon monoxide-inhibited states via attenuated total reflection Fourier-transform infrared (ATR FTIR) difference spectroscopy. Utilizing a novel gas titration setup, we prepared the mixed-valence, CO-inhibited R2CO state as well as the fully-reduced R4 and R4CO states and induced the "active ready" oxidized state OH. These experiments are performed in the dark yielding FTIR difference spectra exclusively triggered by exposure to O2, the natural substrate of CcO. Our data demonstrate that the presence of CO at heme a3 does not impair the catalytic oxidation of CcO when the cycle starts from the fully-reduced states. Interestingly, when starting from the R2CO state, the release of the CO ligand upon purging with inert gas yield a product that is indistinguishable from photolysis-induced states. The observed changes at heme a3 in the catalytic binuclear center (BNC) result from the loss of CO and are unrelated to electronic excitation upon illumination. Based on our experiments, we re-evaluate the assignment of marker bands that appear in time-resolved photolysis and perfusion-induced experiments on CcO.
Collapse
Affiliation(s)
- Federico Baserga
- Freie Universität Berlin, Experimental Molecular Biophysics, Arnimallee 14, D-14195 Berlin, Germany
| | - Julian Storm
- Freie Universität Berlin, Genetic Biophysics, Arnimallee 14, D-14195 Berlin, Germany
| | - Ramona Schlesinger
- Freie Universität Berlin, Genetic Biophysics, Arnimallee 14, D-14195 Berlin, Germany
| | - Joachim Heberle
- Freie Universität Berlin, Experimental Molecular Biophysics, Arnimallee 14, D-14195 Berlin, Germany
| | - Sven T Stripp
- Freie Universität Berlin, Experimental Molecular Biophysics, Arnimallee 14, D-14195 Berlin, Germany; Technische Universität Berlin, Division of Physical Chemistry, Strasse des 17. Juni 115, D-10623 Berlin, Germany.
| |
Collapse
|
4
|
Cao S, Qiu Y, Kalin ML, Huang X. Integrative generalized master equation: A method to study long-timescale biomolecular dynamics via the integrals of memory kernels. J Chem Phys 2023; 159:134106. [PMID: 37787134 PMCID: PMC11005468 DOI: 10.1063/5.0167287] [Citation(s) in RCA: 5] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/11/2023] [Accepted: 09/18/2023] [Indexed: 10/04/2023] Open
Abstract
The generalized master equation (GME) provides a powerful approach to study biomolecular dynamics via non-Markovian dynamic models built from molecular dynamics (MD) simulations. Previously, we have implemented the GME, namely the quasi Markov State Model (qMSM), where we explicitly calculate the memory kernel and propagate dynamics using a discretized GME. qMSM can be constructed with much shorter MD trajectories than the MSM. However, since qMSM needs to explicitly compute the time-dependent memory kernels, it is heavily affected by the numerical fluctuations of simulation data when applied to study biomolecular conformational changes. This can lead to numerical instability of predicted long-time dynamics, greatly limiting the applicability of qMSM in complicated biomolecules. We present a new method, the Integrative GME (IGME), in which we analytically solve the GME under the condition when the memory kernels have decayed to zero. Our IGME overcomes the challenges of the qMSM by using the time integrations of memory kernels, thereby avoiding the numerical instability caused by explicit computation of time-dependent memory kernels. Using our solutions of the GME, we have developed a new approach to compute long-time dynamics based on MD simulations in a numerically stable, accurate and efficient way. To demonstrate its effectiveness, we have applied the IGME in three biomolecules: the alanine dipeptide, FIP35 WW-domain, and Taq RNA polymerase. In each system, the IGME achieves significantly smaller fluctuations for both memory kernels and long-time dynamics compared to the qMSM. We anticipate that the IGME can be widely applied to investigate biomolecular conformational changes.
Collapse
Affiliation(s)
- Siqin Cao
- Department of Chemistry, Theoretical Chemistry Institute, University of Wisconsin-Madison, Madison, Wisconsin 53706, USA
| | - Yunrui Qiu
- Department of Chemistry, Theoretical Chemistry Institute, University of Wisconsin-Madison, Madison, Wisconsin 53706, USA
| | - Michael L. Kalin
- Biophysics Graduate Program, University of Wisconsin-Madison, Madison, Wisconsin 53706, USA
| | - Xuhui Huang
- Department of Chemistry, Theoretical Chemistry Institute, University of Wisconsin-Madison, Madison, Wisconsin 53706, USA
| |
Collapse
|
5
|
Qiu Y, O’Connor MS, Xue M, Liu B, Huang X. An Efficient Path Classification Algorithm Based on Variational Autoencoder to Identify Metastable Path Channels for Complex Conformational Changes. J Chem Theory Comput 2023; 19:4728-4742. [PMID: 37382437 PMCID: PMC11042546 DOI: 10.1021/acs.jctc.3c00318] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/30/2023]
Abstract
Conformational changes (i.e., dynamic transitions between pairs of conformational states) play important roles in many chemical and biological processes. Constructing the Markov state model (MSM) from extensive molecular dynamics (MD) simulations is an effective approach to dissect the mechanism of conformational changes. When combined with transition path theory (TPT), MSM can be applied to elucidate the ensemble of kinetic pathways connecting pairs of conformational states. However, the application of TPT to analyze complex conformational changes often results in a vast number of kinetic pathways with comparable fluxes. This obstacle is particularly pronounced in heterogeneous self-assembly and aggregation processes. The large number of kinetic pathways makes it challenging to comprehend the molecular mechanisms underlying conformational changes of interest. To address this challenge, we have developed a path classification algorithm named latent-space path clustering (LPC) that efficiently lumps parallel kinetic pathways into distinct metastable path channels, making them easier to comprehend. In our algorithm, MD conformations are first projected onto a low-dimensional space containing a small set of collective variables (CVs) by time-structure-based independent component analysis (tICA) with kinetic mapping. Then, MSM and TPT are constructed to obtain the ensemble of pathways, and a deep learning architecture named the variational autoencoder (VAE) is used to learn the spatial distributions of kinetic pathways in the continuous CV space. Based on the trained VAE model, the TPT-generated ensemble of kinetic pathways can be embedded into a latent space, where the classification becomes clear. We show that LPC can efficiently and accurately identify the metastable path channels in three systems: a 2D potential, the aggregation of two hydrophobic particles in water, and the folding of the Fip35 WW domain. Using the 2D potential, we further demonstrate that our LPC algorithm outperforms the previous path-lumping algorithms by making substantially fewer incorrect assignments of individual pathways to four path channels. We expect that LPC can be widely applied to identify the dominant kinetic pathways underlying complex conformational changes.
Collapse
Affiliation(s)
- Yunrui Qiu
- Department of Chemistry, Theoretical Chemistry Institute, University of Wisconsin-Madison, Madison, WI, 53706, USA
| | - Michael S. O’Connor
- Biophysics Graduate Program, University of Wisconsin-Madison, Madison, WI, 53706, USA
| | - Mingyi Xue
- Department of Chemistry, Theoretical Chemistry Institute, University of Wisconsin-Madison, Madison, WI, 53706, USA
| | - Bojun Liu
- Department of Chemistry, Theoretical Chemistry Institute, University of Wisconsin-Madison, Madison, WI, 53706, USA
| | - Xuhui Huang
- Department of Chemistry, Theoretical Chemistry Institute, University of Wisconsin-Madison, Madison, WI, 53706, USA
- Biophysics Graduate Program, University of Wisconsin-Madison, Madison, WI, 53706, USA
| |
Collapse
|
6
|
Kubař T, Elstner M, Cui Q. Hybrid Quantum Mechanical/Molecular Mechanical Methods For Studying Energy Transduction in Biomolecular Machines. Annu Rev Biophys 2023; 52:525-551. [PMID: 36791746 PMCID: PMC10810093 DOI: 10.1146/annurev-biophys-111622-091140] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/17/2023]
Abstract
Hybrid quantum mechanical/molecular mechanical (QM/MM) methods have become indispensable tools for the study of biomolecules. In this article, we briefly review the basic methodological details of QM/MM approaches and discuss their applications to various energy transduction problems in biomolecular machines, such as long-range proton transports, fast electron transfers, and mechanochemical coupling. We highlight the particular importance for these applications of balancing computational efficiency and accuracy. Using several recent examples, we illustrate the value and limitations of QM/MM methodologies for both ground and excited states, as well as strategies for calibrating them in specific applications. We conclude with brief comments on several areas that can benefit from further efforts to make QM/MM analyses more quantitative and applicable to increasingly complex biological problems.
Collapse
Affiliation(s)
- T Kubař
- Institute of Physical Chemistry, Karlsruhe Institute of Technology, Karlsruhe, Germany;
| | - M Elstner
- Institute of Physical Chemistry, Karlsruhe Institute of Technology, Karlsruhe, Germany;
- Institute of Biological Interfaces (IBG-2), Karlsruhe Institute of Technology, Karlsruhe, Germany;
| | - Q Cui
- Department of Chemistry, Boston University, Boston, Massachusetts, USA;
- Department of Physics, Boston University, Boston, Massachusetts, USA
- Department of Biomedical Engineering, Boston University, Boston, Massachusetts, USA
| |
Collapse
|
7
|
Dominic AJ, Cao S, Montoya-Castillo A, Huang X. Memory Unlocks the Future of Biomolecular Dynamics: Transformative Tools to Uncover Physical Insights Accurately and Efficiently. J Am Chem Soc 2023; 145:9916-9927. [PMID: 37104720 DOI: 10.1021/jacs.3c01095] [Citation(s) in RCA: 4] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/29/2023]
Abstract
Conformational changes underpin function and encode complex biomolecular mechanisms. Gaining atomic-level detail of how such changes occur has the potential to reveal these mechanisms and is of critical importance in identifying drug targets, facilitating rational drug design, and enabling bioengineering applications. While the past two decades have brought Markov state model techniques to the point where practitioners can regularly use them to glimpse the long-time dynamics of slow conformations in complex systems, many systems are still beyond their reach. In this Perspective, we discuss how including memory (i.e., non-Markovian effects) can reduce the computational cost to predict the long-time dynamics in these complex systems by orders of magnitude and with greater accuracy and resolution than state-of-the-art Markov state models. We illustrate how memory lies at the heart of successful and promising techniques, ranging from the Fokker-Planck and generalized Langevin equations to deep-learning recurrent neural networks and generalized master equations. We delineate how these techniques work, identify insights that they can offer in biomolecular systems, and discuss their advantages and disadvantages in practical settings. We show how generalized master equations can enable the investigation of, for example, the gate-opening process in RNA polymerase II and demonstrate how our recent advances tame the deleterious influence of statistical underconvergence of the molecular dynamics simulations used to parameterize these techniques. This represents a significant leap forward that will enable our memory-based techniques to interrogate systems that are currently beyond the reach of even the best Markov state models. We conclude by discussing some current challenges and future prospects for how exploiting memory will open the door to many exciting opportunities.
Collapse
Affiliation(s)
- Anthony J Dominic
- Department of Chemistry, University of Colorado Boulder, Boulder, Colorado 80309, USA
| | - Siqin Cao
- Department of Chemistry, Theoretical Chemistry Institute, University of Wisconsin-Madison, Madison, Wisconsin 53706, USA
| | | | - Xuhui Huang
- Department of Chemistry, Theoretical Chemistry Institute, University of Wisconsin-Madison, Madison, Wisconsin 53706, USA
| |
Collapse
|
8
|
Jeong KJ, Jeong S, Lee S, Son CY. Predictive Molecular Models for Charged Materials Systems: From Energy Materials to Biomacromolecules. ADVANCED MATERIALS (DEERFIELD BEACH, FLA.) 2023; 35:e2204272. [PMID: 36373701 DOI: 10.1002/adma.202204272] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 05/11/2022] [Revised: 07/05/2022] [Indexed: 06/16/2023]
Abstract
Electrostatic interactions play a dominant role in charged materials systems. Understanding the complex correlation between macroscopic properties with microscopic structures is of critical importance to develop rational design strategies for advanced materials. But the complexity of this challenging task is augmented by interfaces present in the charged materials systems, such as electrode-electrolyte interfaces or biological membranes. Over the last decades, predictive molecular simulations that are founded in fundamental physics and optimized for charged interfacial systems have proven their value in providing molecular understanding of physicochemical properties and functional mechanisms for diverse materials. Novel design strategies utilizing predictive models have been suggested as promising route for the rational design of materials with tailored properties. Here, an overview of recent advances in the understanding of charged interfacial systems aided by predictive molecular simulations is presented. Focusing on three types of charged interfaces found in energy materials and biomacromolecules, how the molecular models characterize ion structure, charge transport, morphology relation to the environment, and the thermodynamics/kinetics of molecular binding at the interfaces is discussed. The critical analysis brings two prominent field of energy materials and biological science under common perspective, to stimulate crossover in both research field that have been largely separated.
Collapse
Affiliation(s)
- Kyeong-Jun Jeong
- Department of Chemistry, Pohang University of Science and Technology (POSTECH), Pohang, 790-784, South Korea
| | - Seungwon Jeong
- Department of Chemistry, Pohang University of Science and Technology (POSTECH), Pohang, 790-784, South Korea
| | - Sangmin Lee
- Department of Chemistry, Pohang University of Science and Technology (POSTECH), Pohang, 790-784, South Korea
| | - Chang Yun Son
- Department of Chemistry, Pohang University of Science and Technology (POSTECH), Pohang, 790-784, South Korea
| |
Collapse
|
9
|
Gorriz RF, Imhof P. Interplay of Hydration and Protonation Dynamics in the K-Channel of Cytochrome c Oxidase. Biomolecules 2022; 12:biom12111615. [PMID: 36358964 PMCID: PMC9687966 DOI: 10.3390/biom12111615] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/25/2022] [Revised: 10/22/2022] [Accepted: 10/27/2022] [Indexed: 11/06/2022] Open
Abstract
Cytochrome c oxidase is a membrane protein of the respiratory chain that consumes protons and molecular oxygen to produce water and uses the resulting energy to pump protons across the membrane. Our molecular dynamics simulations with an excess proton located at different positions in one of the proton-conducting channels, the K-channel, show a clear dependence of the number of water molecules inside the channel on the proton position. A higher hydration level facilitates the formation of hydrogen-bonded chains along which proton transfer can occur. However, a sufficiently high hydration level for such proton transport is observed only when the excess proton is located above S365, i.e., the lower third of the channel. From the channel entrance up to this point, proton transport is via water molecules as proton carriers. These hydronium ions move with their surrounding water molecules, up to K362, filling and widening the channel. The conformation of K362 depends on its own protonation state and on the hydration level, suggesting its role to be proton transport from a hydronium ion at the height of K362 to the upper part of the channel via a conformational change. The protonation-dependent conformational dynamics of E101 at the bottom of the channel renders proton transfer via E101 unlikely. Instead, its role is rather that of an amplifier of H96’s proton affinity, suggesting H96 as the initial proton acceptor.
Collapse
Affiliation(s)
- Rene F. Gorriz
- Department of Physics, Freie Universität Berlin, Arnimallee 14, 14195 Berlin, Germany
- Computer Chemistry Center, Friedrich-Alexander Universität (FAU) Erlangen-Nürnberg, Nägelsbachstrasse 25, 91052 Erlangen, Germany
| | - Petra Imhof
- Department of Physics, Freie Universität Berlin, Arnimallee 14, 14195 Berlin, Germany
- Computer Chemistry Center, Friedrich-Alexander Universität (FAU) Erlangen-Nürnberg, Nägelsbachstrasse 25, 91052 Erlangen, Germany
- Correspondence:
| |
Collapse
|
10
|
Deng J, Cui Q. Electronic Polarization Is Essential for the Stabilization and Dynamics of Buried Ion Pairs in Staphylococcal Nuclease Mutants. J Am Chem Soc 2022; 144:4594-4610. [PMID: 35239338 PMCID: PMC9616648 DOI: 10.1021/jacs.2c00312] [Citation(s) in RCA: 7] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Buried charged residues play important roles in the modulation of protein stabilities and conformational dynamics and make crucial contributions to protein functions. Considering the generally nonpolar nature of protein interior, a key question concerns the contribution of electronic polarization to the stabilization and properties of buried charges. We answer this question by conducting free energy simulations using the latest polarizable CHARMM force field based on Drude oscillators for a series of Staphylococcal nuclease mutants that involve a buried Glu-Lys pair in different titration states and orientations. While a nonpolarizable model suggests that the ionized form of the buried Glu-Lys pair is more than 40 kcal/mol less stable than the charge-neutral form, the two titration states are comparable in stability when electronic polarization is included explicitly, a result better reconcilable with available experimental data. Analysis of free energy components suggests that additional stabilization of the ionized Glu-Lys pair has contributions from both the enhanced salt-bridge strength and stronger interaction between the ion-pair and surrounding protein residues and penetrated water. Despite the stronger direct interaction between Glu and Lys, the ion-pair exhibits considerably larger and faster structural fluctuations when polarization is included, due to compensation of interactions in the cavity. Collectively, observations from this work provide compelling evidence that electronic polarization is essential to the stability, hydration, dynamics, and therefore function of buried charges in proteins. Therefore, our study advocates for the explicit consideration of electronic polarization for mechanistic and engineering studies that implicate buried charged residues, such as enzymes and ion transporters.
Collapse
Affiliation(s)
- Jiahua Deng
- Department of Chemistry, Boston University, 590 Commonwealth Avenue, Boston, Massachusetts 02215, United States
| | - Qiang Cui
- Department of Chemistry, Boston University, 590 Commonwealth Avenue, Boston, Massachusetts 02215, United States.,Department of Physics, Boston University, 590 Commonwealth Avenue, Boston, Massachusetts 02215, United States.,Department of Biomedical Engineering, Boston University, 44 Cummington Mall, Boston, Massachusetts 02215, United States
| |
Collapse
|
11
|
Maag D, Mast T, Elstner M, Cui Q, Kubař T. O to bR transition in bacteriorhodopsin occurs through a proton hole mechanism. Proc Natl Acad Sci U S A 2021; 118:e2024803118. [PMID: 34561302 PMCID: PMC8488608 DOI: 10.1073/pnas.2024803118] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 08/04/2021] [Indexed: 12/27/2022] Open
Abstract
Extensive classical and quantum mechanical/molecular mechanical (QM/MM) molecular dynamics simulations are used to establish the structural features of the O state in bacteriorhodopsin (bR) and its conversion back to the bR ground state. The computed free energy surface is consistent with available experimental data for the kinetics and thermodynamics of the O to bR transition. The simulation results highlight the importance of the proton release group (PRG, consisting of Glu194/204) and the conserved arginine 82 in modulating the hydration level of the protein cavity. In particular, in the O state, deprotonation of the PRG and downward rotation of Arg82 lead to elevated hydration level and a continuous water network that connects the PRG to the protonated Asp85. Proton exchange through this water network is shown by ∼0.1-μs semiempirical QM/MM free energy simulations to occur through the generation and propagation of a proton hole, which is relayed by Asp212 and stabilized by Arg82. This mechanism provides an explanation for the observation that the D85S mutant of bacteriorhodopsin pumps chloride ions. The electrostatics-hydration coupling mechanism and the involvement of all titration states of water are likely applicable to many biomolecules involved in bioenergetic transduction.
Collapse
Affiliation(s)
- Denis Maag
- Institute of Physical Chemistry, Karlsruhe Institute of Technology, 76131 Karlsruhe, Germany
| | - Thilo Mast
- Institute of Physical Chemistry, Karlsruhe Institute of Technology, 76131 Karlsruhe, Germany
| | - Marcus Elstner
- Institute of Physical Chemistry, Karlsruhe Institute of Technology, 76131 Karlsruhe, Germany
- Institute for Biological Interfaces, Karlsruhe Institute of Technology, 76131 Karlsruhe, Germany
| | - Qiang Cui
- Department of Chemistry, Boston University, Boston, MA 02215
- Department of Physics, Boston University, Boston, MA 02215
- Department of Biomedical Engineering, Boston University, Boston, MA 02215
| | - Tomáš Kubař
- Institute of Physical Chemistry, Karlsruhe Institute of Technology, 76131 Karlsruhe, Germany;
| |
Collapse
|
12
|
Konovalov K, Unarta IC, Cao S, Goonetilleke EC, Huang X. Markov State Models to Study the Functional Dynamics of Proteins in the Wake of Machine Learning. JACS AU 2021; 1:1330-1341. [PMID: 34604842 PMCID: PMC8479766 DOI: 10.1021/jacsau.1c00254] [Citation(s) in RCA: 51] [Impact Index Per Article: 17.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/07/2021] [Indexed: 05/19/2023]
Abstract
Markov state models (MSMs) based on molecular dynamics (MD) simulations are routinely employed to study protein folding, however, their application to functional conformational changes of biomolecules is still limited. In the past few years, the field of computational chemistry has experienced a surge of advancements stemming from machine learning algorithms, and MSMs have not been left out. Unlike global processes, such as protein folding, the application of MSMs to functional conformational changes is challenging because they mostly consist of localized structural transitions. Therefore, it is critical to properly select a subset of structural features that can describe the slowest dynamics of these functional conformational changes. To address this challenge, we recommend several automatic feature selection methods such as Spectral-OASIS. To identify states in MSMs, the chosen features can be subject to dimensionality reduction methods such as TICA or deep learning based VAMPNets to project MD conformations onto a few collective variables for subsequent clustering. Another challenge for the application of MSMs to the study of functional conformational changes is the ability to comprehend their biophysical mechanisms, as MSMs built for these processes often require a large number of states. We recommend the recently developed quasi-MSMs (qMSMs) to address this issue. Compared to MSMs, qMSMs encode the non-Markovian dynamics via the generalized master equation and can significantly reduce the number of states. As a result, qMSMs can be built with a handful of states to facilitate the interpretation of functional conformational changes. In the wake of machine learning, we believe that the rapid advancement in the MSM methodology will lead to their wider application in studying functional conformational changes of biomolecules.
Collapse
Affiliation(s)
- Kirill
A. Konovalov
- Department
of Chemistry, State Key Laboratory of Molecular Neuroscience, The Hong Kong University of Science and Technology, Kowloon, Hong Kong
- Hong
Kong Center for Neurodegenerative Diseases, Hong Kong Science Park, Hong Kong
| | - Ilona Christy Unarta
- Department
of Chemical and Biological Engineering, The Hong Kong University of Science and Technology, Kowloon, Hong Kong
- Hong
Kong Center for Neurodegenerative Diseases, Hong Kong Science Park, Hong Kong
| | - Siqin Cao
- Department
of Chemistry, State Key Laboratory of Molecular Neuroscience, The Hong Kong University of Science and Technology, Kowloon, Hong Kong
- Hong
Kong Center for Neurodegenerative Diseases, Hong Kong Science Park, Hong Kong
| | - Eshani C. Goonetilleke
- Department
of Chemistry, State Key Laboratory of Molecular Neuroscience, The Hong Kong University of Science and Technology, Kowloon, Hong Kong
- Hong
Kong Center for Neurodegenerative Diseases, Hong Kong Science Park, Hong Kong
| | - Xuhui Huang
- Department
of Chemistry, State Key Laboratory of Molecular Neuroscience, The Hong Kong University of Science and Technology, Kowloon, Hong Kong
- Department
of Chemical and Biological Engineering, The Hong Kong University of Science and Technology, Kowloon, Hong Kong
- Hong
Kong Center for Neurodegenerative Diseases, Hong Kong Science Park, Hong Kong
| |
Collapse
|
13
|
Dragelj J, Mroginski MA, Knapp EW. Beating Heart of Cytochrome c Oxidase: The Shared Proton of Heme a3 Propionates. J Phys Chem B 2021; 125:9668-9677. [PMID: 34427096 DOI: 10.1021/acs.jpcb.1c03619] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/02/2023]
Abstract
Cytochrome c oxidase (CcO) pumps protons from the N-side to the P-side and consumes electrons from the P-side of the mitochondrial membrane driven by energy gained from reduction of dioxygen to water. ATP synthesis uses the resulting proton gradient and electrostatic potential difference. Since the distance a proton travels through CcO is too large for a one-step transfer process, proton-loading sites (PLS) that can carry protons transiently are necessary. One specific pump-active PLS couples to the redox reaction, thus energizing the proton to move across the membrane against electric potential and proton gradient. The PLS should also prevent proton backflow. Therefore, the propionates of the two redox-active hemes in CcO were suggested as PLS candidates although, according to CcO crystal structures, none of the four propionates can be protonated on account of strong H-bonds. Here, we show that modeling the local structure around heme a3 propionates enhances significantly their capability of carrying a proton jointly. This was not possible for the propionates of heme a. The modeled structures are stable in molecular dynamics simulations (MDS) and are energetically similar to the crystal structure. Precise electrostatic energy computations of MDS data are used to estimate the pKA values of all titratable residues in CcO. For the modeled structures, the heme a3 propionates have pKA values high enough to host a proton transiently but not too high to fix the proton permanently. The change in pKA throughout the redox reaction is sufficient to push the proton to the P-side of the membrane and to provide the protons with the necessary amount of energy for ATP synthesis.
Collapse
Affiliation(s)
- Jovan Dragelj
- Freie Universität Berlin, Institute for Chemistry and Biochemistry, Fabeckstrasse 36a, 14195 Berlin, Germany.,Department of Chemistry, Technische Universität Berlin, Strasse des 17. Juni 135, 10623 Berlin, Germany
| | - Maria Andrea Mroginski
- Department of Chemistry, Technische Universität Berlin, Strasse des 17. Juni 135, 10623 Berlin, Germany
| | - Ernst Walter Knapp
- Freie Universität Berlin, Institute for Chemistry and Biochemistry, Fabeckstrasse 36a, 14195 Berlin, Germany
| |
Collapse
|
14
|
Deng J, Cui Q. Reverse Protonation of Buried Ion-Pairs in Staphylococcal Nuclease Mutants. J Chem Theory Comput 2021; 17:4550-4563. [PMID: 34143626 DOI: 10.1021/acs.jctc.1c00355] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/24/2023]
Abstract
Although buried titratable residues in protein cavities are often of major functional importance, it is generally challenging to understand their properties such as the ionization state and factors of stabilization based on experimental studies alone. A specific set of examples involve buried Glu-Lys pairs in a series of variants of Staphylococcal nuclease, for which recent structural and thermodynamic studies appeared to suggest that both the stability and the ionization state of the buried Glu-Lys pair are sensitive to its orientation (i.e., Glu23-Lys36 vs Lys23-Glu36). To further clarify the situation, especially ionization states of the buried Glu-Lys pairs, we have conducted extensive molecular dynamics simulations and free energy computations. Microsecond molecular dynamics simulations show that the hydration level of the cavity depends on the orientation of the buried ion-pair therein as well as its ionization state; free energy simulations recapitulate the relative stability of Glu23-Lys36 (EK) vs Lys23-Glu36 (KE) mutants measured experimentally, although the difference is similar in magnitude regardless of the ionization state of the Glu-Lys pair. A complementary set of free energy simulations strongly suggests that, in contrast to the original suggestion in the experimental analysis, the Glu and Lys residues prefer to adopt their charge-neutral rather than the ionized states. This result is consistent with the low dielectric constant computed for water in the cavity, which makes it difficult for the protein cavity to stabilize a pair of charged Glu-Lys residues, even with water penetration. The current study highlights the role of free energy simulations in understanding the ionization state of buried titratable residues and the relevant energetic contributions, forming the basis for the rational design of buried charge networks in proteins.
Collapse
Affiliation(s)
- Jiahua Deng
- Department of Chemistry, Boston University, 590 Commonwealth Avenue, Boston, Massachusetts 02215, United States
| | - Qiang Cui
- Departments of Chemistry, Physics, and Biomedical Engineering, Boston University, 590 Commonwealth Avenue, Boston, Massachusetts 02215, United States
| |
Collapse
|
15
|
Kaur D, Khaniya U, Zhang Y, Gunner MR. Protein Motifs for Proton Transfers That Build the Transmembrane Proton Gradient. Front Chem 2021; 9:660954. [PMID: 34211960 PMCID: PMC8239185 DOI: 10.3389/fchem.2021.660954] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/30/2021] [Accepted: 05/31/2021] [Indexed: 11/13/2022] Open
Abstract
Biological membranes are barriers to polar molecules, so membrane embedded proteins control the transfers between cellular compartments. Protein controlled transport moves substrates and activates cellular signaling cascades. In addition, the electrochemical gradient across mitochondrial, bacterial and chloroplast membranes, is a key source of stored cellular energy. This is generated by electron, proton and ion transfers through proteins. The gradient is used to fuel ATP synthesis and to drive active transport. Here the mechanisms by which protons move into the buried active sites of Photosystem II (PSII), bacterial RCs (bRCs) and through the proton pumps, Bacteriorhodopsin (bR), Complex I and Cytochrome c oxidase (CcO), are reviewed. These proteins all use water filled proton transfer paths. The proton pumps, that move protons uphill from low to high concentration compartments, also utilize Proton Loading Sites (PLS), that transiently load and unload protons and gates, which block backflow of protons. PLS and gates should be synchronized so PLS proton affinity is high when the gate opens to the side with few protons and low when the path is open to the high concentration side. Proton transfer paths in the proteins we describe have different design features. Linear paths are seen with a unique entry and exit and a relatively straight path between them. Alternatively, paths can be complex with a tangle of possible routes. Likewise, PLS can be a single residue that changes protonation state or a cluster of residues with multiple charge and tautomer states.
Collapse
Affiliation(s)
- Divya Kaur
- Department of Chemistry, The Graduate Center, City University of New York, New York, NY, United States.,Department of Physics, City College of New York, New York, NY, United States
| | - Umesh Khaniya
- Department of Physics, City College of New York, New York, NY, United States.,Department of Physics, The Graduate Center, City University of New York, New York, NY, United States
| | - Yingying Zhang
- Department of Physics, City College of New York, New York, NY, United States.,Department of Physics, The Graduate Center, City University of New York, New York, NY, United States
| | - M R Gunner
- Department of Chemistry, The Graduate Center, City University of New York, New York, NY, United States.,Department of Physics, City College of New York, New York, NY, United States.,Department of Physics, The Graduate Center, City University of New York, New York, NY, United States
| |
Collapse
|
16
|
Baserga F, Dragelj J, Kozuch J, Mohrmann H, Knapp EW, Stripp ST, Heberle J. Quantification of Local Electric Field Changes at the Active Site of Cytochrome c Oxidase by Fourier Transform Infrared Spectroelectrochemical Titrations. Front Chem 2021; 9:669452. [PMID: 33987170 PMCID: PMC8111224 DOI: 10.3389/fchem.2021.669452] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/18/2021] [Accepted: 03/24/2021] [Indexed: 11/30/2022] Open
Abstract
Cytochrome c oxidase (CcO) is a transmembrane protein complex that reduces molecular oxygen to water while translocating protons across the mitochondrial membrane. Changes in the redox states of its cofactors trigger both O2 reduction and vectorial proton transfer, which includes a proton-loading site, yet unidentified. In this work, we exploited carbon monoxide (CO) as a vibrational Stark effect (VSE) probe at the binuclear center of CcO from Rhodobacter sphaeroides. The CO stretching frequency was monitored as a function of the electrical potential, using Fourier transform infrared (FTIR) absorption spectroelectrochemistry. We observed three different redox states (R4CO, R2CO, and O), determined their midpoint potential, and compared the resulting electric field to electrostatic calculations. A change in the local electric field strength of +2.9 MV/cm was derived, which was induced by the redox transition from R4CO to R2CO. We performed potential jump experiments to accumulate the R2CO and R4CO species and studied the FTIR difference spectra in the protein fingerprint region. The comparison of the experimental and computational results reveals that the key glutamic acid residue E286 is protonated in the observed states, and that its hydrogen-bonding environment is disturbed upon the redox transition of heme a3. Our experiments also suggest propionate A of heme a3 changing its protonation state in concert with the redox state of a second cofactor, heme a. This supports the role of propionic acid side chains as part of the proton-loading site.
Collapse
Affiliation(s)
- Federico Baserga
- Department of Physics, Experimental Molecular Biophysics, Freie Universität Berlin, Berlin, Germany
| | - Jovan Dragelj
- Macromolecular Modelling Group, Institute of Chemistry and Biochemistry, Freie Universität Berlin, Berlin, Germany.,Modeling of Biomolecular Systems, Technische Universität Berlin, Berlin, Germany
| | - Jacek Kozuch
- Department of Physics, Experimental Molecular Biophysics, Freie Universität Berlin, Berlin, Germany
| | - Hendrik Mohrmann
- Department of Physics, Experimental Molecular Biophysics, Freie Universität Berlin, Berlin, Germany
| | - Ernst-Walter Knapp
- Macromolecular Modelling Group, Institute of Chemistry and Biochemistry, Freie Universität Berlin, Berlin, Germany
| | - Sven T Stripp
- Department of Physics, Experimental Molecular Biophysics, Freie Universität Berlin, Berlin, Germany
| | - Joachim Heberle
- Department of Physics, Experimental Molecular Biophysics, Freie Universität Berlin, Berlin, Germany
| |
Collapse
|
17
|
Abstract
QM/MM simulations have become an indispensable tool in many chemical and biochemical investigations. Considering the tremendous degree of success, including recognition by a 2013 Nobel Prize in Chemistry, are there still "burning challenges" in QM/MM methods, especially for biomolecular systems? In this short Perspective, we discuss several issues that we believe greatly impact the robustness and quantitative applicability of QM/MM simulations to many, if not all, biomolecules. We highlight these issues with observations and relevant advances from recent studies in our group and others in the field. Despite such limited scope, we hope the discussions are of general interest and will stimulate additional developments that help push the field forward in meaningful directions.
Collapse
Affiliation(s)
- Qiang Cui
- Departments of Chemistry, Physics, and Biomedical Engineering, Boston University, 590 Commonwealth Avenue, Boston, Massachusetts 02215, United States
| | - Tanmoy Pal
- Department of Chemistry, Boston University, 590 Commonwealth Avenue, Boston, Massachusetts 02215, United States
| | - Luke Xie
- Department of Chemistry, Boston University, 590 Commonwealth Avenue, Boston, Massachusetts 02215, United States
| |
Collapse
|
18
|
Cai X, Son CY, Mao J, Kaur D, Zhang Y, Khaniya U, Cui Q, Gunner MR. Identifying the proton loading site cluster in the ba 3 cytochrome c oxidase that loads and traps protons. BIOCHIMICA ET BIOPHYSICA ACTA-BIOENERGETICS 2020; 1861:148239. [PMID: 32531221 DOI: 10.1016/j.bbabio.2020.148239] [Citation(s) in RCA: 12] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/09/2020] [Revised: 05/05/2020] [Accepted: 06/04/2020] [Indexed: 12/29/2022]
Abstract
Cytochrome c Oxidase (CcO) is the terminal electron acceptor in aerobic respiratory chain, reducing O2 to water. The released free energy is stored by pumping protons through the protein, maintaining the transmembrane electrochemical gradient. Protons are held transiently in a proton loading site (PLS) that binds and releases protons driven by the electron transfer reaction cycle. Multi-Conformation Continuum Electrostatics (MCCE) was applied to crystal structures and Molecular Dynamics snapshots of the B-type Thermus thermophilus CcO. Six residues are identified as the PLS, binding and releasing protons as the charges on heme b and the binuclear center are changed: the heme a3 propionic acids, Asp287, Asp372, His376 and Glu126B. The unloaded state has one proton and the loaded state two protons on these six residues. Different input structures, modifying the PLS conformation, show different proton distributions and result in different proton pumping behaviors. One loaded and one unloaded protonation states have the loaded/unloaded states close in energy so the PLS binds and releases a proton through the reaction cycle. The alternative proton distributions have state energies too far apart to be shifted by the electron transfers so are locked in loaded or unloaded states. Here the protein can use active states to load and unload protons, but has nearby trapped states, which stabilize PLS protonation state, providing new ideas about the CcO proton pumping mechanism. The distance between the PLS residues Asp287 and His376 correlates with the energy difference between loaded and unloaded states.
Collapse
Affiliation(s)
- Xiuhong Cai
- Department of Physics, City College of New York, 160 Convent Avenue, New York, NY 10031, USA; Department of Physics, Graduate Center, City University of New York, 365 Fifth Avenue, New York, NY 10016, USA
| | - Chang Yun Son
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, CA 91125, USA; Department of Chemistry and Division of Advanced Materials Science, Pohang University of Science and Technology (POSTECH), Pohang 37673, South Korea
| | - Junjun Mao
- Department of Physics, City College of New York, 160 Convent Avenue, New York, NY 10031, USA
| | - Divya Kaur
- Department of Physics, City College of New York, 160 Convent Avenue, New York, NY 10031, USA; Department of Chemistry, Graduate Center, City University of New York, 365 Fifth Avenue, New York, NY 10016, USA
| | - Yingying Zhang
- Department of Physics, City College of New York, 160 Convent Avenue, New York, NY 10031, USA; Department of Physics, Graduate Center, City University of New York, 365 Fifth Avenue, New York, NY 10016, USA
| | - Umesh Khaniya
- Department of Physics, City College of New York, 160 Convent Avenue, New York, NY 10031, USA; Department of Physics, Graduate Center, City University of New York, 365 Fifth Avenue, New York, NY 10016, USA
| | - Qiang Cui
- Department of Chemistry & Department of Biomedical Engineering & Department of Physics, Boston University, 590 Commonwealth Avenue, Boston, MA 02215, USA
| | - M R Gunner
- Department of Physics, City College of New York, 160 Convent Avenue, New York, NY 10031, USA; Department of Physics, Graduate Center, City University of New York, 365 Fifth Avenue, New York, NY 10016, USA; Department of Chemistry, Graduate Center, City University of New York, 365 Fifth Avenue, New York, NY 10016, USA.
| |
Collapse
|
19
|
Palese LL. Oxygen-oxygen distances in protein-bound crystallographic water suggest the presence of protonated clusters. Biochim Biophys Acta Gen Subj 2020; 1864:129480. [DOI: 10.1016/j.bbagen.2019.129480] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/17/2019] [Revised: 10/27/2019] [Accepted: 10/28/2019] [Indexed: 12/11/2022]
|
20
|
Wolf A, Dragelj J, Wonneberg J, Stellmacher J, Balke J, Woelke AL, Hodoscek M, Knapp EW, Alexiev U. The redox-coupled proton-channel opening in cytochrome c oxidase. Chem Sci 2020. [DOI: 10.1039/c9sc06463j] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022] Open
Abstract
The interplay of cytochrome c oxidase's cofactor electrostatics, long-range conformational changes, H-bond rearrangement, and water dynamics enables transient proton-channel activation.
Collapse
Affiliation(s)
- Alexander Wolf
- Physics Department
- Freie Universität Berlin
- 14195 Berlin
- Germany
| | - Jovan Dragelj
- Institute of Chemistry and Biochemistry
- Freie Universität Berlin
- 14195 Berlin
- Germany
| | | | | | - Jens Balke
- Physics Department
- Freie Universität Berlin
- 14195 Berlin
- Germany
| | - Anna Lena Woelke
- Institute of Chemistry and Biochemistry
- Freie Universität Berlin
- 14195 Berlin
- Germany
- Department of Chemistry
| | - Milan Hodoscek
- Institute of Chemistry and Biochemistry
- Freie Universität Berlin
- 14195 Berlin
- Germany
- National Institute of Chemistry
| | - Ernst Walter Knapp
- Institute of Chemistry and Biochemistry
- Freie Universität Berlin
- 14195 Berlin
- Germany
| | - Ulrike Alexiev
- Physics Department
- Freie Universität Berlin
- 14195 Berlin
- Germany
| |
Collapse
|
21
|
Wolf A, Wonneberg J, Balke J, Alexiev U. Electronation-dependent structural change at the proton exit side of cytochrome c oxidase as revealed by site-directed fluorescence labeling. FEBS J 2019; 287:1232-1246. [PMID: 31597007 DOI: 10.1111/febs.15084] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/05/2019] [Revised: 08/17/2019] [Accepted: 10/04/2019] [Indexed: 02/05/2023]
Abstract
Cytochrome c oxidase (CcO), the terminal enzyme of the respiratory chain of mitochondria and many aerobic prokaryotes that function as a redox-coupled proton pump, catalyzes the reduction of molecular oxygen to water. As part of the respiratory chain, CcO contributes to the proton motive force driving ATP synthesis. While many aspects of the enzyme's catalytic mechanisms have been established, a clear picture of the proton exit pathway(s) remains elusive. Here, we aim to gain insight into the molecular mechanisms of CcO through the development of a new homologous mutagenesis/expression system in Paracoccus denitrificans, which allows mutagenesis of CcO subunits 1, 2, and 3. Our system provides true single thiol-reactive CcO variants in a three-subunit base variant with unique labeling sites for the covalent attachment of reporter groups sensitive to nanoenvironmental factors like protonation, polarity, and hydration. To this end, we exchanged six residues on both membrane sides of CcO for cysteines. We show redox-dependent wetting changes at the proton uptake channel and increased polarity at the proton exit side of CcO upon electronation. We suggest an electronation-dependent conformational change to play a role in proton exit from CcO.
Collapse
Affiliation(s)
- Alexander Wolf
- Institute of Experimental Physics, Freie Universität Berlin, Germany
| | - Juliane Wonneberg
- Institute of Experimental Physics, Freie Universität Berlin, Germany
| | - Jens Balke
- Institute of Experimental Physics, Freie Universität Berlin, Germany
| | - Ulrike Alexiev
- Institute of Experimental Physics, Freie Universität Berlin, Germany
| |
Collapse
|
22
|
Abstract
Photosystem II (PSII) uses water as the terminal electron donor, producing oxygen in the Mn4CaO5 oxygen evolving complex (OEC), while cytochrome c oxidase (CcO) reduces O2 to water in its heme–Cu binuclear center (BNC). Each protein is oriented in the membrane to add to the proton gradient. The OEC, which releases protons, is located near the P-side (positive, at low-pH) of the membrane. In contrast, the BNC is in the middle of CcO, so the protons needed for O2 reduction must be transferred from the N-side (negative, at high pH). In addition, CcO pumps protons from N- to P-side, coupled to the O2 reduction chemistry, to store additional energy. Thus, proton transfers are directly coupled to the OEC and BNC redox chemistry, as well as needed for CcO proton pumping. The simulations that study the changes in proton affinity of the redox active sites and the surrounding protein at different states of the reaction cycle, as well as the changes in hydration that modulate proton transfer paths, are described.
Collapse
|
23
|
Malkamäki A, Sharma V. Atomistic insights into cardiolipin binding sites of cytochrome c oxidase. BIOCHIMICA ET BIOPHYSICA ACTA-BIOENERGETICS 2018; 1860:224-232. [PMID: 30414931 DOI: 10.1016/j.bbabio.2018.11.004] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 08/06/2018] [Revised: 09/21/2018] [Accepted: 11/07/2018] [Indexed: 11/19/2022]
Abstract
Mitochondrial cytochrome c oxidase couples the reduction of oxygen to proton pumping. Despite an overall good understanding of its molecular mechanism, the role of cardiolipin in protein function is not understood. Here, we have studied the cardiolipin-protein interactions in a dynamic context by means of atomistic molecular dynamics simulations performed on the entire structure of monomeric and dimeric forms of the enzyme. Several microseconds of simulation data reveal that the crystallographic cardiolipin molecules that glue two monomers together bind weakly in hybrid and single-component lipid bilayers and dissociate rapidly. Atomistic simulations performed in the absence of tightly bound cardiolipin molecules strongly perturb the structural integrity of subunits III and VIIa, thereby highlighting an indispensable nature of lipid-protein interactions in enzyme function such as proton uptake and oxygen channeling. Our results demonstrate the strength of molecular simulations in providing direct atomic description of lipid-protein processes that are difficult to achieve experimentally.
Collapse
Affiliation(s)
- Aapo Malkamäki
- Department of Physics, P. O. Box 64, University of Helsinki, Helsinki 00014, Finland
| | - Vivek Sharma
- Department of Physics, P. O. Box 64, University of Helsinki, Helsinki 00014, Finland; Institute of Biotechnology, P. O. Box 56, University of Helsinki, Helsinki 00014, Finland.
| |
Collapse
|
24
|
Cai X, Haider K, Lu J, Radic S, Son CY, Cui Q, Gunner M. Network analysis of a proposed exit pathway for protons to the P-side of cytochrome c oxidase. BIOCHIMICA ET BIOPHYSICA ACTA-BIOENERGETICS 2018; 1859:997-1005. [DOI: 10.1016/j.bbabio.2018.05.010] [Citation(s) in RCA: 25] [Impact Index Per Article: 4.2] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Received: 03/13/2018] [Revised: 05/11/2018] [Accepted: 05/16/2018] [Indexed: 11/25/2022]
|
25
|
Wikström M, Sharma V. Proton pumping by cytochrome c oxidase – A 40 year anniversary. BIOCHIMICA ET BIOPHYSICA ACTA-BIOENERGETICS 2018; 1859:692-698. [DOI: 10.1016/j.bbabio.2018.03.009] [Citation(s) in RCA: 33] [Impact Index Per Article: 5.5] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Received: 03/08/2018] [Accepted: 03/14/2018] [Indexed: 12/21/2022]
|
26
|
Farahvash A, Stuchebrukhov A. Investigating the Many Roles of Internal Water in Cytochrome c Oxidase. J Phys Chem B 2018; 122:7625-7635. [PMID: 30011995 DOI: 10.1021/acs.jpcb.7b11920] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022]
Abstract
Cytochrome c oxidase (C cO) is the terminal enzyme in the respiratory electron transport chain. As part of its catalytic cycle, C cO transfers protons to its Fe-Cu binuclear center (BNC) to reduce oxygen, and in addition, it pumps protons across the mitochondrial inner, or bacterial, membrane where it is located. It is believed that this proton transport is facilitated by a network of water chains inside the enzyme. Here we present an analysis of the hydration of C cO, including the BNC region, using a semi-empirical hydration program, Dowser++, recently developed in our group. Using high-resolution X-ray data, we show that Dowser++ predictions match very accurately the water molecules seen in the D- and K-channels of C cO, as well as in the vicinity of its BNC. Moreover, Dowser++ predicts many more internal water molecules than is typically seen in the experiment. However, no significant hydration of the catalytic cavity in C cO described recently in the literature is observed. As Dowser++ itself does not account for structural changes of the protein, this result supports the earlier assessment that the proposed wetting transition in the catalytic cavity can only either be due to structural rearrangements of BNC, possibly induced by the charges during the catalytic cycle, or occur transiently, in concert with the proton transfer. Molecular dynamics simulations were performed to investigate the global dynamic nature of Dowser++ waters in C cO, and the results suggest a consistent explanation as to why some predicted water molecules would be missing in the experimental structures. Furthermore, in light of the significant protein hydration predicted by Dowser++, the dielectric constant of the hydrated cavities in C cO was also investigated using the Fröhlich-Kirkwood model; the results indicate that in the cavities where water is packed sufficiently densely the dielectric constant can approach values comparable even to that of bulk water.
Collapse
Affiliation(s)
- Ardavan Farahvash
- Department of Chemistry , University of California-Davis , One Shields Avenue , Davis , California 95616 , United States
| | - Alexei Stuchebrukhov
- Department of Chemistry , University of California-Davis , One Shields Avenue , Davis , California 95616 , United States
| |
Collapse
|
27
|
Supekar S, Kaila VRI. Dewetting transitions coupled to K-channel activation in cytochrome c oxidase. Chem Sci 2018; 9:6703-6710. [PMID: 30310604 PMCID: PMC6115622 DOI: 10.1039/c8sc01587b] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/07/2018] [Accepted: 07/08/2018] [Indexed: 12/20/2022] Open
Abstract
Cytochrome c oxidase (CcO) drives aerobic respiratory chains in all organisms by transducing the free energy from oxygen reduction into an electrochemical proton gradient across a biological membrane.
Cytochrome c oxidase (CcO) drives aerobic respiratory chains in all organisms by transducing the free energy from oxygen reduction into an electrochemical proton gradient across a biological membrane. CcO employs the so-called D- and K-channels for proton uptake, but the molecular mechanism for activation of the K-channel has remained elusive for decades. We show here by combining large-scale atomistic molecular simulations with graph-theoretical water network analysis, and hybrid quantum/classical (QM/MM) free energy calculations, that the K-channel is activated by formation of a reactive oxidized intermediate in the binuclear heme a3/CuB active site. This state induces electrostatic, hydration, and conformational changes that lower the barrier for proton transfer along the K-channel by dewetting pathways that connect the D-channel with the active site. Our combined results reconcile previous experimental findings and indicate that water dynamics plays a decisive role in the proton pumping machinery in CcO.
Collapse
Affiliation(s)
- Shreyas Supekar
- Department Chemie , Technische Universität München , Lichtenbergstraße 4 , D-85748 Garching , Germany .
| | - Ville R I Kaila
- Department Chemie , Technische Universität München , Lichtenbergstraße 4 , D-85748 Garching , Germany .
| |
Collapse
|
28
|
Robinson AC, Schlessman JL, García-Moreno E B. Dielectric Properties of a Protein Probed by Reversal of a Buried Ion Pair. J Phys Chem B 2018; 122:2516-2524. [PMID: 29466010 DOI: 10.1021/acs.jpcb.7b12121] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/19/2023]
Abstract
Thirty years ago, Hwang and Warshel suggested that a microenvironment preorganized to stabilize an ion pair would be incapable of reorganizing to stabilize the reverse ion pair. The implications were that (1) proteins have a limited capacity to reorganize, even under the influence of strong interactions, such as those present when ionizable groups are buried in the hydrophobic interior of a protein, and (2) the inability of proteins to tolerate the reversal of buried ion pairs demonstrates the limitations inherent to continuum electrostatic models of proteins. Previously we showed that when buried individually in the interior of staphylococcal nuclease, Glu23 and Lys36 have p Ka values near pH 7, but when buried simultaneously, they establish a strong interaction of ∼5 kcal/mol and have p Ka values shifted toward more normal values. Here, using equilibrium thermodynamic measurements, crystal structures, and NMR spectroscopy experiments, we show that although the reversed, individual substitutions-Lys23 and Glu36-also have p Ka values near 7, when buried together, they neither establish a strong interaction nor promote reorganization of their microenvironment. These experiments both confirm Warshel's original hypothesis and expand it by showing that it applies to reorganization, as demonstrated by our artificial ion pairs, as well as to preorganization as is commonly argued for motifs that stabilize naturally occurring ion pairs in polar microenvironments. These data constitute a challenging benchmark useful to test the ability of structure-based algorithms to reproduce the compensation between self-energy, Coulomb and polar interactions in hydrophobic environments of proteins.
Collapse
Affiliation(s)
- Aaron C Robinson
- Department of Biophysics , Johns Hopkins University , Baltimore , Maryland 21218 , United States
| | - Jamie L Schlessman
- Chemistry Department , U.S. Naval Academy , Annapolis , Maryland 21402 , United States
| | - Bertrand García-Moreno E
- Department of Biophysics , Johns Hopkins University , Baltimore , Maryland 21218 , United States
| |
Collapse
|
29
|
Ghane T, Gorriz RF, Wrzalek S, Volkenandt S, Dalatieh F, Reidelbach M, Imhof P. Hydrogen-Bonded Network and Water Dynamics in the D-channel of Cytochrome c Oxidase. J Membr Biol 2018; 251:299-314. [PMID: 29435610 DOI: 10.1007/s00232-018-0019-x] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/24/2017] [Accepted: 02/06/2018] [Indexed: 01/09/2023]
Abstract
Proton transfer in cytochrome c oxidase (CcO) from the cellular inside to the binuclear redox centre as well as proton pumping through the membrane takes place through proton entrance via two distinct pathways, the D- and K-channel. Both channels show a dependence of their hydration level on the protonation states of their key residues, K362 for the K-channel, and E286 or D132 for the D-channel. In the oxidative half of CcO's catalytic cycle the D-channel is the proton-conducting path. For this channel, an interplay of protonation state of the D-channel residues with the water and hydrogen-bond dynamics has been observed in molecular dynamics simulations of the CcO protein, embedded in a lipid bi-layer, modelled in different protonation states. Protonation of residue E286 at the end of the D-channel results in a hydrogen-bonded network pointing from E286 to N139, that is against proton transport, and favouring N139 conformations which correspond to a closed asparagine gate (formed by residues N121 and N139). Consequently, the hydration level is lower than with unprotonated E286. In those models, the Asn gate is predominantly open, allowing water molecules to pass and thus increase the hydration level. The hydrogen-bonded network in these states exhibits longer life times of the Asn residues with water than other models and shows the D-channel to be traversable from the entrance, D132, to exit, E286. The D-channel can thus be regarded as auto-regulated with respect to proton transport, allowing proton passage only when required, that is the proton is located at the lower part of the D-channel (D132 to Asn gate) and not at the exit (E286).
Collapse
Affiliation(s)
- Tahereh Ghane
- Institute of Theoretical Physics, Freie Universität Berlin, Arnimallee 14, 14195, Berlin, Germany
| | - Rene F Gorriz
- Institute of Theoretical Physics, Freie Universität Berlin, Arnimallee 14, 14195, Berlin, Germany
| | - Sandro Wrzalek
- Institute of Theoretical Physics, Freie Universität Berlin, Arnimallee 14, 14195, Berlin, Germany
| | - Senta Volkenandt
- Institute of Theoretical Physics, Freie Universität Berlin, Arnimallee 14, 14195, Berlin, Germany
| | - Ferand Dalatieh
- Institute of Theoretical Physics, Freie Universität Berlin, Arnimallee 14, 14195, Berlin, Germany.,R Institute GmbH, Dortustraße 48, 14467, Potsdam, Germany
| | - Marco Reidelbach
- Institute of Theoretical Physics, Freie Universität Berlin, Arnimallee 14, 14195, Berlin, Germany
| | - Petra Imhof
- Institute of Theoretical Physics, Freie Universität Berlin, Arnimallee 14, 14195, Berlin, Germany.
| |
Collapse
|