Ligand Access Channels in Cytochrome P450 Enzymes: A Review

Reinforcing Tunnel Network Exploration in Proteins Using Gaussian Accelerated Molecular Dynamics

Tunnels are structural conduits in biomolecules responsible for transporting chemical compounds and solvent molecules from the active site. They have been shown to be present in a wide variety of enzymes across all functional and structural classes. However, the study of such pathways is experimentally challenging, because they are typically transient. Computational methods, such as molecular dynamics (MD) simulations, have been successfully proposed to explore tunnels. Conventional MD (cMD) provides structural details to characterize tunnels but suffers from sampling limitations to capture rare tunnel openings on longer time scales. Therefore, in this study, we explored the potential of Gaussian accelerated MD (GaMD) simulations to improve the exploration of complex tunnel networks in enzymes. We used the haloalkane dehalogenase LinB and its two variants with engineered transport pathways, which are not only well-known for their application potential but have also been extensively studied experimentally and computationally regarding their tunnel networks and their importance in multistep catalytic reactions. Our study demonstrates that GaMD efficiently improves tunnel sampling and allows the identification of all known tunnels for LinB and its two mutants. Furthermore, the improved sampling provided insight into a previously unknown transient side tunnel (ST). The extensive conformational landscape explored by GaMD simulations allowed us to investigate in detail the mechanism of ST opening. We determined variant-specific dynamic properties of ST opening, which were previously inaccessible due to limited sampling of cMD. Our comprehensive analysis supports multiple indicators of the functional relevance of the ST, emphasizing its potential significance beyond structural considerations. In conclusion, our research proves that the GaMD method can overcome the sampling limitations of cMD for the effective study of tunnels in enzymes, providing further means for identifying rare tunnels in enzymes with the potential for drug development, precision medicine, and rational protein engineering.

Exploring the selectivity of cytochrome P450 for enhanced novel anticancer agent synthesis

Cytochrome P450 monooxygenases are an extensive and unique class of enzymes, which can regio- and stereo-selectively functionalise hydrocarbons by way of oxidation reactions. These enzymes are naturally occurring but have also been extensively applied in a synthesis context, where they are used as efficient biocatalysts. Recently, a biosynthetic pathway where a cytochrome P450 monooxygenase catalyses a critical step of the pathway was uncovered, leading to the production of a number of products that display high antitumour potency. In this work, we use computational techniques to gain insight into the factors that determine the relative yields of the different products. We use conformational search algorithms to understand the substrate stereochemistry. On a machine-learned 3D protein structure, we use molecular docking to obtain a library of favourable poses for substrate-protein interaction. With molecular dynamics, we investigate the most favourable poses for reactivity on a molecular level, allowing us to investigate which protein-substrate interactions favour a given product and thus gain insight into the product selectivity.

Designing cytochrome P450 enzymes for use in cancer gene therapy

Cancer is a significant global socioeconomic burden, as millions of new cases and deaths occur annually. In 2020, almost 10 million cancer deaths were recorded worldwide. Advancements in cancer gene therapy have revolutionized the landscape of cancer treatment. An approach with promising potential for cancer gene therapy is introducing genes to cancer cells that encode for chemotherapy prodrug metabolizing enzymes, such as Cytochrome P450 (CYP) enzymes, which can contribute to the effective elimination of cancer cells. This can be achieved through gene-directed enzyme prodrug therapy (GDEPT). CYP enzymes can be genetically engineered to improve anticancer prodrug conversion to its active metabolites and to minimize chemotherapy side effects by reducing the prodrug dosage. Rational design, directed evolution, and phylogenetic methods are some approaches to developing tailored CYP enzymes for cancer therapy. Here, we provide a compilation of genetic modifications performed on CYP enzymes aiming to build highly efficient therapeutic genes capable of bio-activating different chemotherapeutic prodrugs. Additionally, this review summarizes promising preclinical and clinical trials highlighting engineered CYP enzymes' potential in GDEPT. Finally, the challenges, limitations, and future directions of using CYP enzymes for GDEPT in cancer gene therapy are discussed.

Structure-Function Analysis of the Essential Mycobacterium tuberculosis P450 Drug Target, CYP121A1

Cytochrome P450 CYP121A1 is a well-known drug target against Mycobacterium tuberculosis, the human pathogen that causes the deadly disease tuberculosis (TB). CYP121A1 is a unique P450 enzyme because it uses classical and non-classical P450 catalytic processes and has distinct structural features among P450s. However, a detailed investigation of CYP121A1 protein structures in terms of active site cavity dynamics and key amino acids interacting with bound ligands has yet to be undertaken. To address this research knowledge gap, 53 CYP121A1 crystal structures were investigated in this study. Critical amino acids required for CYP121A1's overall activity were identified and highlighted this enzyme's rigid architecture and substrate selectivity. The CYP121A1-fluconazole crystal structure revealed a novel azole drug-P450 binding mode in which azole heme coordination was facilitated by a water molecule. Fragment-based inhibitor approaches revealed that CYP121A1 can be inhibited by molecules that block the substrate channel or by directly interacting with the P450 heme. This study serves as a reference for the precise understanding of CYP121A1 interactions with different ligands and the structure-function analysis of P450 enzymes in general. Our findings provide critical information for the synthesis of more specific CYP121A1 inhibitors and their development as novel anti-TB drugs.

The Multienzyme Complex Nature of Dehydroepiandrosterone Sulfate Biosynthesis

Dehydroepiandrosterone (DHEA), a precursor of steroid sex hormones, is synthesized by steroid 17-alpha-hydroxylase/17,20-lyase (CYP17A1) with the participation of microsomal cytochrome b5 (CYB5A) and cytochrome P450 reductase (CPR), followed by sulfation by two cytosolic sulfotransferases, SULT1E1 and SULT2A1, for storage and transport to tissues in which its synthesis is not available. The involvement of CYP17A1 and SULTs in these successive reactions led us to consider the possible interaction of SULTs with DHEA-producing CYP17A1 and its redox partners. Text mining analysis, protein-protein network analysis, and gene co-expression analysis were performed to determine the relationships between SULTs and microsomal CYP isoforms. For the first time, using surface plasmon resonance, we detected interactions between CYP17A1 and SULT2A1 or SULT1E1. SULTs also interacted with CYB5A and CPR. The interaction parameters of SULT2A1/CYP17A1 and SULT2A1/CYB5A complexes seemed to be modulated by 3'-phosphoadenosine-5'-phosphosulfate (PAPS). Affinity purification, combined with mass spectrometry (AP-MS), allowed us to identify a spectrum of SULT1E1 potential protein partners, including CYB5A. We showed that the enzymatic activity of SULTs increased in the presence of only CYP17A1 or CYP17A1 and CYB5A mixture. The structures of CYP17A1/SULT1E1 and CYB5A/SULT1E1 complexes were predicted. Our data provide novel fundamental information about the organization of microsomal CYP-dependent macromolecular complexes.

Water Exchange from the Buried Binding Sites of Cytochrome P450 Enzymes 1A2, 2D6, and 3A4 Correlates with Conformational Fluctuations

Human cytochrome P450 enzymes (CYPs) are critical for the metabolism of small-molecule pharmaceuticals (drugs). As such, the prediction of drug metabolism by and drug inhibition of CYP activity is an important component of the drug discovery and design process. Relative to the availability of a wide range of experimental atomic-resolution CYP structures, the development of structure-based CYP activity models has been limited. To better characterize the role of CYP conformational fluctuations in CYP activity, we perform multiple microsecond-scale all-atom explicit-solvent molecular dynamics (MD) simulations on three CYP isoforms, 1A2, 2D6, and 3A4, which together account for the majority of CYP-mediated drug metabolism. The MD simulations employ a variety of positional restraints, ranging from keeping all CYP atoms close to their experimentally determined coordinates to allowing full flexibility. We find that, with full flexibility, large fluctuations in the CYP binding sites correlate with efficient water exchange from these buried binding sites. This is especially true for 1A2, which, when restrained to its crystallographic conformation, is unable to exchange water between the binding site and bulk solvent. These findings imply that, in addition to crystal structures, a representative ensemble of conformational states ought to be included when developing structure-based CYP activity models.

A cytochrome P450 insecticide detoxification mechanism is not conserved across the Megachilidae family of bees

Recent work has demonstrated that many bee species have specific cytochrome P450 enzymes (P450s) that can efficiently detoxify certain insecticides. The presence of these P450s, belonging or closely related to the CYP9Q subfamily (CYP9Q-related), is generally well conserved across the diversity of bees. However, the alfalfa leafcutter bee, Megachile rotundata, lacks CYP9Q-related P450s and is 170-2500 times more sensitive to certain insecticides than bee pollinators with these P450s. The extent to which these findings apply to other Megachilidae bee species remains uncertain. To address this knowledge gap, we sequenced the transcriptomes of four Megachile species and leveraged the data obtained, in combination with publicly available genomic data, to investigate the evolution and function of P450s in the Megachilidae. Our analyses reveal that several Megachilidae species, belonging to the Lithurgini, Megachilini and Anthidini tribes, including all species of the Megachile genus investigated, lack CYP9Q-related genes. In place of these genes Megachile species have evolved phylogenetically distinct CYP9 genes, the CYP9DM lineage. Functional expression of these P450s from M. rotundata reveal they lack the capacity to metabolize the neonicotinoid insecticides thiacloprid and imidacloprid. In contrast, species from the Osmiini and Dioxyini tribes of Megachilidae have CYP9Q-related P450s belonging to the CYP9BU subfamily that are able to detoxify thiacloprid. These findings provide new insight into the evolution of P450s that act as key determinants of insecticide sensitivity in bees and have important applied implications for pesticide risk assessment.

Structural Characterization and Molecular Dynamics Study of the REPI Fusion Protein from Papaver somniferum L

REPI is a pivotal point enzyme in plant benzylisoquinoline alkaloid metabolism as it promotes the evolution of the biosynthetic branch of morphinan alkaloids. Experimental studies of its activity led to the identification of two modules (DRS and DRR) that catalyze two sequential steps of the epimerization of (S)- to (R)-reticuline. Recently, special attention has been paid to its genetic characterization and evolutionary history, but no structural analyses of the REPI protein have been conducted to date. We present here a computational structural characterization of REPI with heme and NADP cofactors in the apo state and in three complexes with substrate (S)-reticuline in DRS and intermediate 1,2-dehydroreticuline in DRS and in DRR. Since no experimental structure exists for REPI, we used its AlphaFold model as a scaffold to build up these four systems, which were submitted to all-atom molecular dynamics (MD) simulations. A comparison of MD results for the four systems revealed key dynamic changes associated with cofactor and ligand binding and provided a dynamic picture of the evolution of their structures and interactions. We also explored the possible dynamic occurrence of tunnels and electrostatic highways potentially involved in alternative mechanisms for channeling the intermediate from DRS to DRR.

Computationally guided bioengineering of the active site, substrate access pathway, and water channels of thermostable cytochrome P450, CYP175A1, for catalyzing the alkane hydroxylation reaction

Understanding structure-function relationships in proteins is pivotal in their development as industrial biocatalysts. In this regard, rational engineering of protein active site access pathways and various tunnels and channels plays a central role in designing competent enzymes with high stability and enhanced efficiency. Here, we report the rational evolution of a thermostable cytochrome P450, CYP175A1, to catalyze the C-H activation reaction of longer-chain alkanes. A strategy combining computational tools with experiments has shown that the substrate scope and enzymatic activity can be enhanced by rational engineering of certain important channels such as the substrate entry and water channels along with the active site of the enzyme. The evolved enzymes showed an improved catalytic rate for hexadecane hydroxylation with high regioselectivity. The Q67L/Y68F mutation showed binding of the substrate in the active site, water channel mutation L80F/V220T showed improved catalytic activity through the peroxide shunt pathway and substrate entry channel mutation W269F/I270A showed better substrate accessibility to the active pocket. All-atom MD simulations provided the rationale for the inactivity of the wild-type CYP175A1 for hexadecane hydroxylation and predicted the above hot-spot residues to enhance the activity. The reaction mechanism was studied by QM/MM calculations for enzyme-substrate complexes and reaction intermediates. Detailed thermal and thermodynamic stability of all the mutants were analyzed and the results showed that the evolved enzymes were thermally stable. The present strategy showed promising results, and insights gained from this work can be applied to the general enzymatic system to expand substrate scope and improve catalytic activity.

Structure-Function Analysis of the Biotechnologically Important Cytochrome P450 107 (CYP107) Enzyme Family

Cytochrome P450 monooxygenases (CYPs; P450s) are a superfamily of heme-containing enzymes that are recognized for their vast substrate range and oxidative multifunctionality. CYP107 family members perform hydroxylation and epoxidation processes, producing a variety of biotechnologically useful secondary metabolites. Despite their biotechnological importance, a thorough examination of CYP107 protein structures regarding active site cavity dynamics and key amino acids interacting with bound ligands has yet to be undertaken. To address this research knowledge gap, 44 CYP107 crystal structures were investigated in this study. We demonstrate that the CYP107 active site cavity is very flexible, with ligand binding reducing the volume of the active site in some situations and increasing volume size in other instances. Polar interactions between the substrate and active site residues result in crucial salt bridges and the formation of proton shuttling pathways. Hydrophobic interactions, however, anchor the substrate within the active site. The amino acid residues within the binding pocket influence substrate orientation and anchoring, determining the position of the hydroxylation site and hence direct CYP107's catalytic activity. Additionally, the amino acid dynamics within and around the binding pocket determine CYP107's multifunctionality. This study serves as a reference for understanding the structure-function analysis of CYP107 family members precisely and the structure-function analysis of P450 enzymes in general. Finally, this work will aid in the genetic engineering of CYP107 enzymes to produce novel molecules of biotechnological interest.

Computer-Aided Lipase Engineering for Improving Their Stability and Activity in the Food Industry: State of the Art

As some of the most widely used biocatalysts, lipases have exhibited extreme advantages in many processes, such as esterification, amidation, and transesterification reactions, which causes them to be widely used in food industrial production. However, natural lipases have drawbacks in terms of organic solvent resistance, thermostability, selectivity, etc., which limits some of their applications in the field of foods. In this systematic review, the application of lipases in various food processes was summarized. Moreover, the general structure of lipases is discussed in-depth, and the engineering strategies that can be used in lipase engineering are also summarized. The protocols of some classical methods are compared and discussed, which can provide some information about how to choose methods of lipase engineering. Thermostability engineering and solvent tolerance engineering are highlighted in this review, and the basic principles for improving thermostability and solvent tolerance are summarized. In the future, comput er-aided technology should be more emphasized in the investigation of the mechanisms of reactions catalyzed by lipases and guide the engineering of lipases. The engineering of lipase tunnels to improve the diffusion of substrates is also a promising prospect for further enhanced lipase activity and selectivity.

Cytochrome P450 Surface Domains Prevent the β-Carotene Monohydroxylase CYP97H1 of Euglena gracilis from Acting as a Dihydroxylase

Molecular biodiversity results from branched metabolic pathways driven by enzymatic regioselectivities. An additional complexity occurs in metabolites with an internal structural symmetry, offering identical extremities to the enzymes. For example, in the terpene family, β-carotene presents two identical terminal closed-ring structures. Theses cycles can be hydroxylated by cytochrome P450s from the CYP97 family. Two sequential hydroxylations lead first to the formation of monohydroxylated β-cryptoxanthin and subsequently to that of dihydroxylated zeaxanthin. Among the CYP97 dihydroxylases, CYP97H1 from Euglena gracilis has been described as the only monohydroxylase. This study aims to determine which enzymatic domains are involved in this regioselectivity, conferring unique monohydroxylase activity on a substrate offering two identical sites for hydroxylation. We explored the effect of truncations, substitutions and domain swapping with other CYP97 members and found that CYP97H1 harbours a unique N-terminal globular domain. This CYP97H1 N-terminal domain harbours a hydrophobic patch at the entrance of the substrate channel, which is involved in the monohydroxylase activity of CYP97H1. This domain, at the surface of the enzyme, highlights the role of distal and non-catalytic domains in regulating enzyme specificity.

Human Cytochrome P450 1, 2, 3 Families as Pharmacogenes with Emphases on Their Antimalarial and Antituberculosis Drugs and Prevalent African Alleles

Precision medicine gives individuals tailored medical treatment, with the genotype determining the therapeutic strategy, the appropriate dosage, and the likelihood of benefit or toxicity. Cytochrome P450 (CYP) enzyme families 1, 2, and 3 play a pivotal role in eliminating most drugs. Factors that affect CYP function and expression have a major impact on treatment outcomes. Therefore, polymorphisms of these enzymes result in alleles with diverse enzymatic activity and drug metabolism phenotypes. Africa has the highest CYP genetic diversity and also the highest burden of malaria and tuberculosis, and this review presents current general information on CYP enzymes together with variation data concerning antimalarial and antituberculosis drugs, while focusing on the first three CYP families. Afrocentric alleles such as CYP2A6*17, CYP2A6*23, CYP2A6*25, CYP2A6*28, CYP2B6*6, CYP2B6*18, CYP2C8*2, CYP2C9*5, CYP2C9*8, CYP2C9*9, CYP2C19*9, CYP2C19*13, CYP2C19*15, CYP2D6*2, CYP2D6*17, CYP2D6*29, and CYP3A4*15 are implicated in diverse metabolic phenotypes of different antimalarials such as artesunate, mefloquine, quinine, primaquine, and chloroquine. Moreover, CYP3A4, CYP1A1, CYP2C8, CYP2C18, CYP2C19, CYP2J2, and CYP1B1 are implicated in the metabolism of some second-line antituberculosis drugs such as bedaquiline and linezolid. Drug-drug interactions, induction/inhibition, and enzyme polymorphisms that influence the metabolism of antituberculosis, antimalarial, and other drugs, are explored. Moreover, a mapping of Afrocentric missense mutations to CYP structures and a documentation of their known effects provided structural insights, as understanding the mechanism of action of these enzymes and how the different alleles influence enzyme function is invaluable to the advancement of precision medicine.

H2O2-Driven Hydroxylation of Steroids Catalyzed by Cytochrome P450 CYP105D18: Exploration of the Substrate Access Channel

CYP105D18 supports H2O2 as an oxygen surrogate for catalysis well and shows high H2O2 resistance capacity. We report the hydroxylation of different steroids using H2O2 as a cosubstrate. Testosterone was regiospecifically hydroxylated to 2β-hydroxytestosterone. Based on the experimental data and molecular docking, we predicted that hydroxylation of methyl testosterone and nandrolone would occur at position 2 in the A-ring, while hydroxylation of androstenedione and adrenosterone was predicted to occur in the B-ring. Further, structure-guided rational design of the substrate access channel was performed with the mutagenesis of residues S63, R82, and F184. Among the mutants, S63A showed a marked decrease in product formation, while F184A showed a significant increase in product formation in testosterone, nandrolone, methyl testosterone, androstenedione, and adrenosterone. The catalytic efficiency (kcat/Km) toward testosterone was increased 1.36-fold in the F184A mutant over that in the wild-type enzyme. These findings might facilitate the potential use of CYP105D18 and further engineering to establish the basis of biotechnological applications. IMPORTANCE The structural modification of steroids is a challenging chemical reaction. Modifying the core ring and the side chain improves the biological activity of steroids. In particular, bacterial cytochrome P450s are used as promiscuous enzymes for the activation of nonreactive carbons of steroids. In the present work, we reported the H2O2-mediated hydroxylation of steroids by CYP105D18, which also overcomes the use of expensive cofactors. Further, exploring the substrate access channel and modifying the bulky amino acid F184A increase substrate conversion while modifying the substrate recognizing amino acid S63 markedly decreases product formation. Exploring the substrate access channel and the rational design of CYP105D18 can improve the substrate conversion, which facilitates the engineering of P450s for industrial application.

Geometry-Based versus Small-Molecule Tracking Method for Tunnel Identification: Benefits and Pitfalls

Different methods for tunnel identification, geometry-based and small-molecule tracking approaches, were compared to provide their benefits and pitfalls. Results obtained for both crystal structures and molecular dynamics (MD) simulations were analyzed to investigate if a more computationally demanding method would be beneficial. Careful examination of the results is essential for the low-diameter tunnel description, and assessment of the tunnel functionality based only on their geometrical parameters is challenging. We showed that the small-molecule tracking approach can provide a detailed description of the system; however, it can also be the most computationally demanding.

Computational Investigation of Ligand Binding of Flavonoids in Cytochrome P450 Receptors

AIM

The cytochrome P450 enzymes play a significant role in regulating cellular and physiological processes by activating endogenous compounds. They also play an essential role in the detoxification process of xenobiotics. Flavonoids belong to a class of polyphenols found in food, such as vegetables, red wine, beer, and fruits, which modulate biological functions in the body.

METHODS

The inhibition of CYP1A1 and CYP1B1 using nutritional sources has been reported as a strategy for cancer prevention. This study investigated the interactions of selected flavonoids binding to the cytochrome P450 enzymes (CYP1A1 and CYP1B1) and their ADMET properties in silico. From docking studies, our findings showed flavonoids, isorhamnetin and pedalitin, to have the strongest binding energies in the crystal structures 6DWM and 6IQ5.

RESULTS

The amino acid residues Asp 313 and Phe 224 in 6DWM interacted with all the ligands investigated, and Ala 330 in 6IQ5 interacted with all the ligands examined. The ligands did not violate any drug-likeness parameters.

CONCLUSION

These data suggest roles for isorhamnetin and pedalitin as potential precursors for natural product- derived therapies.

Synthesis, biological evaluation and computational studies of pyrazole derivatives as Mycobacterium tuberculosis CYP121A1 inhibitors

A series of imidazole and triazole diarylpyrazole derivatives were prepared using an efficient 5-step synthetic scheme and evaluated for binding affinity with Mycobacterium tuberculosis (Mtb) CYP121A1 and antimycobacterial activity against Mtb H37Rv. Antimycobacterial susceptibility was measured using the spot-culture growth inhibition assay (SPOTi): the imidazoles displayed minimum inhibitory concentration (MIC₉₀) in the range of 3.95-12.03 μg mL^-1 (10.07-33.19 μM) with 11f the most active, while the triazoles displayed MIC₉₀ in the range of 4.35-25.63 μg mL^-1 (11.88-70.53 μM) with 12b the most active. Assessment of binding affinity using UV-vis spectroscopy showed that for the imidazole series, the propyloxy (11f) and isopropyloxy (11h) derivatives of the 4-chloroaryl pyrazoles displayed Mtb CYP121A1 type II binding affinity with K _d 11.73 and 17.72 μM respectively compared with the natural substrate cYY (K _d 12.28 μM), while in the triazole series, only the methoxy substitution with the 4-chloroaryl pyrazole (12b) showed good type II Mtb CYP121A1 binding affinity (K _d 5.13 μM). Protein-detected 1D ¹⁹F-NMR spectroscopy as an orthogonal strategy was used to evaluate ligand binding independent of perturbations at the haem. For imidazole and triazole compounds, perturbations were more intense than cYY indicating tighter binding and confirming that ligand coordination occurs in the substrate-binding pocket despite very modest changes in UV-vis absorbance, consistent with computational studies and the demonstrated potential anti-tuberculosis properties of these compounds.

Engineering of Microbial Substrate Promiscuous CYP105A5 for Improving the Flavonoid Hydroxylation

Bacterial cytochrome P450 (CYP) enzymes are versatile biocatalysts that are responsible for the biotransformation of diverse endogenous substances. CYP105A5 from Streptomyces sp. showed substrate flexibility with different flavonoids and was able to catalyze O-demethylation of biochanin A, regioselective C3′-hydroxylation of daidzein, genistein, and naringenin, and additional C8-hydroxylation for daidzein using heterologous redox partners putidaredoxin and putidaredoxin reductase. By rational design of substrate-binding pocket based on experimental data, homology modeling, and molecular docking analysis, we enhanced the product formation rate of flavonoids. The double mutant L100A/I302A and L100A/I408N exhibited greatly enhanced in vivo conversion rates for flavonoid hydroxylation. Particularly, the L100A/I302A mutant’s kcat/Km values and in vivo conversion rate increased by 1.68-fold and 2.57-fold, respectively, for naringenin. Overall, our result might facilitate the potential use of CYP105A5 for future modification and application in whole-cell biocatalysts for the production of valuable polyphenols.

Design principles for site-selective hydroxylation by a Rieske oxygenase

Rieske oxygenases exploit the reactivity of iron to perform chemically challenging C–H bond functionalization reactions. Thus far, only a handful of Rieske oxygenases have been structurally characterized and remarkably little information exists regarding how these enzymes use a common architecture and set of metallocenters to facilitate a diverse range of reactions. Herein, we detail how two Rieske oxygenases SxtT and GxtA use different protein regions to influence the site-selectivity of their catalyzed monohydroxylation reactions. We present high resolution crystal structures of SxtT and GxtA with the native β-saxitoxinol and saxitoxin substrates bound in addition to a Xenon-pressurized structure of GxtA that reveals the location of a substrate access tunnel to the active site. Ultimately, this structural information allowed for the identification of six residues distributed between three regions of SxtT that together control the selectivity of the C–H hydroxylation event. Substitution of these residues produces a SxtT variant that is fully adapted to exhibit the non-native site-selectivity and substrate scope of GxtA. Importantly, we also found that these selectivity regions are conserved in other structurally characterized Rieske oxygenases, providing a framework for predictively repurposing and manipulating Rieske oxygenases as biocatalysts.

SxtT and GxtA are Rieske oxygenases that are involved in paralytic shellfish toxin biosynthesis and catalyze monohydroxylation reactions at different positions on the toxin scaffold. Here, the authors present crystal structures of SxtT and GxtA with the native substrates β-saxitoxinol and saxitoxin as well as a Xenon-pressurized structure of GxtA, which reveal a substrate access tunnel to the active site. Through structure-based mutagenesis studies the authors identify six residues in three different protein regions that determine the substrate specificity and site selectivity of SxtT and GxtA. These findings will aid the rational engineering of other Rieske oxygenases.

A Conserved Allosteric Site on Drug-Metabolizing CYPs: A Systematic Computational Assessment

Cytochrome P450 enzymes (CYPs) are the largest group of enzymes involved in human drug metabolism. Ligand tunnels connect their active site buried at the core of the membrane-anchored protein to the surrounding solvent environment. Recently, evidence of a superficial allosteric site, here denoted as hotspot 1 (H1), involved in the regulation of ligand access in a soluble prokaryotic CYP emerged. Here, we applied multi-scale computational modeling techniques to study the conservation and functionality of this allosteric site in the nine most relevant mammalian CYPs responsible for approximately 70% of drug metabolism. In total, we systematically analyzed over 44 μs of trajectories from conventional MD, cosolvent MD, and metadynamics simulations. Our bioinformatic analysis and simulations with organic probe molecules revealed the site to be well conserved in the CYP2 family with the exception of CYP2E1. In the presence of a ligand bound to the H1 site, we could observe an enlargement of a ligand tunnel in several members of the CYP2 family. Further, we could detect the facilitation of ligand translocation by H1 interactions with statistical significance in CYP2C8 and CYP2D6, even though all other enzymes except for CYP2C19, CYP2E1, and CYP3A4 presented a similar trend. As the detailed comprehension of ligand access and egress phenomena remains one of the most relevant challenges in the field, this work contributes to its elucidation and ultimately helps in estimating the selectivity of metabolic transformations using computational techniques.

Structure-Activity Relationship and Mode-Of-Action Studies Highlight 1-(4-Biphenylylmethyl)-1H-imidazole-Derived Small Molecules as Potent CYP121 Inhibitors

CYP121 of Mycobacterium tuberculosis (Mtb) is an essential target for the development of novel potent drugs against tuberculosis (TB). Besides known antifungal azoles, further compounds of the azole class were recently identified as CYP121 inhibitors with antimycobacterial activity. Herein, we report the screening of a similarity-oriented library based on the former hit compound, the evaluation of affinity toward CYP121, and activity against M. bovis BCG. The results enabled a comprehensive SAR study, which was extended through the synthesis of promising compounds and led to the identification of favorable features for affinity and/or activity and hit compounds with 2.7-fold improved potency. Mode of action studies show that the hit compounds inhibit substrate conversion and highlighted CYP121 as the main antimycobacterial target of our compounds. Exemplified complex crystal structures of CYP121 with three inhibitors reveal a common binding site. Engaging in both hydrophobic interactions as well as hydrogen bonding to the sixth iron ligand, our compounds block a solvent channel leading to the active site heme. Additionally, we report the first CYP inhibitors that are able to reduce the intracellular replication of M. bovis BCG in macrophages, emphasizing their potential as future drug candidates against TB.

Calculation of CYP450 protein-ligand binding and dissociation free energy paths

The function of an enzyme depends on its dynamic structure, and the catalytic mechanism has long been an active focus of research. The principle for interpreting protein selectivity and fidelity stems from optimization of the active site upon protein-substrate complexation, i.e., a lock-and-key configuration, on which most protein-substrate molecule binding recognition, and hence drug discovery, relies. Yet another thought has been to incorporate the protein folding interior tunnels for stereo- and regio-selectivity along the protein-substrate or protein-ligand/inhibitor binding process. Free energy calculations provide valuable information for molecular recognition and protein-ligand binding dynamics and kinetics. In this study, we focused on the kinetics of cytochrome P450 proteins (CYP450s) and the protein interior tunnel structure-dynamics relationship in terms of the substrate binding and leaving mechanism. A case in point is given by the prostaglandin H₂ (PGH₂) homologous isomerase of prostacyclin synthase. To calculate the reactant and product traversing the tunnels to and from the heme site, the free energy paths and tunnel potentials of mean force are constructed from steered molecular dynamics simulations and adaptive basing force umbrella sampling simulations. We explore the binding tunnels and critical residue lining characteristics for the ligand traverse and the underlying mechanism of CYP450 activity. Our theoretical analysis provides insights into the decisive role of the substrate tunnel binding process of the CYP450 mechanism and may be useful in drug design and protein engineering contexts.

F/G Region Rigidity is Inversely Correlated to Substrate Promiscuity of Human CYP Isoforms Involved in Metabolism

Of 57 human cytochrome P450 (CYP) enzymes, 12 metabolize 90% of xenobiotics. To our knowledge, no study has addressed the relation between enzyme dynamics and substrate promiscuity for more than three CYPs. Here, we show by constraint dilution simulations with the Constraint Network Analysis for the 12 isoforms that structural rigidity of the F/G region is significantly inversely correlated to the enzymes' substrate promiscuity. This highlights the functional importance of structural dynamics of the substrate tunnel.

Molecular Dynamics Simulations of a Cytochrome P450 from Tepidiphilus thermophilus (P450-TT) Reveal How Its Substrate-Binding Channel Opens

Cytochrome P450s (P450) are important enzymes in biology with useful biochemical reactions in, for instance, drug and xenobiotics metabolisms, biotechnology, and health. Recently, the crystal structure of a new member of the CYP116B family has been resolved. This enzyme is a cytochrome P450 (CYP116B46) from Tepidiphilus thermophilus (P450-TT) and has potential for the oxy-functionalization of organic molecules such as fatty acids, terpenes, steroids, and statins. However, it was thought that the opening to its hitherto identified substrate channel was too small to allow organic molecules to enter. To investigate this, we performed molecular dynamics simulations on the enzyme. The results suggest that the crystal structure is not relaxed, possibly due to crystal packing effects, and that its tunnel structure is constrained. In addition, the simulations revealed two key amino acid residues at the mouth of the channel; a glutamyl and an arginyl. The glutamyl’s side chain tightens and relaxes the opening to the channel in conjunction with the arginyl’s, though the latter’s side chain is less dramatically changed after the initial relaxation of its conformations. Additionally, it was observed that the effect of increased temperature did not considerably affect the dynamics of the enzyme fold, including the relative solvent accessibility of the amino acid residues that make up the substrate channel wall even as compared to the changes that occurred at room temperature. Interestingly, the substrate channel became distinguishable as a prominent tunnel that is likely to accommodate small- to medium-sized organic molecules for bioconversions. That is, P450-TT has the ability to pass appropriate organic substrates to its active site through its elaborate substrate channel, and notably, is able to control or gate any molecules at the opening to this channel.

Uncovering of cytochrome P450 anatomy by SecStrAnnotator

Protein structural families are groups of homologous proteins defined by the organization of secondary structure elements (SSEs). Nowadays, many families contain vast numbers of structures, and the SSEs can help to orient within them. Communities around specific protein families have even developed specialized SSE annotations, always assigning the same name to the equivalent SSEs in homologous proteins. A detailed analysis of the groups of equivalent SSEs provides an overview of the studied family and enriches the analysis of any particular protein at hand. We developed a workflow for the analysis of the secondary structure anatomy of a protein family. We applied this analysis to the model family of cytochromes P450 (CYPs)-a family of important biotransformation enzymes with a community-wide used SSE annotation. We report the occurrence, typical length and amino acid sequence for the equivalent SSE groups, the conservation/variability of these properties and relationship to the substrate recognition sites. We also suggest a generic residue numbering scheme for the CYP family. Comparing the bacterial and eukaryotic part of the family highlights the significant differences and reveals a well-known anomalous group of bacterial CYPs with some typically eukaryotic features. Our workflow for SSE annotation for CYP and other families can be freely used at address https://sestra.ncbr.muni.cz .

Structural Basis for Selective Oxidation of Phosphorylated Ethylphenols by Cytochrome P450 Monooxygenase CreJ

Selective oxidation of C-H bonds in alkylphenols holds great significance for not only structural derivatization in pharma- and biomanufacturing but also biological degradation of these toxic chemicals in environmental protection. A unique chemomimetic biocatalytic system using enzymes from a p-cresol biodegradation pathway has recently been developed. As the central biocatalyst, the cytochrome P450 monooxygenase CreJ oxidizes diverse p- and m-alkylphenol phosphates with perfect stereoselectivity at different efficiencies. However, the mechanism of regio- and stereoselectivity of this chemomimetic biocatalytic system remained unclear. Here, using p- and m-ethylphenol substrates, we elucidate the CreJ-catalyzed key steps for selective oxidations. The crystal structure of CreJ in complex with m-ethylphenol phosphate was solved and compared with its complex structure with p-ethylphenol phosphate isomer. The results indicate that the conformational changes of substrate-binding residues are slight, while the substrate promiscuity is achieved mainly by the available space in the catalytic cavity. Moreover, the catalytic preferences of regio- and stereoselective hydroxylation for the two ethylphenol substrates is explored by molecular dynamics simulations. The ethyl groups in the complexes display different flexibilities, and the distances of the active oxygen to H _pro-S and H _pro-R of methylene agree with the experimental stereoselectivity. The regioselectivity can be explained by the distances and bond dissociation energy. These results provide not only the mechanistic insights into CreJ regio- and stereoselectivity but also the structural basis for further P450 enzyme design and engineering.IMPORTANCE The key cytochrome P450 monooxygenase CreJ showed excellent regio- and stereoselectivity in the oxidation of various alkylphenol substrates. C-H bond functionalization of these toxic alkylphenols holds great significance for both biological degradation of these environmental chemicals and production of value-added structural derivatives in pharmaceutical and biochemical industries. Our results, combined with in vitro enzymatic assays, crystal structure determination of enzyme-substrate complex, and molecular dynamics simulations, provide not only significant mechanism elucidation of the regio- and stereoselective catalyzation mediated by CreJ but also the promising directions for future engineering efforts of this enzyme toward more useful products. It also has great extendable potential to couple this multifunctional P450 enzyme with other biocatalysts (e.g., hydroxyl-based glycosylase) to access more alkylphenol-derived high-value chemicals through environment-friendly biocatalysis and biotransformation.

Reductive Cytochrome P450 Reactions and Their Potential Role in Bioremediation

Cytochrome P450 enzymes, or P450s, are haem monooxygenases renowned for their ability to insert one atom from molecular oxygen into an exceptionally broad range of substrates while reducing the other atom to water. However, some substrates including many organohalide and nitro compounds present little or no opportunity for oxidation. Under hypoxic conditions P450s can perform reductive reactions, contributing electrons to drive reductive elimination reactions. P450s can catalyse dehalogenation and denitration of a range of environmentally persistent pollutants including halogenated hydrocarbons and nitroamine explosives. P450-mediated reductive dehalogenations were first discovered in the context of human pharmacology but have since been observed in a variety of organisms. Additionally, P450-mediated reductive denitration of synthetic explosives has been discovered in bacteria that inhabit contaminated soils. This review will examine the distribution of P450-mediated reductive dehalogenations and denitrations in nature and discuss synthetic biology approaches to developing P450-based reagents for bioremediation.

ProPores2: Web Service and Stand-Alone Tool for Identifying, Manipulating, and Visualizing Pores in Protein Structures

Surface pockets, cavities, and tunnels in the 3D structures of proteins play integral functional roles such as enabling enzymatic catalysis, ligand binding, or transport of ions or small molecules across biomembranes. ProPores2 facilitates understanding and analysis of these processes by identifying pores and lining residues, determining their axes, and opening closed connections via side-chain rotation. The fast stand-alone tool introduces a novel mode for pore identification, improved axis determination, and additional features such as parallel batch processing and a graphical user interface. The new web service features an integrated and customizable protein viewer with an option to analyze and view more than one structure at once. This feature facilitates side-by-side comparisons of pores in different conformations of the same protein or of identified pores before and after opening gates within the same protein. ProPores2 is freely and publicly available at https://service.bioinformatik.uni-saarland.de/propores.

Tryptophan-75 Is a Low-Energy Channel-Gating Residue that Facilitates Substrate Egress/Access in Cytochrome P450 2D6

CYP2D6 is a major drug metabolizing enzyme with a buried active site. Channels leading to the active site from various enzyme surfaces are believed to facilitate ligand egress and access to the active site. The present study used molecular dynamics (MD) and in vitro studies with CYP2D6*1 and a Trp75-to-Ala mutant to examine channel gating in CYP2D6 by Trp75. MD simulations measured energy landscapes of Trp75 conformations and simulated substrate passage within channel 2b using bufuralol as a model substrate. Trp75 alternated between multiple stable states that supported substrate transport along channel 2b with low-energy barriers between states (∼ -1 kcal/mol). Trp75 conformations were stabilized primarily by hydrogen bonding between Trp75 and Glu222, Asn226, Ala225, or Gln72. Energy barriers were low between Trp75 conformations, allowing Trp75 to easily move between various conformations over time and to function in both binding to and moving substrates in the 2b channel of CYP2D6. Michaelis-Menten kinetic studies completed with purified enzyme in a reconstituted system showed overall reduced enzyme efficiency for metabolism of bufuralol and dextromethorphan by the Trp75Ala mutant compared with CYP2D6*1. In stopped-flow measurements, k off for dextromethorphan was decreased in the absence of Trp75. Our results support a role for Trp75 in substrate shuttling to the active site of CYP2D6. SIGNIFICANCE STATEMENT: Using combined molecular dynamics and in vitro assays, this study shows for the first time a role for Trp75 as a channel entrance gating residue in the mechanism of substrate binding/unbinding in CYP2D6. Energy landscapes derived from molecular dynamics were used to quantitate the strength of gating, and kinetics assays showed the impact on enzyme efficiency and k off of a Trp75Ala mutation.

Paroxetine-Overview of the Molecular Mechanisms of Action

In the 21st century and especially during a pandemic, the diagnosis and treatment of depression is an essential part of the daily practice of many family doctors. It mainly affects patients in the age category 15–44 years, regardless of gender. Anxiety disorders are often diagnosed in children and adolescents. Social phobias can account for up to 13% of these diagnoses. Social anxiety manifests itself in fear of negative social assessment and humiliation, which disrupts the quality of social functioning. Treatment of the above-mentioned disorders is based on psychotherapy and pharmacotherapy. Serious side effects or mortality from antidepressant drug overdose are currently rare. Recent studies indicate that paroxetine (ATC code: N06AB), belonging to the selective serotonin reuptake inhibitors, has promising therapeutic effects and is used off-label in children and adolescents. The purpose of this review is to describe the interaction of paroxetine with several molecular targets in various points of view including the basic chemical and pharmaceutical properties. The central point of the review is focused on the pharmacodynamic analysis based on the molecular mechanism of binding paroxetine to various therapeutic targets.

Advances in enzymatic oxyfunctionalization of aliphatic compounds

Selective oxyfunctionalizations of aliphatic compounds are difficult chemical reactions, where enzymes can play an important role due to their stereo- and regio-selectivity and operation under mild reaction conditions. P450 monooxygenases are well-known biocatalysts that mediate oxyfunctionalization reactions in different living organisms (from bacteria to humans). Unspecific peroxygenases (UPOs), discovered in fungi, have arisen as "dream biocatalysts" of great biotechnological interest because they catalyze the oxyfunctionalization of aliphatic and aromatic compounds, avoiding the necessity of expensive cofactors and regeneration systems, and only depending on H₂O₂ for their catalysis. Here, we summarize recent advances in aliphatic oxyfunctionalization reactions by UPOs, as well as the molecular determinants of the enzyme structures responsible for their activities, emphasizing the differences found between well-known P450s and the novel fungal peroxygenases.

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.

Molecular determinant of substrate binding and specificity of cytochrome P450 2J2

Cytochrome P450 2J2 (CYP2J2) is responsible for the epoxidation of endogenous arachidonic acid, and is involved in the metabolism of exogenous drugs. To date, no crystal structure of CYP2J2 is available, and the proposed structural basis for the substrate recognition and specificity in CYP2J2 varies with the structural models developed using different computational protocols. In this study, we developed a new structural model of CYP2J2, and explored its sensitivity to substrate binding by molecular dynamics simulations of the interactions with chemically similar fluorescent probes. Our results showed that the induced-fit binding of these probes led to the preferred active poses ready for the catalysis by CYP2J2. Divergent conformational dynamics of CYP2J2 due to the binding of each probe were observed. However, a stable hydrophobic clamp composed of residues I127, F310, A311, V380, and I487 was identified to restrict any substrate access to the active site of CYP2J2. Molecular docking of a series of compounds including amiodarone, astemizole, danazol, ebastine, ketoconazole, terfenadine, terfenadone, and arachidonic acid to CYP2J2 confirmed the role of those residues in determining substrate binding and specificity of CYP2J2. In addition to the flexibility of CYP2J2, the present work also identified other factors such as electrostatic potential in the vicinity of the active site, and substrate strain energy and property that have implications for the interpretation of CYP2J2 metabolism.

CYP2J2 Molecular Recognition: A New Axis for Therapeutic Design

Cytochrome P450 (CYP) epoxygenases are a special subset of heme-containing CYP enzymes capable of performing the epoxidation of polyunsaturated fatty acids (PUFA) and the metabolism of xenobiotics. This dual functionality positions epoxygenases along a metabolic crossroad. Therefore, structure-function studies are critical for understanding their role in bioactive oxy-lipid synthesis, drug-PUFA interactions, and for designing therapeutics that directly target the epoxygenases. To better exploit CYP epoxygenases as therapeutic targets, there is a need for improved understanding of epoxygenase structure-function. Of the characterized epoxygenases, human CYP2J2 stands out as a potential target because of its role in cardiovascular physiology. In this review, the early research on the discovery and activity of epoxygenases is contextualized to more recent advances in CYP epoxygenase enzymology with respect to PUFA and drug metabolism. Additionally, this review employs CYP2J2 epoxygenase as a model system to highlight both the seminal works and recent advances in epoxygenase enzymology. Herein we cover CYP2J2's interactions with PUFAs and xenobiotics, its tissue-specific physiological roles in diseased states, and its structural features that enable epoxygenase function. Additionally, the enumeration of research on CYP2J2 identifies the future needs for the molecular characterization of CYP2J2 to enable a new axis of therapeutic design.

Synthesis, Optimization, Antifungal Activity, Selectivity, and CYP51 Binding of New 2-Aryl-3-azolyl-1-indolyl-propan-2-ols

A series of 2-aryl-3-azolyl-1-indolyl-propan-2-ols was designed as new analogs of fluconazole (FLC) by replacing one of its two triazole moieties by an indole scaffold. Two different chemical approaches were then developed. The first one, in seven steps, involved the synthesis of the key intermediate 1-(1H-benzotriazol-1-yl)methyl-1H-indole and the final opening of oxiranes by imidazole or 1H-1,2,4-triazole. The second route allowed access to the target compounds in only three steps, this time with the ring opening by indole and analogs. Twenty azole derivatives were tested against Candida albicans and other Candida species. The enantiomers of the best anti-Candida compound, 2-(2,4-dichlorophenyl)-3-(1H-indol-1-yl)-1-(1H-1,2,4-triazol-1-yl)-propan-2-ol (8g), were analyzed by X-ray diffraction to determine their absolute configuration. The (−)-8g enantiomer (Minimum inhibitory concentration (MIC) = IC₈₀ = 0.000256 µg/mL on C. albicans CA98001) was found with the S-absolute configuration. In contrast the (+)-8g enantiomer was found with the R-absolute configuration (MIC = 0.023 µg/mL on C. albicans CA98001). By comparison, the MIC value for FLC was determined as 0.020 µg/mL for the same clinical isolate. Additionally, molecular docking calculations and molecular dynamics simulations were carried out using a crystal structure of Candida albicans lanosterol 14α-demethylase (CaCYP51). The (−)-(S)-8g enantiomer aligned with the positioning of posaconazole within both the heme and access channel binding sites, which was consistent with its biological results. All target compounds have been also studied against human fetal lung fibroblast (MRC-5) cells. Finally, the selectivity of four compounds on a panel of human P450-dependent enzymes (CYP19, CYP17, CYP26A1, CYP11B1, and CYP11B2) was investigated.

Comparative Dynamics and Functional Mechanisms of the CYP17A1 Tunnels Regulated by Ligand Binding

As an important member of cytochrome P450 (CYP) enzymes, CYP17A1 is a dual-function monooxygenase with a critical role in the synthesis of many human steroid hormones, making it an attractive therapeutic target. The emerging structural information about CYP17A1 and the growing number of inhibitors for these enzymes call for a systematic strategy to delineate and classify mechanisms of ligand transport through tunnels that control catalytic activity. In this work, we applied an integrated computational strategy to different CYP17A1 systems with a panel of ligands to systematically study at the atomic level the mechanism of ligand-binding and tunneling dynamics. Atomistic simulations and binding free energy computations identify the dynamics of dominant tunnels and characterize energetic properties of critical residues responsible for ligand binding. The common transporting pathways including S, 3, and 2c tunnels were identified in CYP17A1 binding systems, while the 2c tunnel is a newly formed pathway upon ligand binding. We employed and integrated several computational approaches including the analysis of functional motions and sequence conservation, atomistic modeling of dynamic residue interaction networks, and perturbation response scanning analysis to dissect ligand tunneling mechanisms. The results revealed the hinge-binding and sliding motions as main functional modes of the tunnel dynamic, and a group of mediating residues as key regulators of tunnel conformational dynamics and allosteric communications. We have also examined and quantified the mutational effects on the tunnel composition, conformational dynamics, and long-range allosteric behavior. The results of this investigation are fully consistent with the experimental data, providing novel rationale to the experiments and offering valuable insights into the relationships between the structure and function of the channel networks and a robust atomistic model of activation mechanisms and allosteric interactions in CYP enzymes.

Carbon-fluorine bond cleavage mediated by metalloenzymes

Fluorochemicals are a widely distributed class of compounds and have been utilized across a wide range of industries for decades. Given the environmental toxicity and adverse health threats of some fluorochemicals, the development of new methods for their decomposition is significant to public health. However, the carbon-fluorine (C-F) bond is among the most chemically robust bonds; consequently, the degradation of fluorinated hydrocarbons is exceptionally difficult. Here, metalloenzymes that catalyze the cleavage of this chemically challenging bond are reviewed. These enzymes include histidine-ligated heme-dependent dehaloperoxidase and tyrosine hydroxylase, thiolate-ligated heme-dependent cytochrome P450, and four nonheme oxygenases, namely, tetrahydrobiopterin-dependent aromatic amino acid hydroxylase, 2-oxoglutarate-dependent hydroxylase, Rieske dioxygenase, and thiol dioxygenase. While much of the literature regarding the aforementioned enzymes highlights their ability to catalyze C-H bond activation and functionalization, in many cases, the C-F bond cleavage has been shown to occur on fluorinated substrates. A copper-dependent laccase-mediated system representing an unnatural radical defluorination approach is also described. Detailed discussions on the structure-function relationships and catalytic mechanisms provide insights into biocatalytic defluorination, which may inspire drug design considerations and environmental remediation of halogenated contaminants.

Capturing Protein-Ligand Recognition Pathways in Coarse-Grained Simulation

Protein-ligand recognition is dynamic and complex. A key approach in deciphering the mechanism underlying the recognition process is to capture the kinetic process of the ligand in its act of binding to its designated protein cavity. Toward this end, ultralong all-atom molecular dynamics simulation has recently emerged as a popular method of choice because of its ability to record these events at high spatial and temporal resolution. However, success via this route comes at an exorbitant computational cost. Herein, we demonstrate that coarse-grained models of the protein, when systematically optimized to maintain its tertiary fold, can capture the complete process of spontaneous protein-ligand binding from bulk media to the cavity at crystallographic precision and within wall clock time that is orders of magnitude shorter than that of all-atom simulations. The exhaustive sampling of ligand exploration in protein and solvent, harnessed by coarse-grained simulation, leads to elucidation of new ligand recognition pathways and discovery of non-native binding poses.

Structural basis for plant lutein biosynthesis from α-carotene

Two cytochrome P450 enzymes, CYP97A3 and CYP97C1, catalyze hydroxylations of the β- and ε-rings of α-carotene to produce lutein. Chirality is introduced at the C-3 atom of both rings, and the reactions are both pro-3R-stereospecific. We determined the crystal structures of CYP97A3 in substrate-free and complex forms with a nonnatural substrate and the structure of CYP97C1 in a detergent-bound form. The structures of CYP97A3 in different states show the substrate channel and the structure of CYP97C1 bound with octylthioglucoside confirms the binding site for the carotenoid substrate. Biochemical assays confirm that the ferredoxin-NADP⁺ reductase (FNR)-ferredoxin pair is used as the redox partner. Details of the pro-3R stereospecificity are revealed in the retinal-bound CYP97A3 structure. Further analysis indicates that the CYP97B clan bears similarity to the β-ring-specific CYP97A clan. Overall, our research describes the molecular basis for the last steps of lutein biosynthesis.

Synthetic Derivatives of (+)-epi-α-Bisabolol Are Formed by Mammalian Cytochromes P450 Expressed in a Yeast Reconstituted Pathway

Identification of the enzyme(s) involved in complex biosynthetic pathways can be challenging. An alternative approach might be to deliberately diverge from the original natural enzyme source and use promiscuous enzymes from other organisms. In this paper, we have tested the ability of a series of human and animal cytochromes P450 involved in xenobiotic detoxification to generate derivatives of (+)-epi-α-bisabolol and attempt to produce the direct precursor of hernandulcin, a sweetener from Lippia dulcis for which the last enzymatic steps are unknown. Screening steps were implemented in vivo in S. cerevisiae optimized for the biosynthesis of oxidized derivatives of (+)-epi-α-bisabolol by coexpressing two key enzymes: the (+)-epi-α-bisabolol synthase and the NADPH cytochrome P450 reductase. Five out of 25 cytochromes P450 were capable of producing new hydroxylated regioisomers of (+)-epi-α-bisabolol. Of the new oxidized bisabolol products, the structure of one compound, 14-hydroxy-(+)-epi-α-bisabolol, was fully elucidated by NMR while the probable structure of the second product was determined. In parallel, the production of (+)-epi-α-bisabolol derivatives was enhanced through the addition of a supplementary genomic copy of (+)-epi-α-bisabolol synthase that augmented the final titer of hydroxylated product to 64 mg/L. We thus demonstrate that promiscuous drug metabolism cytochromes P450 can be used to produce novel compounds from a terpene scaffold.

Functional impact of cytochrome P450 3A (CYP3A) missense variants in cattle

Cytochrome P450 3A is the most important CYP subfamily in humans, and CYP3A4/CYP3A5 genetic variants contribute to inter-individual variability in drug metabolism. However, no information is available for bovine CYP3A (bCYP3A). Here we described bCYP3A missense single nucleotide variants (SNVs) and evaluated their functional effects. CYP3A28, CYP3A38 and CYP3A48 missense SNVs were identified in 300 bulls of Piedmontese breed through targeted sequencing. Wild-type and mutant bCYP3A cDNAs were cloned and expressed in V79 cells. CYP3A-dependent oxidative metabolism of testosterone (TST) and nifedipine (NIF) was assessed by LC-MS/MS. Finally, SNVs functional impact on TST hydroxylation was measured ex vivo in liver microsomes from individually genotyped animals. Thirteen missense SNVs were identified and validated. Five variants showed differences in CYP3A catalytic activity: three CYP3A28 SNVs reduced TST 6β-hydroxylation; one CYP3A38 variant increased TST 16β-hydroxylation, while a CYP3A48 SNV showed enhanced NIF oxidation. Individuals homozygous for rs384467435 SNV showed a reduced TST 6β-hydroxylation. Molecular modelling showed that most of SNVs were distal to CYP3A active site, suggesting indirect effects on the catalytic activity. Collectively, these findings demonstrate the importance of pharmacogenetics studies in veterinary species and suggest bCYP3A genotype variation might affect the fate of xenobiotics in food-producing species such as cattle.

Spontaneous Ligand Access Events to Membrane-Bound Cytochrome P450 2D6 Sampled at Atomic Resolution

The membrane-anchored enzyme Cytochrome P450 2D6 (CYP2D6) is involved in the metabolism of around 25% of marketed drugs and its metabolic performance shows a high interindividual variation. While it was suggested that ligands access the buried active site of the enzyme from the membrane, no proof from unbiased simulations has been provided to support this hypothesis. Laboratory experiments fail to capture the access process which is suspected to influence binding kinetics. Here, we applied unbiased molecular dynamics (MD) simulations to investigate the access of ligands to wild-type CYP2D6, as well as the allelic variant CYP2D6*53. In multiple simulations, substrates accessed the active site of the enzyme from the protein-membrane interface to ultimately adopt a conformation that would allow a metabolic reaction. We propose the necessary steps for ligand access and the results suggest that the increased metabolic activity of CYP2D6*53 might be caused by a facilitated ligand uptake.

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.

Ligand Pathways in Nuclear Receptors

Nuclear receptors (NRs) are ligand-inducible transcription factors that play an essential role in a multitude of physiological processes as well as diseases, rendering them attractive drug targets. Crystal structures revealed the binding site of NRs to be buried in the core of the protein, with no obvious route for ligands to access this cavity. The process of ligand binding is known to be an often-neglected contributor to the efficacy of drug candidates and is thought to influence the selectivity and specificity of NRs. While experimental methods generally fail to highlight the dynamic processes of ligand access or egress on the atomistic scale, computational methods have provided fundamental insight into the pathways connecting the buried binding pocket to the surrounding environment. Methods based on molecular dynamics (MD) and Monte Carlo simulations have been applied to identify pathways and quantify their capability to transport ligands. Here, we systematically review findings of more than 20 years of research in the field, including the applied methodology and controversies. Further, we establish a unified nomenclature to describe the pathways with respect to their location relative to protein secondary structure elements and summarize findings relevant to drug design. Lastly, we discuss the effect of NR interaction partners such as coactivators and corepressors, as well as mutations on the pathways.