Tunneling and Nonadiabatic Effects on a Proton-Coupled Electron Transfer Model for the Qo Site in Cytochrome bc1.

Cytochrome bc1 is a fundamental enzyme for cellular respiration and photosynthesis. This dimeric protein complex catalyzes a proton-coupled electron transfer (PCET) from the reduced coenzyme-Q substrate (Q) to a bimetallic iron-sulfur cluster in the Qo active site. Herein, we combine molecular dynamics simulations of the complete cytochrome bc1 protein with electronic-structure calculations of truncated models and a semiclassical tunneling theory to investigate the electron-proton adiabaticity of the initial reaction catalyzed in the Qo site. After sampling possible orientations between the Q substrate and a histidine side chain that functions as hydrogen acceptor, we find that a truncated model composed by ubiquinol-methyl and imidazole-iron(III)-sulfide captures the expected changes in oxidation and spin states of the electron donor and acceptor. Diabatic electronic surfaces obtained for this model with multiconfigurational wave function calculations demonstrate that this reaction is electronic nonadiabatic, and proton tunneling is faster than mixing of electronic configurations. These results indicate the formalism that should be used to calculate vibronic couplings and kinetic parameters for the initial reaction in the Qo site of cytochrome bc1. This framework for molecular simulation may also be applied to investigate other PCET reactions in the Q-cycle or in various metalloproteins that catalyze proton translocation coupled to redox processes.

Modeling the Hydrolysis of Iron-Sulfur Clusters.

Iron-sulfur (FeS) clusters are essential metal cofactors involved in a wide variety of biological functions. Their catalytic efficiency, biosynthesis, and regulation depend on FeS stability in aqueous solution. Here, molecular modeling is used to investigate the hydrolysis of an oxidized (ferric) mononuclear FeS cluster by bare dissociation and water substitution mechanisms in neutral and acidic solution. First, approximate electronic structure descriptions of FeS reactions by density functional theory are validated against high-level wave function CCSD(T) calculations. Solvation contributions are included by an all-atom model with hybrid quantum chemical/molecular mechanical (QM/MM) potentials and enhanced sampling molecular dynamics simulations. The free energy profile obtained for FeS cluster hydrolysis indicates that the hybrid functional M06 together with an implicit solvent correction capture the most important aspects of FeS cluster reactivity in aqueous solution. Then, 20 reaction channels leading to two consecutive Fe-S bond ruptures were explored with this calibrated model. For all protonation states, nucleophilic substitution with concerted bond breaking and forming to iron is the preferred mechanism, both kinetic and thermodynamically. In neutral solution, proton transfer from water to the sulfur leaving group is also concerted. Dissociative reactions show higher barriers and will not be relevant for FeS reactivity when exposed to solvent. These hydrolysis mechanisms may help to explain the stability and catalytic mechanisms of FeS clusters of multiple sizes and proteins.