Cytochrome P450 structure–function: insights from molecular dynamics simulations

The HNF1α-Regulated LncRNA HNF1α-AS1 Is Involved in the Regulation of Cytochrome P450 Expression in Human Liver Tissues and Huh7 Cells

Expression of cytochrome P450s (P450s) is regulated by epigenetic factors, such as DNA methylation, histone modifications, and noncoding RNAs through different mechanisms. Among these factors, long noncoding RNAs (lncRNAs) have been shown to play important roles in the regulation of gene expression; however, little is known about the effects of lncRNAs on the regulation of P450 expression. The aim of this study was to explore the role of lncRNAs in the regulation of P450 expression by using human liver tissues and hepatoma Huh7 cells. Through lncRNA microarray analysis and quantitative polymerase chain reaction in human liver tissues, we found that the lncRNA hepatocyte nuclear factor 1 alpha antisense 1 (HNF1α-AS1), an antisense RNA of HNF1α, is positively correlated with the mRNA expression of CYP2C8, 2C9, 2C19, 2D6, 2E1, and 3A4 as well as pregnane X receptor (PXR) and constitutive androstane receptor (CAR). Gain- and loss-of-function studies in Huh7 cells transfected with small interfering RNAs or overexpression plasmids showed that HNF1α not only regulated the expression of HNF1α-AS1 and P450s, but also regulated the expression of CAR, PXR, and aryl hydrocarbon receptor (AhR). In turn, HNF1α-AS1 regulated the expression of PXR and most P450s without affecting the expression of HNF1α, AhR, and CAR. Moreover, the rifampicin-induced expression of P450s was also affected by HNF1α and HNF1α-AS1. In summary, the results of this study suggested that HNF1α-AS1 is involved in the HNF1α-mediated regulation of P450s in the liver at both basal and drug-induced levels.

Genetic Polymorphisms and In Silico Mutagenesis Analyses of CYP2C9, CYP2D6, and CYPOR Genes in the Pakistani Population

Diverse distributions of pharmacogenetically relevant variants of highly polymorphic CYP2C9, CYP2D6 and CYPOR genes are responsible for some varied drug responses observed across human populations. There is limited data available regarding the pharmacogenetic polymorphisms and frequency distributions of major allele variants in the Pakistani population. The present in silico mutagenesis study conducted on genotype pharmacogenetic variants and comparative analysis with a global population aims to extend the currently limited pharmacogenetic available evidence for the indigenous Pakistani population. Extracted genomic DNA from 244 healthy individuals' venous blood samples were amplified for distinct variant loci in the CYP2C9, CYP2D6 and CYPOR genes. Two-way sequencing results were compared with standard PubMed data and sequence variant loci confirmed by Chromas. This study revealed significant variations in CYP2C9 (rs1799853, rs1057910 and rs72558189), CYP2D6 (rs16947 and rs1135840), and CYPOR (rs1057868, rs781919285 and rs562750402) variants in intraethnic and interethnic frequency distributions. In silico mutagenesis and three-dimensional protein structural alignment analysis approaches clearly exposed the possible varied impact of rare CYPOR (rs781919285 and rs562750402) single nucleotide polymorphisms (SNPs) and confirmed that the influences of CYP2C9 and CYP2D6 variants are consistent with what was found in earlier studies. This investigation highlighted the need to study pharmacogenetic relevance loci and documentation since evidence could be utilized to elucidate genetic backgrounds of drug metabolism, and provide a basis for future pharmacogenomic studies and adequate dose adjustments in Pakistani and global populations.

Substrate mediated redox partner selectivity of cytochrome P450

Investigating the interplay between cytochrome-P450 and its redox partners (CPR and cytochrome-b5) is vital for understanding the metabolism of most hydrophobic drugs. Dynamic structural interactions with the ternary complex, with and without substrates, captured by NMR reveal a gating mechanism for redox partners to promote P450 function.

Evaluation of modafinil as a perpetrator of metabolic drug-drug interactions using a model informed cocktail reaction phenotyping trial protocol

AIM

To evaluate the capacity for modafinil to be a perpetrator of metabolic drug-drug interactions by altering cytochrome P450 activity following a single dose and dosing to steady state.

METHODS

A single centre, open label, single sequence cocktail drug interaction trial. On days 0, 2 and 8 participants were administered an oral drug cocktail comprising 100 mg caffeine, 30 mg dextromethorphan, 25 mg losartan, 1 mg midazolam and 20 mg enteric-coated omeprazole. Timed blood samples were collected prior to and for up to 6 h post cocktail dosing. Between days 2 and 8 participants orally self-administered 200 mg modafinil each morning.

RESULTS

Following a single 200 mg dose of modafinil mean (± 95% CI) AUC ratios for caffeine, dextromethorphan, losartan, midazolam and omeprazole were 0.95 (± 0.08), 1.01 (± 0.35), 0.97 (± 0.10), 0.98 (± 0.10) and 1.36 (± 0.06), respectively. Following dosing of modafinil to steady state (200 mg for 7 days), AUC ratios for caffeine, dextromethorphan, losartan, midazolam and omeprazole were 0.90 (± 0.16), 0.79 (± 0.09), 0.98 (± 0.11), 0.66 (± 0.12) and 1.90 (± 0.53), respectively.

CONCLUSIONS

These data support consideration of the risk of clinically relevant metabolic drug-drug interactions perpetrated by modafinil when this drug is co-administered with drugs that are primarily cleared by CYP2C19 (single modafinil dose or steady state modafinil dosing) or CYP3A4 (steady state modafinil dosing only) catalysed metabolic pathways.

Cytochrome P450-dependent reactive oxygen species (ROS) production contributes to Mn3O4 nanoparticle-caused liver injury

Mn₃O₄ nanoparticles (NPs) are one of the most important nanomaterials, and have a wide range of applications (i.e., catalysis, solar-electron transformation and molecular adsorption).

Pharmacogenomics of CYP2C9: Functional and Clinical Considerations

CYP2C9 is the most abundant CYP2C subfamily enzyme in human liver and the most important contributor from this subfamily to drug metabolism. Polymorphisms resulting in decreased enzyme activity are common in the CYP2C9 gene and this, combined with narrow therapeutic indices for several key drug substrates, results in some important issues relating to drug safety and efficacy. CYP2C9 substrate selectivity is detailed and, based on crystal structures for the enzyme, we describe how CYP2C9 catalyzes these reactions. Factors relevant to clinical response to CYP2C9 substrates including inhibition, induction and genetic polymorphism are discussed in detail. In particular, we consider the issue of ethnic variation in pattern and frequency of genetic polymorphisms and clinical implications. Warfarin is the most well studied CYP2C9 substrate; recent work on use of dosing algorithms that include CYP2C9 genotype to improve patient safety during initiation of warfarin dosing are reviewed and prospects for their clinical implementation considered. Finally, we discuss a novel approach to cataloging the functional capabilities of rare 'variants of uncertain significance', which are increasingly detected as more exome and genome sequencing of diverse populations is conducted.

Thermodynamics of camphor migration in cytochrome P450cam by atomistic simulations

Understanding the mechanisms of ligand binding to enzymes is of paramount importance for the design of new drugs. Here, we report on the use of a novel biased molecular dynamics (MD) methodology to study the mechanism of camphor binding to cytochrome P450cam. Microsecond-long MD simulations allowed us to observe reaction coordinates characterizing ligand diffusion from the active site of cytochrome P450cam to solvent via three egress routes. These atomistic simulations were used to estimate thermodynamic quantities along the reaction coordinates and indicate diverse binding configurations. The results suggest that the diffusion of camphor along the pathway near the substrate recognition site (SRS) is thermodynamically preferred. In addition, we show that the diffusion near the SRS is triggered by a transition from a heterogeneous collection of closed ligand-bound conformers to the basin comprising the open conformations of cytochrome P450cam. The conformational change accompanying this switch is characterized by the retraction of the F and G helices and the disorder of the B' helix. These results are corroborated by experimental studies and provide detailed insight into ligand binding and conformational behavior of the cytochrome family. The presented methodology is general and can be applied to other ligand-protein systems.

Single or Multiple Access Channels to the CYP450s Active Site? An Answer from Free Energy Simulations of the Human Aromatase Enzyme

Cytochromes P450 (CYP450s), in particular, CYP19A1 and CYP17A1, are key clinical targets of breast and prostate anticancer therapies, critical players in drug metabolism, and their overexpression in tumors is associated with drug resistance. In these enzymes, ligand (substrates, drugs) metabolism occurs in deeply buried active sites accessible only via several grueling channels, whose exact biological role remains unclear. Gaining direct insights on the mechanism by which ligands travel in and out is becoming increasingly important given that channels are involved in the modulation of binding/dissociation kinetics and the specificity of ligands toward a CYP450. This has profound implications for enzymatic efficiency and drug efficacy/toxicity. Here, by applying free energy methods, for a cumulative simulation time of 20 μs, we provide detailed atomistic characterization and free energy profiles of the entry/exit routes preferentially followed by a substrate (androstenedione) and a last-generation inhibitor (letrozole) to/from the catalytic site of CYP19A1 (the human aromatase (HA) enzyme), a key clinical target against breast cancer, studied here as prototypical CYP450. Despite the remarkably different size/shape/hydrophobicity of the ligands, two channels appear accessible to their entrance, while only one exit route appears to be preferential. Our study shows that the preferential paths may be conserved among different CYP450s. Moreover, our results highlight that, at least in the case of HA, ligand channeling is associated with large enzyme structural rearrangements. A wise choice of the computational method and very long simulations are, thus, required to obtain fully converged quantitative free energy profiles, which might be used to design novel biocatalysts or next-generation cytochrome inhibitors with an in silico tuned K_m.

Predicting Regioselectivity and Lability of Cytochrome P450 Metabolism Using Quantum Mechanical Simulations

We describe methods for predicting cytochrome P450 (CYP) metabolism incorporating both pathway-specific reactivity and isoform-specific accessibility considerations. Semiempirical quantum mechanical (QM) simulations, parametrized using experimental data and ab initio calculations, estimate the reactivity of each potential site of metabolism (SOM) in the context of the whole molecule. Ligand-based models, trained using high-quality regioselectivity data, correct for orientation and steric effects of the different CYP isoform binding pockets. The resulting models identify a SOM in the top 2 predictions for between 82% and 91% of compounds in independent test sets across seven CYP isoforms. In addition to predicting the relative proportion of metabolite formation at each site, these methods estimate the activation energy at each site, from which additional information can be derived regarding their lability in absolute terms. We illustrate how this can guide the design of compounds to overcome issues with rapid CYP metabolism.