Search

Photosystem II (PSII) catalyzes light-driven water oxidation that releases dioxygen into our atmosphere and provides the electrons needed for the synthesis of biomass. The catalysis occurs in the oxygen-evolving oxo-manganese-calcium (Mn 4 O 5 Ca) cluster that drives the oxidation and deprotonation of substrate water molecules leading to the O 2 formation. However, despite recent advances, the mechanism of these reactions remains unclear and much debated. Here, we show that the light-driven Tyr161 D1 (Y z ) oxidation adjacent to the Mn 4 O 5 Ca cluster, decreases the barrier for proton transfer from the putative substrate water molecule (W3/W x ) to Glu310 D2 , accessible to the luminal bulk. By combining hybrid quantum/classical (QM/MM) free energy calculations with atomistic molecular dynamics simulations, we probe the energetics of the proton transfer along the Cl1 pathway. We demonstrate that the proton transfer occurs via water molecules and a cluster of conserved carboxylates, driven by redox-triggered electric fields directed along the pathway. Glu65 D1 establishes a local molecular gate that controls the proton transfer to the luminal bulk, while Glu312 D2 acts as a local proton storage site. The identified gating region could be important in preventing backflow of protons to the Mn 4 O 5 Ca cluster. The structural changes, derived here based on the dark-state PSII structure, strongly support recent time-resolved X-ray free electron laser data of the S 3 → S 4 transition (Bhowmick et al . Nature 617 , 2023) and reveal the mechanistic basis underlying deprotonation of the substrate water molecules. Our findings provide insight into the water oxidation mechanism of PSII and show how the interplay between redox-triggered electric fields, ion-pairs, and hydration effects control proton transport reactions., Competing Interests: Competing interests statement:The authors declare no competing interest.

Tuberculosis is one of the most common causes of death worldwide, with a rapid emergence of multi-drug-resistant strains underscoring the need for new antituberculosis drugs. Recent studies indicate that lansoprazole-a known gastric proton pump inhibitor and its intracellular metabolite, lansoprazole sulfide (LPZS)-are potential antituberculosis compounds. Yet, their inhibitory mechanism and site of action still remain unknown. Here, we combine biochemical, computational, and structural approaches to probe the interaction of LPZS with the respiratory chain supercomplex III 2 IV 2 of Mycobacterium smegmatis , a close homolog of Mycobacterium tuberculosis supercomplex. We show that LPZS binds to the Q o cavity of the mycobacterial supercomplex, inhibiting the quinol substrate oxidation process and the activity of the enzyme. We solve high-resolution (2.6 Å) cryo-electron microscopy (cryo-EM) structures of the supercomplex with bound LPZS that together with microsecond molecular dynamics simulations, directed mutagenesis, and functional assays reveal key interactions that stabilize the inhibitor, but also how mutations can lead to the emergence of drug resistance. Our combined findings reveal an inhibitory mechanism of LPZS and provide a structural basis for drug development against tuberculosis., Competing Interests: Competing interests statement:The authors declare no competing interest.

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., (© 2024. The Author(s).)

The alternative oxidase (AOX) is a membrane-bound di-iron enzyme that catalyzes O 2 -driven quinol oxidation in the respiratory chains of plants, fungi, and several pathogenic protists of biomedical and industrial interest. Yet, despite significant biochemical and structural efforts over the last decades, the catalytic principles of AOX remain poorly understood. We develop here multi-scale quantum and classical molecular simulations in combination with biochemical experiments to address the proton-coupled electron transfer (PCET) reactions responsible for catalysis in AOX from Trypanosoma brucei , the causative agent of sleeping sickness. We show that AOX activates and splits dioxygen via a water-mediated PCET reaction, resulting in a high-valent ferryl/ferric species and tyrosyl radical (Tyr220˙) that drives the oxidation of the quinol via electric field effects. We identify conserved carboxylates (Glu215, Asp100) within a buried cluster of ion-pairs that act as a transient proton-loading site in the quinol oxidation process, and validate their function experimentally with point mutations that result in drastic activity reduction and p K a -shifts. Our findings provide a key mechanistic understanding of the catalytic machinery of AOX, as well as a molecular basis for rational drug design against energy transduction chains of parasites. On a general level, our findings illustrate how redox-triggered conformational changes in ion-paired networks control the catalysis via electric field effects., Competing Interests: A. L. M. holds patents and financial interests in the development of phytopathogenic fungicides. The other authors declare no competing financial interest., (This journal is © The Royal Society of Chemistry.)

Hsp90 is a molecular chaperone of central importance for protein homeostasis in the cytosol of eukaryotic cells, with key functional and structural traits conserved from yeast to man. During evolution, Hsp90 has gained additional functional importance, leading to an increased number of interacting co-chaperones and client proteins. Here, we show that the overall conformational transitions coupled to the ATPase cycle of Hsp90 are conserved from yeast to humans, but cycle timing as well as the dynamics are significantly altered. In contrast to yeast Hsp90, the human Hsp90 is characterized by broad ensembles of conformational states, irrespective of the absence or presence of ATP. The differences in the ATPase rate and conformational transitions between yeast and human Hsp90 are based on two residues in otherwise conserved structural elements that are involved in triggering structural changes in response to ATP binding. The exchange of these two mutations allows swapping of the ATPase rate and of the conformational transitions between human and yeast Hsp90. Our combined results show that Hsp90 evolved to a protein with increased conformational dynamics that populates ensembles of different states with strong preferences for the N-terminally open, client-accepting states., (© 2024. The Author(s).)

Botulinum neurotoxins (BoNTs) are the most potent toxins known and are used to treat an increasing number of medical disorders. All BoNTs are naturally co-expressed with a protective partner protein (NTNH) with which they form a 300 kDa complex, to resist acidic and proteolytic attack from the digestive tract. We have previously identified a new botulinum neurotoxin serotype, BoNT/X, that has unique and therapeutically attractive properties. We present the cryo-EM structure of the BoNT/X-NTNH/X complex and the crystal structure of the isolated NTNH protein. Unexpectedly, the BoNT/X complex is stable and protease-resistant at both neutral and acidic pH and disassembles only in alkaline conditions. Using the stabilizing effect of NTNH, we isolated BoNT/X and showed that it has very low potency both in vitro and in vivo. Given the high catalytic activity and translocation efficacy of BoNT/X, low activity of the full toxin is likely due to the receptor-binding domain, which presents very weak ganglioside binding and exposed hydrophobic surfaces., (© 2024. The Author(s).)

Water-mediated proton transfer reactions are central for catalytic processes in a wide range of biochemical systems, ranging from biological energy conversion to chemical transformations in the metabolism. Yet, the accurate computational treatment of such complex biochemical reactions is highly challenging and requires the application of multiscale methods, in particular hybrid quantum/classical (QM/MM) approaches combined with free energy simulations. Here, we combine the unique exploration power of new advanced sampling methods with density functional theory (DFT)-based QM/MM free energy methods for multiscale simulations of long-range protonation dynamics in biological systems. In this regard, we show that combining multiple walkers/well-tempered metadynamics with an extended system adaptive biasing force method (MWE) provides a powerful approach for exploration of water-mediated proton transfer reactions in complex biochemical systems. We compare and combine the MWE method also with QM/MM umbrella sampling and explore the sampling of the free energy landscape with both geometric (linear combination of proton transfer distances) and physical (center of excess charge) reaction coordinates and show how these affect the convergence of the potential of mean force (PMF) and the activation free energy. We find that the QM/MM-MWE method can efficiently explore both direct and water-mediated proton transfer pathways together with forward and reverse hole transfer mechanisms in the highly complex proton channel of respiratory Complex I, while the QM/MM-US approach shows a systematic convergence of selected long-range proton transfer pathways. In this regard, we show that the PMF along multiple proton transfer pathways is recovered by combining the strengths of both approaches in a QM/MM-MWE/focused US (FUS) scheme and reveals new mechanistic insight into the proton transfer principles of Complex I. Our findings provide a promising basis for the quantitative multiscale simulations of long-range proton transfer reactions in biological systems.

The membrane-bound hydrogenase (Mbh) from Pyrococcus furiosus is an archaeal member of the Complex I superfamily. It catalyzes the reduction of protons to H 2 gas powered by a [NiFe] active site and transduces the free energy into proton pumping and Na + /H + exchange across the membrane. Despite recent structural advances, the mechanistic principles of H 2 catalysis and ion transport in Mbh remain elusive. Here, we probe how the redox chemistry drives the reduction of the proton to H 2 and how the catalysis couples to conformational dynamics in the membrane domain of Mbh. By combining large-scale quantum chemical density functional theory (DFT) and correlated ab initio wave function methods with atomistic molecular dynamics simulations, we show that the proton transfer reactions required for the catalysis are gated by electric field effects that direct the protons by water-mediated reactions from Glu21 L toward the [NiFe] site, or alternatively along the nearby His75 L pathway that also becomes energetically feasible in certain reaction steps. These local proton-coupled electron transfer (PCET) reactions induce conformational changes around the active site that provide a key coupling element via conserved loop structures to the ion transport activity. We find that H 2 forms in a heterolytic proton reduction step, with spin crossovers tuning the energetics along key reaction steps. On a general level, our work showcases the role of electric fields in enzyme catalysis and how these effects are employed by the [NiFe] active site of Mbh to drive PCET reactions and ion transport.

Aerobic life is powered by membrane-bound redox enzymes that shuttle electrons to oxygen and transfer protons across a biological membrane. Structural studies suggest that these energy-transducing enzymes operate as higher-order supercomplexes, but their functional role remains poorly understood and highly debated. Here we resolve the functional dynamics of the 0.7 MDa III 2 IV 2 obligate supercomplex from Mycobacterium smegmatis, a close relative of M. tuberculosis, the causative agent of tuberculosis. By combining computational, biochemical, and high-resolution (2.3 Å) cryo-electron microscopy experiments, we show how the mycobacterial supercomplex catalyses long-range charge transport from its menaquinol oxidation site to the binuclear active site for oxygen reduction. Our data reveal proton and electron pathways responsible for the charge transfer reactions, mechanistic principles of the quinone catalysis, and how unique molecular adaptations, water molecules, and lipid interactions enable the proton-coupled electron transfer (PCET) reactions. Our combined findings provide a mechanistic blueprint of mycobacterial supercomplexes and a basis for developing drugs against pathogenic bacteria., (© 2024. The Author(s).)

Aerobic ribonucleotide reductases (RNRs) initiate synthesis of DNA building blocks by generating a free radical within the R2 subunit; the radical is subsequently shuttled to the catalytic R1 subunit through proton-coupled electron transfer (PCET). We present a high-resolution room temperature structure of the class Ie R2 protein radical captured by x-ray free electron laser serial femtosecond crystallography. The structure reveals conformational reorganization to shield the radical and connect it to the translocation path, with structural changes propagating to the surface where the protein interacts with the catalytic R1 subunit. Restructuring of the hydrogen bond network, including a notably short O-O interaction of 2.41 angstroms, likely tunes and gates the radical during PCET. These structural results help explain radical handling and mobilization in RNR and have general implications for radical transfer in proteins.

Directional ion transport across biological membranes plays a central role in many cellular processes. Elucidating the molecular determinants for vectorial ion transport is key to understanding the functional mechanism of membrane-bound ion pumps. The extensive investigation of the light-driven proton pump bacteriorhodopsin from Halobacterium salinarum ( Hs BR) enabled a detailed description of outward proton transport. Although the structure of inward-directed proton pumping rhodopsins is very similar to Hs BR, little is known about their protonation pathway, and hence, the molecular reasons for the vectoriality of proton translocation remain unclear. Here, we employ a combined experimental and theoretical approach to tracking protonation steps in the light-driven inward proton pump xenorhodopsin from Nanosalina sp. ( Ns XeR). Time-resolved infrared spectroscopy reveals the transient deprotonation of D220 concomitantly with deprotonation of the retinal Schiff base. Our molecular dynamics simulations support a proton release pathway from the retinal Schiff base via a hydrogen-bonded water wire leading to D220 that could provide a putative gating point for the proton release and with allosteric interactions to the retinal Schiff base. Our findings support the key role of D220 in mediating proton release to the cytoplasmic side and provide evidence that this residue is not the primary proton acceptor of the proton transiently released by the retinal Schiff base.

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 (TM3 ND6 ) 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.

Recent breakthroughs in neural network-based structure prediction methods, such as AlphaFold2 and RoseTTAFold, have dramatically improved the quality of computational protein structure prediction. These models also provide statistical confidence scores that can estimate uncertainties in the predicted structures, but it remains unclear to what extent these scores are related to the intrinsic conformational dynamics of proteins. Here, we compare AlphaFold2 prediction scores with explicit large-scale molecular dynamics simulations of 28 one- and two-domain proteins with varying degrees of flexibility. We demonstrate a strong correlation between the statistical prediction scores and the explicit motion derived from extensive atomistic molecular dynamics simulations and further derive an elastic network model based on the statistical scores of AlphFold2 (AF-ENM), which we benchmark in combination with coarse-grained molecular dynamics simulations. We show that our AF-ENM method reproduces the global protein dynamics with improved accuracy, providing a powerful way to derive effective molecular dynamics using neural network-based structure prediction models.

Zn 2+ is one of the most versatile biologically available metal ions, but accurate modeling of Zn 2+ -containing metalloproteins at the biomolecular force field level can be challenging. Since most Zn 2+ models are parameterized in bulk solvent, in-depth knowledge about their performance in a protein environment is limited. Thus, we systematically investigate here the behavior of non-polarizable Zn 2+ models for their ability to reproduce experimentally determined metal coordination and ligand binding in metalloproteins. The benchmarking is performed in challenging environments, including mono- (carbonic anhydrase II) and bimetallic (metallo-β-lactamase VIM-2) ligand binding sites. We identify key differences in the performance between the Zn 2+ models with regard to the preferred ligating atoms (charged/non-charged), attraction of water molecules, and the preferred coordination geometry. Based on these results, we suggest suitable simulation conditions for varying Zn 2+ site geometries that could guide the further development of biomolecular Zn 2+ models., (© 2022 The Authors. Journal of Computational Chemistry published by Wiley Periodicals LLC.)

Electron bifurcation is a fundamental energy coupling mechanism widespread in microorganisms that thrive under anoxic conditions. These organisms employ hydrogen to reduce CO 2 , but the molecular mechanisms have remained enigmatic. The key enzyme responsible for powering these thermodynamically challenging reactions is the electron-bifurcating [FeFe]-hydrogenase HydABC that reduces low-potential ferredoxins (Fd) by oxidizing hydrogen gas (H 2 ). By combining single-particle cryo-electron microscopy (cryoEM) under catalytic turnover conditions with site-directed mutagenesis experiments, functional studies, infrared spectroscopy, and molecular simulations, we show that HydABC from the acetogenic bacteria Acetobacterium woodii and Thermoanaerobacter kivui employ a single flavin mononucleotide (FMN) cofactor to establish electron transfer pathways to the NAD(P) + and Fd reduction sites by a mechanism that is fundamentally different from classical flavin-based electron bifurcation enzymes. By modulation of the NAD(P) + binding affinity via reduction of a nearby iron-sulfur cluster, HydABC switches between the exergonic NAD(P) + reduction and endergonic Fd reduction modes. Our combined findings suggest that the conformational dynamics establish a redox-driven kinetic gate that prevents the backflow of the electrons from the Fd reduction branch toward the FMN site, providing a basis for understanding general mechanistic principles of electron-bifurcating hydrogenases.

Aerobic life is powered by membrane-bound enzymes that catalyze the transfer of electrons to oxygen and protons across a biological membrane. Cytochrome c oxidase (C c O) functions as a terminal electron acceptor in mitochondrial and bacterial respiratory chains, driving cellular respiration and transducing the free energy from O 2 reduction into proton pumping. Here we show that C c O creates orientated electric fields around a nonpolar cavity next to the active site, establishing a molecular switch that directs the protons along distinct pathways. By combining large-scale quantum chemical density functional theory (DFT) calculations with hybrid quantum mechanics/molecular mechanics (QM/MM) simulations and atomistic molecular dynamics (MD) explorations, we find that reduction of the electron donor, heme a , leads to dissociation of an arginine (Arg438)-heme a 3 D-propionate ion-pair. This ion-pair dissociation creates a strong electric field of up to 1 V Å -1 along a water-mediated proton array leading to a transient proton loading site (PLS) near the active site. Protonation of the PLS triggers the reduction of the active site, which in turn aligns the electric field vectors along a second, "chemical," proton pathway. We find a linear energy relationship of the proton transfer barrier with the electric field strength that explains the effectivity of the gating process. Our mechanism shows distinct similarities to principles also found in other energy-converting enzymes, suggesting that orientated electric fields generally control enzyme catalysis.

Redox reactions are central to biochemistry and are both controlled by and induce protein structural changes. Here, we describe structural rearrangements and crosstalk within the Bacillus cereus ribonucleotide reductase R2b-NrdI complex, a di-metal carboxylate-flavoprotein system, as part of the mechanism generating the essential catalytic free radical of the enzyme. Femtosecond crystallography at an X-ray free electron laser was utilized to obtain structures at room temperature in defined redox states without suffering photoreduction. Together with density functional theory calculations, we show that the flavin is under steric strain in the R2b-NrdI protein complex, likely tuning its redox properties to promote superoxide generation. Moreover, a binding site in close vicinity to the expected flavin O 2 interaction site is observed to be controlled by the redox state of the flavin and linked to the channel proposed to funnel the produced superoxide species from NrdI to the di-manganese site in protein R2b. These specific features are coupled to further structural changes around the R2b-NrdI interaction surface. The mechanistic implications for the control of reactive oxygen species and radical generation in protein R2b are discussed., Competing Interests: JJ, OA, VS, PS, IK, AB, PS, MD, CP, SG, KS, PA, AB, AO, MC, SO, KT, FF, AB, AB, NS, VY, JY, VK, JK, HL, MH No competing interests declared, (© 2022, John et al.)

Dioxygenases catalyze stereoselective oxygen atom transfer in metabolic pathways of biological, industrial, and pharmaceutical importance, but their precise chemical principles remain controversial. The α-ketoglutarate (αKG)-dependent dioxygenase AsqJ synthesizes biomedically active quinolone alkaloids via desaturation and subsequent epoxidation of a carbon-carbon bond in the cyclopeptin substrate. Here, we combine high-resolution X-ray crystallography with enzyme engineering, quantum-classical (QM/MM) simulations, and biochemical assays to describe a peroxidic intermediate that bridges the substrate and active site metal ion in AsqJ. Homolytic cleavage of this moiety during substrate epoxidation generates an activated high-valent ferryl (Fe IV = O) species that mediates the next catalytic cycle, possibly without the consumption of the metabolically valuable αKG cosubstrate. Our combined findings provide an important understanding of chemical bond activation principles in complex enzymatic reaction networks and molecular mechanisms of dioxygenases.

The mitochondrial electron transport chain maintains the proton motive force that powers adenosine triphosphate (ATP) synthesis. The energy for this process comes from oxidation of reduced nicotinamide adenine dinucleotide (NADH) and succinate, with the electrons from this oxidation passed via intermediate carriers to oxygen. Complex IV (CIV), the terminal oxidase, transfers electrons from the intermediate electron carrier cytochrome c to oxygen, contributing to the proton motive force in the process. Within CIV, protons move through the K and D pathways during turnover. The former is responsible for transferring two protons to the enzyme's catalytic site upon its reduction, where they eventually combine with oxygen and electrons to form water. CIV is the main site for respiratory regulation, and although previous studies showed that steroid binding can regulate CIV activity, little is known about how this regulation occurs. Here, we characterize the interaction between CIV and steroids using a combination of kinetic experiments, structure determination, and molecular simulations. We show that molecules with a sterol moiety, such as glyco-diosgenin and cholesteryl hemisuccinate, reversibly inhibit CIV. Flash photolysis experiments probing the rapid equilibration of electrons within CIV demonstrate that binding of these molecules inhibits proton uptake through the K pathway. Single particle cryogenic electron microscopy (cryo-EM) of CIV with glyco-diosgenin reveals a previously undescribed steroid binding site adjacent to the K pathway, and molecular simulations suggest that the steroid binding modulates the conformational dynamics of key residues and proton transfer kinetics within this pathway. The binding pose of the sterol group sheds light on possible structural gating mechanisms in the CIV catalytic cycle.

The heat shock protein 90 (Hsp90) is a molecular chaperone central to client protein folding and maturation in eukaryotic cells. During its chaperone cycle, Hsp90 undergoes ATPase-coupled large-scale conformational changes between open and closed states, where the N-terminal and middle domains of the protein form a compact dimerized conformation. However, the molecular principles of the switching motion between the open and closed states remain poorly understood. Here we show by integrating atomistic and coarse-grained molecular simulations with small-angle X-ray scattering experiments and NMR spectroscopy data that Hsp90 exhibits rich conformational dynamics modulated by the charged linker, which connects the N-terminal with the middle domain of the protein. We show that the dissociation of these domains is crucial for the conformational flexibility of the open state, with the separation distance controlled by a β-sheet motif next to the linker region. Taken together, our results suggest that the conformational ensemble of Hsp90 comprises highly extended states, which could be functionally crucial for client processing., Competing Interests: Conflict of interest The authors declare that they have no conflicts of interest with the contents of this article., (Copyright © 2022 The Authors. Published by Elsevier Inc. All rights reserved.)

Photosystem II (PSII) catalyzes light-driven water oxidization, releasing O 2 into the atmosphere and transferring the electrons for the synthesis of biomass. However, despite decades of structural and functional studies, the water oxidation mechanism of PSII has remained puzzling and a major challenge for modern chemical research. Here, we show that PSII catalyzes redox-triggered proton transfer between its oxygen-evolving Mn 4 O 5 Ca cluster and a nearby cluster of conserved buried ion-pairs, which are connected to the bulk solvent via a proton pathway. By using multi-scale quantum and classical simulations, we find that oxidation of a redox-active Tyr z (Tyr161) lowers the reaction barrier for the water-mediated proton transfer from a Ca 2+ -bound water molecule (W3) to Asp61 via conformational changes in a nearby ion-pair (Asp61/Lys317). Deprotonation of this W3 substrate water triggers its migration toward Mn1 to a position identified in recent X-ray free-electron laser (XFEL) experiments [Ibrahim et al. Proc. Natl. Acad. Sci. USA 2020, 117, 12,624-12,635]. Further oxidation of the Mn 4 O 5 Ca cluster lowers the proton transfer barrier through the water ligand sphere of the Mn 4 O 5 Ca cluster to Asp61 via a similar ion-pair dissociation process, while the resulting Mn-bound oxo/oxyl species leads to O 2 formation by a radical coupling mechanism. The proposed redox-coupled protonation mechanism shows a striking resemblance to functional motifs in other enzymes involved in biological energy conversion, with an interplay between hydration changes, ion-pair dynamics, and electric fields that modulate the catalytic barriers.

Photosystem II (PSII), the water/plastoquinone photo-oxidoreductase, plays a key energy input role in the biosphere. [Formula: see text], the reduced semiquinone form of the nonexchangeable quinone, is often considered capable of a side reaction with O 2 , forming superoxide, but this reaction has not yet been demonstrated experimentally. Here, using chlorophyll fluorescence in plant PSII membranes, we show that O 2 does oxidize [Formula: see text] at physiological O 2 concentrations with a t 1/2 of 10 s. Superoxide is formed stoichiometrically, and the reaction kinetics are controlled by the accessibility of O 2 to a binding site near [Formula: see text], with an apparent dissociation constant of 70 ± 20 µM. Unexpectedly, [Formula: see text] could only reduce O 2 when bicarbonate was absent from its binding site on the nonheme iron (Fe 2+ ) and the addition of bicarbonate or formate blocked the O 2 -dependant decay of [Formula: see text] These results, together with molecular dynamics simulations and hybrid quantum mechanics/molecular mechanics calculations, indicate that electron transfer from [Formula: see text] to O 2 occurs when the O 2 is bound to the empty bicarbonate site on Fe 2+ A protective role for bicarbonate in PSII was recently reported, involving long-lived [Formula: see text] triggering bicarbonate dissociation from Fe 2+ [Brinkert et al , Proc. Natl. Acad. Sci. U.S.A. 113, 12144-12149 (2016)]. The present findings extend this mechanism by showing that bicarbonate release allows O 2 to bind to Fe 2+ and to oxidize [Formula: see text] This could be beneficial by oxidizing [Formula: see text] and by producing superoxide, a chemical signal for the overreduced state of the electron transfer chain., Competing Interests: The authors declare no competing interest., (Copyright © 2022 the Author(s). Published by PNAS.)

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 1 -10 6 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.

The membrane-bound hydrogenase (Mbh) is a redox-driven Na + /H + transporter that employs the energy from hydrogen gas (H 2 ) 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. Cell 2018 , 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.

Cellular respiration is powered by membrane-bound redox enzymes that convert chemical energy into an electrochemical proton gradient and drive the energy metabolism. By combining large-scale classical and quantum mechanical simulations with cryo-electron microscopy data, we resolve here molecular details of conformational changes linked to proton pumping in the mammalian complex I. Our data suggest that complex I deactivation blocks water-mediated proton transfer between a membrane-bound quinone site and proton-pumping modules, decoupling the energy-transduction machinery. We identify a putative gating region at the interface between membrane domain subunits ND1 and ND3/ND4L/ND6 that modulates the proton transfer by conformational changes in transmembrane helices and bulky residues. The region is perturbed by mutations linked to human mitochondrial disorders and is suggested to also undergo conformational changes during catalysis of simpler complex I variants that lack the "active"-to-"deactive" transition. Our findings suggest that conformational changes in transmembrane helices modulate the proton transfer dynamics by wetting/dewetting transitions and provide important functional insight into the mammalian respiratory complex I., Competing Interests: The authors declare no competing interest., (Copyright © 2021 the Author(s). Published by PNAS.)

Bacteria power their energy metabolism using membrane-bound respiratory enzymes that capture chemical energy and transduce it by pumping protons or Na + ions across their cell membranes. Recent breakthroughs in molecular bioenergetics have elucidated the architecture and function of many bacterial respiratory enzymes, although key mechanistic principles remain debated. In this Review, we present an overview of the structure, function and bioenergetic principles of modular bacterial respiratory chains and discuss their differences from the eukaryotic counterparts. We also discuss bacterial supercomplexes, which provide central energy transduction systems in several bacteria, including important pathogens, and which could open up possible avenues for treatment of disease.

Complex I functions as a primary redox-driven proton pump in aerobic respiratory chains, establishing a proton motive force that powers ATP synthesis and active transport. Recent cryo-electron microscopy (cryo-EM) experiments have resolved the mammalian complex I in the biomedically relevant active (A) and deactive (D) states (Zhu et al., 2016; Fiedorczuk et al., 2016; Agip et al., 2018 [1-3]) that could regulate enzyme turnover, but it still remains unclear how the conformational state and activity are linked. We show here how global motion along the A/D transition accumulates molecular strain at specific coupling regions important for both redox chemistry and proton pumping. Our data suggest that the A/D motion modulates force propagation pathways between the substrate-binding site and the proton pumping machinery that could alter electrostatic and conformational coupling across large distances. Our findings provide a molecular basis to understand how global protein dynamics can modulate the biological activity of large molecular complexes., (Copyright © 2021. Published by Elsevier B.V.)

Photoredox catalysts are integral components of artificial photosystems, and have recently emerged as powerful tools for catalysing numerous organic reactions. However, the development of inexpensive and efficient earth-abundant photoredox catalysts remains a challenge. We here present the photochemical and photophysical properties of a Ni-Mabiq catalyst ([Ni II (Mabiq)]OTf ( 1 ); Mabiq = 2-4:6-8-bis(3,3,4,4-tetramethyldihydropyrrolo)-10-15-(2,2-biquinazolino)-[15]-1,3,5,8,10,14-hexaene1,3,7,9,11,14-N 6 )-and of a Zn-containing analogue ([Zn II (Mabiq)OTf] ( 2 ))-using steady state and time resolved optical spectroscopy, time-dependent density functional theory (TDDFT) calculations, and reactivity studies. The Ni and Zn complexes exhibit similar absorption spectra, but markedly different photochemical properties. These differences arise because the excited states of 2 are ligand-localized, whereas metal-centered states account for the photoactivity of 1 . The distinct properties of the Ni and Zn complexes are manifest in their behavior in the photo-driven aza-Henry reaction and oxidative coupling of methoxybenzylamine., Competing Interests: There are no conflicts to declare., (This journal is © The Royal Society of Chemistry.)

Versatile and ecofriendly methods to perform oxidations at near-neutral pH are of crucial importance for processes aimed at purifying water. Chitosan, a deacetylated form of chitin, is a promising starting material owing to its biocompatibility and ability to form stable films and complexes with metals. Here, we report a novel chitosan-based organometallic complex that was tested both as homogeneous and heterogeneous catalyst in the degradation of contaminants of emerging concern in water. The stoichiometry of the complex was experimentally verified with different metals, namely, Cu(II), Fe(III), Fe(II), Co(II), Pd(II), and Mn(II), and we identified the chitosan-Fe(III) complex as the most efficient catalyst. This complex effectively degraded phenol, triclosan, and 3-chlorophenol in the presence of hydrogen peroxide. A putative ferryl-mediated reaction mechanism is proposed based on experimental data, density functional theory calculations, and kinetic modeling. Finally, a film of the chitosan-Fe(III) complex was synthesized and proven a promising supported heterogeneous catalyst for water purification., (Copyright © 2021 The Authors. Published by Elsevier B.V. All rights reserved.)

Soluble proteins are universally packed with a hydrophobic core and a polar surface that drive the protein folding process. Yet charged networks within the central protein core are often indispensable for the biological function. Here, we show that natural buried ion-pairs are stabilised by amphiphilic residues that electrostatically shield the charged motif from its surroundings to gain structural stability. To explore this effect, we build artificial proteins with buried ion-pairs by combining directed computational design and biophysical experiments. Our findings illustrate how perturbation in charged networks can introduce structural rearrangements to compensate for desolvation effects. We validate the physical principles by resolving high-resolution atomic structures of the artificial proteins that are resistant towards unfolding at extreme temperatures and harsh chemical conditions. Our findings provide a molecular understanding of functional charged networks and how point mutations may alter the protein's conformational landscape.

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.

Rational design of artificial water-splitting catalysts is central for developing new sustainable energy technology. However, the catalytic efficiency of the natural light-driven water-splitting enzyme, photosystem II, has been remarkably difficult to achieve artificially. Here we study the molecular mechanism of ruthenium-based molecular catalysts by integrating quantum chemical calculations with inorganic synthesis and functional studies. By employing correlated ab initio calculations, we show that the thermodynamic driving force for the catalysis is obtained by modulation of π-stacking dispersion interactions within the catalytically active dimer core, supporting recently suggested mechanistic principles of Ru-based water-splitting catalysts. The dioxygen bond forms in a semi-concerted radical coupling mechanism, similar to the suggested water-splitting mechanism in photosystem II. By rationally tuning the dispersion effects, we design a new catalyst with a low activation barrier for the water-splitting. The catalytic principles are probed by synthesis, structural, and electrochemical characterization of the new catalyst, supporting enhanced water-splitting activity under the examined conditions. Our combined findings show that modulation of dispersive interactions provides a rational catalyst design principle for controlling challenging chemistries., Competing Interests: There are no conflicts to declare., (This journal is © The Royal Society of Chemistry.)

In biology, self-assembly of proteins and energy-consuming reaction cycles are intricately coupled. For example, tubulin is activated and deactivated for assembly by a guanosine triphosphate (GTP)-driven reaction cycle, and the emerging microtubules catalyze this reaction cycle by changing the microenvironment of the activated tubulin. Recently, synthetic analogs of chemically fueled assemblies have emerged, but examples in which assembly and reaction cycles are reciprocally coupled remain rare. In this work, we report a peptide that can be activated and deactivated for self-assembly. The emerging assemblies change the microenvironment of their building blocks, which consequently accelerate the rates of building block deactivation and reactivation. We quantitatively understand the mechanisms at play, and we are thus able to tune the catalysis by molecular design of the peptide precursor.

Respiratory complex I (NADH:ubiquinone oxidoreductase) captures the free energy from oxidising NADH and reducing ubiquinone to drive protons across the mitochondrial inner membrane and power oxidative phosphorylation. Recent cryo-EM analyses have produced near-complete models of the mammalian complex, but leave the molecular principles of its long-range energy coupling mechanism open to debate. Here, we describe the 3.0-Å resolution cryo-EM structure of complex I from mouse heart mitochondria with a substrate-like inhibitor, piericidin A, bound in the ubiquinone-binding active site. We combine our structural analyses with both functional and computational studies to demonstrate competitive inhibitor binding poses and provide evidence that two inhibitor molecules bind end-to-end in the long substrate binding channel. Our findings reveal information about the mechanisms of inhibition and substrate reduction that are central for understanding the principles of energy transduction in mammalian 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.

Type II polyketide synthases (PKSs) are multi-enzyme complexes that produce secondary metabolites of medical relevance. Chemical backbones of such polyketides are produced by minimal PKS systems that consist of a malonyl transacylase, an acyl carrier protein and an α/β heterodimeric ketosynthase. Here, we present X-ray structures of all ternary complexes that constitute the minimal PKS system for anthraquinone biosynthesis in Photorhabdus luminescens. In addition, we characterize this invariable core using molecular simulations, mutagenesis experiments and functional assays. We show that malonylation of the acyl carrier protein is accompanied by major structural rearrangements in the transacylase. Principles of an ongoing chain elongation are derived from the ternary complex with a hexaketide covalently linking the heterodimeric ketosynthase with the acyl carrier protein. Our results for the minimal PKS system provide mechanistic understanding of PKSs and a fundamental basis for engineering PKS pathways for future applications.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

While chemical steps involved in bioactive cembranoid biosynthesis have been examined, the corresponding enzymatic mechanisms leading to their formation remain elusive. In the tobacco plant, Nicotiana tabacum, a putative cembratriene-ol synthase (CBTS) initiates the catalytic cascade that lead to the biosynthesis of cembratriene-4,6-diols, which displays antibacterial- and anti-proliferative activities. We report here on structural homology models, functional studies, and mechanistic explorations of this enzyme using a combination of biosynthetic and computational methods. This approach guided us to develop an efficient de novo production of five bioactive non- and monohydroxylated cembranoids. Our homology models in combination with quantum and classical simulations suggested putative principles of the CBTS catalytic cycle, and provided a possible rationale for the formation of premature olefinic side products. Moreover, the functional reconstruction of a N. tabacum -derived class II P450 with a cognate CPR, obtained by transcriptome mining provided for production of bioactive cembratriene-4,6-diols. Our combined findings provide mechanistic insights into cembranoid biosynthesis, and a basis for the sustainable industrial production of highly valuable bioactive cembranoids., Competing Interests: The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (© 2020 The Author(s).)

The heat shock protein 90 (Hsp90) is a molecular chaperone that employs the free energy of ATP hydrolysis to control the folding and activation of several client proteins in the eukaryotic cell. To elucidate how the local ATPase reaction in the active site couples to the global conformational dynamics of Hsp90, we integrate here large-scale molecular simulations with biophysical experiments. We show that the conformational switching of conserved ion pairs between the N-terminal domain, harbouring the active site, and the middle domain strongly modulates the catalytic barrier of the ATP-hydrolysis reaction by electrostatic forces. Our combined findings provide a mechanistic model for the coupling between catalysis and protein dynamics in Hsp90, and show how long-range coupling effects can modulate enzymatic activity.

Methylation of a conserved lysine in C-terminal domain of the molecular chaperone Hsp90 was shown previously to affect its in vivo function. However, the underlying mechanism remained elusive. Through a combined experimental and computational approach, this study shows that this site is very sensitive to sidechain modifications and crucial for Hsp90 activity in vitro and in vivo. Our results demonstrate that this particular lysine serves as a switch point for the regulation of Hsp90 functions by influencing its conformational cycle, ATPase activity, co-chaperone regulation, and client activation of yeast and human Hsp90. Incorporation of the methylated lysine via genetic code expansion specifically shows that upon modification, the conformational cycle of Hsp90 is altered. Molecular dynamics simulations including the methylated lysine suggest specific conformational changes that are propagated through Hsp90. Thus, methylation of the C-terminal lysine allows a precise allosteric tuning of Hsp90 activity via long distances.

Photosynthetic organisms capture light energy to drive their energy metabolism, and employ the chemical reducing power to convert carbon dioxide (CO 2 ) into organic molecules. Photorespiration, however, significantly reduces the photosynthetic yields. To survive under low CO 2 concentrations, cyanobacteria evolved unique carbon-concentration mechanisms that enhance the efficiency of photosynthetic CO 2 fixation, for which the molecular principles have remained unknown. We show here how modular adaptations enabled the cyanobacterial photosynthetic complex I to concentrate CO 2 using a redox-driven proton-pumping machinery. Our cryo-electron microscopy structure at 3.2 Å resolution shows a catalytic carbonic anhydrase module that harbours a Zn 2+ active site, with connectivity to proton-pumping subunits that are activated by electron transfer from photosystem I. Our findings illustrate molecular principles in the photosynthetic complex I machinery that enabled cyanobacteria to survive in drastically changing CO 2 conditions.

Cells employ membrane-embedded antiporter proteins to control their pH, salt concentration, and volume. The large family of cation/proton antiporters is dominated by Na + /H + antiporters that exchange sodium ions against protons, but homologous K + /H + exchangers have recently been characterized. We show experimentally that the electroneutral antiporter NhaP1 of Methanocaldococcus jannaschii (MjNhaP1) is highly selective for Na + ions. We then characterize the ion selectivity in both the inward-open and outward-open states of MjNhaP1 using classical molecular dynamics simulations, free energy calculations, and hybrid quantum/classical (QM/MM) simulations. We show that MjNhaP1 is highly selective for binding of Na + over K + in the inward-open state, yet it is only weakly selective in the outward-open state. These findings are consistent with the function of MjNhaP1 as a sodium-driven deacidifier of the cytosol that maintains a high cytosolic K + concentration in environments of high salinity. By combining experiment and computation, we gain mechanistic insight into the Na + /H + transport mechanism and help elucidate the molecular basis for ion selectivity in cation/proton exchangers.

Quantum chemical calculations are important for elucidating light-capturing mechanisms in photobiological systems. The time-dependent density functional theory (TDDFT) has become a popular methodology because of its balance between accuracy and computational scaling, despite its problems in describing, for example, charge transfer states. As a step toward systematically understanding the performance of TDDFT calculations on biomolecular systems, we study here 17 commonly used density functionals, including seven long-range separated functionals, and compare the obtained results with excitation energies calculated at the approximate second order coupled-cluster theory level (CC2). The benchmarking set includes the first five singlet excited states of 11 chemical analogues of biochromophores from the green fluorescent protein, rhodopsin/bacteriorhodopsin (Rh/bR), and the photoactive yellow protein. We find that commonly used pure density functionals such as BP86, PBE, M11-L, and hybrid functionals with 20-25% of Hartree-Fock (HF) exchange (B3LYP, PBE0) have a tendency to consistently underestimate vertical excitation energies (VEEs) relative to the CC2 values, whereas hybrid density functionals with around 50% HF exchange such as BHLYP, PBE50, and M06-2X and long-range corrected functionals such as CAM-B3LYP, ωPBE, ωPBEh, ωB97X, ωB97XD, BNL, and M11 overestimate the VEEs. We observe that calculations using the CAM-B3LYP and ωPBEh functionals with 65% and 100% long-range HF exchange, respectively, lead to an overestimation of the VEEs by 0.2-0.3 eV for the benchmarking set. To reduce the systematic error, we introduce here two new empirical functionals, CAMh-B3LYP and ωhPBE0, for which we adjusted the long-range HF exchange to 50%. The introduced parameterization reduces the mean signed average (MSA) deviation to 0.07 eV and the root mean square (rms) deviation to 0.17 eV as compared to the CC2 values. In the present study, TDDFT calculations using the aug-def2-TZVP basis sets, the best performing functionals relative to CC2 are ωhPBE0 (rms = 0.17, MSA = 0.06 eV); CAMh-B3LYP (rms = 0.16, MSA = 0.07 eV); and PBE0 (rms = 0.23, MSA = -0.14 eV). For the popular range-separated CAM-B3LYP functional, we obtain an rms value of 0.31 eV and an MSA value of 0.25 eV, which can be compared with the rms and MSA values of 0.37 and -0.31 eV, respectively, as obtained at the B3LYP level.

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 416 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., (© 2019 Boczek et al.)

The recently discovered SAFit class of inhibitors against the Hsp90 co-chaperone FKBP51 show greater than 10 000-fold selectivity over its closely related paralogue FKBP52. However, the mechanism underlying this selectivity remained unknown. By combining NMR spectroscopy, biophysical and computational methods with mutational analysis, we show that the SAFit molecules bind to a transient pocket in FKBP51. This represents a weakly populated conformation resembling the inhibitor-bound state of FKBP51, suggesting conformational selection rather than induced fit as the major binding mechanism. The inhibitor-bound conformation of FKBP51 is stabilized by an allosteric network of residues located away from the inhibitor-binding site. These residues stabilize the Phe67 side chain in a dynamic outward conformation and are distinct in FKBP52, thus rationalizing the basis for the selectivity of SAFit inhibitors. Our results represent a paradigm for the selective inhibition of transient binding pockets., (© 2019 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim.)

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.

Anthraquinones, a widely distributed class of aromatic natural products, are produced by a type II polyketide synthase system in the Gram-negative bacterium Photorhabdus luminescens . Heterologous expression of the antABCDEFGHI anthraquinone biosynthetic gene cluster in Escherichia coli identified AntI as an unusual lyase, catalysing terminal polyketide shortening prior to formation of the third aromatic ring. Functional in vitro and in vivo analysis of AntI using X-ray crystallography, structure-based mutagenesis, and molecular simulations revealed that AntI converts a defined octaketide to the tricyclic anthraquinone ring via retro-Claisen and Dieckmann reactions. Thus, AntI catalyses a so far unobserved multistep reaction in this PKS system.

Complex I functions as an initial electron acceptor in aerobic respiratory chains that reduces quinone and pumps protons across a biological membrane. This remarkable charge transfer process extends ca. 300 Å and it is initiated by a poorly understood proton-coupled electron transfer (PCET) reaction between nicotinamide adenine dinucleotide (NADH) and a protein-bound flavin (FMN) cofactor. We combine here large-scale density functional theory calculations and quantum/classical models with atomistic molecular dynamics simulations to probe the energetics and dynamics of the NADH-driven PCET reaction in complex I. We find that the reaction takes place by concerted hydrogen atom (H • ) transfer that couples to an electron transfer (eT) between the aromatic ring systems of the cofactors and further triggers reduction of the nearby FeS centers. In bacterial, Escherichia coli-like complex I isoforms, reduction of the N1a FeS center increases the binding affinity of the oxidized NAD + that prevents the nucleotide from leaving prematurely. This electrostatic trapping could provide a protective gating mechanism against reactive oxygen species formation. We also find that proton transfer from the transient FMNH • to a nearby conserved glutamate (Glu97) residue favors eT from N1a onward along the FeS chain and modulates the binding of a new NADH molecule. The PCET in complex I isoforms with low-potential N1a centers is also discussed. On the basis of our combined results, we propose a putative mechanistic model for the NADH-driven proton/electron-transfer reaction in complex I.

The light-harvesting chlorophyll (Chl) molecules of photosynthetic systems form the basis for light-driven energy conversion. In biological environments, the Chl chromophores occur in two distinct diastereotopic configurations, where the α and β configurations have a magnesium-ligating histidine residue and a 17-propionic acid moiety on the opposite side or on the same side of the Chl ring, respectively. Although β-ligated Chl dimers occupy conserved positions around the reaction center of photosystem I (PSI), the functional relevance of the α/β configuration of the ligation is poorly understood. We employ here correlated ab initio calculations using the algebraic-diagrammatic construction through second order (ADC(2)) and the approximate second-order coupled cluster (CC2) methods in combination with the reduced virtual space (RVS) approach in studies of the intrinsic excited-state properties of α-ligated and β-ligated Chl dimers of PSI. Our ab initio calculations suggest that the absorption of the α-ligated reaction-center Chl dimer of PSI is redshifted by 0.13-0.14 eV in comparison to the β-ligated dimers due to combined excitonic coupling and strain effects. We also show that time-dependent density functional theory (TDDFT) calculations using range-separated density functionals underestimate the absorption shift between the α- and β-ligated dimers. Our findings may provide a molecular starting point for understanding the energy flow in natural photosynthetic systems, as well as a blueprint for developing new molecules that convert sunlight into other forms of energy.