Conformational plasticity and structure/function relationships in cytochromes P450

Tracking protein-protein interactions by NMR: conformational selection in human steroidogenic cytochrome P450 CYP17A1 induced by cytochrome b5

The human steroidogenic cytochrome P450 CYP17A1 catalyzes two types of reactions in the biosynthetic pathway leading from pregnenolone to testosterone and several other steroid hormones. The first is the hydroxylation of pregnenolone or progesterone to the corresponding 17α-hydroxy steroid, followed by a lyase reaction that converts these 17α-hydroxy intermediates to the androgens dehydroepiandrosterone and androstenedione, respectively. cytochrome b5 (cytb5) is known to act as both an effector and electron donor for the lyase oxidations, markedly stimulating the rate of the lyase reaction in its presence relative to the rate in its absence. Extensive sequential backbone 1H,15N and 13C nuclear magnetic resonance assignments have now been made for oxidized CYP17A1 bound to the prostate cancer drug and inhibitor abiraterone. This is the first eukaryotic P450 for which such assignments are now available. These assignments allow more complete interpretation of the structural perturbations observed upon cytb5 addition. Possible mechanism(s) for the effector activity of cytb5 are discussed in light of this new information.

Small angle X-ray scattering analysis of thermophilic cytochrome P450 CYP119 and the effects of the N-terminal histidine tag

Combining size exclusion chromatography-small angle X-ray scattering (SEC-SAXS) and molecular dynamics (MD) analysis is a promising approach to investigate protein behavior in solution, particularly for understanding conformational changes due to substrate binding in cytochrome P450s (CYPs). This study investigates conformational changes in CYP119, a thermophilic CYP from Sulfolobus acidocaldarius that exhibits structural flexibility similar to mammalian CYPs. Although the crystal structure of ligand-free (open state) and ligand-bound (closed state) forms of CYP119 is known, the overall structure of the enzyme in solution has not been explored until now. It was found that theoretical scattering profiles from the crystal structures of CYP119 did not align with the SAXS data, but conformers from MD simulations, particularly starting from the open state (46 % of all frames), agreed well. Interestingly, a small percentage of closed-state conformers also fit the data (9 %), suggesting ligand-free CYP119 samples ligand-bound conformations. Ab initio SAXS models for N-His tagged CYP119 revealed a tail-like unfolded structure impacting protein flexibility, which was confirmed by in silico modeling. SEC-SAXS analysis of N-His CYP119 indicated pentameric structures in addition to monomers in solution, affecting the stability and activity of the enzyme. This study adds insights into the conformational dynamics of CYP119 in solution.

A Mechanistic Understanding of the Distinct Regio- and Chemoselectivity of Multifunctional P450s by Structural Comparison of IkaD and CftA Complexed with Common Substrates

Regio- and chemoselective C-H activation at multi-positions of a single molecule is fascinating but chemically challenging. The homologous cytochrome P450 enzymes IkaD and CftA catalyze multiple C-H oxidations on the same polycyclic tetramate macrolactam (PoTeM) ikarugamycin, with distinct regio- and chemoselectivity. Herein we provide mechanistic understanding of their functional differences by solving crystal structures of IkaD and CftA in complex with ikarugamycin and unnatural substrates. Distinct conformations of the F/G region in IkaD and CftA are found to differentiate the orientation of PoTeM substrates, by causing different binding patterns with polar moieties to determine site selection, oxidation order, and chemoselectivity. Fine-tuning the polar subpocket altered the regioselectivity of IkaD, indicating that substrate re-orientation by mutating residues distal to the oxidation site could serve as an important method in future engineering of P450 enzymes.

Structural basis of bidirectional allostery across the heme in a cytochrome P450 enzyme

Cytochromes P450 (CYPs) are heme-containing enzymes that are present in all kingdoms of life and share a structurally homologous, globular protein fold. CYPs utilize structures distal to the heme to recognize and coordinate substrates, while the necessary interactions with redox partner proteins are mediated at the opposite, proximal surface. In the current study, we investigated the functional allostery across the heme for the bacterial enzyme CYP121A1, which utilizes a non-polar distal-to-distal dimer interface for specific binding of its dicyclotyrosine substrate. Fluorine-detected Nuclear Magnetic Resonance (19F-NMR) spectroscopy was combined with site-specific labeling of a distal surface residue (S171C of the FG-loop), one residue of the B-helix (N84C), and two proximal surface residues (T103C and T333C) with a thiol-reactive fluorine label. Adrenodoxin was used as a substitute redox protein and was found to promote a closed arrangement of the FG-loop, similar to the addition of substrate alone. Disruption of the protein-protein interface by mutagenesis of two CYP121 basic surface residues removed the allosteric effect. Moreover, 19F-NMR spectra of the proximal surface indicate that ligand-induced allostery modulates the environment at the C-helix but not the meander region of the enzyme. In light of the high degree of structural homology in this family of enzymes, we interpret the findings from this work to represent a conserved allosteric network in CYPs.

Prediction of Sites of Metabolism of CYP3A4 Substrates Utilizing Docking-Derived Geometric Features

Cytochrome P450 3A4 (CYP3A4) is one of the major drug-metabolizing enzymes in the human body and is responsible for the metabolism of ∼50% of clinically used drugs. Therefore, the identification of the compound's sites of metabolism (SOMs) mediated by CYP3A4 is of utmost importance in the early stage of drug discovery and development. Herein, docking-based approaches incorporating geometric features were used for SOMs prediction of CYP3A4 substrates. The cross-docking poses of a relatively large data set containing 474 substrates were analyzed in depth, and a widely observed geometric pattern called the close proximity of SOMs was derived from the poses. On the basis of the close proximity, several structure-based models have been constructed, which demonstrated better performance than those structure-based models using the criterion of Fe-SOM distance. For further improving the prediction performance, the structure-based models were also combined with the well-known ligand-based model SMARTCyp. One combined model exhibited good performance on the SOMs prediction of an external substrate set containing kinase inhibitors, PROTACs, approved drugs, and some lead compounds.

Key regulators in the architecture of substrate access/egress channels in mammalian cytochromes P450 governing flexibility in substrate oxyfunctionalization

Cytochrome P450s (CYP) represent a superfamily of b-type hemoproteins catalyzing oxifunctionalization of a vast array of endogenous and exogenous compounds. The present review focuses on assessment of the topology of prospective determinants in substrate entry and product release channels of mammalian P450s, steering the conformational dynamics of substrate accessibility and productive ligand orientation toward the iron-oxene core. Based on a generalized, CYP3A4-related construct, the sum of critical elements from diverse target enzymes was found to cluster within the known substrate recognition sites. The majority of prevalent substrate access/egress tunnels revealed to be of fairly balanced functional importance. The hydrophobicity profile of the candidates revealed to be the most salient feature in functional interaction throughout the conduits, while bulkiness of the residues imposes steric restrictions on substrate traveling. Thus, small amino acids such as prolines and glycines serve as hinges, driving conformational flexibility in ligand passage. Similarly, bottlenecks in the tunnel architecture, being narrowest encounter points within the CYP3A4 model, have a vital function in substrate selectivity along with clusters of aromatic amino acids acting as gatekeepers. In addition, peripheral patches in conduits may house determinants modulating allosteric cooperativity between remote and central domains in the P450 structure. Remarkably, the bulk critical residues lining tunnels in the various isozymes reside in helices B'/C and F/G inclusive of their interhelical turns as well as in helix I. This suggests these regions to represent hotspots for targeted genetic engineering to tailor more sophisticated mammalian P450s exploitable in industrial, biotechnological and medicinal areas.

Structural comparison of the cytochrome P450 enzymes CYP106A1 and CYP106A2 provides insight into their differences in steroid conversion

Understanding the structural basis of the selectivity of steroid hydroxylation requires detailed structural and functional investigations on various steroid hydroxylases with different selectivities, such as the bacterial cytochrome P450 enzymes. Here, the crystal structure of the cytochrome P450 CYP106A1 from Priestia megaterium was solved. CYP106A1 exhibits a rare additional structural motif of a cytochrome P450, a sixth β-sheet. The protein was found in different unusual conformations corresponding to both open and closed forms even when crystallized without any known substrate. The structural comparison of CYP106A1 with the previously investigated CYP106A2, including docking studies for both isoforms with the substrate cortisol, reveals a completely different orientation of the steroid molecule in the active sites. This distinction convincingly explains the experimentally observed differences in substrate conversion and product formation by the two enzymes.

Mechanism of Melatonin Metabolism by CYP1A1: What Determines the Bifurcation Pathways of Hydroxylation versus Deformylation?

Melatonin, a widely applied cosmetic active ingredient, has a variety of uses as a skin protector through antioxidant and anti-inflammatory functions as well as giving the body UV-induced defenses and immune system support. In the body, melatonin is synthesized from a tryptophan amino acid in a cascade of reactions, but as melatonin is toxic at high concentrations, it is metabolized in the human skin by the cytochrome P450 enzymes. The P450s are diverse heme-based mono-oxygenases that catalyze oxygen atom-transfer processes that trigger metabolism and detoxification reactions in the body. In the catalytic cycle of the P450s, a short-lived high-valent iron(IV)-oxo heme cation radical is formed that has been proposed to be the active oxidant. How and why it activates melatonin in the human body and what the origin of the product distributions is, are unknown. This encouraged us to do a detailed computational study on a typical human P450 isozyme, namely CYP1A1. We initially did a series of molecular dynamics simulations with substrate docked into several orientations. These simulations reveal a number of stable substrate-bound positions in the active site, which may lead to differences in substrate activation channels. Using tunneling analysis on the full protein structures, we show that two of the four binding conformations lead to open substrate-binding pockets. As a result, in these open pockets, the substrate is not tightly bound and can escape back into the solution. In the closed conformations, in contrast, the substrate is mainly oriented with the methoxy group pointing toward the heme, although under a different angle. We then created large quantum cluster models of the enzyme and focused on the chemical reaction mechanisms for melatonin activation, leading to competitive O-demethylation and C⁶-aromatic hydroxylation pathways. The calculations show that active site positioning determines the product distributions, but the bond that is activated is not necessarily closest to the heme in the enzyme-substrate complex. As such, the docking and molecular dynamics positioning of the substrate versus oxidant can give misleading predictions on product distributions. In particular, in quantum mechanics cluster model I, we observe that through a tight hydrogen bonding network, a preferential 6-hydroxylation of melatonin is obtained. However, O-demethylation becomes possible in alternative substrate-binding orientations that have the C⁶-aromatic ring position shielded. Finally, we investigated enzymatic and non-enzymatic O-demethylation processes and show that the hydrogen bonding network in the substrate-binding pocket can assist and perform this step prior to product release from the enzyme.

A unified GCNN model for predicting CYP450 inhibitors by using graph convolutional neural networks with attention mechanism

Undesirable drug-drug interactions (DDIs) may lead to serious adverse side effects when more than two drugs are administered to a patient simultaneously. One of the most common DDIs is caused by unexpected inhibition of a specific human cytochrome P450 (CYP450), which plays a dominant role in the metabolism of the co-administered drugs. Therefore, a unified and reliable method for predicting the potential inhibitors of CYP450 family is extremely important in drug development. In this work, graph convolutional neural network (GCN) with attention mechanism and 1-D convolutional neural network (CNN) were used to extract the features of CYP ligands and the binding sites of CYP450 respectively, which were then combined to establish a unified GCN-CNN (GCNN) model for predicting the inhibitors of 5 dominant CYP isoforms, i.e., 1A2, 2C9, 2C19, 2D6, and 3A4. Overall, the established GCNN model showed good performances on the test samples and achieved better performances than the recently proposed iCYP-MFE model by using the same datasets. Based on the heat-map analysis of the resulting molecular graphs, the key structural determinants of the CYP inhibitors were further explored.

Generation of new inhibitors of selected cytochrome P450 subtypes- In silico study

Physicochemical and pharmacokinetic compound profile has crucial impact on compound potency to become a future drug. Ligands with desired activity profile cannot be used for treatment if they are characterized by unfavourable physicochemical or ADMET properties. In the study, we consider metabolic stability and focus on selected subtypes of cytochrome P450 - proteins, which take part in the first phase of compound transformations in the organism. We develop a protocol for generation of new potential inhibitors of selected cytochrome isoforms. Its subsequent stages are composed of generation and assessment of new derivatives of known cytochrome inhibitors, docking and evaluation of the compound possible inhibition on the basis of the obtained ligand-protein complexes. Besides the library of new potential agents inhibiting particular cytochrome subtypes, we also prepare a graph neural network that predicts the change in activity for all modifications of the starting molecule. In addition, we perform a systematic statistical study on the influence of particular substitutions on the potential inhibition properties of generated compounds (both mono- and di-substitutions are considered), provide explanations of the inhibitory predictions and prepare an on-line visualization platform enabling manual inspection of the results. The developed methodology can greatly support the design of new cytochrome P450 inhibitors with the overarching goal of generation of new metabolically stable compounds. It enables instant evaluation of possible compound-cytochrome interactions and selection of ligands with the highest potential of possessing desired biological activity.

Carving the Active Site of CYP153A7 Monooxygenase for Improving Terminal Hydroxylation of Medium-Chain Fatty Acids

The P450-mediated terminal hydroxylation of non-activated C-H bonds is a chemically challenging reaction. CYP153A7 monooxygenase discovered in Sphingomonas sp. HXN200 belongs to the CYP153A subfamily and shows a pronounced terminal selectivity. Herein, we report the significantly improved terminal hydroxylation activity of CYP153A7 by redesign of the substrate binding pocket based on molecular docking of CYP153A7-C 8:0 and sequence alignments. Some of the resultant single mutants were advantageous over the wild-type enzyme with higher reaction rates, achieving a complete conversion of n- octanoic acid (C 8:0. 1 mM) in a shorter period. Especially, a single-mutation variant, D258E, showed 3.8-fold higher catalytic efficiency than the wild type toward the terminal hydroxylation of medium-chain fatty acid C 8:0 into the high value-added product 8-hydroxyoctanoic acid.

Multifunctional P450 Monooxygenase CftA Diversifies the Clifednamide Pool through Tandem C-H Bond Activations

Polycyclic tetramate macrolactams (PTMs) are a class of structurally complex hybrid polyketide-nonribosomal peptide (PK-NRP) natural products produced by diverse bacteria. Several PTMs display pharmaceutically interesting bioactivities, and the early stages of PTM biosynthesis involving polyketide synthase (PKS) and nonribosomal peptide synthetase (NRPS) enzymology are well studied. However, the timing and mechanisms of post PKS-NRPS oxidations by P450 monooxygenases encoded in PTM biosynthetic gene clusters (BGCs) remain poorly characterized. Here we demonstrate that CftA, encoded in clifednamide-type PTM BGCs, is a multifunctional P450 monooxygenase capable of converting the C29-C30 ethyl side chain of ikarugamycin to either a C29-C30 methyl ketone or a C29-C30 hydroxymethyl ketone through C-H bond activation, resulting in the formation of clifednamide A or clifednamide C, respectively. We also report the complete structure of clifednamide C solved via multidimensional NMR (COSY, HSQC, HMBC, NOESY, and TOCSY) using material purified from an engineered Streptomyces strain optimized for production. Finally, the in vitro reconstitution of recombinant CftA catalytic activity revealed the oxidation cascade for sequential conversion of ikarugamycin to clifednamide A and clifednamide C. Our findings confirm prior genetics-based predictions on the origins of clifednamide complexity via P450s encoded in PTM BGCs and place CftA into a growing group of multifunctional P450s that tailor PTM natural products through late-stage regioselective C-H bond activation.

Point Mutations at a Key Site Alter the Cytochrome P450 OleP Structural Dynamics

Substrate binding to the cytochrome P450 OleP is coupled to a large open-to-closed transition that remodels the active site, minimizing its exposure to the external solvent. When the aglycone substrate binds, a small empty cavity is formed between the I and G helices, the BC loop, and the substrate itself, where solvent molecules accumulate mediating substrate-enzyme interactions. Herein, we analyzed the role of this cavity in substrate binding to OleP by producing three mutants (E89Y, G92W, and S240Y) to decrease its volume. The crystal structures of the OleP mutants in the closed state bound to the aglycone 6DEB showed that G92W and S240Y occupied the cavity, providing additional contact points with the substrate. Conversely, mutation E89Y induces a flipped-out conformation of this amino acid side chain, that points towards the bulk, increasing the empty volume. Equilibrium titrations and molecular dynamic simulations indicate that the presence of a bulky residue within the cavity impacts the binding properties of the enzyme, perturbing the conformational space explored by the complexes. Our data highlight the relevance of this region in OleP substrate binding and suggest that it represents a key substrate-protein contact site to consider in the perspective of redirecting its activity towards alternative compounds.

Assembling the P450 puzzle: on the sources of nonadditivity in drug metabolism

There is an increasing number of indications of an oversimplification in the premise that the cumulative properties of the human drug-metabolizing ensemble represent a simple aggregate of the properties of the constituting enzymes. Recent studies of the functional effects of hetero-association of multiple cytochrome P450 species and their interactions with metabolically related enzymes revealed a tight integration in the drug-metabolizing ensemble. In our opinion, the sources of interindividual variability in drug metabolism can be elucidated only when considering this ensemble as a multienzyme system, the functional parameters of which are determined by interactions between its constituents. In this article, we present a conceptual model providing a mechanistic explanation for the functional effects of the interactions between multiple P450 species and propose a clue to understanding the nonadditive behavior of the drug-metabolizing ensemble.

How can SHAP values help to shape metabolic stability of chemical compounds?

Background

Computational methods support nowadays each stage of drug design campaigns. They assist not only in the process of identification of new active compounds towards particular biological target, but also help in the evaluation and optimization of their physicochemical and pharmacokinetic properties. Such features are not less important in terms of the possible turn of a compound into a future drug than its desired affinity profile towards considered proteins. In the study, we focus on metabolic stability, which determines the time that the compound can act in the organism and play its role as a drug. Due to great complexity of xenobiotic transformation pathways in the living organisms, evaluation and optimization of metabolic stability remains a big challenge.

Results

Here, we present a novel methodology for the evaluation and analysis of structural features influencing metabolic stability. To this end, we use a well-established explainability method called SHAP. We built several predictive models and analyse their predictions with the SHAP values to reveal how particular compound substructures influence the model’s prediction. The method can be widely applied by users thanks to the web service, which accompanies the article. It allows a detailed analysis of SHAP values obtained for compounds from the ChEMBL database, as well as their determination and analysis for any compound submitted by a user. Moreover, the service enables manual analysis of the possible structural modifications via the provision of analogous analysis for the most similar compound from the ChEMBL dataset.

Conclusions

To our knowledge, this is the first attempt to employ SHAP to reveal which substructural features are utilized by machine learning models when evaluating compound metabolic stability. The accompanying web service for metabolic stability evaluation can be of great help for medicinal chemists. Its significant usefulness is related not only to the possibility of assessing compound stability, but also to the provision of information about substructures influencing this parameter. It can assist in the design of new ligands with improved metabolic stability, helping in the detection of privileged and unfavourable chemical moieties during stability optimization. The tool is available at https://metstab-shap.matinf.uj.edu.pl/.

A rotamer relay information system in the epidermal growth factor receptor-drug complexes reveals clues to new paradigm in protein conformational change

Cancer cells can escape the effects of chemotherapy through mutations and upregulation of a tyrosine kinase protein called the epidermal growth factor receptor (EGFR). In the past two decades, four generations of tyrosine kinase inhibitors targeting EGFR have been developed. Using comparative structure analysis of 116 EGFR-drug complex crystal structures, cluster analysis produces two clans of 73 and 43 structures, respectively. The first clan of 73 structures is larger and is comprised mostly of the C-helix-IN conformation while the second clan of 43 structures correlates with the C-helix-OUT conformation. A deep rotamer analysis identifies 43 residues (18%) of the total of 237 residues spanning the kinase structures under investigation with significant rotamer variations between the C-helix-IN and C-helix-OUT clans. The locations of these rotamer variations take on the appearance of side chain conformational relays extending out from points of EGFR mutation to different regions of the EGFR kinase. Accordingly, we propose that key EGFR mutations act singly or together to induce drug resistant conformational changes in EGFR that are communicated via these side chain conformational relays. Accordingly, these side chain conformational relays appear to play a significant role in the development of tumour resistance. This phenomenon also suggests a new paradigm in protein conformational change that is mediated by supportive relays of rotamers on the protein surface, rather than through conventional backbone movements.

Functional Characterization of 21 Rare Allelic CYP1A2 Variants Identified in a Population of 4773 Japanese Individuals by Assessing Phenacetin O-Deethylation

Cytochrome P450 1A2 (CYP1A2), which accounts for approximately 13% of the total hepatic cytochrome content, catalyzes the metabolic reactions of approximately 9% of frequently used drugs, including theophylline and olanzapine. Substantial inter-individual differences in enzymatic activity have been observed among patients, which could be caused by genetic polymorphisms. Therefore, we functionally characterized 21 novel CYP1A2 variants identified in 4773 Japanese individuals by determining the kinetic parameters of phenacetin O-deethylation. Our results showed that most of the evaluated variants exhibited decreased or no enzymatic activity, which may be attributed to potential structural alterations. Notably, the Leu98Gln, Gly233Arg, Ser380del Gly454Asp, and Arg457Trp variants did not exhibit quantifiable enzymatic activity. Additionally, three-dimensional (3D) docking analyses were performed to further understand the underlying mechanisms behind variant pharmacokinetics. Our data further suggest that despite mutations occurring on the protein surface, accumulating interactions could result in the impairment of protein function through the destabilization of binding regions and changes in protein folding. Therefore, our findings provide additional information regarding rare CYP1A2 genetic variants and how their underlying effects could clarify discrepancies noted in previous phenotypical studies. This would allow the improvement of personalized therapeutics and highlight the importance of identifying and characterizing rare variants.

Polymorphism on human aromatase affects protein dynamics and substrate binding: spectroscopic evidence

Human aromatase is a member of the cytochrome P450 superfamily, involved in steroid hormones biosynthesis. In particular, it converts androgen into estrogens being therefore responsible for the correct sex steroids balance. Due to its capacity in producing estrogens it has also been considered as a promising target for breast cancer therapy. Two single-nucleotide polymorphisms (R264C and R264H) have been shown to alter aromatase activity and they have been associated to an increased or decreased risk for estrogen-dependent pathologies. Here, the effect of these mutations on the protein dynamics is investigated by UV/FTIR and time resolved fluorescence spectroscopy. H/D exchange rates were measured by FTIR for the three proteins in the ligand-free, substrate- and inhibitor-bound forms and the data indicate that the wild-type enzyme undergoes a conformational change leading to a more compact tertiary structure upon substrate or inhibitor binding. Indeed, the H/D exchange rates are decreased when a ligand is present. In the variants, the exchange rates in the ligand-free and -bound forms are similar, indicating that a structural change is lacking, despite the single amino acid substitution is located in the peripheral shell of the protein molecule. Moreover, the fluorescence lifetimes data show that the quenching effect on tryptophan-224 observed upon ligand binding in the wild-type, is absent in both variants. Since this residue is located in the catalytic pocket, these findings suggest that substrate entrance and/or retention in the active site is partially compromised in both mutants. A contact network analysis demonstrates that the protein structure is organized in two main clusters, whose connectivity is altered by ligand binding, especially in correspondence of helix-G, where the amino acid substitutions occur. Our findings demonstrate that SNPs resulting in mutations on aromatase surface modify the protein flexibility that is required for substrate binding and catalysis. The cluster analysis provides a rationale for such effect, suggesting helix G as a possible target for aromatase inhibition.

Conformational Plasticity-Rigidity Axis of the Coagulation Factor VII Zymogen Elucidated by Atomistic Simulations of the N-Terminally Truncated Factor VIIa Protease Domain

The vast majority of coagulation factor VII (FVII), a trypsin-like protease, circulates as the inactive zymogen. Activated FVII (FVIIa) is formed upon proteolytic activation of FVII, where it remains in a zymogen-like state and it is fully activated only when bound to tissue factor (TF). The catalytic domains of trypsin-like proteases adopt strikingly similar structures in their fully active forms. However, the dynamics and structures of the available corresponding zymogens reveal remarkable conformational plasticity of the protease domain prior to activation in many cases. Exactly how ligands and cofactors modulate the conformational dynamics and function of these proteases is not entirely understood. Here, we employ atomistic simulations of FVIIa (and variants hereof, including a TF-independent variant and N-terminally truncated variants) to provide fundamental insights with atomistic resolution into the plasticity-rigidity interplay of the protease domain conformations that appears to govern the functional response to proteolytic and allosteric activation. We argue that these findings are relevant to the FVII zymogen, whose structure has remained elusive despite substantial efforts. Our results shed light on the nature of FVII and demonstrate how conformational dynamics has played a crucial role in the evolutionary adaptation of regulatory mechanisms that were not present in the ancestral trypsin. Exploiting this knowledge could lead to engineering of protease variants for use as next-generation hemostatic therapeutics.

Surface hydrophobics mediate functional dimerization of CYP121A1 of Mycobacterium tuberculosis

Tuberculosis is caused by the pathogenic bacterium Mycobacterium tuberculosis (Mtb) and remains the leading cause of death by infection world-wide. The Mtb genome encodes a disproportionate number of twenty cytochrome P450 enzymes, of which the essential enzyme cytochrome P450 121A1 (CYP121A1) remains a target of drug design efforts. CYP121A1 mediates a phenol coupling reaction of the tyrosine dipeptide cyclo-L-Tyr-L-Tyr (cYY). In this work, a structure and function investigation of dimerization was performed as an overlooked feature of CYP121A1 function. This investigation showed that CYP121A1 dimers form via intermolecular contacts on the distal surface and are mediated by a network of solvent-exposed hydrophobic residues. Disruption of CYP121A1 dimers by site-directed mutagenesis leads to a partial loss of specificity for cYY, resulting in an approximate 75% decrease in catalysis. ¹⁹F labeling and nuclear magnetic resonance of the enzyme FG-loop was also combined with protein docking to develop a working model of a functional CYP121A1 dimer. The results obtained suggest that participation of a homodimer interface in substrate selectivity represents a novel paradigm of substrate binding in CYPs, while also providing important mechanistic insight regarding a relevant drug target in the development of novel anti-tuberculosis agents.

A dynamic understanding of cytochrome P450 structure and function through solution NMR

Many economically important biosyntheses incorporate regiospecific and stereospecific oxidations at unactivated carbons. Such oxidations are commonly catalyzed by cytochrome P450 monooxygenases, heme-containing enzymes that activate molecular oxygen while selectively binding and orienting the substrate for reaction. Despite the plethora of P450-catalyzed reactions, the P450 fold is highly conserved, and static structures are often insufficient for characterizing conformational states that contribute to specificity. High-resolution solution nuclear magnetic resonance (NMR) offers insights into dynamic processes and conformational changes that are required of a P450 in order to attain the combination of specificity and efficiency required for these reactions.

Dissecting the Structural Plasticity and Dynamics of Cytochrome P450 2B4 by Molecular Dynamics Simulations

The plasticity of cytochromes P450 (P450s) is known to contribute significantly to their catalytic capacity of metabolizing various substrates. Although numerous studies have been performed, factors governing the plasticity and dynamics of P450s are still not fully understood. In this study, taking CYP2B4 as an example, we dissect the protein plasticity and dynamics in different environments. CYP2B4 is featured by a high degree of plasticity, which exhibits open, closed, and intermediate states. By analyzing the CYP2B4 crystal structures, we identified the structural features for the closed, open, and intermediate states. Interestingly, formation of the dimer structure was found in the open and intermediate states. The subsequent molecular dynamics (MD) simulations of the open structure in water confirmed the importance of the dimer form in stabilizing the open conformations. MD simulations of the closed and open structures in the membrane environment and the free energies for opening the F-G cassette obtained from the umbrella sampling calculations indicate that the membrane environment is important for stabilizing the F-G cassette. The dynamical network analysis indicates that Asp105 on the B-C loop plays an important role in transiting the structure from the open to the intermediate state. Our results thus unveil the mechanisms of dimer formation and open-to-intermediate transition for CYP2B4 in the water and membrane environments.

Structure and Function of the Cytochrome P450 Monooxygenase Cinnamate 4-hydroxylase from Sorghum bicolor

Cinnamate 4-hydroxylase (C4H; CYP73A) is a cytochrome P450 monooxygenase associated externally with the endoplasmic reticulum of plant cells. The enzyme uses NADPH-cytochrome P450 reductase as a donor of electrons and hydroxylates cinnamic acid to form 4-coumaric acid in phenylpropanoid metabolism. In order to better understand the structure and function of this unique class of plant P450 enzymes, we have characterized the enzyme C4H1 from lignifying tissues of sorghum (Sorghum bicolor), encoded by Sobic.002G126600 Here we report the 1.7 Å resolution crystal structure of CYP73A33. The obtained structural information, along with the results of the steady-state kinetic analysis and the absorption spectroscopy titration, displays a high degree of similarity of the structural and functional features of C4H to those of other P450 proteins. Our data also suggest the presence of a putative allosteric substrate-binding site in a hydrophobic pocket on the enzyme surface. In addition, comparing the newly resolved structure with those of well-investigated cytochromes P450 from mammals and bacteria enabled us to identify those residues of critical functional importance and revealed a unique sequence signature that is potentially responsible for substrate specificity and catalytic selectivity of C4H.

Dimension-Tunable Circularly Polarized Luminescent Nanoassemblies with Emerging Selective Chirality and Energy Transfer

The selective interplay between dimensional morphology transition and signal transfer is an important feature for both nanomaterials and biosystems. While most of those reported examples considered either dimensional transition or signal transfer, the integrated interplay or selectivity for these two aspects in single self-assembled system has been rarely studied. Here, we report that a positively charged chiral π-building block could self-assemble into multidimensional nanostructures, which showed tunable circularly polarized luminescence (CPL). Impressively, when these CPL-active multidimensional structures interacted with two achiral dyes (positively charged ThT and negatively charged CNA), 3D nanocubes and 0D nanospheres showed neither chirality transfer nor energy transfer, while 2D nanoplates could successfully trigger a selective chirality or energy transfer depending on the charge type of acceptor dyes, which then emitted an enhanced CPL signal. This work demonstrated rational design of charged π-building block for the construction of dimension controllable and selective signal transfer self-assembly system, which might deepen the understanding the interplay of dimensional structures and signal transfer functions in natural and nano systems.

The Biosynthesis of Enzymatically Oxidized Lipids

Enzymatically oxidized lipids are a specific group of biomolecules that function as key signaling mediators and hormones, regulating various cellular and physiological processes from metabolism and cell death to inflammation and the immune response. They are broadly categorized as either polyunsaturated fatty acid (PUFA) containing (free acid oxygenated PUFA "oxylipins", endocannabinoids, oxidized phospholipids) or cholesterol derivatives (oxysterols, steroid hormones, and bile acids). Their biosynthesis is accomplished by families of enzymes that include lipoxygenases (LOX), cyclooxygenases (COX), cytochrome P450s (CYP), and aldo-keto reductases (AKR). In contrast, non-enzymatically oxidized lipids are produced by uncontrolled oxidation and are broadly considered to be harmful. Here, we provide an overview of the biochemistry and enzymology of LOXs, COXs, CYPs, and AKRs in humans. Next, we present biosynthetic pathways for oxylipins, oxidized phospholipids, oxysterols, bile acids and steroid hormones. Last, we address gaps in knowledge and suggest directions for future work.

Exploring the molecular basis for substrate specificity in homologous macrolide biosynthetic cytochromes P450

Cytochromes P450 (P450s) are nature's catalysts of choice for performing demanding and physiologically vital oxidation reactions. Biochemical characterization of these enzymes over the past decades has provided detailed mechanistic insight and highlighted the diversity of substrates P450s accommodate and the spectrum of oxidative transformations they catalyze. Previously, we discovered that the bacterial P450 MycCI from the mycinamicin biosynthetic pathway in Micromonospora griseorubida possesses an unusually broad substrate scope, whereas the homologous P450 from tylosin-producing Streptomyces fradiae (TylHI) exhibits a high degree of specificity for its native substrate. Here, using biochemical, structural, and computational approaches, we aimed to understand the molecular basis for the disparate reactivity profiles of these two P450s. Turnover and equilibrium binding experiments with substrate analogs revealed that TylHI strictly prefers 16-membered ring macrolides bearing the deoxyamino sugar mycaminose. To help rationalize these results, we solved the X-ray crystal structure of TylHI in complex with its native substrate at 1.99-Å resolution and assayed several site-directed mutants. We also conducted molecular dynamics simulations of TylHI and MycCI and biochemically characterized a third P450 homolog from the chalcomycin biosynthetic pathway in Streptomyces bikiniensis These studies provided a basis for constructing P450 chimeras to gain further insight into the features dictating the differences in reaction profile among these structurally and functionally related enzymes, ultimately unveiling the central roles of key loop regions in influencing substrate binding and turnover. Our work highlights the complex nature of P450/substrate interactions and raises interesting questions regarding the evolution of functional diversity among biosynthetic enzymes.

Directed Evolution of P450 BM3 towards Functionalization of Aromatic O-Heterocycles

The O-heterocycles, benzo-1,4-dioxane, phthalan, isochroman, 2,3-dihydrobenzofuran, benzofuran, and dibenzofuran are important building blocks with considerable medical application for the production of pharmaceuticals. Cytochrome P450 monooxygenase (P450) Bacillus megaterium 3 (BM3) wild type (WT) from Bacillus megaterium has low to no conversion of the six O-heterocycles. Screening of in-house libraries for active variants yielded P450 BM3 CM1 (R255P/P329H), which was subjected to directed evolution and site saturation mutagenesis of four positions. The latter led to the identification of position R255, which when introduced in the P450 BM3 WT, outperformed all other variants. The initial oxidation rate of nicotinamide adenine dinucleotide phosphate (NADPH) consumption increased ≈140-fold (WT: 8.3 ± 1.3 min-1; R255L: 1168 ± 163 min-1), total turnover number (TTN) increased ≈21-fold (WT: 40 ± 3; R255L: 860 ± 15), and coupling efficiency, ≈2.9-fold (WT: 8.8 ± 0.1%; R255L: 25.7 ± 1.0%). Computational analysis showed that substitution R255L (distant from the heme-cofactor) does not have the salt bridge formed with D217 in WT, which introduces flexibility into the I-helix and leads to a heme rearrangement allowing for efficient hydroxylation.

What Your Crystal Structure Will Not Tell You about Enzyme Function

Enzyme function requires that enzyme structures be dynamic. Substrate binding, product release, and transition state stabilization typically involve different enzyme conformers. Furthermore, in multistep enzyme-catalyzed reactions, more than one enzyme conformation may be important for stabilizing different transition states. While X-ray crystallography provides the most detailed structural information of any current methodology, X-ray crystal structures of enzymes capture only those conformations that fit into the crystal lattice, which may or may not be relevant to function. Solution nuclear magnetic resonance (NMR) methods can provide an alternative approach to characterizing enzymes under nonperturbing and controllable conditions, allowing one to identify and localize dynamic processes that are important to function. However, many enzymes are too large for standard approaches to making sequential resonance assignments, a critical first step in analyzing and interpreting the wealth of information inherent in NMR spectra.

This Account describes our long-standing NMR-based research into structural and dynamic aspects of function in the cytochrome P450 monooxygenase superfamily. These heme-containing enzymes typically catalyze the oxidation of unactivated C–H and C=C bonds in a multitude of substrates, often with complete regio- and stereospecificity. Over 600 000 genes in GenBank have been assigned to P450s, yet all known P450 structures exhibit a highly conserved and unique fold. This combination of functional and structural conservation with a vast substrate clientele, each substrate having multiple possible sites for oxidation, makes the P450s a unique target for understanding the role of enzyme structure and dynamics in determining a particular substrate–product combination. P450s are large by solution NMR standards, requiring us to develop specialized approaches for making sequential resonance assignments and interpreting the spectral changes that occur as a function of changing conditions (e.g., oxidation and spin state changes, ligand, substrate or effector binding). Solution conformations are characterized by the fitting of residual dipolar couplings (RDCs) measured for sequence-specifically assigned amide N–H correlations to alignment tensors optimized in the course of restrained molecular dynamics (MD) simulations. The conformational ensembles obtained by such RDC-restrained simulations, which we call “soft annealing”, are then tested by site-directed mutation and spectroscopic and activity assays for relevance. These efforts have gained us insights into cryptic conformational changes associated with substrate and redox partner binding that were not suspected from crystal structures, but were shown by subsequent work to be relevant to function. Furthermore, it appears that many of these changes can be generalized to P450s besides those that we have characterized, providing guidance for enzyme engineering efforts. While past research was primarily directed at the more tractable prokaryotic P450s, our current efforts are aimed at medically relevant human enzymes, including CYP17A1, CYP2D6, and CYP3A4.

The characterisation of two members of the cytochrome P450 CYP150 family: CYP150A5 and CYP150A6 from Mycobacterium marinum

BACKGROUND

Actinobacteria, including the Mycobacteria, have a large component of cytochrome P450 family monooxygenases. This includes Mycobacterium tuberculosis, M. ulcerans and M. marinum, and M. vanbaalenii. These enzymes can abstract CH bonds and have important roles in natural product biosynthesis.

METHODS

Two members of the bacterial CYP150 family, CYP150A5 and CYP150A6 from M. marinum, were produced, purified and characterised. The potential substrate ranges of both enzymes were analysed and the monooxygenase activity of CYP150A5 was reconstituted using a physiological electron transfer partner system. CYP150A6 was structurally characterised by X-ray crystallography.

RESULTS

CYP150A5 was shown to bind various norisoprenoids and terpenoids. It could regioselectively hydroxylate β-ionol. The X-ray crystal structure of substrate-free CYP150A6 was solved to 1.5 Å. This displayed an open conformation with short F and G helices, an unresolved F-G loop region and exposed active site pocket. The active site residues could be identified and important variations were found among the CYP150A enzymes. Haem-binding azole inhibitors were identified for both enzymes.

CONCLUSIONS

The structure of CYP150A6 will facilitate the identification of physiological substrates and the design of better inhibitors for members of this P450 family. Based on the observed differences in substrate binding preference and sequence variations among the active site residues, their roles are predicted to be different.

GENERAL SIGNIFICANCE

Multiple CYP150 family members were found in many bacteria and are prevalent in the Mycobacteria including several human pathogens. Inhibition and structural data are reported here for these enzymes for the first time.

Identification of Genes Encoding Enzymes Catalyzing the Early Steps of Carrot Polyacetylene Biosynthesis

Polyacetylenic lipids accumulate in various Apiaceae species after pathogen attack, suggesting that these compounds are naturally occurring pesticides and potentially valuable resources for crop improvement. These compounds also promote human health and slow tumor growth. Even though polyacetylenic lipids were discovered decades ago, the biosynthetic pathway underlying their production is largely unknown. To begin filling this gap and ultimately enable polyacetylene engineering, we studied polyacetylenes and their biosynthesis in the major Apiaceae crop carrot (Daucus carota subsp. sativus). Using gas chromatography and mass spectrometry, we identified three known polyacetylenes and assigned provisional structures to two novel polyacetylenes. We also quantified these compounds in carrot leaf, petiole, root xylem, root phloem, and root periderm extracts. Falcarindiol and falcarinol predominated and accumulated primarily in the root periderm. Since the multiple double and triple carbon-carbon bonds that distinguish polyacetylenes from ubiquitous fatty acids are often introduced by Δ12 oleic acid desaturase (FAD2)-type enzymes, we mined the carrot genome for FAD2 genes. We identified a FAD2 family with an unprecedented 24 members and analyzed public, tissue-specific carrot RNA-Seq data to identify coexpressed members with root periderm-enhanced expression. Six candidate genes were heterologously expressed individually and in combination in yeast and Arabidopsis (Arabidopsis thaliana), resulting in the identification of one canonical FAD2 that converts oleic to linoleic acid, three divergent FAD2-like acetylenases that convert linoleic into crepenynic acid, and two bifunctional FAD2s with Δ12 and Δ14 desaturase activity that convert crepenynic into the further desaturated dehydrocrepenynic acid, a polyacetylene pathway intermediate. These genes can now be used as a basis for discovering other steps of falcarin-type polyacetylene biosynthesis, to modulate polyacetylene levels in plants, and to test the in planta function of these molecules.

Entropic contribution to enhanced thermal stability in the thermostable P450 CYP119

The enhanced thermostability of thermophilic proteins with respect to their mesophilic counterparts is often attributed to the enthalpy effect, arising from strong interactions between protein residues. Intuitively, these strong interresidue interactions will rigidify the biomolecules. However, the present work utilizing neutron scattering and solution NMR spectroscopy measurements demonstrates a contrary example that the thermophilic cytochrome P450, CYP119, is much more flexible than its mesophilic counterpart, CYP101A1, something which is not apparent just from structural comparison of the two proteins. A mechanism to explain this apparent contradiction is that higher flexibility in the folded state of CYP119 increases its conformational entropy and thereby reduces the entropy gain during denaturation, which will increase the free energy needed for unfolding and thus stabilize the protein. This scenario is supported by thermodynamic data on the temperature dependence of unfolding free energy, which shows a significant entropic contribution to the thermostability of CYP119 and lends an added dimension to enhanced stability, previously attributed only to presence of aromatic stacking interactions and salt bridge networks. Our experimental data also support the notion that highly thermophilic P450s such as CYP119 may use a mechanism that partitions flexibility differently from mesophilic P450s between ligand binding and thermal stability.

Substrate-induced conformational change in cytochrome P450 OleP

The regulation of cytochrome P450 activity is often achieved by structural transitions induced by substrate binding. We describe the conformational transition experienced upon binding by the P450 OleP, an epoxygenase involved in oleandomycin biosynthesis. OleP bound to the substrate analog 6DEB crystallized in 2 forms: one with an ensemble of open and closed conformations in the asymmetric unit and another with only the closed conformation. Characterization of OleP-6DEB binding kinetics, also using the P450 inhibitor clotrimazole, unveiled a complex binding mechanism that involves slow conformational rearrangement with the accumulation of a spectroscopically detectable intermediate where 6DEB is bound to open OleP. Data reported herein provide structural snapshots of key precatalytic steps in the OleP reaction and explain how structural rearrangements induced by substrate binding regulate activity.-Parisi, G., Montemiglio, L. C., Giuffrè, A., Macone, A., Scaglione, A., Cerutti, G., Exertier, C., Savino, C., Vallone, B. Substrate-induced conformational change in cytochrome P450 OleP.

A new approach to understanding structure-function relationships in cytochromes P450 by targeting terpene metabolism in the wild

A strategy for elucidating sequence determinants of function in the class of cytochrome P450 (CYP) enzymes that catalyze the first steps of terpene metabolism in wild microbiomes is described. Wild organisms that can use camphor, terpineol, pinene and limonene were isolated from soils rich in coniferous waste. Cell free extracts and growth beers were analyzed by gas chromatography/mass spectrometry to identify primary oxidative metabolites. For one organism, Pseudomonas nitroreducens TPJM, a cytochrome P450 (CYP108B1) isolated from cell free extracts was demonstrated to catalyze the oxidation of α-terpineol in assays combining the native ferredoxin and putidaredoxin reductase, and the resulting oxidation products identified by gas chromatography/mass spectrometry. Shotgun sequencing of PnTPJM identified four candidate P450 genes, including an apparently fragmentary gene with a high degree of homology with the known enzyme CYP108 (P450_terp).

The crystal structure of P450-TT heme-domain provides the first structural insights into the versatile class VII P450s

The first crystal structure of a class VII P450, CYP116B46 from Tepidiphilus thermophilus, has been solved at 1.9 Å resolution. The structure reveals overall conservation of the P450-fold and a water conduit around the I-helix. Active site residues have been identified and sequence comparisons have been made with other class VII enzymes. A structure similarity search demonstrated that the P450-TT structure is similar to enzymes capable of oxy-functionalization of fatty acids, terpenes, macrolides, steroids and statins. The insight gained from solving this structure will provide a guideline for future engineering and modelling studies on this catalytically promiscuous class of enzymes.

[Methods for determining of cytochrome P450 isozymes functional activity]

The review is dedicated to modern methods and technologies for determining of cytochrome P450 isozymes functional activity, such as absorbance and fluorescent spectroscopy, electron paramagnetic resonance (EPR), nuclear magnetic resonance (NMR), Raman, Mossbauer, and X-ray spectroscopy, surface plasmon resonance (SPR), atomic force microscopy (AFM). Methods of molecular genetic analysis were reviewed from personalized medicine point of view. The use of chromate-mass-spectrometric methods for cytochrome P450-dependent catalytic reactions' products was discussed. The review covers modern electrochemical systems based on cytochrome P450 isozymes for their catalytic activity analysis, their use in practice and further development perspectives for experimental pharmacology, biotechnology and translational medicine.

Some Surprising Implications of NMR-directed Simulations of Substrate Recognition and Binding by Cytochrome P450cam (CYP101A1)

Cytochrome P450_cam (CYP101A1) catalyzes the stereospecific 5-exo hydroxylation of d-camphor by molecular oxygen. Previously, residual dipolar couplings measured for backbone amide ¹H-¹⁵N correlations in both substrate-free and bound forms of CYP101A1 were used as restraints in soft annealing molecular dynamic simulations in order to identify average conformations of the enzyme with and without substrate bound. Multiple substrate-dependent conformational changes remote from the enzyme active site were identified, and site-directed mutagenesis and activity assays confirmed the importance of these changes in substrate recognition. The current work makes use of perturbation response scanning (PRS) and umbrella sampling molecular dynamic of the residual dipolar coupling-derived CYP101A1 structures to probe the roles of remote structural features in enforcing the regio- and stereospecific nature of the hydroxylation reaction catalyzed by CYP101A1. An improper dihedral angle Ψ was defined and used to maintain substrate orientation in the CYP101A1 active site, and it was observed that different values of Ψ result in different PRS response maps. Umbrella sampling methods show that the free energy of the system is sensitive to Ψ, and bound substrate forms an important mechanical link in the transmission of mechanical coupling through the enzyme structure. Finally, a qualitative approach to interpreting PRS maps in terms of the roles of secondary structural features is proposed.

Spectroscopic studies of the cytochrome P450 reaction mechanisms

The cytochrome P450 monooxygenases (P450s) are thiolate heme proteins that can, often under physiological conditions, catalyze many distinct oxidative transformations on a wide variety of molecules, including relatively simple alkanes or fatty acids, as well as more complex compounds such as steroids and exogenous pollutants. They perform such impressive chemistry utilizing a sophisticated catalytic cycle that involves a series of consecutive chemical transformations of heme prosthetic group. Each of these steps provides a unique spectral signature that reflects changes in oxidation or spin states, deformation of the porphyrin ring or alteration of dioxygen moieties. For a long time, the focus of cytochrome P450 research was to understand the underlying reaction mechanism of each enzymatic step, with the biggest challenge being identification and characterization of the powerful oxidizing intermediates. Spectroscopic methods, such as electronic absorption (UV-Vis), electron paramagnetic resonance (EPR), nuclear magnetic resonance (NMR), electron nuclear double resonance (ENDOR), Mössbauer, X-ray absorption (XAS), and resonance Raman (rR), have been useful tools in providing multifaceted and detailed mechanistic insights into the biophysics and biochemistry of these fascinating enzymes. The combination of spectroscopic techniques with novel approaches, such as cryoreduction and Nanodisc technology, allowed for generation, trapping and characterizing long sought transient intermediates, a task that has been difficult to achieve using other methods. Results obtained from the UV-Vis, rR and EPR spectroscopies are the main focus of this review, while the remaining spectroscopic techniques are briefly summarized. This article is part of a Special Issue entitled: Cytochrome P450 biodiversity and biotechnology, edited by Erika Plettner, Gianfranco Gilardi, Luet Wong, Vlada Urlacher, Jared Goldstone.

Substrate recognition by two different P450s: Evidence for conserved roles in a common fold

Cytochrome P450 monooxygenases CYP101A1 and MycG catalyze regio- and stereospecific oxidations of their respective substrates, d-camphor and mycinamicin IV. Despite the low sequence homology between the two enzymes (29% identity) and differences in size and hydrophobicity of their substrates, the conformational changes that occur upon substrate binding in both enzymes as determined by solution NMR methods show some striking similarities. Many of the same secondary structural features in both enzymes are perturbed, suggesting the existence of a common mechanism for substrate binding and recognition in the P450 superfamily.

Molecular cloning and 3D model of first cytochrome P450 from CYP3A subfamily in saltwater crocodile (Crocodylus porosus)

Cytochrome P450s (CYPs) play critical role in oxidative metabolism of numerous xenobiotics and endogenous compounds. The first CYP3A subfamily member in saltwater crocodile has been cloned and modelled for three-dimensional (3D) structure. The full-length cDNA was obtained employing reverse transcription polymerase chain reaction (RT-PCR) strategy and rapid amplification of cDNA ends (RACE). The cDNA sequence of 1659 nucleotides includes 132 nucleotides from 5' untranslated region (UTR), an open reading frame of 1527 nucleotides encoding 509 amino acids designated as CYP3A163. The alignment of CYP3A163 sequence with CYP3A subfamily across the lineages exhibit the loss of 1 residue in birds and 7 residues in mammals in comparison to reptiles suggesting the adaptation processes during evolution. The amino acid identity of CYP3A163 with Alligator mississippiensis CYP3A77 and Homo sapiens CYP3A4 is 91% and 62% respectively. The 3D structure of CYP3A163 modelled using human CYP3A4 structure as a template with Phyre² software, represents high similarity with its functionally important motifs and catalytic domain. Both sequence and structure of CYP3A163 display the common and conserved features of CYP3A subfamily. Overall, this study provides primary molecular and structural data of CYP3A163 required to investigate the xenobiotic metabolism in saltwater crocodiles.

NADH reduction of nitroaromatics as a probe for residual ferric form high-spin in a cytochrome P450

The existence of a substrate-sensitive equilibrium between high spin (S=5/2) and low spin (S=1/2) ferric iron is a well-established phenomenon in the cytochrome P450 (CYP) superfamily, although its origins are still a subject of discussion. A series of mutations that strongly perturb the spin state equilibrium in the camphor hydroxylase CYP101A1 were recently described (Colthart et al., Sci. Rep. 6, 22035 (2016)). Wild type CYP101A1 as well as some CYP101A1 mutants are herein shown to be capable of catalyzing the reduction of nitroacetophenones by NADH to the corresponding anilino compounds (nitroreductase or NRase activity). The distinguishing characteristic between those mutants that catalyze the reduction and those that cannot appears to be the extent to which residual high spin form exists in the absence of the native substrate d-camphor, with those showing the largest spin state shifts upon camphor binding also exhibiting NRase activity. Optical and EPR spectroscopy was used to further examine these phenomena. These results suggest that reduction of nitroaromatics may provide a useful probe of residual high spin states in the CYP superfamily. This article is part of a Special Issue entitled: Cytochrome P450 biodiversity and biotechnology, edited by Erika Plettner, Gianfranco Gilardi, Luet Wong, Vlada Urlacher, Jared Goldstone.

Structural insights into the function of steroidogenic cytochrome P450 17A1

Cytochrome P450 17A1 (CYP17A1) operates at the core of human steroidogenesis, directing precursors into mineralocorticoids, glucocorticoids, or sex steroids. Although the 17α-hydroxylase and 17,20-lyase activities of this dual function enzyme have been investigated extensively, until recently no CYP17A1 structures were available to inform our understanding. Structures of CYP17A1 with a range of steroidal inhibitors and substrates are now available. This review relates functional knowledge of this enzyme to structural features defining the selective differentiation between its various substrates. While both hydroxylase and lyase substrates have similar orientations with respect to the heme, subtle differences in hydrogen bonding between CYP17A1 and the C3 substituent at the opposite end of ligands appear to correlate with differential substrate utilization and product formation. Complementary structural information from solution NMR supports cytochrome b₅ allosteric modulation of the lyase reaction, implicating regions involved in ligand access to the otherwise buried active site.

Physical Studies of P450-P450 Interactions: Predicting Quaternary Structures of P450 Complexes in Membranes from Their X-ray Crystal Structures

Cytochrome P450 enzymes, which catalyze oxygenation reactions of both exogenous and endogenous chemicals, are membrane bound proteins that require interaction with their redox partners in order to function. Those responsible for drug and foreign compound metabolism are localized primarily in the endoplasmic reticulum of liver, lung, intestine, and other tissues. More recently, the potential for P450 enzymes to exist as supramolecular complexes has been shown by the demonstration of both homomeric and heteromeric complexes. The P450 units in these complexes are heterogeneous with respect to their distribution and function, and the interaction of different P450s can influence P450-specific metabolism. The goal of this review is to examine the evidence supporting the existence of physical complexes among P450 enzymes. Additionally, the review examines the crystal lattices of different P450 enzymes derived from X-ray diffraction data to make assumptions regarding possible quaternary structures in membranes and in turn, to predict how the quaternary structures could influence metabolism and explain the functional effects of specific P450-P450 interactions.

Conformational Mobility in Cytochrome P450 3A4 Explored by Pressure-Perturbation EPR Spectroscopy

We used high hydrostatic pressure as a tool for exploring the conformational landscape of human cytochrome P450 3A4 (CYP3A4) by electron paramagnetic resonance and fluorescence spectroscopy. Site-directed incorporation of a luminescence resonance energy transfer donor-acceptor pair allowed us to identify a pressure-dependent equilibrium between two states of the enzyme, where an increase in pressure increased the spatial separation between the two distantly located fluorophores. This transition is characterized by volume change (ΔV°) and P1/2 values of -36.8 ± 5.0 mL/mol and 1.45 ± 0.33 kbar, respectively, which corresponds to a Keq° of 0.13 ± 0.06, so that only 15% of the enzyme adopts the pressure-promoted conformation at ambient pressure. This pressure-promoted displacement of the equilibrium is eliminated by the addition of testosterone, an allosteric activator. Using site-directed spin labeling, we demonstrated that the pressure- and testosterone-sensitive transition is also revealed by pressure-induced changes in the electron paramagnetic resonance spectra of a nitroxide side chain placed at position 85 or 409 of the enzyme. Furthermore, we observed a pressure-induced displacement of the emission maxima of a solvatochromic fluorophore (7-diethylamino-3-((((2-maleimidyl)ethyl)amino)carbonyl) coumarin) placed at the same positions, which suggests a relocation to a more polar environment. Taken together, the results reveal an effector-dependent conformational equilibrium between open and closed states of CYP3A4 that involves a pronounced change at the interface between the region of α-helices A/A' and the meander loop of the enzyme, where residues 85 and 409 are located. Our study demonstrates the high potential of pressure-perturbation strategies for studying protein conformational landscapes.

Whole-Genome Identification, Phylogeny, and Evolution of the Cytochrome P450 Family 2 (CYP2) Subfamilies in Birds

The cytochrome P450 (CYP) superfamily defends organisms from endogenous and noxious environmental compounds, and thus is crucial for survival. However, beyond mammals the molecular evolution of CYP2 subfamilies is poorly understood. Here, we characterized the CYP2 family across 48 avian whole genomes representing all major extant bird clades. Overall, 12 CYP2 subfamilies were identified, including the first description of the CYP2F, CYP2G, and several CYP2AF genes in avian genomes. Some of the CYP2 genes previously described as being lineage-specific, such as CYP2K and CYP2W, are ubiquitous to all avian groups. Furthermore, we identified a large number of CYP2J copies, which have been associated previously with water reabsorption. We detected positive selection in the avian CYP2C, CYP2D, CYP2H, CYP2J, CYP2K, and CYP2AC subfamilies. Moreover, we identified new substrate recognition sites (SRS0, SRS2_SRS3, and SRS3.1) and heme binding areas that influence CYP2 structure and function of functional importance as under significant positive selection. Some of the positively selected sites in avian CYP2D are located within the same SRS1 region that was previously linked with the metabolism of plant toxins. Additionally, we find that selective constraint variations in some avian CYP2 subfamilies are consistently associated with different feeding habits (CYP2H and CYP2J), habitats (CYP2D, CYP2H, CYP2J, and CYP2K), and migratory behaviors (CYP2D, CYP2H, and CYP2J). Overall, our findings indicate that there has been active enzyme site selection on CYP2 subfamilies and differential selection associated with different life history traits among birds.

The Role of Protein-Protein and Protein-Membrane Interactions on P450 Function

This symposium summary, sponsored by the ASPET, was held at Experimental Biology 2015 on March 29, 2015, in Boston, Massachusetts. The symposium focused on: 1) the interactions of cytochrome P450s (P450s) with their redox partners; and 2) the role of the lipid membrane in their orientation and stabilization. Two presentations discussed the interactions of P450s with NADPH-P450 reductase (CPR) and cytochrome b5. First, solution nuclear magnetic resonance was used to compare the protein interactions that facilitated either the hydroxylase or lyase activities of CYP17A1. The lyase interaction was stimulated by the presence of b5 and 17α-hydroxypregnenolone, whereas the hydroxylase reaction was predominant in the absence of b5. The role of b5 was also shown in vivo by selective hepatic knockout of b5 from mice expressing CYP3A4 and CYP2D6; the lack of b5 caused a decrease in the clearance of several substrates. The role of the membrane on P450 orientation was examined using computational methods, showing that the proximal region of the P450 molecule faced the aqueous phase. The distal region, containing the substrate-access channel, was associated with the membrane. The interaction of NADPH-P450 reductase (CPR) with the membrane was also described, showing the ability of CPR to "helicopter" above the membrane. Finally, the endoplasmic reticulum (ER) was shown to be heterogeneous, having ordered membrane regions containing cholesterol and more disordered regions. Interestingly, two closely related P450s, CYP1A1 and CYP1A2, resided in different regions of the ER. The structural characteristics of their localization were examined. These studies emphasize the importance of P450 protein organization to their function.