Conformational landscape and the selectivity of cytochrome P450cam

Emerging approaches to investigating functional protein dynamics in modular redox enzymes: Nitric oxide synthase as a model system

Approximately 80% of eukaryotic and 65% of prokaryotic proteins are composed of multiple folding units (i.e., domains) connected by flexible linkers. These dynamic protein architectures enable diverse, essential functions such as electron transfer, respiration, and biosynthesis. This review critically assesses recent advancements in methods for studying protein dynamics, with a particular focus on modular, multidomain nitric oxide synthase (NOS) enzymes. Moving beyond traditional static "snapshots" of protein structures, current research emphasizes the dynamic nature of proteins, viewing them as flexible architectures modulated by conformational changes and interactions. In this context, the review discusses key developments in the integration of quantitative crosslinking mass spectrometry (qXL MS) with AlphaFold 2 predictions, which provides a powerful approach to disentangling NOS structural dynamics and understanding their modulation by external regulatory cues. Additionally, advances in site-specific infrared (IR) spectroscopy offer exciting potential in providing rich details about the conformational dynamics of NOSs in docked states. Moreover, optimization of genetic code expansion machinery enables the generation of genuine phosphorylated NOS enzymes, paving the way for detailed biophysical and functional analyses of phosphorylation's role in shaping NOS activity and structural flexibility; notably, this approach also empowers site-specific IR probe labeling with cyano groups. By embracing and leveraging AI-driven tools like AlphaFold 2 for structural and conformational modeling, alongside solution-based biophysical methods such as qXL MS and site-specific IR spectroscopy, researchers will gain integrative insights into functional protein dynamics. Collectively, these breakthroughs highlight the transformative potential of modern approaches in driving fundamental biological chemistry research.

Facile Generation of Cyanoselenocysteine as a Vibrational Label for Measuring Protein Dynamics on Longer Time Scales by 2D IR Spectroscopy

Two-dimensional infrared (2D IR) spectroscopy is a powerful technique for measuring molecular heterogeneity and dynamics with a high spatiotemporal resolution. The methods can be applied to characterize specific residues of proteins by incorporating frequency-resolved vibrational labels. However, the time scale of dynamics that 2D IR spectroscopy can measure is limited by the vibrational label's excited-state lifetime due to the decay of 2D IR absorption bands. To extend this time scale, vibrational labels with longer lifetimes are sought. An effective approach to inhibiting intramolecular energy relaxation is to isolate the vibration from the rest of the molecule by inserting a heavy atom bridge. Although this strategy has been demonstrated through the generation of functionalized amino acids, a straightforward route to their selective incorporation into proteins is often unclear. A facile approach for the attachment of a cyano group at cysteine to generate a thiocyanate has contributed to its adoption as a vibrational label of proteins. We demonstrate that an analogous route can be used for introducing cyanoselenocysteine to generate a selenocyanate vibrational label containing a heavier bridge atom. We confirm by infrared pump-probe and 2D IR spectroscopy longer vibrational lifetimes of 100-250 ps, depending on the solvent, which enable the collection of 2D IR spectra to measure frequency dynamics on longer time scales.

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.

Interdomain Interactions Modulate the Active Site Dynamics of Human Inducible Nitric Oxide Synthase

Nitric oxide synthase (NOS) is a homodimeric flavohemoprotein responsible for catalyzing the oxidation of l-arginine (l-Arg) to citrulline and nitric oxide. Electrons are supplied for the reaction via interdomain electron transfer between an N-terminal heme-containing oxygenase domain and a FMN-containing (sub)domain of a C-terminal reductase domain. Extensive attention has focused on elucidating how conformational dynamics regulate electron transfer between the domains. Here we investigate the impact of the interdomain FMN-heme interaction on the heme active site dynamics of inducible NOS (iNOS). Steady state linear and time-resolved two-dimensional infrared (2D IR) spectroscopy was applied to probe a CO ligand at the heme within the oxygenase domain for full-length and truncated or mutated constructs of human iNOS. Whereas the linear IR spectra of the CO ligand were identical among the constructs, 2D IR spectroscopy revealed variation in the frequency dynamics. The wild-type constructs that can properly form the FMN/oxygenase docked state due to the presence of both the FMN and oxygenase domains showed slower dynamics than the oxygenase domain alone. Introduction of the mutation (E546N) predicted to perturb electrostatic interactions between the domains resulted in measured dynamics intermediate between those for the full-length and individual oxygenase domain, consistent with perturbation to the docked/undocked equilibrium. These results indicate that docking of the FMN domain to the oxygenase domain not only brings the FMN cofactor within electron transfer distance of the heme domain but also modulates the dynamics sensed by the CO ligand within the active site in a way expected to promote efficient electron transfer.

Transparent window 2D IR spectroscopy of proteins

Proteins are complex, heterogeneous macromolecules that exist as ensembles of interconverting states on a complex energy landscape. A complete, molecular-level understanding of their function requires experimental tools to characterize them with high spatial and temporal precision. Infrared (IR) spectroscopy has an inherently fast time scale that can capture all states and their dynamics with, in principle, bond-specific spatial resolution. Two-dimensional (2D) IR methods that provide richer information are becoming more routine but remain challenging to apply to proteins. Spectral congestion typically prevents selective investigation of native vibrations; however, the problem can be overcome by site-specific introduction of amino acid side chains that have vibrational groups with frequencies in the "transparent window" of protein spectra. This Perspective provides an overview of the history and recent progress in the development of transparent window 2D IR of proteins.

Protein Dynamics by Two-Dimensional Infrared Spectroscopy

Proteins function as ensembles of interconverting structures. The motions span from picosecond bond rotations to millisecond and longer subunit displacements. Characterization of functional dynamics on all spatial and temporal scales remains challenging experimentally. Two-dimensional infrared spectroscopy (2D IR) is maturing as a powerful approach for investigating proteins and their dynamics. We outline the advantages of IR spectroscopy, describe 2D IR and the information it provides, and introduce vibrational groups for protein analysis. We highlight example studies that illustrate the power and versatility of 2D IR for characterizing protein dynamics and conclude with a brief discussion of the outlook for biomolecular 2D IR.

Dynamics underlying hydroxylation selectivity of cytochrome P450cam

Structural heterogeneity and the dynamics of the complexes of enzymes with substrates can determine the selectivity of catalysis; however, fully characterizing how remains challenging as heterogeneity and dynamics can vary at the spatial level of an amino acid residue and involve rapid timescales. We demonstrate the nascent approach of site-specific two-dimensional infrared (IR) spectroscopy to investigate the archetypical cytochrome P450, P450cam, to better delineate the mechanism of the lower regioselectivity of hydroxylation of the substrate norcamphor in comparison to the native substrate camphor. Specific locations are targeted throughout the enzyme by selectively introducing cyano groups that have frequencies in a spectrally isolated region of the protein IR spectrum as local vibrational probes. Linear and two-dimensional IR spectroscopy were applied to measure the heterogeneity and dynamics at each probe and investigate how they differentiate camphor and norcamphor recognition. The IR data indicate that the norcamphor complex does not fully induce a large-scale conformational change to a closed state of the enzyme adopted in the camphor complex. Additionally, a probe directed at the bound substrate experiences rapidly interconverting states in the norcamphor complex that explain the hydroxylation product distribution. Altogether, the study reveals large- and small-scale structural heterogeneity and dynamics that could contribute to selectivity of a cytochrome P450 and illustrates the approach of site-selective IR spectroscopy to elucidate protein dynamics.

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.

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.

Site-Specific 1D and 2D IR Spectroscopy to Characterize the Conformations and Dynamics of Protein Molecular Recognition

Proteins exist as ensembles of interconverting states on a complex energy landscape. A complete, molecular-level understanding of their function requires knowledge of the populated states and thus the experimental tools to characterize them. Infrared (IR) spectroscopy has an inherently fast time scale that can capture all states and their dynamics with, in principle, bond-specific spatial resolution, and 2D IR methods that provide richer information are becoming more routine. Although application of IR spectroscopy for investigation of proteins is challenged by spectral congestion, the issue can be overcome by site-specific introduction of amino acid side chains that have IR probe groups with frequency-resolved absorptions, which furthermore enables selective characterization of different locations in proteins. Here, we briefly introduce the biophysical methods and summarize the current progress toward the study of proteins. We then describe our efforts to apply site-specific 1D and 2D IR spectroscopy toward elucidation of protein conformations and dynamics to investigate their involvement in protein molecular recognition, in particular mediated by dynamic complexes: plastocyanin and its binding partner cytochrome f, cytochrome P450s and substrates or redox partners, and Src homology 3 domains and proline-rich peptide motifs. We highlight the advantages of frequency-resolved probes to characterize specific, local sites in proteins and uncover variation among different locations, as well as the advantage of the fast time scale of IR spectroscopy to detect rapidly interconverting states. In addition, we illustrate the greater insight provided by 2D methods and discuss potential routes for further advancement of the field of biomolecular 2D IR spectroscopy.

Heterogeneous and Highly Dynamic Interface in Plastocyanin-Cytochrome f Complex Revealed by Site-Specific 2D-IR Spectroscopy

Transient protein complexes are crucial for sustaining dynamic cellular processes. The complexes of electron-transfer proteins are a notable example, such as those formed by plastocyanin (Pc) and cytochrome f (cyt f) in the photosynthetic apparatus. The dynamic and heterogeneous nature of these complexes, however, makes their study challenging. To better elucidate the complex of Nostoc Pc and cyt f, 2D-IR spectroscopy coupled to site-specific labeling with cyanophenylalanine infrared (IR) probes was employed to characterize how the local environments at sites along the surface of Pc were impacted by cyt f binding. The results indicate that Pc most substantially engages with cyt f via the hydrophobic patch around the copper redox site. Complexation with cyt f led to an increase in inhomogeneous broadening of the probe absorptions, reflective of increased heterogeneity of interactions with their environment. Notably, most of the underlying states interconverted very rapidly (1 to 2 ps), suggesting a complex with a highly mobile interface. The data support a model of the complex consisting of a large population of an encounter complex. Additionally, the study demonstrates the application of 2D-IR spectroscopy with site-specifically introduced probes to reveal new quantitative insight about dynamic biochemical systems.

Conformational Change Induced by Putidaredoxin Binding to Ferrous CO-ligated Cytochrome P450cam Characterized by 2D IR Spectroscopy

The importance of conformational dynamics to protein function is now well-appreciated. An outstanding question is whether they are involved in the effector role played by putidaredoxin (Pdx) in its reduction of the O₂ complex of cytochrome P450cam (P450cam), an archetypical member of the cytochrome P450 superfamily. Recent studies have reported that binding of Pdx induces a conformational change from a closed to an open state of ferric P450cam, but a similar conformational change does not appear to occur for the ferrous, CO-ligated enzyme. To better understand the effector role of Pdx when binding the ferrous, CO-ligated P450cam, we applied 2D IR spectroscopy to compare the conformations and dynamics of the wild-type (wt) enzyme in the absence and presence of Pdx, as well as of L358P P450cam (L358P), which has served as a putative model for the Pdx complex. The CO vibrations of the Pdx complex and L358P report population of two conformational states in which the CO experiences distinct environments. The dynamics among the CO frequencies indicate that the energy landscape of substates within one conformation are reflective of the closed state of P450cam, and for the other conformation, differ from the free wt enzyme, but are equivalent between the Pdx complex and L358P. The two states co-populated by the Pdx complex are postulated to reflect a loosely bound encounter complex and a more tightly bound state, as is commonly observed for the dynamic complexes of redox partners. Significantly, this study shows that the binding of Pdx to ferrous, CO-ligated P450cam does perturb the conformational ensemble in a way that might underlie the effector role of Pdx.

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

The broad and variable substrate specificity of cytochrome P450 enzymes makes them a model system for studying the determinants of protein molecular recognition. The archetypal cytochrome P450cam (P450cam) is a relatively specific P450, a feature once attributed to the high rigidity of its active site. However, increasingly studies have provided evidence of the importance of conformational changes to P450cam activity. Here we used infrared (IR) spectroscopy to investigate the molecular recognition of P450cam. Toward this goal, and to assess the influence of a hydrogen bond (H-bond) between active site residue Y96 and substrates, two variants in which Y96 is replaced by a cyanophenyl (Y96CNF) or phenyl (Y96F) group were characterized in complexes with the substrates camphor, isoborneol, and camphane. These combinations allow for a comparison of complexes in which the moieties on both the protein and substrate can serve as a H-bond donor, acceptor, or neither. The IR spectra of heme-bound CO and the site-specifically incorporated CN of Y96CNF were analyzed to characterize the number and nature of environments in each protein, both in the free and bound states. Although the IR spectra do not support the idea that protein-substrate H-bonding is central to P450cam recognition, the data altogether suggest that the differing conformational heterogeneity in the active site of the P450cam variants and changes in heterogeneity upon binding of different substrates likely contribute to their variable affinities via a conformational selection mechanism. This study further extends our understanding of the molecular recognition of archetypal P450cam and demonstrates the application of IR spectroscopy combined with selective protein modification to delineate protein-ligand interactions.

Extended timescale 2D IR probes of proteins: p-cyanoselenophenylalanine

The importance of dynamics to the function of proteins is well appreciated, but the difficulty in their measurement impedes investigation into their precise role(s). 2D IR spectroscopy is a developing approach for the study of dynamics and has motivated efforts to develop spectrally resolved IR probe groups that enable its application for measuring the dynamics at specific sites in a protein. A challenge with this approach is that the timescales accessible are limited by the vibrational lifetimes of the probes. Toward development of better probes for 2D IR spectroscopy of protein dynamics, we report the characterization of p-cyano-seleno-phenylalanine (CNSePhe), a derivative of the well established IR probe p-cyano-phenylalanine (CNPhe), by FT IR, pump-probe, and 2D IR spectroscopy. The incorporation of the heavy Se atom decouples the CN vibration from the rest in the molecule. Although this leads to a reduction of the transition dipole strength, and thus a reduction in signal intensity, it also dramatically increases the vibrational lifetime, enabling collection of 2D IR spectra for analysis of molecular dynamics on much longer timescales. Interestingly, we also find that the lifetime for CNSePhe shows increased sensitivity to the presence of hydrogen bonding interactions with the CN, suggesting that the probe should be useful for interpretation of CN spectra and possibly for the study of solvation.

Site-Specific Characterization of Cytochrome P450cam Conformations by Infrared Spectroscopy

Conformational changes are central to protein function but challenging to characterize with both high spatial and temporal precision. The inherently fast time scale and small chromophores of infrared (IR) spectroscopy are well-suited for characterization of potentially rapidly fluctuating environments, and when frequency-resolved probes are incorporated to overcome spectral congestion, enable characterization of specific sites in proteins. We selectively incorporated p-cyanophenylalanine (CNF) as a vibrational probe at five distinct locations in the enzyme cytochrome P450cam and used IR spectroscopy to characterize the environments in substrate and/or ligand complexes reflecting those in the catalytic cycle. Molecular dynamics (MD) simulations were performed to provide a structural basis for spectral interpretation. Together the experimental and simulation data suggest that the CN frequencies are sensitive to both long-range influences, resulting from the particular location of a residue within the enzyme, as well as short-range influences from hydrogen bonding and packing interactions. The IR spectra demonstrate that the environments and effects of substrate and/or ligand binding are different at each position probed and also provide evidence that a single site can experience multiple environments. This study illustrates how IR spectroscopy, when combined with the spectral decongestion and spatial selectivity afforded by CNF incorporation, provides detailed information about protein structural changes that underlie function.

Resolution of Site-Specific Conformational Heterogeneity in Proline-Rich Molecular Recognition by Src Homology 3 Domains

Conformational heterogeneity and dynamics are increasingly evoked in models of protein molecular recognition but are challenging to experimentally characterize. Here we combine the inherent temporal resolution of infrared (IR) spectroscopy with the spatial resolution afforded by selective incorporation of carbon-deuterium (C-D) bonds, which provide frequency-resolved absorptions within a protein IR spectrum, to characterize the molecular recognition of the Src homology 3 (SH3) domain of the yeast protein Sho1 with its cognate proline-rich (PR) sequence of Pbs2. The IR absorptions of C-D bonds introduced at residues along a peptide of the Pbs2 PR sequence report on the changes in the local environments upon binding to the SH3 domain. Interestingly, upon forming the complex the IR spectra of the peptides labeled with C-D bonds at either of the two conserved prolines of the PXXP consensus recognition sequence show more absorptions than there are C-D bonds, providing evidence for the population of multiple states. In contrast, the NMR spectra of the peptides labeled with (13)C at the same residues show only single resonances, indicating rapid interconversion on the NMR time scale. Thus, the data suggest that the SH3 domain recognizes its cognate peptide with a component of induced fit molecular recognition involving the adoption of multiples states, which have previously gone undetected due to interconversion between the populated states that is too fast to resolve using conventional methods.

Whole-Cell Biocatalysts for Stereoselective C-H Amination Reactions

Enantiomerically pure chiral amines are ubiquitous chemical building blocks in bioactive pharmaceutical products and their synthesis from simple starting materials is of great interest. One of the most attractive strategies is the stereoselective installation of a chiral amine through C-H amination, which is a challenging chemical transformation. Herein we report the application of a multienzyme cascade, generated in a single bacterial whole-cell system, which is able to catalyze stereoselective benzylic aminations with ee values of 97.5%. The cascade uses four heterologously expressed recombinant enzymes with cofactors provided by the host cell and isopropyl amine added as the amine donor. The cascade presents the first example of the successful de novo design of a single whole-cell biocatalyst for formal stereoselective C-H amination.

Considerations of Protein Subpockets in Fragment-Based Drug Design

While the fragment-based drug design approach continues to gain importance, gaps in the tools and methods available in the identification and accurate utilization of protein subpockets have limited the scope. The importance of these features of small molecule-protein recognition is highlighted with several examples. A generalized solution for the identification of subpockets and corresponding chemical fragments remains elusive, but there are numerous advancements in methods that can be used in combination to address subpockets. Finally, additional examples of approaches that consider the relative importance of small-molecule co-dependence of protein conformations are highlighted to emphasize an increased significance of subpockets, especially at protein interfaces.