Molecular mechanisms for generating transmembrane proton gradients

Redox-Activated Proton Transfer through a Redundant Network in the Qo Site of Cytochrome bc1

Proton translocation catalyzed by cytochrome bc1 (respiratory complex III) during coenzyme-Q redox cycling is a critical bioenergetic process, yet its detailed molecular mechanism remains incompletely understood. In this study, the energetics of the long-range proton transfers through multiple proton-conducting wires in the Qo site of the bc1 complex was investigated computationally using hybrid QM/MM simulations and a specialized reaction coordinate. Key reactive groups and proton transfer mechanisms were characterized, confirming the propionate-A group of heme bL as a plausible proton acceptor. Upon coenzyme-Q oxidation, a Grotthuss hopping mechanism is activated, facilitating proton transfer along three distinct pathways with comparable barriers and stability. These pathways operate redundantly, forming a robust proton-conducting network, and account for the unusual experimental behavior observed in single-point mutations. Energetic analyses exclude charged closed-shell species as likely intermediates and propose a reaction sequence for coenzyme-Q oxidation proceeding as QH2 → QH• → Q0, either via coupled proton-electron transfers or stepwise mechanisms involving open-shell intermediates. These findings elucidate mechanistic details of the Q-cycle and improve our understanding of the catalytic reactions supporting redox-activated proton transfer in respiratory enzymes.

Finding the E-channel proton loading sites by calculating the ensemble of protonation microstates

The aerobic electron transfer chain builds a proton gradient by proton coupled electron transfer reactions through a series of proteins. Complex I is the first enzyme in the sequence. Here transfer of two electrons from NADH to quinone yields four protons pumped from the membrane N- (negative, higher pH) side to the P- (positive, lower pH) side. Protons move through three linear antiporter paths, with a few amino acids and waters providing the route; and through the E-channel, a complex of competing paths, with clusters of interconnected protonatable residues. Proton loading sites (PLS) transiently bind protons as they are transported from N- to P-compartments. PLS can be individual residues or extended clusters of residues. The program MCCE uses Monte Carlos sampling to analyze the E-channel proton binding in equilibrium with individual Molecular Dynamics snapshots from trajectories of Thermus thermuphillus Complex I in the apo, quinone and quinol bound states. At pH 7, the five E-channel subunits (Nqo4, Nqo7, Nqo8, Nqo10, and Nqo11) take >25,000 protonation microstates, each with different residues protonated. The microstate explosion is tamed by analyzing interconnected clusters of residues along the proton transfer paths. A proton is bound and released from a cluster of five coupled residues on the protein N-side and to six coupled residues in the protein center. Loaded microstates bind protons to sites closer to the P-side in the forward pumping direction. MCCE microstate analysis identifies strongly coupled proton binding amongst individual residues in the two PLS clusters.

Similarity between oxygen evolution in photosystem II and oxygen reduction in cytochrome c oxidase via proton coupled electron transfers. A unified view of the oxygenic life from four electron oxidation-reduction reactions

Basic concepts and theoretical foundations of broken symmetry (BS) and post BS methods for strongly correlated electron systems (SCES) such as electron-transfer (ET) diradical, multi-center polyradicals with spin frustration are described systematically to elucidate structures, bonding and reactivity of the high-valent transition metal oxo bonds in metalloenzymes: photosystem II (PSII) and cytochrome c oxidase (CcO). BS hybrid DFT (HDFT) and DLPNO coupled-cluster (CC) SD(T0) computations are performed to elucidate electronic and spin states of CaMn4Ox cluster in the key step for oxygen evolution, namely S4 [S3 with Mn(IV) = O + Tyr161-O radical] state of PSII and PM [Fe(IV) = O + HO-Cu(II) + Tyr161-O radical] step for oxygen reduction in CcO. The cycle of water oxidation catalyzed by the CaMn4Ox cluster in PSII and the cycle of oxygen reduction catalyzed by the CuA-Fea-Fea3-CuB cluster in CcO are examined on the theoretical grounds, elucidating similar concerted and/or stepwise proton transfer coupled electron transfer (PT-ET) processes for the four-electron oxidation in PSII and four-electron reduction in CcO. Interplay between theory and experiments have revealed that three electrons in the metal sites and one electron in tyrosine radical site are characteristic for PT-ET in these biological redox reaction systems, indicating no necessity of harmful Mn(V) = O and Fe(V) = O bonds with strong oxyl-radical character. Implications of the computational results are discussed in relation to design of artificial systems consisted of earth abundant transition metals for water oxidation.

Hybrid Quantum Mechanical/Molecular Mechanical Methods For Studying Energy Transduction in Biomolecular Machines

Hybrid quantum mechanical/molecular mechanical (QM/MM) methods have become indispensable tools for the study of biomolecules. In this article, we briefly review the basic methodological details of QM/MM approaches and discuss their applications to various energy transduction problems in biomolecular machines, such as long-range proton transports, fast electron transfers, and mechanochemical coupling. We highlight the particular importance for these applications of balancing computational efficiency and accuracy. Using several recent examples, we illustrate the value and limitations of QM/MM methodologies for both ground and excited states, as well as strategies for calibrating them in specific applications. We conclude with brief comments on several areas that can benefit from further efforts to make QM/MM analyses more quantitative and applicable to increasingly complex biological problems.

Tools for analyzing protonation states and for tracing proton transfer pathways with examples from the Rb. sphaeroides photosynthetic reaction centers

Protons participate in many reactions. In proteins, protons need paths to move in and out of buried active sites. The vectorial movement of protons coupled to electron transfer reactions establishes the transmembrane electrochemical gradient used for many reactions, including ATP synthesis. Protons move through hydrogen bonded chains of waters and hydroxy side chains via the Grotthuss mechanism and by proton binding and release from acidic and basic residues. MCCE analysis shows that proteins exist in a large number of protonation states. Knowledge of the equilibrium ensemble can provide a rational basis for setting protonation states in simulations that fix them, such as molecular dynamics (MD). The proton path into the Q_B site in the bacterial reaction centers (RCs) of Rb. sphaeroides is analyzed by MD to provide an example of the benefits of using protonation states found by the MCCE program. A tangled web of side chains and waters link the cytoplasm to Q_B. MCCE analysis of snapshots from multiple trajectories shows that changing the input protonation state of a residue in MD biases the trajectory shifting the proton affinity of that residue. However, the proton affinity of some residues is more sensitive to the input structure. The proton transfer networks derived from different trajectories are quite robust. There are some changes in connectivity that are largely restricted to the specific residues whose protonation state is changed. Trajectories with Q_B^•- are compared with earlier results obtained with Q_B [Wei et. al Photosynthesis Research volume 152, pages153-165 (2022)] showing only modest changes. While introducing new methods the study highlights the difficulty of establishing the connections between protein conformation.

Identification of amino acid residues in a proton release pathway near the bacteriochlorophyll dimer in reaction centers from Rhodobacter sphaeroides

Insight into control of proton transfer, a crucial attribute of cellular functions, can be gained from investigations of bacterial reaction centers. While the uptake of protons associated with the reduction of the quinone is well characterized, the release of protons associated with the oxidized bacteriochlorophyll dimer has been poorly understood. Optical spectroscopy and proton release/uptake measurements were used to examine the proton release characteristics of twelve mutant reaction centers, each containing a change in an amino acid residue near the bacteriochlorophyll dimer. The mutant reaction centers had optical spectra similar to wild-type and were capable of transferring electrons to the quinones after light excitation of the bacteriochlorophyll dimer. They exhibited a large range in the extent of proton release and in the slow recovery of the optical signal for the oxidized dimer upon continuous illumination. Key roles were indicated for six amino acid residues, Thr L130, Asp L155, Ser L244, Arg M164, Ser M190, and His M193. Analysis of the results points to a hydrogen-bond network that contains these residues, with several additional residues and bound water molecules, forming a proton transfer pathway. In addition to proton transfer, the properties of the pathway are proposed to be responsible for the very slow charge recombination kinetics observed after continuous illumination. The characteristics of this pathway are compared to proton transfer pathways near the secondary quinone as well as those found in photosystem II and cytochrome c oxidase.

Characterizing Protein Protonation Microstates Using Monte Carlo Sampling

Proteins are polyelectrolytes with acidic and basic amino acids Asp, Glu, Arg, Lys, and His, making up ≈25% of the residues. The protonation state of residues, cofactors, and ligands defines a "protonation microstate". In an ensemble of proteins some residues will be ionized and others neutral, leading to a mixture of protonation microstates rather than in a single one as is often assumed. The microstate distribution changes with pH. The protein environment also modifies residue proton affinity so microstate distributions change in different reaction intermediates or as ligands are bound. Particular protonation microstates may be required for function, while others exist simply because there are many states with similar energy. Here, the protonation microstates generated in Monte Carlo sampling in MCCE are characterized in HEW lysozyme as a function of pH and bacterial photosynthetic reaction centers (RCs) in different reaction intermediates. The lowest energy and highest probability microstates are compared. The ΔG, ΔH, and ΔS between the four protonation states of Glu35 and Asp52 in lysozyme are shown to be calculated with reasonable precision. At pH 7 the lysozyme charge ranges from 6 to 10, with 24 accepted protonation microstates, while RCs have ≈50,000. A weighted Pearson correlation analysis shows coupling between residue protonation states in RCs and how they change when the quinone in the Q_B site is reduced. Protonation microstates can be used to define input MD parameters and provide insight into the motion of protons coupled to reactions.

O to bR transition in bacteriorhodopsin occurs through a proton hole mechanism

Extensive classical and quantum mechanical/molecular mechanical (QM/MM) molecular dynamics simulations are used to establish the structural features of the O state in bacteriorhodopsin (bR) and its conversion back to the bR ground state. The computed free energy surface is consistent with available experimental data for the kinetics and thermodynamics of the O to bR transition. The simulation results highlight the importance of the proton release group (PRG, consisting of Glu194/204) and the conserved arginine 82 in modulating the hydration level of the protein cavity. In particular, in the O state, deprotonation of the PRG and downward rotation of Arg82 lead to elevated hydration level and a continuous water network that connects the PRG to the protonated Asp85. Proton exchange through this water network is shown by ∼0.1-μs semiempirical QM/MM free energy simulations to occur through the generation and propagation of a proton hole, which is relayed by Asp212 and stabilized by Arg82. This mechanism provides an explanation for the observation that the D85S mutant of bacteriorhodopsin pumps chloride ions. The electrostatics-hydration coupling mechanism and the involvement of all titration states of water are likely applicable to many biomolecules involved in bioenergetic transduction.

Activating the Cpx response induces tolerance to antisense PNA delivered by an arginine-rich peptide in Escherichia coli

Cell-penetrating peptides (CPPs) are increasingly used for cellular drug delivery in both pro- and eukaryotic cells, and oligoarginines have attracted special attention. How arginine-rich CPPs translocate across the cell envelope, particularly for prokaryotes, is still unknown. Arginine-rich CPPs efficiently deliver antimicrobial peptide nucleic acid (PNA) to its intracellular mRNA target in bacteria. We show that resistance to PNA conjugated to an arginine-rich CPP in Escherichia coli requires multiple genetic modifications and is specific for the CPP part and not to the PNA part. An integral part of the resistance was the constitutively activated Cpx-envelope stress response system (cpx∗), which decreased the cytoplasmic membrane potential. This indicates an indirect energy-dependent uptake mechanism for antimicrobials conjugated to arginine-rich CPPs. In agreement, cpx∗ mutants showed low-level resistance to aminoglycosides and an arginine-rich CPP conjugated to a peptide targeting the DNA sliding clamp, i.e., similar uptake in E. coli for these antimicrobial compounds.

"Ion Pool" Structural Ion Storage Device: A New Strategy to Collect Ions by Nanoconfinement Effects

Ion storage structure widely exists in organisms, which is used to harvesting energy in environment and converting it into ion concentration gradient to maintain complex life activities. The construction of ion storage structures relies on isolating the biological body fluids by biofilm systems, which can also be regarded as local ions confinement. Mimicking this ions storage process, an "ion pool" structural ion storage device is proposed in this research by artificial ion nanochannels, which can transform the electric power into ion concentration gradient. It is consisted of micrometer-sized ions reservoir and nanosized ions filters. Ions can be isolated within the "pool" and performed ultrahigh ions enrichment or depletion behavior deviated from bulk. Through numerical simulation by Poisson-Nernst-Planck equations, "ion pool" structural device achieves nearly 20 000 rectification ratio with low surface charge. An "ion pool" structural ions storage device is also constructed with block copolymer and polyethylene terephthalate composite membranes, a super high rectification ratio of 3184.0 is achieved from the experiment, which is the highest reported so far. The ion storage efficiency of the device reaches 14.90%, which is an order of magnitude better than non-"ion pool" structural nanofluid devices.

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.

Proton exit pathways surrounding the oxygen evolving complex of photosystem II

Photosystem II allows water to be the primary electron source for the photosynthetic electron transfer chain. Water is oxidized to dioxygen at the Oxygen Evolving Complex (OEC), a Mn₄CaO₅ inorganic core embedded on the lumenal side of PSII. Water-filled channels surrounding the OEC must bring in substrate water molecules, remove the product protons to the lumen, and may transport the product oxygen. Three water-filled channels, denoted large, narrow, and broad, extend from the OEC towards the aqueous surface more than 15 Å away. However, the role of each pathway in the transport in and out of the OEC is yet to be established. Here, we combine Molecular Dynamics (MD), Multi Conformation Continuum Electrostatics (MCCE) and Network Analysis to compare and contrast the three potential proton transfer paths. Hydrogen bond network analysis shows that near the OEC the waters are highly interconnected with similar free energy for hydronium at all locations. The paths diverge as they move towards the lumen. The water chain in the broad channel is better connected than in the narrow and large channels, where disruptions in the network are observed approximately 10 Å from the OEC. In addition, the barrier for hydronium translocation is lower in the broad channel. Thus, a proton released from any location on the OEC can access all paths, but the likely exit to the lumen passes through PsbO via the broad channel.

Identifying the proton loading site cluster in the ba3 cytochrome c oxidase that loads and traps protons

Cytochrome c Oxidase (CcO) is the terminal electron acceptor in aerobic respiratory chain, reducing O₂ to water. The released free energy is stored by pumping protons through the protein, maintaining the transmembrane electrochemical gradient. Protons are held transiently in a proton loading site (PLS) that binds and releases protons driven by the electron transfer reaction cycle. Multi-Conformation Continuum Electrostatics (MCCE) was applied to crystal structures and Molecular Dynamics snapshots of the B-type Thermus thermophilus CcO. Six residues are identified as the PLS, binding and releasing protons as the charges on heme b and the binuclear center are changed: the heme a₃ propionic acids, Asp287, Asp372, His376 and Glu126B. The unloaded state has one proton and the loaded state two protons on these six residues. Different input structures, modifying the PLS conformation, show different proton distributions and result in different proton pumping behaviors. One loaded and one unloaded protonation states have the loaded/unloaded states close in energy so the PLS binds and releases a proton through the reaction cycle. The alternative proton distributions have state energies too far apart to be shifted by the electron transfers so are locked in loaded or unloaded states. Here the protein can use active states to load and unload protons, but has nearby trapped states, which stabilize PLS protonation state, providing new ideas about the CcO proton pumping mechanism. The distance between the PLS residues Asp287 and His376 correlates with the energy difference between loaded and unloaded states.

Hydrogen bond network analysis reveals the pathway for the proton transfer in the E-channel of T. thermophilus Complex I

Complex I, NADH-ubiquinone oxidoreductase, is the first enzyme in the mitochondrial and bacterial aerobic respiratory chain. It pumps four protons through four transiently open pathways from the high pH, negative, N-side of the membrane to the positive, P-side driven by the exergonic transfer of electrons from NADH to a quinone. Three protons transfer through subunits descended from antiporters, while the fourth, E-channel is unique. The path through the E-channel is determined by a network analysis of hydrogen bonded pathways obtained by Monte Carlo sampling of protonation states, polar hydrogen orientation and water occupancy. Input coordinates are derived from molecular dynamics trajectories comparing oxidized, reduced (dihydro) and no menaquinone-8 (MQ). A complex proton transfer path from the N- to the P-side is found consisting of six clusters of highly connected hydrogen-bonded residues. The network connectivity depends on the presence of quinone and its redox state, supporting a role for this cofactor in coupling electron and proton transfers. The N-side is more organized with MQ-bound complex I facilitating proton entry, while the P-side is more connected in the apo-protein, facilitating proton exit. Subunit Nqo8 forms the core of the E channel; Nqo4 provides the N-side entry, Nqo7 and then Nqo10 join the pathway in the middle, while Nqo11 contributes to the P-side exit.

Potential of Engineered Bacteriorhodopsins as Photoactivated Biomaterials in Modulating Neural Stem Cell Behavior

Bacteriorhodopsin (BR), a light-sensitive bacterial proton pump, has been demonstrated the capacity for regulating the neural activity in mammalian cells. Because of the difficulty in production and purification in large quantities, the BR proteins have neither been directly employed to biomedical applications nor verified the functionality by protein administration. Previously, we have invented a highly expressible bacteriorhodopsin (HEBR) and established the massive production protocol. In the current study, we mass-produced the two types of HEBR proteins that have normal or abnormal activity on the proton pumping, and then we treated murine neural stem cells (NSCs) with these HEBR proteins. We discovered that the cell behaviors including growth, metabolism, mitochondrial inner membrane potential, and differentiation were obviously affected in NSCs after the treatment of HEBR proteins. Particularly, these effects induced by HEBR proteins were correlated to their proton pump activity and could be altered by cell culture substrate materials. Current findings suggest that the engineered light-sensitive HEBR protein can serve as a biological material to directly influence the multiple behaviors of mammalian cells, which is further modified by the cell culture substrate material, revealing the versatile potential of HEBR protein in biomaterial applications.

Tracing the Pathways of Waters and Protons in Photosystem II and Cytochrome c Oxidase

Photosystem II (PSII) uses water as the terminal electron donor, producing oxygen in the Mn4CaO5 oxygen evolving complex (OEC), while cytochrome c oxidase (CcO) reduces O2 to water in its heme–Cu binuclear center (BNC). Each protein is oriented in the membrane to add to the proton gradient. The OEC, which releases protons, is located near the P-side (positive, at low-pH) of the membrane. In contrast, the BNC is in the middle of CcO, so the protons needed for O2 reduction must be transferred from the N-side (negative, at high pH). In addition, CcO pumps protons from N- to P-side, coupled to the O2 reduction chemistry, to store additional energy. Thus, proton transfers are directly coupled to the OEC and BNC redox chemistry, as well as needed for CcO proton pumping. The simulations that study the changes in proton affinity of the redox active sites and the surrounding protein at different states of the reaction cycle, as well as the changes in hydration that modulate proton transfer paths, are described.

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.

Effect of a pH Gradient on the Protonation States of Cytochrome c Oxidase: A Continuum Electrostatics Study

Cytochrome c oxidase (CcO) couples the reduction of dioxygen to water with transmembrane proton pumping, which leads to the generation of an electrochemical gradient. In this study we analyze how one of the components of the electrochemical gradient, the difference in pH across the membrane, or ΔpH, influences the protonation states of residues in CcO. We modified our continuum electrostatics/Monte Carlo (CE/MC) method in order to include the ΔpH and applied it to the study of CcO, in what is, to our best knowledge, the first CE/MC study of CcO in the presence of a pH gradient. The inclusion of a transmembrane pH gradient allows for the identification of residues whose titration behavior depends on the pH on both sides of the membrane. Among the several residues with unusual titration profiles, three are well-known key residues in the proton transfer process of CcO: E286_I, Y288_I, and K362_I. All three residues have been previously identified as being critical for the catalytic or proton pumping functions of CcO. Our results suggest that when the pH gradient increases, these residues may be part of a regulatory mechanism to stem the proton flow.

Design and fine-tuning redox potentials of metalloproteins involved in electron transfer in bioenergetics

Redox potentials are a major contributor in controlling the electron transfer (ET) rates and thus regulating the ET processes in the bioenergetics. To maximize the efficiency of the ET process, one needs to master the art of tuning the redox potential, especially in metalloproteins, as they represent major classes of ET proteins. In this review, we first describe the importance of tuning the redox potential of ET centers and its role in regulating the ET in bioenergetic processes including photosynthesis and respiration. The main focus of this review is to summarize recent work in designing the ET centers, namely cupredoxins, cytochromes, and iron-sulfur proteins, and examples in design of protein networks involved these ET centers. We then discuss the factors that affect redox potentials of these ET centers including metal ion, the ligands to metal center and interactions beyond the primary ligand, especially non-covalent secondary coordination sphere interactions. We provide examples of strategies to fine-tune the redox potential using both natural and unnatural amino acids and native and nonnative cofactors. Several case studies are used to illustrate recent successes in this area. Outlooks for future endeavors are also provided. This article is part of a Special Issue entitled Biodesign for Bioenergetics--the design and engineering of electronic transfer cofactors, proteins and protein networks, edited by Ronald L. Koder and J.L. Ross Anderson.

Structural Changes and Proton Transfer in Cytochrome c Oxidase

In cytochrome c oxidase electron transfer from cytochrome c to O2 is linked to transmembrane proton pumping, which contributes to maintaining a proton electrochemical gradient across the membrane. The mechanism by which cytochrome c oxidase couples the exergonic electron transfer to the endergonic proton translocation is not known, but it presumably involves local structural changes that control the alternating proton access to the two sides of the membrane. Such redox-induced structural changes have been observed in X-ray crystallographic studies at residues 423-425 (in the R. sphaeroides oxidase), located near heme a. The aim of the present study is to investigate the functional effects of these structural changes on reaction steps associated with proton pumping. Residue Ser425 was modified using site-directed mutagenesis and time-resolved spectroscopy was used to investigate coupled electron-proton transfer upon reaction of the oxidase with O2. The data indicate that the structural change at position 425 propagates to the D proton pathway, which suggests a link between redox changes at heme a and modulation of intramolecular proton-transfer rates.

Proton-Coupled Electron Transfer During the S-State Transitions of the Oxygen-Evolving Complex of Photosystem II

The oxygen-evolving complex (OEC) of photosystem II (PSII) is a unique Mn4O5Ca cluster that catalyzes water oxidation via four photoactivated electron transfer steps. As the protein influence on the redox and protonation chemistry of the OEC remains an open question, we present a classical valence model of the OEC that allows the redox state of each Mn and the protonation state of bridging μ-oxos and terminal waters to remain in equilibrium with the PSII protein throughout the redox cycle. We find that the last bridging oxygen loses its proton during the transition from S0 to S1. Two possible S2 states are found depending on the OEC geometry: S2 has Mn4(IV) with a proton lost from a terminal water (W1) trapped by the nearby D1-D61 if O5 is closer to Mn4, or Mn1(IV), with partial deprotonation of D1-H337 and D1-E329 if O5 is closer to Mn1. In S3, the OEC is Mn4(IV) with W2 deprotonated. The estimated OEC Em's range from +0.7 to +1.3 V, enabling oxidation by P680(+), the primary electron donor in PSII. In chloride-depleted PSII, the proton release increases during the S1 to S2 transition, leaving the OEC unable to properly advance through the water-splitting cycle.

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.

Halorhodopsin pumps Cl- and bacteriorhodopsin pumps protons by a common mechanism that uses conserved electrostatic interactions

Key mutations differentiate the functions of homologous proteins. One example compares the inward ion pump halorhodopsin (HR) and the outward proton pump bacteriorhodopsin (BR). Of the nine essential buried ionizable residues in BR, six are conserved in HR. However, HR changes three BR acids, D85 in a central cluster of ionizable residues, D96, nearer the intracellular, and E204, nearer the extracellular side of the membrane to the small, neutral amino acids T111, V122, and T230, respectively. In BR, acidic amino acids are stationary anions whose proton affinity is modulated by conformational changes, establishing a sequence of directed binding and release of protons. Multiconformation continuum electrostatics calculations of chloride affinity and residue protonation show that, in reaction intermediates where an acid is ionized in BR, a Cl(-) is bound to HR in a position near the deleted acid. In the HR ground state, Cl(-) binds tightly to the central cluster T111 site and weakly to the extracellular T230 site, recovering the charges on ionized BR-D85 and neutral E204 in BR. Imposing key conformational changes from the BR M intermediate into the HR structure results in the loss of Cl(-) from the central T111 site and the tight binding of Cl(-) to the extracellular T230 site, mirroring the changes that protonate BR-D85 and ionize E204 in BR. The use of a mobile chloride in place of D85 and E204 makes HR more susceptible to the environmental pH and salt concentrations than BR. These studies shed light on how ion transfer mechanisms are controlled through the interplay of protein and ion electrostatics.

Characterizing the proton loading site in cytochrome c oxidase

Cytochrome c oxidase (CcO) uses the energy released by reduction of O2 to H2O to drive eight charges from the high pH to low pH side of the membrane, increasing the electrochemical gradient. Four electrons and protons are used for chemistry, while four more protons are pumped. Proton pumping requires that residues on a pathway change proton affinity through the reaction cycle to load and then release protons. The protonation states of all residues in CcO are determined in MultiConformational Continuum Electrostatics simulations with the protonation and redox states of heme a, a3, Cu(B), Y288, and E286 used to define the catalytic cycle. One proton is found to be loaded and released from residues identified as the proton loading site (PLS) on the P-side of the protein in each of the four CcO redox states. Thus, the same proton pumping mechanism can be used each time CcO is reduced. Calculations with structures of Rhodobacter sphaeroides, Paracoccus denitrificans, and bovine CcO derived by crystallography and molecular dynamics show the PLS functions similarly in different CcO species. The PLS is a cluster rather than a single residue, as different structures show 1-4 residues load and release protons. However, the proton affinity of the heme a3 propionic acids primarily determines the number of protons loaded into the PLS; if their proton affinity is too low, less than one proton is loaded.

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.

Mutational control of bioenergetics of bacterial reaction center probed by delayed fluorescence

The free energy gap between the metastable charge separated state P(+)QA(-) and the excited bacteriochlorophyll dimer P* was measured by delayed fluorescence of the dimer in mutant reaction center proteins of the photosynthetic bacterium Rhodobacter sphaeroides. The mutations were engineered both at the donor (L131L, M160L, M197F and M202H) and acceptor (M265I and M234E) sides. While the donor side mutations changed systematically the number of H-bonds to P, the acceptor side mutations modified the energetics of QA by altering the van-der-Waals and electronic interactions (M265IT) and H-bond network to the acidic cluster around QB (M234EH, M234EL, M234EA and M234ER). All mutants decreased the free energy gap of the wild type RC (~890meV), i.e. destabilized the P(+)QA(-) charge pair by 60-110meV at pH8. Multiple modifications in the hydrogen bonding pattern to P resulted in systematic changes of the free energy gap. The destabilization showed no pH-dependence (M234 mutants) or slight increase (WT, donor-side mutants and M265IT above pH8) with average slope of 10-15meV/pH unit over the 6-10.5pH range. In wild type and donor-side mutants, the free energy change of the charge separation consisted of mainly enthalpic term but the acceptor side mutants showed increased entropic (even above that of enthalpic) contributions. This could include softening the structure of the iron ligand (M234EH) and the QA binding pocket (M265IT) and/or increase of the multiplicity of the electron transfer of charge separation in the acceptor side upon mutation.