Complex I functions as a redox-driven proton pump in aerobic respiratory chains. By reducing quinone (Q), complex I employs the free energy released in the process to thermodynamically drive proton pumping across its membrane domain. The initial Q reduction step plays a central role in activating the proton pumping machinery. In order to probe the energetics, dynamics, and molecular mechanism for the proton-coupled electron transfer process linked to the Q reduction, we employ here multiscale quantum and classical molecular simulations. We identify that both ubiquinone (UQ) and menaquinone (MQ) can form stacking and hydrogen-bonded interactions with the conserved Q-binding-site residue His-38 and that conformational changes between these binding modes modulate the Q redox potentials and the rate of electron transfer (eT) from the terminal N2 iron–sulfur center. We further observe that, while the transient formation of semiquinone is not proton-coupled, the second eT process couples to a semiconcerted proton uptake from conserved tyrosine (Tyr-87) and histidine (His-38) residues within the active site. Our calculations indicate that both UQ and MQ have low redox potentials around −260 and −230 mV, respectively, in the Q-binding site, respectively, suggesting that release of the Q toward the membrane is coupled to an energy transduction step that could thermodynamically drive proton pumping in complex I.

The conversion of light energy into ion gradients across biological membranes is one of the most fundamental reactions in primary biological energy transduction. Recently, the structure of the first light-activated Na+ pump, Krokinobacter eikastus rhodopsin 2 (KR2), was resolved at atomic resolution [Kato HE, et al. (2015) Nature 521:48–53]. To elucidate its molecular mechanism for Na+ pumping, we perform here extensive classical and quantum molecular dynamics (MD) simulations of transient photocycle states. Our simulations show how the dynamics of key residues regulate water and ion access between the bulk and the buried light-triggered retinal site. We identify putative Na+ binding sites and show how protonation and conformational changes gate the ion through these sites toward the extracellular side. We further show by correlated ab initio quantum chemical calculations that the obtained putative photocycle intermediates are in close agreement with experimental transient optical spectroscopic data. The combined results of the ion translocation and gating mechanisms in KR2 may provide a basis for the rational design of novel light-driven ion pumps with optogenetic applications.

Chlorophylls are light-capturing units found in photosynthetic proteins. We study here the ground and excited state properties of monomeric, dimeric, and tetrameric models of the special chlorophyll/bacteriochlorophyll (Chl/BChl) pigment (P) centers P700 and P680/P870 of type I and type II photosystems, respectively. In the excited state calculations, we study the performance of the algebraic diagrammatic construction through second-order (ADC(2)) method in combination with the reduced virtual space (RVS) approach and the recently developed Laplace-transformed scaled-opposite-spin (LT-SOS) algorithm, which allows us, for the first time, to address multimeric effects at correlated ab initio levels using large basis sets. At the LT-SOS-RVS-ADC(2)/def2-TZVP level, we obtain vertical excitation energies (VEEs) of 2.00–2.07 and 1.52–1.62 eV for the P680/P700 and the P870 pigment models, respectively, which agree well with the experimental absorption maxima of 1.82, 1.77, and 1.43 eV for P680, P700, and P870, respectively. In the P680/P870 models, we find that the photoexcitation leads to a π → π* transition in which the exciton is delocalized between the adjacent Chl/BChl molecules of the central pair, whereas the exciton is localized to a single chlorophyll molecule in the P700 model. Consistent with experiments, the calculated excitonic splittings between the central pairs of P680, P700, and P870 models are 80, 200, and 400 cm–1, respectively. The calculations show that the electron affinity of the radical cation of the P680 model is 0.4 V larger than for the P870 model and 0.2 V larger than for P700. The chromophore stacking interaction is found to strongly influence the electron localization properties of the light-absorbing pigments, which may help to elucidate mechanistic details of the charge separation process in type I and type II photosystems.

The Photoactive Yellow Protein (PYP) is a light-driven photoreceptor, responsible for the phototaxis of halophilic bacteria. Recently, a new short-lived intermediate (pR0) was characterized in the PYP photocycle using combined time-resolved X-ray crystallography and density functional theory calculations. The pR0 species was identified as a highly contorted cis-intermediate, which is stabilized by hydrogen bonds with protein residues. Here we show by hybrid quantum mechanics/classical mechanics (QM/MM) molecular dynamics simulations, and first-principles calculations of optical properties, that the optical shifts in the early steps of the PYP photocycle originate from the conversion of light energy into molecular strain, stored in the pR0 state, and its relaxation in subsequent reaction steps. Our calculations quantitatively reproduce experimental data, which enables us to identify molecular origins of the optical shifts. Our combined approach suggests that the short-lived pR0 intermediate stores ∼1/3 of the photon energy as molecular strain, thus providing the thermodynamic driving force for later conformational changes in the protein.

Proton-coupled electron transfer (PCET) processes are elementary chemical reactions involved in a broad range of radical and redox reactions. Elucidating fundamental PCET reaction mechanisms are thus of central importance for chemical and biochemical research. Here we use quantum chemical density functional theory (DFT), time-dependent density functional theory (TDDFT), and the algebraic diagrammatic-construction through second-order (ADC(2)) to study the mechanism, thermodynamic driving force effects, and reaction barriers of both ground state proton transfer (pT) and photoinduced proton-coupled electron transfer (PCET) between nitrosylated phenyl-phenol compounds and hydrogen-bonded t-butylamine as an external base. We show that the obtained reaction barriers for the ground state pT reactions depend linearly on the thermodynamic driving force, with a Brønsted slope of 1 or 0. Photoexcitation leads to a PCET reaction, for which we find that the excited state reaction barrier depends on the thermodynamic driving force with a Brønsted slope of 1/2. To support the mechanistic picture arising from the static potential energy surfaces, we perform additional molecular dynamics simulations on the excited state energy surface, in which we observe a spontaneous PCET between the donor and the acceptor groups. Our findings suggest that a Brønsted analysis may distinguish the ground state pT and excited state PCET processes.

Accurate predictions of photoexcitation properties are a major challenge for modern methods of theoretical chemistry. We show here how approximate coupled-cluster singles and doubles (CC2) calculations in combination with the reduced virtual space (RVS) approach can be employed in studies of excited states of large biomolecular systems. The RVS-CC2 approach is used for accurately predicting optical properties of the p-hydroxybenzylidene-dihydroimidazolinone (p-HBDI) chromophore embedded in green fluorescent protein (GFP) models using quantum mechanical calculations in combination with large basis sets. We study the lowest excited states for the isolated and protein-embedded chromophore in two different protonation states, and show how omitting high-lying virtual orbitals in the RVS calculation of excitation energies renders large-scale CC2 studies computationally feasible. We also discuss how the error introduced by the RVS approach can be systematically estimated and controlled. The obtained CC2 excitation energies of 3.13–3.27 and 2.69–2.77 eV for the two protonation states of different protein models are in excellent agreement with the maxima of the experimental absorption spectra of 3.12–3.14 and 2.61–2.64 eV, respectively. Thus, the calculated energy splitting between the excited states of the two protonation states is 0.44–0.52 eV, which agrees very well with the experimental value of 0.48–0.51 eV. The calculations at the RVS-CC2 level on the protein models show the importance of using large QM regions in studies of biochromophores embedded in proteins.

Complex I functions as a redox-linked proton pump in the respiratory chains of mitochondria and bacteria, driven by the reduction of quinone (Q) by NADH. Remarkably, the distance between the Q reduction site and the most distant proton channels extends nearly 200 Å. To elucidate the molecular origin of this long-range coupling, we apply a combination of large-scale molecular simulations and a site-directed mutagenesis experiment of a key residue. In hybrid quantum mechanics/molecular mechanics simulations, we observe that reduction of Q is coupled to its local protonation by the His-38/Asp-139 ion pair and Tyr-87 of subunit Nqo4. Atomistic classical molecular dynamics simulations further suggest that formation of quinol (QH2) triggers rapid dissociation of the anionic Asp-139 toward the membrane domain that couples to conformational changes in a network of conserved charged residues. Site-directed mutagenesis data confirm the importance of Asp-139; upon mutation to asparagine the Q reductase activity is inhibited by 75%. The current results, together with earlier biochemical data, suggest that the proton pumping in complex I is activated by a unique combination of electrostatic and conformational transitions.

Retinal is the light-absorbing biochromophore responsible for the activation of vision pigments and light-driven ion pumps. Nature has evolved molecular tuning mechanisms that significantly shift the optical properties of the retinal pigments to enable their absorption of visible light. Using large-scale quantum chemical calculations at the density functional theory level combined with frozen density embedding theory, we show here how the protein environment of vision pigments tunes the absorption of retinal by electrostatically dominated interactions between the chromophore and the surrounding protein residues. The calculations accurately reproduce the experimental absorption maxima of rhodopsin and the red, green, and blue color pigments. We further identify key interactions responsible for the color-shifting effects by mutating the rhodopsin structure in silico, and we find that deprotonation of the retinyl is likely to be responsible for the blue-shifted absorption in the blue cone vision pigment.

Complex I serves as the primary electron entry point into the mitochondrial and bacterial respiratory chains. It catalyzes the reduction of quinones by electron transfer from NADH, and couples this exergonic reaction to the translocation of protons against an electrochemical proton gradient. The membrane domain of the enzyme extends ∼180 Å from the site of quinone reduction to the most distant proton pathway. To elucidate possible mechanisms of the long-range proton-coupled electron transfer process, we perform large-scale atomistic molecular dynamics simulations of the membrane domain of complex I from Escherichia coli. We observe spontaneous hydration of a putative proton entry channel at the NuoN/K interface, which is sensitive to the protonation state of buried glutamic acid residues. In hybrid quantum mechanics/classical mechanics simulations, we find that the observed water wires support rapid proton transfer from the protein surface to the center of the membrane domain. To explore the functional relevance of the pseudosymmetric inverted-repeat structures of the antiporter-like subunits NuoL/M/N, we constructed a symmetry-related structure of a possible alternate-access state. In molecular dynamics simulations, we find the resulting structural changes to be metastable and reversible at the protein backbone level. However, the increased hydration induced by the conformational change persists, with water molecules establishing enhanced lateral connectivity and pathways for proton transfer between conserved ionizable residues along the center of the membrane domain. Overall, the observed water-gated transitions establish conduits for the unidirectional proton translocation processes, and provide a possible coupling mechanism for the energy transduction in complex I.

The Photoactive Yellow Protein (PYP) is a light-driven photoreceptor, responsible for the phototaxis of halophilic bacteria. Recently, a new short-lived intermediate (pR0) was characterized in the PYP photocycle using combined time-resolved X-ray crystallography and density functional theory calculations. The pR0 species was identified as a highly contorted cis-intermediate, which is stabilized by hydrogen bonds with protein residues. Here we show by hybrid quantum mechanics/classical mechanics (QM/MM) molecular dynamics simulations, and first-principles calculations of optical properties, that the optical shifts in the early steps of the PYP photocycle originate from the conversion of light energy into molecular strain, stored in the pR0 state, and its relaxation in subsequent reaction steps. Our calculations quantitatively reproduce experimental data, which enables us to identify molecular origins of the optical shifts. Our combined approach suggests that the short-lived pR0 intermediate stores ∼1/3 of the photon energy as molecular strain, thus providing the thermodynamic driving force for later conformational changes in the protein.