The ins and outs of cytochrome P450s

Ranking Single Fluorescent Protein-Based Calcium Biosensor Performance by Molecular Dynamics Simulations

Genetically encoded fluorescent biosensors (GEFBs) have become indispensable tools for visualizing biological processes in vivo. A typical GEFB is composed of a sensory domain (SD) that undergoes a conformational change upon ligand binding or enzymatic reaction; the SD is genetically fused with a fluorescent protein (FP). The changes in the SD allosterically modulate the chromophore environment whose spectral properties are changed. Single fluorescent (FP)-based biosensors, a subclass of GEFBs, offer a simple experimental setup; they are easy to produce in living cells, structurally stable, and simple to use due to their single-wavelength operation. However, they pose a significant challenge for structure optimization, especially concerning the length and residue content of linkers between the FP and SD, which affect how well the chromophore responds to conformational change in the SD. In this work, we use all-atom molecular dynamics simulations to analyze the dynamic properties of a series of calmodulin-based calcium biosensors, all with different FP-SD interaction interfaces and varying degrees of calcium binding-dependent fluorescence change. Our results indicate that biosensor performance can be predicted based on distribution of water molecules around the chromophore and shifts in hydrogen bond occupancies between the ligand-bound and ligand-free sensor structures.

Functional significance of CYP2B6 gene rare allelic variants identified in Japanese individuals

Cytochrome P450 2B6 (CYP2B6) catalyzes the metabolism of many drugs, including efavirenz and propofol. Genetic polymorphisms in CYP2B6 alter its enzymatic activity and substantially affect its pharmacokinetics. High-frequency variants, such as CYP2B6*6, are associated with the risk of developing side effects due to reduced CYP2B6 activity. However, the impact of rare alterations on enzyme function remains unknown, and some of these variants may significantly decrease the CYP2B6 activity. Therefore, in this study, we evaluated in vitro the functional alterations in 29 missense variants of the CYP2B6 gene identified in 8,380 Japanese individuals. Wild-type CYP2B6 and 29 rare CYP2B6 variants were transiently expressed in mammalian cells. The expression levels of variant CYP2B6 proteins in the microsomal fractions extracted from 293FT cells were assessed using western blotting and reduced-carbon monoxide difference spectroscopy, and a specific peak at 450 nm was detected in the wild-type and 19 variants. Furthermore, kinetic parameters were determined by assaying the reactions with efavirenz and propofol and quantifying the metabolite concentrations. We found that 12 variants had significantly lower or abolished enzymatic activity with both the substrates. In silico three-dimensional docking and molecular-dynamics simulations suggested that these functional changes were due to conformational changes in essential regions, such as the heme-binding site and ligand channels involved in transporting substrates to the active site. These findings have implications for predicting the plasma concentrations of CYP2B6 substrates and controlling their side effects.

Cytochrome P450 monooxygenase systems: Diversity and plasticity for adaptive stress response

Superfamily of cytochromes P450 (CYPs) is composed of heme-thiolate-containing monooxygenase enzymes, which play crucial roles in the biosynthesis, bioactivation, and detoxification of a variety of organic compounds, both endogenic and exogenic. Majority of CYP monooxygenase systems are multi-component and contain various redox partners, cofactors and auxiliary proteins, which contribute to their diversity in both prokaryotes and eukaryotes. Recent progress in bioinformatics and computational biology approaches make it possible to undertake whole-genome and phylogenetic analyses of CYPomes of a variety of organisms. Considerable variations in sequences within and between CYP families and high similarity in secondary and tertiary structures between all CYPs along with dramatic conformational changes in secondary structure elements of a substrate binding site during catalysis have been reported. This provides structural plasticity and substrate promiscuity, which underlie functional diversity of CYPs. Gene duplication and mutation events underlie CYP evolutionary diversity and emergence of novel selectable functions, which provide the involvement of CYPs in high adaptability to changing environmental conditions and dietary restrictions. In our review, we discuss the recent advancements and challenges in the elucidating the evolutionary origin and mechanisms underlying the CYP monooxygenase system diversity and plasticity. Our review is in the view of hypothesis that diversity of CYP monooxygenase systems is translated into the broad metabolic profiles, and this has been acquired during the long evolutionary time to provide structural plasticity leading to high adaptative capabilities to environmental stress conditions.

A Computational Pipeline Observes the Flexibility and Dynamics of Plant Cytochrome P450 Binding Sites

Binding site flexibility and dynamics strongly affect the ability of proteins to accommodate substrates and inhibitors. The significance of these properties is particularly pronounced for proteins that are inherently flexible, such as cytochrome P450 enzymes (CYPs). While the research on human CYPs provides detailed knowledge on both structural and functional level, such analyses are still lacking for their plant counterparts. This study aims to bridge this gap. We developed a novel computational pipeline consisting of two steps. Firstly, we use molecular dynamics (MD) simulations to capture the full conformational ensemble for a certain plant CYP. Subsequently, we developed and applied a comprehensive methodology to analyze a number of binding site properties-size, flexibility, shape, hydrophobicity, and accessibility-using the fpocket and mdpocket packages on MD-generated trajectories. The workflow was validated on human CYPs 1A2, 2A6, and 3A4, as their binding site characteristics are well known. Not only could we confirm known binding site properties, but we also identified and named previously unseen binding site channels for CYPs 1A2 and 2A6. The pipeline was then applied to plant CYPs, leading to the first categorization of 15 chosen plant CYPs based on their binding site's (dis)similarities. This study provides a foundation for the largely uncharted fields of plant CYP substrate specificity and facilitates a more precise understanding of their largely unknown specific biological functions. It offers new insights into the structural and functional dynamics of plant CYPs, which may facilitate a more accurate understanding of the fate of agrochemicals or the biotechnological design and exploitation of enzymes with specific functions. Additionally, it serves as a reference for future structural-functional analyses of CYP enzymes across various biological kingdoms.

Aromatase inhibitors for the treatment of breast cancer: An overview (2019-2023)

Aromatase inhibition is considered a legitimate approach for the treatment of ER-positive (ER+) breast cancer as it accounts for more than 70% of breast cancer cases. Aromatase inhibitor therapy has been demonstrated to be highly effective in decreasing tumour size, increasing survival rates, and lowering the chance of cancer recurrence. The present review deliberates the pathophysiology and the role of aromatase in estrogen biosynthesis. Estrogen biosynthesis, various androgens, and their function in the human body have also been discussed. The salient aspects of the aromatase active site, its mode of action, and AIs, along with their intended interactions with presently FDA-approved inhibitors, have been briefly discussed. It has been detailed how different reported AIs were designed, their SAR investigations, in silico analysis, and biological evaluations. Various AIs from multiple origins, such as synthetic and semi-synthetic, have also been discussed.

Multiresolution molecular dynamics simulations reveal the interplay between conformational variability and functional interactions in membrane-bound cytochrome P450 2B4

Cytochrome P450 2B4 (CYP 2B4) is one of the best-characterized CYPs and serves as a key model system for understanding the mechanisms of microsomal class II CYPs, which metabolize most known drugs. The highly flexible nature of CYP 2B4 is apparent from crystal structures that show the active site with either a wide open or a closed heme binding cavity. Here, we investigated the conformational ensemble of the full-length CYP 2B4 in a phospholipid bilayer, using multiresolution molecular dynamics (MD) simulations. Coarse-grained MD simulations revealed two predominant orientations of CYP 2B4's globular domain with respect to the bilayer. Their refinement by atomistic resolution MD showed adaptation of the enzyme's interaction with the lipid bilayer, leading to open configurations that facilitate ligand access to the heme binding cavity. CAVER analysis of enzyme tunnels, AquaDuct analysis of water routes, and Random Acceleration Molecular Dynamics simulations of ligand dissociation support the conformation-dependent passage of molecules between the active site and the protein surroundings. Furthermore, simulation of the re-entry of the inhibitor bifonazole into the open conformation of CYP 2B4 resulted in binding at a transient hydrophobic pocket within the active site cavity that may play a role in substrate binding or allosteric regulation. Together, these results show how the open conformation of CYP 2B4 facilitates the binding of substrates from and release of products to the membrane, whereas the closed conformation prolongs the residence time of substrates or inhibitors and selectively allows the passage of smaller reactants via the solvent and water channels.

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.

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.

Molecular Mechanisms Underlying Single Nucleotide Polymorphism-Induced Reactivity Decrease in CYP2D6

Cytochrome P450 2D6 (CYP2D6) is one of the most important enzymes involved in drug metabolism. Genetic polymorphism can influence drug metabolism by CYP2D6 such that a therapy is seriously affected by under- or overdosing of drugs. However, a general explanation at the atomistic level for poor activity is missing so far. Here we show for the 20 most common single nucleotide polymorphisms (SNPs) of CYP2D6 that poor metabolism is driven by four mechanisms. We found in extensive all-atom molecular dynamics simulations that the rigidity of the I-helix (central helix), distance between central phenylalanines (stabilizing bound substrate), availability of basic residues on the surface of CYP2D6 (binding of cytochrome P450 reductase), and position of arginine 132 (electron transfer to heme) are essential for an extensive function of the enzyme. These results were applied to SNPs with unknown effects, and potential SNPs that may lead to poor drug metabolism were identified. The revealed molecular mechanisms might be important for other drug-metabolizing cytochrome P450 enzymes.

Sequence-Guided Redesign of an Omega-Transaminase from Bacillus megaterium for the Asymmetric Synthesis of Chiral Amines

ω-Transaminases (ω-TAs) are attractive biocatalysts asymmetrically catalyzing ketones to chiral amines. However, poor non-native catalytic activity and substrate promiscuity severely hamper its wide application in industrial production. Protein engineering efforts have generally focused on reshaping the substrate-binding pockets of ω-TAs. However, hotspots around the substrate tunnel as well as distant sites outside the pockets may also affect its activity. In this study, the ω-TA from Bacillus megaterium (BmeTA) was selected for engineering. The tunnel mutation Y164F synergy with distant mutation A245T which was acquired through a multiple sequence alignment showed improved soluble expression, a 3.7-fold higher specific activity and a 19.9-fold longer half-life at 45 °C. Molecule Dynamics simulation explains the mechanism of improved catalytic activity, enhanced thermostability and improved soluble expression of BmeTAY164F/A245T(2 M). Finally, the resting cells of 2 M were used for biocatalytic processes. 450 mM of S-methoxyisopropylamine (S-MOIPA) was obtained with an ee value of 97.3 % and a conversion rate of 90 %, laying the foundation for its industrial production. Mutant 2 M was also found to be more advantageous in catalyzing the transamination of various ketones. These results demonstrated that sites that are far away from the active center also play an important role in the redesign of ω-TAs.

Studying the Collective Functional Response of a Receptor in Alchemical Ligand Binding Free Energy Simulations with Accelerated Solvation Layer Dynamics

Ligand binding free energy simulations (LB-FES) that involve sampling of protein functional conformations have been longstanding challenges in research on molecular recognition. Particularly, modeling of the conformational transition pathway and design of the heuristic biasing mechanism are severe bottlenecks for the existing enhanced configurational sampling (ECS) methods. Inspired by the key role of hydration in regulating conformational dynamics of macromolecules, this report proposes a novel ECS approach that facilitates binding-associated structural dynamics by accelerated hydration transitions in combination with the λ-exchange of free energy perturbation (FEP). Two challenging protein-ligand binding processes involving large configurational transitions of the receptor are studied, with hydration transitions at binding sites accelerated by Hamiltonian-simulated annealing of the hydration layer. Without the need for pathway analysis or ad hoc barrier flattening potential, LB-FES were performed with FEP/λ-exchange molecular dynamics simulation at a minor overhead for annealing of the hydration layer. The LB-FES studies showed that the accelerated rehydration significantly enhances the collective conformational transitions of the receptor, and convergence of binding affinity calculations is obtained at a sweet-spot simulation time scale. Alchemical LB-FES with the proposed ECS strategy is free from the effort of trial and error for the setup and realizes efficient on-the-fly sampling for the collective functional response of the receptor and bound water and therefore presents a practical approach to high-throughput screening in drug discovery.

Unraveling the Ligand-Binding Sites of CYP3A4 by Molecular Dynamics Simulations with Solvent Probes

Cytochrome P450 3A4 (CYP3A4) is one of the most important drug-metabolizing enzymes in the human body and is well known for its complicated, atypical kinetic characteristics. The existence of multiple ligand-binding sites in CYP3A4 has been widely recognized as being capable of interfering with the active pocket through allosteric effects. The identification of ligand-binding sites other than the canonical active site above the heme is especially important for understanding the atypical kinetic characteristics of CYP3A4 and the intriguing association between the ligand and the receptor. In this study, we first employed mixed-solvent molecular dynamics (MixMD) simulations coupled with the online computational predictive tools to explore potential ligand-binding sites in CYP3A4. The MixMD approach demonstrates better performance in dealing with the receptor flexibility compared with other computational tools. From the sites identified by MixMD, we then picked out multiple sites for further exploration using ensemble docking and conventional molecular dynamics (cMD) simulations. Our results indicate that three extra sites are suitable for ligand binding in CYP3A4, including one experimentally confirmed site and two novel sites.

MD investigation on the differences in the dynamic interactions between the specific ligand azamulin and two CYP3A isoforms, 3A4 and 3A5

The unmarked potential drug molecule azamulin has been found to be a specific inhibitor of CYP3A4 and CYP3A5 in recent years, but this molecule also shows different binding ability and affinity to the two CYP3A isoforms. In order to explore the microscopic mechanism, conventional molecular dynamics (MD) simulation methods were performed to study the dynamic interactions between two isoforms and azamulin. The simulation results show that the binding of the ligand leads to different structural properties of two CYP3A proteins. First of all, compared with apo-CYP3A4, the binding of the ligand azamulin can lead to changes in the structural flexibility of CYP3A4, i.e., holo-CYP3A4 is more flexible than apo-CYP3A4. The structural changes of CYP3A5 are just the opposite. The ligand binding increases the rigidity of CYP3A5. Furthermore, the representative structures of the production phase in the MD simulation were in details analyzed to obtain the microscopic interactions between the ligand azamulin and two CYP3A isoforms at the atomic level. It is speculated that the difference of composition and interaction of the active sites is the fundamental cause of the change of structural properties of the two proteins.Communicated by Ramaswamy H. Sarma.

Investigating the association between CYP2J2 inhibitors and QT prolongation: a literature review

Drug withdrawal post-marketing due to cardiotoxicity is a major concern for drug developers, regulatory agencies, and patients. One common mechanism of cardiotoxicity is through inhibition of cardiac ion channels, leading to prolongation of the QT interval and sometimes fatal arrythmias. Recently, oxylipin signaling compounds have been shown to bind to and alter ion channel function, and disruption in their cardiac levels may contribute to QT prolongation. Cytochrome P450 2J2 (CYP2J2) is the predominant CYP isoform expressed in cardiomyocytes, where it oxidizes arachidonic acid to cardioprotective epoxyeicosatrienoic acids (EETs). In addition to roles in vasodilation and angiogenesis, EETs bind to and activate various ion channels. CYP2J2 inhibition can lower EET levels and decrease their ability to preserve cardiac rhythm. In this review, we investigated the ability of known CYP inhibitors to cause QT prolongation using Certara's Drug Interaction Database. We discovered that among the multiple CYP isozymes, CYP2J2 inhibitors were more likely to also be QT-prolonging drugs (by approximately 2-fold). We explored potential binding interactions between these inhibitors and CYP2J2 using molecular docking and identified four amino acid residues (Phe61, Ala223, Asn231, and Leu402) predicted to interact with QT-prolonging drugs. The four residues are located near the opening of egress channel 2, highlighting the potential importance of this channel in CYP2J2 binding and inhibition. These findings suggest that if a drug inhibits CYP2J2 and interacts with one of these four residues, then it may have a higher risk of QT prolongation and more preclinical studies are warranted to assess cardiovascular safety.

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.

ChannelsDB 2.0: a comprehensive database of protein tunnels and pores in AlphaFold era

ChannelsDB 2.0 is an updated database providing structural information about the position, geometry and physicochemical properties of protein channels-tunnels and pores-within deposited biomacromolecular structures from PDB and AlphaFoldDB databases. The newly deposited information originated from several sources. Firstly, we included data calculated using a popular CAVER tool to complement the data obtained using original MOLE tool for detection and analysis of protein tunnels and pores. Secondly, we added tunnels starting from cofactors within the AlphaFill database to enlarge the scope of the database to protein models based on Uniprot. This has enlarged available channel annotations ∼4.6 times as of 1 September 2023. The database stores information about geometrical features, e.g. length and radius, and physico-chemical properties based on channel-lining amino acids. The stored data are interlinked with the available UniProt mutation annotation data. ChannelsDB 2.0 provides an excellent resource for deep analysis of the role of biomacromolecular tunnels and pores. The database is available free of charge: https://channelsdb2.biodata.ceitec.cz.

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.

Graphene quantum dots blocking the channel egresses of cytochrome P450 enzyme (CYP3A4) reveals potential toxicity

Graphene quantum dots (GQDs) have garnered significant attention, particularly in the biomedical domain. However, extensive research reveals a dichotomy concerning the potential toxicity of GQDs, presenting contrasting outcomes. Therefore, a comprehensive understanding of GQD biosafety necessitates a detailed supplementation of their toxicity profile. In this study, employing a molecular dynamics (MD) simulation approach, we systematically investigate the potential toxicity of GQDs on the CYP3A4 enzyme. We construct two distinct simulation systems, wherein a CYP3A4 protein is enveloped by either GQDs or GOQDs (graphene oxide quantum dots). Our results elucidate that GQDs come into direct contact with the bottleneck residues of Channels 2a and 2b of CYP3A4. Furthermore, GQDs entirely cover the exits of Channels 2a and 2b, implying a significant hindrance posed by GQDs to these channels and consequently leading to toxicity towards CYP3A4. In-depth analysis reveals that the adsorption of GQDs to the exits of Channels 2a and 2b is driven by a synergistic interplay of hydrophobic and van der Waals (vdW) interactions. In contrast, GOQDs only partially obstruct Channel 1 of CYP3A4, indicating a weaker influence on CYP3A4 compared to GQDs. Our findings underscore the potential deleterious impact of GQDs on the CYP3A4 enzyme, providing crucial molecular insights into GQD toxicology.

3site Multisubstrate-Bound State of Cytochrome P450cam

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

Hydration and Structural Adaptations of the Human CYP1A1, CYP1A2, and CYP1B1 Active Sites by Molecular Dynamics Simulations

Cytochromes CYP1A1, CYP1A2, and CYP1B1, the members of the cytochrome P450 family 1, catalyze the metabolism of endogenous compounds, drugs, and non-drug xenobiotics which include substances involved in the process of carcinogenesis, cancer chemoprevention, and therapy. In the present study, the interactions of three selected polymethoxy-trans-stilbenes, analogs of a bioactive polyphenol trans-resveratrol (3,5,4'-trihydroxy-trans-stilbene) with the binding sites of CYP1 isozymes were investigated with molecular dynamics (MD) simulations. The most pronounced structural changes in the CYP1 binding sites were observed in two substrate recognition sites (SRS): SRS2 (helix F) and SRS3 (helix G). MD simulations show that the number and position of water molecules occurring in CYP1 APO and in the structures complexed with ligands are diverse. The presence of water in binding sites results in the formation of water-protein, water-ligand, and bridging ligand-water-protein hydrogen bonds. Analysis of the solvent and substrate channels opening during the MD simulation showed significant differences between cytochromes in relation to the solvent channel and the substrate channels 2c, 2ac, and 2f. The results of this investigation lead to a deeper understanding of the molecular processes that occur in the CYP1 binding sites and may be useful for further molecular studies of CYP1 functions.

Pharmacokinetics during therapeutic hypothermia in neonates: from pathophysiology to translational knowledge and physiologically-based pharmacokinetic (PBPK) modeling

INTRODUCTION

Perinatal asphyxia (PA) still causes significant morbidity and mortality. Therapeutic hypothermia (TH) is the only effective therapy for neonates with moderate to severe hypoxic-ischemic encephalopathy after PA. These neonates need additional pharmacotherapy, and both PA and TH may impact physiology and, consequently, pharmacokinetics (PK) and pharmacodynamics (PD).

AREAS COVERED

This review provides an overview of the available knowledge in PubMed (until November 2022) on the pathophysiology of neonates with PA/TH. In vivo pig models for this setting enable distinguishing the effect of PA versus TH on PK and translating this effect to human neonates. Available asphyxia pig models and methodological considerations are described. A summary of human neonatal PK of supportive pharmacotherapy to improve neurodevelopmental outcomes is provided.

EXPERT OPINION

To support drug development for this population, knowledge from clinical observations (PK data, real-world data on physiology), preclinical (in vitro and in vivo (minipig)) data, and molecular and cellular biology insights can be integrated into a predictive physiologically-based PK (PBPK) framework, as illustrated by the I-PREDICT project (Innovative physiology-based pharmacokinetic model to predict drug exposure in neonates undergoing cooling therapy). Current knowledge, challenges, and expert opinion on the future directions of this research topic are provided.

Nanodisc-embedded cytochrome P450 P3A4 binds diverse ligands by distributing conformational dynamics to its flexible elements

Cytochrome P450 3A4 (CYP3A4) metabolizes a wide range of drugs and toxins. Interactions of CYP3A4 with ligands are difficult to predict due to promiscuity and conformational flexibility. To better understand CYP3A4 conformational responses to ligands we use hydrogen deuterium exchange mass spectrometry (HDX-MS) to investigate the effect of ligands on nanodisc-embedded CYP3A4. For a subset of CYP3A4-ligand complexes, differences in the low-frequency modes derived by principal component analyses of molecular dynamics trajectories mirrored the HDX-MS results. The effects of ligands are distributed to flexible elements of CYP3A4 between stretches of secondary structure. The largest effects occur in the F- and G-helices, where most ligands increase the flexibility of the F-helix and connecting loops and decrease the flexibility of the C-term of the G-helix. Most ligands affect the E-F-G, CD and HI regions of the protein. Ligand-dependent differences are observed in the A"-A' loop, BC region, E-helix, K-β1 region, proximal loop, and C-term loop. Correlated HDX responses were observed in the CD region and the C-term of the G-helix that were most pronounced for Type II ligands. Collectively, the HDX and molecular dynamics results suggest that CYP3A4 accommodates diverse binding partners by propagating local backbone fluctuations from the binding site onto the flexible regions of the enzyme via long-range interactions that are differentially modulated by ligands. In contrast to the paradigm wherein ligands decrease protein dynamics at their binding site, a wide range of ligands modestly increase CYP3A4 dynamics throughout the protein including effects remote from the active site.

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

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

Improved 2α-Hydroxylation Efficiency of Steroids by CYP154C2 Using Structure-Guided Rational Design

Cytochrome P450 enzymes are promising biocatalysts for industrial use because they catalyze site-selective C-H oxidation and have diverse catalytic reactions and a broad substrate range. In this study, the 2α-hydroxylation activity of CYP154C2 from Streptomyces avermitilis MA-4680T toward androstenedione (ASD) was identified by an in vitro conversion assay. The testosterone (TES)-bound structure of CYP154C2 was solved at 1.42 Å, and this structure was used to design eight mutants, including single, double, and triple mutants, to improve the conversion efficiency. Mutants L88F/M191F and M191F/V285L were found to enhance the conversion rates significantly (i.e., 8.9-fold and 7.4-fold for TES, 46.5-fold and 19.5-fold for ASD, respectively) compared with the wild-type (WT) enzyme while retaining high 2α-position selectivity. The substrate binding affinity of the L88F/M191F mutant toward TES and ASD was enhanced compared with that of WT CYP154C2, supporting the measured increase in the conversion efficiencies. Moreover, the total turnover number and kcat/Km of the L88F/M191F and M191F/V285L mutants increased significantly. Interestingly, all mutants containing L88F generated 16α-hydroxylation products, suggesting that L88 in CYP154C2 plays a vital role in substrate selectivity and that the amino acid corresponding to L88 in the 154C subfamily affects the orientation of steroid binding and substrate selectivity. IMPORTANCE Hydroxylated derivatives of steroids play essential roles in medicine. Cytochrome P450 enzymes selectively hydroxylate methyne groups on steroids, which can dramatically change their polarity, biological activity and toxicity. There is a paucity of reports on the 2α-hydroxylation of steroids, and documented 2α-hydroxylate P450s show extremely low conversion efficiency and/or low regio- and stereoselectivity. This study conducted crystal structure analysis and structure-guided rational engineering of CYP154C2 and efficiently enhanced the conversion efficiency of TES and ASD with high regio- and stereoselectivity. Our results provide an effective strategy and theoretical basis for the 2α-hydroxylation of steroids, and the structure-guided rational design of P450s should facilitate P450 applications in the biosynthesis of steroid drugs.

Multidisciplinary Insights into the Structure-Function Relationship of the CYP2B6 Active Site

Cytochrome P450 2B6 (CYP2B6) is a highly polymorphic human enzyme involved in the metabolism of many clinically relevant drugs, environmental toxins, and endogenous molecules with disparate structures. Over the last 20-plus years, in silico and in vitro studies of CYP2B6 using various ligands have provided foundational information regarding the substrate specificity and structure-function relationship of this enzyme. Approaches such as homology modeling, X-ray crystallography, molecular docking, and kinetic activity assays coupled with CYP2B6 mutagenesis have done much to characterize this originally neglected monooxygenase. However, a complete understanding of the structural details that make new chemical entities substrates of CYP2B6 is still lacking. Surprisingly little in vitro data has been obtained about the structure-function relationship of amino acids identified to be in the CYP2B6 active site. Since much attention has already been devoted to elucidating the function of CYP2B6 allelic variants, here we review the salient findings of in silico and in vitro studies of the CYP2B6 structure-function relationship with a deliberate focus on the active site. In addition to summarizing these complementary approaches to studying structure-function relationships, we note gaps/challenges in existing data such as the need for more CYP2B6 crystal structures, molecular docking results with various ligands, and data coupling CYP2B6 active site mutagenesis with kinetic parameter measurement under standard expression conditions. Harnessing in silico and in vitro techniques in tandem to understand the CYP2B6 structure-function relationship will likely offer further insights into CYP2B6-mediated metabolism. SIGNIFICANCE STATEMENT: The apparent importance of cytochrome P450 2B6 (CYP2B6) in the metabolism of various xenobiotics and endogenous molecules has grown since its discovery with many in silico and in vitro studies offering a partial description of its structure-function relationship. Determining the structure-function relationship of CYP2B6 is difficult but may be aided by thorough biochemical investigations of the CYP2B6 active site that provide a more complete pharmacological understanding of this important enzyme.

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.

Multiple-P450 Gene Co-Up-Regulation in the Development of Permethrin Resistance in the House Fly, Musca domestica

This paper reports a study conducted at the whole transcriptome level to characterize the P450 genes involved in the development of pyrethroid resistance, utilizing expression profile analyses of 86 cytochrome P450 genes in house fly strains with different levels of resistance to pyrethroids/permethrin. Interactions among the up-regulated P450 genes and possible regulatory factors in different autosomes were examined in house fly lines with different combinations of autosomes from a resistant house fly strain, ALHF. Eleven P450 genes that were significantly up-regulated, with levels > 2-fold those in the resistant ALHF house flies, were in CYP families 4 and 6 and located on autosomes 1, 3 and 5. The expression of these P450 genes was regulated by trans- and/or cis-acting factors, especially on autosomes 1 and 2. An in vivo functional study indicated that the up-regulated P450 genes also conferred permethrin resistance in Drosophila melanogaster transgenic lines. An in vitro functional study confirmed that the up-regulated P450 genes are able to metabolize not only cis- and trans-permethrin, but also two metabolites of permethrin, PBalc and PBald. In silico homology modeling and the molecular docking methodology further support the metabolic capacity of these P450s for permethrin and substrates. Taken together, the findings of this study highlight the important function of multi-up-regulated P450 genes in the development of insecticide resistance in house flies.

Computational method for aromatase-related proteins using machine learning approach

Human aromatase enzyme is a microsomal cytochrome P450 and catalyzes aromatization of androgens into estrogens during steroidogenesis. For breast cancer therapy, third-generation aromatase inhibitors (AIs) have proven to be effective; however patients acquire resistance to current AIs. Thus there is a need to predict aromatase-related proteins to develop efficacious AIs. A machine learning method was established to identify aromatase-related proteins using a five-fold cross validation technique. In this study, different SVM approach-based models were built using the following approaches like amino acid, dipeptide composition, hybrid and evolutionary profiles in the form of position-specific scoring matrix (PSSM); with maximum accuracy of 87.42%, 84.05%, 85.12%, and 92.02% respectively. Based on the primary sequence, the developed method is highly accurate to predict the aromatase-related proteins. Prediction scores graphs were developed using the known dataset to check the performance of the method. Based on the approach described above, a webserver for predicting aromatase-related proteins from primary sequence data was developed and implemented at https://bioinfo.imtech.res.in/servers/muthu/aromatase/home.html. We hope that the developed method will be useful for aromatase protein related research.

Identification and drug metabolic characterization of four new CYP2C9 variants CYP2C9*72-*75 in the Chinese Han population

Cytochrome 2C9 (CYP2C9), one of the most important drug metabolic enzymes in the human hepatic P450 superfamily, is required for the metabolism of 15% of clinical drugs. Similar to other CYP2C family members, CYP2C9 gene has a high genetic polymorphism which can cause significant racial and inter-individual differences in drug metabolic activity. To better understand the genetic distribution pattern of CYP2C9 in the Chinese Han population, 931 individuals were recruited and used for the genotyping in this study. As a result, seven synonymous and 14 non-synonymous variations were identified, of which 4 missense variants were designated as new alleles CYP2C9*72, *73, *74 and *75, resulting in the amino acid substitutions of A149V, R150C, Q214H and N418T, respectively. When expressed in insect cell microsomes, all four variants exhibited comparable protein expression levels to that of the wild-type CYP2C9 enzyme. However, drug metabolic activity analysis revealed that these variants exhibited significantly decreased catalytic activities toward three CYP2C9 specific probe drugs, as compared with that of the wild-type enzyme. These data indicate that the amino acid substitution in newly designated variants can cause reduced function of the enzyme and its clinical significance still needs further investigation in the future.

Unraveling the Abnormal Molecular Mechanism of Suicide Inhibition of Cytochrome P450 3A4

Suicide inhibition of the CYP3A4 enzyme by a drug inactivates the enzyme in the drug biotransformation process and often shows safety concerns about the drug. Despite extensive experimental studies, the abnormal molecular mechanism of a suicide inhibitor that forms a covalent bond with the residue far away from the catalytically active center of CYP3A4 inactivating the enzyme remains elusive. Here, the authors used molecular simulation approaches to study in detail how diquinone methide (DQR), the metabolite product of raloxifene, unbinds from CYP3A4 and inactivates the enzyme at the atomistic level. The results clearly indicate that in one of the intermediate states formed in its unbinding process, DQR covalently binds to Cys239, a residue far away from the catalytically active center of CYP3A4, and hinders the substrate from entering or leaving the enzyme. This work therefore provides an unprecedented way of clarifying the abnormal mechanism of suicide inhibition of the CYP3A4 enzyme.

Human Orphan Cytochrome P450 2U1 Catalyzes the ω-Hydroxylation of Leukotriene B4

Cytochrome P450 2U1 (CYP2U1) identified from the human genome remains poorly known since few data are presently available on its physiological function(s) and substrate(s) specificity. CYP2U1 mutations are associated with complicated forms of hereditary spastic paraplegia, alterations of mitochondrial architecture and bioenergetics. In order to better know the biological roles of CYP2U1, we used a bioinformatics approach. The analysis of the data invited us to focus on leukotriene B₄ (LTB₄), an important inflammatory mediator. Here, we show that CYP2U1 efficiently catalyzes the hydroxylation of LTB₄ predominantly on its ω-position. We also report docking experiments of LTB₄ in a 3D model of truncated CYP2U1 that are in agreement with this hydroxylation regioselectivity. The involvement of CYP2U1 in the metabolism of LTB₄ could have strong physiological consequences in cerebral pathologies including ischemic stroke because CYP2U1 is predominantly expressed in the brain.

Cholesterol's inhibition effect on entering of chlorzoxazone into phosphatidylethanolamine bilayer: Relevance to cytochrome P450 drug metabolism at endoplasmic reticulum membranes

Many drugs are metabolized by cytochrome P450 (CYP) in the endoplasmic reticulum (ER) membrane. Recent studies have shown that CYP-substrate drugs reach the CYP active site after entering the lipid hydrophobic part of the ER membrane. To clarify the role of cholesterol (Chol) in the CYP-related drug metabolic process, we investigated the lipid bilayer entry of CYP-substrate drugs using a model membrane system as follows. The model membrane system comprised palmitoyl-oleoyl-phosphatidylethanolamine (POPE) and Chol. Phosphatidylethanolamine is the second major phospholipid component of ER membranes. Chlorzoxazone (CZX) was used as the CYP-substrate drug. Calorimetric measurements showed that the addition of CZX to POPE bilayers decreased the gel-liquid crystal phase transition temperature; X-ray diffraction indicated that CZX distributes into the liquid crystal phase bilayers but not practically the gel phase POPE bilayers. In the presence of Chol, dialysis and X-ray structural analyses showed that Chol inhibited CZX entry into the bilayer with an increase in Chol concentration. The Chol concentration in the ER membrane (5-10 mol%) is much lower than that in the plasma membrane (approximately 30 mol%). This fact may allow CYP-substrate drugs to enter the hydrophobic portion of the ER membrane more easily than other organelle membranes, yielding efficient drug metabolism.

Site directed mutagenesis of Catharanthus roseus (+)-vincadifformine 19-hydroxylase (CYP71BY3) results in two distinct enzymatic functions

The most abundant monoterpenoid indole alkaloids (MIAs) in Catharanthus roseus roots include lochnericine and (+)-echitovenine. The formation of (+)-echitovenine involves a 3-step pathway including (+)-vincadifformine-19-hydroxylase (V19H) that differentiates it from a parallel pathway involved in the formation of lochnericine, hörhammericine and its O-acetylated derivative. Homology based modeling and docking experiments in the present study show that (+) and (-) vincadifformine can occupy the V19H active site and is proven experimentally by showing that (-)-vincadifformine is a competitive inhibitor of V19H. Comparative modeling of V19H with tabersonine 3-oxidase (T3O) and tabersonine 19-hydroxylase (T19H) that accept (-)-aspidosperma MIAs identified four conserved amino acid residues in T3O and T19H that were different in the V19H binding site and were used to generate a series of single-, double-, or four-point mutations in V19H. While all mutants retained their ability to convert (+)-vincadifformine to (+)-minovincinine only the four-point mutant gained T3O activity enabling it to convert (-)-tabersonine to tabersonine 2,3-epoxide. The gain of T3O-like activity following mutagenesis without the loss of V19H activity supports the hypothesis that V19H shares a common ancestor to T3O which is involved in vindoline biosynthesis in C. roseus leaves.

The Structures of the Steroid Binding CYP142 Cytochrome P450 Enzymes from Mycobacterium ulcerans and Mycobacterium marinum

The steroid binding CYP142 cytochrome P450 enzymes of Mycobacterium species are involved in the metabolism of cholesterol and its derivatives. The equivalent enzyme from Mycobacterium ulcerans was studied to compare the degree of functional conservation between members of this CYP family. We compared substrate binding of the CYP142A3 enzymes of M. ulcerans and M. marinum and CYP142A1 from M. tuberculosis using UV-vis spectroscopy. The catalytic oxidation of cholesterol derivatives by all three enzymes was undertaken. Both CYP142A3 enzymes were structurally characterized by X-ray crystallography. The amino acid sequences of the CYP142A3 enzymes are more similar to CYP142A1 from M. tuberculosis than CYP142A2 from Mycolicibacterium smegmatis. Both CYP142A3 enzymes have substrate binding properties, which are more resemblant to CYP142A1 than CYP142A2. The cholest-4-en-3-one-bound X-ray crystal structure of both CYP142A3 enzymes were determined at a resolution of <1.8 Å, revealing the substrate binding mode at a high level of detail. The structures of the cholest-4-en-3-one binding CYP142 enzymes from M. ulcerans and M. marinum demonstrate how the steroid binds in the active site of these enzymes. They provide an explanation for the high selectivity of the enzyme for terminal methyl C-H bond oxidation to form 26-hydroxy derivatives. These enzymes in pathogenic Mycobacterium species are candidates for inhibition. The work here demonstrates that similar drug molecules could target these CYP142 enzymes from different species in order to combat Buruli ulcer or tuberculosis.

"Commandeuring" Xenobiotic Metabolism: Advances in Understanding Xenobiotic Metabolism

The understanding of how exogenous chemicals (xenobiotics) are metabolized, distributed, and eliminated is critical to determine the impact of the chemical and its metabolites to the (human) organism. This is part of the research and educational discipline ADMET (absorption, distribution, metabolism, elimination, and toxicity). Here, we review the work of Jan Commandeur and colleagues who have not only made a significant impact in understanding of phase I and phase II metabolism of several important compounds but also contributed greatly to the development of experimental techniques for the study of xenobiotic metabolism. Jan Commandeur’s work has covered a broad area of research, such as the development of online screening methodologies, the use of a combination of enzyme mutagenesis and molecular modeling for structure–activity relationship (SAR) studies, and the development of novel probe substrates. This work is the bedrock of current activities and brings the field closer to personalized (cohort-based) pharmacology, toxicology, and hazard/risk assessment.

Chemical manipulation of mitochondrial function affects metabolism of red carotenoids in a marine copepod (Tigriopus californicus)

The shared-pathway hypothesis offers a cellular explanation for the connection between ketocarotenoid pigmentation and individual quality. Under this hypothesis, ketocarotenoid metabolism shares cellular pathways with mitochondrial oxidative phosphorylation such that red carotenoid-based coloration is inextricably linked mitochondrial function. To test this hypothesis, we exposed Tigriopus californicus copepods to a mitochondrially targeted protonophore, 2,4-dinitrophenol (DNP), to induce proton leak in the inner mitochondrial membranes. We then measured whole-animal metabolic rate and ketocarotenoid accumulation. As observed in prior studies of vertebrates, we observed that DNP treatment of copepods significantly increased respiration and that DNP-treated copepods accumulated more ketocarotenoid than control animals. Moreover, we observed a relationship between ketocarotenoid concentration and metabolic rate, and this association was strongest in DNP-treated copepods. These data support the hypothesis that ketocarotenoid and mitochondrial metabolism are biochemically intertwined. Moreover, these results corroborate observations in vertebrates, perhaps suggesting a fundamental connection between ketocarotenoid pigmentation and mitochondrial function that should be explored further.

The involvement of CYP1A2 in biodegradation of dioxins in pigs

2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD) is one of the most harmful chemicals showing resistance to biodegradation. The majority of TCDD effects is mediated by the aryl hydrocarbon receptor (AhR) pathway. TCDD binding to AhR results in the activation of cytochrome P450 enzymes (CYP1A1, CYP1A2, CYP1B1) involved in dioxin biodegradation. The goal of the study was to explore the potentialrole of CYP1A2 in the metabolism of TCDD. We investigated a molecular structure of CYP1A2 and the binding selectivity and affinity between the pig CYP1A2 and: 1/ DiCDD or TCDD (dioxins differing in toxicity and biodegradability) or 2/ their selected metabolites. pCYP1A2 demonstrated higher affinity towards DiCDD and TCDD than other pCYP1 enzymes. All dioxin-pCYP1A2 complexes were found to be stabilized by hydrophobic interactions. The calculated distances between the heme oxygen and the dioxin carbon nearest to the oxygen, reflecting the hydroxylating potential of CYP1A2, were higher than in other pCYP1 enzymes. The distances between the heme iron and the nearest dioxin carbon exceeded 5 Å, a distance sufficient to allow the metabolites to leave the active site. However, the molecular dynamics simulations revealed that two access channels of CYP1A2 were closed upon binding the majority of the examined dioxins. Moreover, the binding of dioxin metabolites did not promote opening of channel S–an exit for hydroxylated products. It appears that the undesired changes in the behavior of access channels prevail over the hydroxylating potential of CYP1A2 towards TCDD and the favorable distances, ultimately trapping the metabolites at the enzyme’s active site.

Molecular Basis of the Recognition of Cholesterol by Cytochrome P450 46A1 along the Major Access Tunnel

CYP46A1 is an important potential target for the treatment of Alzheimer's disease (AD), which is the most common neurodegenerative disease among older individuals. However, the binding mechanism between CYP46A1 and substrate cholesterol (CH) has not been clarified and will not be conducive to the research of relevant drug molecules. In this study, we integrated molecular docking, molecular dynamics (MD) simulations, and adaptive steered MD simulations to explore the recognition and binding mechanism of CYP46A1 with CH. Two key factors affecting the interaction between CH and CYP46A1 are determined: one is a hydrophobic cavity formed by seven hydrophobic residues (F80, Y109, L112, I222, W368, F371, and T475), which provides nonpolar interactions to stabilize CH, and the other is a hydrogen bond formed by H81 and CH, which ensures the binding direction of CH. In addition, the tunnel analysis results show that tunnel 2a is identified as the primary pathway of CH. The entry of CH induces tunnel 2e to close and tunnel w to open. Our results may provide effective clues for the design of drugs based on the substrate for AD and improve our understanding of the structure-function of CYP46A1.

Modulating the Coupling Efficiency of P450 BM3 by Controlling Water Diffusion through Access Tunnel Engineering

Cytochromes P450 have gained much interest for their broad substrate scope in the catalysis of oxidation reactions for pharmaceuticals, plastics, and hormones. However, achieving high coupling efficiency by the engineering of P450s is still a big challenge. The presence of extra water around the active site is deemed to be related to uncoupling. In this study, the access tunnels of P450 BM3 from Bacillus megaterium are engineered to control water access from bulk solvent to the active site. Nine residues located in tunnels are investigated by site-saturation mutagenesis to reduce water diffusion, thereby improving the coupling efficiency. The recombined variant N319L/T411V/T436A shows improved coupling efficiency (from 31.2 % to 52.6 %). Tunnel polarity analysis and molecular dynamics simulation further indicate that reduced water molecules around the active site lead to higher coupling efficiency. Overall, this study provides valuable insight on improving coupling efficiency by controlling water diffusion through tunnel engineering.

Determination of Multidirectional Pathways for Ligand Release from the Receptor: A New Approach Based on Differential Evolution

Steered molecular dynamics (SMD) simulation is a powerful method in computer-aided drug design as it can be used to access the relative binding affinity with high precision but with low computational cost. The success of SMD depends on the choice of the direction along which the ligand is pulled from the receptor-binding site. In most simulations, the unidirectional pathway was used, but in some cases, this choice resulted in the ligand colliding with the complex surface of the exit tunnel. To overcome this difficulty, several variants of SMD with multidirectional pulling have been proposed, but they are not completely devoid of disadvantages. Here, we have proposed to determine the direction of pulling with a simple scoring function that minimizes the receptor–ligand interaction, and an optimization algorithm called differential evolution is used for energy minimization. The effectiveness of our protocol was demonstrated by finding expulsion pathways of Huperzine A and camphor from the binding site of Torpedo California acetylcholinesterase and P450cam proteins, respectively, and comparing them with the previous results obtained using memetic sampling and random acceleration molecular dynamics. In addition, by applying this protocol to a set of ligands bound with LSD1 (lysine specific demethylase 1), we obtained a much higher correlation between the work of pulling force and experimental data on the inhibition constant IC50 compared to that obtained using the unidirectional approach based on minimal steric hindrance.

Structural characterization of the homotropic cooperative binding of azamulin to human cytochrome P450 3A5

Cytochrome P450 3A4 and 3A5 catalyze the metabolic clearance of a large portion of therapeutic drugs. Azamulin is used as a selective inhibitor for 3A4 and 3A5 to define their roles in metabolism of new chemical entities during drug development. In contrast to 3A4, 3A5 exhibits homotropic cooperativity for the sequential binding of two azamulin molecules at concentrations used for inhibition. To define the underlying sites and mechanisms for cooperativity, an X-ray crystal structure of 3A5 was determined with two azamulin molecules in the active site that are stacked in an antiparallel orientation. One azamulin resides proximal to the heme in a pose similar to the 3A4-azamulin complex. Comparison to the 3A5 apo structure indicates that the distal azamulin in 3A5 ternary complex causes a significant induced fit that excludes water from the hydrophobic surfaces of binding cavity and the distal azamulin, which is augmented by the stacking interaction with the proximal azamulin. Homotropic cooperativity was not observed for the binding of related pleuromutilin antibiotics, tiamulin, retapamulin, and lefamulin, to 3A5, which are larger and unlikely to bind in the distal site in a stacked orientation. Formation of the 3A5 complex with two azamulin molecules may prevent time-dependent inhibition that is seen for 3A4 by restricting alternate product formation and/or access of reactive intermediates to vulnerable protein sites. These results also contribute to a better understanding of sites for cooperative binding and the differential structural plasticity of 3A5 and 3A4 that contribute to differential substrate and inhibitor binding.

X-ray Crystal Structures of Methane Monooxygenase Hydroxylase Complexes with Variants of Its Regulatory Component: Correlations with Altered Reaction Cycle Dynamics

Full activity of soluble methane monooxygenase (sMMO) depends upon the formation of a 1:1 complex of the regulatory protein MMOB with each alpha subunit of the (αβγ)₂ hydroxylase, sMMOH. Previous studies have shown that mutations in the core region of MMOB and in the N- and C-termini cause dramatic changes in the rate constants for steps in the sMMOH reaction cycle. Here, X-ray crystal structures are reported for the sMMOH complex with two double variants within the core region of MMOB, DBL1 (N107G/S110A), and DBL2 (S109A/T111A), as well as two variants in the MMOB N-terminal region, H33A and H5A. DBL1 causes a 150-fold decrease in the formation rate constant of the reaction cycle intermediate P, whereas DBL2 accelerates the reaction of the dinuclear Fe(IV) intermediate Q with substrates larger than methane by three- to fourfold. H33A also greatly slows P formation, while H5A modestly slows both formation of Q and its reactions with substrates. Complexation with DBL1 or H33A alters the position of sMMOH residue R245, which is part of a conserved hydrogen-bonding network encompassing the active site diiron cluster where P is formed. Accordingly, electron paramagnetic resonance spectra of sMMOH:DBL1 and sMMOH:H33A complexes differ markedly from that of sMMOH:MMOB, showing an altered electronic environment. In the sMMOH:DBL2 complex, the position of M247 in sMMOH is altered such that it enlarges a molecular tunnel associated with substrate entry into the active site. The H5A variant causes only subtle structural changes despite its kinetic effects, emphasizing the precise alignment of sMMOH and MMOB required for efficient catalysis.

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.

Relaxed Substrate Requirements of Sterol 14α-Demethylase from Naegleria fowleri Are Accompanied by Resistance to Inhibition

Naegleria fowleri is the protozoan pathogen that causes primary amoebic meningoencephalitis (PAM), with the death rate exceeding 97%. The amoeba makes sterols and can be targeted by sterol biosynthesis inhibitors. Here, we characterized N. fowleri sterol 14-demethylase, including catalytic properties and inhibition by clinical antifungal drugs and experimental substituted azoles with favorable pharmacokinetics and low toxicity. None of them inhibited the enzyme stoichiometrically. The highest potencies were displayed by posaconazole (IC50 = 0.69 μM) and two of our compounds (IC50 = 1.3 and 0.35 μM). Because both these compounds penetrate the brain with concentrations reaching minimal inhibitory concentration (MIC) values in an N. fowleri cellular assay, we report them as potential drug candidates for PAM. The 2.1 Å crystal structure, in complex with the strongest inhibitor, provides an explanation connecting the enzyme weaker substrate specificity with lower sensitivity to inhibition. It also provides insight into the enzyme/ligand molecular recognition process and suggests directions for the design of more potent inhibitors.

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.