Conformational Heterogeneity and the Affinity of Substrate Molecular Recognition by Cytochrome P450cam

Unnatural Amino Acids for Biological Spectroscopy and Microscopy

Due to advances in methods for site-specific incorporation of unnatural amino acids (UAAs) into proteins, a large number of UAAs with tailored chemical and/or physical properties have been developed and used in a wide array of biological applications. In particular, UAAs with specific spectroscopic characteristics can be used as external reporters to produce additional signals, hence increasing the information content obtainable in protein spectroscopic and/or imaging measurements. In this Review, we summarize the progress in the past two decades in the development of such UAAs and their applications in biological spectroscopy and microscopy, with a focus on UAAs that can be used as site-specific vibrational, fluorescence, electron paramagnetic resonance (EPR), or nuclear magnetic resonance (NMR) probes. Wherever applicable, we also discuss future directions.

The Role H-Bonding and Supramolecular Structures in Homogeneous and Enzymatic Catalysis

The article analyzes the role of hydrogen bonds and supramolecular structures in enzyme catalysis and model systems. Hydrogen bonds play a crucial role in many enzymatic reactions. However, scientists have only recently attempted to harness the power of hydrogen bonds in homogeneous catalytic systems. One of the newest directions is associated with attempts to control the properties of catalysts by influencing the "second coordination sphere" of metal complexes. The role H-bonding, and the building of stable supramolecular nanostructures due to intermolecular H-bonds, based on catalytic active heteroligand iron (Fe) or nickel (Ni) complexes formed during hydrocarbon oxidations were assessed via the AFM (Atomic-force microscopy) method, which was proposed and applied by authors of this manuscript. Th is article also discusses the roles of hydrogen bonds and supramolecular structures in oxidation reactions catalyzed by heteroligand Ni and Fe complexes, which are not only effective homogeneous catalysts but also structural and functional models of Oxygenases.

3site Multisubstrate-Bound State of Cytochrome P450cam

We identified a multisubstrate-bound state, hereby referred as a 3site state, in cytochrome P450cam via integrating molecular dynamics simulation with nuclear magnetic resonance (NMR) pseudocontact shift measurements. The 3site state is a result of simultaneous binding of three camphor molecules in three locations around P450cam: (a) in a well-established "catalytic" site near heme, (b) in a kink-separated "waiting" site along channel-1, and (c) in a previously reported "allosteric" site at E, F, G, and H helical junctions. These three spatially distinct binding modes in the 3site state mutually communicate with each other via homotropic allostery and act cooperatively to render P450cam functional. The 3site state shows a significantly superior fit with NMR pseudo contact shift (PCS) data with a Q-score of 0.045 than previously known bound states and consists of D251 free of salt-bridges with K178 and R186, rendering the enzyme functionally primed. To date, none of the reported cocomplex of P450cam with its redox partner putidaredoxin (pdx) has been able to match solution NMR data and controversial pdx-induced opening of P450cam's channel-1 remains a matter of recurrent discourse. In this regard, inclusion of pdx to the 3site state is able to perfectly fit the NMR PCS measurement with a Q-score of 0.08 and disfavors the pdx-induced opening of channel-1, reconciling previously unexplained remarkably fast hydroxylation kinetics with a koff of 10.2 s-1. Together, our findings hint that previous experimental observations may have inadvertently captured the 3site state as an in vitro solution state, instead of the catalytic state alone, and provided a distinct departure from the conventional understanding of cytochrome P450.

The Dual Function of PhOH Included in the Coordination Sphere of the Nickel Complexes in the Processes of Oxidation with Dioxygen

The role of ligands in the regulation of the catalytic activity of Ni-complexes (Ni(acac)₂) in green process-selective ethylbenzene oxidation with O₂ into α-phenyl ethyl hydroperoxide is considered in this article. The dual function of phenol (PhOH) included in the coordination sphere of the nickel complex as an antioxidant or catalyst depends on the ligand environment of the metal. The role of intermolecular H-bonds and supramolecular structures (AFM method) in the mechanisms of selective catalysis by nickel complexes in chemical and biological oxidation reactions is analyzed.

A deep autoencoder framework for discovery of metastable ensembles in biomacromolecules

Biomacromolecules manifest dynamic conformational fluctuation and involve mutual interconversion among metastable states. A robust mapping of their conformational landscape often requires the low-dimensional projection of the conformational ensemble along optimized collective variables (CVs). However, the traditional choice for the CV is often limited by user-intuition and prior knowledge about the system, and this lacks a rigorous assessment of their optimality over other candidate CVs. To address this issue, we propose an approach in which we first choose the possible combinations of inter-residue Cα-distances within a given macromolecule as a set of input CVs. Subsequently, we derive a non-linear combination of latent space embedded CVs via auto-encoding the unbiased molecular dynamics simulation trajectories within the framework of the feed-forward neural network. We demonstrate the ability of the derived latent space variables in elucidating the conformational landscape in four hierarchically complex systems. The latent space CVs identify key metastable states of a bead-in-a-spring polymer. The combination of the adopted dimensional reduction technique with a Markov state model, built on the derived latent space, reveals multiple spatially and kinetically well-resolved metastable conformations for GB1 β-hairpin. A quantitative comparison based on the variational approach-based scoring of the auto-encoder-derived latent space CVs with the ones obtained via independent component analysis (principal component analysis or time-structured independent component analysis) confirms the optimality of the former. As a practical application, the auto-encoder-derived CVs were found to predict the reinforced folding of a Trp-cage mini-protein in aqueous osmolyte solution. Finally, the protocol was able to decipher the conformational heterogeneities involved in a complex metalloenzyme, namely, cytochrome P450.

Role of PhOH and Tyrosine in Selective Oxidation of Hydrocarbons

Earlier, we established that nickel or iron heteroligand complexes, which include PhOH (nickel complexes) or tyrosine residue (nickel or iron complexes), are not only hydrocarbon oxidation catalysts (in the case of PhOH), but also simulate the active centers of enzymes (PhOH, tyrosine). The AFM method established the self-organization of nickel or iron heteroligand complexes, which included tyrosine residue or PhOH, into supramolecular structures on a modified silicon surface. Supramolecular structures were formed as a result of H-bonds and other non-covalent intermolecular interactions and, to a certain extent, reflected the structures involved in the mechanisms of reactions of homogeneous and enzymatic catalysis. Using the AFM method, we obtained evidence at the model level in favor of the involvement of the tyrosine fragment as one of the possible regulatory factors in the functioning of Ni(Fe)ARD dioxygenases or monooxygenases of the family of cytochrome P450. The principles of actions of these oxygenases were used to create highly efficient catalytic systems for the oxidation of hydrocarbons.

Reconciling conformational heterogeneity and substrate recognition in cytochrome P450

Cytochrome P450, the ubiquitous metalloenzyme involved in detoxification of foreign components, has remained one of the most popular systems for substrate-recognition process. However, despite being known for its high substrate specificity, the mechanistic basis of substrate-binding by archetypal system cytochrome P450cam has remained at odds with the contrasting reports of multiple diverse crystallographic structures of its substrate-free form. Here, we address this issue by elucidating the probability of mutual dynamical transition to the other crystallographic pose of cytochrome P450cam and vice versa via unbiased all-atom computer simulation. A robust Markov state model, constructed using adaptively sampled 84-μs-long molecular dynamics simulation trajectories, maps the broad and heterogenous P450cam conformational landscape into five key substates. In particular, the Markov state model identifies an intermediate-assisted dynamic equilibrium between a pair of conformations of P450cam, in which the substrate-recognition sites remain "closed" and "open," respectively. However, the estimate of a significantly higher stationary population of closed conformation, coupled with faster rate of open → closed transition than its reverse process, dictates that the net conformational equilibrium would be swayed in favor of "closed" conformation. Together, the investigation quantitatively infers that although a potential substrate of cytochrome P450cam would, in principle, explore a diverse array of conformations of substrate-free protein, it would mostly encounter a "closed" or solvent-occluded conformation and hence would follow an induced-fit-based recognition process. Overall, the work reconciles multiple precedent crystallographic, spectroscopic investigations and establishes how a statistical elucidation of conformational heterogeneity in protein would provide crucial insights in the mechanism of potential substrate-recognition process.

Kinetic Evidence for an Induced Fit Mechanism in the Binding of the Substrate Camphor by Cytochrome P450cam

Bacterial cytochrome P450 (P450) 101A1 (P450cam) has served as a prototype among the P450 enzymes and has high catalytic activity towards its cognate substrate, camphor. X-ray crystallography and NMR and IR spectroscopy have demonstrated the existence of multiple conformations of many P450s, including P450cam. Kinetic studies have indicated that substrate binding to several P450s is dominated by a conformational selection process, in which the substrate binds an individual conformer(s) of the unliganded enzyme. P450cam was found to differ in that binding of the substrate camphor is dominated by an induced fit mechanism, in which the enzyme binds camphor and then changes conformation, as evidenced by the equivalence of binding eigenvalues observed when varying both camphor and P450cam concentrations. The accessory protein putidaredoxin had no effect on substrate binding. Estimation of the rate of dissociation of the P450cam·camphor complex (15 s-1) and fitting of the data yield a minimal kinetic mechanism in which camphor binds (1.5 × 107 M-1 s-1) and the initial P450cam•camphor complex undergoes a reversible equilibrium (k forward 112 s-1, k reverse 28 s-1) to a final complex. This induced fit mechanism differs from those reported for several mammalian P450s and bacterial P450BM-3, indicative of the diversity of how P450s recognize multiple substrates. However, similar behavior was not observed with the alternate substrates (+)-α-pinene and 2-adamantanone, which probably utilize a conformational selection process.

Low-Frequency Protein Motions Coupled to Catalytic Sites

This review examines low-frequency vibrational modes of proteins and their coupling to enzyme catalytic sites. That protein motions are critical to enzyme function is clear, but the kinds of motions present in proteins and how they are involved in function remain unclear. Several models of enzyme-catalyzed reaction suggest that protein dynamics may be involved in the chemical step of the catalyzed reaction, but the evidence in support of such models is indirect. Spectroscopic studies of low-frequency protein vibrations consistently show that there are underdamped modes of the protein with frequencies in the tens of wavenumbers where overdamped behavior would be expected. Recent studies even show that such underdamped vibrations modulate enzyme active sites. These observations suggest that increasingly sophisticated spectroscopic methods will be able to unravel the link between low-frequency protein vibrations and enzyme function.

AFM Research in Catalysis and Medicine

Background:

In this study, we show that the AFM method not only allows monitoring the morphological changes in biological structures fixed on the surface due to H-bonds, but also makes it possible to study the self-organization of metal complexes by simulating the active center of enzymes due to intermolecular H-bonds into stable nanostructures; the sizes of which are much smaller than the studied biological objects. The possible role of intermolecular hydrogen bonds in the formation of stable supramolecular metal complexes, which are effective catalysts for the oxidation of alkyl arenes to hydroperoxides by molecular oxygen and mimic the selective active sites of enzymes, was first studied by AFM.

Methods and Results:

The formation of supramolecular structures due to intermolecular hydrogen bonds and, possibly, other non-covalent interactions, based on homogenous catalysts and models of active centers enzymes, heteroligand nickel and iron complexes, was proven by AFM-technique. AFM studies of supramolecular structures were carried out using NSG30 cantilever with a radius of curvature of 2 nm, in the tapping mode. To form nanostructures on the surface of a hydrophobic, chemically modified silicon surface as a substrate, the sample was prepared using a spin-coating process from solutions of the nickel and iron complexes. The composition and the structure of the complex Ni2(acac)(OAc)3·NMP·2H2O were determined in earlier works using various methods: mass spectrometry, UV- and IR-spectroscopy, elemental analysis, and polarography. Self-assembly of supramolecular structures is due to intermolecular interactions with a certain coordination of these interactions, which may be a consequence of the properties of the components themselves, the participation of hydrogen bonds and other non-covalent interactions, as well as the balance of the interaction of these components with the surface. Using AFM, approaches have been developed for fixing on the surface and quantifying parameters of cells.

Conclusion:

This study summarizes the authors' achievements in using the atomic force microscopy (AFM) method to study the role of intermolecular hydrogen bonds (and other non-covalent interactions) and supramolecular structures in the mechanisms of catalysis. The data obtained from AFM based on nickel and iron complexes, which are effective catalysts and models of active sites of enzymes, indicate a high probability of the formation of supramolecular structures in real conditions of catalytic oxidation, and can bring us closer to understanding enzymes activity. With a sensitive AFM method, it is possible to observe the self-organization of model systems into stable nanostructures due to H-bonds and possibly other non-covalent interactions, which can be considered as a step towards modeling the active sites of enzymes. Methodical approaches of atomic force microscopy for the study of morphological changes of cells have been developed.

Horizontal 'gene drives' harness indigenous bacteria for bioremediation

Engineering bacteria to clean-up oil spills is rapidly advancing but faces regulatory hurdles and environmental concerns. Here, we develop a new technology to harness indigenous soil microbial communities for bioremediation by flooding local populations with catabolic genes for petroleum hydrocarbon degradation. Overexpressing three enzymes (almA, xylE, p450cam) in Escherichia coli led to degradation of 60-99% of target hydrocarbon substrates. Mating experiments, fluorescence microscopy and TEM revealed indigenous bacteria could obtain these vectors from E. coli through several mechanisms of horizontal gene transfer (HGT), including conjugation and cytoplasmic exchange through nanotubes. Inoculating petroleum-polluted sediments with E. coli carrying the vector pSF-OXB15-p450camfusion showed that the E. coli cells died after five days but a variety of bacteria received and carried the vector for over 60 days after inoculation. Within 60 days, the total petroleum hydrocarbon content of the polluted soil was reduced by 46%. Pilot experiments show that vectors only persist in indigenous populations when under selection pressure, disappearing when this carbon source is removed. This approach to remediation could prime indigenous bacteria for degrading pollutants while providing minimal ecosystem disturbance.

Proton relay network in P450cam formed upon docking of putidaredoxin

Cytochromes P450 are versatile heme-based enzymes responsible for vital life processes. Of these, P450cam (substrate camphor) has been most studied. Despite this, precise mechanisms of the key O─O cleavage step remain partly elusive to date; effects observed in various enzyme mutants remain partly unexplained. We have carried out extended (to 1000 ns) MM-MD and follow-on quantum mechanics/molecular mechanics computations, both on the well-studied FeOO state and on Cpd(0) (compound 0). Our simulations include (all camphor-bound): (a) WT (wild type), FeOO state. (b) WT, Cpd(0). (c) Pdx (Putidaredoxin, redox partner of P450)-docked-WT, FeOO state. (d) Pdx-docked WT, Cpd(0). (e) Pdx-docked T252A mutant, Cpd(0). Among our key findings: (a) Effect of Pdx docking appears to go far beyond that indicated in prior studies: it leads to specific alterations in secondary structure that create the crucial proton relay network. (b) Specific proton relay networks we identify are: FeOO(H)⋯T252⋯nH ₂ O⋯D251 in WT; FeOO(H)⋯nH ₂ O⋯D251 in T252A mutant; both occur with Pdx docking. (c) Direct interaction of D251 with -FeOOH is, respectively, rare/frequent in WT/T252A mutant. (d) In WT, T252 is in the proton relay network. (e) Positioning of camphor appears significant: when camphor is part of H-bonding network, second protonation appears to be facilitated.

Mapping the Substrate Recognition Pathway in Cytochrome P450

Cytochrome P450s are ubiquitous metalloenzymes involved in the metabolism and detoxification of foreign components via catalysis of the hydroxylation reactions of a vast array of organic substrates. However, the mechanism underlying the pharmaceutically critical process of substrate access to the catalytic center of cytochrome P450 is a long-standing puzzle, further complicated by the crystallographic evidence of a closed catalytic center in both substrate-free and substrate-bound cytochrome P450. Here, we address a crucial question whether the conformational heterogeneity prevalent in cytochrome P450 translates to heterogeneous pathways for substrate access to the catalytic center of these metalloenzymes. By atomistically capturing the full process of spontaneous substrate association from bulk solvent to the occluded catalytic center of an archetypal system P450cam in multi-microsecond-long continuous unbiased molecular dynamics simulations, we here demonstrate that the substrate recognition in P450cam always occurs through a single well-defined dominant pathway. The simulated final bound pose resulting from these unguided simulations is in striking resemblance with the crystallographic bound pose. Each individual binding trajectory reveals that the substrate, initially placed at random locations in bulk solvent, spontaneously lands on a single key channel on the protein-surface of P450cam and resides there for an uncharacteristically long period, before correctly identifying the occluded target-binding cavity. Surprisingly, the passage of substrate to the closed catalytic center is not accompanied by any large-scale opening in protein. Rather, the unbiased simulated trajectories (∼57 μs) and underlying Markov state model, in combination with free-energy analysis, unequivocally show that the substrate recognition process in P450cam needs a substrate-induced side-chain displacement coupled with a complex array of dynamical interconversions of multiple metastable substrate conformations. Further, the work reconciles multiple precedent experimental and theoretical observations on P450cam and establishes a comprehensive view of substrate-recognition in cytochrome P450 that only occurs via substrate-induced structural rearrangements.

Long-Range Modulations of Electric Fields in Proteins

Electrostatic interactions are essential for controlling the protein structure and function. Whereas so far experimental and theoretical efforts focused on the effect of local electrostatics, this work aims at elucidating the long-range modulation of electric fields in proteins upon binding to charged surfaces. The study is based on cytochrome c (Cytc) variants carrying nitrile reporters for the vibrational Stark effect that are incorporated into the protein via genetic engineering and chemical modification. The Cytc variants were thoroughly characterized with respect to possible structural perturbations due to labeling. For the proteins in solution, the relative hydrogen bond occupancy and the calculated electric fields, both obtained from molecular dynamics (MD) simulations, and the experimental nitrile stretching frequencies were used to develop a relationship for separating hydrogen-bonding and non-hydrogen-bonding electric field effects. This relationship provides an excellent description for the stable Cytc variants in solution. For the proteins bound to Au electrodes coated with charged self-assembled monolayers (SAMs), the underlying MD simulations can only account for the electric field changes Δ E_ads due to the formation of the electrostatic SAM-Cytc complexes but not for the additional contribution, Δ E_int, representing the consequences of the potential drops over the electrode/SAM/protein interfaces. Both Δ E_ads and Δ E_int, determined at distances between 20 and 30 Å with respect to the SAM surface, are comparable in magnitude to the non-hydrogen-bonding electric field in the unbound protein. This long-range modulation of the internal electric field may be of functional relevance for proteins in complexes with partner proteins (Δ E_ads) and attached to membranes (Δ E_ads + Δ E_int).