Intricate role of water in proton transport through cytochrome c oxidase

Vitamin C as Scavenger of Reactive Oxygen Species during Healing after Myocardial Infarction

Currently, coronary artery bypass and reperfusion therapies are considered the gold standard in long-term treatments to restore heart function after acute myocardial infarction. As a drawback of these restoring strategies, reperfusion after an ischemic insult and sudden oxygen exposure lead to the exacerbated synthesis of additional reactive oxidative species and the persistence of increased oxidation levels. Attempts based on antioxidant treatment have failed to achieve an effective therapy for cardiovascular disease patients. The controversial use of vitamin C as an antioxidant in clinical practice is comprehensively systematized and discussed in this review. The dose-dependent adsorption and release kinetics mechanism of vitamin C is complex; however, this review may provide a holistic perspective on its potential as a preventive supplement and/or for combined precise and targeted therapeutics in cardiovascular management therapy.

Role of water in cyclooxygenase catalysis and design of anti-inflammatory agents targeting two sites of the enzyme

While designing the anti-inflammatory agents targeting cyclooxygenase-2 (COX-2), we first identified a water loop around the heme playing critical role in the enzyme catalysis. The results of molecular dynamic studies supported by the strong hydrogen-bonding equilibria of the participating atoms, radical stabilization energies, the pK_a of the H-donor/acceptor sites and the cyclooxygenase activity of pertinent muCOX-2 ravelled the working of the water–peptide channel for coordinating the flow of H·/electron between the heme and Y385. Based on the working of H·/electron transfer channel between the 12.5 Å distant radical generation and the radical disposal sites, a series of molecules was designed and synthesized. Among this category of compounds, an appreciably potent anti-inflammatory agent exhibiting IC₅₀ 0.06 μM against COX-2 and reversing the formalin induced analgesia and carageenan induced inflammation in mice by 90% was identified. Further it was revealed that, justifying its bidentate design, the compound targets water loop (heme bound site) and the arachidonic acid binding pockets of COX-2.

Ultrafast Proton Transport between a Hydroxy Acid and a Nitrogen Base along Solvent Bridges Governed by the Hydroxide/Methoxide Transfer Mechanism

Aqueous proton transport plays a key role in acid–base neutralization and energy transport through biological membranes and hydrogen fuel cells. Extensive experimental and theoretical studies have resulted in a highly detailed elucidation of one of the underlying microscopic mechanisms for aqueous excess proton transport, known as the von Grotthuss mechanism, involving different hydrated proton configurations with associated high fluxional structural dynamics. Hydroxide transport, with approximately 2-fold-lower bulk diffusion rates compared to those of excess protons, has received much less attention. We present femtosecond UV/IR pump–probe experiments and ab initio molecular dynamics simulations of different proton transport pathways of bifunctional photoacid 7-hydroxyquinoline (7HQ) in water/methanol mixtures. For 7HQ solvent-dependent photoacidity, free-energy–reactivity correlation behavior and quantum mechanics/molecular mechanics (QM/MM) trajectories point to a dominant OH^–/CH₃O^– transport pathway for all water/methanol mixing ratios investigated. Our joint ultrafast infrared spectroscopic and ab initio molecular dynamics study provides conclusive evidence for the hydrolysis/methanolysis acid–base neutralization pathway, as formulated by Manfred Eigen half a century ago. Our findings on the distinctly different acid–base reactivities for aromatic hydroxyl and aromatic nitrogen functionalities suggest the usefulness of further exploration of these free-energy–reactivity correlations as a function of solvent polarity. Ultimately the determination of solvent-dependent acidities will contribute to a better understanding of proton-transport mechanisms at weakly polar surfaces and near polar or ionic regions in transmembrane proton pump proteins or hydrogen fuel cell materials.

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.

Multiscale modelization in a small virus: Mechanism of proton channeling and its role in triggering capsid disassembly

In this work, we assess a previously advanced hypothesis that predicts the existence of ion channels in the capsid of small and non-enveloped icosahedral viruses. With this purpose we examine Triatoma Virus (TrV) as a case study. This virus has a stable capsid under highly acidic conditions but disassembles and releases the genome in alkaline environments. Our calculations range from a subtle sub-atomic proton interchange to the dismantling of a large-scale system representing several million of atoms. Our results provide structure-based explanations for the three roles played by the capsid to enable genome release. First, we observe, for the first time, the formation of a hydrophobic gate in the cavity along the five-fold axis of the wild-type virus capsid, which can be disrupted by an ion located in the pore. Second, the channel enables protons to permeate the capsid through a unidirectional Grotthuss-like mechanism, which is the most likely process through which the capsid senses pH. Finally, assuming that the proton leak promotes a charge imbalance in the interior of the capsid, we model an internal pressure that forces shell cracking using coarse-grained simulations. Although qualitatively, this last step could represent the mechanism of capsid opening that allows RNA release. All of our calculations are in agreement with current experimental data obtained using TrV and describe a cascade of events that could explain the destabilization and disassembly of similar icosahedral viruses.

Plant and animal small non-enveloped viruses are composed of a capsid shell that encloses the genome. One of the multiple functions played by the capsid is to protect the genome against host defenses and to withstand environmental aggressions, such as dehydration. This highly specialized capsule selectively recognizes and binds to the target tissue infected by the virus. In the viral cycle, the ultimate function of the capsid is to release the genome. Observations of many viruses demonstrate that the pH of the medium can trigger genome release. Nevertheless, the mechanism underlying this process at the atomic level is poorly understood. In this work, we computationally modeled the mechanism by which the capsid senses environmental pH and the destabilization process that permits genome release. Our calculations predict that a cavity that traverses the capsid functions as a hydrophobic gate, a feature already observed in membrane ion channels. Moreover, our results predict that this cavity behaves as a proton diode because the proton transit can only occur from the capsid interior to the exterior. In turn, our calculations describe a cascade of events that could explain the destabilization and dismantling of an insect virus, but this description could also apply to many vertebrate viruses.

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

Validating the Water Flooding Approach by Comparing It to Grand Canonical Monte Carlo Simulations

The study of the function of proteins on a quantitative level requires consideration of the water molecules in and around the protein. This requirement presents a major computational challenge due to the fact that the insertion of water molecules can have a very high activation barrier and would require a long simulation time. Recently, we developed a water flooding (WF) approach which is based on a postprocessing Monte Carlo ranking of possible water configurations. This approach appears to provide a very effective way for assessing the insertion free energies and determining the most likely configurations of the internal water molecules. Although the WF approach was used effectively in modeling challenging systems that have not been addressed reliably by other microscopic approaches, it was not validated by a comparison to the more rigorous grand canonical Monte Carlo (GCMC) method. Here we validate the WF approach by comparing its performance to that of the GCMC method. It is found that the WF approach reproduces the GCMC results in well-defined test cases but does so much faster. This established the WF approach as a useful strategy for finding correct water configurations in proteins and thus to provide a powerful way for studies of the functions of proteins.

Cavity hydration dynamics in cytochrome c oxidase and functional implications

Cytochrome c oxidase (CcO) is a transmembrane protein that uses the free energy of O₂ reduction to generate the proton concentration gradient across the membrane. The regulation of competitive proton transfer pathways has been established to be essential to the vectorial transport efficiency of CcO, yet the underlying mechanism at the molecular level remains lacking. Recent studies have highlighted the potential importance of hydration-level change in an internal cavity that connects the proton entrance channel, the site of O₂ reduction, and the putative proton exit route. In this work, we use atomistic molecular dynamics simulations to investigate the energetics and timescales associated with the volume fluctuation and hydration-level change in this central cavity. Extensive unrestrained molecular dynamics simulations (accumulatively [Formula: see text]4 [Formula: see text]s) and free energy computations for different chemical states of CcO support a model in which the volume and hydration level of the cavity are regulated by the protonation state of a propionate group of heme a₃ and, to a lesser degree, the redox state of heme a and protonation state of Glu286. Markov-state model analysis of [Formula: see text]2-[Formula: see text]s trajectories suggests that hydration-level change occurs on the timescale of 100-200 ns before the proton-loading site is protonated. The computed energetic and kinetic features for the cavity wetting transition suggest that reversible hydration-level change of the cavity can indeed be a key factor that regulates the branching of proton transfer events and therefore contributes to the vectorial efficiency of proton transport.

Protonation-State-Dependent Communication in Cytochrome c Oxidase

Proton transfer in cytochrome c oxidase from the cellular inside to the binuclear redox center (BNC) can occur through two distinct pathways, the D- and K-channels. For the protein to function as both a redox enzyme and a proton pump, proton transfer into the protein toward the BNC or toward a proton loading site (and ultimately through the membrane) must be highly regulated. The P_R → F transition is the first step in a catalytic cycle that requires proton transfer from the bulk at the N-side to the BNC. Molecular dynamics simulations of the P_R → F intermediate of this transition, with 16 different combinations of protonation states of key residues in the D- and K-channel, show the impact of the K-channel on the D-channel to be protonation-state dependent. Strength as well as means of communication, correlations in positions, or communication along the hydrogen-bonded network depends on the protonation state of the K-channel residue K362. The conformational and hydrogen-bond dynamics of the D-channel residue N139 is regulated by an interplay of protonation in the D-channel and K362. N139 thus assumes a gating function by which proton passage through the D-channel toward E286 is likely facilitated for states with protonated K362 and unprotonated E286. In contrast, proton passage through the D-channel is hindered by N139's preference for a closed conformation in situations with protonated E286.

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

Cytochrome c oxidase (CcO) catalyzes the reduction of oxygen to water and uses the released free energy to pump protons against the transmembrane proton gradient. To better understand the proton-pumping mechanism of the wild-type (WT) CcO, much attention has been given to the mutation of amino acid residues along the proton translocating D-channel that impair, and sometimes decouple, proton pumping from the chemical catalysis. Although their influence has been clearly demonstrated experimentally, the underlying molecular mechanisms of these mutants remain unknown. In this work, we report multiscale reactive molecular dynamics simulations that characterize the free-energy profiles of explicit proton transport through several important D-channel mutants. Our results elucidate the mechanisms by which proton pumping is impaired, thus revealing key kinetic gating features in CcO. In the N139T and N139C mutants, proton back leakage through the D-channel is kinetically favored over proton pumping due to the loss of a kinetic gate in the N139 region. In the N139L mutant, the bulky L139 side chain inhibits timely reprotonation of E286 through the D-channel, which impairs both proton pumping and the chemical reaction. In the S200V/S201V double mutant, the proton affinity of E286 is increased, which slows down both proton pumping and the chemical catalysis. This work thus not only provides insight into the decoupling mechanisms of CcO mutants, but also explains how kinetic gating in the D-channel is imperative to achieving high proton-pumping efficiency in the WT CcO.

Insights into proton translocation in cbb3 oxidase from MD simulations

Heme-copper oxidases are membrane protein complexes that catalyse the final step of the aerobic respiration, namely the reduction of oxygen to water. The energy released during catalysis is coupled to the active translocation of protons across the membrane, which contributes to the establishment of an electrochemical gradient that is used for ATP synthesis. The distinctive C-type (or cbb₃) cytochrome c oxidases, which are mostly present in proteobacteria, exhibit a number of unique structural and functional features, including high catalytic activity at low oxygen concentrations. At the moment, the functioning mechanism of C-type oxidases, in particular the proton transfer/pumping mechanism presumably via a single proton channel, is still poorly understood. In this work we used all-atom molecular dynamics simulations and continuum electrostatics calculations to obtain atomic-level insights into the hydration and dynamics of a cbb₃ oxidase. We provide the details of the water dynamics and proton transfer pathways for both the "chemical" and "pumped" protons, and show that formation of protonic connections is strongly affected by the protonation state of key residues, namely H243, E323 and H337.

Monte Carlo Simulations of Glu-242 in Cytochrome c Oxidase

Monte Carlo (MC) simulations of conformational changes and protonation of Glu-242, a key residue that shuttles protons in cytochrome c oxidase (CcO), are reported. Previous studies suggest that this residue may play a role of the valve of the enzyme proton pump. Here we examine how sensitive the results of simulations are to the computational method used. We applied both molecular mechanic (MM) and hybrid quantum mechanic:molecular mechanic (QM:MM) methods and find that the results are qualitatively different. The results indicate that the mechanism for proton gating in CcO is still an open issue.

Computationally Efficient Multiscale Reactive Molecular Dynamics to Describe Amino Acid Deprotonation in Proteins

An important challenge in the simulation of biomolecular systems is a quantitative description of the protonation and deprotonation process of amino acid residues. Despite the seeming simplicity of adding or removing a positively charged hydrogen nucleus, simulating the actual protonation/deprotonation process is inherently difficult. It requires both the explicit treatment of the excess proton, including its charge defect delocalization and Grotthuss shuttling through inhomogeneous moieties (water and amino residues), and extensive sampling of coupled condensed phase motions. In a recent paper (J. Chem. Theory Comput. 2014, 10, 2729-2737), a multiscale approach was developed to map high-level quantum mechanics/molecular mechanics (QM/MM) data into a multiscale reactive molecular dynamics (MS-RMD) model in order to describe amino acid deprotonation in bulk water. In this article, we extend the fitting approach (called FitRMD) to create MS-RMD models for ionizable amino acids within proteins. The resulting models are shown to faithfully reproduce the free energy profiles of the reference QM/MM Hamiltonian for PT inside an example protein, the ClC-ec1 H(+)/Cl(-) antiporter. Moreover, we show that the resulting MS-RMD models are computationally efficient enough to then characterize more complex 2-dimensional free energy surfaces due to slow degrees of freedom such as water hydration of internal protein cavities that can be inherently coupled to the excess proton charge translocation. The FitRMD method is thus shown to be an effective way to map ab initio level accuracy into a much more computationally efficient reactive MD method in order to explicitly simulate and quantitatively describe amino acid protonation/deprotonation in proteins.

Spectroscopic and Crystallographic Evidence for the Role of a Water-Containing H-Bond Network in Oxidase Activity of an Engineered Myoglobin

Heme-copper oxidases (HCOs) catalyze efficient reduction of oxygen to water in biological respiration. Despite progress in studying native enzymes and their models, the roles of non-covalent interactions in promoting this activity are still not well understood. Here we report EPR spectroscopic studies of cryoreduced oxy-F33Y-CuBMb, a functional model of HCOs engineered in myoglobin (Mb). We find that cryoreduction at 77 K of the O2-bound form, trapped in the conformation of the parent oxyferrous form, displays a ferric-hydroperoxo EPR signal, in contrast to the cryoreduced oxy-wild-type (WT) Mb, which is unable to deliver a proton and shows a signal from the peroxo-ferric state. Crystallography of oxy-F33Y-CuBMb reveals an extensive H-bond network involving H2O molecules, which is absent from oxy-WTMb. This H-bonding proton-delivery network is the key structural feature that transforms the reversible oxygen-binding protein, WTMb, into F33Y-CuBMb, an oxygen-activating enzyme that reduces O2 to H2O. These results provide direct evidence of the importance of H-bond networks involving H2O in conferring enzymatic activity to a designed protein. Incorporating such extended H-bond networks in designing other metalloenzymes may allow us to confer and fine-tune their enzymatic activities.

Stay Wet, Stay Stable? How Internal Water Helps the Stability of Thermophilic Proteins

We present a systematic computational investigation of the internal hydration of a set of homologous proteins of different stability content and molecular complexities. The goal of the study is to verify whether structural water can be part of the molecular mechanisms ensuring enhanced stability in thermophilic enzymes. Our free-energy calculations show that internal hydration in the thermophilic variants is generally more favorable, and that the cumulated effect of wetting multiple sites results in a meaningful contribution to stability. Moreover, thanks to a more effective capability to retain internal water, some thermophilic proteins benefit by a systematic gain from internal wetting up to their optimal working temperature. Our work supports the idea that internal wetting can be viewed as an alternative molecular variable to be tuned for increasing protein stability.

Microscopic basis for kinetic gating in Cytochrome c oxidase: insights from QM/MM analysis

Understanding the mechanism of vectorial proton pumping in biomolecules requires establishing the microscopic basis for the regulation of both thermodynamic and kinetic features of the relevant proton transfer steps.

Understanding the mechanism of vectorial proton pumping in biomolecules requires establishing the microscopic basis for the regulation of both thermodynamic and kinetic features of the relevant proton transfer steps. For the proton pump cytochrome c oxidase, while the regulation of thermodynamic driving force for key proton transfers has been discussed in great detail, the microscopic basis for the control of proton transfer kinetics has been poorly understood. Here we carry out extensive QM/MM free energy simulations to probe the kinetics of relevant proton transfer steps and analyze the effects of local structure and hydration level. We show that protonation of the proton loading site (PLS, taken to be a propionate of heme a₃) requires a concerted process in which a key glutamic acid (Glu286H) delivers the proton to the PLS while being reprotonated by an excess proton coming from the D-channel. The concerted nature of the mechanism is a crucial feature that enables the loading of the PLS before the cavity containing Glu286 is better hydrated to lower its pK_a to experimentally measured range; the charged rather than dipolar nature of the process also ensures a tight coupling with heme a reduction, as emphasized by Siegbahn and Blomberg. In addition, we find that rotational flexibility of the PLS allows its protonation before that of the binuclear center (the site where oxygen gets reduced to water). Together with our recent study (P. Goyal, et al., Proc. Natl. Acad. Sci. U. S. A., 2013, 110, 18886–18891) that focused on the modulation of Glu286 pK_a, the current work suggests a mechanism that builds in a natural sequence for the protonation of the PLS prior to that of the binuclear center. This provides microscopic support to the kinetic constraints revealed by kinetic network analysis as essential elements that ensure an efficient vectorial proton transport in cytochrome c oxidase.

Exploring O2 diffusion in A-type cytochrome c oxidases: molecular dynamics simulations uncover two alternative channels towards the binuclear site

Cytochrome c oxidases (Ccoxs) are the terminal enzymes of the respiratory chain in mitochondria and most bacteria. These enzymes couple dioxygen (O₂) reduction to the generation of a transmembrane electrochemical proton gradient. Despite decades of research and the availability of a large amount of structural and biochemical data available for the A-type Ccox family, little is known about the channel(s) used by O₂ to travel from the solvent/membrane to the heme a₃-Cu_B binuclear center (BNC). Moreover, the identification of all possible O₂ channels as well as the atomic details of O₂ diffusion is essential for the understanding of the working mechanisms of the A-type Ccox. In this work, we determined the O₂ distribution within Ccox from Rhodobacter sphaeroides, in the fully reduced state, in order to identify and characterize all the putative O₂ channels leading towards the BNC. For that, we use an integrated strategy combining atomistic molecular dynamics (MD) simulations (with and without explicit O₂ molecules) and implicit ligand sampling (ILS) calculations. Based on the 3D free energy map for O₂ inside Ccox, three channels were identified, all starting in the membrane hydrophobic region and connecting the surface of the protein to the BNC. One of these channels corresponds to the pathway inferred from the X-ray data available, whereas the other two are alternative routes for O₂ to reach the BNC. Both alternative O₂ channels start in the membrane spanning region and terminate close to Y288_I. These channels are a combination of multiple transiently interconnected hydrophobic cavities, whose opening and closure is regulated by the thermal fluctuations of the lining residues. Furthermore, our results show that, in this Ccox, the most likely (energetically preferred) routes for O₂ to reach the BNC are the alternative channels, rather than the X-ray inferred pathway.

Cytochrome c oxidases (Ccoxs), the terminal enzymes of the respiratory electron transport chain in eukaryotes and many prokaryotes, are key enzymes in aerobic respiration. These proteins couple the reduction of molecular dioxygen to water with the creation of a transmembrane electrochemical proton gradient. Over the last decades, most of the Ccoxs research focused on the mechanisms and energetics of reduction and/or proton pumping, and little emphasis has been given to the pathways used by dioxygen to reach the binuclear center, where dioxygen reduction takes place. In particular, the existence and the characteristics of the channel(s) used by O₂ to travel from the solvent/membrane to the binuclear site are still unclear. In this work, we combine all-atom molecular dynamics simulations and implicit ligand sampling calculations in order to identify and characterize the O₂ delivery channels in the Ccox from Rhodobacter sphaeroides. Altogether, our results suggest that, in this Ccox, O₂ can diffuse via three well-defined channels that start in membrane region (where O₂ solubility is higher than in the water). One of these channels corresponds to the pathway inferred from the X-ray data available, whereas the other two are alternative routes for O₂ to reach the binuclear center.

Role of Internal Water on Protein Thermal Stability: The Case of Homologous G Domains

In this work, we address the question of whether the enhanced stability of thermophilic proteins has a direct connection with internal hydration. Our model systems are two homologous G domains of different stability: the mesophilic G domain of the elongation factor thermal unstable protein from E. coli and the hyperthermophilic G domain of the EF-1α protein from S. solfataricus. Using molecular dynamics simulation at the microsecond time scale, we show that both proteins host water molecules in internal cavities and that these molecules exchange with the external solution in the nanosecond time scale. The hydration free energy of these sites evaluated via extensive calculations is found to be favorable for both systems, with the hyperthermophilic protein offering a slightly more favorable environment to host water molecules. We estimate that, under ambient conditions, the free energy gain due to internal hydration is about 1.3 kcal/mol in favor of the hyperthermophilic variant. However, we also find that, at the high working temperature of the hyperthermophile, the cavities are rather dehydrated, meaning that under extreme conditions other molecular factors secure the stability of the protein. Interestingly, we detect a clear correlation between the hydration of internal cavities and the protein conformational landscape. The emerging picture is that internal hydration is an effective observable to probe the conformational landscape of proteins. In the specific context of our investigation, the analysis confirms that the hyperthermophilic G domain is characterized by multiple states and it has a more flexible structure than its mesophilic homologue.

Multiscale Reactive Molecular Dynamics for Absolute pKa Predictions and Amino Acid Deprotonation

Accurately calculating a weak acid's pK_a from simulations remains a challenging task. We report a multiscale theoretical approach to calculate the free energy profile for acid ionization, resulting in accurate absolute pK_a values in addition to insights into the underlying mechanism. Importantly, our approach minimizes empiricism by mapping electronic structure data (QM/MM forces) into a reactive molecular dynamics model capable of extensive sampling. Consequently, the bulk property of interest (the absolute pK_a) is the natural consequence of the model, not a parameter used to fit it. This approach is applied to create reactive models of aspartic and glutamic acids. We show that these models predict the correct pK_a values and provide ample statistics to probe the molecular mechanism of dissociation. This analysis shows changes in the solvation structure and Zundel-dominated transitions between the protonated acid, contact ion pair, and bulk solvated excess proton.

Changing hydration level in an internal cavity modulates the proton affinity of a key glutamate in cytochrome c oxidase

Cytochrome c oxidase contributes to the transmembrane proton gradient by removing two protons from the high-pH side of the membrane each time the binuclear center active site is reduced. One proton goes to the binuclear center, whereas the other is pumped to the low-pH periplasmic space. Glutamate 286 (Glu286) has been proposed to serve as a transiently deprotonated proton donor. Using unrestrained atomistic molecular dynamics simulations, we show that the size of and water distribution in the hydrophobic cavity that holds Glu286 is controlled by the protonation state of the propionic acid of heme a3, a group on the proton outlet pathway. Protonation of the propionate disrupts hydrogen bonding to two side chains, allowing a loop to swing open. Continuum electrostatics and atomistic free-energy perturbation calculations show that the resultant changes in hydration and electrostatic interactions lower the Glu proton affinity by at least 5 kcal/mol. These changes in the internal hydration level occur in the absence of major conformational transitions and serve to stabilize needed transient intermediates in proton transport. The trigger is not the protonation of the Glu of interest, but rather the protonation of a residue ∼10 Å away. Thus, unlike local water penetration to stabilize a new charge, this finding represents a specific role for water molecules in the protein interior, mediating proton transfers and facilitating ion transport.

Functions of the hydrophilic channels in protonmotive cytochrome c oxidase

The structures and functions of hydrophilic channels in electron-transferring membrane proteins are discussed. A distinction is made between proton channels that can conduct protons and dielectric channels that are non-conducting but can dielectrically polarize in response to the introduction of charge changes in buried functional centres. Functions of the K, D and H channels found in A1-type cytochrome c oxidases are reviewed in relation to these ideas. Possible control of function by dielectric channels and their evolutionary relation to proton channels is explored.

Molecular mechanisms for generating transmembrane proton gradients

Membrane proteins use the energy of light or high energy substrates to build a transmembrane proton gradient through a series of reactions leading to proton release into the lower pH compartment (P-side) and proton uptake from the higher pH compartment (N-side). This review considers how the proton affinity of the substrates, cofactors and amino acids are modified in four proteins to drive proton transfers. Bacterial reaction centers (RCs) and photosystem II (PSII) carry out redox chemistry with the species to be oxidized on the P-side while reduction occurs on the N-side of the membrane. Terminal redox cofactors are used which have pKas that are strongly dependent on their redox state, so that protons are lost on oxidation and gained on reduction. Bacteriorhodopsin is a true proton pump. Light activation triggers trans to cis isomerization of a bound retinal. Strong electrostatic interactions within clusters of amino acids are modified by the conformational changes initiated by retinal motion leading to changes in proton affinity, driving transmembrane proton transfer. Cytochrome c oxidase (CcO) catalyzes the reduction of O2 to water. The protons needed for chemistry are bound from the N-side. The reduction chemistry also drives proton pumping from N- to P-side. Overall, in CcO the uptake of 4 electrons to reduce O2 transports 8 charges across the membrane, with each reduction fully coupled to removal of two protons from the N-side, the delivery of one for chemistry and transport of the other to the P-side.

Evolution of the genetic code by incorporation of amino acids that improved or changed protein function

Fifty years have passed since the genetic code was deciphered, but how the genetic code came into being has not been satisfactorily addressed. It is now widely accepted that the earliest genetic code did not encode all 20 amino acids found in the universal genetic code as some amino acids have complex biosynthetic pathways and likely were not available from the environment. Therefore, the genetic code evolved as pathways for synthesis of new amino acids became available. One hypothesis proposes that early in the evolution of the genetic code four amino acids-valine, alanine, aspartic acid, and glycine-were coded by GNC codons (N = any base) with the remaining codons being nonsense codons. The other sixteen amino acids were subsequently added to the genetic code by changing nonsense codons into sense codons for these amino acids. Improvement in protein function is presumed to be the driving force behind the evolution of the code, but how improved function was achieved by adding amino acids has not been examined. Based on an analysis of amino acid function in proteins, an evolutionary mechanism for expansion of the genetic code is described in which individual coded amino acids were replaced by new amino acids that used nonsense codons differing by one base change from the sense codons previously used. The improved or altered protein function afforded by the changes in amino acid function provided the selective advantage underlying the expansion of the genetic code. Analysis of amino acid properties and functions explains why amino acids are found in their respective positions in the genetic code.

Capturing the energetics of water insertion in biological systems: the water flooding approach

Consistent description of the effect of internal water in proteins has been a major challenge for both simulation and experimental studies. Describing this effect has been particularly important and elusive in cases of charges in protein interiors. Here, we present a new microscopic method that provides an efficient way for simulating the energetics of water insertion. Instead of performing explicit Monte Carlo (MC) moves on the insertion process, which generally involves an enormous number of rejected attempts, our method is based on generating trial configurations with excess amount of internal water, estimating the relevant free energy by the linear response approximation, and then using a postprocessing MC treatment to filter out a limited number of configurations from a large possible set. Our approach is validated on particularly challenging test cases including the pK(a) of the V66D mutation in Staphylococcal nuclease, Glu286 in cytochrome c oxidase (CcO) and the energetics of a protonated water molecule in the D channel of CcO. The new postprocessing method allows us to reproduce the relevant energetics of highly unstable charges in protein interiors using fully microscopic calculations and provides a substantial improvement over regular microscopic free energy estimates. This advance established the effectiveness of our water insertion strategy in challenging cases that have not been addressed successfully by other microscopic methods. Furthermore, our study provides a new exciting view on the crucial effect of water penetration in key biological systems as well as a new view on the nature of the dielectric in protein interiors.

Molecular dynamics simulations reveal proton transfer pathways in cytochrome C-dependent nitric oxide reductase

Nitric oxide reductases (NORs) are membrane proteins that catalyze the reduction of nitric oxide (NO) to nitrous oxide (N₂O), which is a critical step of the nitrate respiration process in denitrifying bacteria. Using the recently determined first crystal structure of the cytochrome c-dependent NOR (cNOR) [Hino T, Matsumoto Y, Nagano S, Sugimoto H, Fukumori Y, et al. (2010) Structural basis of biological N2O generation by bacterial nitric oxide reductase. Science 330: 1666–70.], we performed extensive all-atom molecular dynamics (MD) simulations of cNOR within an explicit membrane/solvent environment to fully characterize water distribution and dynamics as well as hydrogen-bonded networks inside the protein, yielding the atomic details of functionally important proton channels. Simulations reveal two possible proton transfer pathways leading from the periplasm to the active site, while no pathways from the cytoplasmic side were found, consistently with the experimental observations that cNOR is not a proton pump. One of the pathways, which was newly identified in the MD simulation, is blocked in the crystal structure and requires small structural rearrangements to allow for water channel formation. That pathway is equivalent to the functional periplasmic cavity postulated in cbb₃ oxidase, which illustrates that the two enzymes share some elements of the proton transfer mechanisms and confirms a close evolutionary relation between NORs and C-type oxidases. Several mechanisms of the critical proton transfer steps near the catalytic center are proposed.

Denitrification is an anaerobic process performed by several bacteria as an alternative to aerobic respiration. A key intermediate step is catalyzed by the nitric oxide reductase (NOR) enzyme, which is situated in the cytoplasmic membrane. Proton delivery to the catalytic site inside NOR is an important part of its functioning. In this work we use molecular dynamics simulations to describe water distribution and to identify proton transfer pathways in cNOR. Our results reveal two channels from the periplasmic side of the membrane and none from the cytoplasmic side, indicating that cNOR is not a proton pump. It is our hope that these results will provide a basis for further experimental and computational studies aimed to understand details of the NOR mechanism. Furthermore, this work sheds light on the molecular evolution of respiratory enzymes.

Expanding the view of proton pumping in cytochrome c oxidase through computer simulation

In cytochrome c oxidase (CcO), a redox-driven proton pump, protons are transported by the Grotthuss shuttling via hydrogen-bonded water molecules and protonatable residues. Proton transport through the D-pathway is a complicated process that is highly sensitive to alterations in the amino acids or the solvation structure in the channel, both of which can inhibit proton pumping and enzymatic activity. Simulations of proton transport in the hydrophobic cavity showed a clear redox state dependence. To study the mechanism of proton pumping in CcO, multi-state empirical valence bond (MS-EVB) simulations have been conducted, focusing on the proton transport through the D-pathway and the hydrophobic cavity next to the binuclear center. The hydration structures, transport pathways, effects of residues, and free energy surfaces of proton transport were revealed in these MS-EVB simulations. The mechanistic insight gained from them is herein reviewed and placed in context for future studies.

Gating and regulation of the cytochrome c oxidase proton pump

As a consumer of 95% of the oxygen we breathe, cytochrome c oxidase plays a major role in the energy balance of the cell. Regulation of its oxygen reduction and proton pumping activity is therefore critical to physiological function in health and disease. The location and structure of pathways for protons that are required to support cytochrome c oxidase activity are still under debate, with respect to their requirements for key residues and fixed waters, and how they are gated to prevent (or allow) proton backflow. Recent high resolution structures of bacterial and mammalian forms reveal conserved lipid and steroid binding sites as well as redox-linked conformational changes that provide new insights into potential regulatory ligands and gating modes. Mechanistic interpretation of these findings and their significance for understanding energy regulation is discussed.

Glu-286 rotation and water wire reorientation are unlikely the gating elements for proton pumping in cytochrome C oxidase

One of the key unresolved issues regarding proton pumping in cytochrome c oxidase (CcO) is the identity of the gating element that prevents the backflow of protons. In this study, we analyze two popular proposals for this element: isomerization of the key branching residue (Glu-286) and (re)orientation of water molecules in the hydrophobic cavity. Using a multifaceted set of computational analyses that involve CcO embedded in either an implicit or explicit treatment of lipid membrane, we show that neither Glu-286 nor active-site water likely constitutes the gating element. Detailed energetic and structural analyses of the simulation results indicate that the gating-relevant properties of these structural motifs observed in previous work are likely a result of the simplified computational models employed in those studies.

High resolution structure of the ba3 cytochrome c oxidase from Thermus thermophilus in a lipidic environment

The fundamental chemistry underpinning aerobic life on Earth involves reduction of dioxygen to water with concomitant proton translocation. This process is catalyzed by members of the heme-copper oxidase (HCO) superfamily. Despite the availability of crystal structures for all types of HCO, the mode of action for this enzyme is not understood at the atomic level, namely how vectorial H⁺ and e^- transport are coupled. Toward addressing this problem, we report wild type and A120F mutant structures of the ba₃-type cytochrome c oxidase from Thermus thermophilus at 1.8 Å resolution. The enzyme has been crystallized from the lipidic cubic phase, which mimics the biological membrane environment. The structures reveal 20 ordered lipid molecules that occupy binding sites on the protein surface or mediate crystal packing interfaces. The interior of the protein encloses 53 water molecules, including 3 trapped in the designated K-path of proton transfer and 8 in a cluster seen also in A-type enzymes that likely functions in egress of product water and proton translocation. The hydrophobic O₂-uptake channel, connecting the active site to the lipid bilayer, contains a single water molecule nearest the Cu_B atom but otherwise exhibits no residual electron density. The active site contains strong electron density for a pair of bonded atoms bridging the heme Fe_a3 and Cu_B atoms that is best modeled as peroxide. The structure of ba₃-oxidase reveals new information about the positioning of the enzyme within the membrane and the nature of its interactions with lipid molecules. The atomic resolution details provide insight into the mechanisms of electron transfer, oxygen diffusion into the active site, reduction of oxygen to water, and pumping of protons across the membrane. The development of a robust system for production of ba₃-oxidase crystals diffracting to high resolution, together with an established expression system for generating mutants, opens the door for systematic structure-function studies.

Molecular basis of proton uptake in single and double mutants of cytochrome c oxidase

Cytochrome c oxidase, the terminal enzyme of the respiratory chain, utilizes the reduction of dioxygen into water to pump protons across the mitochondrial inner membrane. The principal pathway of proton uptake into the enzyme, the D channel, is a 2.5 nm long channel-like cavity named after a conserved, negatively charged aspartic acid (D) residue thought to help recruiting protons to its entrance (D132 in the first subunit of the S. sphaeroides enzyme). The single-point mutation of D132 to asparagine (N), a neutral residue, abolishes enzyme activity. Conversely, replacing conserved N139, one-third into the D channel, by D, induces a decoupled phenotype, whereby oxygen reduction proceeds but not proton pumping. Intriguingly, the double mutant D132N/N139D, which conserves the charge of the D channel, restores the wild-type phenotype. We use molecular dynamics simulations and electrostatic calculations to examine the structural and physical basis for the coupling of proton pumping and oxygen chemistry in single and double N139D mutants. The potential of mean force for the conformational isomerization of N139 and N139D side chains reveals the presence of three rotamers, one of which faces the channel entrance. This out-facing conformer is metastable in the wild-type and in the N139D single mutant, but predominant in the double mutant thanks to the loss of electrostatic repulsion with the carboxylate group of D132. The effects of mutations and conformational isomerization on the pKa of E286, an essential proton-shuttling residue located at the top of the D channel, are shown to be consistent with the electrostatic control of proton pumping proposed recently (Fadda et al 2008 Biochim. Biophys. Acta 1777 277-84). Taken together, these results suggest that preserving the spatial distribution of charges at the entrance of the D channel is necessary to guarantee both the uptake and the relay of protons to the active site of the enzyme. These findings highlight the interplay of long-range electrostatic forces and local structural fluctuations in the control of proton movement and provide a physical explanation for the restoration of proton pumping activity in the double mutant.

The D-channel of cytochrome oxidase: an alternative view

The D-pathway in A-type cytochrome c oxidases conducts protons from a conserved aspartate on the negatively charged N-side of the membrane to a conserved glutamic acid at about the middle of the membrane dielectric. Extensive work in the past has indicated that all four protons pumped across the membrane on reduction of O(2) to water are transferred via the D-pathway, and that it is also responsible for transfer of two out of the four "chemical protons" from the N-side to the binuclear oxygen reduction site to form product water. The function of the D-pathway has been discussed in terms of an apparent pK(a) of the glutamic acid. After reacting fully reduced enzyme with O(2), the rate of formation of the F state of the binuclear heme-copper active site was found to be independent of pH up to pH~9, but to drop off at higher pH with an apparent pK(a) of 9.4, which was attributed to the glutamic acid. Here, we present an alternative view, according to which the pH-dependence is controlled by proton transfer into the aspartate residue at the N-side orifice of the D-pathway. We summarise experimental evidence that favours a proton pump mechanism in which the proton to be pumped is transferred from the glutamic acid to a proton-loading site prior to proton transfer for completion of oxygen reduction chemistry. The mechanism is discussed by which the proton-pumping activity is decoupled from electron transfer by structural alterations of the D-pathway. This article is part of a Special Issue entitled: Allosteric cooperativity in respiratory proteins.

Modeling and simulation of ion-coupled and ATP-driven membrane proteins

The molecular mechanisms of membrane proteins that are activated either by ions or by ATP are just beginning to come into focus, as long-awaited structural data are revealed. This information is being leveraged and supplemented to great effect by molecular modeling and computer simulation studies. Important examples include the homology modeling of eukaryotic protein structures based on distantly related templates, as well as the use of internal structural symmetry for modeling different states in conformational cycles. Molecular simulation studies have elucidated the location and coordination structure of ion binding sites, and explained their selectivity, while also providing tantalizing insights into the mechanisms that couple conformational change to ion translocation or ATP hydrolysis.