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.
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.
Retinal is the light-absorbing biochromophore responsible for the activation of vision pigments and light-driven ion pumps. Nature has evolved molecular tuning mechanisms that significantly shift the ...optical properties of the retinal pigments to enable their absorption of visible light. Using large-scale quantum chemical calculations at the density functional theory level combined with frozen density embedding theory, we show here how the protein environment of vision pigments tunes the absorption of retinal by electrostatically dominated interactions between the chromophore and the surrounding protein residues. The calculations accurately reproduce the experimental absorption maxima of rhodopsin and the red, green, and blue color pigments. We further identify key interactions responsible for the color-shifting effects by mutating the rhodopsin structure in silico, and we find that deprotonation of the retinyl is likely to be responsible for the blue-shifted absorption in the blue cone vision pigment.
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.
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.
The FKBP51 inhibitor SAFit1 binds a minor conformation of FKBP51 with an outward conformation of Phe67 in the unbound state, suggesting SAFit1 binds to FKBP51 using conformational selection, rather than induced fit. Stabilizing Phe67 in this dynamic outward conformation involves interactions with FKBP51‐specific residues, explaining its selectivity for FKBP51 vs. FKBP52. This shows the importance of protein dynamics in designing selective inhibitors.
Proton transfer (pT) reactions in biochemical processes are often mediated by chains of hydrogen-bonded water molecules. We use hybrid density functional calculations to study pT along quasi ...one-dimensional water arrays that connect an imidazolium-imidazole proton donor-acceptor pair. We characterize the structures of intermediates and transition states, the energetics, and the dynamics of the pT reactions, including vibrational contributions to kinetic isotope effects. In molecular dynamics simulations of pT transition paths, we find that for short water chains with four water molecules, the pT reactions are semi-concerted. The formation of a high-energy hydronium intermediate next to the proton-donating group is avoided by a simultaneous transfer of a proton from the donor to the first water molecule, and from the first water molecule into the water chain. Lowering the dielectric constant of the environment and increasing the water chain length both reduce the barrier for pT. We study the effect of the driving force on the energetics of the pT reaction by changing the proton affinity of the donor and acceptor groups through halogen and methyl substitutions. We find that the barrier of the pT reaction depends linearly on the proton affinity of the donor but is nearly independent of the proton affinity of the acceptor, corresponding to Brønsted slopes of one and zero, respectively.
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.
This benchmarking study reveals key performance differences and limitations of widely‐used non‐polarizable Zn2+ models in metalloproteins. The characterization of these models regarding their preferred coordination geometry, preferred ligating atoms, and ligand binding stability during molecular dynamics simulations allowed us to suggest suitable simulation conditions for varying Zn2+ site geometries, which can guide further development of these models.
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.
•Global motions of mammalian complex I are analyzed with network models.•Twisting and scissoring motions of membrane and hydrophilic domains are identified.•The global modes have a high overlap with cryo-EM derived structural classes.•The bacterial complex I have motions similar to the mammalian enzyme.•Supernumerary subunits 42 kDa and B13 modulate the global motions of complex I.
Proton-coupled electron transfer (PCET) processes are elementary chemical reactions involved in a broad range of radical and redox reactions. Elucidating fundamental PCET reaction mechanisms are thus ...of central importance for chemical and biochemical research. Here we use quantum chemical density functional theory (DFT), time-dependent density functional theory (TDDFT), and the algebraic diagrammatic-construction through second-order (ADC(2)) to study the mechanism, thermodynamic driving force effects, and reaction barriers of both ground state proton transfer (pT) and photoinduced proton-coupled electron transfer (PCET) between nitrosylated phenyl-phenol compounds and hydrogen-bonded t-butylamine as an external base. We show that the obtained reaction barriers for the ground state pT reactions depend linearly on the thermodynamic driving force, with a Brønsted slope of 1 or 0. Photoexcitation leads to a PCET reaction, for which we find that the excited state reaction barrier depends on the thermodynamic driving force with a Brønsted slope of 1/2. To support the mechanistic picture arising from the static potential energy surfaces, we perform additional molecular dynamics simulations on the excited state energy surface, in which we observe a spontaneous PCET between the donor and the acceptor groups. Our findings suggest that a Brønsted analysis may distinguish the ground state pT and excited state PCET processes.
Mitochondrial proton-pumping NADH:ubiquinone oxidoreductase (respiratory complex I) comprises more than 40 polypeptides and contains eight canonical FeS clusters. The integration of subunits and ...insertion of cofactors into the nascent complex is a complicated multistep process that is aided by assembly factors. We show that the accessory NUMM subunit of complex I (human NDUFS6) harbors a Zn-binding site and resolve its position by X-ray crystallography. Chromosomal deletion of the NUMM gene or mutation of Zn-binding residues blocked a late step of complex I assembly. An accumulating assembly intermediate lacked accessory subunit N7BM (NDUFA12), whereas a paralog of this subunit, the assembly factor N7BML (NDUFAF2), was found firmly bound instead. EPR spectroscopic analysis and metal content determination after chromatographic purification of the assembly intermediate showed that NUMM is required for insertion or stabilization of FeS cluster N4.
Significance Respiratory complex I is the largest membrane protein complex in mitochondria and has a central function in energy metabolism. Numerous human diseases are linked with complex I dysfunction or assembly defects. The concerted assembly of more than 40 subunits and the insertion of cofactors is aided by specific chaperones. In addition to eight FeS clusters, complex I comprises a Zn-binding site of unknown function. Combining X-ray structural analysis of complex I crystals with quantum chemical modeling and proteomic and spectroscopic analysis of a purified assembly intermediate, we show that accessory subunit NUMM (human ortholog NDUFS6) binds Zn at the interface of two functional modules of the enzyme complex and is required for a specific step of complex I biogenesis.