Reactive Metabolites from Thiazole-Containing Drugs: Quantum Chemical Insights into Biotransformation and Toxicity

Drugs containing thiazole and aminothiazole groups are known to generate reactive metabolites (RMs) catalyzed by cytochrome P450s (CYPs). These RMs can covalently modify essential cellular macromolecules and lead to toxicity and induce idiosyncratic adverse drug reactions. Molecular docking and quantum chemical hybrid DFT study were carried out to explore the molecular mechanisms involved in the biotransformation of thiazole (TZ) and aminothiazole (ATZ) groups leading to RM epoxide, S-oxide, N-oxide, and oxaziridine. The energy barrier required for the epoxidation is 13.63 kcal/mol, that is lower than that of S-oxidation, N-oxidation, and oxaziridine formation (14.56, 17.90, and 20.20, kcal/mol respectively). The presence of the amino group in ATZ further facilitates all the metabolic pathways, for example, the barrier for the epoxidation reaction is reduced by ∼2.5 kcal/mol. Some of the RMs/their isomers are highly electrophilic and tend to form covalent bonds with nucleophilic amino acids, finally leading to the formation of metabolic intermediate complexes (MICs). The energy profiles of these competitive pathways have also been explored.

Virtual screening of potentially endocrine-disrupting chemicals against nuclear receptors and its application to identify PPARγ-bound fatty acids

Nuclear receptors (NRs) are key regulators of energy homeostasis, body development, and sexual reproduction. Xenobiotics binding to NRs may disrupt natural hormonal systems and induce undesired adverse effects in the body. However, many chemicals of concerns have limited or no experimental data on their potential or lack-of-potential endocrine-disrupting effects. Here, we propose a virtual screening method based on molecular docking for predicting potential endocrine-disrupting chemicals (EDCs) that bind to NRs. For 12 NRs, we systematically analyzed how multiple crystal structures can be used to distinguish actives and inactives found in previous high-throughput experiments. Our method is based on (i) consensus docking scores from multiple structures at a single functional state (agonist-bound or antagonist-bound), (ii) multiple functional states (agonist-bound and antagonist-bound), and (iii) multiple pockets (orthosteric site and alternative sites) of these NRs. We found that the consensus enrichment from multiple structures is better than or comparable to the best enrichment from a single structure. The discriminating power of this consensus strategy was further enhanced by a chemical similarity-weighted scoring scheme, yielding better or comparable enrichment for all studied NRs. Applying this optimized method, we screened 252 fatty acids against peroxisome proliferator-activated receptor gamma (PPARγ) and successfully identified 3 previously unknown fatty acids with Kd = 100-250 μM including two furan fatty acids: furannonanoic acid (FNA) and furanundecanoic acid (FUA), and one cyclopropane fatty acid: phytomonic acid (PTA). These results suggested that the proposed method can be used to rapidly screen and prioritize potential EDCs for further experimental evaluations.

Modular mimicry and engagement of the Hippo pathway by Marburg virus VP40: Implications for filovirus biology and budding

Ebola (EBOV) and Marburg (MARV) are members of the Filoviridae family, which continue to emerge and cause sporadic outbreaks of hemorrhagic fever with high mortality rates. Filoviruses utilize their VP40 matrix protein to drive virion assembly and budding, in part, by recruitment of specific WW-domain-bearing host proteins via its conserved PPxY Late (L) domain motif. Here, we screened an array of 115 mammalian, bacterially expressed and purified WW-domains using a PPxY-containing peptide from MARV VP40 (mVP40) to identify novel host interactors. Using this unbiased approach, we identified Yes Associated Protein (YAP) and Transcriptional co-Activator with PDZ-binding motif (TAZ) as novel mVP40 PPxY interactors. YAP and TAZ function as downstream transcriptional effectors of the Hippo signaling pathway that regulates cell proliferation, migration and apoptosis. We demonstrate that ectopic expression of YAP or TAZ along with mVP40 leads to significant inhibition of budding of mVP40 VLPs in a WW-domain/PPxY dependent manner. Moreover, YAP colocalized with mVP40 in the cytoplasm, and inhibition of mVP40 VLP budding was more pronounced when YAP was localized predominantly in the cytoplasm rather than in the nucleus. A key regulator of YAP nuclear/cytoplasmic localization and function is angiomotin (Amot); a multi-PPxY containing protein that strongly interacts with YAP WW-domains. Interestingly, we found that expression of PPxY-containing Amot rescued mVP40 VLP egress from either YAP- or TAZ-mediated inhibition in a PPxY-dependent manner. Importantly, using a stable Amot-knockdown cell line, we found that expression of Amot was critical for efficient egress of mVP40 VLPs as well as egress and spread of authentic MARV in infected cell cultures. In sum, we identified novel negative (YAP/TAZ) and positive (Amot) regulators of MARV VP40-mediated egress, that likely function in part, via competition between host and viral PPxY motifs binding to modular host WW-domains. These findings not only impact our mechanistic understanding of virus budding and spread, but also may impact the development of new antiviral strategies.

By screening an array of 115 mammalian WW-domains with the PPxY motif from MARV VP40 (mVP40), we identified YAP1 and TAZ, transcriptional effectors of the Hippo pathway, as mVP40 interactors, and demonstrated that ectopically expressed YAP1 or TAZ inhibited budding of mVP40 virus-like particles (VLPs) in a WW-domain/PPxY dependent manner. Angiomotin (Amot), a multi-PPxY containing regulator of YAP1 nuclear/cytoplasmic localization and function, rescued mVP40 VLP egress from either YAP1- or TAZ-mediated inhibition in a PPxY-dependent manner. Indeed, endogenous Amot expression was critical for egress of mVP40 VLPs and authentic MARV. In sum, we have revealed a link between the Hippo pathway and filovirus egress by identifying negative (YAP/TAZ) and positive (Amot) regulators of MARV VP40-mediated egress.

Exploring PfDHFR reaction surface: A combined molecular dynamics and QM/MM analysis

The substrate to the enzyme PfDHFR (Plasmodium falciparum Dihydrofolate Reductase) is a small molecule dihydrofolate (DHF), it gets converted to tetrahydrofolate (THF) in the active site of the enzyme. The PfDHFR reaction surface involves the protonation of DHF to DHFP as an initial step before the catalytic conversion. The binding affinities of all these species (DHF, DHFP and THF) contribute to the mechanism of DHFR catalytic action. Molecular dynamics (MD) simulations and Quantum Mechanics/Molecular Mechanics (QM/MM) analysis were performed to evaluate the binding affinity and molecular recognition interactions of the substrate DHF/DHFP and the product THF, in the active site of wild-type PfDHFR (wtPfDHFR). The binding affinities of the cofactor NADPH/NADP⁺ were also estimated in all the three complexes. The molecular dynamics (MD) simulations of the substrate, product and cofactor in the cavities of wtPfDHFR revealed the variation of the atomic level interactions during the course of the catalytic conversion. It was found that the DHFP binds very strongly to the PfDHFR active site and pulls the cofactor NADPH closer to itself. The QM/MM analysis revealed that the binding energy of DHFP (-59.82 kcal/mol) and NADPH (-100.24 kcal/mol) in DHFP-wtPfDHFR complex, is higher in comparison to the binding energy of DHF (-38.67 kcal/mol) and NADPH (-77.53 kcal/mol) in DHF-wtPfDHFR complex and the binding energy of THF (-30.72 kcal/mol) and NADP⁺ (-73.72 kcal/mol) in THF-wtPfDHFR complex. The hydride ion donor-acceptor distance (DAD) analysis was also carried out. This combined MD and QM/MM analysis revealed that the protonation of DHF increases the proximity between the substrate and the cofactor, thus facilitates the reaction profile of PfDHFR.

Copper(ii)-catalyzed Chan–Lam cross-coupling: chemoselective N-arylation of aminophenols

Copper(ii)-catalyzed boronic acid promoted chemoselective N-arylation of unprotected aminophenols has been developed.

Toxicity Originating from Thiophene Containing Drugs: Exploring the Mechanism using Quantum Chemical Methods

Drug metabolism of thiophene containing substrates by cytochrome P450s (CYP450) leads to toxic side effects, for example, nephrotoxicity (suprofen, ticlopidine), hepatotoxicity (tienilic acid), thrombotic thrombocytopenic purpura (clopidogrel), and aplastic anemia (ticlopidine). The origin of toxicity in these cases has been attributed to two different CYP450 mediated metabolic reactions: S-oxidation and epoxidation. In this work, the molecular level details of the bioinorganic chemistry associated with the generation of these competitive reactions are reported. Density functional theory was utilized (i) to explore the molecular mechanism for S-oxidation and epoxidation using the radical cationic center Cpd I [(iron(IV)-oxo-heme porphine system with SH(-) as the axial ligand, to mimic CYP450s] as the model oxidant, (ii) to establish the 3D structures of the reactants, transition states, and products on both the metabolic pathways, and (iii) to examine the potential energy (PE) profile for both the pathways to determine the energetically preferred toxic metabolite formation. The energy barrier required for S-oxidation was observed to be 14.75 kcal/mol as compared to that of the epoxidation reaction (13.23 kcal/mol) on the doublet PE surface of Cpd I. The formation of the epoxide metabolite was found to be highly exothermic (-23.24 kcal/mol), as compared to S-oxidation (-8.08 kcal/mol). Hence, on a relative scale the epoxidation process was observed to be thermodynamically and kinetically more favorable. The energy profiles associated with the reactions of the S-oxide and epoxide toxic metabolites were also explored. This study helps in understanding the CYP450-catalyzed toxic reactions of drugs containing the thiophene ring at the atomic level.