Electric fields control water-gated proton transfer in cytochrome c oxidase.

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 (CcO) functions as a terminal electron acceptor in mitochondrial and bacterial respiratory chains, driving cellular respiration and transducing the free energy from O2 reduction into proton pumping. Here we show that CcO 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 a3 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.

Structural basis of mammalian complex IV inhibition by steroids.

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.

Bicarbonate-controlled reduction of oxygen by the QA semiquinone in Photosystem II in membranes.

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 O2, forming superoxide, but this reaction has not yet been demonstrated experimentally. Here, using chlorophyll fluorescence in plant PSII membranes, we show that O2 does oxidize [Formula: see text] at physiological O2 concentrations with a t1/2 of 10 s. Superoxide is formed stoichiometrically, and the reaction kinetics are controlled by the accessibility of O2 to a binding site near [Formula: see text], with an apparent dissociation constant of 70 ± 20 µM. Unexpectedly, [Formula: see text] could only reduce O2 when bicarbonate was absent from its binding site on the nonheme iron (Fe2+) and the addition of bicarbonate or formate blocked the O2-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 O2 occurs when the O2 is bound to the empty bicarbonate site on Fe2+ A protective role for bicarbonate in PSII was recently reported, involving long-lived [Formula: see text] triggering bicarbonate dissociation from Fe2+ [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 O2 to bind to Fe2+ 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.

Mechanistic Principles of Hydrogen Evolution in the Membrane-Bound Hydrogenase.

The membrane-bound hydrogenase (Mbh) from Pyrococcus furiosus is an archaeal member of the Complex I superfamily. It catalyzes the reduction of protons to H2 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 H2 catalysis and ion transport in Mbh remain elusive. Here, we probe how the redox chemistry drives the reduction of the proton to H2 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 Glu21L toward the [NiFe] site, or alternatively along the nearby His75L 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 H2 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.

Deactivation blocks proton pathways in the mitochondrial complex I.

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.

Extended conformational states dominate the Hsp90 chaperone dynamics.

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.

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.

Molecular Basis of the Electron Bifurcation Mechanism in the [FeFe]-Hydrogenase Complex HydABC.

Electron bifurcation is a fundamental energy coupling mechanism widespread in microorganisms that thrive under anoxic conditions. These organisms employ hydrogen to reduce CO2, 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 (H2). 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.

Benchmarking biomolecular force field-based Zn2+ for mono- and bimetallic ligand binding sites.

Zn2+ is one of the most versatile biologically available metal ions, but accurate modeling of Zn2+ -containing metalloproteins at the biomolecular force field level can be challenging. Since most Zn2+ 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 Zn2+ 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 Zn2+ 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 Zn2+ site geometries that could guide the further development of biomolecular Zn2+ models.

Molecular Principles of Redox-Coupled Protonation Dynamics in Photosystem II.

Photosystem II (PSII) catalyzes light-driven water oxidization, releasing O2 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 Mn4O5Ca 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 Tyrz (Tyr161) lowers the reaction barrier for the water-mediated proton transfer from a Ca2+-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 Mn4O5Ca cluster lowers the proton transfer barrier through the water ligand sphere of the Mn4O5Ca cluster to Asp61 via a similar ion-pair dissociation process, while the resulting Mn-bound oxo/oxyl species leads to O2 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.

Peroxy Intermediate Drives Carbon Bond Activation in the Dioxygenase AsqJ.

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 (FeIV = 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.

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 (101-106 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 (H2) 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.

Site-specific ubiquitylation and SUMOylation using genetic-code expansion and sortase.

Post-translational modification of proteins with ubiquitin and ubiquitin-like proteins (Ubls) is central to the regulation of eukaryotic cellular processes. Our ability to study the effects of ubiquitylation, however, is limited by the difficulty to prepare homogenously modified proteins in vitro and by the impossibility to selectively trigger specific ubiquitylation events in living cells. Here we combine genetic-code expansion, bioorthogonal Staudinger reduction and sortase-mediated transpeptidation to develop a general tool to ubiquitylate proteins in an inducible fashion. The generated ubiquitin conjugates display a native isopeptide bond and bear two point mutations in the ubiquitin C terminus that confer resistance toward deubiquitinases. Nevertheless, physiological integrity of sortase-generated diubiquitins in decoding cellular functions via recognition by ubiquitin-binding domains is retained. Our approach allows the site-specific attachment of Ubls to nonrefoldable, multidomain proteins and enables inducible and ubiquitin-ligase-independent ubiquitylation of proteins in mammalian cells, providing a powerful tool to dissect the biological functions of ubiquitylation with temporal control.

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, QH2, 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 QH2 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.

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 Tyr416 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.

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.

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.

Reciprocal Coupling in Chemically Fueled Assembly: A Reaction Cycle Regulates Self-Assembly and Vice Versa.

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.

Symmetry-related proton transfer pathways in respiratory complex I.

Complex I functions as the initial electron acceptor in aerobic respiratory chains of most organisms. This gigantic redox-driven enzyme employs the energy from quinone reduction to pump protons across its complete approximately 200-Å membrane domain, thermodynamically driving synthesis of ATP. Despite recently resolved structures from several species, the molecular mechanism by which complex I catalyzes this long-range proton-coupled electron transfer process, however, still remains unclear. We perform here large-scale classical and quantum molecular simulations to study the function of the proton pump in complex I from Thermus thermophilus The simulations suggest that proton channels are established at symmetry-related locations in four subunits of the membrane domain. The channels open up by formation of quasi one-dimensional water chains that are sensitive to the protonation states of buried residues at structurally conserved broken helix elements. Our combined data provide mechanistic insight into long-range coupling effects and predictions for site-directed mutagenesis experiments.