Molecular simulation and modeling of complex I

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.

From the 'black box' to 'domino effect' mechanism: what have we learned from the structures of respiratory complex I

My group and myself have studied respiratory complex I for almost 30 years, starting in 1994 when it was known as a L-shaped giant 'black box' of bioenergetics. First breakthrough was the X-ray structure of the peripheral arm, followed by structures of the membrane arm and finally the entire complex from Thermus thermophilus. The developments in cryo-EM technology allowed us to solve the first complete structure of the twice larger, ∼1 MDa mammalian enzyme in 2016. However, the mechanism coupling, over large distances, the transfer of two electrons to pumping of four protons across the membrane remained an enigma. Recently we have solved high-resolution structures of mammalian and bacterial complex I under a range of redox conditions, including catalytic turnover. This allowed us to propose a robust and universal mechanism for complex I and related protein families. Redox reactions initially drive conformational changes around the quinone cavity and a long-distance transfer of substrate protons. These set up a stage for a series of electrostatically driven proton transfers along the membrane arm ('domino effect'), eventually resulting in proton expulsion from the distal antiporter-like subunit. The mechanism radically differs from previous suggestions, however, it naturally explains all the unusual structural features of complex I. In this review I discuss the state of knowledge on complex I, including the current most controversial issues.

On the role of ubiquinone in the proton translocation mechanism of respiratory complex I

Complex I converts oxidoreduction energy into a proton electrochemical gradient across the inner mitochondrial or bacterial cell membrane. This gradient is the primary source of energy for aerobic synthesis of ATP. Oxidation of reduced nicotinamide adenine dinucleotide (NADH) by ubiquinone (Q) yields NAD⁺ and ubiquinol (QH₂ ), which is tightly coupled to translocation of four protons from the negatively to the positively charged side of the membrane. Electrons from NADH oxidation reach the iron-sulfur centre N2 positioned near the bottom of a tunnel that extends circa 30 Å from the membrane domain into the hydrophilic domain of the complex. The tunnel is occupied by ubiquinone, which can take a distal position near the N2 centre or proximal positions closer to the membrane. Here, we review important structural, kinetic and thermodynamic properties of ubiquinone that define its role in complex I function. We suggest that this function exceeds that of a mere substrate or electron acceptor and propose that ubiquinone may be the redox element of complex I coupling electron transfer to proton translocation.

High-resolution structure and dynamics of mitochondrial complex I-Insights into the proton pumping mechanism

Mitochondrial NADH:ubiquinone oxidoreductase (complex I) is a 1-MDa membrane protein complex with a central role in energy metabolism. Redox-driven proton translocation by complex I contributes substantially to the proton motive force that drives ATP synthase. Several structures of complex I from bacteria and mitochondria have been determined, but its catalytic mechanism has remained controversial. We here present the cryo-EM structure of complex I from Yarrowia lipolytica at 2.1-Å resolution, which reveals the positions of more than 1600 protein-bound water molecules, of which ~100 are located in putative proton translocation pathways. Another structure of the same complex under steady-state activity conditions at 3.4-Å resolution indicates conformational transitions that we associate with proton injection into the central hydrophilic axis. By combining high-resolution structural data with site-directed mutagenesis and large-scale molecular dynamic simulations, we define details of the proton translocation pathways and offer insights into the redox-coupled proton pumping mechanism of complex I.

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.

Ubiquinone Binding and Reduction by Complex I-Open Questions and Mechanistic Implications

NADH: ubiquinone oxidoreductase (complex I) is the first enzyme complex of the respiratory chain. Complex I is a redox-driven proton pump that contributes to the proton motive force that drives ATP synthase. The structure of complex I has been analyzed by x-ray crystallography and electron cryo-microscopy and is now well-described. The ubiquinone (Q) reduction site of complex I is buried in the peripheral arm and a tunnel-like structure is thought to provide access for the hydrophobic substrate from the membrane. Several intermediate binding positions for Q in the tunnel were identified in molecular simulations. Structural data showed the binding of native Q molecules and short chain analogs and inhibitors in the access pathway and in the Q reduction site, respectively. We here review the current knowledge on the interaction of complex I with Q and discuss recent hypothetical models for the coupling mechanism.

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.

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.

Respiratory complex I - Mechanistic insights and advances in structure determination

Complex I is the largest and most intricate redox-driven proton pump of the respiratory chain. The structure of bacterial and mitochondrial complex I has been determined by X-ray crystallography and cryo-EM at increasing resolution. The recent cryo-EM structures of the complex I-like NDH complex and membrane bound hydrogenase open a new and more comprehensive perspective on the complex I superfamily. Functional studies and molecular modeling approaches have greatly advanced our understanding of the catalytic cycle of complex I. However, the molecular mechanism by which energy is extracted from the redox reaction and utilized to drive proton translocation is unresolved and a matter of ongoing debate. Here, we review progress in structure determination and functional characterization of complex I and discuss current mechanistic models.

Long-range proton-coupled electron transfer in biological energy conversion: towards mechanistic understanding of respiratory complex I

Biological energy conversion is driven by efficient enzymes that capture, store and transfer protons and electrons across large distances. Recent advances in structural biology have provided atomic-scale blueprints of these types of remarkable molecular machinery, which together with biochemical, biophysical and computational experiments allow us to derive detailed energy transduction mechanisms for the first time. Here, I present one of the most intricate and least understood types of biological energy conversion machinery, the respiratory complex I, and how its redox-driven proton-pump catalyses charge transfer across approximately 300 Å distances. After discussing the functional elements of complex I, a putative mechanistic model for its action-at-a-distance effect is presented, and functional parallels are drawn to other redox- and light-driven ion pumps.

Single-Molecule Imaging and Computational Microscopy Approaches Clarify the Mechanism of the Dimerization and Membrane Interactions of Green Fluorescent Protein

Green fluorescent protein (GFP) is widely used as a biomarker in living systems; however, GFP and its variants are prone to forming low-affinity dimers under physiological conditions. This undesirable tendency is exacerbated when fluorescent proteins (FP) are confined to membranes, fused to naturally-oligomeric proteins, or expressed at high levels in cells. Oligomerization of FPs introduces artifacts into the measurement of subunit stoichiometry, as well as interactions between proteins fused to FPs. Introduction of a single mutation, A206K, has been shown to disrupt hydrophobic interactions in the region responsible for GFP dimerization, thereby contributing to its monomerization. Nevertheless, a detailed understanding of how this single amino acid-dependent inhibition of dimerization in GFP occurs at the atomic level is still lacking. Single-molecule experiments combined with computational microscopy (atomistic molecular dynamics) revealed that the amino group of A206 contributes to GFP dimer formation via a multivalent electrostatic interaction. We further showed that myristoyl modification is an efficient mechanism to promote membrane attachment of GFP. Molecular dynamics-based site-directed mutagenesis has been used to identify the key functional residues in FPs. The data presented here have been utilized as a monomeric control in downstream single-molecule studies, facilitating more accurate stoichiometry quantification of functional protein complexes in living cells.

Insights into the respiratory chain and oxidative stress

Reactive oxygen species (ROS) are highly reactive reduced oxygen molecules that result from aerobic metabolism. The common forms are the superoxide anion (O2∙-) and hydrogen peroxide (H2O2) and their derived forms, hydroxyl radical (HO∙) and hydroperoxyl radical (HOO∙). Their production sites in mitochondria are reviewed. Even though being highly toxic products, ROS seem important in transducing information from dysfunctional mitochondria. Evidences of signal transduction mediated by ROS in mitochondrial deficiency contexts are then presented in different organisms such as yeast, mammals or photosynthetic organisms.

Redox-coupled quinone dynamics in the respiratory complex I

Complex I couples the free energy released from quinone (Q) reduction to pump protons across the biological membrane in the respiratory chains of mitochondria and many bacteria. The Q reduction site is separated by a large distance from the proton-pumping membrane domain. To address the molecular mechanism of this long-range proton-electron coupling, we perform here full atomistic molecular dynamics simulations, free energy calculations, and continuum electrostatics calculations on complex I from Thermus thermophilus We show that the dynamics of Q is redox-state-dependent, and that quinol, QH₂, moves out of its reduction site and into a site in the Q tunnel that is occupied by a Q analog in a crystal structure of Yarrowia lipolytica We also identify a second Q-binding site near the opening of the Q tunnel in the membrane domain, where the Q headgroup forms strong interactions with a cluster of aromatic and charged residues, while the Q tail resides in the lipid membrane. We estimate the effective diffusion coefficient of Q in the tunnel, and in turn the characteristic time for Q to reach the active site and for QH₂ to escape to the membrane. Our simulations show that Q moves along the Q tunnel in a redox-state-dependent manner, with distinct binding sites formed by conserved residue clusters. The motion of Q to these binding sites is proposed to be coupled to the proton-pumping machinery in complex I.

Role of water and protein dynamics in proton pumping by respiratory complex I

Membrane bound respiratory complex I is the key enzyme in the respiratory chains of bacteria and mitochondria, and couples the reduction of quinone to the pumping of protons across the membrane. Recently solved crystal or electron microscopy structures of bacterial and mitochondrial complexes have provided significant insights into the electron and proton transfer pathways. However, due to large spatial separation between the electron and proton transfer routes, the molecular mechanism of coupling remains unclear. Here, based on atomistic molecular dynamics simulations performed on the entire structure of complex I from Thermus thermophilus, we studied the hydration of the quinone-binding site and the membrane-bound subunits. The data from simulations show rapid diffusion of water molecules in the protein interior, and formation of hydrated regions in the three antiporter-type subunits. An unexpected water-protein based connectivity between the middle of the Q-tunnel and the fourth proton channel is also observed. The protonation-state dependent dynamics of key acidic residues in the Nqo8 subunit suggest that the latter may be linked to redox-coupled proton pumping in complex I. We propose that in complex I the proton and electron transfer paths are not entirely separate, instead the nature of coupling may in part be ‘direct’.