Membrane-embedded substrate recognition by cytochrome P450 3A4

Small angle scattering reveals the orientation of cytochrome P450 19A1 in lipoprotein nanodiscs

Human aromatase (CYP19A1), the cytochrome P450 enzyme responsible for conversion of androgens to estrogens, was incorporated into lipoprotein nanodiscs (NDs) and interrogated by small angle X-ray and neutron scattering (SAXS/SANS). CYP19A1 was associated with the surface and centered at the edge of the long axis of the ND membrane. In the absence of the N-terminal anchor, the amphipathic A'- and G'-helices were predominately buried in the lipid head groups, with the possibly that their hydrophobic side chains protrude into the hydrophobic, aliphatic tails. The prediction is like that for CYP3A4 based on SAXS employing a similar modeling approach. The orientation of CYP19A1 in a ND is consistent with our previous predictions based on molecular dynamics simulations and lends additional credibility to the notion that CYP19A1 captures substrates from the membrane.

[Metabolic Activities Catalyzed by Human Cytochrome P450 (CYP) 2D6 and CYP3A Subfamily Members and Effect of Various Compounds, Including Endogenous Steroid Hormones, on These Activities]

My research focused on the effects of various drugs on (1) dopamine formation from p-tyramine catalyzed by polymorphic cytochrome P450 (CYP or P450) 2D6 variants and (2) endogenous steroid hormone hydroxylation catalyzed by CYP3A subfamily members (CYP3A4, CYP3A5, CYP3A7). The activation (cooperativity) of metabolic reactions catalyzed by P450s was especially emphasized. The effects of various psychotropic agents on dopamine formation from p-tyramine, catalyzed by wild-type CYP2D6.1 and CYP2D6 variants, including CYP2D6.2 (Arg296Cys;Ser486Thr), CYP2D6.10 (Pro34Ser;Ser486Thr), and CYP2D6.39 (Ser486Thr) were compared. Michaelis (Km) and inhibition (Ki) constants of the psychotropic agents in the presence of CYP2D6.10 were higher than those observed in the presence of other CYP2D6 variants. Fluvoxamine, fluoxetine, milnacipran, and haloperidol activated CYP2D6-catalyzed dopamine formation [decreasing the Km and/or increasing the maximal velocity (kcat)], and this activation was CYP2D6 variant-dependent. Regarding the CYP3A subfamily, the effects of various compounds including endogenous steroid hormones on the 6β-hydroxylation of steroid hormones, such as testosterone, progesterone, and cortisol, were determined; it was found that testosterone, dehydroepiandrosterone, and/or α-naphthoflavone activated 6β-hydroxylation of cortisol and/or progesterone, but the effects varied in the presence of different CYP3A subfamily members. Further studies are required to confirm the mechanisms and therapeutic relevance of these activation phenomena.

Differential Effects of Clotrimazole on X-Ray Crystal Structures of Human Cytochromes P450 3A5 and 3A4

Cytochromes P450 CYP3A5 and CYP3A4 exhibit differential plasticity that underlies differences in drug metabolism and drug-drug interactions. To extend previous studies, CYP3A4 and CYP3A5 were cocrystallized with clotrimazole, a compact ligand that binds to the heme iron in the catalytic center of the active site. Binding studies indicate that clotrimazole exhibits tight binding to CYP3A5 with a binding affinity (Kd) of <0.01 μM like that of CYP3A4. A single clotrimazole is bound to the heme iron in CYP3A4 that triggers expansion of active site cavity that reflects a loss of aromatic interactions between phenylalanine sidechains in the distal active site and increased conformational entropy for the F-F' connector due to reorientation of Phe-304 to accommodate clotrimazole. In contrast to CYP3A4, the CYP3A5 Phe-304 exhibits an induced fit along with Phe-213 to form edge-to-face aromatic interactions with heme-bound clotrimazole. These aromatic interactions between aromatic amino acids propagate by induced fits with a second clotrimazole residing in the distal active site and a third clotrimazole bound in an expanded entrance channel as well as between the three clotrimazoles. The large, expanded entrance channel surrounded by the C-terminal loop and the F' and A' helices in CYP3A5 suggests conformational selection for the binding of clotrimazole due to its large girth, which may also cause the entrance channel to remain open after the binding of the first clotrimazole to the heme iron. The additional binding sites suggest a path for sequential binding of one molecule to reach and bind to the heme iron. SIGNIFICANCE STATEMENT: Clotrimazole binds to the heme iron of CYP3A5 and CYP3A4. In CYP3A5, two clotrimazoles also bind in the distal active site and in an expanded entrance channel. Aromatic interactions between clotrimazoles and phenylalanine sidechains including Phe-304 indicate induced fits for each clotrimazole. In contrast to CYP3A5, displacement of the CYP3A4 Phe-304 rotamer by clotrimazole leads to extensive disruption of phenylalanine interactions that limit the space above the heme, to an expanded active site cavity, and to increased CYP3A4 conformational heterogeneity.

Induced Fit Describes Ligand Binding to Membrane-Associated Cytochrome P450 3A4

Cytochrome P450 3A4 (CYP3A4) is the dominant P450 involved in human xenobiotic metabolism. Competition for CYP3A4 therefore underlies several adverse drug-drug interactions. Despite its clinical significance, the mechanisms CYP3A4 uses to bind diverse ligands remain poorly understood. Highly monodisperse CYP3A4 embedded in anionic lipoprotein nanodiscs containing an equal mixture of 1-palmitoyl-2-oleoyl-glycero-3-phosphocholine and 1-palmitoyl-2-oleoyl-sn-glycero-3-phospho-(1'-rac-glycerol) were used to determine which of the limiting kinetic schemes that include protein conformational change, conformational selection (CS) or induced fit (IF), best described the binding of four known irreversible inhibitors. Azamulin, retapamulin, pleuromutilin, and mibrefadil binding to CYP3A4 nanodiscs conformed to a single-site binding model. Exponential fits of stopped-flow UV-visible absorption spectroscopy data supported multiple-step binding mechanisms. Trends in the rates of relaxation to equilibrium with increasing ligand concentrations were ambiguous as to whether IF or CS was involved; however, global fitting and consideration of the rate constants favored an IF mechanism. In the case of mibrefadil, a transient complex was observed in the stopped-flow UV-visible experiment, definitively assigning the presence of IF in ligand binding. While these studies only consider a small region of CYP3A4's vast ligand space, they provide kinetic evidence that CYP3A4 can use an IF mechanism. SIGNIFICANCE STATEMENT: CYP3A4 is capable of oxidizing numerous xenobiotics, including many drugs. Such promiscuity could not be achieved without conformational changes to accommodate diverse substrates. It is unknown whether conformational heterogeneity is present before (conformational selection) or after (induced fit) ligand binding. Stopped-flow measurements of suicide inhibitors binding to nanodisc-embedded CYP3A4 combined with sophisticated numerical analyses support that induced fit better describes ligand binding to this important enzyme.

Interaction of CYP3A4 with caffeine: First insights into multiple substrate binding

Human cytochrome P450 3A4 (CYP3A4) is a major drug-metabolizing enzyme that shows extreme substrate promiscuity. Moreover, its large and malleable active site can simultaneously accommodate several substrate molecules of the same or different nature, which may lead to cooperative binding and allosteric behavior. Due to difficulty of crystallization of CYP3A4-substrate complexes, it remains unknown how multiple substrates can arrange in the active site. We determined crystal structures of CYP3A4 bound to three and six molecules of caffeine, a psychoactive alkaloid serving as a substrate and modulator of CYP3A4. In the ternary complex, one caffeine binds to the active site suitably for C8-hydroxylation, most preferable for CYP3A4. In the senary complex, three caffeine molecules stack parallel to the heme with the proximal ligand poised for 3-N-demethylation. However, the caffeine stack forms extensive hydrophobic interactions that could preclude product dissociation and multiple turnovers. In both complexes, caffeine is also bound in the substrate channel and on the outer surface known as a peripheral site. At all sites, aromatic stacking with the caffeine ring(s) is likely a dominant interaction, while direct and water-mediated polar contacts provide additional stabilization for the substrate-bound complexes. Protein-ligand interactions via the active site R212, intrachannel T224, and peripheral F219 were experimentally confirmed, and the latter two residues were identified as important for caffeine association. Collectively, the structural, spectral, and mutagenesis data provide valuable insights on the ligand binding mechanism and help better understand how purine-based pharmaceuticals and other aromatic compounds could interact with CYP3A4 and mediate drug-drug interactions.

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.

Crystal Structure of CYP3A4 Complexed with Fluorol Identifies the Substrate Access Channel as a High-Affinity Ligand Binding Site

Cytochrome P450 3A4 (CYP3A4) is a major human drug-metabolizing enzyme, notoriously known for its extreme substrate promiscuity, allosteric behavior, and implications in drug–drug interactions. Despite extensive investigations, the mechanism of ligand binding to CYP3A4 is not fully understood. We determined the crystal structure of CYP3A4 complexed with fluorol, a small fluorescent dye that can undergo hydroxylation. In the structure, fluorol associates to the substrate channel, well suited for the binding of planar polyaromatic molecules bearing polar groups, through which stabilizing H-bonds with the polar channel residues, such as Thr224 and Arg372, can be established. Mutagenesis, spectral, kinetic, and functional data confirmed the involvement but not strict requirement of Thr224 for the association of fluorol. Collectively, our data identify the substrate channel as a high-affinity ligand binding site and support the notion that hydrophobic ligands first dock to the nearby peripheral surface, before migrating to the channel and, subsequently, into the active site.

Effect of Human Cytochrome P450 2D6 Polymorphism on Progesterone Hydroxylation

BACKGROUND AND OBJECTIVES

Herein, hydroxylation activities at the 6β-position and 21-position of progesterone mediated by human cytochrome P450 (CYP) 2D6 and its variants and the effects of psychotropic drugs on these hydroxylation activities were compared to clarify whether CYP2D6 polymorphisms and psychotropic drugs impact neurosteroid levels in the brain.

METHODS

Progesterone was incubated with CYP2D6.1, CYP2D6.2 (Arg296Cys, Ser486Thr), CYP2D6.10 (Pro34Ser, Ser486Thr), and CYP2D6.39 (Ser486Thr) in the absence or presence of typical psychotropic drugs (fluvoxamine, fluoxetine, paroxetine, fluphenazine, and milnacipran) and endogenous steroids (testosterone and cortisol). Then, 6β- and 21-hydroxyprogesterone levels were determined by high-performance liquid chromatography.

RESULTS

Although the Michaelis-Menten constants (Km) for progesterone 6β- and 21-hydroxylation reactions mediated by the different CYP2D6 variants were similar, the maximal velocity (Vmax) values of the reactions mediated by CYP2D6.1 and CYP2D6.2 were the highest, followed by those mediated by CYP2D6.39 and CYP2D6.10. Thus, the of progesterone 6β- and/or 21-hydroxylation reactions mediated by CYP2D6.1 and CYP2D6.2 showed the highest Vmax/Km values, followed by the reactions mediated by CYP2D6.39. All investigated compounds inhibited progesterone 21-hydroxylation mediated by CYP2D6 variants at high concentrations. Interestingly, at low concentrations, fluoxetine increased progesterone 21-hydroxylation mediated by CYP2D6.1, but not that mediated by CYP2D6.2 or CYP2D6.10. In addition, the Km value for CYP2D6.2 was elevated in the presence of fluoxetine, whereas the value for CYP2D6.1 was unaltered; however, Vmax values of both CYP2D6.1 and CYP2D6.2 were increased. Paroxetine competitively inhibited CYP2D6.1- and CYP2D6.2-mediated progesterone 21-hydroxylation.

CONCLUSIONS

These results suggest that CYP2D6 polymorphism can affect the stimulation/inhibition of progesterone 21-hydroxylation.

Silibinin Suppresses the Hyperlipidemic Effects of the ALK-Tyrosine Kinase Inhibitor Lorlatinib in Hepatic Cells

The third-generation anaplastic lymphoma tyrosine kinase inhibitor (ALK-TKI) lorlatinib has a unique side effect profile that includes hypercholesteremia and hypertriglyceridemia in >80% of lung cancer patients. Here, we tested the hypothesis that lorlatinib might directly promote the accumulation of cholesterol and/or triglycerides in human hepatic cells. We investigated the capacity of the hepatoprotectant silibinin to modify the lipid-modifying activity of lorlatinib. To predict clinically relevant drug−drug interactions if silibinin were used to clinically manage lorlatinib-induced hyperlipidemic effects in hepatic cells, we also explored the capacity of silibinin to interact with and block CYP3A4 activity using in silico computational descriptions and in vitro biochemical assays. A semi-targeted ultrahigh pressure liquid chromatography accurate mass quadrupole time-of-flight mass spectrometry with electrospray ionization (UHPLC-ESI-QTOF-MS/MS)-based lipidomic approach revealed that short-term treatment of hepatic cells with lorlatinib promotes the accumulation of numerous molecular species of cholesteryl esters and triglycerides. Silibinin treatment significantly protected the steady-state lipidome of hepatocytes against the hyperlipidemic actions of lorlatinib. Lipid staining confirmed the ability of lorlatinib to promote neutral lipid overload in hepatocytes upon long-term exposure, which was prevented by co-treatment with silibinin. Computational analyses and cell-free biochemical assays predicted a weak to moderate inhibitory activity of clinically relevant concentrations of silibinin against CYP3A4 when compared with recommended (rosuvastatin) and non-recommended (simvastatin) statins for lorlatinib-associated dyslipidemia. The elevated plasma cholesterol and triglyceride levels in lorlatinib-treated lung cancer patients might involve primary alterations in the hepatic accumulation of lipid intermediates. Silibinin could be clinically explored to reduce the undesirable hyperlipidemic activity of lorlatinib in lung cancer patients.

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.

Investigation of the molecular and mechanistic basis for the regioselective metabolism of midazolam by cytochrome P450 3A4

Cytochrome P450 3A4 (CYP3A4) is the most important P450 enzyme for drug metabolism and drug-drug interaction, due to it being responsible for the biotransformation of approximately 50% of clinically used drugs. Advance knowledge of the molecular and mechanistic basis of CYP3A4 regioselective metabolism is beneficial for understanding the production of metabolites, and may allow personalized metabolic pathways or designing pathway-specific therapeutics. In this work, we focus on investigating the ligand-receptor interactions, substrate conformational transition, and key factors regulating the specificity of metabolic pathways using midazolam (MDZ) as a probe. Here, three types of substrate-binding conformations related to the diversity of MDZ metabolites are identified. The results also suggest that an allosteric site for MDZ is located near the F'-helix, A-anchor, and C-terminal loop of CYP3A4. The presence of an effector in the allosteric site can accelerate the conformational transition of the substrate via modulating a "sandwich" structure, and may affect the proportion of metabolites at high substrate concentration. We hope that the results can improve the understanding of the CYP3A4 structure and function, and provide a new perspective for drug development.

Preparative Production of Functionalized (N- and O-Heterocyclic) Polycyclic Aromatic Hydrocarbons by Human Cytochrome P450 3A4 in a Bioreactor

Polycyclic aromatic hydrocarbons (PAHs) and their N- and O-containing derivatives (N-/O-PAHs) are environmental pollutants and synthetically attractive building blocks in pharmaceuticals. Functionalization of PAHs can be achieved via C-H activation by cytochrome P450 enzymes (e.g., P450 CYP3A4) in an environmentally friendly manner. Despite its broad substrate scope, the contribution of CYP3A4 to metabolize common PAHs in humans was found to be small. We recently showcased the potential of CYP3A4 in whole-cell biocatalysis with recombinant yeast Komagataella phaffii (Pichia pastoris) catalysts for the preparative-scale synthesis of naturally occurring metabolites in humans. In this study, we aimed at exploring the substrate scope of CYP3A4 towards (N-/O)-PAHs and conducted a bioconversion experiment at 10 L scale to validate the synthetic potential of CYP3A4 for the preparative-scale production of functionalized PAH metabolites. Hydroxylated products were purified and characterized using HPLC and NMR analysis. In total, 237 mg of fluorenol and 48 mg of fluorenone were produced from 498 mg of fluorene, with peak productivities of 27.7 μmol/L/h for fluorenol and 5.9 μmol/L/h for fluorenone; the latter confirmed that CYP3A4 is an excellent whole-cell biocatalyst for producing authentic human metabolites.

Mechanistic Understanding from Molecular Dynamics in Pharmaceutical Research 2: Lipid Membrane in Drug Design

We review the use of molecular dynamics (MD) simulation as a drug design tool in the context of the role that the lipid membrane can play in drug action, i.e., the interaction between candidate drug molecules and lipid membranes. In the standard "lock and key" paradigm, only the interaction between the drug and a specific active site of a specific protein is considered; the environment in which the drug acts is, from a biophysical perspective, far more complex than this. The possible mechanisms though which a drug can be designed to tinker with physiological processes are significantly broader than merely fitting to a single active site of a single protein. In this paper, we focus on the role of the lipid membrane, arguably the most important element outside the proteins themselves, as a case study. We discuss work that has been carried out, using MD simulation, concerning the transfection of drugs through membranes that act as biological barriers in the path of the drugs, the behavior of drug molecules within membranes, how their collective behavior can affect the structure and properties of the membrane and, finally, the role lipid membranes, to which the vast majority of drug target proteins are associated, can play in mediating the interaction between drug and target protein. This review paper is the second in a two-part series covering MD simulation as a tool in pharmaceutical research; both are designed as pedagogical review papers aimed at both pharmaceutical scientists interested in exploring how the tool of MD simulation can be applied to their research and computational scientists interested in exploring the possibility of a pharmaceutical context for their research.

Role of amino acids at positions 34, 296, and 486 of cytochrome P450 2D6 in the stimulatory and inhibitory effects of psychotropic agents on dopamine formation from p-tyramine

The effects of psychotropic agents such as fluvoxamine, fluoxetine, paroxetine, milnacipran, and fluphenazine on dopamine formation from p-tyramine catalysed by cytochrome P450 (CYP) 2D6.2 (Arg296Cys;Ser486Thr), CYP2D6.10 (Pro34Ser;Ser486Thr), and CYP2D6.39 (Ser486Thr) were compared with the effects on dopamine formation from p-tyramine by CYP2D6.1. Michaelis constants (K_m) and maximal velocity (k_cat) values for dopamine formation and inhibition constants (K_i) of the psychotropic agents were determined.For CYP2D6.39, the k_cat values for fluvoxamine, fluoxetine, and milnacipran, but not for paroxetine and fluphenazine, gradually increased with increasing concentrations, indicating activation of the catalysed reaction.Fluphenazine competitively inhibited dopamine formation catalysed by all variants, with a higher K_i value for CYP2D6.10. Among the three compounds that have a trifluoromethyl group in their chemical structure, only fluvoxamine and fluoxetine, as well as milnacipran that does not have this group, decreased K_m values and/or increased k_cat values for dopamine formation, suggesting that the group may not be essential for the activation.These findings indicate that substitution of amino acids at positions 34 and 486 can affect the affinity (K_m) and enzymatic activity (k_cat), respectively, for milnacipran and that the effect of substitution of arginine to cysteine at the 296^th position on the activation would be effector dependent.

Structural Dynamics of Cytochrome P450 3A4 in the Presence of Substrates and Cytochrome P450 Reductase

Cytochrome P450 3A4 (CYP3A4) is the most important drug-metabolizing enzyme in humans and has been associated with harmful drug interactions. The activity of CYP3A4 is known to be modulated by several compounds and by the electron transfer partner, cytochrome P450 reductase (CPR). The underlying mechanism of these effects, however, is poorly understood. We have used hydrogen-deuterium exchange mass spectrometry to investigate the impact of binding of CPR and of three different substrates (7-benzyloxy-4-trifluoromethyl-coumarin, testosterone, and progesterone) on the conformational dynamics of CYP3A4. Here, we report that interaction of CYP3A4 with substrates or with the oxidized or reduced forms of CPR leads to a global rigidification of the CYP3A4 structure. This was evident from the suppression of deuterium exchange in several regions of CYP3A4, including regions known to be involved in protein-protein interactions (helix C) and substrate binding and specificity (helices B' and E, and loop K/β1). Furthermore, the bimodal isotopic distributions observed for some CYP3A4-derived peptides were drastically impacted upon binding to CPR and/or substrates, suggesting the existence of stable CYP3A4 conformational populations that are perturbed by ligand/CPR binding. The results have implications for understanding the mechanisms of ligand binding, allostery, and catalysis in CYP enzymes.

Comparison of the Stimulatory and Inhibitory Effects of Steroid Hormones and α-Naphthoflavone on Steroid Hormone Hydroxylation Catalyzed by Human Cytochrome P450 3A Subfamilies

The inhibitory and stimulatory effects of steroid hormones and related compounds on the hydroxylation activity at the 6β-position of two steroid hormones, progesterone and testosterone, by CYP3A4, polymorphically expressed CYP3A5, and fetal CYP3A7 were compared to clarify the catalytic properties of the predominant forms of the human CYP3A subfamily. Hydroxylation activities of progesterone and testosterone by CYP3A4, CYP3A5, and CYP3A7 were estimated using HPLC. The Michaelis constants (K_m) for progesterone 6β-hydroxylation by CYP3A5 were markedly decreased in the presence of dehydroepiandrosterone (DHEA) and α-naphthoflavone (ANF), whereas progesterone and DHEA competitively inhibited testosterone 6β-hydroxylation mediated by CYP3A4, and progesterone competitively inhibited CYP3A5-mediated activity, which was weaker than that for CYP3A4. ANF noncompetitively inhibited testosterone 6β-hydroxylation mediated by both CYP3A4 and CYP3A5. Progesterone and testosterone 6β-hydroxylation mediated by CYP3A7 was inhibited or unaffected by DHEA, pregnenolone, and ANF. These results suggested that DHEA and ANF stimulated progesterone 6β-hydroxylation by CYP3A5 but not by CYP3A4 and CYP3A7; however, progesterone, DHEA, and ANF inhibited testosterone 6β-hydroxylation mediated by all CYP3A subfamily members. The inhibitory/stimulatory pattern of steroid-steroid interactions is different among CYP3A subfamily members and CYP3A5 is the most sensitive in terms of activation among the CYP3A subfamily members investigated.

Comparison of the inhibitory effects of azole antifungals on cytochrome P450 3A4 genetic variants

Cytochrome P450 (CYP) 3A4 is one of the major drug-metabolizing enzymes. Genetic variants of CYP3A4 with altered activity are one of the factors responsible for interindividual differences in drug metabolism. Azole antifungals inhibit CYP3A4 to cause clinically significant drug-drug interactions. In the present quantitative study, we investigated the inhibitory effects of three azole antifungals (ketoconazole, voriconazole, and fluconazole) on testosterone metabolism by recombinant CYP3A4 genetic variants (CYP3A4.1 (WT), CYP3A4.2, CYP3A4.7, CYP3A4.16, and CYP3A4.18) and compared them with those previously reported for itraconazole. The inhibition constants (K_i) of ketoconazole, voriconazole, and fluconazole for rCYP3A4.1 were 3.6 nM, 3.2 μM, and 16.1 μM, respectively. The K_i values of these azoles for rCYP3A4.16 were 13.9-, 13.6-, and 6.2-fold higher than those for rCYP3A4.1, respectively, whereas the K_i value of itraconazole for rCYP3A4.16 was 0.54-fold of that for rCYP3A4.1. The other genetic variants had similar effects on the K_i values of the three azoles, whereas a very different pattern was seen for itraconazole. In conclusion, itraconazole has unique characteristics that are distinct from those shared by the other azole anti-fungal drugs ketoconazole, voriconazole, and fluconazole with regard to the influence of genetic variations on the inhibition of CYP3A4.

The Effect of Force-Field Parameters on Cytochrome P450-Membrane Interactions: Structure and Dynamics

The simulation of membrane proteins requires compatible protein and lipid force fields that reproduce the properties of both the protein and the lipid bilayer. Cytochrome P450 enzymes are bitopic membrane proteins with a transmembrane helical anchor and a large cytosolic globular domain that dips into the membrane. As such, they are representative and challenging examples of membrane proteins for simulations, displaying features of both peripheral and integral membrane proteins. We performed molecular dynamics simulations of three cytochrome P450 isoforms (2C9, 2C19 and 1A1) in a 2-oleoyl-1-palmitoyl-sn-glycerol-3-phosphocholine bilayer using two AMBER force field combinations: GAFF-LIPID with ff99SB for the protein, and LIPID14 with ff14SB for the protein. Comparison of the structural and dynamic properties of the proteins, the lipids and the protein-membrane interactions shows differing sensitivity of the cytochrome P450 isoforms to the choice of force field, with generally better agreement with experiment for the LIPID14 + ff14SB combination.

Comparison of steroid hormone hydroxylation mediated by cytochrome P450 3A subfamilies

Hydroxylation activity at the 6β-position of steroid hormones (testosterone, progesterone, and cortisol) by human cytochromes P450 (CYP) 3A4, polymorphic CYP3A5, and fetal CYP3A7 were compared to understand the catalytic properties of the major forms of human CYP3A subfamily. Testosterone, progesterone, and cortisol 6β-hydroxylation activities of recombinant CYP3A4, CYP3A5, and CYP3A7 were determined by liquid chromatography. Michaelis constants (K_m) for CYP3A7-mediated 6β-hydroxylation of testosterone, progesterone, and cortisol were similar to those of CYP3A4 and CYP3A5. The maximal velocity (k_cat) and k_cat/K_m values for CYP3A4 were the highest, followed by CYP3A5 and those for CYP3A7 were the lowest among three CYP3A subfamily members. A decrease in K_m values for progesterone 6β-hydroxylation by CYP3A4, CYP3A5, and CYP3A7 in the presence of testosterone was observed, and the k_cat values for CYP3A5 gradually increased with increasing testosterone. This indicated that testosterone stimulated progesterone 6β-hydroxylation by all three CYP3A subfamily members. However, progesterone inhibited testosterone 6β-hydroxylation mediated by CYP3A4, CYP3A5, and CYP3A7. In conclusion, the k_cat values, rather than K_m values, for 6β-hydroxylation of three steroid hormones mediated by CYP3A7 were different from those for CYP3A4 and CYP3A5. In addition, the inhibitory/stimulatory pattern of steroid-steroid interactions would be different among CYP3A subfamily members.

Dynamics and Location of the Allosteric Midazolam Site in Cytochrome P4503A4 in Lipid Nanodiscs

Promiscuous and allosteric drug interactions with cytochrome P450 3A4 (CYP3A4) are ubiquitous but incompletely understood at the molecular level. A classic allosteric CYP3A4 drug interaction includes the benzodiazepine midazolam (MDZ). MDZ exhibits homotropic and heterotropic allostery when metabolized to 1'-hydroxy and 4-hydroxy metabolites in varying ratios. The combination of hydrogen-deuterium exchange mass spectrometry (HDX-MS) and Gaussian accelerated molecular dynamics (GaMD) simulations of CYP3A4 in lipid nanodiscs and in a lipid bilayer, respectively, reveals MDZ-dependent changes in dynamics in a membrane environment. The F-, G-, and intervening helices, as well as the loop preceding the β1-sheets, display the largest observed changes in HDX. The GaMD suggests a potential allosteric binding site for MDZ in the F'- and G'-regions, which undergo significant increases in HDX at near-saturating MDZ concentrations. The HDX-MS and GaMD results confirm that changes in dynamics are most significant near the developing consensus allosteric site, and these changes are distinct from those observed previously with the nonallosteric inhibitor ketoconazole. The results suggest that the allosteric MDZ remains mobile in its binding site at the Phe-cluster. The results further suggest that this binding site remains dynamic or changes the depth of insertion in the membrane.

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.

Differing Membrane Interactions of Two Highly Similar Drug-Metabolizing Cytochrome P450 Isoforms: CYP 2C9 and CYP 2C19

The human cytochrome P450 (CYP) 2C9 and 2C19 enzymes are two highly similar isoforms with key roles in drug metabolism. They are anchored to the endoplasmic reticulum membrane by their N-terminal transmembrane helix and interactions of their cytoplasmic globular domain with the membrane. However, their crystal structures were determined after N-terminal truncation and mutating residues in the globular domain that contact the membrane. Therefore, the CYP-membrane interactions are not structurally well-characterized and their dynamics and the influence of membrane interactions on CYP function are not well understood. We describe herein the modeling and simulation of CYP 2C9 and CYP 2C19 in a phospholipid bilayer. The simulations revealed that, despite high sequence conservation, the small sequence and structural differences between the two isoforms altered the interactions and orientations of the CYPs in the membrane bilayer. We identified residues (including K72, P73, and I99 in CYP 2C9 and E72, R73, and H99 in CYP 2C19) at the protein-membrane interface that contribute not only to the differing orientations adopted by the two isoforms in the membrane, but also to their differing substrate specificities by affecting the substrate access tunnels. Our findings provide a mechanistic interpretation of experimentally observed effects of mutagenesis on substrate selectivity.

Multiscale Simulations of Biological Membranes: The Challenge To Understand Biological Phenomena in a Living Substance

Biological membranes are tricky to investigate. They are complex in terms of molecular composition and structure, functional over a wide range of time scales, and characterized by nonequilibrium conditions. Because of all of these features, simulations are a great technique to study biomembrane behavior. A significant part of the functional processes in biological membranes takes place at the molecular level; thus computer simulations are the method of choice to explore how their properties emerge from specific molecular features and how the interplay among the numerous molecules gives rise to function over spatial and time scales larger than the molecular ones. In this review, we focus on this broad theme. We discuss the current state-of-the-art of biomembrane simulations that, until now, have largely focused on a rather narrow picture of the complexity of the membranes. Given this, we also discuss the challenges that we should unravel in the foreseeable future. Numerous features such as the actin-cytoskeleton network, the glycocalyx network, and nonequilibrium transport under ATP-driven conditions have so far received very little attention; however, the potential of simulations to solve them would be exceptionally high. A major milestone for this research would be that one day we could say that computer simulations genuinely research biological membranes, not just lipid bilayers.

Ilexgenin A enhances the effects of simvastatin on non-alcoholic fatty liver disease without changes in simvastatin pharmacokinetics

Cardiovascular disease (CVD) is the most common cause of death in patients with non-alcoholic fatty liver disease (NAFLD). New therapeutic strategies which have the potential for slowing down the evolution of NAFLD and reducing CVD-related mortality are urgently needed. Statins are well recognized in the treatment of dyslipidemia, but their use in the treatment of NAFLD is limited due to the safety concerns. Ilexgenin A (IA) is one of the main bioactive compounds in 'Shan-lv-cha', an herbal tea commonly used in China. In the present study, we investigated the possible synergistic therapeutic effects of IA and simvastatin (SV) on NAFLD. IA or SV showed beneficial effects on the rats with NAFLD by lowering the liver weight, liver index and plasma levels of alanine aminotransferase and aspartate aminotransferase, regulating abnormal metabolism of lipids and ameliorating steatosis in liver. IA significantly enhanced the hypolipidemic and anti-inflammation effects of SV. Furthermore, a sensitive, accurate, convenient and reproducible LC-MS method was developed to investigate the effects of IA on the pharmacokinetics of SV. No significant changes were observed in pharmacokinetic parameters of SV and simvastatin hydroxy acid in the IA plus SV co-treated group in comparison with those in the group treated with SV alone. The mRNA levels and activity of CYP3A1 were not altered by IA. In conclusion, the results obtained from the present study should be helpful for further clinical application of SV and IA alone or in combination.