Electrostatics, hydration, and proton transfer dynamics in the membrane domain of respiratory complex I

Coupling between electrons' spin and proton transfer in chiral biological crystals

Proton transport plays a fundamental role in many biological and chemical systems. In life, proton transport is crucial for biochemical and physiological functions. It is usually accepted that the main mechanism of proton transfer is a result of hopping between neighboring water molecules and amino acid side chains. It was recently suggested that the proton transfer can be simultaneously coupled with electron transfer. As life is homochiral, proton transfer in biology is occurring in a chiral environment. In this environment, the chiral-induced spin selectivity effect relating to electron transfer and chirality is expected to occur. The present work establishes that the proton transfer is coupled to a specific electron spin polarization in lysozyme crystals, associating proton transfer to electron movement and polarization. To preserve total angular momentum, this motion may be coupled to chiral phonons that propagate in the crystal. Our work shows that the interaction of the electrons' spin and phonons is very significant in proton transfer through lysosome crystals. Injecting the opposite electron spin into the lysosome crystal results in a significant change in proton transfer impedance. This study presents the support for the proton-coupled electron transfer mechanism and indicates the importance of spin polarization in the process.

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.

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.

Rotary mechanism of the prokaryotic Vo motor driven by proton motive force

ATP synthases play a crucial role in energy production by utilizing the proton motive force (pmf) across the membrane to rotate their membrane-embedded rotor c-ring, and thus driving ATP synthesis in the hydrophilic catalytic hexamer. However, the mechanism of how pmf converts into c-ring rotation remains unclear. This study presents a 2.8 Å cryo-EM structure of the Vo domain of V/A-ATPase from Thermus thermophilus, revealing precise orientations of glutamate (Glu) residues in the c12-ring. Three Glu residues face a water channel, with one forming a salt bridge with the Arginine in the stator (a/Arg). Molecular dynamics (MD) simulations show that protonation of specific Glu residues triggers unidirectional Brownian motion of the c12-ring towards ATP synthesis. When the key Glu remains unprotonated, the salt bridge persists, blocking rotation. These findings suggest that asymmetry in the protonation of c/Glu residues biases c12-ring movement, facilitating rotation and ATP synthesis.

Dissected antiporter modules establish minimal proton-conduction elements of the respiratory complex I

The respiratory Complex I is a highly intricate redox-driven proton pump that powers oxidative phosphorylation across all domains of life. Yet, despite major efforts in recent decades, its long-range energy transduction principles remain highly debated. We create here minimal proton-conducting membrane modules by engineering and dissecting the key elements of the bacterial Complex I. By combining biophysical, biochemical, and computational experiments, we show that the isolated antiporter-like modules of Complex I comprise all functional elements required for conducting protons across proteoliposome membranes. We find that the rate of proton conduction is controlled by conformational changes of buried ion-pairs that modulate the reaction barriers by electric field effects. The proton conduction is also modulated by bulky residues along the proton channels that are key for establishing a tightly coupled proton pumping machinery in Complex I. Our findings provide direct experimental evidence that the individual antiporter modules are responsible for the proton transport activity of Complex I. On a general level, our findings highlight electrostatic and conformational coupling mechanisms in the modular energy-transduction machinery of Complex I with distinct similarities to other enzymes.

A computational study to assess the pathogenicity of single or combinations of missense variants on respiratory complex I

Variants found in the respiratory complex I (CI) subunit genes encoded by mitochondrial DNA can cause severe genetic diseases. However, it is difficult to establish a priori whether a single or a combination of CI variants may impact oxidative phosphorylation. Here we propose a computational approach based on coarse-grained molecular dynamics simulations aimed at investigating new CI variants. One of the primary CI variants associated with the Leber hereditary optic neuropathy (m.14484T>C/MT-ND6) was used as a test case and was investigated alone or in combination with two additional rare CI variants whose role remains uncertain. We found that the primary variant positioned in the E-channel region, which is fundamental for CI function, stiffens the enzyme dynamics. Moreover, a new mechanism for the transition between π- and α-conformation in the helix carrying the primary variant is proposed. This may have implications for the E-channel opening/closing mechanism. Finally, our findings show that one of the rare variants, located next to the primary one, further worsens the stiffening, while the other rare variant does not affect CI function. This approach may be extended to other variants candidate to exert a pathogenic impact on CI dynamics, or to investigate the interaction of multiple variants.

Perspective on Integrative Simulations of Bioenergetic Domains

Bioenergetic processes in cells, such as photosynthesis or respiration, integrate many time and length scales, which makes the simulation of energy conversion with a mere single level of theory impossible. Just like the myriad of experimental techniques required to examine each level of organization, an array of overlapping computational techniques is necessary to model energy conversion. Here, a perspective is presented on recent efforts for modeling bioenergetic phenomena with a focus on molecular dynamics simulations and its variants as a primary method. An overview of the various classical, quantum mechanical, enhanced sampling, coarse-grained, Brownian dynamics, and Monte Carlo methods is presented. Example applications discussed include multiscale simulations of membrane-wide electron transport, rate kinetics of ATP turnover from electrochemical gradients, and finally, integrative modeling of the chromatophore, a photosynthetic pseudo-organelle.

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.

Mutation at the entrance of the quinone cavity severely disrupts quinone binding in respiratory complex I

In all resolved structures of complex I, there exists a tunnel-like Q-chamber for ubiquinone binding and reduction. The entrance to the Q-chamber in ND1 subunit forms a narrow bottleneck, which is rather tight and requires thermal conformational changes for ubiquinone to get in and out of the binding chamber. The substitution of alanine with threonine at the bottleneck (AlaThr MUT), associated with 3460/ND1 mtDNA mutation in human complex I, is implicated in Leber's Hereditary Optic Neuropathy (LHON). Here, we show the AlaThr MUT further narrows the Q-chamber entrance cross-section area by almost 30%, increasing the activation free energy barrier of quinone passage by approximately 5 kJ mol-1. This severely disrupts quinone binding and reduction as quinone passage through the bottleneck is slowed down almost tenfold. Our estimate of the increase in free energy barrier is entirely due to the bottleneck narrowing, leading to a reduction of the transition state entropy between WT and MUT, and thus more difficult quinone passage. Additionally, we investigate details of possible water exchange between the Q-chamber and membrane. We find water exchange is dynamic in WT but may be severely slowed in MUT. We propose that LHON symptoms caused by 3460/ND1 mtDNA mutation are due to slowed quinone binding. This leads to an increased production of reactive oxidative species due to upstream electron backup at the FMN site of complex I, thus resulting in a mt bioenergetic defect.

Quinone Catalysis Modulates Proton Transfer Reactions in the Membrane Domain of Respiratory Complex I

Complex I is a redox-driven proton pump that drives electron transport chains and powers oxidative phosphorylation across all domains of life. Yet, despite recently resolved structures from multiple organisms, it still remains unclear how the redox reactions in Complex I trigger proton pumping up to 200 Å away from the active site. Here, we show that the proton-coupled electron transfer reactions during quinone reduction drive long-range conformational changes of conserved loops and trans-membrane (TM) helices in the membrane domain of Complex I from Yarrowia lipolytica. We find that the conformational switching triggers a π → α transition in a TM helix (TM3ND6) and establishes a proton pathway between the quinone chamber and the antiporter-like subunits, responsible for proton pumping. Our large-scale (>20 μs) atomistic molecular dynamics (MD) simulations in combination with quantum/classical (QM/MM) free energy calculations show that the helix transition controls the barrier for proton transfer reactions by wetting transitions and electrostatic effects. The conformational switching is enabled by re-arrangements of ion pairs that propagate from the quinone binding site to the membrane domain via an extended network of conserved residues. We find that these redox-driven changes create a conserved coupling network within the Complex I superfamily, with point mutations leading to drastic activity changes and mitochondrial disorders. On a general level, our findings illustrate how catalysis controls large-scale protein conformational changes and enables ion transport across biological membranes.

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.

Single protonation of the reduced quinone in respiratory complex I drives four-proton pumping

Complex I is a key proton-pumping enzyme in bacterial and mitochondrial respiratory electron transport chains. Using quantum chemistry and electrostatic calculations, we have examined the pKa of the reduced quinone QH-/QH2 in the catalytic cavity of complex I. We find that pKa (QH-/QH2) is very high, above 20. This means that the energy of a single protonation reaction of the doubly reduced quinone (i.e. the reduced semiquinone QH-) is sufficient to drive four protons across the membrane with a potential of 180 mV. Based on these calculations, we propose a possible scheme of redox-linked proton pumping by complex I. The model explains how the energy of the protonation reaction can be divided equally among four pumping units of the pump, and how a single proton can drive translocation of four additional protons in multiple pumping blocks.

Mechanical Allosteric Couplings of Redox-Induced Conformational Changes in Respiratory Complex I

We apply linear response theory to calculate mechanical allosteric couplings in respiratory complex I between the iron sulfur cluster N2, located in the catalytic cavity, and the membrane part of the enzyme, separated from it by more than 50 Å. According to our hypothesis, the redox reaction of ubiquinone in the catalytic cavity of the enzyme generates an unbalanced charge that via repulsion of the charged redox center N2 produces local mechanical stress that transmits into the membrane part of the enzyme where it induces proton pumping. Using coarse-grained simulations of the enzyme, we calculated mechanistic allosteric couplings that reveal the pathways of the mechanical transmission of the stress along the enzyme. The results shed light on the recent experimental studies where a stabilization of the enzyme with an introduced disulfide bridge resulted in the abolishing of proton pumping. Simulation of the disulfide bond action indicates a dramatic change of the mechanistic coupling pathways in line with experimental findings.

Human clinical mutations in mitochondrially encoded subunits of Complex I can be successfully modeled in E. coli

Respiratory Complex I is the site of a large fraction of the mutations that appear to cause mitochondrial disease. Seven of its subunits are mitochondrially encoded, and therefore, such mutants are particularly difficult to construct in cell-culture model systems. We have selected 13 human clinical mutations found in ND2, ND3, ND4, ND4L, ND5 and ND6 that are generally found at subunit interfaces, and not in critical residues. These mutations have been modeled in E. coli subunits of Complex I, nuoN, nuoA, nuoM, nuoK, nuoL, and nuoJ, respectively. All mutants were expressed from a plasmid encoding the entire nuo operon, and membrane vesicles were analyzed for deamino-NADH oxidase activity, and proton translocation activity. ND5 mutants were also analyzed using a time-delayed expression system, recently described by this lab. Other mutants were analyzed for the ability to associate in subcomplexes, after expression of subsets of the genes. For most mutants there was a positive correlation between those that were previously determined to be pathogenic, or likely to be pathogenic, and those that we found with compromised Complex I activity or subunit interactions in E. coli. In conclusion, this approach provides another way to explore the deleterious effects of human mitochondrial mutations, and it can contribute to molecular understanding of such mutations.

Resolving Chemical Dynamics in Biological Energy Conversion: Long-Range Proton-Coupled Electron Transfer in Respiratory Complex I

Biological energy conversion is catalyzed by membrane-bound proteins that transduce chemical or light energy into energy forms that power endergonic processes in the cell. At a molecular level, these catalytic processes involve elementary electron-, proton-, charge-, and energy-transfer reactions that take place in the intricate molecular machineries of cell respiration and photosynthesis. Recent developments in structural biology, particularly cryo-electron microscopy (cryoEM), have resolved the molecular architecture of several energy transducing proteins, but detailed mechanistic principles of their charge transfer reactions still remain poorly understood and a major challenge for modern biochemical research. To this end, multiscale molecular simulations provide a powerful approach to probe mechanistic principles on a broad range of time scales (femtoseconds to milliseconds) and spatial resolutions (10¹–10⁶ atoms), although technical challenges also require balancing between the computational accuracy, cost, and approximations introduced within the model. Here we discuss how the combination of atomistic (aMD) and hybrid quantum/classical molecular dynamics (QM/MM MD) simulations with free energy (FE) sampling methods can be used to probe mechanistic principles of enzymes responsible for biological energy conversion. We present mechanistic explorations of long-range proton-coupled electron transfer (PCET) dynamics in the highly intricate respiratory chain enzyme Complex I, which functions as a redox-driven proton pump in bacterial and mitochondrial respiratory chains by catalyzing a 300 Å fully reversible PCET process. This process is initiated by a hydride (H^–) transfer between NADH and FMN, followed by long-range (>100 Å) electron transfer along a wire of 8 FeS centers leading to a quinone biding site. The reduction of the quinone to quinol initiates dissociation of the latter to a second membrane-bound binding site, and triggers proton pumping across the membrane domain of complex I, in subunits up to 200 Å away from the active site. Our simulations across different size and time scales suggest that transient charge transfer reactions lead to changes in the internal hydration state of key regions, local electric fields, and the conformation of conserved ion pairs, which in turn modulate the dynamics of functional steps along the reaction cycle. Similar functional principles, which operate on much shorter length scales, are also found in some unrelated proteins, suggesting that enzymes may employ conserved principles in the catalysis of biological energy transduction processes.

Functional Dynamics of an Ancient Membrane-Bound Hydrogenase

The membrane-bound hydrogenase (Mbh) is a redox-driven Na⁺/H⁺ transporter that employs the energy from hydrogen gas (H₂) production to catalyze proton pumping and Na⁺/H⁺ exchange across cytoplasmic membranes of archaea. Despite a recently resolved structure of this ancient energy-transducing enzyme [Yu et al. Cell2018, 173, 1636–1649], the molecular principles of its redox-driven ion-transport mechanism remain puzzling and of major interest for understanding bioenergetic principles of early cells. Here we use atomistic molecular dynamics (MD) simulations in combination with data clustering methods and quantum chemical calculations to probe principles underlying proton reduction as well as proton and sodium transport in Mbh from the hyperthermophilic archaeon Pyrococcus furiosus. We identify putative Na⁺ binding sites and proton pathways leading across the membrane and to the NiFe-active center as well as conformational changes that regulate ion uptake. We suggest that Na⁺ binding and protonation changes at a putative ion-binding site couple to proton transfer across the antiporter-like MbhH subunit by modulating the conformational state of a conserved ion pair at the subunit interface. Our findings illustrate conserved coupling principles within the complex I superfamily and provide functional insight into archaeal energy transduction mechanisms.

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.

Reverse Protonation of Buried Ion-Pairs in Staphylococcal Nuclease Mutants

Although buried titratable residues in protein cavities are often of major functional importance, it is generally challenging to understand their properties such as the ionization state and factors of stabilization based on experimental studies alone. A specific set of examples involve buried Glu-Lys pairs in a series of variants of Staphylococcal nuclease, for which recent structural and thermodynamic studies appeared to suggest that both the stability and the ionization state of the buried Glu-Lys pair are sensitive to its orientation (i.e., Glu23-Lys36 vs Lys23-Glu36). To further clarify the situation, especially ionization states of the buried Glu-Lys pairs, we have conducted extensive molecular dynamics simulations and free energy computations. Microsecond molecular dynamics simulations show that the hydration level of the cavity depends on the orientation of the buried ion-pair therein as well as its ionization state; free energy simulations recapitulate the relative stability of Glu23-Lys36 (EK) vs Lys23-Glu36 (KE) mutants measured experimentally, although the difference is similar in magnitude regardless of the ionization state of the Glu-Lys pair. A complementary set of free energy simulations strongly suggests that, in contrast to the original suggestion in the experimental analysis, the Glu and Lys residues prefer to adopt their charge-neutral rather than the ionized states. This result is consistent with the low dielectric constant computed for water in the cavity, which makes it difficult for the protein cavity to stabilize a pair of charged Glu-Lys residues, even with water penetration. The current study highlights the role of free energy simulations in understanding the ionization state of buried titratable residues and the relevant energetic contributions, forming the basis for the rational design of buried charge networks in proteins.

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.

Evolution of complex I-like respiratory complexes

The modern-day respiratory complex I shares a common ancestor with the membrane-bound hydrogenase (MBH) and membrane-bound sulfane sulfur reductase (MBS). MBH and MBS use protons and sulfur as their respective electron sinks, which helped to conserve energy during early life in the Proterozoic era when the Earth's atmosphere was low in oxygen. MBH and MBS likely evolved from an integration of an ancestral, membrane-embedded, multiple resistance and pH antiporter and a soluble redox-active module encompassing a [NiFe] hydrogenase. In this review, we discuss how the structures of MBH, MBS, multiple resistance and pH, photosynthetic NADH dehydrogenase-like complex type-1, and complex I, which have been determined recently, thanks to the advent of high-resolution cryo-EM, have significantly improved our understanding of the catalytic reaction mechanisms and the evolutionary relationships of the respiratory complexes.

Mitochondrial Structure and Bioenergetics in Normal and Disease Conditions

Mitochondria are ubiquitous intracellular organelles found in almost all eukaryotes and involved in various aspects of cellular life, with a primary role in energy production. The interest in this organelle has grown stronger with the discovery of their link to various pathologies, including cancer, aging and neurodegenerative diseases. Indeed, dysfunctional mitochondria cannot provide the required energy to tissues with a high-energy demand, such as heart, brain and muscles, leading to a large spectrum of clinical phenotypes. Mitochondrial defects are at the origin of a group of clinically heterogeneous pathologies, called mitochondrial diseases, with an incidence of 1 in 5000 live births. Primary mitochondrial diseases are associated with genetic mutations both in nuclear and mitochondrial DNA (mtDNA), affecting genes involved in every aspect of the organelle function. As a consequence, it is difficult to find a common cause for mitochondrial diseases and, subsequently, to offer a precise clinical definition of the pathology. Moreover, the complexity of this condition makes it challenging to identify possible therapies or drug targets.

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.

Molecular Dynamics Simulation of T10609C and C10676G Mutations of Mitochondrial ND4L Gene Associated With Proton Translocation in Type 2 Diabetes Mellitus and Cataract Patients

The mutation rate of mitochondrial DNA (mtDNA) is 17 times higher than nuclear DNA, and these mutations can cause mitochondrial disease in 1 of 10.000 people. The T10609C mutation was identified in type 2 diabetes mellitus (T2DM) patients and the C10676G mutation in cataract patients, with both mutations occurring in the ND4L gene of mtDNA that encodes ND4L protein. ND4L protein, a subunit of complex I in the respiratory complex, has been shown to play a role in the proton translocation process. The purpose of this study was to investigate the effect of both mutations on the proton translocation mechanism. Mutation mapping showed changes in amino acids M47T (T10609C) and C69W (C10676G). The 100 ns molecular dynamics (MD) simulations performed on native and mutants of ND4L-ND6 subunits. It is revealed that the native model had a similar proton translocation pathway to that of complex I from other organisms. Interestingly, the mutant M47T and C69W showed the interruption of the translocation pathway by a hydrogen bond formation between Glu34 and Tyr157. It is observed that the mutations were restricting the passage of water molecules through the transmembrane region. These results could help to develop the computational assay for the validation of a specific genetic biomarker for T2DM and cataracts.

The coupling mechanism of mammalian respiratory complex I

Mitochondrial complex I couples NADH:ubiquinone oxidoreduction to proton pumping by an unknown mechanism. Here, we present cryo-electron microscopy structures of ovine complex I in five different conditions, including turnover, at resolutions up to 2.3 to 2.5 angstroms. Resolved water molecules allowed us to experimentally define the proton translocation pathways. Quinone binds at three positions along the quinone cavity, as does the inhibitor rotenone that also binds within subunit ND4. Dramatic conformational changes around the quinone cavity couple the redox reaction to proton translocation during open-to-closed state transitions of the enzyme. In the induced deactive state, the open conformation is arrested by the ND6 subunit. We propose a detailed molecular coupling mechanism of complex I, which is an unexpected combination of conformational changes and electrostatic interactions.

Water-Gated Proton Transfer Dynamics in Respiratory Complex I

The respiratory complex I transduces redox energy into an electrochemical proton gradient in aerobic respiratory chains, powering energy-requiring processes in the cell. However, despite recently resolved molecular structures, the mechanism of this gigantic enzyme remains poorly understood. By combining large-scale quantum and classical simulations with site-directed mutagenesis and biophysical experiments, we show here how the conformational state of buried ion-pairs and water molecules control the protonation dynamics in the membrane domain of complex I and establish evolutionary conserved long-range coupling elements. We suggest that an electrostatic wave propagates in forward and reverse directions across the 200 Å long membrane domain during enzyme turnover, without significant dissipation of energy. Our findings demonstrate molecular principles that enable efficient long-range proton–electron coupling (PCET) and how perturbation of this PCET machinery may lead to development of mitochondrial disease.

Scalable molecular dynamics on CPU and GPU architectures with NAMD

NAMDis a molecular dynamics program designed for high-performance simulations of very large biological objects on CPU- and GPU-based architectures. NAMD offers scalable performance on petascale parallel supercomputers consisting of hundreds of thousands of cores, as well as on inexpensive commodity clusters commonly found in academic environments. It is written in C++ and leans on Charm++ parallel objects for optimal performance on low-latency architectures. NAMD is a versatile, multipurpose code that gathers state-of-the-art algorithms to carry out simulations in apt thermodynamic ensembles, using the widely popular CHARMM, AMBER, OPLS, and GROMOS biomolecular force fields. Here, we review the main features of NAMD that allow both equilibrium and enhanced-sampling molecular dynamics simulations with numerical efficiency. We describe the underlying concepts utilized by NAMD and their implementation, most notably for handling long-range electrostatics; controlling the temperature, pressure, and pH; applying external potentials on tailored grids; leveraging massively parallel resources in multiple-copy simulations; and hybrid quantum-mechanical/molecular-mechanical descriptions. We detail the variety of options offered by NAMD for enhanced-sampling simulations aimed at determining free-energy differences of either alchemical or geometrical transformations and outline their applicability to specific problems. Last, we discuss the roadmap for the development of NAMD and our current efforts toward achieving optimal performance on GPU-based architectures, for pushing back the limitations that have prevented biologically realistic billion-atom objects to be fruitfully simulated, and for making large-scale simulations less expensive and easier to set up, run, and analyze. NAMD is distributed free of charge with its source code at www.ks.uiuc.edu.

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.

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

BACKGROUND

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

METHODS

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

RESULTS

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

CONCLUSIONS

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

GENERAL SIGNIFICANCE

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

Proton motive function of the terminal antiporter-like subunit in respiratory complex I

In the aerobic respiratory chains of many organisms, complex I functions as the first electron input. By reducing ubiquinone (Q) to ubiquinol, it catalyzes the translocation of protons across the membrane as far as ~200 Å from the site of redox reactions. Despite significant amount of structural and biochemical data, the details of redox coupled proton pumping in complex I are poorly understood. In particular, the proton transfer pathways are extremely difficult to characterize with the current structural and biochemical techniques. Here, we applied multiscale computational approaches to identify the proton transfer paths in the terminal antiporter-like subunit of complex I. Data from combined classical and quantum chemical simulations reveal for the first time structural elements that are exclusive to the subunit, and enables the enzyme to achieve coupling between the spatially separated Q redox reactions and proton pumping. By studying long time scale protonation and hydration dependent conformational dynamics of key amino acid residues, we provide novel insights into the proton pumping mechanism of complex I.

Respiratory complex I - structure, mechanism and evolution

Respiratory complex I is an intricate multi-subunit membrane protein with a central function in aerobic energy metabolism. During the last years, structures of mitochondrial complex I and respiratory supercomplexes were determined by cryo-EM at increasing resolution. Structural and computational studies have shed light on the dynamics of proton translocation pathways, the interaction of complex I with lipids and the unusual access pathway of ubiquinone to the active site. Recent advances in understanding complex I function include characterization of specific conformational changes that are critical for proton pumping. Cryo-EM structures of the NADH dehydrogenase-like (NDH) complex of photosynthesis and a bacterial membrane bound hydrogenase (MBH) have provided a broader perspective on the complex I superfamily.

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.

Autophosphorylation activates c-Src kinase through global structural rearrangements

The prototypical kinase c-Src plays an important role in numerous signal transduction pathways, where its activity is tightly regulated by two phosphorylation events. Phosphorylation at a specific tyrosine by C-terminal Src kinase inactivates c-Src, whereas autophosphorylation is essential for the c-Src activation process. However, the structural consequences of the autophosphorylation process still remain elusive. Here we investigate how the structural landscape of c-Src is shaped by nucleotide binding and phosphorylation of Tyr⁴¹⁶ using biochemical experiments, hydrogen/deuterium exchange MS, and atomistic molecular simulations. We show that the initial steps of kinase activation involve large rearrangements in domain orientation. The kinase domain is highly dynamic and has strong cross-talk with the regulatory domains, which are displaced by autophosphorylation. Although the regulatory domains become more flexible and detach from the kinase domain because of autophosphorylation, the kinase domain gains rigidity, leading to stabilization of the ATP binding site and a 4-fold increase in enzymatic activity. Our combined results provide a molecular framework of the central steps in c-Src kinase regulation process with possible implications for understanding general kinase activation mechanisms.

Electric field modulated redox-driven protonation and hydration energetics in energy converting enzymes

Biological energy conversion is catalysed by proton-coupled electron transfer (PCET) reactions that form the chemical basis of respiratory and photosynthetic enzymes. Despite recent advances in structural, biophysical, and computational experiments, the mechanistic principles of these reactions still remain elusive. Based on common functional features observed in redox enzymes, we study here generic mechanistic models for water-mediated long-range PCET reactions. We show how a redox reaction within a buried protein environment creates an electric field that induces hydration changes between the proton acceptor and donor groups, and in turn, lowers the reaction barrier and increases the thermodynamic driving forces for the water-mediated PCET process. We predict linear free energy relationships, and discuss the proposed mechanism in context of PCET in cytochrome c oxidase.

Molecular dynamics and structural models of the cyanobacterial NDH-1 complex

NDH-1 is a gigantic redox-driven proton pump linked with respiration and cyclic electron flow in cyanobacterial cells. Based on experimentally resolved X-ray and cryo-EM structures of the respiratory complex I, we derive here molecular models of two isoforms of the cyanobacterial NDH-1 complex involved in redox-driven proton pumping (NDH-1L) and CO2-fixation (NDH-1MS). Our models show distinct structural and dynamic similarities to the core architecture of the bacterial and mammalian respiratory complex I. We identify putative plastoquinone-binding sites that are coupled by an electrostatic wire to the proton pumping elements in the membrane domain of the enzyme. Molecular simulations suggest that the NDH-1L isoform undergoes large-scale hydration changes that support proton-pumping within antiporter-like subunits, whereas the terminal subunit of the NDH-1MS isoform lacks such structural motifs. Our work provides a putative molecular blueprint for the complex I-analogue in the photosynthetic energy transduction machinery and demonstrates that general mechanistic features of the long-range proton-pumping machinery are evolutionary conserved in the complex I-superfamily.

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.

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

Cytochrome c oxidase (CcO) drives aerobic respiratory chains in all organisms by transducing the free energy from oxygen reduction into an electrochemical proton gradient across a biological membrane.

Cytochrome c oxidase (CcO) drives aerobic respiratory chains in all organisms by transducing the free energy from oxygen reduction into an electrochemical proton gradient across a biological membrane. CcO employs the so-called D- and K-channels for proton uptake, but the molecular mechanism for activation of the K-channel has remained elusive for decades. We show here by combining large-scale atomistic molecular simulations with graph-theoretical water network analysis, and hybrid quantum/classical (QM/MM) free energy calculations, that the K-channel is activated by formation of a reactive oxidized intermediate in the binuclear heme a₃/Cu_B active site. This state induces electrostatic, hydration, and conformational changes that lower the barrier for proton transfer along the K-channel by dewetting pathways that connect the D-channel with the active site. Our combined results reconcile previous experimental findings and indicate that water dynamics plays a decisive role in the proton pumping machinery in CcO.

Cryo-EM structures of complex I from mouse heart mitochondria in two biochemically defined states

Complex I (NADH:ubiquinone oxidoreductase) uses the reducing potential of NADH to drive protons across the energy-transducing inner membrane and power oxidative phosphorylation in mammalian mitochondria. Recent cryo-EM analyses have produced near-complete models of all 45 subunits in the bovine, ovine and porcine complexes and have identified two states relevant to complex I in ischemia-reperfusion injury. Here, we describe the 3.3-Å structure of complex I from mouse heart mitochondria, a biomedically relevant model system, in the 'active' state. We reveal a nucleotide bound in subunit NDUFA10, a nucleoside kinase homolog, and define mechanistically critical elements in the mammalian enzyme. By comparisons with a 3.9-Å structure of the 'deactive' state and with known bacterial structures, we identify differences in helical geometry in the membrane domain that occur upon activation or that alter the positions of catalytically important charged residues. Our results demonstrate the capability of cryo-EM analyses to challenge and develop mechanistic models for mammalian complex I.

Structural Stability of Peptidic His-Containing Proton Wire in Solution and in the Adsorbed State

Molecular wires are functional molecules applicable in the field of transfer processes in technological and biochemical applications. Besides molecular wires with the ability to transfer electrons, research is currently focused on molecular wires with high proton affinity and proton transfer ability. Recently, proposed peptidic proton wires (H wires) are one example. Their ability to mediate the transport of protons from aqueous solutions onto the surface of a Hg electrode in a catalytic hydrogen evolution reaction was investigated by constant-current chronopotentiometric stripping. However, elucidating the structure of H wires and rationalizing their stability are key requirements for their further research and application. In this article, we focus on the His (H) and Ala (A)-containing peptidic H wire A₃-(H-A₂)₆ in solution and after its immobilization onto the electrode surface in the presence of the secondary structure stabilizer 2,2,2-trifluoroethanol (TFE). We found that the solvent containing more than 25% of TFE stabilizes the helical structure of A₃-(H-A₂)₆ not only in solution but also in the adsorbed state. The TFE efficacy to stabilize α-helical structure was confirmed using high-resolution nuclear magnetic resonance, circular dichroism, and molecular dynamics simulation. Experimental and theoretical results indicated A₃-(H-A₂)₆ to be a high proton-affinity peptidic H wire with an α-helical structure stabilized by TFE, which was confirmed in a comparative study with hexahistidine as an example of a peptide with a definitely disordered and random coil structure. The results presented here could be used for further investigation of the peptidic H wires and for the application of electrochemical methods in the research of proton transfer phenomena in general.

How inter-subunit contacts in the membrane domain of complex I affect proton transfer energetics

The respiratory complex I is a redox-driven proton pump that employs the free energy released from quinone reduction to pump protons across its complete ca. 200 Å wide membrane domain. Despite recently resolved structures and molecular simulations, the exact mechanism for the proton transport process remains unclear. Here we combine large-scale molecular simulations with quantum chemical density functional theory (DFT) models to study how contacts between neighboring antiporter-like subunits in the membrane domain of complex I affect the proton transfer energetics. Our combined results suggest that opening of conserved Lys/Glu ion pairs within each antiporter-like subunit modulates the barrier for the lateral proton transfer reactions. Our work provides a mechanistic suggestion for key coupling effects in the long-range force propagation process of complex I.

Structure of an Ancient Respiratory System

Hydrogen gas-evolving membrane-bound hydrogenase (MBH) and quinone-reducing complex I are homologous respiratory complexes with a common ancestor, but a structural basis for their evolutionary relationship is lacking. Here, we report the cryo-EM structure of a 14-subunit MBH from the hyperthermophile Pyrococcus furiosus. MBH contains a membrane-anchored hydrogenase module that is highly similar structurally to the quinone-binding Q-module of complex I while its membrane-embedded ion-translocation module can be divided into a H⁺- and a Na⁺-translocating unit. The H⁺-translocating unit is rotated 180° in-membrane with respect to its counterpart in complex I, leading to distinctive architectures for the two respiratory systems despite their largely conserved proton-pumping mechanisms. The Na⁺-translocating unit, absent in complex I, resembles that found in the Mrp H⁺/Na⁺ antiporter and enables hydrogen gas evolution by MBH to establish a Na⁺ gradient for ATP synthesis near 100°C. MBH also provides insights into Mrp structure and evolution of MBH-based respiratory enzymes.

Peroxiredoxin-3 attenuates traumatic neuronal injury through preservation of mitochondrial function

Peroxiredoxins (PRDXs) are a highly conserved family of thiol peroxidases that scavenge peroxides in cells. PRDX3 is one member of PRDXs localized in the mitochondria, and has been shown to be involved in antioxidant defense and redox signaling. In this study, we investigated the role of PRDX3 in neuronal trauma using a traumatic neuronal injury (TNI) model in primary cultured cortical neurons. We found that TNI significantly decreased the expression of PRDX3 at both mRNA and protein levels. Overexpression of PRDX3 by lentivirus (LV-PRDX3) transfection attenuated lactate dehydrogenase (LDH) release and neuronal apoptosis after TNI. The results of immunostaining showed that LV-PRDX3 transfection markedly reduced TNI-induced intracellular ROS production, protein radical formation and lipid peroxidation. In addition, overexpression of PRDX3 preserved mitochondrial membrane potential (MMP) levels and ATP generation, and inhibited mitochondrial cytochrome c release in TNI-injured neurons. The results of polymerase chain reaction (PCR) showed that PRDX3 overexpression also increased mitochondrial DNA (mtDNA) content and upregulated the expression of mitochondrial biogenesis-related factors. Taken together, our data demonstrate that PRDX3 protects against TNI insult by preserving mitochondrial function and mitochondrial biogenesis, and may have potential therapeutic value for traumatic brain injury (TBI).

Classical Molecular Dynamics with Mobile Protons

An important limitation of standard classical molecular dynamics simulations is the inability to make or break chemical bonds. This restricts severely our ability to study processes that involve even the simplest of chemical reactions, the transfer of a proton. Existing approaches for allowing proton transfer in the context of classical mechanics are rather cumbersome and have not achieved widespread use and routine status. Here we reconsider the combination of molecular dynamics with periodic stochastic proton hops. To ensure computational efficiency, we propose a non-Boltzmann acceptance criterion that is heuristically adjusted to maintain the correct or desirable thermodynamic equilibria between different protonation states and proton transfer rates. Parameters are proposed for hydronium, Asp, Glu, and His. The algorithm is implemented in the program CHARMM and tested on proton diffusion in bulk water and carbon nanotubes and on proton conductance in the gramicidin A channel. Using hopping parameters determined from proton diffusion in bulk water, the model reproduces the enhanced proton diffusivity in carbon nanotubes and gives a reasonable estimate of the proton conductance in gramicidin A.

Implications of short time scale dynamics on long time processes

This review provides a comprehensive overview of the structural dynamics in topical gas- and condensed-phase systems on multiple length and time scales. Starting from vibrationally induced dissociation of small molecules in the gas phase, the question of vibrational and internal energy redistribution through conformational dynamics is further developed by considering coupled electron/proton transfer in a model peptide over many orders of magnitude. The influence of the surrounding solvent is probed for electron transfer to the solvent in hydrated I-. Next, the dynamics of a modified PDZ domain over many time scales is analyzed following activation of a photoswitch. The hydration dynamics around halogenated amino acid side chains and their structural dynamics in proteins are relevant for iodinated TyrB26 insulin. Binding of nitric oxide to myoglobin is a process for which experimental and computational analyses have converged to a common view which connects rebinding time scales and the underlying dynamics. Finally, rhodopsin is a paradigmatic system for multiple length- and time-scale processes for which experimental and computational methods provide valuable insights into the functional dynamics. The systems discussed here highlight that for a comprehensive understanding of how structure, flexibility, energetics, and dynamics contribute to functional dynamics, experimental studies in multiple wavelength regions and computational studies including quantum, classical, and more coarse grained levels are required.