New enantioenriched β-indolyl ketones as aromatase inhibitors: Unraveling heme-ligand interactions by MD simulation and MMPBSA analysis

A series of enantioenriched β-indolyl ketones as aromatase inhibitors (AI) is synthesized through the Michael-type Friedel-Crafts alkylation of indole. A highly efficient bifunctionalized amino catalyst is developed to access structurally diverse β-indolyl ketones in high yields (up to 91%) and excellent enantioselectivity (enantiomeric ratio up to 98:2). All the synthesized compounds demonstrated promising aromatase inhibitory potential, where ortho-substituted analogs (3c and 3e) were found most active with IC50 values of 0.68 and 0.90 µM, respectively. Both of these compounds exhibited significant cytotoxicity (IC50 = 0.34 and 0.37 µM) against the MCF-7 breast cancer cell line in the (3-[4,5-dimethylthiazol-2-yl]-2,5-diphenyl tetrazolium bromide) assay. Molecular docking studies of the synthesized compounds demonstrate favorable binding interactions with the estrogens controlling CYP19A1 (3EQM) and metabolizing CYP3A4 (5VCC) enzymes. Molecular dynamic (MD) simulation analysis revealed the essentiality of heme-ligand interactions to build a stable protein-ligand complex. An average root mean square deviation of 0.35 nm observed during a 100-ns MD simulation and binding free energy in the range of -190 to -227 kJ/mol calculated by g_mmpbsa analysis authenticated the stability of the 3c-3EQM complex. ADMET and drug-likeness parameters supported the suitability of these indole derivatives as the drug lead to develop potent inhibitors for estrogen-dependent breast cancer.

Computational Predictive and Electrochemical Detection of Metabolites (CP-EDM) of Piperine

In this article, we introduce a proof-of-concept strategy, Computational Predictive and Electrochemical Detection of Metabolites (CP-EDM), to expedite the discovery of drug metabolites. The use of a bioactive natural product, piperine, that has a well-curated metabolite profile but an unpredictable computational metabolism (Biotransformer v3.0) was selected. We developed an electrochemical reaction to oxidize piperine into a range of metabolites, which were detected by LC-MS. A series of chemically plausible metabolites were predicted based on ion fragmentation patterns. These metabolites were docked into the active site of CYP3A4 using Autodock4.2. From the clustered low-energy profile of piperine in the active site, it can be inferred that the most likely metabolic position of piperine (based on intermolecular distances to the Fe-oxo active site) is the benzo[d][1,3]dioxole motif. The metabolic profile was confirmed by comparison with the literature, and the electrochemical reaction delivered plausible metabolites, vide infra, thus, demonstrating the power of the hyphenated technique of tandem electrochemical detection and computational evaluation of binding poses. Taken together, we outline a novel approach where diverse data sources are combined to predict and confirm a metabolic outcome for a bioactive structure.

Prediction of Sites of Metabolism of CYP3A4 Substrates Utilizing Docking-Derived Geometric Features

Cytochrome P450 3A4 (CYP3A4) is one of the major drug-metabolizing enzymes in the human body and is responsible for the metabolism of ∼50% of clinically used drugs. Therefore, the identification of the compound's sites of metabolism (SOMs) mediated by CYP3A4 is of utmost importance in the early stage of drug discovery and development. Herein, docking-based approaches incorporating geometric features were used for SOMs prediction of CYP3A4 substrates. The cross-docking poses of a relatively large data set containing 474 substrates were analyzed in depth, and a widely observed geometric pattern called the close proximity of SOMs was derived from the poses. On the basis of the close proximity, several structure-based models have been constructed, which demonstrated better performance than those structure-based models using the criterion of Fe-SOM distance. For further improving the prediction performance, the structure-based models were also combined with the well-known ligand-based model SMARTCyp. One combined model exhibited good performance on the SOMs prediction of an external substrate set containing kinase inhibitors, PROTACs, approved drugs, and some lead compounds.

A multicentric consortium study demonstrates that dimethylarginine dimethylaminohydrolase 2 is not a dimethylarginine dimethylaminohydrolase

Dimethylarginine dimethylaminohydrolase 1 (DDAH1) protects against cardiovascular disease by metabolising the risk factor asymmetric dimethylarginine (ADMA). However, the question whether the second DDAH isoform, DDAH2, directly metabolises ADMA has remained unanswered. Consequently, it is still unclear if DDAH2 may be a potential target for ADMA-lowering therapies or if drug development efforts should focus on DDAH2's known physiological functions in mitochondrial fission, angiogenesis, vascular remodelling, insulin secretion, and immune responses. Here, an international consortium of research groups set out to address this question using in silico, in vitro, cell culture, and murine models. The findings uniformly demonstrate that DDAH2 is incapable of metabolising ADMA, thus resolving a 20-year controversy and providing a starting point for the investigation of alternative, ADMA-independent functions of DDAH2.

Drug-drug interactions that alter the exposure of glucuronidated drugs: Scope, UDP-glucuronosyltransferase (UGT) enzyme selectivity, mechanisms (inhibition and induction), and clinical significance

Drug-drug interactions (DDIs) arising from the perturbation of drug metabolising enzyme activities represent both a clinical problem and a potential economic loss for the pharmaceutical industry. DDIs involving glucuronidated drugs have historically attracted little attention and there is a perception that interactions are of minor clinical relevance. This review critically examines the scope and aetiology of DDIs that result in altered exposure of glucuronidated drugs. Interaction mechanisms, namely inhibition and induction of UDP-glucuronosyltransferase (UGT) enzymes and the potential interplay with drug transporters, are reviewed in detail, as is the clinical significance of known DDIs. Altered victim drug exposure arising from modulation of UGT enzyme activities is relatively common and, notably, the incidence and importance of UGT induction as a DDI mechanism is greater than generally believed. Numerous DDIs are clinically relevant, resulting in either loss of efficacy or an increased risk of adverse effects, necessitating dose individualisation. Several generalisations relating to the likelihood of DDIs can be drawn from the known substrate and inhibitor selectivities of UGT enzymes, highlighting the importance of comprehensive reaction phenotyping studies at an early stage of drug development. Further, rigorous assessment of the DDI liability of new chemical entities that undergo glucuronidation to a significant extent has been recommended recently by regulatory guidance. Although evidence-based approaches exist for the in vitro characterisation of UGT enzyme inhibition and induction, the availability of drugs considered appropriate for use as 'probe' substrates in clinical DDI studies is limited and this should be research priority.

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.

A Comprehensive Analysis of Human CYP3A4 Crystal Structures as a Potential Tool for Molecular Docking-Based Site of Metabolism and Enzyme Inhibition Studies

The notable ability of human liver cytochrome P450 3A4 (CYP3A4) to metabolize diverse xenobiotics encourages researchers to explore in-depth the mechanism of enzyme action. Numerous CYP3A4 protein crystal structures have been deposited in protein data bank (PDB) and are majorly used in molecular docking analysis. The quality of the molecular docking results depends on the three-dimensional CYP3A4 protein crystal structures from the PDB. Present review endeavors to provide a brief outline of some technical parameters of CYP3A4 PDB entries as valuable information for molecular docking research. PDB entries between 22 April 2004 and 2 June 2021 were compiled and the active sites were thoroughly observed. The present review identified 76 deposited PDB entries and described basic information that includes CYP3A4 from human genetic, Escherichia coli (E. coli) use for protein expression, crystal structure obtained from X-ray diffraction method, taxonomy ID 9606, Uniprot ID P08684, ligand–protein structure description, co-crystal ligand, protein site deposit and resolution ranges between 1.7[Formula: see text]Å and 2.95[Formula: see text]Å. The observation of protein–ligand interactions showed the various residues on the active site depending on the ligand. The residues Ala305, Ser119, Ala370, Phe304, Phe108, Phe213 and Phe215 have been found to frequently interact with ligands from CYP3A4 PDB. Literature surveys of 17 co-crystal ligands reveal multiple mechanisms that include competitive inhibition, noncompetitive inhibition, mixed-mode inhibition, mechanism-based inhibition, substrate with metabolite, inducer, or combination modes of action. This overview may help researchers choose a trustworthy CYP3A4 protein structure from the PDB database to apply the protein in molecular docking analysis for drug discovery.

Assessment of CYP2C9 Structural Models for Site of Metabolism Prediction

Structure-based prediction of a compound's potential sites of metabolism (SOMs) mediated by cytochromes P450 (CYPs) is highly advantageous in the early stage of drug discovery. However, the accuracy of the SOMs prediction can be influenced by several factors. CYP2C9 is one of the major drug-metabolizing enzymes in humans and is responsible for the metabolism of ∼13 % of clinically used drugs. In this study, we systematically evaluated the effects of protein crystal structure models, scoring functions, heme forms, conserved active-site water molecules, and protein flexibility on SOMs prediction of CYP2C9 substrates. Our results demonstrated that, on average, ChemScore and GlideScore outperformed four other scoring functions: Vina, GoldScore, ChemPLP, and ASP. The performance of the crystal structure models with pentacoordinated heme was generally superior to that of the hexacoordinated iron-oxo heme (referred to as Compound I) models. Inclusion of the conserved active-site water molecule improved the prediction accuracy of GlideScore, but reduced the accuracy of ChemScore. In addition, the effect of the conserved water on SOMs prediction was found to be dependent on the receptor model and the substrate. We further found that one of snapshots from molecular dynamics simulations on the apo form can improve the prediction accuracy when compared to the crystal structural model.

Differential inhibition of human CYP2C8 and molecular docking interactions elicited by sorafenib and its major N-oxide metabolite

The tyrosine kinase inhibitor sorafenib (SOR) is being used increasingly in combination with other anticancer agents like paclitaxel, but this increases the potential for drug toxicity. SOR inhibits several human CYPs, including CYP2C8, which is a major enzyme in the elimination of oncology drugs like paclitaxel and imatinib. It has been reported that CYP2C8 inhibition by SOR in human liver microsomes is potentiated by NADPH-dependent biotransformation. This implicates a SOR metabolite in enhanced inhibition, although the identity of that metabolite is presently unclear. The present study evaluated the capacity of the major N-oxide metabolite of SOR (SNO) to inhibit CYP2C8-dependent paclitaxel 6α-hydroxylation. The IC₅₀ of SNO against CYP2C8 activity was found to be 3.7-fold lower than that for the parent drug (14 μM versus 51 μM). In molecular docking studies, both SOR and SNO interacted with active site residues in CYP2C8, but four additional major hydrogen and halogen bonding interactions were identified between SNO and amino acids in the B-B' loop region and helixes F' and I that comprise the catalytic region of the enzyme. In contrast, the binding of both SOR and SNO to active site residues in the closely related human CYP2C9 enzyme was similar, as were the IC₅₀s determined against CYP2C9-mediated losartan oxidation. These findings suggest that the active metabolite SNO could impair the elimination of coadministered drugs that are substrates for CYP2C8, and mediate toxic adverse events, perhaps in those individuals in whom SNO is formed extensively.

Inhibition of Hepatic CYP2D6 by the Active N-Oxide Metabolite of Sorafenib

The multikinase inhibitor sorafenib (SOR) is used to treat patients with hepatocellular and renal carcinomas. SOR undergoes CYP-mediated biotransformation to a pharmacologically active N-oxide metabolite (SNO) that has been shown to accumulate to varying extents in individuals. Kinase inhibitors like SOR are frequently coadministered with a range of other drugs to improve the efficacy of anticancer drug therapy and to treat comorbidities. Recent evidence has suggested that SNO is more effective than SOR as an inhibitor of CYP3A4-mediated midazolam 1'-hydroxylation. CYP2D6 is also reportedly inhibited by SOR. The present study assessed the possibility that SNO might contribute to CYP2D6 inhibition. The inhibition kinetics of CYP2D6-mediated dextromethorphan O-demethylation were analyzed in human hepatic microsomes, with SNO found to be ~ 19-fold more active than SOR (K_is 1.8 ± 0.3 μM and 34 ± 11 μM, respectively). Molecular docking studies of SOR and SNO were undertaken using multiple crystal structures of CYP2D6. Both molecules mediated interactions with key amino acid residues in putative substrate recognition sites of CYP2D6. However, a larger number of H-bonding interactions was noted between the N-oxide moiety of SNO and active site residues that account for its greater inhibition potency. These findings suggest that SNO has the potential to contribute to pharmacokinetic interactions involving SOR, perhaps in those individuals in whom SNO accumulates.