Predicting routes of phase I and II metabolism based on quantum mechanics and machine learning.

1. Unexpected metabolism could lead to the failure of many late-stage drug candidates or even the withdrawal of approved drugs. Thus, it is critical to predict and study the dominant routes of metabolism in the early stages of research. In this study, we describe the development and validation of a 'WhichEnzyme' model that accurately predicts the enzyme families most likely to be responsible for a drug-like molecule's metabolism. Furthermore, we combine this model with our previously published regioselectivity models for Cytochromes P450, Aldehyde Oxidases, Flavin-containing Monooxygenases, UDP-glucuronosyltransferases and Sulfotransferases - the most important Phase I and Phase II drug metabolising enzymes - and a 'WhichP450' model that predicts the Cytochrome P450 isoform(s) responsible for a compound's metabolism. The regioselectivity models are based on a mechanistic understanding of these enzymes' actions, and use quantum mechanical simulations with machine learning methods to accurately predict sites of metabolism and the resulting metabolites. We train heuristic based on the outputs of the 'WhichEnzyme', 'WhichP450', and regioselectivity models to determine the most likely routes of metabolism and metabolites to be observed experimentally. Finally, we demonstrate that this combination delivers high sensitivity in identifying experimentally reported metabolites and higher precision than other methods for predicting in vivo metabolite profiles.

Predicting Regioselectivity of Cytosolic Sulfotransferase Metabolism for Drugs.

Cytosolic sulfotransferases (SULTs) are a family of enzymes responsible for the sulfation of small endogenous and exogenous compounds. SULTs contribute to the conjugation phase of metabolism and share substrates with the uridine 5'-diphospho-glucuronosyltransferase (UGT) family of enzymes. UGTs are considered to be the most important enzymes in the conjugation phase, and SULTs are an auxiliary enzyme system to them. Understanding how the regioselectivity of SULTs differs from that of UGTs is essential from the perspective of developing novel drug candidates. We present a general ligand-based SULT model trained and tested using high-quality experimental regioselectivity data. The current study suggests that, unlike other metabolic enzymes in the modification and conjugation phases, the SULT regioselectivity is not strongly influenced by the activation energy of the rate-limiting step of the catalysis. Instead, the prominent role is played by the substrate binding site of SULT. Thus, the model is trained only on steric and orientation descriptors, which mimic the binding pocket of SULT. The resulting classification model, which predicts whether a site is metabolized, achieved a Cohen's kappa of 0.71.

Predicting Regioselectivity of AO, CYP, FMO, and UGT Metabolism Using Quantum Mechanical Simulations and Machine Learning.

Unexpected metabolism in modification and conjugation phases can lead to the failure of many late-stage drug candidates or even withdrawal of approved drugs. Thus, it is critical to predict the sites of metabolism (SoM) for enzymes, which interact with drug-like molecules, in the early stages of the research. This study presents methods for predicting the isoform-specific metabolism for human AOs, FMOs, and UGTs and general CYP metabolism for preclinical species. The models use semi-empirical quantum mechanical simulations, validated using experimentally obtained data and DFT calculations, to estimate the reactivity of each SoM in the context of the whole molecule. Ligand-based models, trained and tested using high-quality regioselectivity data, combine the reactivity of the potential SoM with the orientation and steric effects of the binding pockets of the different enzyme isoforms. The resulting models achieve κ values of up to 0.94 and AUC of up to 0.92.

Predicting reactivity to drug metabolism: beyond P450s-modelling FMOs and UGTs.

We present a study based on density functional theory calculations to explore the rate limiting steps of product formation for oxidation by Flavin-containing Monooxygenase (FMO) and glucuronidation by the UDP-glucuronosyltransferase (UGT) family of enzymes. FMOs are responsible for the modification phase of metabolism of a wide diversity of drugs, working in conjunction with Cytochrome P450 (CYP) family of enzymes, and UGTs are the most important class of drug conjugation enzymes. Reactivity calculations are important for prediction of metabolism by CYPs and reactivity alone explains around 70-85% of the experimentally observed sites of metabolism within CYP substrates. In the current work we extend this approach to propose model systems which can be used to calculate the activation energies, i.e. reactivity, for the rate-limiting steps for both FMO oxidation and glucuronidation of potential sites of metabolism. These results are validated by comparison with the experimentally observed reaction rates and sites of metabolism, indicating that the presented models are suitable to provide the basis of a reactivity component within generalizable models to predict either FMO or UGT metabolism.

WhichP450: a multi-class categorical model to predict the major metabolising CYP450 isoform for a compound.

In the development of novel pharmaceuticals, the knowledge of how many, and which, Cytochrome P450 isoforms are involved in the phase I metabolism of a compound is important. Potential problems can arise if a compound is metabolised predominantly by a single isoform in terms of drug-drug interactions or genetic polymorphisms that would lead to variations in exposure in the general population. Combined with models of regioselectivities of metabolism by each isoform, such a model would also aid in the prediction of the metabolites likely to be formed by P450-mediated metabolism. We describe the generation of a multi-class random forest model to predict which, out of a list of the seven leading Cytochrome P450 isoforms, would be the major metabolising isoforms for a novel compound. The model has a 76% success rate with a top-1 criterion and an 88% success rate for a top-2 criterion and shows significant enrichment over randomised models.

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.

Aminoalkyl phenyl sulfones--a novel series of 5-HT7 receptor ligands.

A novel series of 5-HT(7) receptor ligands has been identified and evaluated, providing compounds showing a broad spectrum of functional activities and good selectivity over selected receptors and ion channels.

Multiple structural features of steroids mediate subtype-selective effects on human alpha4beta3delta GABAA receptors.

Neurosteroids have been shown to mediate some of their physiological effects via a modulatory site on type A inhibitory gamma-aminobutyric acid (GABAA) receptors. In particular, recent evidence has implicated selective potentiation of the delta subunit of GABAA receptors as an important mediator of in vitro and in vivo neurosteroid activity. However, this has been demonstrated for only a very small number of steroids, so both the generality of this finding, and the structural features of steroids which mediate functional delta-selectivity, are unclear. We have used a potentiometric assay based on fluorescence resonance energy transfer to measure GABA-activated responses in L(tk-) cells stably transfected with human GABAA receptor alpha4beta3delta and alpha4beta3gamma2 receptor subtypes. A set of 28 steroids were evaluated on these subtypes to characterise their functional potency and efficacy in modulating GABA responses. For most compounds there was a clear separation of their efficacy profiles between the receptor subtypes, with a substantially larger maximal response at the alpha4beta3delta receptor. 5beta-Pregnan-3beta-ol-20-one, 5beta-pregnane-3alpha,20beta-diol and 5beta-pregnane-3alpha,17alpha-diol-11,20-dione showed particularly high efficacy for alpha4beta3delta. No compounds were identified that simply inhibited responses at delta-containing receptors. However, 5beta-pregnane-3alpha,17alpha,20beta-triol, prednisolone 21-acetate, 4-pregnene-17alpha,20alpha-diol-3-one-20-acetate, 4-pregnen-20alpha-ol-3-one, and 5beta-pregnane-3alpha,17alpha,21-triol-20-one inhibited, though did not abolish, GABA responses at the alpha4beta3gamma2 subtype, while evoking modest-amplitude potentiation of alpha4beta3delta responses. Molecular modelling on this compound series using principal components analysis indicates that several structural features of steroids underlie their relative functional selectivity for potentiation of delta-containing GABAA receptors.