π-π Stacking mediated drug-drug interactions in human CYP2E1

Automated Analysis of Soft Matter Interfaces, Interactions, and Self-Assembly with PySoftK

Molecular dynamics simulations have become essential tools in the study of soft matter and biological macromolecules. The large amount of high-dimensional data associated with such simulations does not straightforwardly elucidate the atomistic mechanisms that underlie complex materials and molecular processes. Analysis of these simulations is complicated: the dynamics intrinsic to soft matter simulations necessitates careful application of specific, and often complex, algorithms to extract meaningful molecular scale understanding. There is an ongoing need for high-quality automated computational workflows to facilitate this analysis in a reproducible manner with minimal user input. In this work, we introduce a series of molecular simulation analysis tools for investigating interfaces, molecular interactions (including ring-ring stacking), and self-assembly. In addition, we include a number of auxiliary tools, including a useful function to unwrap molecular structures that are greater than half the length of their corresponding simulation box. These tools are contained in the PySoftK software package, making the application of these algorithms straightforward for the user. These new simulation analysis tools within PySoftK will support high-quality, reproducible analysis of soft matter and biomolecular simulations to bring about new predictive understanding in nano- and biotechnology.

Inhibitory effects of the main metabolites of Apatinib on CYP450 isozymes in human and rat liver microsomes

PURPOSE

The inhibitory effect of Apatinib on cytochrome P450 (CYP450) enzymes has been studied. However, it is unknown whether the inhibition is related to the major metabolites, M1-1, M1-2 and M1-6.

METHODS

A 5-in-1 cocktail system composed of CYP2B6/Cyp2b1, CYP2C9/Cyp2c11, CYP2E1/Cyp2e1, CYP2D6/Cyp2d1 and CYP3A/Cyp3a2 was used in this study. Firstly, the effects of APA and its main metabolites on the activities of HLMs, RLMs and recombinant isoforms were examined. The reaction mixture included HLMs, RLMs or recombinant isoforms (CYP3A4.1, CYP2D6.1, CYP2D6.10 or CYP2C9.1), analyte (APA, M1-1, M1-2 or M1-6), probe substrates. The reactions were pre-incubated for 5 min at 37 °C, followed by the addition of NAPDH to initiate the reactions, which continued for 40 min. Secondly, IC50 experiments were conducted to determine if the inhibitions were reversible. The reaction mixture of the "+ NADPH Group" included HLMs or RLMs, 0 to 100 of μM M1-1 or M1-2, probe substrates. The reactions were pre-incubated for 5 min at 37 °C, and then NAPDH was added to initiate reactions, which proceeded for 40 min. The reaction mixture of the "- NADPH Group" included HLMs or RLMs, probe substrates, NAPDH. The reactions were pre-incubated for 30 min at 37 °C, and then 0 to 100 μM of M1-1 or M1-2 was added to initiate the reactions, which proceeded for 40 min. Finally, the reversible inhibition of M1-1 and M1-2 on isozymes was determined. The reaction mixture included HLMs or RLMs, 0 to 10 μM of M1-1 or M1-2, probe substrates with concentrations ranging from 0.25Km to 2Km.

RESULTS

Under the influence of M1-6, the activity of CYP2B6, 2C9, 2E1 and 3A4/5 was increased to 193.92%, 210.82%, 235.67% and 380.12% respectively; the activity of CYP2D6 was reduced to 92.61%. The inhibitory effects of M1-1 on CYP3A4/5 in HLMs and on Cyp2d1 in RLMs, as well as the effect of M1-2 on CYP3A in HLMs, were determined to be noncompetitive inhibition, with the Ki values equal to 1.340 μM, 1.151 μM and 1.829 μM, respectively. The inhibitory effect of M1-1 on CYP2B6 and CYP2D6 in HLMs, as well as the effect of M1-2 on CYP2C9 and CYP2D6 in HLMs, were determined to be competitive inhibition, with the Ki values equal to 12.280 μM, 2.046 μM, 0.560 μM and 4.377 μM, respectively. The inhibitory effects of M1-1 on CYP2C9 in HLMs and M1-2 on Cyp2d1 in RLMs were determined to be mixed-type, with the Ki values equal to 0.998 μM and 0.884 μM. The parameters could not be obtained due to the atypical kinetics of CYP2E1 in HLMs under the impact of M1-2.

CONCLUSIONS

M1-1 and M1-2 exhibited inhibition for several CYP450 isozymes, especially CYP2B6, 2C9, 2D6 and 3A4/5. This observation may uncover potential drug-drug interactions and provide valuable insights for the clinical application of APA.

Phenylalanine Residues in the Active Site of CYP2E1 Participate in Determining the Binding Orientation and Metabolism-Dependent Genotoxicity of Aromatic Compounds

The composition of amino acids forming the active site of a CYP enzyme is impactful in its substrate selectivity. For CYP2E1, the role of PHE residues in the formation of effective binding orientations for its aromatic substrates remains unclear. In this study, molecular docking and molecular dynamics analysis were performed to reflect the interactions between PHEs in the active site of human CYP2E1 and various aromatic compounds known as its substrates. The results indicated that the orientation of 1-methylpyrene (1-MP) in the active site was highly determined by the presence of PHEs, PHE478 contributing to the binding free energy most significantly. Moreover, by building a random forest model the relationship between each of 19 molecular descriptors of polychlorinated biphenyl (PCB) compounds (from molecular docking, quantum mechanics, and physicochemical properties) and their human CYP2E1-dependent mutagenicityas established mostly in our lab, was investigated. The presence of PHEs did not appear to significantly modify the electronic or structural feature of each bound ligand (PCB), instead, the flexibility of the conformation of PHEs contributed substantially to the effective binding energy and orientation. It is supposed that PHE residues adjust their own conformation to permit a suitablly shaped cavity for holding the ligand and forming its orientation as favorable for a biochemical reaction. This study has provided some insights into the role of PHEs in guiding the interactive adaptation of the active site of human CYP2E1 for the binding and metabolism of aromatic substrates.

Effects of resveratrol on drug- and carcinogen-metabolizing enzymes, implications for cancer prevention

Resveratrol is a polyphenol found in grape skins and peanuts that has demonstrated many health benefits including protection against aging, cardiovascular and metabolic disease, neurological decline, and cancer. The anticancer properties of resveratrol have been attributed to a variety of mechanisms, including its general inhibition of phase I metabolism and induction of phase II metabolism. The effects of resveratrol on these enzymes, however, are still unclear, as in vitro evidence often contrasts with animal studies and clinical trials. Reasons for these variances could include the low bioavailability of resveratrol and the effects of resveratrol metabolites. Due to resveratrol's interactions with drug-metabolizing enzymes and drug transporters, individuals concurrently taking pharmacological doses of resveratrol with other supplements or medications could potentially experience nutrient-drug interactions. This review summarizes the known effects of resveratrol and its main metabolites on drug metabolism in order to help characterize which populations might benefit from resveratrol for the prevention of cancer, as well as those that may need to avoid supplementation due to potential drug interactions.

Pi-pi Stacking Mediated Cooperative Mechanism for Human Cytochrome P450 3A4

Human Cytochrome P450 3A4 (CYP3A4) is an important member of the cytochrome P450 superfamily with responsibility for metabolizing ~50% of clinical drugs. Experimental evidence showed that CYP3A4 can adopt multiple substrates in its active site to form a cooperative binding model, accelerating substrate metabolism efficiency. In the current study, we constructed both normal and cooperative binding models of human CYP3A4 with antifungal drug ketoconazoles (KLN). Molecular dynamics simulation and free energy calculation were then carried out to study the cooperative binding mechanism. Our simulation showed that the second KLN in the cooperative binding model had a positive impact on the first one binding in the active site by two significant pi-pi stacking interactions. The first one was formed by Phe215, functioning to position the first KLN in a favorable orientation in the active site for further metabolism reactions. The second one was contributed by Phe304. This pi-pi stacking was enhanced in the cooperative binding model by the parallel conformation between the aromatic rings in Phe304 and the dioxolan moiety of the first KLN. These findings can provide an atomic insight into the cooperative binding in CYP3A4, revealing a novel pi-pi stacking mechanism for drug-drug interactions.

Structural basis for cooperative binding of azoles to CYP2E1 as interpreted through guided molecular dynamics simulations

CYP2E1 metabolizes a wide array of small, hydrophobic molecules, resulting in their detoxification or activation into carcinogens through Michaelis-Menten as well as cooperative mechanisms. Nevertheless, the molecular determinants for CYP2E1 specificity and metabolic efficiency toward these compounds are still unknown. Herein, we employed computational docking studies coupled to molecular dynamics simulations to provide a critical perspective for understanding a structural basis for cooperativity observed for an array of azoles from our previous binding and catalytic studies (Hartman et al., 2014). The resulting 28 CYP2E1 complexes in this study revealed a common passageway for azoles that included a hydrophobic steric barrier causing a pause in movement toward the active site. The entrance to the active site acted like a second sieve to restrict access to the inner chamber. Collectively, these interactions impacted the final orientation of azoles reaching the active site and hence could explain differences in their biochemical properties observed in our previous studies, such as the consequences of methylation at position 5 of the azole ring. The association of a second azole demonstrated significant differences in interactions stabilizing the bound complex than observed for the first binding event. Intermolecular interactions occurred between the two azoles as well as CYP2E1 residue side chains and backbone and involved both hydrophobic contacts and hydrogen bonds. The relative importance of these interactions depended on the structure of the respective azoles indicating the absence of specific defining criteria for binding unlike the well-characterized dominant role of hydrophobicity in active site binding. Consequently, the structure activity relationships described here and elsewhere are necessary to more accurately identify factors impacting the observation and significance of cooperativity in CYP2E1 binding and catalysis toward drugs, dietary compounds, and pollutants.

Structural and energetic analysis to provide insight residues of CYP2C9, 2C11 and 2E1 involved in valproic acid dehydrogenation selectivity

Docking and molecular dynamics (MD) simulation have been two computational techniques used to gain insight about the substrate orientation within protein active sites, allowing to identify potential residues involved in the binding and catalytic mechanisms. In this study, both methods were combined to predict the regioselectivity in the binding mode of valproic acid (VPA) on three cytochrome P-450 (CYP) isoforms CYP2C9, CYP2C11, and CYP2E1, which are involved in the biotransformation of VPA yielding reactive hepatotoxic intermediate 2-n-propyl-4-pentenoic acid (4nVPA). There are experimental data about hydrogen atom abstraction of the C4-position of VPA to yield 4nVPA, however, there are not structural evidence about the binding mode of VPA and 4nVPA on CYPs. Therefore, the complexes between these CYP isoforms and VPA or 4nVPA were studied to explore their differences in binding and energetic stabilization. Docking results showed that VPA and 4nVPA are coupled into CYPs binding site in a similar conformation, but it does not explain the VPA hydrogen atom abstraction. On the other hand, MD simulations showed a set of energetic states that reorient VPA at the first ns, then making it susceptible to a dehydrogenation reaction. For 4nVPA, multiple binding modes were observed in which the different states could favor either undergo other reaction mechanism or ligand expulsion from the binding site. Otherwise, the energetic and entropic contribution point out a similar behavior for the three CYP complexes, showing as expected a more energetically favorable binding free energy for the complexes between CYPs and VPA than with 4nVPA.