Major cytochrome inhibitors may result in drug–drug interactions.
1. Introduction
Roughly seventy-four percent (74%) of all physician office visits involve drug therapy, with 2.9 billion drugs prescribed annually
[1]. These prescriptions accounted for USD 335 billion in patient costs in 2020
[2]. Half of Americans take at least one prescription drug and nearly a quarter of persons use three or more medications concomitantly
[1]. An estimated 6.7% of hospitalized patients experience an adverse drug reaction with a fatality rate of 0.32%, responsible for 106,000 deaths annually
[3]. These data suggest adverse drug reactions are the fourth leading cause of death in the United States
[3][4][5]. Furthermore, these statistics exclude ambulatory adverse drug reactions and are therefore likely an underrepresentation of the morbidity and mortality burden of adverse drug reactions. Nearly 350,000 adverse drug reactions occur in nursing homes annually, where polypharmacy is commonplace
[6]. Many of these drug–drug interactions are avoidable with proper patient medication management and by deprescribing. While efforts to reduce drug–drug interactions are important, a reduction in the number of medications cannot always be achieved without impacting therapeutic outcomes.
Adverse drug reactions commonly occur due to drug–drug interactions as a result of major cytochrome P450 (CYP) inhibition
[7]. CYPs are heme-containing membrane-bound enzymes that predominantly reside in the smooth endoplasmic reticulum and mitochondria of hepatocytes and in the intestines
[8][9]. In mammals, 57 CYP isoforms have been identified with the function of performing the oxidative metabolism of xenobiotics and endogenous compounds
[8]. Of the 57 CYP isoforms, 5 CYPs (CYPs 3A4, 2D6, 2C19, 2C9, and 1A2) are responsible for the metabolism of more than 80% of clinically used drugs
[7] (Figure 1). Enzyme inhibition occurs when two co-administered drugs share an identical mechanism of biotransformation and compete for metabolism in the same enzyme receptor site
[9]. As such, the more potent inhibitor will prevail, interrupting biotranformation, resulting in the diminished metabolism of the competing drug. This results in elevated serum drug levels of the unmetabolized drug, increasing the probability of adverse toxicological outcomes. This is especially important for drugs with a narrow therapeutic index, such as chemotherapeutics. It is important to note that inhibitors of major cytochromes can promote adverse drug reactions when administered as monotherapy. Inhibitors of major cytochromes may alter their own metabolism at high doses, leading to toxic plasma drug accumulation. These interactions pose a significant safety risk and should be examined closer during the initial drug design process. The CYP system and drug biotransformation has been reviewed extensively elsewhere
[10].
Figure 1. Cartoon depiction of cytochrome inhibition and pie chart. (
A) Orally bioavailable drugs are absorbed from the intestinal lumen into the mesenteric capillaries and transported to the liver via the portal vein. Many drugs undergo phase I biotransformation by major cytochromes within hepatocytes. This process enzymatically converts lipid-soluble compounds to more water-soluble compounds to facilitate the excretion. In the event of polypharmacy and/or the co-administration of CYP inhibitors, serum accumulation of unmetabolized drugs may occur, leading to untoward toxicities. (
B) Pie chart demonstrating the breakdown of the five major cytochromes and their contributions to the oxidative metabolism of xenobiotics and organic endogenous molecules.
Traditionally, pharmacokinetics has been neglected in preclinical drug design. Roughly four in five investigational new drugs fail as a result of absorption, distribution, metabolism, excretion and toxicity (ADMET) issues
[11]. Tools that assist in ADMET prediction may be advantageous in early drug design by reducing the drug attrition rate in advanced preclinical- and early clinical-stage studies. In particular, the design of drugs devoid of major cytochrome inhibition would reduce the burden of drug–drug interactions in patients where polypharmacy is unavoidable. Several reviews exist that focus on platforms that utilize a diverse range of prediction methods
[12][13][14][15][16][17][18][19][20]. Many of these reviews focus on outdated tools, as well as platforms that are exclusively available through a fee-based licensing model. More recently, numerous machine learning-based platforms have emerged that are user-friendly, are not as computationally expensive as pre-existing methods, and demonstrate improved prediction accuracy.
2. Machine Learning-Based Methods
Current in vitro and in vivo experimental approaches to investigate drug ADMET properties are highly effective in assessing pharmacokinetics; however, such approaches are extraordinarily expensive and are commonly performed late in drug development during IND-enabling studies
[13]. In silico investigation of CYP interactions early in drug development allows investigators to select for compounds that are less likely to fail in late-stage preclinical or early clinical trials due to undesirable pharmacokinetics (Figure 2). Such tools exist and are effective at predicting CYP inhibitors, common sites of metabolism (SoMs), and prospective drug metabolites. Hereby we highlight computational approaches available for ADMET predictions, focusing primarily on CYP inhibition and metabolism. The aim of this discussion will be to highlight emerging platforms that utilize machine learning-based methods, with a special emphasis placed on applications that are free to the user and available through a public webserver. A complete overview of platforms available through a fee-based licensing model, as well as older free prediction methods and methods for analyzing CYP molecular dynamics, are reviewed elsewhere
[13][14][15][16][17][18][19][20].
Figure 2. Workflow for incorporating ADMET predictions into early-stage drug design and development. Drug candidates can be readily screened using high-throughput computational approaches. Drugs with acceptable pharmacokinetic parameters as defined by the investigator will proceed to in vitro screening. Drugs with unacceptable pharmacokinetic properties will be triaged early in drug development. The following approach will result in a reduction in ADMET-related drug failures, reducing the cost of drug discovery efforts.
Various platforms capable of making enzyme-inhibitor predictions exist. These techniques use an array of methods, ranging from molecular docking with flexible multi-dimensional QSAR to machine learning based approaches. It should be noted that flexible docking-based methods, such as VirtualToxLab, are computationally expensive, take longer to receive results, and commonly require a software license
[13][21]. Machine learning-based tools are capable of producing accurate results in a timely fashion and can be hosted via a publicly available web server. These platforms tend to be more user-friendly, and are commonly able to screen a large number of structures. These platforms each use unique descriptors and utilize their own methods of data generation. As such, the accuracy of said platforms is subject to variability, contingent on the descriptors and molecular fingerprints used, the number of input structures used to build and train the model, and how the data are interpreted. In this section, we review the best publicly available machine learning-based platforms free to the scientific community.
3. Conclusions
More than half of drug failures are a result of poor ADMET. More specifically, drug–drug interactions as a result of major cytochrome inhibition are commonplace and can lead to life-threatening adverse drug reactions. A majority of pre-clinical-stage drug discovery programs focus solely on binding affinity and selectivity
[22]. Although binding affinity is of utmost importance, ADMET properties ensure that a drug is properly absorbed, distributed, and cleared; features that equally contribute to the success of a drug in the clinic. As such, in silico approaches have been developed to be introduced during early drug screening, with the goal of reducing the rate of drug attrition (Table 1).
Table 1. Summary of publicly available machine learning-based methods for cytochrome inhibition and metabolism.
Name |
CYPs |
Prediction Type |
ML Method |
No. Structures |
Avg. Accuracy 1 |
Additional Features |
pkCSM |
3A4, 2D6, 2C19, 2C9, 1A2 |
Inhibition |
Graph-based signatures |
18,000 |
0.810 (0.780–0.853) |
Comprehensive ADMET predictions (23 total) |
DeepCYP |
3A4, 2D6, 2C19, 2C9, 1A2 |
Inhibition |
Multitask autoencoder deep neural network |
13,000 |
0.864 (0.809–0.968) |
Assigns probabilities for CYP inhibition |
SuperCYPs-Pred |
3A4, 2D6, 2C19, 2C9, 1A2 |
Inhibition |
Random forests |
41,963 2 |
0.930 (0.840–0.970) |
Assigns probabilities for CYP inhibition |
vNN-ADMET |
3A4, 2D6, 2C19, 2C9, 1A2 |
Inhibition |
Variable nearest neighbors |
6261 |
0.890 (0.870–0.910) |
|
AdmetSAR 2.0 |
3A4, 2D6, 2C19, 2C9, 1A2 |
Inhibition |
Random forests, support vector machines, k-nearest neighbors |
96,000 3 |
0.784 (0.645–0.855) |
Comprehensive ADMET predictions (47 total); ADMETopt for lead optimization |
SwissADME |
3A4, 2D6l 2C19, 2C9, 1A2 |
Inhibition |
Support vector machines |
16,561 4 |
0.794 (0.720–0.800) |
Predictions of physicochemical properties, pharmacokinetics, and drug likeness; high throughput |
CypRules |
3A4, 2D6, 2C19, 2C9, 1A2 |
Inhibition |
Decision trees |
16,561 |
0.812 (0.730–0.900) |
High throughput |
CypReact |
3A4, 2E1, 2D6,2C19, 2C9, 2C8, 2B6, 2A6, 1A2 |
Sites of metabolism |
LBM learning algorithm |
2685 |
Unavailable |
Metabolite predictions; Additional CYPs |