Search

The electron-conducting circuitry of life represents an as-yet untapped resource of exquisite, nanoscale biomolecular engineering. Here, we report the characterization and structure of a de novo diheme "maquette" protein, 4D2, which we subsequently use to create an expanded, modular platform for heme protein design. A well-folded monoheme variant was created by computational redesign, which was then utilized for the experimental validation of continuum electrostatic redox potential calculations. This demonstrates how fundamental biophysical properties can be predicted and fine-tuned. 4D2 was then extended into a tetraheme helical bundle, representing a 7 nm molecular wire. Despite a molecular weight of only 24 kDa, electron cryomicroscopy illustrated a remarkable level of detail, indicating the positioning of the secondary structure and the heme cofactors. This robust, expressible, highly thermostable and readily designable modular platform presents a valuable resource for redox protein design and the future construction of artificial electron-conducting circuitry.

In redox metalloenzymes, the process of electron transfer often involves the concerted movement of a proton. These processes are referred to as proton-coupled electron transfer, and they underpin a wide variety of biological processes, including respiration, energy conversion, photosynthesis, and metalloenzyme catalysis. The mechanisms of proton delivery are incompletely understood, in part due to an absence of information on exact proton locations and hydrogen bonding structures in a bona fide metalloenzyme proton pathway. Here, we present a 2.1-Å neutron crystal structure of the complex formed between a redox metalloenzyme (ascorbate peroxidase) and its reducing substrate (ascorbate). In the neutron structure of the complex, the protonation states of the electron/proton donor (ascorbate) and all of the residues involved in the electron/proton transfer pathway are directly observed. This information sheds light on possible proton movements during heme-catalyzed oxygen activation, as well as on ascorbate oxidation., Competing Interests: The authors declare no competing interest.

Protein dynamics have controversially been proposed to be at the heart of enzyme catalysis, but identification and analysis of dynamical effects in enzyme-catalyzed reactions have proved very challenging. Here, we tackle this question by comparing an enzyme with its heavy ((15)N, (13)C, (2)H substituted) counterpart, providing a subtle probe of dynamics. The crucial hydride transfer step of the reaction (the chemical step) occurs more slowly in the heavy enzyme. A combination of experimental results, quantum mechanics/molecular mechanics simulations, and theoretical analyses identify the origins of the observed differences in reactivity. The generally slightly slower reaction in the heavy enzyme reflects differences in environmental coupling to the hydride transfer step. Importantly, the barrier and contribution of quantum tunneling are not affected, indicating no significant role for "promoting motions" in driving tunneling or modulating the barrier. The chemical step is slower in the heavy enzyme because protein motions coupled to the reaction coordinate are slower. The fact that the heavy enzyme is only slightly less active than its light counterpart shows that protein dynamics have a small, but measurable, effect on the chemical reaction rate.

Cytochrome P450 enzymes play key roles in the metabolism of the majority of drugs. Improved models for prediction of likely metabolites will contribute to drug development. In this work, two possible metabolic routes (aromatic carbon oxidation and O-demethylation) of dextromethorphan are compared using molecular dynamics (MD) simulations and density functional theory (DFT). The DFT results on a small active site model suggest that both reactions might occur competitively. Docking and MD studies of dextromethorphan in the active site of P450 2D6 show that the dextromethorphan is located close to heme oxygen in a geometry apparently consistent with competitive metabolism. In contrast, calculations of the reaction path in a large protein model [using a hybrid quantum mechanical-molecular mechanics (QM/MM) method] show a very strong preference for O-demethylation, in accordance with experimental results. The aromatic carbon oxidation reaction is predicted to have a high activation energy, due to the active site preventing formation of a favorable transition-state structure. Hence, the QM/MM calculations demonstrate a crucial role of many active site residues in determining reactivity of dextromethorphan in P450 2D6. Beyond substrate binding orientation and reactivity of Compound I, successful metabolite predictions must take into account the detailed mechanism of oxidation in the protein. These results demonstrate the potential of QM/MM methods to investigate specificity in drug metabolism.