Cytochrome P450 flexibility

Biocatalytic, Enantioenriched Primary Amination of Tertiary C-H Bonds

Intermolecular functionalization of tertiary C-H bonds to construct fully substituted stereogenic carbon centers represents a formidable challenge: without the assistance of directing groups, state-of-the-art catalysts struggle to introduce chirality to racemic tertiary sp 3 -carbon centers. Direct asymmetric functionalization of such centers is a worthy reactivity and selectivity goal for modern biocatalysis. Here we present an engineered nitrene transferase (P411-TEA-5274), derived from a bacterial cytochrome P450, that is capable of aminating tertiary C-H bonds to provide chiral α-tertiary primary amines with high efficiency (up to 2300 total turnovers) and selectivity (up to >99% enantiomeric excess (e.e.)). The construction of fully substituted stereocenters with methyl and ethyl groups underscores the enzyme's remarkable selectivity. A comprehensive substrate scope study demonstrates the biocatalyst's compatibility with diverse functional groups and tertiary C-H bonds. Mechanistic studies elucidate how active-site residues distinguish between the enantiomers and enable the enzyme to perform this transformation with excellent enantioselectivity.

Substrate Conformational Switch Enables the Stereoselective Dimerization in P450 NascB: Insights from Molecular Dynamics Simulations and Quantum Mechanical/Molecular Mechanical Calculations

P450 NascB catalyzes the coupling of cyclo-(l-tryptophan-l-proline) (1) to generate (-)-naseseazine C (2) through intramolecular C-N bond formation and intermolecular C-C coupling. A thorough understanding of its catalytic mechanism is crucial for the engineering or design of P450-catalyzed C-N dimerization reactions. By employing MD simulations, QM/MM calculations, and enhanced sampling, we assessed various mechanisms from recent works. Our study demonstrates that the most favorable pathway entails the transfer of a hydrogen atom from N7-H to Cpd I. Subsequently, there is a conformational change in the substrate radical, shifting it from the Re-face to the Si-face of N7 in Substrate 1. The Si-face conformation of Substrate 1 is stabilized by the protein environment and the π-π stacking interaction between the indole ring and heme porphyrin. The subsequent intermolecular C3-C6' bond formation between Substrate 1 radical and Substrate 2 occurs via a radical attack mechanism. The conformational switch of the Substrate 1 radical not only lowers the barrier of the intermolecular C3-C6' bond formation but also yields the correct stereoselectivity observed in experiments. In addition, we evaluated the reactivity of the ferric-superoxide species, showing it is not reactive enough to initiate the hydrogen atom abstraction from the indole NH group of the substrate. Our simulation provides a comprehensive mechanistic insight into how the P450 enzyme precisely controls both the intramolecular C-N cyclization and intermolecular C-C coupling. The current findings align with the available experimental data, emphasizing the pivotal role of substrate dynamics in governing P450 catalysis.

Multilevel superposition for deciphering the conformational variability of protein ensembles

The dynamics and variability of protein conformations are directly linked to their functions. Many comparative studies of X-ray protein structures have been conducted to elucidate the relevant conformational changes, dynamics and heterogeneity. The rapid increase in the number of experimentally determined structures has made comparison an effective tool for investigating protein structures. For example, it is now possible to compare structural ensembles formed by enzyme species, variants or the type of ligands bound to them. In this study, the author developed a multilevel model for estimating two covariance matrices that represent inter- and intra-ensemble variability in the Cartesian coordinate space. Principal component analysis using the two estimated covariance matrices identified the inter-/intra-enzyme variabilities, which seemed to be important for the enzyme functions, with the illustrative examples of cytochrome P450 family 2 enzymes and class A $\beta$-lactamases. In P450, in which each enzyme has its own active site of a distinct size, an active-site motion shared universally between the enzymes was captured as the first principal mode of the intra-enzyme covariance matrix. In this case, the method was useful for understanding the conformational variability after adjusting for the differences between enzyme sizes. The developed method is advantageous in small ensemble-size problems and hence promising for use in comparative studies on experimentally determined structures where ensemble sizes are smaller than those generated, for example, by molecular dynamics simulations.

Asparagine-85 Stabilizes a Structural Active Site Water Network in CYP121A1 of Mycobacterium tuberculosis

The cytochrome P450 enzyme CYP121A1 endogenously catalyzes the formation of a carbon-carbon bond between the two phenol groups of dicyclotyrosine (cYY) in Mycobacterium tuberculosis (Mtb). One of 20 CYP enzymes in Mtb, CYP121A1 continues to garner significant interest as a potential drug target. The accompanying reports the use of 19F NMR spectroscopy, reconstituted activity assays, and molecular dynamics simulations to investigate the significance of hydrogen bonding interactions that were theorized to stabilize a static active site water network. The active site residue Asn-85, whose hydrogen bonds with the diketopiperazine ring of cYY contributes to a contiguous active site water network in the absence of cYY, was mutated to a serine (N85S) and to a glutamine (N85Q). These conservative changes in the hydrogen bond donor side chain result in inactivation of the enzyme. Moreover, the N85S mutation induces reverse type-I binding as measured by absorbance difference spectra. NMR spectra monitoring the ligand-adaptive FG-loop and the active site Trp-182 side chain confirm that disruption of the active site water network also significantly alters the structure of the active site. These data were consistent with dynamics simulations of N85S and N85Q that demonstrate that a compromised water network is responsible for remodeling of the active site B-helix and a repositioning of cYY toward the heme. These findings implicate a slowly exchanging water network as a critical factor in CYP121A1 function and a likely contributor to the unusual rigidity of the structure.

Enzymatic Assembly of Diverse Lactone Structures: An Intramolecular C-H Functionalization Strategy

Lactones are cyclic esters with extensive applications in materials science, medicinal chemistry, and the food and perfume industries. Nature's strategy for the synthesis of many lactones found in natural products always relies on a single type of retrosynthetic strategy, a C-O bond disconnection. Here, we describe a set of laboratory-engineered enzymes that use a new-to-nature C-C bond-forming strategy to assemble diverse lactone structures. These engineered "carbene transferases" catalyze intramolecular carbene insertions into benzylic or allylic C-H bonds, which allow for the synthesis of lactones with different ring sizes and ring scaffolds from simple starting materials. Starting from a serine-ligated cytochrome P450 variant previously engineered for other carbene-transfer activities, directed evolution generated a variant P411-LAS-5247, which exhibits a high activity for constructing a five-membered ε-lactone, lactam, and cyclic ketone products (up to 5600 total turnovers (TTN) and >99% enantiomeric excess (ee)). Further engineering led to variants P411-LAS-5249 and P411-LAS-5264, which deliver six-membered δ-lactones and seven-membered ε-lactones, respectively, overcoming the thermodynamically unfavorable ring strain associated with these products compared to the γ-lactones. This new carbene-transfer activity was further extended to the synthesis of complex lactone scaffolds based on fused, bridged, and spiro rings. The enzymatic platform developed here complements natural biosynthetic strategies for lactone assembly and expands the structural diversity of lactones accessible through C-H functionalization.

Enantio- and Diastereoenriched Enzymatic Synthesis of 1,2,3-Polysubstituted Cyclopropanes from (Z/E)-Trisubstituted Enol Acetates

In nature and synthetic chemistry, stereoselective [2 + 1] cyclopropanation is the most prevalent strategy for the synthesis of chiral cyclopropanes, a class of key pharmacophores in pharmaceuticals and bioactive natural products. One of the most extensively studied reactions in the organic chemist's arsenal, stereoselective [2 + 1] cyclopropanation, largely relies on the use of stereodefined olefins, which can require elaborate laboratory synthesis or tedious separation to ensure high stereoselectivity. Here, we report engineered hemoproteins derived from a bacterial cytochrome P450 that catalyze the synthesis of chiral 1,2,3-polysubstituted cyclopropanes, regardless of the stereopurity of the olefin substrates used. Cytochrome P450BM3 variant P411-INC-5185 exclusively converts (Z)-enol acetates to enantio- and diastereoenriched cyclopropanes and in the model reaction delivers a leftover (E)-enol acetate with 98% stereopurity, using whole Escherichia coli cells. P411-INC-5185 was further engineered with a single mutation to enable the biotransformation of (E)-enol acetates to α-branched ketones with high levels of enantioselectivity while simultaneously catalyzing the cyclopropanation of (Z)-enol acetates with excellent activities and selectivities. We conducted docking studies and molecular dynamics simulations to understand how active-site residues distinguish between the substrate isomers and enable the enzyme to perform these distinct transformations with such high selectivities. Computational studies suggest the observed enantio- and diastereoselectivities are achieved through a stepwise pathway. These biotransformations streamline the synthesis of chiral 1,2,3-polysubstituted cyclopropanes from readily available mixtures of (Z/E)-olefins, adding a new dimension to classical cyclopropanation methods.

Computation-Aided Engineering of Cytochrome P450 for the Production of Pravastatin

CYP105AS1 is a cytochrome P450 from Amycolatopsis orientalis that catalyzes monooxygenation of compactin to 6-epi-pravastatin. For fermentative production of the cholesterol-lowering drug pravastatin, the stereoselectivity of the enzyme needs to be inverted, which has been partially achieved by error-prone PCR mutagenesis and screening. In the current study, we report further optimization of the stereoselectivity by a computationally aided approach. Using the CoupledMoves protocol of Rosetta, a virtual library of mutants was designed to bind compactin in a pro-pravastatin orientation. By examining the frequency of occurrence of beneficial substitutions and rational inspection of their interactions, a small set of eight mutants was predicted to show the desired selectivity and these variants were tested experimentally. The best CYP105AS1 variant gave >99% stereoselective hydroxylation of compactin to pravastatin, with complete elimination of the unwanted 6-epi-pravastatin diastereomer. The enzyme-substrate complexes were also examined by ultrashort molecular dynamics simulations of 50 × 100 ps and 5 × 22 ns, which revealed that the frequency of occurrence of near-attack conformations agreed with the experimentally observed stereoselectivity. These results show that a combination of computational methods and rational inspection could improve CYP105AS1 stereoselectivity beyond what was obtained by directed evolution. Moreover, the work lays out a general in silico framework for specificity engineering of enzymes of known structure.

Role of P34S, G169R, R296C, and S486T Substitutions in Ligand Access and Catalysis for Cytochrome P450 2D6 Allelic Variants CYP2D6*14A and CYP2D6*14B

BACKGROUND

Genetic polymorphism of cytochrome P450 (CYP) contributes to variability in drug metabolism, clearance, and response. This study aimed to investigate the functional and molecular basis for altered ligand binding and catalysis in CYP2D6*14A and CYP2D6*14B, two unique alleles common in the Asian population.

METHODS

CYP proteins expressed in Escherichia coli were studied using the substrate 3-cyano-7- ethoxycoumarin (CEC) and inhibitor probes (quinidine, fluoxetine, paroxetine, terbinafine) in the enzyme assay. Computer modelling was additionally used to create three-dimensional structures of the CYP2D6*14 variants.

RESULTS

Kinetics data indicated significantly reduced intrinsic clearance in CYP2D6*14 variants, suggesting that P34S, G169R, R296C, and S486T substitutions worked cooperatively to alter the conformation of the active site that negatively impacted the deethylase activity of CYP2D6. For the inhibition studies, IC50 values decreased in quinidine, paroxetine, and terbinafine but increased in fluoxetine, suggesting a varied ligand-specific susceptibility to inhibition. Molecular docking further demonstrated the role of P34S and R296C in altering access channel dimensions, thereby affecting ligand access and binding and subsequently resulting in varied inhibition potencies.

CONCLUSION

In summary, the differential selectivity of CYP2D6*14 variants for the ligands (substrate and inhibitor) was governed by the alteration of the active site and access channel architecture induced by the natural mutations found in the alleles.

Mechanistic Insights into the P450 TleB-Catalyzed Unusual Intramolecular C-N Bond Formation Involved in the Biosynthesis of Indolactam V

Indolactam V, a known biosynthetic precursor of indolactam alkaloids, is the main pharmacophore of this family and exhibits potential protein kinase C activation. A key step in the biosynthesis of indolactam V is the formation of an indole-fused nine-membered lactam core by intramolecular C-N bond formation. In this work, we report a computational study of the unique cytochrome P450 TleB enzyme-catalyzed direct and selective C-H bond amination reaction that can generate indolactam V from the dipeptide N-methylvalyl-tryptophanol. By performing molecular dynamics simulations and quantum-mechanical/molecular-mechanical calculations, we revealed that the C-H bond amination involves one step of proton transfer from N1-H of the indole ring to the Fe^IV═O unit, one step of hydrogen abstraction of N13-H in the side chain of the substrate by the Fe^IV-OH unit, and diradical coupling, in which two conformational changes of the side chain of the substrate are necessary. In the enzyme-substrate complex of TleB, the N-H bond of the indole ring of the substrate forms a strong hydrogen bond with the Fe^IV═O unit in compound I, and the porphyrin radical cation accepts an electron from the substrate to form the closed-shell electronic configuration. Thus, compound I in the enzyme-substrate complex cannot be described as Fe^IV═O coupled to a porphyrin radical cation, which is different from those of other P450 enzymes. Besides, two stages of conformational changes of the side chains of the substrate may increase the relative energies of reaction intermediates by 10-12 kcal/mol. From the structure point of view, it is the rotatable long side chain of the substrate and the large flexible active pocket of TleB that make the intramolecular diradical coupling feasible. Our findings may provide useful information to further understand the Tleb-catalyzed intramolecular C-H bond amination and the other bio-catalyzed intramolecular diradical coupling.

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.

Navigating the Unnatural Reaction Space: Directed Evolution of Heme Proteins for Selective Carbene and Nitrene Transfer

Despite the astonishing diversity of naturally occurring biocatalytic processes, enzymes do not catalyze many of the transformations favored by synthetic chemists. Either nature does not care about the specific products, or if she does, she has adopted a different synthetic strategy. In many cases, the appropriate reagents used by synthetic chemists are not readily accessible to biological systems. Here, we discuss our efforts to expand the catalytic repertoire of enzymes to encompass powerful reactions previously known only in small-molecule catalysis: formation and transfer of reactive carbene and nitrene intermediates leading to a broad range of products, including products with bonds not known in biology. In light of the structural similarity of iron carbene (Fe=C(R¹)(R²)) and iron nitrene (Fe=NR) to the iron oxo (Fe=O) intermediate involved in cytochrome P450-catalyzed oxidation, we have used synthetic carbene and nitrene precursors that biological systems have not encountered and repurposed P450s to catalyze reactions that are not known in the natural world. The resulting protein catalysts are fully genetically encoded and function in intact microbial cells or cell-free lysates, where their performance can be improved and optimized by directed evolution. By leveraging the catalytic promiscuity of P450 enzymes, we evolved a range of carbene and nitrene transferases exhibiting excellent activity toward these new-to-nature reactions. Since our initial report in 2012, a number of other heme proteins including myoglobins, protoglobins, and cytochromes c have also been found and engineered to promote unnatural carbene and nitrene transfer. Due to the altered active-site environments, these heme proteins often displayed complementary activities and selectivities to P450s.

Using wild-type and engineered heme proteins, we and others have described a range of selective carbene transfer reactions, including cyclopropanation, cyclopropenation, Si–H insertion, B–H insertion, and C–H insertion. Similarly, a variety of asymmetric nitrene transfer processes including aziridination, sulfide imidation, C–H amidation, and, most recently, C–H amination have been demonstrated. The scopes of these biocatalytic carbene and nitrene transfer reactions are often complementary to the state-of-the-art processes based on small-molecule transition-metal catalysts, making engineered biocatalysts a valuable addition to the synthetic chemist’s toolbox. Moreover, enabled by the exquisite regio- and stereocontrol imposed by the enzyme catalyst, this biocatalytic platform provides an exciting opportunity to address challenging problems in modern synthetic chemistry and selective catalysis, including ones that have eluded synthetic chemists for decades.

Towards resolving the complex paramagnetic nuclear magnetic resonance (NMR) spectrum of small laccase: assignments of resonances to residue-specific nuclei

Laccases efficiently reduce dioxygen to water in an active site containing a tri-nuclear copper centre (TNC). The dynamics of the protein matrix is a determining factor in the efficiency in catalysis. To probe mobility, nuclear magnetic resonance (NMR) spectroscopy is highly suitable. However, several factors complicate the assignment of resonances to active site nuclei in laccases. The paramagnetic nature causes large shifts and line broadening. Furthermore, the presence of slow chemical exchange processes of the imidazole rings of copper ligand results in peak doubling. A third complicating factor is that the enzyme occurs in two states, the native intermediate (NI) and resting oxidized (RO) states, with different paramagnetic properties. The present study aims at resolving the complex paramagnetic NMR spectra of the TNC of Streptomyces coelicolor small laccase (SLAC). With a combination of paramagnetically tailored NMR experiments, all eight His Nδ 1 and Hδ 1 resonances for the NI state are identified, as well as His Hβ protons for the RO state. With the help of second-shell mutagenesis, selective resonances are tentatively assigned to the histidine ligands of the copper in the type-2 site. This study demonstrates the utility of the approaches used for the sequence-specific assignment of the paramagnetic NMR spectra of ligands in the TNC that ultimately may lead to a description of the underlying motion.

Molecular and Structural Evolution of Cytochrome P450 Aromatase

Aromatase is the cytochrome P450 enzyme converting androgens into estrogen in the last phase of steroidogenesis. As estrogens are crucial in reproductive biology, aromatase is found in vertebrates and the invertebrates of the genus Branchiostoma, where it carries out the aromatization reaction of the A-ring of androgens that produces estrogens. Here, we investigate the molecular evolution of this unique and highly substrate-selective enzyme by means of structural, sequence alignment, and homology modeling, shedding light on its key role in species conservation. The alignments led to the identification of a core structure that, together with key and unique amino acids located in the active site and the substrate recognition sites, has been well conserved during evolution. Structural analysis shows what their roles are and the reason why they have been preserved. Moreover, the residues involved in the interaction with the redox partner and some phosphorylation sites appeared late during evolution. These data reveal how highly substrate-selective cytochrome P450 has evolved, indicating that the driving forces for evolution have been the optimization of the interaction with the redox partner and the introduction of phosphorylation sites that give the possibility of modulating its activity in a rapid way.

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.

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.

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.

Microsecond MD simulations of human CYP2D6 wild-type and five allelic variants reveal mechanistic insights on the function

Characterization of cytochrome P450 2D6 (CYP2D6) and the impact of the major identified allelic variants on the activity of one of the most dominating drug-metabolising enzymes is essential to increase drug safety and avoid adverse reactions. Microsecond molecular dynamics simulations have been performed to capture the dynamic signatures of this complex enzyme and five allelic variants with diverse enzymatic activity. In addition to the apo simulations, three substrates (bufuralol, veliparib and tamoxifen) and two inhibitors (prinomastat and quinidine) were included to explore their influence on the structure and dynamical features of the enzyme. Our results indicate that the altered enzyme activity can be attributed to changes in the hydrogen bonding network within the active site, and local structural differences in flexibility, position and shape of the binding pocket. In particular, the increased (CYP2D6*53) or the decreased (CYP2D6*17) activity seems to be related to a change in dynamics of mainly the BC loop due to a modified hydrogen bonding network around this region. In addition, the smallest active site volume was found for CYP2D6*4 (no activity). CYP2D6*2 (normal activity) showed no major differences in dynamic behaviour compared to the wild-type.

Catalytic modulation of human cytochromes P450 17A1 and P450 11B2 by phospholipid

Unlike most of the drug-metabolizing cytochrome P450s, microsomal P450 17A1 and mitochondrial P450 11B2 catalyze sequential multi-step reactions in steroid biosynthesis. The membrane phospholipid composition might be one parameter that modulates the efficiency and processivity of specific pathways. Here we systematically examined the effects of physiologically relevant phospholipids on the catalysis of purified P450 17A1, P450 11B2, and P450 11B1 in reconstituted assay systems. Both dioleoylphosphatidylcholine (DOPC, 18:1) and dilauroylphosphatidylcholine (DLPC, 12:0) were found to be very efficient in reconstituting 17-hydroxylase and 1720-lyase reactions of P450 17A1. Phosphatidylethanolamine (PE) specifically enhanced 1720-lyase activity up to 2.4-fold in the presence of phosphatidylcholine. On the other hand, P450 11B2-catalyzed production of aldosterone from 11-deoxycorticosterone was very low and from 18-hydroxycorticosterone nil, implying low processivity. DOPC or cardiolipin, which is exclusively located in the inner mitochondrial membrane, maximized aldosterone yield. In sharp contrast, reconstitution of homologous P450 11B1 with DOPC significantly decreased corticosterone formation without affecting the synthesis of 18-hydroxycorticosterone. The intrinsic fluorescence of P450 17A1 and 11B2 increased in the presence of DOPC, DLPC and PE. Acrylamide quenching studies showed that PE decreased solvent accessibility for tryptophan in P450 17A1, as did 20:4 PC or 18:2 PC for P450 11B2. A moderately positive correlation between the proportion of high-spin substrate-bound species and catalytic activity was only observed in the presence of phosphatidylcholines with low-temperature phase transition. These results demonstrate the potential for phospholipids to regulate the activity of steroidogenic P450 activities and thereby steroid hormone biosynthetic pathways.

Ligand Access Channels in Cytochrome P450 Enzymes: A Review

Quantitative structure-activity relationships may bring invaluable information on structural elements of both enzymes and substrates that, together, govern substrate specificity. Buried active sites in cytochrome P450 enzymes are connected to the solvent by a network of channels exiting at the distal surface of the protein. This review presents different in silico tools that were developed to uncover such channels in P450 crystal structures. It also lists some of the experimental evidence that actually suggest that these predicted channels might indeed play a critical role in modulating P450 functions. Amino acid residues at the entrance of the channels may participate to a first global ligand recognition of ligands by P450 enzymes before they reach the buried active site. Moreover, different P450 enzymes show different networks of predicted channels. The plasticity of P450 structures is also important to take into account when looking at how channels might play their role.

Homology Modeling and Molecular Dynamics Simulation Combined with X-ray Solution Scattering Defining Protein Structures of Thromboxane and Prostacyclin Synthases

A combination of molecular dynamics (MD) simulations and X-ray scattering (SAXS) has emerged as the approach of choice for studying protein structures and dynamics in solution. This approach has potential applications for membrane proteins that neither are soluble nor form crystals easily. We explore the water-coupled dynamic structures of thromboxane synthase (TXAS) and prostacyclin synthase (PGIS) from scanning HPLC-SAXS measurements combined with MD ensemble analyses. Both proteins are heme-containing enzymes in the cytochrome P450 family, known as prostaglandin H₂ (PGH₂) isomerase, with counter-functions in regulation of platelet aggregation. Currently, the X-ray crystallographic structures of PGIS are available, but those for TXAS are not. The use of homology modeling of the TXAS structure with ns-μs explicit water solvation MD simulations allows much more accurate estimation of the configuration space with loop motion and origin of the protein behaviors in solution. In contrast to the stability of the conserved PGIS structure in solution, the pronounced TXAS flexibility has been revealed to have unstructured loop regions in connection with the characteristic P450 structural elements. The MD-derived and experimental-solution SAXS results are in excellent agreement. The significant protein internal motions, whole-molecule structures, and potential problems with protein folding, crystallization, and functionality are examined.

Structural insights into the mechanism of the drastic changes in enzymatic activity of the cytochrome P450 vitamin D3 hydroxylase (CYP107BR1) caused by a mutation distant from the active site

Cytochromes P450 (P450s) are haem-containing enzymes that catalyze medically and industrially important oxidative reactions, and many P450s have been subjected to directed evolution and site-directed mutagenesis to improve their activity and substrate specificity. Nonetheless, in most cases the mechanism that leads to drastic changes in specific activity after the introduction of an amino-acid substitution distant from the active-site pocket is unclear. Here, two crystal structures of inactive mutants of the P450 vitamin D3 hydroxylase (Vdh), Vdh-F106V and Vdh-L348M, which were obtained in the course of protein-engineering experiments on Vdh, are reported. The overall structures of these mutants show an open conformation similar to that of wild-type Vdh (Vdh-WT), whereas a rearrangement of the common main-chain hydrogen bonds is observed in the CD-loop (residues 102-106), resulting in a more compactly folded CD-loop relative to that of Vdh-WT. The previously reported structures of Vdh-WT and of the highly active Vdh-T107A and Vdh-K1 mutants have a more stretched CD-loop, with partial formation of 310-helix-type hydrogen bonds, both in the open and closed states. Molecular-dynamics simulations also showed that the frequency of the 310-helix is significantly reduced in Vdh-F106V and Vdh-L348M. The closed conformation is crucial for substrate and ferredoxin binding to initiate the catalytic reaction of Vdh. Therefore, it is implied that the small local structural changes observed in this study might disrupt the conformational transition from the open to the closed state, thereby leading to a complete loss of vitamin D3 hydroxylase activity.

Quasielastic neutron scattering in biology: Theory and applications

Neutrons scatter quasielastically from stochastic, diffusive processes, such as overdamped vibrations, localized diffusion and transitions between energy minima. In biological systems, such as proteins and membranes, these relaxation processes are of considerable physical interest. We review here recent methodological advances and applications of quasielastic neutron scattering (QENS) in biology, concentrating on the role of molecular dynamics simulation in generating data with which neutron profiles can be unambiguously interpreted. We examine the use of massively-parallel computers in calculating scattering functions, and the application of Markov state modeling. The decomposition of MD-derived neutron dynamic susceptibilities is described, and the use of this in combination with NMR spectroscopy. We discuss dynamics at very long times, including approximations to the infinite time mean-square displacement and nonequilibrium aspects of single-protein dynamics. Finally, we examine how neutron scattering and MD can be combined to provide information on lipid nanodomains. This article is part of a Special Issue entitled "Science for Life" Guest Editor: Dr. Austen Angell, Dr. Salvatore Magazù and Dr. Federica Migliardo.

Discovery of a regioselectivity switch in nitrating P450s guided by molecular dynamics simulations and Markov models

The dynamic motions of protein structural elements, particularly flexible loops, are intimately linked with diverse aspects of enzyme catalysis. Engineering of these loop regions can alter protein stability, substrate binding, and even dramatically impact enzyme function. When these flexible regions are structurally unresolvable, computational reconstruction in combination with large-scale molecular dynamics simulations can be used to guide the engineering strategy. Here, we present a collaborative approach consisting of both experiment and computation that led to the discovery of a single mutation in the F/G loop of the nitrating cytochrome P450 TxtE that simultaneously controls loop dynamics and completely shifts the enzyme's regioselectivity from the C4 to the C5 position of L-tryptophan. Furthermore, we find that this loop mutation is naturally present in a subset of homologous nitrating P450s and confirm that these uncharacterized enzymes exclusively produce 5-nitro-L-tryptophan, a previously unknown biosynthetic intermediate.

Cytochrome P450 17A1 Interactions with the FMN Domain of Its Reductase as Characterized by NMR

To accomplish key physiological processes ranging from drug metabolism to steroidogenesis, human microsomal cytochrome P450 enzymes require the sequential input of two electrons delivered by the FMN domain of NADPH-cytochrome P450 reductase. Although some human microsomal P450 enzymes can instead accept the second electron from cytochrome b5, for human steroidogenic CYP17A1, the cytochrome P450 reductase FMN domain delivers both electrons, and b5 is an allosteric modulator. The structural basis of these key but poorly understood protein interactions was probed by solution NMR using the catalytically competent soluble domains of each protein. Formation of the CYP17A1·FMN domain complex induced differential line broadening of the NMR signal for each protein. Alterations in the exchange dynamics generally occurred for residues near the surface of the flavin mononucleotide, including 87-90 (loop 1), and for key CYP17A1 active site residues. These interactions were modulated by the identity of the substrate in the buried CYP17A1 active site and by b5. The FMN domain outcompetes b5 for binding to CYP17A1 in the three-component system. These results and comparison with previous NMR studies of the CYP17A1·b5 complex suggest a model of CYP17A1 enzyme regulation.

F429 Regulation of Tunnels in Cytochrome P450 2B4: A Top Down Study of Multiple Molecular Dynamics Simulations

The root causes of the outcomes of the single-site mutation in enzymes remain by and large not well understood. This is the case of the F429H mutant of the cytochrome P450 (CYP) 2B4 enzyme where the substitution, on the proximal surface of the active site, of a conserved phenylalanine 429 residue with histidine seems to hamper the formation of the active species, Compound I (porphyrin cation radical-Fe^(IV) = O, Cpd I) from the ferric hydroperoxo (Fe^(III)OOH^-, Cpd 0) precursor. Here we report a study based on extensive molecular dynamic (MD) simulations of 4 CYP-2B4 point mutations compared to the WT enzyme, having the goal of better clarifying the importance of the proximal Phe429 residue on CYP 2B4 catalytic properties. To consolidate the huge amount of data coming from five simulations and extract the most distinct structural features of the five species studied we made an extensive use of cluster analysis. The results show that all studied single polymorphisms of F429, with different side chain properties: i) drastically alter the reservoir of conformations accessible by the protein, perturbing global dynamics ii) expose the thiolate group of residue Cys436 to the solvent, altering the electronic properties of Cpd0 and iii) affect the various ingress and egress channels connecting the distal sites with the bulk environment, altering the reversibility of these channels. In particular, it was observed that the wild type enzyme exhibits unique structural features as compared to all mutant species in terms of weak interactions (hydrogen bonds) that generate a completely different dynamical behavior of the complete system. Albeit not conclusive, the current computational investigation sheds some light on the subtle and critical effects that proximal single-site mutations can exert on the functional mechanisms of human microsomal CYPs which should go rather far beyond local structure characterization.

Incorporation of charged residues in the CYP2J2 F-G loop disrupts CYP2J2-lipid bilayer interactions

CYP2J2 epoxygenase is an extrahepatic, membrane bound cytochrome P450 (CYP) that is primarily found in the heart and mediates endogenous fatty acid metabolism. CYP2J2 interacts with membranes through an N-terminal anchor and various non-contiguous hydrophobic residues. The molecular details of the motifs that mediate membrane interactions are complex and not fully understood. To gain better insights of these complex protein-lipid interactions, we employed molecular dynamics (MD) simulations using a highly mobile membrane mimetic (HMMM) model that enabled multiple independent spontaneous membrane binding events to be captured. Simulations revealed that CYP2J2 engages with the membrane at the F-G loop through hydrophobic residues Trp-235, Ille-236, and Phe-239. To explore the role of these residues, three F-G loop mutants were modeled from the truncated CYP2J2 construct (Δ34) which included Δ34-I236D, Δ34-F239H and Δ34-I236D/F239H. Using the HMMM coordinates of CYP2J2, the simulations were extended to a full POPC membrane which showed a significant decrease in the depth of insertion for each of the F-G loop mutants. The CYP2J2 F-G loop mutants were expressed in E. coli and were shown to be localized to the cytosolic fraction at a greater percentage relative to construct Δ34. Notably, the functional data demonstrated that the double mutant, Δ34-I236D/F239H, maintained native-like enzymatic activity. The membrane insertion characteristics were examined by monitoring CYP2J2 Trp-quenching fluorescence spectroscopy upon binding nanodiscs containing pyrene phospholipids. Relative to the Δ34 construct, the F-G loop mutants exhibited lower Trp quenching and membrane insertion. Taken together, the results suggest that the mutants exhibit a different membrane topology in agreement with the MD simulations and provide important evidence towards the involvement of key residues in the F-G loop of CYP2J2.

Comparison of intrinsic dynamics of cytochrome p450 proteins using normal mode analysis

Cytochrome P450 enzymes are hemeproteins that catalyze the monooxygenation of a wide-range of structurally diverse substrates of endogenous and exogenous origin. These heme monooxygenases receive electrons from NADH/NADPH via electron transfer proteins. The cytochrome P450 enzymes, which constitute a diverse superfamily of more than 8,700 proteins, share a common tertiary fold but < 25% sequence identity. Based on their electron transfer protein partner, cytochrome P450 proteins are classified into six broad classes. Traditional methods of pro are based on the canonical paradigm that attributes proteins' function to their three-dimensional structure, which is determined by their primary structure that is the amino acid sequence. It is increasingly recognized that protein dynamics play an important role in molecular recognition and catalytic activity. As the mobility of a protein is an intrinsic property that is encrypted in its primary structure, we examined if different classes of cytochrome P450 enzymes display any unique patterns of intrinsic mobility. Normal mode analysis was performed to characterize the intrinsic dynamics of five classes of cytochrome P450 proteins. The present study revealed that cytochrome P450 enzymes share a strong dynamic similarity (root mean squared inner product > 55% and Bhattacharyya coefficient > 80%), despite the low sequence identity (< 25%) and sequence similarity (< 50%) across the cytochrome P450 superfamily. Noticeable differences in Cα atom fluctuations of structural elements responsible for substrate binding were noticed. These differences in residue fluctuations might be crucial for substrate selectivity in these enzymes.

Structure, dynamics, and function of the monooxygenase P450 BM-3: insights from computer simulations studies

The monooxygenase P450 BM-3 is a NADPH-dependent fatty acid hydroxylase enzyme isolated from soil bacterium Bacillus megaterium. As a pivotal member of cytochrome P450 superfamily, it has been intensely studied for the comprehension of structure-dynamics-function relationships in this class of enzymes. In addition, due to its peculiar properties, it is also a promising enzyme for biochemical and biomedical applications. However, despite the efforts, the full understanding of the enzyme structure and dynamics is not yet achieved. Computational studies, particularly molecular dynamics (MD) simulations, have importantly contributed to this endeavor by providing new insights at an atomic level regarding the correlations between structure, dynamics, and function of the protein. This topical review summarizes computational studies based on MD simulations of the cytochrome P450 BM-3 and gives an outlook on future directions.

Structural, functional, and spectroscopic characterization of the substrate scope of the novel nitrating cytochrome P450 TxtE

A novel cytochrome P450 enzyme, TxtE, was recently shown to catalyze the direct aromatic nitration of L-tryptophan. This unique chemistry inspired us to ask whether TxtE could serve as a platform for engineering new nitration biocatalysts to replace current harsh synthetic methods. As a first step toward this goal, and to better understand the wild-type enzyme, we obtained high-resolution structures of TxtE in its substrate-free and substrate-bound forms. We also screened a library of substrate analogues for spectroscopic indicators of binding and for production of nitrated products. From these results, we found that the wild-type enzyme accepts moderate decoration of the indole ring, but the amino acid moiety is crucial for binding and correct positioning of the substrate and therefore less amenable to modification. A nitrogen atom is essential for catalysis, and a carbonyl must be present to recruit the αB'1 helix of the protein to seal the binding pocket.

Characterization of the active site properties of CYP4F12

Cytochrome P450 4F12 is a drug-metabolizing enzyme that is primarily expressed in the liver, kidney, colon, small intestine, and heart. The properties of CYP4F12 that may impart an increased catalytic selectivity (decreased promiscuity) were explored through in vitro metabolite elucidation, kinetic isotope effect experiments, and computational modeling of the CYP4F12 active site. By using astemizole as a probe substrate for CYP4F12 and CYP3A4, it was observed that although CYP4F12 favored astemizole O-demethylation as the primary route of metabolism, CYP3A4 was capable of metabolizing astemizole at multiple sites on the molecule. Deuteration of astemizole at the site of O-demethylation resulted in an isotope effect of 7.1 as well as an 8.3-fold decrease in the rate of clearance for astemizole by CYP4F12. Conversely, although an isotope effect of 3.8 was observed for the formation of the O-desmethyl metabolite when deuterated astemizole was metabolized by CYP3A4, there was no decrease in the clearance of astemizole. Development of a homology model of CYP4F12 based on the crystal structure of cytochrome P450 BM3 predicted an active site volume for CYP4F12 that was approximately 76% of the active site volume of CYP3A4. As predicted, multiple favorable binding orientations were available for astemizole docked into the active site of CYP3A4, but only a single binding orientation with the site of O-demethylation oriented toward the heme was identified for CYP4F12. Overall, it appears that although CYP4F12 may be capable of binding similar ligands to other cytochrome P450 enzymes such as CYP3A4, the ability to achieve catalytically favorable orientations may be inherently more difficult because of the increased steric constraints of the CYP4F12 active site.

Dynamics and flexibility of human aromatase probed by FTIR and time resolved fluorescence spectroscopy

Human aromatase (CYP19A1) is a steroidogenic cytochrome P450 converting androgens into estrogens. No ligand-free crystal structure of the enzyme is available to date. The crystal structure in complex with the substrate androstenedione and the steroidal inhibitor exemestane shows a very compact conformation of the enzyme, leaving unanswered questions on the conformational changes that must occur to allow access of the ligand to the active site. As H/D exchange kinetics followed by FTIR spectroscopy can provide information on the conformational changes in proteins where solvent accessibility is affected, here the amide I region was used to measure the exchange rates of the different elements of the secondary structure for aromatase in the ligand-free form and in the presence of the substrate androstenedione and the inhibitor anastrozole. Biphasic exponential functions were found to fit the H/D exchange data collected as a function of time. Two exchange rates were assigned to two populations of protons present in different flexible regions of the protein. The addition of the substrate androstenedione and the inhibitor anastrozole lowers the H/D exchange rates of the α-helices of the enzyme when compared to the ligand-free form. Furthermore, the presence of the inhibitor anastrozole lowers exchange rate constant (k₁) for β-sheets from 0.22±0.06 min⁻¹ for the inhibitor-bound enzyme to 0.12±0.02 min⁻¹ for the free protein. Dynamics effects localised in helix F were studied by time resolved fluorescence. The data demonstrate that the fluorescence lifetime component associated to Trp224 emission undergoes a shift toward longer lifetimes (from ≈5.0 to ≈5.5 ns) when the substrate or the inhibitor are present, suggesting slower dynamics in the presence of ligands. Together the results are consistent with different degrees of flexibility of the access channel and therefore different conformations adopted by the enzyme in the free, substrate- and inhibitor-bound forms.

Coupled flexibility change in cytochrome P450cam substrate binding determined by neutron scattering, NMR, and molecular dynamics simulation

Neutron scattering and nuclear magnetic resonance relaxation experiments are combined with molecular dynamics (MD) simulations in a novel, to our knowledge, approach to investigate the change in internal dynamics on substrate (camphor) binding to a protein (cytochrome P450cam). The MD simulations agree well with both the neutron scattering, which furnishes information on global flexibility, and the nuclear magnetic resonance data, which provides residue-specific order parameters. Decreased fluctuations are seen in the camphor-bound form using all three techniques, dominated by changes in specific regions of the protein. The combined experimental and simulation results permit a detailed description of the dynamical change, which involves modifications in the coupling between the dominant regions and concomitant substrate access channel closing, via specific salt-bridge, hydrogen-bonding, and hydrophobic interactions. The work demonstrates how the combination of complementary experimental spectroscopies with MD simulation can provide an in-depth description of functional dynamical protein changes.

Diversity of P450 enzymes in the biosynthesis of natural products

Diverse oxygenation patterns of natural products generated by secondary metabolic pathways in microorganisms and plants are largely achieved through the tailoring reactions catalysed by cytochrome P450 enzymes (P450s). P450s are a large family of oxidative hemoproteins found in all life forms from prokaryotes to humans. Understanding the reactivity and selectivity of these fascinating C-H bond-activating catalysts will advance their use in generating valuable pharmaceuticals and products for medicine, agriculture and industry. A major strength of this P450 group is its set of established enzyme-substrate relationships, the source of the most detailed knowledge on how P450 enzymes work. Engineering microbial-derived P450 enzymes to accommodate alternative substrates and add new functions continues to be an important near- and long-term practical goal driving the structural characterization of these molecules. Understanding the natural evolution of P450 structure-function should accelerate metabolic engineering and directed evolutionary approaches to enhance diversification of natural product structures and other biosynthetic applications.

CYP51 structures and structure-based development of novel, pathogen-specific inhibitory scaffolds

CYP51 (sterol 14α-demethylase) is a cytochrome P450 enzyme essential for sterol biosynthesis and the primary target for clinical and agricultural antifungal azoles. The azoles that are currently in clinical use for systemic fungal infections represent modifications of two basic scaffolds, ketoconazole and fluconazole, all of them being selected based on their antiparasitic activity in cellular experiments. By studying direct inhibition of CYP51 activity across phylogeny including human pathogens Trypanosoma brucei, Trypanosoma cruzi and Leishmania infantum, we identified three novel protozoa-specific inhibitory scaffolds, their inhibitory potency correlating well with antiprotozoan activity. VNI scaffold (carboxamide containing β-phenyl-imidazoles) is the most promising among them: killing T. cruzi amastigotes at low nanomolar concentration, it is also easy to synthesize and nontoxic. Oral administration of VNI (up to 400 mg/kg) neither leads to mortality nor reveals significant side effects up to 48 h post treatment using an experimental mouse model of acute toxicity. Trypanosomatidae CYP51 crystal structures determined in the ligand-free state and complexed with several azole inhibitors as well as a substrate analog revealed high rigidity of the CYP51 substrate binding cavity, which must be essential for the enzyme strict substrate specificity and functional conservation. Explaining profound potency of the VNI inhibitory scaffold, the structures also outline guidelines for its further development. First steps of the VNI scaffold optimization have been undertaken; the results presented here support the notion that CYP51 structure-based rational design of more efficient, pathogen-specific inhibitors represents a highly promising direction.

Plasticity of CYP2B enzymes: structural and solution biophysical methods

In the past three years, major advances in understanding cytochrome P450 2B (CYP2B) structure-function relationships have been made through determination of multiple ligand-bound and one ligand-free X-ray crystal structure of CYP2B4 and one ligand-bound X-ray crystal structure of CYP2B6. These structures have provided insight into the features that provide the high degree of plasticity of the enzymes. A combination of a phenylalanine cluster that allows for concerted movement of helices F through G and a conserved set of electrostatic interactions involving Arg(262) facilitates movement of this region to accommodate binding of ligands of various sizes without perturbing most of the P450 fold. Integrating solution based techniques such as NMR or deuterium exchange mass spectrometry (DXMS) with computational methods including molecular docking has provided further insight into enzyme behavior upon ligand binding. In addition, extended molecular dynamics simulations have provided a link between an open and a closed conformation of ligand-free CYP2B4 found in crystal structures. Other studies revealed the utility of rational engineering in improving stability of P450s to facilitate structural studies. The solution and computational results combined with the X-ray crystal structures yield a comprehensive picture of how these enzymes adopt different conformations to bind various ligands.

Investigation of indazole unbinding pathways in CYP2E1 by molecular dynamics simulations

Human microsomal cytochrome P450 2E1 (CYP2E1) can oxidize not only low molecular weight xenobiotic compounds such as ethanol, but also many endogenous fatty acids. The crystal structure of CYP2E1 in complex with indazole reveals that the active site is deeply buried into the protein center. Thus, the unbinding pathways and associated unbinding mechanisms remain elusive. In this study, random acceleration molecular dynamics simulations combined with steered molecular dynamics and potential of mean force calculations were performed to identify the possible unbinding pathways in CYP2E1. The results show that channel 2c and 2a are most likely the unbinding channels of CYP2E1. The former channel is located between helices G and I and the B-C loop, and the latter resides between the region formed by the F-G loop, the B-C loop and the β1 sheet. Phe298 and Phe478 act as the gate keeper during indazole unbinding along channel 2c and 2a, respectively. Previous site-directed mutagenesis experiments also supported these findings.

Motion and flexibility in human cytochrome p450 aromatase

The crystal structures of human placental aromatase in complex with the substrate androstenedione and exemestane have revealed an androgen-specific active site and the structural basis for higher order organization. However, X-ray structures do not provide accounts of movements due to short-range fluctuations, ligand binding and protein-protein association. In this work, we conduct normal mode analysis (NMA) revealing the intrinsic fluctuations of aromatase, deduce the internal modes in membrane-free and membrane-integrated monomers as well as the intermolecular modes in oligomers, and propose a quaternary organization for the endoplasmic reticulum (ER) membrane integration. Dynamics of the crystallographic oligomers from NMA is found to be in agreement with the isotropic thermal factors from the X-ray analysis. Calculations of the root mean square fluctuations of the C-alpha atoms from their equilibrium positions confirm that the rigid-core structure of aromatase is intrinsic regardless of the changes in steroid binding interactions, and that aromatase self-association does not deteriorate the rigidity of the catalytic cleft. Furthermore, NMA on membrane-integrated aromatase shows that the internal modes in all likelihood contribute to breathing of the active site access channel. The collective intermolecular hinge bending and twisting modes provide the flexibility in the quaternary association necessary for membrane integration of the aromatase oligomers. Taken together, fluctuations of the active site, the access channel, and the heme-proximal cavity, and a dynamic quaternary organization could all be essential components of the functional aromatase in its role as an ER membrane-embedded steroidogenic enzyme.

Tailoring reactions catalyzed by heme-dependent enzymes: spectroscopic characterization of the L-tryptophan-nitrating cytochrome P450 TxtE

There is a truly vast quantity of research articles and textbooks, aimed at a variety of audiences, on cytochromes P450. However, a large amount of specialized terminology has become associated with these enzymes, which can be daunting to those new to the field. The aim of this chapter is to give a brief overview of the functions and importance of cytochromes P450 with particular emphasis on their roles as tailoring enzymes in natural product biosynthetic pathways. Differences between the biosynthetic enzymes and their catabolic counterparts are highlighted. Assays used to investigate substrate binding to cytochromes P450 are described using TxtE, a recently discovered unique nitrating enzyme involved in thaxtomin A biosynthesis, as an example.

Active-site hydration and water diffusion in cytochrome P450cam: a highly dynamic process

Long-timescale molecular dynamics simulations (300 ns) are performed on both the apo- (i.e., camphor-free) and camphor-bound cytochrome P450cam (CYP101). Water diffusion into and out of the protein active site is observed without biased sampling methods. During the course of the molecular dynamics simulation, an average of 6.4 water molecules is observed in the camphor-binding site of the apo form, compared to zero water molecules in the binding site of the substrate-bound form, in agreement with the number of water molecules observed in crystal structures of the same species. However, as many as 12 water molecules can be present at a given time in the camphor-binding region of the active site in the case of apo-P450cam, revealing a highly dynamic process for hydration of the protein active site, with water molecules exchanging rapidly with the bulk solvent. Water molecules are also found to exchange locations frequently inside the active site, preferentially clustering in regions surrounding the water molecules observed in the crystal structure. Potential-of-mean-force calculations identify thermodynamically favored trans-protein pathways for the diffusion of water molecules between the protein active site and the bulk solvent. Binding of camphor in the active site modifies the free-energy landscape of P450cam channels toward favoring the diffusion of water molecules out of the protein active site.

Substrate preferences and catalytic parameters determined by structural characteristics of sterol 14alpha-demethylase (CYP51) from Leishmania infantum

Leishmaniasis is a major health problem that affects populations of ∼90 countries worldwide, with no vaccine and only a few moderately effective drugs. Here we report the structure/function characterization of sterol 14α-demethylase (CYP51) from Leishmania infantum. The enzyme catalyzes removal of the 14α-methyl group from sterol precursors. The reaction is essential for membrane biogenesis and therefore has great potential to become a target for antileishmanial chemotherapy. Although L. infantum CYP51 prefers C4-monomethylated sterol substrates such as C4-norlanosterol and obtusifoliol (V(max) of ∼10 and 8 min(-1), respectively), it is also found to 14α-demethylate C4-dimethylated lanosterol (V(max) = 0.9 min(-1)) and C4-desmethylated 14α-methylzymosterol (V(max) = 1.9 min(-1)). Binding parameters with six sterols were tested, with K(d) values ranging from 0.25 to 1.4 μM. Thus, L. infantum CYP51 is the first example of a plant-like sterol 14α-demethylase, where requirements toward the composition of the C4 atom substituents are not strict, indicative of possible branching in the postsqualene portion of sterol biosynthesis in the parasite. Comparative analysis of three CYP51 substrate binding cavities (Trypanosoma brucei, Trypanosoma cruzi, and L. infantum) suggests that substrate preferences of plant- and fungal-like protozoan CYP51s largely depend on the differences in the enzyme active site topology. These minor structural differences are also likely to underlie CYP51 catalytic rates and drug susceptibility and can be used to design potent and specific inhibitors.

Homology modelling of Drosophila cytochrome P450 enzymes associated with insecticide resistance

BACKGROUND

Overexpression of the cytochrome P450 gene Cyp6g1 confers resistance against DDT and a broad range of other insecticides in Drosophila melanogaster Meig. In the absence of crystal structures of CYP6G1 or complexes with its substrates, structural studies rely on homology modelling and ligand docking to understand P450-substrate interactions.

RESULTS

Homology models are presented for CYP6G1, a P450 associated with resistance to DDT and neonicotinoids, and two other enzymes associated with insecticide resistance in D. melanogaster, CYP12D1 and CYP6A2. The models are based on a template of the X-ray structure of the phylogenetically related human CYP3A4, which is known for its broad substrate specificity. The model of CYP6G1 has a much smaller active site cavity than the template. The cavity is also 'V'-shaped and is lined with hydrophobic residues, showing high shape and chemical complementarity with the molecular characteristics of DDT. Comparison of the DDT-CYP6G1 complex and a non-resistant CYP6A2 homology model implies that tight-fit recognition of this insecticide is important in CYP6G1. The active site can accommodate differently shaped substrates ranging from imidacloprid to malathion but not the pyrethroids permethrin and cyfluthrin.

CONCLUSION

The CYP6G1, CYP12D1 and CYP6A2 homology models can provide a structural insight into insecticide resistance in flies overexpressing P450 enzymes with broad substrate specificities.

Structural evidence for enhancement of sequential vitamin D3 hydroxylation activities by directed evolution of cytochrome P450 vitamin D3 hydroxylase

Vitamin D(3) hydroxylase (Vdh) isolated from actinomycete Pseudonocardia autotrophica is a cytochrome P450 (CYP) responsible for the biocatalytic conversion of vitamin D(3) (VD(3)) to 1α,25-dihydroxyvitamin D(3) (1α,25(OH)(2)VD(3)) by P. autotrophica. Although its biological function is unclear, Vdh is capable of catalyzing the two-step hydroxylation of VD(3), i.e. the conversion of VD(3) to 25-hydroxyvitamin D(3) (25(OH)VD(3)) and then of 25(OH)VD(3) to 1α,25(OH)(2)VD(3), a hormonal form of VD(3). Here we describe the crystal structures of wild-type Vdh (Vdh-WT) in the substrate-free form and of the highly active quadruple mutant (Vdh-K1) generated by directed evolution in the substrate-free, VD(3)-bound, and 25(OH)VD(3)-bound forms. Vdh-WT exhibits an open conformation with the distal heme pocket exposed to the solvent both in the presence and absence of a substrate, whereas Vdh-K1 exhibits a closed conformation in both the substrate-free and substrate-bound forms. The results suggest that the conformational equilibrium was largely shifted toward the closed conformation by four amino acid substitutions scattered throughout the molecule. The substrate-bound structure of Vdh-K1 accommodates both VD(3) and 25(OH)VD(3) but in an anti-parallel orientation. The occurrence of the two secosteroid binding modes accounts for the regioselective sequential VD(3) hydroxylation activities. Moreover, these structures determined before and after directed evolution, together with biochemical and spectroscopic data, provide insights into how directed evolution has worked for significant enhancement of both the VD(3) 25-hydroxylase and 25(OH)VD(3) 1α-hydroxylase activities.

How does the reductase help to regulate the catalytic cycle of cytochrome P450 3A4 using the conserved water channel?

Water molecules play a major role in the P450 catalytic cycle. Here, we locate the preferred water pathways and their gating mechanisms for the human cytochrome P450 3A4 (CYP3A4) and elucidate the role of the cytochrome P450 reductase (CPR) in turning on and activating these water channels. We perform explicit solvent molecular dynamic simulations of CYP3A4, unbound and bound to two substrates, and with and without the flavin mononucleotide (FMN)-binding domain of CPR. We observe in/out passage of water molecules via a water-specific and conserved channel (aqueduct) located between the active site and the heme proximal side. We find that the aqueduct gating mechanism is mediated by R375, the conserved arginine that salt bridges with the heme 7-propionate. When R375 rotates, it opens the aqueduct and establishes a connection between a cluster of active site water molecules network and the bulk solvent. The aqueduct region overlaps with the CPR binding-site to CYP3A4. Indeed, we find that when the FMN domain of CPR binds to CYP3A4, the aqueduct fully opens up, thereby allowing a flow of water molecules. The aqueduct's opening can permit proton transfer, shuttling the protons to the active site through ordered water molecules. In addition, the expulsion of water molecules via the aqueduct contributes to substrate binding. As such, the CPR binding has a function: it triggers the aqueduct's opening and thereby enables a proton shuttle pathway, which is needed for the dioxygen activation. This mechanism could be a general paradigm in P450s.

Regio- and stereospecificity of filipin hydroxylation sites revealed by crystal structures of cytochrome P450 105P1 and 105D6 from Streptomyces avermitilis

The polyene macrolide antibiotic filipin is widely used as a probe for cholesterol and a diagnostic tool for type C Niemann-Pick disease. Two position-specific P450 enzymes are involved in the post-polyketide modification of filipin during its biosynthesis, thereby providing molecular diversity to the "filipin complex." CYP105P1 and CYP105D6 from Streptomyces avermitilis, despite their high sequence similarities, catalyze filipin hydroxylation at different positions, C26 and C1', respectively. Here, we determined the crystal structure of the CYP105P1-filipin I complex. The distal pocket of CYP105P1 has the second largest size among P450 hydroxylases that act on macrolide substrates. Compared with previously determined substrate-free structures, the FG helices showed significant closing motion on substrate binding. The long BC loop region adopts a unique extended conformation without a B' helix. The binding site is essentially hydrophobic, but numerous water molecules are involved in recognizing the polyol side of the substrate. Therefore, the distal pocket of CYP105P1 provides a specific environment for the large filipin substrate to bind with its pro-S side of position C26 directed toward the heme iron. The ligand-free CYP105D6 structure was also determined. A small sub-pocket accommodating the long alkyl side chain of filipin I was observed in the CYP105P1 structure but was absent in the CYP105D6 structure, indicating that filipin cannot bind to CYP105D6 with a similar orientation due to steric hindrance. This observation can explain the strict regiospecificity of these enzymes.

Theoretical characterization of substrate access/exit channels in the human cytochrome P450 3A4 enzyme: involvement of phenylalanine residues in the gating mechanism

The human cytochrome P450 3A4 mono-oxygenates ∼50% of all drugs. Its substrates/products enter/leave the active site by access/exit channels. Here, we perform steered molecular dynamics simulations, pulling the products temazepam and testosterone-6βOH out of the P450 3A4 enzyme in order to identify the preferred substrate/product pathways and their gating mechanism. We locate six different egress pathways of products from the active site with different exit preferences for the two products and find that there is more than just one access/exit channel in CYP3A4. The so-called solvent channel manifests the largest opening for both tested products, thereby identifying this channel as a putative substrate channel. Most channels consist of one or two π-stacked phenylalanine residues that serve as gate keepers. The oxidized drug breaks the hydrophobic interactions of the gating residues and forms mainly hydrophobic contacts with the gate. We argue that product exit preferences in P450s are regulated by protein−substrate specificity.