Toxic Metabolite Formation from Troglitazone (TGZ): New Insights from a DFT Study

Machine Learning Enables Accurate Prediction of Quinone Formation during Drug Metabolism

Metabolism helps in the elimination of drugs from the human body by making them more hydrophilic. Sometimes, drugs can be bioactivated to highly reactive metabolites or intermediates during metabolism. These reactive metabolites are often responsible for the toxicities associated with the drugs. Identification of reactive metabolites of drug candidates can be very helpful in the initial stages of drug discovery. Quinones are soft electrophiles that are generated as reactive intermediates during metabolism. Quinones make up more than 40% of the reactive metabolites. In this work, a reliable data set of 510 molecules was used to develop machine learning and deep learning-based predictive models to predict the formation of quinone-type metabolites. For representing molecules, two-dimensional (2D) descriptors, PubChem fingerprints, electro-topological state (E-state) fingerprints, and metabolic reactivity-based descriptors were used. Developed models were compared to the existing Xenosite web server using the untouched test set of 102 molecules. The best model achieved an accuracy of 86.27%, while the Xenosite server could achieve an accuracy of only 52.94% on the test set. Descriptor analysis revealed that the presence of greater numbers of polar moieties in a molecule can prevent the formation of quinone-type metabolites. In addition, the presence of a nitrogen atom in an aromatic ring and the presence of metabolophores V51, V52, and V53 (SMARTCyp descriptors) decrease the probability of quinone formation. Finally, a tool based on the best machine learning models was developed, which is accessible at http://14.139.57.41/quinonepred/.

Remdesivir: Mechanism of Metabolic Conversion from Prodrug to Drug

BACKGROUND

Remdesivir (GS-5734) has emerged as a promising drug during the challenging times of COVID-19 pandemic. Being a prodrug, it undergoes several metabolic reactions before converting to its active triphosphate metabolite. It is important to establish the atomic level details and explore the energy profile of the prodrug to drug conversion process.

METHODS

In this work, Density Functional Theory (DFT) calculations were performed to explore the entire metabolic path. Further, the potential energy surface (PES) diagram for the conversion of prodrug remdesivir to its active metabolite was established. The role of catalytic triad of Hint1 phosphoramidase enzyme in P-N bond hydrolysis was also studied on a model system using combined molecular docking and quantum mechanics approach.

RESULTS

The overall energy of reaction is 11.47 kcal/mol exergonic and the reaction proceeds through many steps requiring high activation energies. In the absence of a catalyst, the P-N bond breaking step requires 41.78 kcal/mol, which is reduced to 14.26 kcal/mol in a catalytic environment.

CONCLUSION

The metabolic pathways of model system of remdesivir (MSR) were completely explored completely and potential energy surface diagrams at two levels of theory, B3LYP/6-311++G(d, p) and B3LYP/6-31+G(d), were established and compared. The results highlight the importance of an additional water molecule in the metabolic reaction. The P-N bond cleavage step of the metabolic process requires the presence of an enzymatic environment.

Computer-Aided Estimation of Biological Activity Profiles of Drug-Like Compounds Taking into Account Their Metabolism in Human Body

Most pharmaceutical substances interact with several or even many molecular targets in the organism, determining the complex profiles of their biological activity. Moreover, due to biotransformation in the human body, they form one or several metabolites with different biological activity profiles. Therefore, the development and rational use of novel drugs requires the analysis of their biological activity profiles, taking into account metabolism in the human body. In silico methods are currently widely used for estimating new drug-like compounds' interactions with pharmacological targets and predicting their metabolic transformations. In this study, we consider the estimation of the biological activity profiles of organic compounds, taking into account the action of both the parent molecule and its metabolites in the human body. We used an external dataset that consists of 864 parent compounds with known metabolites. It is shown that the complex assessment of active pharmaceutical ingredients' interactions with the human organism increases the quality of computer-aided estimates. The toxic and adverse effects showed the most significant difference: reaching 0.16 for recall and 0.14 for precision.

A simple model to solve a complex drug toxicity problem

Linear drug toxicity models like therapeutic index (TI), physicochemical rules (rule of five, 3/75), ligand efficiency indices (LEI), ideal pharmacokinetic (PK) and pharmacodynamic (PD) profiles are widely used in drug discovery and development. In spite of this, predicting drug toxicity at various stages remains challenging and the overall productivity (<20%) and ultimate benefit to the patients remain low. A simple drug toxicity model, "Drug Toxicity Index" (DTI), is developed here using 711 oral drugs. DTI redefines drug toxicity as scaled biphasic and exponential functions of PD, PK and physicochemical parameters. PD, PK and physicochemical toxicity contributions were estimated from the on and off target IC₅₀, maximum unbound plasma drug concentration (free C _max), and log D values, respectively. These contributions are then scaled by molar dose and oral bioavailability and the logarithm of the sum of scaled contributions is DTI. Drugs with DTI above the WHO ATC drug category specific average values consistently have toxic profiles, while drugs with DTI below this average are relatively safe. DTI performs better than standard rules for lead optimization, LEI and exposure based TIs in identifying safe and toxic drugs. DTI classifies 392 drugs reported in the US-FDA's Liver Toxicity Knowledge Base (LTKB) with an AUC for ROC curves of 0.91-0.64 for different WHO ATC categories. DTI has been used to predict network meta-analysis results on relative toxicity within/across eight different therapeutic areas. It is useful in understanding PD, PK and physicochemical toxicity contributions and identifying potentially toxic drugs and the toxicity of recently approved drugs. Decision trees are proposed for applying the DTI concept in preclinical drug discovery and clinical trial settings. DTI can potentially reduce failure in drug discovery and might be useful in therapeutic drug monitoring and in xenobiotic and environmental toxicity studies.

Natural 6-hydroxy-chromanols and -chromenols: structural diversity, biosynthetic pathways and health implications

We present the first comprehensive and systematic review on the structurally diverse toco-chromanols and -chromenols found in photosynthetic organisms, including marine organisms, and as metabolic intermediates in animals. The focus of this work is on the structural diversity of chromanols and chromenols that result from various side chain modifications. We describe more than 230 structures that derive from a 6-hydroxy-chromanol- and 6-hydroxy-chromenol core, respectively, and comprise di-, sesqui-, mono- and hemiterpenes. We assort the compounds into a structure-activity relationship with special emphasis on anti-inflammatory and anti-carcinogenic activities of the congeners. This review covers the literature published from 1970 to 2017.

CC/B–N substitution in five membered heterocycles. A computational analysis

Novel five-membered azaboroles are aromatic, stable under neutral conditions, isomer stabilization energy is explained using σ-bond and aromatic stabilization energies.

3-Methylene-2,4-chromandione in situ trapping: introducing molecular diversity on 4-hydroxycoumarin

3-Methylene-2,4-chromandione trapped in a solid-state stable Mannich adduct, is released in solution.

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.

DFT study on the oxygen transfer mechanism in nitroethenediamine based H2-receptor antagonists using the bis-dithiolene complex as the model catalyst for N-oxide reductase enzyme

Nitroethenediamine is an important functional unit, which is present in H2-receptor antagonists. These drugs show low bioavailability due to the bacterial degradation caused by the N-oxide reductase type of enzymes present in the human colon. Quantum chemical studies have been carried out to elucidate the mechanism of metabolic degradation of nitroethenediamine in the active site of N-oxide reductase. Three different pathways have been explored for the N-oxide bond cleavage by the model system, Mo(IV) bis-dithiolene complex [Mo(OMe)(mdt)2](-), (where mdt=1,2-dimethyl-ethene-1,2-dithiolate) using B3LYP/6-311+G(d,p) and M06/6-311+G(d,p) Density Functional Theory methods. The oxygen atom transfer from the nitrogen atom of nitroethenediamine to the Mo(IV) complex, involves simultaneous weakening of the N-oxide bond and the formation of Mo-O bond through a least motion path. During this transfer, Mo center is converted from a square pyramidal geometry to a distorted octahedral geometry, to facilitate the process of oxygen atom transfer. The energy barrier for the oxygen atom transfer from the imine tautomer has been estimated to be 25.9kcal/mol however, the overall reaction has been found to be endothermic. On the other hand, oxygen transfer reaction from the nitronic acid tautomer requires 30.5kcal/mol energy leading to a highly exothermic metabolite (M-1) directly hence, this path can be considered thermodynamically favorable for this metabolite. The alternative path involving the oxygen atom transfer from the enamine tautomer requires comparatively a higher energy barrier (32.6kcal/mol) and leads to a slightly endothermic metabolite. This study established the structural and energetic details associated with the Mo(IV) bis-dithiolene complex that catalyzes the degradation of nitroethenediamine based drug molecules.

Optimization of troglitazone derivatives as potent anti-proliferative agents: towards more active and less toxic compounds

Δ2-Troglitazone derivatives were shown to exhibit anti-proliferative activity in a PPARγ-independent manner. We prepared various compounds in order to increase their potency and decrease their toxicity towards non-malignant primary cultured hepatocytes. Many compounds induced viabilities less than 20% at 10 μM on various cancer cell lines. Furthermore, five of them showed hepatocyte viability of 80% or more at 200 μM. In addition, compounds 17 and 18 exhibited promising maximum tolerated doses on a murine model, enabling future investigations.

3D organotypic cultures of human HepaRG cells: a tool for in vitro toxicity studies

Drug-induced human hepatotoxicity is difficult to predict using the current in vitro systems. In this study, long-term 3D organotypic cultures of the human hepatoma HepaRG cell line were prepared using a high-throughput hanging drop method. The organotypic cultures were maintained for 3 weeks and assessed for (1) liver specific functions, including phase I enzyme and transporter activities, (2) expression of liver-specific proteins, and (3) responses to three drugs (acetaminophen, troglitazone, and rosiglitazone). Our results show that the organotypic cultures maintain high liver-specific functionality during 3 weeks of culture. The immunohistochemistry analyses illustrate that the organotypic cultures express liver-specific markers such as albumin, CYP3A4, CYP2E1, and MRP-2 throughout the cultivation period. Accordingly, the production rates of albumin and glucose, as well as CYP2E1 activity, were significantly higher in the 3D versus the 2D cultures. Toxicity studies show that the organotypic cultures are more sensitive to acetaminophen- and rosiglitazone-induced toxicity but less sensitive to troglitazone-induced toxicity than the 2D cultures. Furthermore, the EC50 value (2.7mM) for acetaminophen on the 3D cultures was similar to in vivo toxicity. In summary, the results from our study suggest that the 3D organotypic HepaRG culture is a promising in vitro tool for more accurate assessment of acute and also possibly for chronic drug-induced hepatotoxicity.

Density functional study on the cytochrome-mediated S-oxidation: identification of crucial reactive intermediate on the metabolic path of thiazolidinediones

S-Oxidation is an important cytochrome P450 (CYP450)-catalyzed reaction, and the structural and energetic details of this process can only be studied by using quantum chemical methods. Thiazolidinedione (TZD) ring metabolism involving initial S-oxidation leads to the generation of reactive metabolites (RMs) and subsequent toxicity forcing the withdrawal of the glitazone class of drugs, thus, the study of the biochemical pathway of TZD ring metabolism is a subject of interest. The S-oxidation of the TZD ring and the formation of the isocyanate intermediate (ISC) was implicated as a possible pathway; however, there are several questions still unanswered in this biochemical pathway. The current study focuses on the CYP450-mediated S-oxidation, fate of the sulfoxide product (TZDSO), ring cleavage to ISC, and formation of nucleophilic adducts. The process of S-oxidation was explored by using Cpd I (iron(IV)-oxo porphyrin, to mimic CYP450) at TZVP/6-311+G(d) basis set. The barriers were calculated after incorporating dispersion and solvent corrections. The metabolic conversion from TZDSO to ISC (studied at B3LYP/6-311++G(2df,3pd)//B3LYP/6-31+G(d)) required a novel protonated intermediate, TZDSOH(+). The effect of higher basis sets (6-311+G(d,p), aug-cc-pvqz) on this conversion was studied. TZDSOH(+) was observed to be more reactive and thermodynamically accessible than ISC, indicating that TZDSOH(+) is the actual reactive intermediate leading to toxicity of the TZD class of compounds.

Synthesis of new troglitazone derivatives: anti-proliferative activity in breast cancer cell lines and preliminary toxicological study

Breast cancer is the most prevalent cancer in women. The development of resistances to therapeutic agents and the absence of targeted therapy for triple negative breast cancer motivate the search for alternative treatments. With this aim in mind, we synthesised new derivatives of troglitazone, a compound which was formerly used as an anti-diabetic agent and which exhibits anti-proliferative activity on various cancer cell lines. Among the compounds prepared, some displayed micromolar activity against hormone-dependent and hormone-independent breast cancer cells. Furthermore, the influence of the compounds on the viability of primary cultures of human hepatocytes was evaluated. This enabled us to obtain for the first time interesting structure-toxicity relationships in this family of compounds, resulting in 6b and 8b, which show good anti-proliferative activities and poor toxicity towards hepatocytes, compared to troglitazone.