Understanding the essential proton-pumping kinetic gates and decoupling mutations in cytochrome c oxidase

QM/CG-MM: Systematic Embedding of Quantum Mechanical Systems in a Coarse-Grained Environment with Accurate Electrostatics

Quantum Mechanics/Molecular Mechanics (QM/MM) can describe chemical reactions in molecular dynamics (MD) simulations at a much lower cost than ab initio MD. Still, it is prohibitively expensive for many systems of interest because such systems usually require long simulations for sufficient statistical sampling. Additional MM degrees of freedom are often slow and numerous but secondary in interest. Coarse-graining (CG) is well-known to be able to speed up sampling through both reduction in simulation cost and the ability to accelerate the dynamics. Therefore, embedding a QM system in a CG environment can be a promising way of expediting sampling without compromising the information about the QM subsystem. Sinitskiy and Voth first proposed the theory of Quantum Mechanics/Coarse-grained Molecular Mechanics (QM/CG-MM) with a bottom-up CG mapping. Mironenko and Voth subsequently introduced the DFT-QM/CG-MM formalism to couple a Density Functional Theory (DFT) treated QM system and to an apolar environment. Here, we present a more complete theory that addresses MM environments with significant polarity by explicitly accounting for the electrostatic coupling. We demonstrate our QM/CG-MM method with a chloride-methyl chloride SN2 reaction system in acetone, which is sensitive to solvent polarity. The method accurately recapitulates the potential of mean force for the substitution reaction, and the reaction barrier from the best model agrees with the atomistic simulations within sampling error. These models also have generalizability. In two other reactive systems that they have not been trained on, the QM/CG-MM model still achieves the same level of agreement with the atomistic QM/MM models. Finally, we show that in these examples the speed-up in the sampling is proportional to the acceleration of the rotational dynamics of the solvent in the CG system.

Molecular Dynamics Simulation of Complex Reactivity with the Rapid Approach for Proton Transport and Other Reactions (RAPTOR) Software Package

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.

Computational Characterization of the Reactivity of Compound I in Unspecific Peroxygenases

Unspecific peroxygenases (UPOs) are emerging as promising biocatalysts for selective oxyfunctionalization of unactivated C-H bonds. However, their potential in large-scale synthesis is currently constrained by suboptimal chemical selectivity. Improving the selectivity of UPOs requires a deep understanding of the molecular basis of their catalysis. Recent molecular simulations have sought to unravel UPO's selectivity and inform their design principles. However, most of these studies focused on substrate-binding poses. Few researchers have investigated how the reactivity of CpdI, the principal oxidizing intermediate in the catalytic cycle, influences selectivity in a realistic protein environment. Moreover, the influence of protein electrostatics on the reaction kinetics of CpdI has also been largely overlooked. To bridge this gap, we used multiscale simulations to interpret the regio- and enantioselective hydroxylation of the n-heptane substrate catalyzed by Agrocybe aegerita UPO (AaeUPO). We comprehensively characterized the energetics and kinetics of the hydrogen atom-transfer (HAT) step, initiated by CpdI, and the subsequent oxygen rebound step forming the product. Notably, our approach involved both free energy and potential energy evaluations in a quantum mechanics/molecular mechanics (QM/MM) setting, mitigating the dependence of results on the choice of initial conditions. These calculations illuminate the thermodynamics and kinetics of the HAT and oxygen rebound steps. Our findings highlight that both the conformational selection and the distinct chemical reactivity of different substrate hydrogen atoms together dictate the regio- and enantio-selectivity. Building on our previous study of CpdI's formation in AaeUPO, our results indicate that the HAT step is the rate-limiting step in the overall catalytic cycle. The subsequent oxygen rebound step is swift and retains the selectivity determined by the HAT step. We also pinpointed several polar and charged amino acid residues whose electrostatic potentials considerably influence the reaction barrier of the HAT step. Notably, the Glu196 residue is pivotal for both the CpdI's formation and participation in the HAT step. Our research offers in-depth insights into the catalytic cycle of AaeUPO, which will be instrumental in the rational design of UPOs with enhanced properties.

Crystallographic cyanide-probing for cytochrome c oxidase reveals structural bases suggesting that a putative proton transfer H-pathway pumps protons

Cytochrome c oxidase (CcO) reduces O2 in the O2-reduction site by sequential four-electron donations through the low-potential metal sites (CuA and Fea). Redox-coupled X-ray crystal structural changes have been identified at five distinct sites including Asp51, Arg438, Glu198, the hydroxyfarnesyl ethyl group of heme a, and Ser382, respectively. These sites interact with the putative proton-pumping H-pathway. However, the metal sites responsible for each structural change have not been identified, since these changes were detected as structural differences between the fully reduced and fully oxidized CcOs. Thus, the roles of these structural changes in the CcO function are yet to be revealed. X-ray crystal structures of cyanide-bound CcOs under various oxidation states showed that the O2-reduction site controlled only the Ser382-including site, while the low-potential metal sites induced the other changes. This finding indicates that these low-potential site-inducible structural changes are triggered by sequential electron-extraction from the low-potential sites by the O2-reduction site and that each structural change is insensitive to the oxidation and ligand-binding states of the O2-reduction site. Because the proton/electron coupling efficiency is constant (1:1), regardless of the reaction progress in the O2-reduction site, the structural changes induced by the low-potential sites are assignable to those critically involved in the proton pumping, suggesting that the H-pathway, facilitating these low-potential site-inducible structural changes, pumps protons. Furthermore, a cyanide-bound CcO structure suggests that a hypoxia-inducible activator, Higd1a, activates the O2-reduction site without influencing the electron transfer mechanism through the low-potential sites, kinetically confirming that the low-potential sites facilitate proton pump.

Understanding the Multifaceted Mechanism of Compound I Formation in Unspecific Peroxygenases through Multiscale Simulations

Unspecific peroxygenases (UPOs) can selectively oxyfunctionalize unactivated hydrocarbons by using peroxides under mild conditions. They circumvent the oxygen dilemma faced by cytochrome P450s and exhibit greater stability than the latter. As such, they hold great potential for industrial applications. A thorough understanding of their catalysis is needed to improve their catalytic performance. However, it remains elusive how UPOs effectively convert peroxide to Compound I (CpdI), the principal oxidizing intermediate in the catalytic cycle. Previous computational studies of this process primarily focused on heme peroxidases and P450s, which have significant differences in the active site from UPOs. Additionally, the roles of peroxide unbinding in the kinetics of CpdI formation, which is essential for interpreting existing experiments, have been understudied. Moreover, there has been a lack of free energy characterizations with explicit sampling of protein and hydration dynamics, which is critical for understanding the thermodynamics of the proton transport (PT) events involved in CpdI formation. To bridge these gaps, we employed multiscale simulations to comprehensively characterize the CpdI formation in wild-type UPO from Agrocybe aegerita (AaeUPO). Extensive free energy and potential energy calculations were performed in a quantum mechanics/molecular mechanics setting. Our results indicate that substrate-binding dehydrates the active site, impeding the PT from H2O2 to a nearby catalytic base (Glu196). Furthermore, the PT is coupled with considerable hydrogen bond network rearrangements near the active site, facilitating subsequent O-O bond cleavage. Finally, large unbinding free energy barriers kinetically stabilize H2O2 at the active site. These findings reveal a delicate balance among PT, hydration dynamics, hydrogen bond rearrangement, and cosubstrate unbinding, which collectively enable efficient CpdI formation. Our simulation results are consistent with kinetic measurements and offer new insights into the CpdI formation mechanism at atomic-level details, which can potentially aid the design of next-generation biocatalysts for sustainable chemical transformations of feedstocks.

Transient water wires mediate selective proton transport in designed channel proteins

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.

Protonation-State Dependence of Hydration and Interactions in the Two Proton-Conducting Channels of Cytochrome c Oxidase

Cytochrome c Oxidase (CcO), a membrane protein of the respiratory chain, pumps protons against an electrochemical gradient by using the energy of oxygen reduction to water. The ("chemical") protons required for this reaction and those pumped are taken up via two distinct channels, named D-channel and K-channel, in a step-wise and highly regulated fashion. In the reductive phase of the catalytic cycle, both channels transport protons so that the pumped proton passes the D-channel before the "chemical" proton has crossed the K-channel. By performing molecular dynamics simulations of CcO in the O→E redox state (after the arrival of the first reducing electron) with various combinations of protonation states of the D- and K-channels, we analysed the effect of protonation on the two channels. In agreement with previous work, the amount of water observed in the D-channel was significantly higher when the terminal residue E286 was not (yet) protonated than when the proton arrived at this end of the D-channel and E286 was neutral. Since a sufficient number of water molecules in the channel is necessary for proton transport, this can be understood as E286 facilitating its own protonation. K-channel hydration shows an even higher dependence on the location of the excess proton in the K-channel. Also in agreement with previous work, the K-channel exhibits a very low hydration level that likely hinders proton transfer when the excess proton is located in the lower part of the K-channel, that is, on the N-side of S365. Once the proton has passed S365 (towards the reaction site, the bi-nuclear centre (BNC)), the amount of water in the K-channel provides hydrogen-bond connectivity that renders proton transfer up to Y288 at the BNC feasible. No significant direct effect of the protonation state of one channel on the hydration level, hydrogen-bond connectivity, or interactions between protein residues in the other channel could be observed, rendering proton conductivity in the two channels independent of each other. Regulation of the order of proton uptake and proton passage in the two channels such that the "chemical" proton leaves its channel last must, therefore, be achieved by other means of communication, such as the location of the reducing electron.

Recent progress in experimental studies on the catalytic mechanism of cytochrome c oxidase

Cytochrome c oxidase (CcO) reduces molecular oxygen (O₂) to water, coupled with a proton pump from the N-side to the P-side, by receiving four electrons sequentially from the P-side to the O₂-reduction site-including Fe_a3 and Cu_B-via the two low potential metal sites; Cu_A and Fe_a. The catalytic cycle includes six intermediates as follows, R (Fe_a3 ²⁺, Cu_B ¹⁺, Tyr244OH), A (Fe_a3 ²⁺-O₂, Cu_B ¹⁺, Tyr244OH), P_m (Fe_a3 ⁴⁺ = O^2-, Cu_B ²⁺-OH^-, Tyr244O•), F (Fe_a3 ⁴⁺ = O^2-, Cu_B ²⁺-OH^-, Tyr244OH), O (Fe_a3 ³⁺-OH^-, Cu_B ²⁺-OH^-, Tyr244OH), and E (Fe_a3 ³⁺-OH^-, Cu_B ¹⁺-H₂O, Tyr244OH). CcO has three proton conducting pathways, D, K, and H. The D and K pathways connect the N-side surface with the O₂-reduction site, while the H-pathway is located across the protein from the N-side to the P-side. The proton pump is driven by electrostatic interactions between the protons to be pumped and the net positive charges created during the O₂ reduction. Two different proton pump proposals, each including either the D-pathway or H-pathway as the proton pumping site, were proposed approximately 30 years ago and continue to be under serious debate. In our view, the progress in understanding the reaction mechanism of CcO has been critically rate-limited by the resolution of its X-ray crystallographic structure. The improvement of the resolutions of the oxidized/reduced bovine CcO up to 1.5/1.6 Å resolution in 2016 provided a breakthrough in the understanding of the reaction mechanism of CcO. In this review, experimental studies on the reaction mechanism of CcO before the appearance of the 1.5/1.6 Å resolution X-ray structures are summarized as a background description. Following the summary, we will review the recent (since 2016) experimental findings which have significantly improved our understanding of the reaction mechanism of CcO including: 1) redox coupled structural changes of bovine CcO; 2) X-ray structures of all six intermediates; 3) spectroscopic findings on the intermediate species including the Tyr244 radical in the P_m form, a peroxide-bound form between the A and Pm forms, and F_r, a one-electron reduced F-form; 4) time resolved X-ray structural changes during the photolysis of CO-bound fully reduced CcO using XFEL; 5) a simulation analysis for the Pm→Pr→F transition.

Generalized Transition State Theory Treatment of Water-Assisted Proton Transport Processes in Proteins

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.

Multiscale simulation unravels the light-regulated reversible inhibition of dihydrofolate reductase by phototrexate

Molecular photoswitches are widely used in photopharmacology, where the biomolecular functions are photo-controlled reversibly with high spatiotemporal precision. Despite the success of this field, it remains elusive how the protein environment modulates the photochemical properties of photoswitches. Understanding this fundamental question is critical for designing more effective light-regulated drugs with mitigated side effects. In our recent work, we employed first-principles non-adiabatic dynamics simulations to probe the effects of protein on the trans to cis photoisomerization of phototrexate (PTX), a photochromic analog of the anticancer therapeutic methotrexate that inhibits the target enzyme dihydrofolate reductase (DHFR). Building upon this study, in this work, we employ multiscale simulations to unravel the full photocycle underlying the light-regulated reversible inhibition of DHFR by PTX, which remains elusive until now. First-principles non-adiabatic dynamics simulations reveal that the cis to trans photoisomerization quantum yield is hindered in the protein due to backward isomerization on the ground-state following non-adiabatic transition, which arises from the favorable binding of the cis isomer with the protein. However, free energy simulations indicate that cis to trans photoisomerization significantly decreases the binding affinity of the PTX. Thus, the cis to trans photoisomerization most likely precedes the ligand unbinding from the protein. We propose the most probable photocycle of the PTX-DHFR system. Our comprehensive simulations highlight the trade-offs among the binding affinity, photoisomerization quantum yield, and the thermal stability of the ligand's different isomeric forms. As such, our work reveals new design principles of light-regulated drugs in photopharmacology.

Recent Advances in Structural Studies of Cytochrome bd and Its Potential Application as a Drug Target

Cytochrome bd is a triheme copper-free terminal oxidase in membrane respiratory chains of prokaryotes. This unique molecular machine couples electron transfer from quinol to O₂ with the generation of a proton motive force without proton pumping. Apart from energy conservation, the bd enzyme plays an additional key role in the microbial cell, being involved in the response to different environmental stressors. Cytochrome bd promotes virulence in a number of pathogenic species that makes it a suitable molecular drug target candidate. This review focuses on recent advances in understanding the structure of cytochrome bd and the development of its selective inhibitors.

Multiscale Simulation of an Influenza A M2 Channel Mutant Reveals Key Features of Its Markedly Different Proton Transport Behavior

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.

A quantitative paradigm for water-assisted proton transport through proteins and other confined spaces

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.

Accurate and Transferable Reactive Molecular Dynamics Models from Constrained Density Functional Theory

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 pK_a values in water as well as their significantly shifted pK_a 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 pK_a 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.

Protein Motifs for Proton Transfers That Build the Transmembrane Proton Gradient

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.

Influenza A M2 Inhibitor Binding Understood through Mechanisms of Excess Proton Stabilization and Channel Dynamics

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.

Multiscale Simulation Reveals Passive Proton Transport Through SERCA on the Microsecond Timescale

The sarcoplasmic reticulum Ca²⁺-ATPase (SERCA) transports two Ca²⁺ ions from the cytoplasm to the reticulum lumen at the expense of ATP hydrolysis. In addition to transporting Ca²⁺, 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 Ca²⁺. Previous studies have shown the existence of a transient water-filled pore in SERCA that connects the Ca²⁺ 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 Ca²⁺ transport. We propose that this proton transport mechanism is operational and serves as a functional conduit for passive proton transport across the sarcoplasmic reticulum.

Understanding and Tracking the Excess Proton in Ab Initio Simulations; Insights from IR Spectra

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.

An Unusual Amino Acid Substitution Within Hummingbird Cytochrome c Oxidase Alters a Key Proton-Conducting Channel

Hummingbirds in flight exhibit the highest mass-specific metabolic rate of all vertebrates. The bioenergetic requirements associated with sustained hovering flight raise the possibility of unique amino acid substitutions that would enhance aerobic metabolism. Here, we have identified a non-conservative substitution within the mitochondria-encoded cytochrome c oxidase subunit I (COI) that is fixed within hummingbirds, but not among other vertebrates. This unusual change is also rare among metazoans, but can be identified in several clades with diverse life histories. We performed atomistic molecular dynamics simulations using bovine and hummingbird COI models, thereby bypassing experimental limitations imposed by the inability to modify mtDNA in a site-specific manner. Intriguingly, our findings suggest that COI amino acid position 153 (bovine numbering convention) provides control over the hydration and activity of a key proton channel in COX. We discuss potential phenotypic outcomes linked to this alteration encoded by hummingbird mitochondrial genomes.

What Coordinate Best Describes the Affinity of the Hydrated Excess Proton for the Air-Water Interface?

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.

Interfacial solvation and slow transport of hydrated excess protons in non-ionic reverse micelles

This work employs molecular dynamics simulations to investigate the solvation and transport properties of hydrated excess protons (with a hydronium-like core structure) in non-ionic Igepal CO-520 reverse micelles of various sizes in a non-polar solvent. Multiscale Reactive Molecular Dynamics (MS-RMD) simulations were used to describe vehicular and hopping diffusion during the proton transport process. As detailed herein, an excess proton shows a marked tendency to localize in the interfacial region of micellar water pools. Slow proton transport was observed which becomes faster with increasing micellar size. Further analysis reveals that the slow diffusion of an excess proton is a combined result of slow water diffusion and the low proton hopping rate. This study also confirms that a low proton hopping rate in reverse micelles stems from the interfacial solvation of hydrated excess protons and the immobilization of interfacial water. The low water density in the interfacial region makes it difficult to form a complete hydrogen bond network near the hydrated excess proton, and therefore locks in the orientation of hydrated proton cations. The immobilization of the interfacial water also slows the relaxation of the overall hydrogen bond network.

Proton transfer in the D-channel of cytochrome c oxidase modeled by a transition network approach

BACKGROUND

Determination of proton uptake pathways in Cytochrome c Oxidase is difficult due to the complexity of the system. The transition networks approach allows sampling of proton transfer pathways without predefined reaction coordinate.

METHODS

Computation of the proton transfer pathways in a model of the D-channel of cytochrome c oxidase has been performed by a transition network approach that combines discrete, optimisation based and molecular dynamics based sampling.

RESULTS

The optimal pathway involves an opening of the so-called asparagine gate, hydration of the asparagine region, the formation of a hydrogen-bonded chain, and finally concerted proton hole transport along this chain. The optimal pathway finds the protonation of residue H26 close to the channel entrance favourable for lowering the transition energies of subsequent steps, in particular, opening of the Asn gate and formation of a hydrogen-bonded chain. Residue Y33 plays an important role in shuttling the transferred proton hole.

CONCLUSIONS

The optimal pathway found by the transition network approach shows the same important characteristics as pathways determined earlier by other methods. The computed barrier and reaction energies are also in good agreement with previous studies. The transition network approach provides an alternative to explore pathways in complex systems.

GENERAL SIGNIFICANCE

The correct function of the enzyme as oxidase and proton pump depends on the interplay of several redox and proton transport steps. Understanding the proton transport mechanism is therefore key to understanding the protein's function. The complex nature of long- distances proton transfer through a protein requires a non-trivial simulation strategy.

Proton-Induced Conformational and Hydration Dynamics in the Influenza A M2 Channel

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.

Cryo-EM structure of Neurospora crassa respiratory complex IV

In fungi, the mitochondrial respiratory chain complexes (complexes I-IV) are responsible for oxidative phosphorylation, as in higher eukaryotes. Cryo-EM was used to identify a 200 kDa membrane protein from Neurospora crassa in lipid nanodiscs as cytochrome c oxidase (complex IV) and its structure was determined at 5.5 Å resolution. The map closely resembles the cryo-EM structure of complex IV from Saccharomyces cerevisiae. Its ten subunits are conserved in S. cerevisiae and Bos taurus, but other transmembrane subunits are missing. The different structure of the Cox5a subunit is typical for fungal complex IV and may affect the interaction with complex III in a respiratory supercomplex. Additional density was found between the matrix domains of the Cox4 and Cox5a subunits that appears to be specific to N. crassa.

The yeast mitochondrial proteins Rcf1 and Rcf2 support the enzymology of the cytochrome c oxidase complex and generation of the proton motive force

The yeast mitochondrial proteins Rcf1 and Rcf2 are associated with a subpopulation of the cytochrome bc ₁-cytochrome c oxidase supercomplex and have been proposed to play a role in the assembly and/or modulation of the activity of the cytochrome c oxidase (complex IV, CIV). Yeast mutants deficient in either Rcf1 or Rcf2 proteins can use aerobic respiration-based metabolism for growth, but the absence of both proteins results in a strong growth defect. In this study, using assorted biochemical and biophysical analyses of Rcf1/Rcf2 single and double null-mutant yeast cells and mitochondria, we further explored how Rcf1 and Rcf2 support aerobic respiration and growth. We show that the absence of Rcf1 physically reduces the levels of CIV and diminishes the ability of the CIV that is present to maintain a normal mitochondrial proton motive force (PMF). Although the absence of Rcf2 did not noticeably affect the physical content of CIV, the PMF generated by CIV was also lower than normal. Our results indicate that the detrimental effects of the absence of Rcf1 and Rcf2 proteins on the CIV complex are distinct in terms of CIV assembly/accumulation and additive in terms of the ability of CIV to generate PMF. Thus, the combined absence of Rcf1 and Rcf2 alters both CIV physiology and assembly. We conclude that the slow aerobic growth of the Rcf1/Rcf2 double null mutant results from diminished generation of mitochondrial PMF by CIV and limits the level of CIV activity required for maintenance of the PMF and growth under aerobic conditions.

Prediction of perturbed proton transfer networks

The transfer of protons through proton translocating channels is a complex process, for which direct samplings of different protonation states and side chain conformations in a transition network calculation provide an efficient, bias-free description. In principle, a new transition network calculation is required for every unsampled change in the system of interest, e.g. an unsampled protonation state change, which is associated with significant computational costs. Transition networks void of or including an unsampled change are termed unperturbed or perturbed, respectively. Here, we present a prediction method, which is based on an extensive coarse-graining of the underlying transition networks to speed up the calculations. It uses the minimum spanning tree and a corresponding sensitivity analysis of an unperturbed transition network as initial guess and refinement parameter for the determination of an unknown, perturbed transition network. Thereby, the minimum spanning tree defines a sub-network connecting all nodes without cycles and minimal edge weight sum, while the sensitivity analysis analyzes the stability of the minimum spanning tree towards individual edge weight reductions. Using the prediction method, we are able to reduce the calculation costs in a model system by up to 80%, while important network properties are maintained in most predictions.

Synthetic Fe/Cu Complexes: Toward Understanding Heme-Copper Oxidase Structure and Function

Heme-copper oxidases (HCOs) are terminal enzymes on the mitochondrial or bacterial respiratory electron transport chain, which utilize a unique heterobinuclear active site to catalyze the 4H+/4e- reduction of dioxygen to water. This process involves a proton-coupled electron transfer (PCET) from a tyrosine (phenolic) residue and additional redox events coupled to transmembrane proton pumping and ATP synthesis. Given that HCOs are large, complex, membrane-bound enzymes, bioinspired synthetic model chemistry is a promising approach to better understand heme-Cu-mediated dioxygen reduction, including the details of proton and electron movements. This review encompasses important aspects of heme-O2 and copper-O2 (bio)chemistries as they relate to the design and interpretation of small molecule model systems and provides perspectives from fundamental coordination chemistry, which can be applied to the understanding of HCO activity. We focus on recent advancements from studies of heme-Cu models, evaluating experimental and computational results, which highlight important fundamental structure-function relationships. Finally, we provide an outlook for future potential contributions from synthetic inorganic chemistry and discuss their implications with relevance to biological O2-reduction.

Atomistic insights into cardiolipin binding sites of cytochrome c oxidase

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.

Investigating the Many Roles of Internal Water in Cytochrome c Oxidase

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.

Dewetting transitions coupled to K-channel activation in cytochrome c oxidase

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 a₃/Cu_B 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.

Hydrogen-Bonded Network and Water Dynamics in the D-channel of Cytochrome c Oxidase

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).

Oxygen Activation and Energy Conservation by Cytochrome c Oxidase

This review focuses on the type A cytochrome c oxidases (CcO), which are found in all mitochondria and also in several aerobic bacteria. CcO catalyzes the respiratory reduction of dioxygen (O₂) to water by an intriguing mechanism, the details of which are fairly well understood today as a result of research for over four decades. Perhaps even more intriguingly, the membrane-bound CcO couples the O₂ reduction chemistry to translocation of protons across the membrane, thus contributing to generation of the electrochemical proton gradient that is used to drive the synthesis of ATP as catalyzed by the rotary ATP synthase in the same membrane. After reviewing the structure of the core subunits of CcO, the active site, and the transfer paths of electrons, protons, oxygen, and water, we describe the states of the catalytic cycle and point out the few remaining uncertainties. Finally, we discuss the mechanism of proton translocation and the controversies in that area that still prevail.

Insights into functions of the H channel of cytochrome c oxidase from atomistic molecular dynamics simulations

Proton pumping A-type cytochrome c oxidase (CcO) terminates the respiratory chains of mitochondria and many bacteria. Three possible proton transfer pathways (D, K, and H channels) have been identified based on structural, functional, and mutational data. Whereas the D channel provides the route for all pumped protons in bacterial A-type CcOs, studies of bovine mitochondrial CcO have led to suggestions that its H channel instead provides this route. Here, we have studied H-channel function by performing atomistic molecular dynamics simulations on the entire, as well as core, structure of bovine CcO in a lipid-solvent environment. The majority of residues in the H channel do not undergo large conformational fluctuations. Its upper and middle regions have adequate hydration and H-bonding residues to form potential proton-conducting channels, and Asp51 exhibits conformational fluctuations that have been observed crystallographically. In contrast, throughout the simulations, we do not observe transient water networks that could support proton transfer from the N phase toward heme a via neutral His413, regardless of a labile H bond between Ser382 and the hydroxyethylfarnesyl group of heme a In fact, the region around His413 only became sufficiently hydrated when His413 was fixed in its protonated imidazolium state, but its calculated pK_a is too low for this to provide the means to create a proton transfer pathway. Our simulations show that the electric dipole moment of residues around heme a changes with the redox state, hence suggesting that the H channel could play a more general role as a dielectric well.