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

Searching for proton transfer channels in respiratory complex I

We have explored a strategy to identify potential proton transfer channels using computational analysis of a protein structure based on Voronoi partitioning and applied it for the analysis of proton transfer pathways in redox-driven proton-pumping respiratory complex I. The analysis results in a network of connected voids/channels, which represent the dual structure of the protein; we then hydrated the identified channels using our water placement program Dowser++. Many theoretical water molecules found in the channels perfectly match the observed experimental water molecules in the structure; some other predicted water molecules have not been resolved in the experiments. The channels are of varying cross sections. Some channels are big enough to accommodate water molecules that are suitable to conduct protons; others are too narrow to hold water but require only minor conformational changes to accommodate proton transfer. We provide a preliminary analysis of the proton conductivity of the network channels, classifying the proton transfer channels as open, closed, and partially open, and discuss possible conformational changes that can modulate, i.e., open and close, the channels.

Identification and characterization of a new pathologic mutation in a large Leber hereditary optic neuropathy pedigree

BACKGROUND

Most patients suffering from Leber hereditary optic neuropathy carry one of the three classic pathologic mutations, but not all individuals with these genetic alterations develop the disease. There are different risk factors that modify the penetrance of these mutations. The remaining patients carry one of a set of very rare genetic variants and, it appears that, some of the risk factors that modify the penetrance of the classical pathologic mutations may also affect the phenotype of these other rare mutations.

RESULTS

We describe a large family including 95 maternally related individuals, showing 30 patients with Leber hereditary optic neuropathy. The mutation responsible for the phenotype is a novel transition, m.3734A > G, in the mitochondrial gene encoding the ND1 subunit of respiratory complex I. Molecular-genetic, biochemical and cellular studies corroborate the pathogenicity of this genetic change.

CONCLUSIONS

With the study of this family, we confirm that, also for this very rare mutation, sex and age are important factors modifying penetrance. Moreover, this pedigree offers an excellent opportunity to search for other genetic or environmental factors that additionally contribute to modify penetrance.

Role of Protonation States in the Stability of Molecular Dynamics Simulations of High-Resolution Membrane Protein Structures

Classical molecular dynamics (MD) simulations provide unmatched spatial and time resolution of protein structure and function. However, the accuracy of MD simulations often depends on the quality of force field parameters and the time scale of sampling. Another limitation of conventional MD simulations is that the protonation states of titratable amino acid residues remain fixed during simulations, even though protonation state changes coupled to conformational dynamics are central to protein function. Due to the uncertainty in selecting protonation states, classical MD simulations are sometimes performed with all amino acids modeled in their standard charged states at pH 7. Here, we performed and analyzed classical MD simulations on high-resolution cryo-EM structures of two large membrane proteins that transfer protons by catalyzing protonation/deprotonation reactions. In simulations performed with titratable amino acids modeled in their standard protonation (charged) states, the structure diverges far from its starting conformation. In comparison, MD simulations performed with predetermined protonation states of amino acid residues reproduce the structural conformation, protein hydration, and protein-water and protein-protein interactions of the structure much better. The results support the notion that it is crucial to perform basic protonation state calculations, especially on structures where protonation changes play an important functional role, prior to the launch of any conventional MD simulations. Furthermore, the combined approach of fast protonation state prediction and MD simulations can provide valuable information about the charge states of amino acids in the cryo-EM sample. Even though accurate prediction of protonation states in proteinaceous environments currently remains a challenge, we introduce an approach of combining pKa prediction with cryo-EM density map analysis that helps in improving not only the protonation state predictions but also the atomic modeling of density data.

Enticing a Proton using Single Ammonia Molecule as Bait

In microhydrated acid-solvent clusters, deprotonation of an acid is assisted by a critical number of solvent molecules and a solvent electric field. Born-Oppenheimer molecular dynamics simulations reveal that trifluoroacetic acid undergoes spontaneous proton transfer in water clusters, with the critical number being five. Acetic acid and phenol, on the other hand, do not dissociate even in the presence of a large number of water molecules (in excess of 40). The addition of a single ammonia molecule to the water cluster, which interacts directly with the protic group, lowers the critical number of solvent water molecules required for proton transfer to three and seven in the case of acetic acid and phenol, respectively. The population of the undissociated and the proton-transferred structures get dispersed to form separate islands on the electric field versus the O-H distance representation with the cusp representing the critical values. The critical electric fields for the spontaneous proton transfer are around 254, 237, and 318 MV cm-1 for trifluoroacetic acid, acetic acid, and phenol, respectively. In the case of phenol, the free energy profiles suggest that proton transfer to the ammonia moiety embedded in water promotes proton transfer efficiently due to the higher basicity of ammonia and enhanced hydrogen bonding network of solvent water, vis-à-vis phenol-ammonia clusters.

Investigation of hydrated channels and proton pathways in a high-resolution cryo-EM structure of mammalian complex I

Respiratory complex I, a key enzyme in mammalian metabolism, captures the energy released by reduction of ubiquinone by NADH to drive protons across the inner mitochondrial membrane, generating the proton-motive force for ATP synthesis. Despite remarkable advances in structural knowledge of this complicated membrane-bound enzyme, its mechanism of catalysis remains controversial. In particular, how ubiquinone reduction is coupled to proton pumping and the pathways and mechanisms of proton translocation are contested. We present a 2.4-Å resolution cryo-EM structure of complex I from mouse heart mitochondria in the closed, active (ready-to-go) resting state, with 2945 water molecules modeled. By analyzing the networks of charged and polar residues and water molecules present, we evaluate candidate pathways for proton transfer through the enzyme, for the chemical protons for ubiquinone reduction, and for the protons transported across the membrane. Last, we compare our data to the predictions of extant mechanistic models, and identify key questions to answer in future work to test them.

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.

Tunnel dynamics of quinone derivatives and its coupling to protein conformational rearrangements in respiratory complex I

Respiratory complex I in mitochondria and bacteria catalyzes the transfer of electrons from NADH to quinone (Q). The free energy available from the reaction is used to pump protons and to establish a membrane proton electrochemical gradient, which drives ATP synthesis. Even though several high-resolution structures of complex I have been resolved, how Q reduction is linked with proton pumping, remains unknown. Here, microsecond long molecular dynamics (MD) simulations were performed on Yarrowia lipolytica complex I structures where Q molecules have been resolved in the ~30 Å long Q tunnel. MD simulations of several different redox/protonation states of Q reveal the coupling between the Q dynamics and the restructuring of conserved loops and ion pairs. Oxidized quinone stabilizes towards the N2 FeS cluster, a binding mode not previously described in Yarrowia lipolytica complex I structures. On the other hand, reduced (and protonated) species tend to diffuse towards the Q binding sites closer to the tunnel entrance. Mechanistic and physiological relevance of these results are discussed.

Subunits E-F-G of E. coli Complex I can form an active complex when expressed alone, but in time-delayed assembly co-expression of B-CD-E-F-G is optimal

Respiratory Complex I from E. coli is a proto-type of the mitochondrial enzyme, consisting of a 6-subunit peripheral arm (B-CD-E-F-G-I) and a 7-subunit membrane arm. When subunits E-F-G (N-module), were expressed alone they formed an active complex as determined by co-immunoprecipitation and native gel electrophoresis. When co-expressed with subunits B and CD, only a complex of E-F-G was found. When these five subunits were co-expressed with subunit I and two membrane subunits, A and H, a complex of B-CD-E-F-G-I was membrane-bound, constituting the N- and Q-modules. Assembly of Complex I was also followed by splitting the genes between two plasmids, in three different groupings, and expressing them simultaneously, or with time-delay of expression from one plasmid. When the B-CD-E-F-G genes were co-expressed after a time-delay, assembly was over 90 % of that when the whole operon was expressed together. In summary, E-F-G was the only soluble subcomplex detected in these studies, but assembly was not optimal when these subunits were expressed either first or last. Co-expression of subunits B and CD with E-F-G provided a higher level of assembly, indicating that integrated assembly of N- and Q-modules provides a more efficient pathway.

Rotational Mechanism of FO Motor in the F-Type ATP Synthase Driven by the Proton Motive Force

In FOF1 ATP synthase, driven by the proton motive force across the membrane, the FO motor rotates the central rotor and induces conformational changes in the F1 motor, resulting in ATP synthesis. Recently, many near-atomic resolution structural models have been obtained using cryo-electron microscopy. Despite high resolution, however, static information alone cannot elucidate how and where the protons pass through the FO and how proton passage is coupled to FO rotation. Here, we review theoretical and computational studies based on FO structure models. All-atom molecular dynamics (MD) simulations elucidated changes in the protonation/deprotonation of glutamate—the protein-carrier residue—during rotation and revealed the protonation states that form the “water wire” required for long-range proton hopping. Coarse-grained MD simulations unveiled a free energy surface based on the protonation state and rotational angle of the rotor. Hybrid Monte Carlo and MD simulations showed how proton transfer is coupled to rotation.

Comparison of proton transfer paths to the QA and QB sites of the Rb. sphaeroides photosynthetic reaction centers

The photosynthetic bacterial reaction centers from purple non-sulfur bacteria use light energy to drive the transfer of electrons from cytochrome c to ubiquinone. Ubiquinone bound in the Q_A site cycles between quinone, Q_A, and anionic semiquinone, Q_A^·-, being reduced once and never binding protons. In the Q_B site, ubiquinone is reduced twice by Q_A^·-, binds two protons and is released into the membrane as the quinol, QH₂. The network of hydrogen bonds formed in a molecular dynamics trajectory was drawn to investigate proton transfer pathways from the cytoplasm to each quinone binding site. Q_A is isolated with no path for protons to enter from the surface. In contrast, there is a complex and tangled network requiring residues and waters that can bring protons to Q_B. There are three entries from clusters of surface residues centered around HisH126, GluH224, and HisH68. The network is in good agreement with earlier studies, Mutation of key nodes in the network, such as SerL223, were previously shown to slow proton delivery. Mutational studies had also shown that double mutations of residues such as AspM17 and AspL210 along multiple paths in the network presented here slow the reaction, while single mutations do not. Likewise, mutation of both HisH126 and HisH128, which are at the entry to two paths reduce the rate of proton uptake.

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.

Poor Person's pH Simulation of Membrane Proteins

pH conditions are central to the functioning of all biomolecules. However, implications of pH changes are nontrivial on a molecular scale. Though a rigorous microscopic definition of pH exists, its implementation in classical molecular dynamics (MD) simulations is cumbersome, and more so in large integral membrane systems. In this chapter, an integrative pipeline is described that combines Multi-Conformation Continuum Electrostatics (MCCE) computations with MD simulations to capture the effect of transient protonation states on the coupled conformational changes in transmembrane proteins. The core methodologies are explained, and all the software required to set up this pipeline are outlined with their key parameters. All associated analyses of structure and function are provided using two case studies, namely those of bioenergetic complexes: NADH dehydrogenase (complex I) and V_o domain of V-type ATPase. The hybrid MCCE-MD pipeline has allowed the discovery of hydrogen bond networks, ligand binding pathways, and disease-causing mutations.

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.

Functional Water Wires Catalyze Long-Range Proton Pumping in the Mammalian Respiratory Complex I

The respiratory complex I is a gigantic (1 MDa) redox-driven proton pump that reduces the ubiquinone pool and generates proton motive force to power ATP synthesis in mitochondria. Despite resolved molecular structures and biochemical characterization of the enzyme from multiple organisms, its long-range (∼300 Å) proton-coupled electron transfer (PCET) mechanism remains unsolved. We employ here microsecond molecular dynamics simulations to probe the dynamics of the mammalian complex I in combination with hybrid quantum/classical (QM/MM) free energy calculations to explore how proton pumping reactions are triggered within its 200 Å wide membrane domain. Our simulations predict extensive hydration dynamics of the antiporter-like subunits in complex I that enable lateral proton transfer reactions on a microsecond time scale. We further show how the coupling between conserved ion pairs and charged residues modulate the proton transfer dynamics, and how transmembrane helices and gating residues control the hydration process. Our findings suggest that the mammalian complex I pumps protons by tightly linked conformational and electrostatic coupling principles.

Key role of quinone in the mechanism of respiratory complex I

Complex I is the first and the largest enzyme of respiratory chains in bacteria and mitochondria. The mechanism which couples spatially separated transfer of electrons to proton translocation in complex I is not known. Here we report five crystal structures of T. thermophilus enzyme in complex with NADH or quinone-like compounds. We also determined cryo-EM structures of major and minor native states of the complex, differing in the position of the peripheral arm. Crystal structures show that binding of quinone-like compounds (but not of NADH) leads to a related global conformational change, accompanied by local re-arrangements propagating from the quinone site to the nearest proton channel. Normal mode and molecular dynamics analyses indicate that these are likely to represent the first steps in the proton translocation mechanism. Our results suggest that quinone binding and chemistry play a key role in the coupling mechanism of complex I.

Complex I (NADH:ubiquinone oxidoreductase) is the first enzyme of the respiratory chain in bacteria and mitochondria. Here, the authors present cryo-EM and crystal structures of T. thermophilus complex I in different conformational states and further analyse them by Normal Mode Analysis and molecular dynamics simulations and conclude that quinone redox reactions are important for the coupling mechanism of complex I.

Coupling of quinone dynamics to proton pumping in respiratory complex I

Respiratory complex I (NADH:quinone oxidoreductase) plays a central role in generating the proton electrochemical gradient in mitochondrial and bacterial membranes, which is needed to generate ATP. Several high-resolution structures of complex I have been determined, revealing its intricate architecture and complementing the biochemical and biophysical studies. However, the molecular mechanism of long-range coupling between ubiquinone (Q) reduction and proton pumping is not known. Computer simulations have been applied to decipher the dynamics of Q molecule in the ~30 Å long Q tunnel. In this short report, we discuss the binding and dynamics of Q at computationally predicted Q binding sites, many of which are supported by structural data on complex I. We suggest that the binding of Q at these sites is coupled to proton pumping by means of conformational rearrangements in the conserved loops of core subunits.

Structure and mechanism of the Mrp complex, an ancient cation/proton antiporter

Multiple resistance and pH adaptation (Mrp) antiporters are multi-subunit Na⁺ (or K⁺)/H⁺ exchangers representing an ancestor of many essential redox-driven proton pumps, such as respiratory complex I. The mechanism of coupling between ion or electron transfer and proton translocation in this large protein family is unknown. Here, we present the structure of the Mrp complex from Anoxybacillus flavithermus solved by cryo-EM at 3.0 Å resolution. It is a dimer of seven-subunit protomers with 50 trans-membrane helices each. Surface charge distribution within each monomer is remarkably asymmetric, revealing probable proton and sodium translocation pathways. On the basis of the structure we propose a mechanism where the coupling between sodium and proton translocation is facilitated by a series of electrostatic interactions between a cation and key charged residues. This mechanism is likely to be applicable to the entire family of redox proton pumps, where electron transfer to substrates replaces cation movements.