Search

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.

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

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.

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

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.

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.

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

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.

Quantum chemical (QC) calculations provide a basis for deriving a microscopic understanding of enzymes and photobiological systems. Here we describe how QC models can be used to explore the electronic structure, dynamics, and energetics of biomolecules. We introduce the hybrid quantum mechanics/classical mechanics (QM/MM) approach, where a quantum mechanically described system of interest is embedded in a classically described force field representation of the biochemical surroundings. We also discuss the QM cluster model approach, as well as embedding theories, that provide complementary methodologies to model quantum mechanical effects in biomolecules. The chapter also provides some practical guides for building quantum biochemical models using the quinone reduction catalysis in respiratory complex I and a model reaction in solution as examples.

Complex I couples the free energy released from quinone (Q) reduction to pump protons across the biological membrane in the respiratory chains of mitochondria and many bacteria. The Q reduction site is separated by a large distance from the proton-pumping membrane domain. To address the molecular mechanism of this long-range proton-electron coupling, we perform here full atomistic molecular dynamics simulations, free energy calculations, and continuum electrostatics calculations on complex I from Thermus thermophilus We show that the dynamics of Q is redox-state-dependent, and that quinol, QH 2 , moves out of its reduction site and into a site in the Q tunnel that is occupied by a Q analog in a crystal structure of Yarrowia lipolytica We also identify a second Q-binding site near the opening of the Q tunnel in the membrane domain, where the Q headgroup forms strong interactions with a cluster of aromatic and charged residues, while the Q tail resides in the lipid membrane. We estimate the effective diffusion coefficient of Q in the tunnel, and in turn the characteristic time for Q to reach the active site and for QH 2 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., Competing Interests: The authors declare no conflict of interest., (Copyright © 2018 the Author(s). Published by PNAS.)

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

In Photosystem II (PSII), the Mn 4 CaO 5 -cluster of the active site advances through five sequential oxidation states (S 0 to S 4 ) before water is oxidized and O 2 is generated. Here, we have studied the transition between the low spin (LS) and high spin (HS) configurations of S 2 using EPR spectroscopy, quantum chemical calculations using Density Functional Theory (DFT), and time-resolved UV-visible absorption spectroscopy. The EPR experiments show that the equilibrium between S 2 LS and S 2 HS is pH dependent, with a pK a  ≈ 8.3 (n ≈ 4) for the native Mn 4 CaO 5 and pK a  ≈ 7.5 (n ≈ 1) for Mn 4 SrO 5 . The DFT results suggest that exchanging Ca with Sr modifies the electronic structure of several titratable groups within the active site, including groups that are not direct ligands to Ca/Sr, e.g., W1/W2, Asp61, His332 and His337. This is consistent with the complex modification of the pK a upon the Ca/Sr exchange. EPR also showed that NH 3 addition reversed the effect of high pH, NH 3 -S 2 LS being present at all pH values studied. Absorption spectroscopy indicates that NH 3 is no longer bound in the S 3 Tyr Z state, consistent with EPR data showing minor or no NH 3 -induced modification of S 3 and S 0 . In both Ca-PSII and Sr-PSII, S 2 HS was capable of advancing to S 3 at low temperature (198 K). This is an experimental demonstration that the S 2 LS is formed first and advances to S 3 via the S 2 HS state without detectable intermediates. We discuss the nature of the changes occurring in the S 2 LS to S 2 HS transition which allow the S 2 HS to S 3 transition to occur below 200 K. This work also provides a protocol for generating S 3 in concentrated samples without the need for saturating flashes., (Copyright © 2018 Elsevier B.V. All rights reserved.)

The respiratory complex I is an enzyme responsible for the conversion of chemical energy into an electrochemical proton motive force across the membrane. Despite extensive studies, the mechanism by which the activity of this enormous, ca. 1 MDa, redox-coupled proton pump is regulated still remains unclear. Recent structural studies (Zhu et al., Nature 2016; Fiedorczuk et al., Nature 2016) resolved complex I in different conformations connected to the active-to-deactive (A/D) transition that regulate complex I activity in several species. Based on anisotropic network models (ANM) and principal component analysis (PCA), we identify here transitions between experimentally resolved structures of the mammalian complex I as low-frequency collective motions of the enzyme, highlighting similarities and differences between the bacterial and mammalian enzymes. Despite the reduced complexity of the smaller bacterial enzyme, our results suggest that the global dynamics of complex I is overall conserved. We further probe how the supernumerary subunits could be involved in the modulation of the A/D-transition, and show that in particular the 42 kDa and B13 subunits affect the global motions of the mammalian enzyme., (Copyright © 2018. Published by Elsevier B.V.)

Respiratory complex I (NADH:ubiquinone oxidoreductase), one of the largest membrane-bound enzymes in mammalian cells, powers ATP synthesis by using the energy from electron transfer from NADH to ubiquinone-10 to drive protons across the energy-transducing mitochondrial inner membrane. Ubiquinone-10 is extremely hydrophobic, but in complex I the binding site for its redox-active quinone headgroup is ∼20 Å above the membrane surface. Structural data suggest it accesses the site by a narrow channel, long enough to accommodate almost all of its ∼50-Å isoprenoid chain. However, how ubiquinone/ubiquinol exchange occurs on catalytically relevant timescales, and whether binding/dissociation events are involved in coupling electron transfer to proton translocation, are unknown. Here, we use proteoliposomes containing complex I, together with a quinol oxidase, to determine the kinetics of complex I catalysis with ubiquinones of varying isoprenoid chain length, from 1 to 10 units. We interpret our results using structural data, which show the hydrophobic channel is interrupted by a highly charged region at isoprenoids 4-7. We demonstrate that ubiquinol-10 dissociation is not rate determining and deduce that ubiquinone-10 has both the highest binding affinity and the fastest binding rate. We propose that the charged region and chain directionality assist product dissociation, and that isoprenoid stepping ensures short transit times. These properties of the channel do not benefit the exhange of short-chain quinones, for which product dissociation may become rate limiting. Thus, we discuss how the long channel does not hinder catalysis under physiological conditions and the possible roles of ubiquinone/ubiquinol binding/dissociation in energy conversion., Competing Interests: The authors declare no conflict of interest.

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., Competing Interests: The authors declare no conflict of interest.

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

Cytochrome c oxidase (CcO) is a redox-driven proton pump that powers aerobic respiratory chains. We show here by multi-scale molecular simulations that a protonated water cluster near the active site is likely to serve as the transient proton-loading site (PLS) that stores a proton during the pumping process. The pKa of this water cluster is sensitive to the redox states of the enzyme, showing distinct similarities to other energy converting proton pumps., (© 2016 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim.)

Rational design of light-capturing properties requires understanding the molecular and electronic structure of chromophores in their native chemical or biological environment. We employ here large-scale quantum chemical calculations to study the light-capturing properties of retinal in recently designed human cellular retinol binding protein II (hCRBPII) variants (Wang et al. Science, 2012, 338, 1340-1343). Our calculations show that these proteins absorb across a large part of the visible spectrum by combined polarization and electrostatic effects. These effects stabilize the ground or excited state energy levels of the retinal by perturbing the Schiff-base or β-ionone moieties of the chromophore, which in turn modulates the amount of charge transfer within the molecule. Based on the predicted tuning principles, we design putative in silico mutations that further shift the absorption properties of retinal in hCRBPII towards the ultraviolet and infrared regions of the spectrum., (© 2016 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim.)

Photosystem II (PSII) catalyzes light-driven water splitting in nature and is the key enzyme for energy input into the biosphere. Important details of its mechanism are not well understood. In order to understand the mechanism of water splitting, we perform here large-scale density functional theory (DFT) calculations on the active site of PSII in different oxidation, spin and ligand states. Prior to formation of the O-O bond, we find that all manganese atoms are oxidized to Mn(IV) in the S3 state, consistent with earlier studies. We find here, however, that the formation of the S3 state is coupled to the movement of a calcium-bound hydroxide (W3) from the Ca to a Mn (Mn1 or Mn4) in a process that is triggered by the formation of a tyrosyl radical (Tyr-161) and its protonated base, His-190. We find that subsequent oxidation and deprotonation of this hydroxide on Mn1 result in formation of an oxyl-radical that can exergonically couple with one of the oxo-bridges (O5), forming an O-O bond. When O(2) leaves the active site, a second Ca-bound water molecule reorients to bridge the gap between the manganese ions Mn1 and Mn4, forming a new oxo-bridge for the next reaction cycle. Our findings are consistent with experimental data, and suggest that the calcium ion may control substrate water access to the water oxidation sites., (Copyright © 2016. Published by Elsevier B.V.)

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

Hsp90 is a molecular chaperone involved in the activation of numerous client proteins, including many kinases. The most stringent kinase client is the oncogenic kinase v-Src. To elucidate how Hsp90 chaperones kinases, we reconstituted v-Src kinase chaperoning in vitro and show that its activation is ATP-dependent, with the cochaperone Cdc37 increasing the efficiency. Consistent with in vivo results, we find that Hsp90 does not influence the almost identical c-Src kinase. To explain these findings, we designed Src kinase chimeras that gradually transform c-Src into v-Src and show that their Hsp90 dependence correlates with compactness and folding cooperativity. Molecular dynamics simulations and hydrogen/deuterium exchange of Hsp90-dependent Src kinase variants further reveal increased transitions between inactive and active states and exposure of specific kinase regions. Thus, Hsp90 shifts an ensemble of conformations of v-Src toward high activity states that would otherwise be metastable and poorly populated.

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

This chapter introduces the theory and applications of commonly used methods of electronic structure calculation, with particular emphasis on methods applicable for modelling biomolecular systems. This chapter is sectioned as follows. We start by presenting ab initio methods, followed by a treatment of density functional theory (DFT) and some recent advances in semi-empirical methods. Treatment of excited states as well as basis sets are also presented.

To understand how signaling proteins function, it is crucial to know the time-ordered sequence of events that lead to the signaling state. We recently developed on the BioCARS 14-IDB beamline at the Advanced Photon Source the infrastructure required to characterize structural changes in protein crystals with near-atomic spatial resolution and 150-ps time resolution, and have used this capability to track the reversible photocycle of photoactive yellow protein (PYP) following trans-to-cis photoisomerization of its p-coumaric acid (pCA) chromophore over 10 decades of time. The first of four major intermediates characterized in this study is highly contorted, with the pCA carbonyl rotated nearly 90° out of the plane of the phenolate. A hydrogen bond between the pCA carbonyl and the Cys69 backbone constrains the chromophore in this unusual twisted conformation. Density functional theory calculations confirm that this structure is chemically plausible and corresponds to a strained cis intermediate. This unique structure is short-lived (∼600 ps), has not been observed in prior cryocrystallography experiments, and is the progenitor of intermediates characterized in previous nanosecond time-resolved Laue crystallography studies. The structural transitions unveiled during the PYP photocycle include trans/cis isomerization, the breaking and making of hydrogen bonds, formation/relaxation of strain, and gated water penetration into the interior of the protein. This mechanistically detailed, near-atomic resolution description of the complete PYP photocycle provides a framework for understanding signal transduction in proteins, and for assessing and validating theoretical/computational approaches in protein biophysics.

Heme-copper oxidases (HCOs) are terminal electron acceptors in aerobic respiration. They catalyze the reduction of molecular oxygen to water with concurrent pumping of protons across the mitochondrial and bacterial membranes. Protons required for oxygen reduction chemistry and pumping are transferred through proton uptake channels. Recently, the crystal structure of the first C-type member of the HCO superfamily was resolved [Buschmann et al. Science 329 (2010) 327-330], but crystallographic water molecules could not be identified. Here we have used molecular dynamics (MD) simulations, continuum electrostatic approaches, and quantum chemical cluster calculations to identify proton transfer pathways in cytochrome cbb(3). In MD simulations we observe formation of stable water chains that connect the highly conserved Glu323 residue on the proximal side of heme b(3) both with the N- and the P-sides of the membrane. We propose that such pathways could be utilized for redox-coupled proton pumping in the C-type oxidases. Electrostatics and quantum chemical calculations suggest an increased proton affinity of Glu323 upon reduction of high-spin heme b(3). Protonation of Glu323 provides a mechanism to tune the redox potential of heme b(3) with possible implications for proton pumping., (Copyright © 2011. Published by Elsevier B.V.)

Cytochrome c oxidase (CcO) is the terminal enzyme of the respiratory chain. By reducing oxygen to water, it generates a proton gradient across the mitochondrial or bacterial membrane. Recently, two independent X-ray crystallographic studies ((Aoyama et al. Proc. Natl. Acad. Sci. USA 106 (2009) 2165-2169) and (Koepke et al. Biochim. Biophys. Acta 1787 (2009) 635-645)), suggested that a peroxide dianion might be bound to the active site of oxidized CcO. We have investigated this hypothesis by combining quantum chemical calculations with a re-refinement of the X-ray crystallographic data and optical spectroscopic measurements. Our data suggest that dianionic peroxide, superoxide, and dioxygen all form a similar superoxide species when inserted into a fully oxidized ferric/cupric binuclear site (BNC). We argue that stable peroxides are unlikely to be confined within the oxidized BNC since that would be expected to lead to bond splitting and formation of the catalytic P intermediate. Somewhat surprisingly, we find that binding of dioxygen to the oxidized binuclear site is weakly exergonic, and hence, the observed structure might have resulted from dioxygen itself or from superoxide generated from O(2) by the X-ray beam. We show that the presence of O(2) is consistent with the X-ray data. We also discuss how other structures, such as a mixture of the aqueous species (H(2)O+OH(-) and H(2)O) and chloride fit the experimental data., (Copyright © 2011 Elsevier B.V. All rights reserved.)

The proton-pumping cbb(3)-type cytochrome c oxidases catalyze cell respiration in many pathogenic bacteria. For reasons not yet understood, the apparent dioxygen (O(2)) affinity in these enzymes is very high relative to other members of the heme-copper oxidase (HCO) superfamily. Based on density functional theory (DFT) calculations on intermediates of the oxygen scission reaction in active-site models of cbb(3)- and aa(3)-type oxidases, we find that a transient peroxy intermediate (I(P), Fe[III]-OOH(-)) is ~6kcal/mol more stable in the former case, resulting in more efficient kinetic trapping of dioxygen and hence in a higher apparent oxygen affinity. The major molecular basis for this stabilization is a glutamate residue, polarizing the proximal histidine ligand of heme b(3) in the active site., (Published by Elsevier B.V.)

Cellular respiration is driven by cytochrome c oxidase (CcO), which reduces oxygen to water and couples the released energy to proton pumping across the mitochondrial or bacterial membrane. Proton pumping in CcO involves proton transfer from the negatively charged side of the membrane to a transient proton-loading or pump site (PLS), before it is ejected to the opposite side. Although many details of the reaction mechanism are known, the exact location of the PLS has remained elusive. We report here results from combined classical molecular dynamics simulations and continuum electrostatic calculations, which show that the hydrogen-bonded system around the A-propionate of heme a₃ dissociates reversibly upon reduction of heme a. The dissociation increases the pK(a) value of the propionate to a value above ~9, making it accessible for redox-state dependent protonation. The redox state of heme a is of key importance in controlling proton leaks by polarizing the PLS both statically and dynamically. These findings suggest that the propionate region of heme a₃ fulfills the criteria of the pump site in the proton translocation mechanism of CcO., (Copyright © 2010 Elsevier B.V. All rights reserved.)

Cytochrome c oxidase (CcO) is the terminal enzyme of the respiratory chain that catalyzes respiratory reduction of dioxygen (O(2)) to water in all eukaryotes and many aerobic bacteria. CcO, and its homologs among the heme-copper oxidases, has an active site composed of an oxygen-binding heme and a copper center in the vicinity, plus another heme group that donates electrons to this site. In most oxidoreduction enzymes, electron transfer (eT) takes place by quantum-mechanical electron tunneling. Here we show by independent molecular dynamics and quantum-chemical methods that the heme-heme eT in CcO differs from the majority of cases in having an exceptionally low reorganization energy. We show that the rate of interheme eT in CcO may nevertheless be predicted by the Moser-Dutton equation if reinterpreted as the average of the eT rates between all individual atoms of the donor and acceptor weighed by the respective packing densities between them. We argue that this modification may be necessary at short donor/acceptor distances comparable to the donor/acceptor radii.

Cytochrome cbb3 is a distinct member of the superfamily of respiratory heme-copper oxidases, and is responsible for driving the respiratory chain in many pathogenic bacteria. Like the canonical heme-copper oxidases, cytochrome cbb3 reduces oxygen to water and couples the released energy to pump protons across the bacterial membrane. Homology modeling and recent electron paramagnetic resonance (EPR) studies on wild type and a mutant cbb3 enzyme [V. Rauhamäki et al. J. Biol. Chem. 284 (2009) 11301-11308] have led us to perform high-level quantum chemical calculations on the active site. These calculations bring molecular insight into the unique hydrogen bonding between the proximal histidine ligand of heme b3 and a conserved glutamate, and indicate that the catalytic mechanism involves redox-coupled proton transfer between these residues. The calculated spin densities give insight in the difference in EPR spectra for the wild type and a recently studied E383Q-mutant cbb3-enzyme. Furthermore, we show that the redox-coupled proton movement in the proximal cavity of cbb3-enzymes contributes to the low redox potential of heme b3, and suggest its potential implications for the high apparent oxygen affinity of these enzymes., (Copyright (c) 2010 Elsevier B.V. All rights reserved.)

Cytochrome c oxidase (CcO) is the terminal enzyme of aerobic respiration. The energy released from the reduction of molecular oxygen to water is used to pump protons across the mitochondrial or bacterial membrane. The pump function introduces a mechanistic requirement of a valve that prevents protons from flowing backwards during the process. It was recently found that Glu-242, a key amino acid in transferring protons to be pumped across the membrane and to the site of oxygen reduction, fulfils the function of such a valve by preventing simultaneous contact to the pump site and to the proton-conducting D-channel (Kaila V.R.I. et al. Proc. Natl. Acad. Sci. USA 105, 2008). Here we have incorporated the valve model into the framework of the reaction mechanism. The function of the Glu valve is studied by exploring how the redox state of the surrounding metal centers, dielectric effects, and membrane potential, affects the energetics and leaks of this valve. Parallels are drawn between the dynamics of Glu-242 and the long-standing obscure difference between the metastable O(H) and stable O states of the binuclear center. Our model provides a suggestion for why reduction of the former state is coupled to proton translocation while reduction of the latter is not.

The CuB metal center is at the core of the active site of the heme-copper oxidases, comprising a copper atom ligating three histidine residues one of which is covalently bonded to a tyrosine residue. Using quantum chemical methodology, we have studied the CuB site in several redox and ligand states proposed to be intermediates of the catalytic cycle. The importance of the His-Tyr crosslink was investigated by comparing energetics, charge, and spin distributions between systems with and without the crosslink. The His-Tyr bond was shown to decrease the proton affinity and increase the electron affinity of both Tyr-244 and the copper. A previously unnoticed internal electronic equilibrium between the copper atom and the tyrosine was observed, which seems to be coupled to the unique structure of the system. In certain states the copper and Tyr-244 compete for the unpaired electron, the localization of which is determined by the oxygenous ligand of the copper. This electronic equilibrium was found to be sensitive to the presence of a positive charge 10 A away from the center, simulating the effect of Lys-319 in the K-pathway of proton transfer. The combined results provide an explanation for why the heme-copper oxidases need two pathways of proton uptake, and why the K-pathway is active only in the second half of the reaction cycle.

The cytochrome c oxidases (CcO), which are responsible for most O(2) consumption in biology, are also redox-linked proton pumps that effectively convert the free energy of O(2) reduction to an electrochemical proton gradient across mitochondrial and bacterial membranes. Recently, time-resolved measurements have elucidated the sequence of events in proton translocation, and shed light on the underlying molecular mechanisms. One crucial property of the proton pump mechanism has received less attention, viz. how proton leaks are avoided. Here, we will analyse this topic and demonstrate how the key proton-carrying residue Glu-242 (numbering according to the sequence of subunit I of bovine heart CcO) functions as a valve that has the effect of minimising back-leakage of the pumped proton.

Aerobic life is based on a molecular machinery that utilizes oxygen as a terminal electron sink. The membrane-bound cytochrome c oxidase (CcO) catalyzes the reduction of oxygen to water in mitochondria and many bacteria. The energy released in this reaction is conserved by pumping protons across the mitochondrial or bacterial membrane, creating an electrochemical proton gradient that drives production of ATP. A crucial question is how the protons pumped by CcO are prevented from flowing backwards during the process. Here, we show by molecular dynamics simulations that the conserved glutamic acid 242 near the active site of CcO undergoes a protonation state-dependent conformational change, which provides a valve in the pumping mechanism. The valve ensures that at any point in time, the proton pathway across the membrane is effectively discontinuous, thereby preventing thermodynamically favorable proton back-leakage while maintaining an overall high efficiency of proton translocation. Suppression of proton leakage is particularly important in mitochondria under physiological conditions, where production of ATP takes place in the presence of a high electrochemical proton gradient.

In many cytochrome c oxidases glutamic acid 242 is required for proton transfer to the binuclear heme a(3)/Cu(B) site, and for proton pumping. When present, the side chain of Glu-242 is orientated "down" towards the proton-transferring D-pathway in all available crystal structures. A nonpolar cavity "above" Glu-242 is empty in these structures. Yet, proton transfer from Glu-242 to the binuclear site, and for proton-pumping, is well established, and the cavity has been proposed to at least transiently contain water molecules that would mediate proton transfer. Such proton transfer has been proposed to require isomerisation of the Glu-242 side chain into an "up" position pointing towards the cavity. Here, we have explored the molecular dynamics of the protonated Glu-242 side chain. We find that the "up" position is preferred energetically when the cavity contains four water molecules, but the "down" position is favoured with less water. We conclude that the cavity might be deficient in water in the crystal structures, possibly reflecting the "resting" state of the enzyme, and that the "up/down" equilibrium of Glu-242 may be coupled to the presence of active-site water molecules produced by O(2) reduction.