REVIEW article

Front. Chem., 11 September 2020

Sec. Medicinal and Pharmaceutical Chemistry

Volume 8 - 2020 | https://doi.org/10.3389/fchem.2020.00726

Computational Approaches in Preclinical Studies on Drug Discovery and Development

  • 1. Key Laboratory of Big Data Mining and Precision Drug Design of Guangdong Medical University, Research Platform Service Management Center, Dongguan, China

  • 2. Key Laboratory of Pesticide & Chemical Biology, Ministry of Education, College of Chemistry, Central China Normal University, Wuhan, China

  • 3. The Second School of Clinical Medicine, Guangdong Medical University, Dongguan, China

  • 4. Key Laboratory for Research and Development of Natural Drugs of Guangdong Province, School of Pharmacy, Guangdong Medical University, Dongguan, China

  • 5. Marine Biomedical Research Institute of Guangdong Zhanjiang, Zhanjiang, China

Article metrics

View details

379

Citations

40,6k

Views

11,9k

Downloads

Abstract

Because undesirable pharmacokinetics and toxicity are significant reasons for the failure of drug development in the costly late stage, it has been widely recognized that drug ADMET properties should be considered as early as possible to reduce failure rates in the clinical phase of drug discovery. Concurrently, drug recalls have become increasingly common in recent years, prompting pharmaceutical companies to increase attention toward the safety evaluation of preclinical drugs. In vitro and in vivo drug evaluation techniques are currently more mature in preclinical applications, but these technologies are costly. In recent years, with the rapid development of computer science, in silico technology has been widely used to evaluate the relevant properties of drugs in the preclinical stage and has produced many software programs and in silico models, further promoting the study of ADMET in vitro. In this review, we first introduce the two ADMET prediction categories (molecular modeling and data modeling). Then, we perform a systematic classification and description of the databases and software commonly used for ADMET prediction. We focus on some widely studied ADMT properties as well as PBPK simulation, and we list some applications that are related to the prediction categories and web tools. Finally, we discuss challenges and limitations in the preclinical area and propose some suggestions and prospects for the future.

Introduction

Drug development is a complicated, risky, and time-consuming process that can be divided into several stages, including disease-related genomics, target identification and validation, lead discovery and optimization, preclinical studies, and clinical trials (Tang et al., 2006) (Figure 1). During early drug discovery, the activities and specificities of candidate drugs are usually assessed at an early stage, and pharmacokinetics and toxicities are evaluated at a relatively late stage (Selick et al., 2002). However, the undesirable efficacy and safety, mainly caused by absorption, distribution, metabolism, excretion, and toxicity (ADMET) characteristics, resulted in the failure of many candidate drugs in the final stage (Caldwell et al., 2009). Cook et al. (2014) comprehensively reviewed the results of AstraZeneca's small-molecule drug projects from 2005 to 2010 based on a longitudinal study. They found that unacceptable safety and toxicity were the most important reasons for the failure of more than half of all project closures. As with the development of drug discovery, it was realized that it is important to filter and optimize the ADMET properties for drugs at an early stage, which has been accepted and widely used to reduce the attrition rate in drug research and development. A “fail early, fail cheap” strategy has been employed by many pharmaceutical companies (Yu and Adedoyin, 2003). Pharmacokinetics and toxicity assessments of preclinical drugs are of great value in reducing the failure rate of new chemical entities (NCEs) in clinical trials (Kola and Landis, 2004; Yang et al., 2018b; Ferreira and Andricopulo, 2019). In recent years, in vitro and in vivo ADMET prediction methods have been widely used, but it is impractical to perform complex and expensive ADMET experiments on a large number of compounds (Cheng et al., 2013; Patel C. N. et al., 2020). Thus, an in silico strategy to predict ADMET properties has become very attractive as a cost-saving and high-throughput alternative to experimental measurement methods.

Figure 1

With the rapid development of computer technologies, the high-throughput screening of compounds, application of combinatorial chemistry, and ability of compound synthesis have increased dramatically. The early demands for ADMET data on lead compounds have also significantly increased, and methods for evaluating ADMET in vitro are gradually increasing. Many in silico methods have been successfully applied to the in vitro prediction of ADMET, and in silico models have also been developed to replace in vivo models for the prediction of pharmacokinetics, toxicity, and other parameters (Zhu et al., 2011; Wang et al., 2015; Alqahtani, 2017). In silico ADMET prediction has progressed with the continuous development of cheminformatics and has entered the era of big data (Ferreira and Andricopulo, 2019). Two in silico approach categories can be used for ADMET prediction: molecular modeling and data modeling. Molecular modeling is based on the three-dimensional structures of proteins. It includes multiple methods such as molecular docking, molecular dynamics (MD) simulation, and quantum mechanics (QM) calculation (Bowen and Guener, 2013; Cheng et al., 2013; Silva-Junior et al., 2017). Data modeling includes quantitative structure–activity relationship (QSAR) (Cumming et al., 2013) and physiologically-based pharmacokinetic (PBPK) modeling (Fan and de Lannoy, 2014). Due to the increase in number of properties that need to be predicted, a series of ADMET software programs capable of comprehensive property prediction have been developed. The development from in silico approaches to ADMET software has undergone a long process of predicting property parameters from less to more at early to late timepoints (Figure 2). This review first provides a detailed introduction to the two in silico approaches of ADMET prediction. Then, we summarize the widely used databases and software related to ADMET prediction. Finally, we analyze the problems and challenges faced by computer model prediction methods as well as the tools, and we propose some of our own prospects for future development in this area.

Figure 2

In silico Approaches

Molecular Modeling

Molecular modeling, based on the three-dimensional structures of proteins, is an important category in predicting ADMET properties and includes methods such as pharmacophore modeling, molecular docking, MD simulations, and QM calculations (Figure 3). As more and more three-dimensional structures of ADMET proteins become available, molecular modeling can complement or even surpass QSAR studies (Moroy et al., 2012). Applying molecular modeling to perform ADMET prediction is a challenge because the ADMET proteins usually have flexible and large binding cavities. Many promising results of molecular modeling in predicting compound metabolism have been reported. The methods in these cases can be generally divided into ligand-based and structure-based and help not only to analyze metabolic properties but also to further optimize compound toxicity, bioavailability, and other parameters (Lin et al., 2003).

Figure 3

Ligand-Based Methods

Ligand-based methods derive information on proteins' active sites based on the shapes, electronic properties, and conformations of inhibitors, substrates or metabolites; this information depends on the assumption that the metabolic properties of compounds are entirely the result of their chemical structures and characteristics (de Groot et al., 2004; Andrade et al., 2014). In this category, pharmacophore modeling is one of the most widely used methods. The interactions between ligands and receptors can be predicted by constructing a pharmacophore model to cover the structures or properties of ligands in three-dimensional space and then to simulate the spatial and chemical properties of binding sites (de Groot, 2006). Therefore, the availability of ligand data is essential to the construction of pharmacophore models. In recent years, there have been many cases of using pharmacophore models to screen promising compounds with outstanding ADMET properties (Nandekar et al., 2016; El-Zahabi et al., 2019; Mohan et al., 2020; Patel D. B. et al., 2020; Rawat and Verma, 2020). For example, Nandekar et al. (2016) generated and validated a pharmacophore model to screen anticancer compounds acting via cytochrome P450 1A1 (CYP1A1). Nine compounds that have preferred pharmacophore characteristics and are capable of generating reactive metabolites were finally selected for further study. Rawat and Verma (2020) developed a pharmacophore model to discover new dual target inhibitors of Plasmodium falciparum dihydroorotate dehydrogenase (PfDHODH) and cytochrome bc1 complex (PfCytbc1) to treat malaria. The molecule MMV007571, which has been validated as an efficient multi-target inhibitor, was used to extract features from the binding information for the model construction. The model was used to screen a library including more than 40,000 molecules. After a series of experiments, two compounds were developed with the desired properties in binding potential and pharmacokinetic characters.

The shape-focused method is another category of ligand-based methods. This method is based on the fact that the shapes of a ligand and the binding site of its receptor should be complementary. Thus, molecules that have a comparable shape should be able to bind to the same receptor. It is more possible for a ligand to bind with the same target if it has greater similarity with the reference molecule (Putta and Beroza, 2007). This method requires only one reference molecule (shape) to perform a screening. However, more models should be constructed to cover more diverse chemical space if different shapes are available (Perez-Nueno and Ritchie, 2011), particularly for highly flexible proteins. Some studies using the shape-focused method have been reported in recent years (Reddy et al., 2013; Chen et al., 2015; Kumar et al., 2015; Prabhu and Singh, 2019). For example, Chen et al. (2015) presented a shape-based virtual screening to find new cores for the design of acetylcholinesterase (AChE) inhibitors. The shape of commercial inhibitor tacrine was used to search for new potential inhibitors. Two hit compounds were finally identified with good ADMET properties and better activity than tacrine.

With the improvement of computer hardware performance, the time-consuming QM calculation in ADMET prediction has become possible and popular. The QM calculation can be used to evaluate the bond break, which is a step required for metabolic transformation (Andrade et al., 2014). Moreover, this calculation uses an accurate means of describing electrons in atoms and molecules (Modi, 2003). Hence, QM calculation is very helpful in ADMET prediction. Extensively increasing studies involving the application of QM methods have been conducted to describe ADMET properties of new compounds (Li et al., 2012; Taxak et al., 2013; Kavitha et al., 2015; Sasahara et al., 2015; Evangelista et al., 2016; Mondal et al., 2017). The ab initio (Hartree-Fock), semiempirical (AM1 and PM3), and density functional theory (DFT) approaches are most commonly used in these studies (Silva-Junior et al., 2017). For example, Mondal et al. (2017) used the DFT method to study the absorption profile and antimicrobial activity for five sulfonamide Schiff bases. The geometries compared well with the experimental value. Sasahara et al. (2015) used the DFT method to evaluate the metabolic selectivity of antipsychotic thioridazine by CYP450 2D6. This study revealed the importance of the substrate orientation in the reaction center of this enzyme for the metabolic reaction.

All ligand-based methods need to address the problem of the uncertainty of metabolic enzyme binding sites. If reliable structural data for a metabolic enzyme are lacking, the properties of ligand binding to the enzyme can only be speculative, and minor modifications to a ligand may cause a significant decrease in ligand-protein affinity (Kirchmair et al., 2012). Therefore, it is difficult to predict ADMET using ligand-based methods without reliable protein structure data. However, the use of ligand-based methods in metabolism prediction can easily eliminate inappropriate compounds and reduce the number of compounds that fail during the synthetic evaluation cycle and more expensive late stages.

Structure-Based Methods

Structure-based methods can be used not only to predict the ADMET properties of compounds, but also to study specific interactions between small molecules and ADMET proteins. In general, they focus on obtaining binding modes from the static structures of protein-ligand complexes, regardless of time-dependent conformational fluctuations (Cheng T. et al., 2012). For structure-based methods, changes in receptor conformation must be considered mainly because of the possible interaction between the existing large and flexible binding cavities with diverse ligands. By performing MD calculation to simulate the dynamic changes of spatial shape, we may obtain an adequate conformation sampling to search a stable and reliable binding mode for ADMET prediction. In recent years, structure-based methods have been widely used, for instance, to predict the binding patterns of substrates (Macalino et al., 2015), conformational changes in enzymes (Cheng T. et al., 2012), and their catalytic effects on physiological systems (Cui and Karplus, 2003), to evaluate substrate affinity, instability and metabolic pathways (Sun and Scott, 2010), and to assess the relationship between metabolism and carcinogenicity (Fratev and Benfenati, 2008).

The molecular modeling strategy makes important contributions to the rationalization of metabolic reactions of compounds, allowing the simulation of binding modes between drugs and macromolecules in the ADMET process at atomic or molecular levels. With the rapid development of structural elucidation and pharmacokinetic calculation techniques, structure-based methods are becoming increasingly predictive and accurate. However, the molecular modeling strategy is still limited by drawbacks such as a requirement to accurately analyze the structural flexibility of proteins (Kazmi et al., 2019). Additionally, high-resolution experimental structural data of the target will be more conducive to our accurate prediction of the drug metabolic fate (Doss et al., 2014). Surely, combining improved structure-based and ligand-based methods can create synergistic effects in metabolic prediction, enabling more comprehensive descriptions of metabolic reactions (Issa et al., 2017; Kar and Leszczynski, 2017).

Data Modeling

There are two widely used data modeling methods to predict ADMET-related properties, QSAR, and PBPK modeling (Cheng et al., 2013). ADMET analysis and prediction in QSAR mainly depends on many molecular descriptors, including topological, geometrical physicochemical, or electronic descriptors. Many properties, such as blood-brain barrier (BBB), clinical adverse effects, percent protein binding (%PPB), lipophilicity (logP), preclinical toxicological endpoints, and metabolism of pharmaceutical substances, can be predicted using the QSAR method. The PBPK modeling always predicts parameters concerning the dose size and dose frequency, such as the volume of distribution at steady-state (Vss), total drug clearance (CL), and fraction of dose that reaches the portal vein (Fabs), because most drugs are taken orally.

QSAR

QSAR, which employs mathematical models to describe relationships between molecular structures and their biological activities, has been used in pharmaceutical chemistry since the 1960s (Hansch, 1981). The classical QSAR developed by Hansch was used to predict ADMET (Hansch, 1981). It is mainly based on the hypothesis that similar molecules exhibit similar properties (Patterson et al., 1996). Thus, two prediction methods are primarily considered: (1) prediction based on molecular similarity (pharmacophore- and molecular fragment-based methods) and (2) prediction based on property similarity (log P, log D, and others) (Yongye and Medina-Franco, 2013). In these two methods, the accuracy of prediction depends on the attribute characteristics in the applicability domain contained in the training set. Thus, when using a model constructed with a specific training set, prediction should be performed using compounds that have a similar structure space to those in the training set to improve the prediction accuracy of the QSAR model, since compounds with similar distribution in chemical space are more likely to exhibit similar biological activities (Huang and Fan, 2011).

QSAR is a method for using various biochemical and physical data to construct models. In QSAR studies, compounds can be mathematically codified as molecular descriptors, and the relationship between molecular descriptors and defined properties is constructed by statistical methods, after which a generated model is used to predict the corresponding properties of new compounds (Michielan and Moro, 2010). This method first transforms molecular structure into molecular descriptors, which are then used to establish prediction models by using statistical approaches or machine learning techniques such as support vector machine (SVM) and K Nearest Neighbor (kNN) (Wang S. et al., 2016; Wu et al., 2019; Yang et al., 2019; Fu et al., 2020). For example, Schyman et al. (2017) used the variable nearest neighbor (vNN) method to develop 15 ADMET prediction models and to use them to quickly assess some potential drug candidates, including toxicity, microsomal stability, mutagenicity, and likelihood of causing drug-induced liver injury. Belekar et al. (2015) developed a computational model to identify compounds as breast cancer resistance protein (BCRP) inhibitors or not by using various machine learning approaches like SVM, k-NN, and the artificial neural network (ANN). The prediction accuracy of all three approaches was over 85%. Finally, internal and external cross-validation were performed to confirm the reliability of the model before it is used on new predictions to find molecules outside the training set (Figure 4).

Figure 4

QSAR modeling uses a large number of descriptors that allow lookups, enable structure/response associations, and help with similarity and substructure searches (Khan, 2010). Most of the available descriptors can be divided into three categories: (1) two-dimensional molecular topology information; (2) three-dimensional molecular structure; and (3) physicochemical and electronic descriptors, which are commonly used to predict ADMET-related properties (Danishuddin and Khan, 2016; Tabeshpour et al., 2018; Zhang et al., 2018). Jiang et al. (2020) developed a series of QSAR models by using 379 molecular descriptors to discriminate BCRP inhibitors. The descriptors characterized the physicochemical, two-dimensional substructures, and drug-like properties of the studied compounds. Lapins et al. (2018) constructed a QSAR model to predict the lipophilicity of compounds by using a signature molecular descriptor, which is related to the molecular two-dimensional topology information from 1.6 million compounds. Xu et al. (2017) developed three deep learning-based QSAR models to evaluate the acute oral toxicity (AOT) of compounds. The atom and bond information extracted from over 2,000 two-dimensional molecule structures were used as descriptors to construct the models. The best model achieved an external prediction accuracy over 94%, which is more efficient than traditional fingerprints or descriptors. Bujak et al. (2015) predicted the permeability of the BBB of chemical compounds using molecular energy-related descriptors in combination with the well-known lipophilicity descriptors. The data indicate that the QSAR model has important information value, and these descriptors may have supportive value in predicting the blood brain distribution (Bujak et al., 2015). Therefore, accurate prediction of ADMET parameters mainly depends on the selection of a suitable modeling method, molecular descriptors of specific ADMET endpoints, and large experimental data sets related to these endpoints. Only in this way can the ADMET properties of the candidate compounds be predicted precisely.

At present, many tools used for ADMET prediction have been developed based on QSAR methods. These tools utilize different descriptors to define the collected data, and then the mathematical model fitted from the training set is used to predict the properties. We listed three widely used QSAR-based ADEMT prediction tools and related studies herein. (1) The Danish QSAR Database (http://qsar.food.dtu.dk/) collected estimates from over 200 QSAR models from free and commercial platforms, including descriptors like ecotoxicity, environmental fate, physicochemical properties, and ADMET. Trivedi et al. (2019, 2020) used the online tool Danish QSAR database to determine the ADMET properties for potential hits for Dengue fever and H1N1 flu, respectively. Multicase acute aquatic toxicity, carcinogenicity, arylhydroxylase activity, lethal body burden, bioconcentration, mutagenicity, biodegradation, environmental partitioning, and general properties are included in the ADMET properties. Finally, 12 compounds were identified as potential leads against dengue fever and 18 compounds against H1N1 flu. (2) The OCED Toolbox (https://www.oecd.org/chemicalsafety/risk-assessment/oecd-qsar-toolbox.htm), a package for toxicity prediction, was also developed based on QSAR. Han et al. (2019) used OECD QSAR Toolbox 4.1 to predict the genotoxicity for ceftazidime (CAZ) and its impurities to improve quality control of drugs. (3) ADMET PredictorTM (https://www.simulations-plus.com/software/admetpredictor/) is another tool utilizing QSAR to predict ADMET parameters of compounds. Takac et al. (2019) used ADMET PredictorTM to investigate the potential impact and safety profile with respect to the environment and health for 25 selected entactogen molecules. The chemical structure (including 1D and 2D) information was used as the input for ADMET PredictorTM. Lipophilicity parameters, volume of distribution, jejunal permeability, solubility, and logarithm of the brain/blood partition coefficient were predicted in this case. Alarn and Khan (2019) used the ADMET PredictorTM to predict pharmacokinetics, pharmacodynamics, and toxicity parameters of flavone analogs to reveal their anticancer activity. Different physicochemical properties were calculated as descriptors to build the model, and then numerous properties were predicted, such as solubility, lipophilicity, permeability, absorption, bioavailability, BBB, transporters, plasma-protein binding, and volume of distribution.

Although the use of QSAR models has made considerable progress in ADMET prediction, these models cannot yet be used to replace in vitro or in vivo studies for all endpoints. The QSAR method is always limited by its model expansion capability, and large experimental data are always needed for model construction. The narrow data distribution may induce over fitting and lead to inaccurate prediction results. For example, Verheyen et al. (2017) estimated the QSAR models used for the prediction of eye and skin irritation/corrosion in Derek Nexus, Toxtree and Case Ultra. They found that the prediction results were unsatisfactory because of the narrow application range and low accuracy. Thus, validation and documentation for a constructed model is important prior to use.

PBPK Modeling

Most traditional models for predicting drug pharmacokinetics are empirical models. With a deeper understanding of the pharmacokinetic mechanism of drugs, PBPK models have been developed to predict PK properties (De Buck and Mackie, 2007). The PBPK model is an arithmetical model that combines drug data (e.g., drug concentration and clearance rate) and species physiology parameters to replicate the PK profile of a drug in plasma and tissues, aiming to describe in vivo drug pharmacokinetics that are related to tissue volume, administration routes, blood flow, biotransformation pathways, and interactions with tissues or organs in the body (Espie et al., 2009; Jones et al., 2015). The origin of the PBPK models can be traced back to Teorell's work in 1935. Teorell introduced a multicompartment model to simulate pharmacokinetics, which organically combined physiology and biology for the first time (Teorell, 1935; Zhao et al., 2011). Teorell's work has since attracted serious attention to the PBPK model.

In recent years, PBPK modeling has been substantially improved, making it more widely applicable for the research and development of drugs (Edginton et al., 2008; Rowland et al., 2011; Zhuang and Lu, 2016). Moreover, the increase in preclinical data, especially in vitro data, has promoted the development of PBPK models and simulations in drug discovery (Zhuang and Lu, 2016). PBPK modeling describes the physical and biological disposition of each compartment by dividing organisms into individual organs, and the most common processes are related to blood transportation and penetration, distribution between blood and organ tissue, and metabolic excretion, among others (Schmitt and Willmann, 2004). Since PBPK integrates large amounts of drug-specific data, parameters, and species physiology (systematic data), there are two kinds of parameters in PBPK models, which use the concentration-time curves of all organs and blood as output information (Nestorov, 2003). The first type consists of physiological parameters, such as tissue volume, blood flow, and cardiac output. Recently, due to the extensive application of in vitro-in vivo extrapolation (IVIVE) in PBPK, many researchers have predicted the disposal of drugs in vivo through in vitro metabolism and transport data, indicating that metabolic enzymes and transporter expression data have become important physiological parameters (Rostami-Hodjegan, 2012). The second type consists of drug-related parameters, such as the plasma ratio, organ/blood partition coefficient, and permeability (Nestorov, 2003). Recently, PBPK models have been widely constructed to predict drug-related parameters (Chow et al., 2016; Pathak et al., 2019; Song et al., 2020). For example, Chow et al. used a physiologically based model to predict drug solubility and effective permeability (Chow et al., 2016) to examine the potential impact of excipients on oral drug absorption.

Databases

In the past 10 years, with rapid development, a number of related databases storing pharmacokinetic parameters have emerged. We collected some of the most commonly used databases and classified them as ADMET-related databases and auxiliary databases. A brief introduction to these databases, including website links, data scales, and descriptions, is provided in Tables 1, 2. For the ADMET-related databases, users can submit information on the compounds they want to query through the corresponding modules. Then, shape screening or pharmacophore screening will be performed to obtain additional targets or bioactivity information on similar ligands that match the query molecule. The ADMET-related properties can also be obtained from the query result. The auxiliary databases mainly focus on providing structural information about compounds. Although some ADMET-related information is provided in the search results, it is not complete, and not every compound is associated with such information.

Table 1

Database nameAvailabilityDescriptionURL
ACD/PerceptaCommercialUsed to calculate PhysChem properties, ADME properties, and toxicity endpointshttps://www.acdlabs.com/products/percepta/index.php
ADME DatabaseCommercialA commercial database for studying drug interactions and ADME, updated quarterlyhttps://www.fujitsu.com/jp/group/kyushu/en/solutions/industry/lifescience/admedatabase/
AurSCOPE ADMECommercialA complete annotated, structured knowledge base can be used to design predictive models and identify potential drug interactionshttp://www.aureus-sciences.com/
KnowitallCommercialThe world's largest mass spectrometry library, can be used to quickly and accurately identify the spectrumhttp://www.bio-rad.com/
ToxBankCommercialA cross-cluster project can be used for comprehensive data analysis of toxicology and alternative detection of repeated dose toxicity testshttp://toxbank.net/
ACToRFreeExploring and visualizing complex computational toxicology informationhttps://actor.epa.gov/actor/home.xhtml
BRENDAFreeDatabase of enzyme and enzyme-ligand informationhttp://www.brenda-enzymes.org/
CEBSFreeA toxicological resource that can be used to model, predict analyze and assess effects of time and dose on responses to experimental conditionshttps://manticore.niehs.nih.gov/cebssearch
ChemTunesFreeA unique cheminformatics platform and expert QC'ed database that facilitate and support the safety and risk assessment process for chemical substanceshttps://www.mn-am.com/products/chemtunes
CTDFreeA public resource storing scientific data about relationships between genes, chemicals, and human diseaseshttp://ctdbase.org/
DIDBFreeAllows researchers to perform the assessment of human PK-based drug interactions and drug safetyhttps://www.druginteractioninfo.org/
DSSToxFreeA public resource of chemical structure information with existing toxicity data for supporting better predictive toxicologyhttps://www.epa.gov/chemical-research/distributed-structure-searchable-toxicity-dsstox-database
eChemPortalFreeChemical substance search, chemical property data search and GHS search are availablehttp://www.oecd.org/chemicalsafety/risk-assessment/echemportalglobalportaltoinformationchemicalsubstances.htm
IDA2PMFreeAn integrated database for data access transformation and analysis visualization and prognostic modelinghttp://idaapm.helsinki.fi/
Liceptor DatabaseFreeA ligand database containing 2D structures, related molecular descriptors and bioactivity data consisting of assays, functions and therapeutic hintshttp://www.evolvus.com/
Metabolism and Transport DatabaseFreeA database related to small molecule transport metabolism that can be used for computational analysis and modelinghttp://www-metrabase.ch.cam.ac.uk/
MetaCycFreeA database of metabolic pathways that can predict the metabolic pathways of sequenced genomeshttps://metacyc.org/
NTPFreeNTP studied at any chemical with the potential to impact healthhttps://sandbox.ntp.niehs.nih.gov/neurotox/
Pharmaco Kinetics Knowledge Base (PKKB)FreeA free database that collects ADMET data and can be used for ADMET modelinghttp://cadd.zju.edu.cn/pkkb/
RepdoseFreeA toxicity database storing in vitro and in vivo data and a prediction system for the safety and risk evaluationhttps://repdose.item.fraunhofer.de/about_repdose.html
RTECSFreeContains additional information related to the chemical industry and occupational safety and health that can be used to assess workers' exposure to chemicalshttps://www.atsdr.cdc.gov/substances/index.asp
SuperToxicFreeCan perform similarity screening, risk assessment, and link to other databaseshttp://bioinformatics.charite.de/supertoxic/
T3DBFreeDescribe the relationship between toxins and targets and the mechanism of toxic action, which can be used for prediction of toxicity, prediction of toxic targets, etc.http://www.t3db.ca/
TOXNETFreeIt is composed of a set of databases dealing with environmental health of toxicological hazardous chemicals and related fieldshttps://toxnet.nlm.nih.gov/
ToxRefDBFreeStore data from in vivo animal toxicity tests and provide toxicity endpoints for predictive modelinghttps://catalog.data.gov/dataset/toxcast-toxrefdb

Some widely used ADMET-related databases.

Table 2

Database NameAvailabilityDescriptionScaleURL
PubChemFreeA public database storing biological properties of small moleculesContains biological test results for more than 700,000 compoundshttp://pubchem.ncbi.nlm.nih.gov/
DrugBankFreeStoring information about drugs and related targetsContains 13,490 drug entrieshttps://www.drugbank.ca/
STITCHFreeA search tool and resource for interactions of chemicals and proteinsInteraction data for more than 68, 000 different chemicalshttp://stitch.embl.de/
ChEMBLdbFreeAn open large bioactive database that can display metabolite pathways and link to metabolite metabolizing enzymes and information document source dataContaining more than 1.6 million individual compound structures denoted in the database, with 14 million activity values from more than 1.2 million assayshttps://www.ebi.ac.uk/chembldb/
BindingDBFreeA web-based resource containing experimental binding affinities, focusing mainly on the interactions between potential drug-targets and drug-like molecules1,794,819 binding data, for 7,438 protein targets and 796,104 small moleculeshttp://www.bindingdb.org/bind/index.jsp
ChemProtFreeA resource to perform in silico estimation of small molecules with the integration of molecular, cellular, and disease-related proteins complexesover 1.7 million compounds with 7.8 million bioactivity measurements for 19,504 proteinshttp://potentia.cbs.dtu.dk/ChemProt/
SIDERFreeContains information about marketed drugs and their adverse reactions, which can be used to quickly track trace the origin of an extracted side effect5,868 side effect and 1,430 drugs and 139756 drug-SE pairshttp://sideeffects.embl.de/drugs/
MetaADEDBFreeA comprehensive database of adverse drug reactions (adr), which is used to predict clinical adr and the source of adverse drug side effects of personalized drugs, and to predict drug interactions with targets and drug action patterns3,060 chemicals (including more than 1,300 FDA approved and experimental drugs) and 13, 256 ADEshttp://lmmd.ecust.edu.cn/online_services/metaadedb/
TTDFreeProvides data about known and studied therapeutic macromolecules, the targeted disease, pathway information, and the associated drugs directed at each of these targets2,589 targets and 31,614 drugshttp://bidd.nus.edu.sg/group/cjttd/
KEGGFreeA reference repository that is widely used to integrate and interpret large-scale datasets obtained by genome sequencing and other high-throughput experimental technologies4 categories (systems, genomic, chemical and health information) from 18 databaseshttps://www.kegg.jp/

Some auxiliary databases for ADMET prediction.

ADMET-Related Databases

At present, many in silico models are used to predict ADMET, but massive amounts of data are needed to build them. The quality and quantity of the data are closely related to the accuracy of model prediction, so reliable experimental data are the key to successful prediction (Dearden, 2007; Alqahtani, 2017). Currently, there are some databases that can help ADMET prediction, such as the ADME Database (Shang et al., 2017), SuperToxic (Schmidt et al., 2009), PKKB (Cao et al., 2012), and DSSTox (Williams et al., 2017). By using these databases, users can obtain helpful data sets for use in external algorithms to generate prediction models. The databases can also be used directly to perform prediction through search functions, such as the similarity search or prediction. In addition, these models can be updated as new experimental data are added to the database.

The ADME Database (https://www.fujitsu.com/jp/group/kyushu/en/solutions/industry/lifescience/admedatabase/), developed by Zagreb University and Fujitsu in 2004, is a commercial database that specializes in pharmacokinetics information. It provides comprehensive data on drug-metabolizing enzymes and drug transporters that are specific to humans. The data have been widely used in drug research and development, such as ADME prediction and drug-drug interactions. Users can search for classification, metabolic reactions, and kinetics-related information about compounds by structure or substructure. However, the database currently limits large-scale downloads of user data, as well as public dissemination of some models.

SuperToxic collects toxins from different sources (animals, plants, synthetic, etc.), compiles ~60 000 compounds with their structures, and integrates some chemical properties and commercial availability information (Schmidt et al., 2009). These compounds are classified based on their toxicity, which derives from more than 2 million measurements. The values can be used to study the relationship between chemical structures and functions of toxins for evaluating the risk of their use. Users can easily query the structure and toxicity information of all compounds with corresponding properties through a structure search, name search, or property search. SuperToxic also allows users to browse the data by choosing an alphabetic character or numbers to present all entries starting with the selection (Schmidt et al., 2009). The available CASRN or NSC numbers in the database can also be recorded. The toxicity information retrieved from the database includes the dosage, type of test (toxicity measurement, such as LD50), and cell lines or organisms that determine the toxicity (Schmidt et al., 2009). The database also integrates software packages that are widely used in modern composite database construction, such as Marvin Sketch (molecule drawing), JMol (visual inspection), and MyChem/OpenBabel (property calculation). Furthermore, SuperToxic was connected to the Protein Data Bank, UniProt, and KEGG databases to identify potential targets in biochemical pathways to search for compounds (Schmidt et al., 2009).

The EPA Distributed Structure-Searchable Toxicity (DSSTox) database, which provides a series of documented, standardized and complete structure annotated toxicity information files, can be very useful for SAR model development (Richard and Williams, 2002; Williams et al., 2017). To allow wider use of the database, DSSTox was designed to use a structure data file (SDF), a public and industry-standard import/export file format storing chemical structures and property information that can be used as input for any chemical relational database (CRD) application or converted to data tables. It is one of the best-curated public datasets available at present, and the data stored in it are regarded as a standard reference for publicly available structural toxicity-based data.

Except for the databases introduced above, some newly constructed ADMET-related databases should also be of concern, for example, Comparative Toxicogenomics Database (CTD) (Davis et al., 2019), The Toxicity Reference Database (ToxRefDB) (Watford et al., 2019), and The Chemical Effects in Biological Systems database (CEBS) (Lea et al., 2017). CTD is a powerful and public database designed to enhance understanding of how environmental exposures influence human health. It provides data on chemical–gene/protein interactions, chemical–disease and gene–disease relationships that are combined with pathway and function data to help develop hypotheses about the inherent mechanisms of diseases affected by the environment (Davis et al., 2019). ToxRefDB collects data from in vivo animal toxicity tests and provides toxicity endpoints for predictive modeling. Approximately 28,000 datasets from nearly 400 endpoints have been generated and stored. The recent update of ToxRefDB has added connections to other resources and significantly enhanced the utility of predicting toxicology (Watford et al., 2019). CEBS offers a toxicology resource that compiles individual and concise animal data from 11,000 test articles and more than 8,000 studies encompassing all available National Toxicology Program (NTP) carcinogenicity, genetic toxicity, and short-term toxicity studies. The high-quality data in CEBS is very useful for constructing a more accurate model for toxicity prediction (Lea et al., 2017). We may infer from the recent constructed databases that the quantity and quality of the data determine the quality of the model and will be our focus in the future.

Auxiliary Databases

In addition to ADMET-related databases, databases of biological activity, pathways, and side effects are important for ADMET prediction. Most of these databases are free and open to visitation, such as DrugBank (Wishart et al., 2018), PubChem (Kim et al., 2019), and ChEMBL (Gaulton et al., 2017). Although they are rarely used to predict ADMET-related properties directly, they can provide structural information to build models or be queried for information about compounds. Users can also download the predicted compound structure and use it as an input file for other software.

DrugBank, a comprehensive database, integrated thousands of well-studied drugs and drug targets with their physical, chemical, biological, and pharmaceutical data (Wishart et al., 2006, 2008). DrugBank 4.0 was further expanded to contain data on ADMET and other kinds of QSAR information (Law et al., 2013). DrugBank 5.0, the latest version, has further updated this information (Wishart et al., 2018). Users can use chemical shifts or mass-to-charge ratio (m/z) lists to search DrugBank's spectral library for approximate or exact matches. DrugBank also systematically classifies compounds into different types based on structural features and structural similarities and allows users to query it by using a simple text (Law et al., 2013).

Pubchem, a public database of small molecules with their biological properties, consists of three interconnected parts: (1) Compound, storing over 102 million unique chemical structures provided by various depositors; (2) Substance, containing more than 251 million records including complexes, extracts, mixtures, and non-characteristics; and (3) BioAssay, containing more than 1,067,000 bioassays, providing composite adjacent structures, substructures, similar structures, bioactivity data, and other search functions (Kim et al., 2019).

ChEMBL is an open data database containing two-dimensional structures, calculated properties (molecular weight, lipophilicity, etc.) and abstract biological activities (pharmacology and ADMET data) for numerous drug-like bioactive compounds (Gaulton et al., 2017). It is composed of three different datasets that were originally developed by Inpharmatica, including StARlite, CandiStore, and DrugStore (Overington, 2009). The data in ChEMBL were extracted from the scientific literature and designed to meet the needs of users to intelligently cluster relevant information and integrate data across therapeutic studies and areas.

Software

Favorable ADMET characteristics are important as early requirements for drug candidates to reduce late failure and cost. However, many ADMET properties are highly dependent on each other, so they need to be optimized simultaneously in preclinical studies on drug discovery and development. Nevertheless, concurrent optimization of multiparameter ADMET is the most difficult and least attractive stage. As a result, early prediction of ADMET involved only some simple properties, such as logP, logD, and logS. With increasing experimental data, an increasing number of in silico models were developed to predict more complex ADMET parameters, such as the human intestinal absorption rate, oral bioavailability, blood-brain barrier permeation rate, Caco-2 permeability, human intestinal absorption, drug interactions, P-glycoprotein, plasma protein binding rate, CYP metabolic enzymes, and kidney clearance (Pires et al., 2015; Dong et al., 2018) (Figure 5). People have also attempted to integrate these models to predict ADMET parameters concurrently, and many studies have described these in silico models and their predicted properties (Dickins and van de Waterbeemd, 2004; Wang et al., 2015). In addition, software integrating these models to predict ADMET parameters concurrently has been developed. Tables 3, 4 list some of these software packages (free and commercial) with their functions. We also compared five commonly used software in Table 5 with each other to visualize their detailed functions.

Figure 5

Table 3

SoftwareURLLogSLogPLogDSolTPSAHIACaco-2PPBBBBVdMetCLToxP-gpperpKaMDCK
ACD/I-labhttp://www.acdlabs.com/home/
ADMETlabhttp://admet.scbdd.com/
admetSARhttp://lmmd.ecust.edu.cn/admetsar1/
FAF-Drug4http://fafdrugs4.mti.univ-paris-diderot.fr.
Lazarhttps://www.in-silico.de
OCHEMhttp://ochem.eu
OECD Toolboxhttp://toolbox.oasis-lmc.org/
OSIRIS property explorerhttp://www.organic-chemistry.org/prog/peo/
pkCSMhttps://smartcyp.sund.ku.dk
SMARTCyphttps://nodepit.com/node/org.lhasalimited.knime.metabolism.encapsulated.smartcyp.SMARTCypNodeFactory
SwissADMEhttp://www.swissadme.ch
ToxCreatehttps://github.com/opentox/toxcreate
ToxTreehttp://toxtree.sourceforge.net/#carousel0
VCCLAB (ALOGPS 2.1)http://www.vcclab.org
vNN-ADMEThttps://vnnadmet.bhsai.org/vnnadmet/login.xhtml

Some free ADMET software and related properties in their prediction.

logS, aqueous solubility; LogP, octanol-water partition coefficient; LogD, octanol-water distribution coefficient; Sol, solubility; TPSA, topological polar surface area; HIA, human intestinal absorption; PPB, plasma protein binding; BBB, blood-brain barrier; Vd, volume of distribution; Met, metabolism; CL, clearance; Tox, toxicity; P-gp, P glycoprotein; Per, Permeability; pKa, acidity coefficient; MDCK, madin Darby canine kidney cell line.

Table 4

SoftwareURLLogSLogPLogDSolTPSAHIACaco-2PPBBBBVdMetCLToxP-gpPerpKaMDCK
ACD/Percepta Platformhttps://www.acdlabs.com/products/percepta/
ADMEWORKShttps://www.fujitsu.com/jp/group/kyushu/en/solutions/industry/lifescience/admeworks/
CompuDrug's Pallas Systemhttp://www.compudrug.com/pallas_system
Derek Nexushttps://www.lhasalimited.org/products/derek-nexus.htm
MCASE, CASE, CASETOXhttps://www.multicase.com/
MetaSitehttp://www.moldiscovery.com/soft_metasite.php
PASShttp://genexplain.com/pass/
Schrodinger QikProphttps://www.schrodinger.com/qikprop/
Simulations Plus ADMET Predictorhttps://www.simulations-plus.com/software/admetpredictor/
StarDrophttps://www.optibrium.com/stardrop/stardrop-p450-models.php
TIMEShttp://oasis-lmc.org/products/software/times.aspx
TOPKAThttp://www.moldiscovery.com/software/vsplus/
VolSurf+http://www.moldiscovery.com/software/vsplus/

Some commercial ADMET software and related properties in their prediction.

logS, aqueous solubility; LogP, octanol-water partition coefficient; LogD, octanol-water distribution coefficient; Sol, solubility; TPSA, topological polar surface area; HIA, human intestinal absorption; PPB, plasma protein binding; BBB, blood-brain barrier; Vd, volume of distribution; Met, metabolism; CL, clearance; Tox, toxicity; P-gp, P glycoprotein; Per, Permeability; pKa, acidity coefficient; MDCK, madin Darby canine kidney cell line.

Table 5

ToolsAvailabilityBatch computationEndpointsDatabaseDruglikeness rulesDruglikeness modelSystematic evaluationMedicinal chemistry friendliness propertiesPhysico- chemistry propertiesSimilarityQSAR modelAlgorithmsTraining setsPattern recognition
SwissADMEFreeNumber, 19 Contents: B, A, D, M××××××
ADMETlabFreeNumber: 31 Contents: B, A, D, M, E, T√ (288,967 entries; 5 similarity searching strategies)√ (5 rules)××
admetSAR 2.0Free×Number: 47 Contents: B, A, D, M, E, T√ (210,000 entries)××××
LazarFree×Number: 11 Contents: T: Acute toxicity; BBB; Carcinogenicity, LOAEL, Maximum Recommended Daily Dose, Mutagenicity××××××
ToxTreeFree×Number: 6 Contents: M, T×××××××××

Comparison of five commonly used ADMET software programs.

B, basic physicochemical property; A, absorption; D, distribution; M, metabolism; E, excretion; T, toxicity.

SwissADME is a hybrid web server that was developed by the Swiss Institute of Bioinformatics (Daina et al., 2017). It supports diverse input formats and can predict and analyze the ADME properties of numerous compounds in batches submitted from all over the world. This software outputs different types of physicochemical properties of drugs: water solubility, lipophilicity, physicochemical properties, druglikeness, pharmacokinetics, and medicinal chemistry, which can be directly exported and saved as a data file in CSV (comma-separated values) format and read by programs such as WordPad and Excel (Daina et al., 2017). In addition, it supplies a bioavailability radar map to quickly and intuitively evaluate the druglikeness of small molecules, facilitating its use for non-experts without professional knowledge (Daina et al., 2017). The server uses a variety of rules to evaluate the same property and provides the evaluation criteria and basis for most of the predicted data. However, in the prediction of whether the compound is a substrate or an inhibitor of the CYP enzyme, only a propensity is output, rather than a probability output similar to that of admetSAR. SwissADME also supplies a link for the one-click submission of the queried molecules to other servers in the Swiss series for further analysis (Daina et al., 2017).

Although in vivo toxicology is still the gold standard for identifying drug side effects, it is now believed that this method will not help reduce the large consumption rate in late clinical development (Merlot, 2010). Many computational tools have been developed to predict drug toxicity, helping to decrease the attrition rate of molecular compounds in drug discovery and reduce drug development time and cost. In recent years, the predictive power of these toxicology prediction systems has tremendously improved, covering more complex toxicological endpoints, such as hepatotoxicity, teratogenicity, nephrotoxicity, and carcinogenicity (Muster et al., 2008; Yang et al., 2018b). Currently, many commercial and free web-based toxicity predictors are available, such as Lazar (Maunz et al., 2013) and Toxtree (Patlewicz et al., 2008; Bhhatarai et al., 2016).

Lazar, developed by in silico toxicology GMBH, is a tool based on OpenTox (an integrated interface for an interoperable prognostic toxicology framework) to predict toxicological endpoints such as carcinogenicity, reproductive toxicity, and long-term toxicity (Hardy et al., 2010). It uses data mining algorithms to predict the toxicity of new compounds based on the experimental training data. Data sets with chemical structure and biological activity can be used as training data. Thus, Lazar can serve as a universal predictive algorithm for any biological endpoint if sufficient experimental data are available, so users no longer need to consider chemical, biological, or toxicological expertise but derive in silico models from statistical standards (Maunz et al., 2013). Users need only to input the structure of the compound, and Lazar will search the database for a series of similar compounds and corresponding experimental data, which will be used to construct a local QSAR model. The prediction results using the model will be displayed in a graphical interface, which provides structural features and compounds similar to the query compounds and toxicity properties for each fragment (Helma, 2006; Maunz et al., 2013).

Toxtree is a free software program that was commissioned by the European Chemicals Agency (ECB), the joint research center of the European Commission (Bhatia et al., 2015). It was originally designed to enable effective development of the Cramer decision tree. The latest version of Toxtree included additional projects, such as corrosion rules, BfR/SICRET skin irritation, and the Verhaar scheme, with a total of 14 functional modules (Bhatia et al., 2015). It includes physiochemical exclusion rules and structural alert inclusion rules, which are used to categorize compounds (Bhhatarai et al., 2016). Unlike Lazar, it has no training set. Its prediction is based on structural filters, so there is no applicability domain. It handles molecular structure information by using a decision tree model for risk assessment (Bhhatarai et al., 2016). Users can access it at http://toxtree.sourceforge.net/ to predict the toxicity of structures of interest.

ADMETlab is a platform for systematic ADMET estimation based on a comprehensive collection of ADMET databases (Dong et al., 2018). The platform includes four functional modules, which are used for drug similarity assessment (based on Lipinski's rule of five and the Druglikeness model), ADMET endpoint prediction, system evaluation, and database/similarity search (Dong et al., 2018). Among them, “ADMET prediction” is the main module used; the other three are auxiliary modules. Users can query one or more compounds with the platform by using the SMILES, uploading an SDF format file, or drawing the chemical structure using the embedded JME editor. After the compound is uploaded, the platform will use multiple pharmacokinetics models built by the different integrated data sets to make extensive predictions of the ADMET properties. Model prediction results are output in an interactive data table containing predicted values and structures. The software allows batch prediction, and users can apply the “drug similarity assessment” module to filter out compounds that are unlikely to be lead compounds or drugs, achieving the purpose of preliminary screening (Dong et al., 2018). Users must select one model to acquire results for one or multiple molecules, which is proper for screening compounds at specific endpoints, and the results will provide reasonable ADMET recommendations for each endpoint. Therefore, users can perform rapid screening of ADMET properties based on these independent specific prediction models and even further deliberately optimize the chemical structures of compounds, making them more likely to become drugs (Dong et al., 2018). Considering the very large amount of collected data and numerous constructed QSPR models, ADMETlab is currently one of the most comprehensive platforms used in ADMET prediction.

AdmetSAR is a free and comprehensive tool for ADMET property prediction (Cheng F. et al., 2012). The ADMET-related property data stored in AdmetSAR were collected from published literature. AdmetSAR also includes a searchable tool called ADMET-Simulator, which combines predictive and superior QSAR models in a toolbox based on chemical informatics and can predict ~50 ADMET endpoints. AdmetSAR enables users to easily search for ADMET properties by querying CASRN, the common name, or the structure (Yang et al., 2018a). The new version of admetSAR (version 2.0) mainly focuses on in silico prediction of chemical ADMET properties (Yang et al., 2018a). More than 40 predictive models trained by state-of-the-art machine learning methods were implemented in admetSAR. Four functions were developed: (1) customizable ADMET risk filters, (2) QSAR-based ADMET property prediction, (3) toxicity prediction for environmental chemicals, and (4) environmental hazard assessment. ADMETopt (Yang et al., 2018c) is a new module added in version 2.0 for lead compound optimization according to the predicted ADMET properties.

Except for ADMET software, an increasing number of PBPK software have been developed to perform systematics of the drug process in human body. The development of PBPK software has further promoted the use of PBPK modeling methods (Bouzom and Walther, 2008; Edginton et al., 2008; Perdaems et al., 2010). PBPK software was used to build and use the PBPK models, which can be useful for the estimation of pharmacokinetic parameters during the drug development process. At present, PBPK software can be roughly divided into two categories, the user customized software and user-friendly software. A brief introduction to the two types of software is provided in Table 6, including the company/institute and website links. In Table 7, we compare the functions of the commonly used software packages for both categories, such as WinNonlin (https://www.certara.com/) and GastroPlus (https://www.simulations-plus.com/).

Table 6

SoftwareCompany/InstituteURL
USER CUSTOMIZED SOFTWARE
KineticaThermo Fisher Scientific Inc.http://kineticadownload.com/Kinetica5.1-SP1/Default.asp
MATLAB-simulinkThe MathWorks Inc.https://www.mathworks.com/
NONMENICONhttps://www.iconplc.com/innovation/nonmem/
SAAM IIWashington Universityhttps://tegvirginia.com/software/saam-ii/
WinNonlinCertara USA Inc.https://www.certara.com/software/phoenix-winnonlin/
USER-FRIENDLY SOFTWARE
Cloe® PKCyprotexhttps://www.cyprotex.com/insilico/
GastroPlusSimulations Plus Inc.https://www.simulations-plus.com/
Medici-PKComputing in Technologyhttp://www.cit-wulkow.de/
PK—SimBayer Technology Serviceshttp://www.systems-biology.com/sb/
SimCYPCertara USA Inc.https://www.certara.com/software/simcyp-pbpk/

Two types of several PBPK software programs.

Table 7

User customized softwareLatest versionModeling languageAuxiliary toolOperation methodModel
WinNonlinWinNonlin v6.4RPhoenix® NLME™Built-in options or a combination of graphics and text commandsNon-compartmental analysis (NCA)
NONMEMNONMEM v7.3FORTRANWings for NONMEM, prianaFully based on textTwo-compartment or three-compartment model
User-friendly softwareCompanySpeciesRoutes of administrationFeaturesApplication
GastroPlusSimulations Plus Inc.Human, rat, dog, mouse, monkey, user definedi.v., p.o., ocular, pulmonary, lingual, sublingual, buccalACAT model PEAR functionTen modules for PBPK modeling, built-in multi-person physiological treatment model, IVIVC; single simulation, batch simulation, parameter sensitivity analysis (PSA), population simulation
SimCYPCertara USA Inc.Human, rat, mouse, dogi.v., p.o., pulmonary, skinADAM model, Extensive data base with physiological informationPBPK modeling function, IVIVE; DDI simulation and population differential prediction; time course of simulated metabolites
PK-SimBayer Technology ServicesHuman, rat, mouse, dog, monkey, minipigi.v., p.o., subcutaneous, dermal, user definedWB-PBPK model, various PBPK calculation methodsHigh-flexibility-all model parameters accessible for specific investigations; Simulating physiological variability in reactions

Comparison of commonly used PBPK software programs for two categories.

It was initially used in engineering and mathematics of the user customized software because the essence of PBPK modeling is mathematical modeling (Bouzom et al., 2012). Hence, the application of these software programs in PBPK modeling is a natural transformation. The user customized software requires users to write their own programs to build the PK model at the beginning stage of development. This procedure requires users to have proficient programming skills as well as expertise in the field. Recently, specific PK or PBPK modules and equation libraries, as well as visual graphical interfaces, have been added to some of these software programs (Bouzom et al., 2012). By using these software programs, users can quickly generate standard PBPK models by following existing templates that already contain standard codes and equations, greatly facilitating user operation.

The user-friendly type, customized for PBPK modeling, has a graphical interface. It requires no modeling language and programming, so it is relatively simple to operate. Originally, some such software programs were specifically modeled for predicting a specific property of the ADME process, such as absorption (GastroPlus, Tian et al., 2011) or metabolism (SimCYP, Jamei et al., 2009). These software programs gradually evolved into complex PBPK modeling tools for the entire body. Recently, the function of updated versions has become increasingly sophisticated. Now, they not only model their specific areas but also simulate the whole-body pharmacokinetic process, which is absorption, metabolism, and excretion, etc. (Li M. et al., 2017; Byun et al., 2020). These software programs can perform various tasks, such as simulation, parameter evaluation, and sensitivity analysis, simply by inputting specific drug parameters and choosing certain model options (Bouzom et al., 2012).

Applications

The in silico applications of predicting ADMET profiles in 2016-2018 were collected by searching PubMed. We analyzed the search results and briefly introduce how software programs and methods predict the properties of ADMETs.

Molecular Modeling

Most applications of molecular modeling focus on predicting the strength of the interaction between a molecule and a metabolic enzyme or transporter. For instance, Niu et al. (2016) performed docking studies of flavokawain A (FKA) and its target CYP450. FKA shows obvious inhibition of different CYP isoforms, and subsequent inhibition experiments showed that CYP3A2 was the primary isoform contributing to the metabolism of FKA. Gong et al. (2018) performed a molecular docking experiment to study the binding mode between sauchinone and CYPs. The results showed the interactions of sauchinone in the active site of CYP2B6, 2C19, 2E1, and 3A4. In addition to the above examples, the details of 22 representative studies are included in Table 8.

Table 8

YearCompoundMetabolizing enzymeToolReferences
20164 Steroid derivativesCYP1B1 (3PM0)Gold-5.4Poirier et al., 2016
20168 omeprazole-based analogsCYP2C19 (4GQS)GOLD 5.2Li et al., 2016
2016Flavokawain ACYP2D6 (3QM4)SYBYL-X2.0Niu et al., 2016
2016Ketoconazole, Resveratrol, MAR, DAR and TARCYP3A4 (2V0M)Discovery Studio 4.0Basheer et al., 2016
2016Naringenin, 6',7'-dihydrox-ybergamottinCYP1A1 (4I8V)AutodockSantes-Palacios et al., 2016
2016Pyridine, Piperidine and Azol scaffoldsUBE2D4Discovery Studio 3.5Ramatenki et al., 2016
2016Rab38 inhibitorsRab38Discovery Studio 4.0Abdelmonsef et al., 2016
20171′-S-1′-acetoxychavicol acetateCYP1A2 (2HIF), CYP2D6 (3QM4), CYP3A4 (4D6Z)AutoDock Vina version 1.1.2Haque et al., 2017
201712 estrone (E1), 17 β-estradiol (E2) derivativesCYP1B1(3PM0)GOLD 5.4Dutour et al., 2017
2017Progesterone (PGS)CYP3A4 (1W0F)Schrödinger Suite, 2012Du H. et al., 2017
2017Quinoxaline, Diazepine, PiperazineUBE2NL enzymeLigPrep version 5.6, SchrödingerRamatenki et al., 2017
2017Resveratrol, Nitrostilbene, Dimethoxy-nitrostilbene, KetoconazoleCYP3A4(2V0M)Discovery Studio 4.0 (CDOCKER), Schrodinger Suite 2016 platform (Glide docking)Basheer et al., 2017
2017Sulfonyl hydrazonesMAO-A and BAutoDock VinaAbid et al., 2017
2017Wilfortrine, Wilforine, Wilfordine, Euonymine, WilforgineCYP3A4 (1W0F)Discovery Studio (CDOCKER)Wang L. et al., 2017
2017XIAPCaspase-3AutoDock VinaProkop et al., 2017
2017α-Naphthoflavone (ANF) 7-ethoxyresorufin (7ER)CYP1A2 (2HI4)GOLD 5.2.2Watanabe et al., 2017
201815 Vinca derivativesCYP3A4 (3NXU), CYP3A5GOLD 5.2Saba and Seal, 2018
20183 CompoundsCYP3A4 (4NY4)Schrödinger Release 201702Vaz et al., 2018
2018Bavachina, Neobavaisoflavoneb, Corylifol AcCYP1A2 (2HI4)a, b, c, CYP2C9 (1R9O)a, b, c, CYP2C19 (4GQS)a, b, c, CYP2D6 (3TGB)a, b, c, CYP3A4 (1WOF)a, b, c, CYP2E1 (3T3Z)a, b, cDiscovery Studio 4.1 (CDOCKER)Wang L. et al., 2018
2018FAK and Triazinic inhibitorsFAKSchrödinger 9.0Cheng P. et al., 2018
2018Metconazole (MEZ) isomersCYP3A4 (2V0M)AutoDock VinaZhuang S. et al., 2018
2018Paracetamola, PilocarpinebCYP2E1 (3T3Z)a,b, CYP3A4aSYBYL-X 1.3Wang Y. et al., 2018
2018SauchinoneCYP3A4 (3UA1), CYP2B6 (3IBD), CYP2C19 (4GQS), CYP2E1 (3GPH)AutoDock 4.2.6Gong et al., 2018
2018Sulfaphenazole, Chondroitin disaccharide Δdi-4S (C4S), Glucosamine 3-sulfate, Glucosamine 6-sulfate, Diacerein, RheinCYP2C9 (1R9O)Discovery Studio 4.0 (CDOCKER)Tan et al., 2018

Applications of molecular modeling in predicting metabolic properties from 2016 to 2018.

Superscript values denote that the metabolizing enzymes in the second column correspond to the query molecules in the first column, respectively.

QSAR

Prediction of the pharmacokinetic properties of compounds using QSAR relies mainly on traditional models or software developed based on constructed data sets. Table 9 lists seven typical applications of the QSAR method. For example, Khan et al. (2016) utilized the QSAR model in ACD/I-lab to determine multiple ADMET properties (such as logS, logP, logD, BBB) for 6 compounds targeting heat-shock protein 90 (Hsp90). Then, six compounds were designed according to BBB and antiangiogenic properties. One molecule (compound 6) was observed to inhibit Hsp90 with a predicted efficiency of BBB permeation of 0.55 kcal/mol in comparison to the experimental value of 0.625 kcal/mol. Ajay Kumar et al. (2018) performed 3D-QSAR studies to filter compounds based on ADME properties by using Schrodinger. Fifty hit compounds targeted to transforming growth factor-β (TGF-β) type I were screened based on predicted ADME properties (such as BBB, logS, and Lipinski's rule of five). Seven molecules were finally selected as the lead compound for subsequent research.

Table 9

YearTool/MethodCompoundPropertiesReferences
2016ACD/I-labAngiogenic inhibitor for brain tumorlogS, logP, logD, logBB, hERG inhibition, HBA, HBD, MWKhan et al., 2016
2016QikProp 4.6N-pyridyl, Pyrimidine benzamideslogS, logBB, MDCK, logKP, metab, CNS, loghERG, HOAMalik et al., 2016
2017ACD ChemSketchFlavonoidsMW, IC50, Index of refraction, Surface tension, Density, Polarity, logPDas et al., 2017
2017CoMFA modelCITCO, α-naphtholphthalein, diethylstilbestrol, TPP, phenytoin, (R)-ethotoin, (S)-ethotoinDrug-drug interactions, logPKato et al., 2017
2017Discovery Studio v3.5Novel dibenzofuran derivativesPSA, Solubility, HIA, Cytochrome P450 2D6, BBB, PPB, HepatotoxicityMa et al., 2017
2017Self-organizing molecular field analysis (SOMFA)Curcumin analogsMW, IC50, logPVerma and Thareja, 2017
2018CoMFA and CoMSIAScopoletin Phenolic Ether DerivativesLC50, TPSALuo et al., 2018
2018MLR, SYBYLX v1.3, CoMFA and CoMSIAAmyloid β aggregation inhibitorsIC50Aswathy et al., 2018
2018Quantum Mechanics/ Molecular Mechanics (QM/MM)DCHA inhibitorsIC50Kollar and Frecer, 2018
2018QikProp 4.6TGF-type I inhibitorsIC50, BBB, logSAjay Kumar et al., 2018

Applications of QSAR in predicting ADMET properties from 2016 to 2018.

LogP, octanol-water partition coefficient; LogD, octanol-water distribution coefficient; logS, aqueous solubility; logBB, blood-brain barrier; hERG, human ether-a-go-go related gene block; HOA, human oral absorption; PPB, plasma protein binding; MW, Molecular weight; LC50, median lethal concentration; TPSA, topological polar surface area; PSA, polar surface area; IC50, half maximal inhibitory concentration; HIA, human intestinal absorption; MDCK, madin Darby canine kidney cell line.

PBPK

PBPK modeling has been consistently performed to predict pharmacokinetics with the help of some widely used software programs (Table 7). Seventy-four applications using three software packages are listed in Table 10. WinNonlin is one of the most widely used software programs in the prediction of pharmacokinetics. For example, Gestrich et al. (2018) used WinNonlin (v6.4) to analyze compartmental and non-compartmental gentamicin plasma concentrations vs. time. The peak drug concentrations and AUCs in young adults and older alpacas were compared, and both were significantly lower in young adults than in geriatric alpacas. The increased drug exposure and decreased clearance in geriatric alpacas created a greater risk of ADRs and/or therapeutic failure. Another software, NONMEM, is also widely used to predict the impact of drugs on the target population. It is the “gold standard” software package for analysis of population PK/PD data. For example, Polepally et al. (2018) used non-linear mixed-effects modeling in NONMEM (version 7.3) to analyze concentration-time data to estimate the effect of age on the pharmacokinetic parameters of lamotrigine (LTG). By comparing the pharmacokinetic characteristics of young adult and elderly epilepsy patients, it was concluded that the bioavailability of LTG was not affected by age (Polepally et al., 2018). However, LTG CL in the elderly was 27.2% lower than in young epilepsy patients. These findings are very useful for clinicians to offer optimal epilepsy care and support to elderly patients starting low-dose treatment (Polepally et al., 2018).

Table 10

YearToolCompound /PreparationsPropertiesReferences
2016GastroPlusAlectinibOral bioavailability, Aqueous solubility, Cmax, AUCParrott et al., 2016
2016GastroPlusBisoprolol, Nifedipine, Cimetidine, FurosemideOral absorptionHansmann et al., 2016
2016GastroPlusCarvedilol loaded nanocapsules (CLN)AUC, CmaxGeorge et al., 2016
2016GastroPlusKetoconazole, ErythromycinAUC, CmaxBoetsch et al., 2016
2016GastroPlusLevofloxacinlogP, Plasma protein binding, AUC, CmaxZhu et al., 2016
2016GastroPlusMet XR 1000 mg tabletsDose, logD, AUC, CmaxChen W. et al., 2016
2016PKSimAzathioprineHepatotoxicityThiel et al., 2016
2016PKSimMorphine and FurosemideVss valueSchlender et al., 2016
2016WinNonlin6-chloro-9-nitro-5-oxo-5H-benzo-(a)-phenoxazine (CNOB)AUC, Vd, CL, half-lifeWang J. H. et al., 2016
2016WinNonlinArbekacinVdHagihara et al., 2016
2016WinNonlinBusulfan (BU)AUC, apparent clearancede Castro et al., 2016
2016WinNonlinCetuximab (CTX), Capecitabine (CCB)Cmax, AUC, VdRachar et al., 2016
2016WinNonlinGlibenclamideCL, Vd, half-lifeRambiritch et al., 2016
2016WinNonlinImatinib mesylateBioavailability, Tmax, CmaxArora et al., 2016
2016WinNonlinMidazolamCmax, VdVuu et al., 2016
2016WinNonlinMidazolam, IrinotecanClearance rateLee et al., 2016
2016WinNonlinPerfluorooctanoic acid (PFOA), Perfluorooctanesulfonic acid (PFOS), Perfluorohexane sulfonic acid (PFHxS)AUC, renal clearance, Cmax, VdKim et al., 2016
2016WinNonlinPropofolPlasma concentrations, VdChen J. Y. et al., 2016
2016WinNonlinSodium succinate, Polysorbate, Arginine, Phosphate-buffered saline (PBS)CL, Half-life, CmaxGupta et al., 2016
2016WinNonlinTacrine hydrochlorideSkin penetrationPatel et al., 2016
2016WinNonlinTilmicosinAUC, CmaxZhang et al., 2016
2016WinNonlinTreosulfan (TREO)Cmax, AUCRomanski et al., 2016
2017WinNonlinTacrolimusPlasma clearanceDavid-Neto et al., 2017
2017GastroPlusBasmisanillogD, SolubilityYang et al., 2017
2017GastroPlusBuagafuranPlasma Protein Binding, logPYang et al., 2017
2017GastroPlusCompound A (CPD A)Bioavailability, Cmax, AUCStillhart et al., 2017
2017GastroPlusMangiferinAqueous solubility, logD, logP, PermeabilityKhurana et al., 2017
2017PKSimCefazolin, Cefuroxime, CefradineCLDallmann et al., 2017
2017PKSimEndogenous IgGCL, VdNiederalt et al., 2017
2017PKSimFentanyl, Alfentanil, Thiopental, Omadacycline, Amiodarone, PropylthiouracilPlasma concentrationPilari et al., 2017
2017PKSimVorinostatDoseMoj et al., 2017
2017WinNonlinAcetylkitasamycinCmax, Tmax, AUCNan et al., 2017
2017WinNonlinBenznidazoleCmax, AUCMolina et al., 2017
2017WinNonlinCeftiofurCL, VdWang J. et al., 2017
2017WinNonlinCloxacillinDoseBurmanczuk et al., 2017
2017WinNonlinDanofloxacinAUCZhang N. et al., 2017
2017WinNonlinDiaveridineOral bioavailability, Cmax, AUCLi Y. F. et al., 2017
2017WinNonlinEnrofloxacinPlasma concentrationShan et al., 2017
2017WinNonlinIohexolCLZhang C. et al., 2017
2017WinNonlinLurasidoneCmax, AUCHu et al., 2017
2017WinNonlinMeropenemPlasma concentrationKong et al., 2017
2017WinNonlinMetolazoneAUC, Cmax, TmaxLi X. et al., 2017
2017WinNonlinPsilocybinNephrotoxicityBrown et al., 2017
2017WinNonlinPyrazinamideCmaxVoelkner et al., 2017
2017WinNonlinSarafloxacinDoseYu et al., 2017
2017WinNonlinSildenafilCL, Tmax, VdOlguin et al., 2017
2017WinNonlinTenofovirCL, Plasma concentrationDu X. et al., 2017
2017WinNonlinTilmicosinAUC, Tmax, half-lifeZhang L. et al., 2017
2017WinNonlinTreosulfanLiver, Brain toxicityRomanski et al., 2017
2017WinNonlinTulathromycinDoseZhou et al., 2017
2018GastroPlusCefadroxilPermeability, logP, Aqueous solubility, Distribution volumeHu and Smith, 2018
2018GastroPlusCompound-ABioavailability, logP, Permeability, Aqueous solubilityKou et al., 2018
2018GastroPlusDasatinibAUC, Cmax, TmaxVaidhyanathan et al., 2018
2018GastroPlusDPP-4 inhibitorsBioavailabilityDaga et al., 2018
2018GastroPlusLanabecestat (AZD3293)BioavailabilityYe et al., 2018
2018GastroPlusTramadolRenal clearanceT'Jollyn et al., 2018
2018PKSimCarvedilolPlasma concentrationIbarra et al., 2018
2018PKSimEscitalopramPlasma exposureDelaney et al., 2018
2018PKSimIndomethacin, FelodipinePlasma exposureKeemink et al., 2018
2018PKSimPregabalinPlasma concentrationsIdkaidek et al., 2018b
2018WinNonlinCefquinomeCmax, AUCShan et al., 2018
2018WinNonlinDexmedetomidineCL, VdSong et al., 2018
2018WinNonlinGentamicinAUCGestrich et al., 2018
2018WinNonlinGliclazideVd, CLShaik et al., 2018
2018WinNonlinLetrozoleDoseArora et al., 2018
2018WinNonlinMoxidectinCL, Tmax, CmaxXiao et al., 2018
2018WinNonlinPenicillin GVd, CL, half-lifePadari et al., 2018
2018WinNonlinSeroquel XR, Quesero XRAUC, Cmax, CLHuang et al., 2018
2018WinNonlinSitagliptinCmax, Tmax, AUC, T1/2Sangle et al., 2018
2018WinNonlinVitacoxibPlasma concentrations, Cmax, AUCWang et al., 2018a
2018WinNonlinAmoxicillin/clavulanic acid tabletsCmax, Tmax, AUC, T1/2De Velde et al., 2018
2018WinNonlinPhenylbutyric acid (PBA), phenylacetic acid (PAA), and phenylacetylglutamine (PAGN), UPAGNPlasma concentrations, AUC, TmaxBerry et al., 2018
2018WinNonlinSiponimodT1/2Jin et al., 2018
2018WinNonlinLacosamideT1/2, Vd, Cmax, CL, AUCFranquiz et al., 2018
2018WinNonlinTildipirosin (TD)Cmax, AUC, TmaxWang et al., 2018b
2018WinNonlinStreptomycinT1/2, AUC, Tmax, CmaxChen and Gao, 2018
2018WinNonlinPantoprazoleT1/2, AUC, Vd/F, CL/FShakhnovich et al., 2018
2018WinNonlinFVIII concentratesT1/2Cheng X. et al., 2018
2018WinNonlinAdherence to tenofovir disoproxil fumarate/emtricitabine (TDF/FTC)Cmax, AUCIbrahim et al., 2018
2018WinNonlinValsartan, Hydrochlorothiazide (HCT)Cmax, Tmax, AUC, T1/2, KelIdkaidek et al., 2018a
2018WinNonlinIndividual total serum cortisol, unbound serum cortisol and salivary cortisoneCmax, AUC, Tmax, BioavailabilityJohnson et al., 2018

Applications of PBPK in predicting ADMET properties from 2016 to 2018.

LogP, octanol-water partition coefficient; LogD, octanol-water distribution coefficient; AUC, area under the curve; Cmax, peak concentration; Vss, steady-state; Vd, volume of distribution; T1/2, half-life; logP, octanol-water partition coefficient; logD, octanol-water distribution coefficient; CL, clearance; CL/F, body clearance corrected for bioavailability; Vd/F, volume of distribution corrected for bioavailability. Kel, elimination rate constant; F, bioavailability; Tmax, time to maximum concentration.

GastroPlus and SimCYP are the mainstream PBPK emulation software programs. In recent years, applications have shown a growing trend in the use of these software programs along with the improvement of software functions. For example, Ye et al. (2018) constructed an absorption model using GastroPlus to predict the potential effects of different gastric pH levels on the pharmacokinetics of lanabecestat and found that changes in gastric pH had a minimal influence on clinical exposure to lanabecestat. They also compared the bioavailability of two tablet formulations and an oral solution. The results showed that the 90% confidence intervals for geometric mean ratios were within standard bioequivalence boundaries for all other pharmacokinetic parameters, indicating that both tablet formulations were located in the accepted bioequivalence criteria compared with the oral solution (Ye et al., 2018). Boland et al. (2018) used SIMCYP to generate a dose-concentration model by using data from different genders, ages, and oral morphine formulations. The model was then validated against clinical pharmacokinetics data and used to calculate the association of the morphine dose with the plasma concentration. Finally, the analysis showed that older age, female sex, modified-release formulation, and inferior renal function were related to higher plasma concentrations (Boland et al., 2018). This result can help clinicians provide personalized prescription decisions.

ADMET

We listed ADMET prediction applications by using AdmetSAR, SwissADME, Lazar, or Toxtree in Table 11. For example, Petrescu et al. (2019) used the AdmetSAR computational program to study the cytotoxicity of 15 phenolic compounds. The results showed that these compounds were much less toxic to aquatic life than synthetic pesticides. Roman et al. (2018) used SwissADME, FAFDrugs4, and admetSAR to predict the ADMET profiles and pharmacokinetics of 31 anabolic and androgen steroids in humans. The results revealed that the investigated steroids showed high gastrointestinal absorption and good oral bioavailability, which may be useful in the inhibition of human cytochromes associated with the metabolism of xenobiotics. In addition, the side effects of the studied steroids in humans were also predicted. Silva et al. (2019) predicted the theoretical toxicity of fluconazole (FNZ) by using Lazar to study the toxicity profile of FNZ toward human peripheral blood mononuclear cells (PBMCs) cultured in vitro. The results showed that FNZ had potential mutagenic, tumorigenic, stimulating, and carcinogenic effects (Silva et al., 2019). Zhuang J. et al. (2018) used Toxtree (version 2.6.13) to evaluate the toxicity of extractables from multilayer coextrusion bags, and their prediction results revealed the types of extractables as well as the bioaccumulation factor and mutagenicity.

Table 11

YearToolCompound/PreparationsPropertiesReferences
2016admetSARGSK-3 targeting ligandsMutagenicity, BBB, HIA, Caco-2, MDCK, PPB, AMES test, carcinogenicity, rat acute toxicity, P-gp substrate/inhibitor probabilityNisha et al., 2016
2016LazarA dataset of air toxins (332 chemicals), A subset of the gold carcinogenic potency database (480 chemicals)Carcinogenicity, mutagenicityPradeep et al., 2016
2016LazarThe Schiff bases of Benzothiazol-2-ylamine, Thiazolo [5, 4-b] pyridin-2-ylamineMax, daily dose, acute toxicity, LC50Shukla et al., 2016
2016ToxTreeCinnamaldehyde, EugenolBiodegradability, genotoxicity, carcinogenicity, bioaccumulation, developmental toxicity, mutagenicity, LD50, LC50Absalan et al., 2016
2016ToxTreePC-replacement products—the 48 substances chosen based on the publications of (Simoneau et al., 2012) and (Onghena et al., 2014, 2015)Genotoxic carcinogenicityMertens et al., 2016
2016ToxTreeVaccine constituentsCarcinogenicity, mutagenicity, genotoxicity, LD50White et al., 2016
2017admetSAR3-bromopyruvate, Dibromopyruvate (DBPA), Propionic acid (PA)BBB, Human intestinal absorption, Caco-2 permeability, P-glp substrate, AMES toxicity, Acute oral toxicity, Acute toxicity, CYP450 substrate and inhibitor, hERGYadav et al., 2017
2017admetSARHistone deacetylase (HDACs) inhibitorsLogS, Caco-2 permeabilityUba and Yelekci, 2017
2017admetSARVernonia anthelmintica (L.)Bioavailability, HIA, Caco-2, metabolism CYPWang J. Y. et al., 2017
2017Lazar2-amino6-methylpyridine, 6-heptenoic acid, 2-methylphenolCarcinogenicity, mutagenicityFrenzel et al., 2017
2017Lazar, ToxTree8 volatile organic compounds (VOC)Carcinogenicity, mutagenicityGuerra et al., 2017
2017SwissADME3,7-dimethyl-2,6-octadienal, 2-pentene-2-methylPhysicochemical properties, lipophilicity, hydrophilicitySimhadri Vsdna et al., 2017
2017SwissADME338 different chemical pesticidesLipophilicity, TPSA, molar refractivity, BBB permeant, GI absorptionChedik et al., 2017
2017SwissADMEGingerGI absorption, BBB, skin permeability, P-gp substrateSanni and Fatoki, 2017
2017SwissADMEPolyphenolsLipophilicity, water solubilityYugandhar et al., 2017
2017SwissADMETributyltin (IV) complex carboxylic acid derivativeGI absorption, BBB, LogP, water solubility, GI, Caco-2 Cells, Ames TestWaseem et al., 2017
2017SwissADMEXeronineLipophilicity, GI absorption, solubility, bioavailabilitySanni et al., 2017
2017ToxTree400 compoundsskin/eye irritation, corrosionVerheyen et al., 2017
2017ToxTree80 commercially available chemicals (38 liquids and 42 solids)Eye irritation, corrosionGeerts et al., 2017
2018Lazar(-)-Asimilobine, Aloin, Annoretine, Chrysothrone, Coptisine, Elymoclavine, ThalicminineGenotoxicity, carcinogenicity, mutagenicityGluck et al., 2018
2018LazarNewly proposed heterocyclic derivativesCarcinogenicity, mutagenicityAzad et al., 2018
2018LazarThe synthesized 1,3,5-trisubstituted−2-pyrazoline derivatives (5a-5t)Maximum recommended daily dose, reproductive toxicity, carcinogenicity, Mutagenicity, Acute toxicity, LC50Tripathi et al., 2018
2018SwissADME21 Organosilicone compoundsTPSA, logP, GI absorption, BBBShaaban et al., 2018
2018SwissADMEMycotoxins (DON,3-AcDON,15-AcDON)HIA, BBB penetration, Mutagenicity, Carcinogenicity, Acute toxicityTaroncher et al., 2018
2018SwissADMENAZ2329Lipophilicity, Water Solubility, GI absorption, BBB permeantAgoni et al., 2018
2018SwissADME, admetSAR, ToxTree31 anabolic hormones and androgen hormonesGastrointestinal absorption, Blood brain barrier, P-gp substrate, Skin permeability, Carcinogenicity, hERG, Ames toxicityRoman et al., 2018
2018ToxTree48 selected sensitizing and non-sensitizing ASSkin sensitizationBraeuning et al., 2018
2018ToxTreeBis [2,4-bis(2-methyl-2-propanyl) phenyl] hydrogen phosphateLD50, Bioaccumulation factor, MutagenicityZhuang J. et al., 2018
2018ToxTreePrinted paper and board FCM substancesAmes/bacterial mutagenicityVan Bossuyt et al., 2018
2018SwissADME, ToxTreeDominant phytochemicals from Rheum palmatum, Rubus coreanus and Sanguisorba officinalislogP, logD, TPSA, log S, GI absorption, CYP450 isoforms inhibitor probability, GenotoxicityNosrati et al., 2018
2018SwissADME8 mPGES-1 bindersPAINSLauro et al., 2018
2018SwissADMEA new series of synthesized quinazoline derivativesMolecular weight, logP, HBA, HBD, TPSA, Lipinski's RO5, LeadlikenessNasab et al., 2018
2018SwissADME107 Compounds containing biaryl scaffoldMolecular weight, logP, HBA, HBD, Lipinski's RO5, TPSA, BBB, GI absorptionKhalid et al., 2018
2018SwissADMEGenistein, Daidzein and GlyciteinLipinski's RO5, TPSA, Num. rotatable bondsShaji, 2018
2018SwissADMEFour series of diphyllin-related compoundsLogP, PAINSLindstrom et al., 2018

Applications of predicting ADMET properties using software from 2016 to 2018.

BBB, blood brain barrier; P-gp, P glycoprotein; HIA, human intestinal absorption; MDCK, madin Darby canine kidney cell line; PPB, plasma protein binding; GI absorption, gastrointestinal absorption; TPSA, topological polar surface area; LogP, octanol-water partition coefficient; LogD, octanol-water distribution coefficient; logS, aqeous solubility; F, Bioavailability; hERG, human ether-a-go-go related gene block; PAINS, pan-assay interference structures; LC50, median lethal concentration; LD50, median lethal dose; RO5, rule of five.

Discussion

Deficiencies in Current in silico Methods

Each in silico method has its own characteristics and application scope. Hence, we need to select the most proper method for more accurate prediction. However, some methods have obvious deficiencies that may affect the prediction results. For example, molecular modeling plays a major role only in predicting metabolism and can assess only the possible interactions between compounds and metabolic enzymes; it cannot explicitly evaluate the ADMET risks of candidate compounds. The scoring function in molecular docking also affects the accuracy of ADMET prediction (de Graaf et al., 2006; Zhou et al., 2006). For instance, Kemp et al. (2004) and de Graaf et al. (2006) used different scoring functions to evaluate the binding affinities between cytochrome P450 and its inhibitors or substrates, respectively. Kemp et al. (2004) docked 33 compounds to P450, and the results revealed a correlation coefficient of R2 = 0.61 between the docking scores and active compounds. The docking scores were only able to identify several compounds as CYP2D6 inhibitors. Although de Graaf et al. (2006) integrated six scoring functions to identify the substrates of P450, the highest predicted accuracy (GOLD-Chemscore) identified 60% known substrates in the top 5% results during the virtual screening. Only high-affinity CYP2D6 ligands could be predicted. Therefore, docking methods with scoring functions are mostly applied for coarse screening of a series of compounds.

Compared with molecular modeling, data modeling can predict more properties, but its prediction accuracy depends on the quality and quantity of data. The QSAR method, the main strategy in data modeling, has limited value without an estimated model applicability domain for predicting biological or physicochemical properties (Sushko et al., 2010). The predictive ability of models will be limited if the predicted chemical is outside the chemical space where the models were developed (Sheridan et al., 2004). Furthermore, the descriptors used in the model construction for structural transformation are too simplistic and inadequate to predict the behavior of a drug in a whole organism. Therefore, there is a need to develop molecular descriptors containing more information. Many developed QSAR models have been validated only by internal validation without sufficient external validation, which is considered a necessary factor to build a reliable QSAR model (Roy et al., 2012). QSAR prediction is based on the principle that similar molecules have similar properties (Patterson et al., 1996), but in some special cases, such as CYP metabolism, similar molecules may have different activities, which are known as activity cliffs (Guha, 2011). Therefore, the ADMET properties of the compounds in the human body are not independent but are also affected by other factors. The PBPK method can predict multiple properties, but only provides common information about the biological behavior of organs or tissues. It is also limited by the mathematical form of the PBPK equations, which ignores the structural and physical properties of drug compounds (Huynh et al., 2009). Moreover, a large amount of experimental data is required when constructing models. Due to the lack of proper and easily accessible databases related to physiological properties, the data used to build models can be obtained only from the literature (Rowland et al., 2004). However, the obtained data are relatively limited, which reduces the predictability of the models. Finally, although some PBPK software programs have been developed, most of them are commercial, and users must participate in the training of software companies to make them more useful (Lave et al., 2007).

The existing ADMET software can perform faster and more convenient predictions of multiple properties to obtain more comprehensive prediction results. However, we can see in practical applications that the software is applicable only for qualitative analysis of compounds and cannot accurately predict the quantitative values of some properties. Moreover, issues are observed in the data quality and quantity of these software. While more experimental data are needed to further optimize the software, integrating unconfirmed data into the software to predict new compounds will decrease the prediction accuracy.

Comparison of the Applications of Three Methods in Predicting ADMET Properties

To compare the application trends of each method more intuitively, we counted the number of applications of these three methods from 2016 to 2018 (Figure 6). Applications using PBPK modeling software have exhibited an upward trend in recent years, which means that this approach will be the mainstream pharmacokinetic evaluation in the future. We believe that more user-friendly software will be developed, making the prediction process more convenient for users. Applications using ADMET software for forecasting are also on the rise. The main reason for these changes may be the demand for multi-property prediction and the phenomenon of drug recall, leading to the hope of predicting more pharmacokinetic properties in advance to reduce drug development costs. This increase also confirms that researchers have gradually integrated software predictions into the early stages of drug discovery to improve the success rate of drug development through complementarity and collaboration. In addition, we found that the QSAR method is rarely used alone, potentially due to the limitations of the QSAR method prediction. Therefore, we can increase or integrate the data volume of some specific QSAR models to generate a more comprehensive prediction model. Researchers can also combine different modeling methods and then perform relevant predictions. If a combined model can predict the properties well, then it can be used as a consensus approach to improve ADMET prediction accuracy.

Figure 6

Previous Review and Prospective Studies of ADMET and PBPK Simulations

To date, there have been many reviews on ADMET and PBPK, which provide summaries from different perspectives. We found 11 articles that are closely related to this review after careful retrieval of the published literature. We have classified these articles into three categories based on the perspective of the descriptions: (1) directly related to our review, (2) machine learning (ML) methods, and (3) advances in in silico ADMET modeling.

Five reviews are directly related to this article. Yamashita and Hashida (2004) reviewed the application of structure-based methods and QSAR methods in predicting ADMET in detail. Cheng et al. (2013) introduced the recent progress and current challenges of QSAR in predicting ADMET and then discussed several new promising research directions that could be employed for systemic in silico ADMET prediction. Alqahtani (2017) reviewed the in silico models for predicting the ADMET properties of compounds and provided a comprehensive overview of the latest modeling methods and algorithms, as well as the application prospects of PBPK in predicting pharmacokinetics. Wang et al. (2015) briefly introduced the development of in silico models for ADMET prediction. They also focused on the modeling approaches, related applications, and potential advantages or disadvantages of these models used in drug discovery. Wishart (2007) highlighted ADMET property prediction, as well as ADME-related databases and software, and briefly introduced the application of PBPK and related software in ADMET.

Three out of 11 reviews are about machine learning methods for predicting ADMET. Ferreira and Andricopulo (2019) provided a detailed description of current machine learning approaches to ADMET modeling, focusing on key advances from 2017 to 2018. Tao et al. (2015) reviewed the progress of machine learning methods in ADMET prediction and discussed the performance, applications, and challenges in developing machine learning methods. Maltarollo et al. (2015) described the applications of ML methods in ADMET prediction.

The remaining three reviews addressed the progress of ADMET modeling. Specifically, Lin et al. (2003) introduced the desirable properties of new chemical entities (NCEs) from an ADMET perspective and discussed basic concepts, important tools, reagents, and experimental approaches used by researchers in predicting human pharmacokinetics. van de Waterbeemd and Gifford (2003) summarized the endpoints of pharmacokinetics, metabolism, and toxicity. Wang and Hou (2009) introduced the properties of ADMET, discussed the latest corresponding in silico models, and provided a brief summary of some software and databases. These articles introduced ADMET, PBPK, and related research progress from different perspectives.

Conclusion

In this review, we provide a detailed and comprehensive introduction to currently used approaches or tools in predicting ADMET properties, including the basic principles, classification, and applications. In addition, we collect related applications from published articles over the past 3 years and analyze the trends in these applications. The purpose of this review is to help readers quickly understand these approaches and the characteristics of the related tools (databases and software). It may also provide readers with a better understanding of how existing tools can be applied to pharmacokinetic predictions. We are convinced that more accurate predictions due to users' familiarity with existing online services will increase the importance of in silico ADMET prediction in pharmacokinetics. In the future, we expect not only a reduced failure rate in drug development and drug recalls but also a faster timeline from R&D to market, as well as decreased costs during the late stage of development.

Statements

Author contributions

ZH, FW, and YZ contributed to the design and conception of the study. LL, YZ, FW, XS, GC, XW, XL, and MT performed information retrieval and analysis. FW, YZ, LL, and XS wrote the manuscript. YZ, FW, LL, XS, GC, XW, XL, and MT created the tables and figures. ZH and FW guided the manuscript writing and revised the manuscript. ZH provided financial support. All authors contributed to manuscript revision and have read and approved the submitted version.

Funding

This work was supported by the National Natural Science Foundation of China (31770774), the Provincial Major Project of Basic or Applied Research in Natural Science, the Guangdong Provincial Education Department (2016KZDXM038), and the Higher Education Reform Project of Guangdong Province (2019268).

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

References

  • 1

    AbdelmonsefA. H.DulapalliR.DasariT.PadmaraoL. S.MukkeraT.VuruputuriU. (2016). Identification of novel antagonists for rab38 protein by homology modeling and virtual screening. Comb. Chem. High Throughput Screen.19, 875892. 10.2174/1386207319666161026153237

  • 2

    AbidS. M. A.YounusH. A.Al-RashidaM.ArshadZ.MaryumT.GilaniM. A.et al. (2017). Sulfonyl hydrazones derived from 3-formylchromone as non-selective inhibitors of MAO-A and MAO-B: synthesis, molecular modelling and in-silico ADME evaluation. Bioorg. Chem.75, 291302. 10.1016/j.bioorg.2017.10.001

  • 3

    AbsalanA.Mesbah-NaminS. A.TiraihiT.TaheriT. (2016). The effects of cinnamaldehyde and eugenol on human adipose-derived mesenchymal stem cells viability, growth and differentiation: a cheminformatics and in vitro study. Avicenna J. Phytomed.6, 643657.

  • 4

    AgoniC.RamharackP.SolimanM. E. S. (2018). Allosteric inhibition induces an open WPD-loop: a new avenue towards glioblastoma therapy. RSC Adv.8, 4018740197. 10.1039/C8RA08427K

  • 5

    Ajay KumarT. V.AthavanA. A. S.LoganathanC.SaravananK.KabilanS.ParthasarathyV. (2018). Design, 3D QSAR modeling and docking of TGF-beta type I inhibitors to target cancer. Comput. Biol. Chem.76, 232244. 10.1016/j.compbiolchem.2018.07.011

  • 6

    AlarnS.KhanF. (2019). 3D-QSAR, docking, ADME/Tox studies on flavone analogs reveal anticancer activity through Tankyrase inhibition. Sci. Rep.9:5414. 10.1038/s41598-019-41984-7

  • 7

    AlqahtaniS. (2017). In silico ADME-Tox modeling: progress and prospects. Expert Opin. Drug Metab. Toxicol.13, 11471158. 10.1080/17425255.2017.1389897

  • 8

    AndradeC. H.SilvaD. C.BragaR. C. (2014). In silico prediction of drug metabolism by P450. Curr. Drug Metab.15, 514525. 10.2174/1389200215666140908102530

  • 9

    AroraP.AdamsC. H.GudelskyG.DasGuptaB.DesaiP. B. (2018). Plasma and brain pharmacokinetics of letrozole and drug interaction studies with temozolomide in NOD-scid gamma mice and sprague dawley rats. Cancer Chemother. Pharmacol.83, 8189. 10.1007/s00280-018-3705-6

  • 10

    AroraR.SharmaM.MonifT.IyerS. (2016). A multi-centric bioequivalence trial in Ph+ chronic myeloid leukemia patients to assess bioequivalence and safety evaluation of generic imatinib mesylate 400 mg tablets. Cancer Res. Treat.48, 11201129. 10.4143/crt.2015.436

  • 11

    AswathyL.JishaR. S.MasandV. H.GajbhiyeJ. M.ShibiI. G. (2018). Design of novel amyloid beta aggregation inhibitors using QSAR, pharmacophore modeling, molecular docking and ADME prediction. In Silico Pharmacol.6:12. 10.1007/s40203-018-0049-1

  • 12

    AzadI.NasibullahM.KhanT.HassanF.AkhterY. (2018). Exploring the novel heterocyclic derivatives as lead molecules for design and development of potent anticancer agents. J. Mol. Graph Model.81, 211228. 10.1016/j.jmgm.2018.02.013

  • 13

    BasheerL.SchultzK.GuttmanY.KeremZ. (2017). In silico and in vitro inhibition of cytochrome P450 3A by synthetic stilbenoids. Food Chem.237, 895903. 10.1016/j.foodchem.2017.06.040

  • 14

    BasheerL.SchultzK.KeremZ. (2016). Inhibition of cytochrome P450 3A by acetoxylated analogues of resveratrol in in vitro and in silico models. Sci. Rep.6:31557. 10.1038/srep31557

  • 15

    BelekarV.LingineniK.GargP. (2015). Classification of breast cancer resistant protein (BCRP) inhibitors and non-inhibitors using machine learning approaches. Comb. Chem. High Throughput Screen.18, 476485. 10.2174/1386207318666150525094503

  • 16

    BerryS. A.VockleyJ.VinksA. A.DongM.DiazG. A.McCandlessS. E.et al. (2018). Pharmacokinetics of glycerol phenylbutyrate in pediatric patients 2 months to 2 years of age with urea cycle disorders. Mol. Genet. Metab.125, 251257. 10.1016/j.ymgme.2018.09.001

  • 17

    BhatiaS.SchultzT.RobertsD.ShenJ.KromidasL.Marie ApiA. (2015). Comparison of cramer classification between toxtree, the OECD QSAR toolbox and expert judgment. Regul. Toxicol. Pharmacol.71, 5262. 10.1016/j.yrtph.2014.11.005

  • 18

    BhhataraiB.WilsonD. M.ParksA. K.CarneyE. W.SpencerP. J. (2016). Evaluation of TOPKAT, toxtree, and derek nexus in silico models for ocular irritation and development of a knowledge-based framework to improve the prediction of severe irritation. Chem. Res. Toxicol.29, 81022. 10.1021/acs.chemrestox.5b00531

  • 19

    BoetschC.ParrottN.FowlerS.PoirierA.HainzlD.BankenL.et al. (2016). Effects of cytochrome P450 3A4 inhibitors-ketoconazole and erythromycin-on bitopertin pharmacokinetics and comparison with physiologically based modelling predictions. Clin. Pharmacokinet.55, 237247. 10.1007/s40262-015-0312-0

  • 20

    BolandJ. W.JohnsonM.FerreiraD.BerryD. J. (2018). In silico (computed) modelling of doses and dosing regimens associated with morphine levels above international legal driving limits. Palliat. Med.32, 12221232. 10.1177/0269216318773956

  • 21

    BouzomF.BallK.PerdaemsN.WaltherB. (2012). Physiologically based pharmacokinetic (PBPK) modelling tools: how to fit with our needs?Biopharm. Drug Dispos.33, 5571. 10.1002/bdd.1767

  • 22

    BouzomF.WaltherB. (2008). Pharmacokinetic predictions in children by using the physiologically based pharmacokinetic modelling. Fundam. Clin. Pharmacol.22, 579587. 10.1111/j.1472-8206.2008.00648.x

  • 23

    BowenJ. P.GuenerO. F. (2013). A perspective on quantum mechanics calculations in ADMET predictions. Curr. Top. Med. Chem.13, 12571272. 10.2174/15680266113139990032

  • 24

    BraeuningC.BraeuningA.MielkeH.HolzwarthA.PeiserM. (2018). Evaluation and improvement of QSAR predictions of skin sensitization for pesticides. SAR QSAR Environ. Res.29, 823846. 10.1080/1062936X.2018.1518261

  • 25

    BrownR. T.NicholasC. R.CozziN. V.GassmanM. C.CooperK. M.MullerD.et al. (2017). Pharmacokinetics of escalating doses of oral psilocybin in healthy adults. Clin. Pharmacokinet.56, 15431554. 10.1007/s40262-017-0540-6

  • 26

    BujakR.Struck-LewickaW.KaliszanM.KaliszanR.MarkuszewskiM. J. (2015). Blood-brain barrier permeability mechanisms in view of quantitative structure-activity relationships (QSAR). J. Pharm. Biomed. Anal.108, 2937. 10.1016/j.jpba.2015.01.046

  • 27

    BurmanczukA.GrabowskiT.OsypiukM.PolskaB.KowalskiC. (2017). Determination of cloxacillin residues in dairy cows after intramammary administration. J. Vet. Pharmacol. Ther.40, 552560. 10.1111/jvp.12391

  • 28

    ByunJ. H.HanD. G.ChoH. J.YoonI. S.JungI. H. (2020). Recent advances in physiologically based pharmacokinetic and pharmacodynamic models for anticancer nanomedicines. Arch. Pharm. Res.43, 8099. 10.1007/s12272-020-01209-2

  • 29

    CaldwellG. W.YanZ.TangW.DasguptaM.HastingB. (2009). ADME optimization and toxicity assessment in early- and late-phase drug discovery. Curr. Top. Med. Chem.9, 965980. 10.2174/156802609789630929

  • 30

    CaoD.WangJ.ZhouR.LiY.YuH.HouT. (2012). ADMET evaluation in drug discovery. 11. Pharmacokinetics knowledge base (PKKB): a comprehensive database of pharmacokinetic and toxic properties for drugs. J. Chem. Inf. Model.52, 11321137. 10.1021/ci300112j

  • 31

    ChedikL.Mias-LucquinD.BruyereA.FardelO. (2017). In silico prediction for intestinal absorption and brain penetration of chemical pesticides in humans. Int. J. Environ. Res. Public Health14:708. 10.3390/ijerph14070708

  • 32

    ChenJ. Y.YiM.YaoS. L.ZhangX. P. (2016). Propofol target-controlled infusion modeling in rabbits: pharmacokinetic and pharmacodynamic analysis. J. Huazhong Univ. Sci. Technol. Med. Sci.36, 428433. 10.1007/s11596-016-1604-9

  • 33

    ChenL.GaoZ. (2018). Continuous hypoxia reduces the concentration of streptomycin in the blood. BMC Infect. Dis.18:120. 10.1186/s12879-018-3027-7

  • 34

    ChenW.DesaiD.GoodD.CrisonJ.TimminsP.ParuchuriS.et al. (2016). Mathematical model-based accelerated development of extended-release metformin hydrochloride tablet formulation. AAPS PharmSciTech.17, 10071013. 10.1208/s12249-015-0423-9

  • 35

    ChenY.LiuZ. L.FuT. M.LiW.XuX. L.SunH. P. (2015). Discovery of new acetylcholinesterase inhibitors with small core structures through shape-based virtual screening. Bioorg. Med. Chem. Lett.25, 34423446. 10.1016/j.bmcl.2015.07.026

  • 36

    ChengF.LiW.LiuG.TangY. (2013). In silico ADMET prediction: recent advances, current challenges and future trends. Curr. Top. Med. Chem.13, 12731289. 10.2174/15680266113139990033

  • 37

    ChengF.LiW.ZhouY.ShenJ.WuZ.LiuG.et al. (2012). admetSAR: a comprehensive source and free tool for assessment of chemical ADMET properties. J. Chem. Inf. Model.52, 30993105. 10.1021/ci300367a

  • 38

    ChengP.LiJ.WangJ.ZhangX.ZhaiH. (2018). Investigations of FAK inhibitors: a combination of 3D-QSAR, docking, and molecular dynamics simulations studies. J. Biomol. Struct. Dyn.36, 15291549. 10.1080/07391102.2017.1329095

  • 39

    ChengT.LiQ.ZhouZ.WangY.BryantS. H. (2012). Structure-based virtual screening for drug discovery: a problem-centric review. AAPS J.14, 133141. 10.1208/s12248-012-9322-0

  • 40

    ChengX.LiP.ChenZ.ZhangN.ZhenY.ZhaoL.et al. (2018). Break-through bleeding in relation to pharmacokinetics of factor VIII in paediatric patients with severe haemophilia A. Haemophilia24, 120125. 10.1111/hae.13373

  • 41

    ChowE. C. Y.TalattofA.TsakalozouE.FanJ.ZhaoL.ZhangX. (2016). Using physiologically based pharmacokinetic (PBPK) modeling to evaluate the impact of pharmaceutical excipients on oral drug absorption: sensitivity analyses. AAPS J.18, 15001511. 10.1208/s12248-016-9964-4

  • 42

    CookD.BrownD.AlexanderR.MarchR.MorganP.SatterthwaiteG.et al. (2014). Lessons learned from the fate of AstraZeneca's drug pipeline: a five-dimensional framework. Nat. Rev. Drug Discov.13, 419431. 10.1038/nrd4309

  • 43

    CuiQ.KarplusM. (2003). Catalysis and specificity in enzymes: a study of triosephosphate isomerase and comparison with methyl glyoxal synthase. Adv. Protein Chem.66, 315372. 10.1016/S0065-3233(03)66008-0

  • 44

    CummingJ. G.DavisA. M.MuresanS.HaeberleinM.ChenH. (2013). Chemical predictive modelling to improve compound quality. Nat. Rev. Drug Discov.12, 948962. 10.1038/nrd4128

  • 45

    DagaP. R.BolgerM. B.HaworthI. S.ClarkR. D.MartinE. J. (2018). Physiologically based pharmacokinetic modeling in lead optimization. 1. Evaluation and adaptation of gastroplus to predict bioavailability of medchem series. Mol. Pharm.15, 821830. 10.1021/acs.molpharmaceut.7b00972

  • 46

    DainaA.MichielinO.ZoeteV. (2017). SwissADME: a free web tool to evaluate pharmacokinetics, drug-likeness and medicinal chemistry friendliness of small molecules. Sci. Rep.7:42717. 10.1038/srep42717

  • 47

    DallmannA.InceI.SolodenkoJ.MeyerM.WillmannS.EissingT.et al. (2017). Physiologically based pharmacokinetic modeling of renally cleared drugs in pregnant women. Clin. Pharmacokinet.56, 15251541. 10.1007/s40262-017-0538-0

  • 48

    Danishuddin KhanA. U. (2016). Descriptors and their selection methods in QSAR analysis: paradigm for drug design. Drug Discov. Today21, 12911302. 10.1016/j.drudis.2016.06.013

  • 49

    DasS.LaskarM. A.SarkerS. D.ChoudhuryM. D.ChoudhuryP. R.MitraA.et al. (2017). Prediction of anti-alzheimer's activity of flavonoids targeting acetylcholinesterase in silico. Phytochem. Anal.28, 324331. 10.1002/pca.2679

  • 50

    David-NetoE.RomanoP.Kamada TriboniA. H.RamosF.AgenaF.Almeida Rezende EbnerP.et al. (2017). Longitudinal pharmacokinetics of tacrolimus in elderly compared with younger recipients in the first 6 months after renal transplantation. Transplantation101, 13651372. 10.1097/TP.0000000000001369

  • 51

    DavisA. P.GrondinC. J.JohnsonR. J.SciakyD.McMorranR.WiegersJ.et al. (2019). The comparative toxicogenomics database: update 2019. Nucleic Acids Res.47, D948D954. 10.1093/nar/gky868

  • 52

    De BuckS. S.MackieC. E. (2007). Physiologically based approaches towards the prediction of pharmacokinetics: in vitro-in vivo extrapolation. Expert Opin. Drug Metab. Toxicol.3, 865878. 10.1517/17425255.3.6.865

  • 53

    de CastroF. A.SimoesB. P.GodoyA. L.Bertagnoli TrigoF. M.CoelhoE. B.LanchoteV. L. (2016). Use of an oral busulfan test dose in patients undergoing hematopoietic stem cell transplantation treated with or without fludarabine. J. Clin. Pharmacol.56, 15551562. 10.1002/jcph.758

  • 54

    de GraafC.OostenbrinkC.KeizersP. H.van der WijstT.JongejanA.VermeulenN. P. (2006). Catalytic site prediction and virtual screening of cytochrome P450 2D6 substrates by consideration of water and rescoring in automated docking. J. Med. Chem.49, 24172430. 10.1021/jm0508538

  • 55

    de GrootM. J. (2006). Designing better drugs: predicting cytochrome P450 metabolism. Drug Discov. Today11, 601606. 10.1016/j.drudis.2006.05.001

  • 56

    de GrootM. J.KirtonS. B.SutcliffeM. J. (2004). In silico methods for predicting ligand binding determinants of cytochromes P450. Curr. Top. Med. Chem.4, 18031824. 10.2174/1568026043387061

  • 57

    De VeldeF.De WinterB. C. M.KochB. C. P.Van GelderT.MoutonJ. W.COMBACTE-NET consortium (2018). Highly variable absorption of clavulanic acid during the day: a population pharmacokinetic analysis. J. Antimicrob. Chemother.73, 469476. 10.1093/jac/dkx376

  • 58

    DeardenJ. C. (2007). In silico prediction of ADMET properties: how far have we come?Expert Opin. Drug Metab. Toxicol.3, 635639. 10.1517/17425255.3.5.635

  • 59

    DelaneyS. R.MalikP. R. V.StefanC.EdgintonA. N.ColantonioD. A.ItoS. (2018). Predicting escitalopram exposure to breastfeeding infants: integrating analytical and in silico techniques. Clin. Pharmacokinet.57, 16031611. 10.1007/s40262-018-0657-2

  • 60

    DickinsM.van de WaterbeemdH. (2004). Simulation models for drug disposition and drug interactions. Drug Discov. Today BIOSILICO2, 3845. 10.1016/S1741-8364(04)02388-1

  • 61

    DongJ.WangN. N.YaoZ. J.ZhangL.ChengY.OuyangD.et al. (2018). ADMETlab: a platform for systematic ADMET evaluation based on a comprehensively collected ADMET database. J. Cheminform.10:29. 10.1186/s13321-018-0283-x

  • 62

    DossC. G. P.ChakrabortyC.NarayanV.KumarD. T. (2014). Computational approaches and resources in single amino acid substitutions analysis toward clinical research, in Advances in Protein Chemistry and Structural Biology, Vol. 94, ed DonevR. (San Diego, CA: Elsevier Academic Press), 365423. 10.1016/B978-0-12-800168-4.00010-X

  • 63

    DuH.LiJ.CaiY.ZhangH.LiuG.TangY.et al. (2017). Computational investigation of ligand binding to the peripheral site in CYP3A4: conformational dynamics and inhibitor discovery. J. Chem. Inf. Model.57, 616626. 10.1021/acs.jcim.7b00012

  • 64

    DuX.KouH.FuQ.LiY.ZhuZ.LiT. (2017). Steady-state pharmacokinetics of tenofovir disoproxil fumarate in human immunodeficiency virus-infected Chinese patients. Expert Rev. Clin. Pharmacol.10, 783788. 10.1080/17512433.2017.1321480

  • 65

    DutourR.Cortes-BenitezF.RoyJ.PoirierD. (2017). Structure-based design and synthesis of new estrane-pyridine derivatives as cytochrome P450 (CYP) 1B1 inhibitors. ACS Med. Chem. Lett.8, 11591164. 10.1021/acsmedchemlett.7b00265

  • 66

    EdgintonA. N.TheilF. P.SchmittW.WillmannS. (2008). Whole body physiologically-based pharmacokinetic models: their use in clinical drug development. Expert Opin. Drug Metab. Toxicol.4, 11431152. 10.1517/17425255.4.9.1143

  • 67

    El-ZahabiH. S. A.KhalifaM. M. A.GadoY. M. H.FarragA. M.ElaasserM. M.SafwatN. A.et al. (2019). New thiobarbituric acid scaffold-based small molecules: synthesis, cytotoxicity, 2D-QSAR, pharmacophore modelling and in-silico ADME screening. Eur. J. Pharm. Sci.130, 124136. 10.1016/j.ejps.2019.01.023

  • 68

    EspieP.TytgatD.Sargentini-MaierM. L.PoggesiI.WateletJ. B. (2009). Physiologically based pharmacokinetics (PBPK). Drug Metab. Rev.41, 391407. 10.1080/10837450902891360

  • 69

    EvangelistaW.WeirR. L.EllingsonS. R.HarrisJ. B.KapoorK.SmithJ. C.et al. (2016). Ensemble-based docking: from hit discovery to metabolism and toxicity predictions. Bioorg. Med. Chem.24, 49284935. 10.1016/j.bmc.2016.07.064

  • 70

    FanJ.de LannoyI. A. (2014). Pharmacokinetics. Biochem. Pharmacol.87, 93120. 10.1016/j.bcp.2013.09.007

  • 71

    FerreiraL. L. G.AndricopuloA. D. (2019). ADMET modeling approaches in drug discovery. Drug Discov. Today24, 11571165. 10.1016/j.drudis.2019.03.015

  • 72

    FranquizM. J.KalariaS. N.ArmahizerM. J.GopalakrishnanM.McCarthyP. J.BadjatiaN. (2018). Lacosamide pharmacokinetics in a critically Ill patient receiving continuous venovenous hemofiltration. Pharmacotherapy38, E17E21. 10.1002/phar.2063

  • 73

    FratevF.BenfenatiE. (2008). A combination of 3D-QSAR, docking, local-binding energy (LBE) and GRID study of the species differences in the carcinogenicity of benzene derivatives chemicals. J. Mol. Graph Model.27, 147160. 10.1016/j.jmgm.2008.04.004

  • 74

    FrenzelF.BuhrkeT.WenzelI.AndrackJ.HielscherJ.LampenA. (2017). Use of in silico models for prioritization of heat-induced food contaminants in mutagenicity and carcinogenicity testing. Arch. Toxicol.91, 31573174. 10.1007/s00204-016-1924-3

  • 75

    FuL.LiuL.YangZ. J.LiP.DingJ. J.YunY. H.et al. (2020). Systematic modeling of logD7.4 based on ensemble machine learning, group contribution, and matched molecular pair analysis. J. Chem. Inf. Model.60, 6376. 10.1021/acs.jcim.9b00718

  • 76

    GaultonA.HerseyA.NowotkaM.BentoA. P.ChambersJ.MendezD.et al. (2017). The ChEMBL database in 2017. Nucleic Acids Res.45, D945D954. 10.1093/nar/gkw1074

  • 77

    GeertsL.AdriaensE.AlepeeN.GuestR.WilloughbyJ. A.Sr.KandarovaH.et al. (2017). CON4EI: evaluation of QSAR models for hazard identification and labelling of eye irritating chemicals. Toxicol. In Vitro49, 9098. 10.1016/j.tiv.2017.09.004

  • 78

    GeorgeJ. K.SinghS. K.VermaP. R. (2016). In vivo in silico pharmacokinetic simulation studies of carvedilol-loaded nanocapsules using GastroPlus. Ther. Deliv.7, 305318. 10.4155/tde-2015-0004

  • 79

    GestrichA.BedeniceD.CeresiaM.ZaghloulI. (2018). Pharmacokinetics of intravenous gentamicin in healthy young-adult compared to aged alpacas. J. Vet. Pharmacol. Ther.41, 581587. 10.1111/jvp.12506

  • 80

    GluckJ.BuhrkeT.FrenzelF.BraeuningA.LampenA. (2018). In silico genotoxicity and carcinogenicity prediction for food-relevant secondary plant metabolites. Food Chem. Toxicol.116, 298306. 10.1016/j.fct.2018.04.024

  • 81

    GongE. C.CheaS.BalupuriA.KangN. S.ChinY. W.ChoiY. H. (2018). Enzyme kinetics and molecular docking studies on cytochrome 2B6, 2C19, 2E1, and 3A4 activities by sauchinone. Molecules23:555. 10.3390/molecules23030555

  • 82

    GuerraL. R.de SouzaA. M. T.CortesJ. A.LioneV. O. F.CastroH. C.AlvesG. G. (2017). Assessment of predictivity of volatile organic compounds carcinogenicity and mutagenicity by freeware in silico models. Regul. Toxicol. Pharmacol.91, 18. 10.1016/j.yrtph.2017.09.030

  • 83

    GuhaR. (2011). The ups and downs of structure-activity landscapes. Methods Mol. Biol.672, 101117. 10.1007/978-1-60761-839-3_3

  • 84

    GuptaP.KamathA. V.ParkS.ChiuH.LutmanJ.MaiaM.et al. (2016). Preclinical pharmacokinetics of MHAA4549A, a human monoclonal antibody to influenza A virus, and the prediction of its efficacious clinical dose for the treatment of patients hospitalized with influenza A. MAbs8, 991997. 10.1080/19420862.2016.1167294

  • 85

    HagiharaM.KatoH.HamadaY.HiraiJ.SakanashiD.SuematsuH.et al. (2016). Population pharmacokinetics of arbekacin in different infectious disease settings and evaluation of dosing regimens. J. Infect. Chemother.22, 436443. 10.1016/j.jiac.2016.03.008

  • 86

    HanY.ZhangJ. P.HuC. Q.ZhangX.MaB. F.ZhangP. P. (2019). In silico ADME and toxicity prediction of ceftazidime and its impurities. Front. Pharmacol.10:434. 10.3389/fphar.2019.00434

  • 87

    HanschC. (1981). The physicochemical approach to drug design and discovery (QSAR). Drug Dev. Res.1, 267309. 10.1002/ddr.430010403

  • 88

    HansmannS.DarwichA.MargolskeeA.AaronsL.DressmanJ. (2016). Forecasting oral absorption across biopharmaceutics classification system classes with physiologically based pharmacokinetic models. J. Pharm. Pharmacol.68, 15011515. 10.1111/jphp.12618

  • 89

    HaqueA.LeongK. H.LoY. L.AwangK.NagoorN. H. (2017). In vitro inhibitory mechanisms and molecular docking of 1'-S-1'-acetoxychavicol acetate on human cytochrome P450 enzymes. Phytomedicine31, 19. 10.1016/j.phymed.2017.05.002

  • 90

    HardyB.DouglasN.HelmaC.RautenbergM.JeliazkovaN.JeliazkovV.et al. (2010). Collaborative development of predictive toxicology applications. J. Cheminform.2:7. 10.1186/1758-2946-2-7

  • 91

    HelmaC. (2006). Lazy structure-activity relationships (lazar) for the prediction of rodent carcinogenicity and salmonella mutagenicity. Mol. Divers.10, 147158. 10.1007/s11030-005-9001-5

  • 92

    HuC.WangY.SongR.YuC.LuoX.JiaJ. (2017). Single- and multiple-dose pharmacokinetics, safety and tolerability of lurasidone in healthy Chinese subjects. Clin. Drug Investig.37, 861871. 10.1007/s40261-017-0546-8

  • 93

    HuY.SmithD. E. (2018). In silico prediction of the absorption and disposition of cefadroxil in humans using an intestinal permeability method scaled from humanized PepT1 mice. Drug Metab. Dispos.47, 173183. 10.1124/dmd.118.084236

  • 94

    HuangJ.FanX. (2011). Why QSAR fails: an empirical evaluation using conventional computational approach. Mol. Pharm.8, 600608. 10.1021/mp100423u

  • 95

    HuangX.ZhangS.MaY.YangH.HeC.TianR.et al. (2018). Bioequivalence of two quetiapine extended release tablets in Chinese healthy volunteers under fasting and fed conditions and effects of food on pharmacokinetic profiles. Drug Des. Devel. Ther.13, 255264. 10.2147/DDDT.S182965

  • 96

    HuynhL.MasereeuwR.FriedbergT.Ingelman-SundbergM.ManivetP. (2009). In silico platform for xenobiotics ADME-T pharmacological properties modeling and prediction. Part I: beyond the reduction of animal model use. Drug Discov. Today14, 401405. 10.1016/j.drudis.2009.01.009

  • 97

    IbarraM.ValianteC.SopenaP.SchiavoA.LorierM.VazquezM.et al. (2018). Integration of in vitro biorelevant dissolution and in silico PBPK model of carvedilol to predict bioequivalence of oral drug products. Eur. J. Pharm. Sci.118, 176182. 10.1016/j.ejps.2018.03.032

  • 98

    IbrahimM. E.BrooksK. M.Castillo-MancillaJ. R.McHughC.MorrowM.BrothersJ.et al. (2018). Short communication: bioequivalence of tenofovir and emtricitabine after coencapsulation with the proteus ingestible sensor. AIDS Res. Hum. Retroviruses34, 835837. 10.1089/aid.2018.0081

  • 99

    IdkaidekN.AghaH.ArafatT. (2018a). Saliva versus plasma bioequivalence of valsartan/hydrochlorothiazide in humans: validation of classes II and IV drugs of the salivary excretion classification system. Drug Res.68, 5459. 10.1055/s-0043-117775

  • 100

    IdkaidekN.HamadiS.El-AssiM.Al-ShalalfehA.Al-GhazawiA. (2018b). Saliva versus plasma therapeutic drug monitoring of pregabalin in jordanian patients. Drug Res.68, 596600. 10.1055/a-0600-2113

  • 101

    IssaN. T.WathieuH.OjoA.ByersS. W.DakshanamurthyS. (2017). Drug metabolism in preclinical drug development: a survey of the discovery process, toxicology, and computational tools. Curr. Drug Metab.18, 556565. 10.2174/1389200218666170316093301

  • 102

    JameiM.MarciniakS.FengK.BarnettA.TuckerG.Rostami-HodjeganA. (2009). The simcyp population-based ADME simulator. Expert Opin. Drug Metab. Toxicol.5, 211223. 10.1517/17425250802691074

  • 103

    JiangD.LeiT.WangZ.ShenC.CaoD.HouT. (2020). ADMET evaluation in drug discovery. 20. Prediction of breast cancer resistance protein inhibition through machine learning. J. Cheminform.12:16. 10.1186/s13321-020-00421-y

  • 104

    JinY.BorellH.GardinA.UferM.HuthF.CamenischG. (2018). In vitro studies and in silico predictions of fluconazole and CYP2C9 genetic polymorphism impact on siponimod metabolism and pharmacokinetics. Eur. J. Clin. Pharmacol.74, 455464. 10.1007/s00228-017-2404-2

  • 105

    JohnsonT. N.WhitakerM. J.KeevilB.RossR. J. (2018). Bioavailability of oral hydrocortisone corrected for binding proteins and measured by LC-MS/MS using serum cortisol and salivary cortisone. J. Bioequivalence Bioavailab.10, 001003. 10.4172/jbb.1000365

  • 106

    JonesH.ChenY.GibsonC.HeimbachT.ParrottN.PetersS.et al. (2015). Physiologically based pharmacokinetic modeling in drug discovery and development: a pharmaceutical industry perspective. Clin. Pharmacol. Ther.97, 247262. 10.1002/cpt.37

  • 107

    KarS.LeszczynskiJ. (2017). Recent advances of computational modeling for predicting drug metabolism: a perspective. Curr. Drug Metab.18, 11061122. 10.2174/1389200218666170607102104

  • 108

    KatoH.YamaotsuN.IwazakiN.OkamuraS.KumeT.HironoS. (2017). Precise prediction of activators for the human constitutive androstane receptor using structure-based three-dimensional quantitative structure-activity relationship methods. Drug Metab. Pharmacokinet.32, 179188. 10.1016/j.dmpk.2017.02.001

  • 109

    KavithaR.KarunagaranS.ChandraboseS. S.LeeK. W.MeganathanC. (2015). Pharmacophore modeling, virtual screening, molecular docking studies and density functional theory approaches to identify novel ketohexokinase (KHK) inhibitors. Biosystems138, 3952. 10.1016/j.biosystems.2015.10.005

  • 110

    KazmiS. R.JunR.YuM.-S.JungC.NaD. (2019). In silico approaches and tools for the prediction of drug metabolism and fate: a review. Comput. Biol. Med.106, 5464. 10.1016/j.compbiomed.2019.01.008

  • 111

    KeeminkJ.SjogrenE.HolmR.BergstromC. A. S. (2018). Does the intake of ethanol affect oral absorption of poorly soluble drugs?J. Pharm. Sci.108, 17651771. 10.1016/j.xphs.2018.12.006

  • 112

    KempC. A.FlanaganJ. U.van EldikA. J.MarechalJ. D.WolfC. R.RobertsG. C.et al. (2004). Validation of model of cytochrome P450 2D6: an in silico tool for predicting metabolism and inhibition. J. Med. Chem.47, 53405346. 10.1021/jm049934e

  • 113

    KhalidS.ZahidM. A.AliH.KimY. S.KhanS. (2018). Biaryl scaffold-focused virtual screening for anti-aggregatory and neuroprotective effects in alzheimer's disease. BMC Neurosci.19:74. 10.1186/s12868-018-0472-6

  • 114

    KhanM. S.MajidA. M.IqbalM. A.MajidA. S.Al-MansoubM.HaqueR. S. (2016). Designing the angiogenic inhibitor for brain tumor via disruption of VEGF and IL17A expression. Eur. J. Pharm. Sci.93, 304318. 10.1016/j.ejps.2016.08.032

  • 115

    KhanM. T. (2010). Predictions of the ADMET properties of candidate drug molecules utilizing different QSAR/QSPR modelling approaches. Curr. Drug Metab.11, 285295. 10.2174/138920010791514306

  • 116

    KhuranaR. K.KaurR.KaurM.KaurR.KaurJ.KaurH.et al. (2017). Exploring and validating physicochemical properties of mangiferin through GastroPlus((R)) software. Future Sci. OA3:FSO167. 10.4155/fsoa-2016-0055

  • 117

    KimS.ChenJ.ChengT.GindulyteA.HeJ.HeS.et al. (2019). PubChem 2019 update: improved access to chemical data. Nucleic Acids Res.47, D1102D1109. 10.1093/nar/gky1033

  • 118

    KimS. J.HeoS. H.LeeD. S.HwangI. G.LeeY. B.ChoH. Y. (2016). Gender differences in pharmacokinetics and tissue distribution of 3 perfluoroalkyl and polyfluoroalkyl substances in rats. Food Chem. Toxicol.97, 243255. 10.1016/j.fct.2016.09.017

  • 119

    KirchmairJ.WilliamsonM. J.TyzackJ. D.TanL.BondP. J.BenderA.et al. (2012). Computational prediction of metabolism: sites, products, SAR, P450 enzyme dynamics, and mechanisms. J. Chem. Inf. Model.52, 617648. 10.1021/ci200542m

  • 120

    KolaI.LandisJ. (2004). Can the pharmaceutical industry reduce attrition rates?Nat. Rev. Drug Discov.3, 711715. 10.1038/nrd1470

  • 121

    KollarJ.FrecerV. (2018). Diarylcyclopropane hydroxamic acid inhibitors of histone deacetylase 4 designed by combinatorial approach and QM/MM calculations. J. Mol. Graph Model.85, 97110. 10.1016/j.jmgm.2018.08.008

  • 122

    KongL.TangY.ZhangX.LuG.YuM.ShiQ.et al. (2017). Pharmacokinetic/pharmacodynamic analysis of meropenem for the treatment of nosocomial pneumonia in intracerebral hemorrhage patients by monte carlo simulation. Ann. Pharmacother.51, 970975. 10.1177/1060028017719715

  • 123

    KouD.ZhangC.YiuH.NgT.LubachJ. W.JansonM.et al. (2018). In vitro, in silico, and in vivo assessments of intestinal precipitation and its impact on bioavailability of a BCS class 2 basic compound. Mol. Pharm.15, 16071617. 10.1021/acs.molpharmaceut.7b01143

  • 124

    KumarR.GargP.BharatamP. V. (2015). Shape-based virtual screening, docking, and molecular dynamics simulations to identify Mtb-ASADH inhibitors. J. Biomol. Struct. Dyn.33, 10821093. 10.1080/07391102.2014.929535

  • 125

    LapinsM.ArvidssonS.LampaS.BergA.SchaalW.AlvarssonJ.et al. (2018). A confidence predictor for logD using conformal regression and a support-vector machine. J. Cheminform.10:17. 10.1186/s13321-018-0271-1

  • 126

    LauroG.CantoneV.PotenzaM.FischerK.KoeberleA.WerzO.et al. (2018). Discovery of 3-hydroxy-3-pyrrolin-2-one-based mPGES-1 inhibitors using a multi-step virtual screening protocol. Medchemcomm9, 20282036. 10.1039/C8MD00497H

  • 127

    LaveT.ParrottN.GrimmH. P.FleuryA.ReddyM. (2007). Challenges and opportunities with modelling and simulation in drug discovery and drug development. Xenobiotica37, 12951310. 10.1080/00498250701534885

  • 128

    LawV.KnoxC.DjoumbouY.JewisonT.GuoA. C.LiuY.et al. (2013). DrugBank 4.0: shedding new light on drug metabolism. Nucleic Acids Res.42, D1091D1097. 10.1093/nar/gkt1068

  • 129

    LeaI. A.GongH.PalejaA.RashidA.FostelJ. (2017). CEBS: a comprehensive annotated database of toxicological data. Nucleic Acids Res.45, D964D971. 10.1093/nar/gkw1077

  • 130

    LeeL. S.SengK. Y.WangL. Z.YongW. P.HeeK. H.SohT. I.et al. (2016). Phenotyping of UGT1A1 activity using raltegravir predicts pharmacokinetics and toxicity of irinotecan in FOLFIRI. PLoS ONE11:e0147681. 10.1371/journal.pone.0147681

  • 131

    LiD.WangY.HanK. (2012). Recent density functional theory model calculations of drug metabolism by cytochrome P450. Coord Chem. Rev.256, 11371150. 10.1016/j.ccr.2012.01.016

  • 132

    LiJ.DuH.WuZ.SuH.LiuG.TangY.et al. (2016). Interactions of omeprazole-based analogues with cytochrome P450 2C19: a computational study. Mol. Biosyst.12, 19131921. 10.1039/C6MB00139D

  • 133

    LiM.ZouP.TynerK.LeeS. (2017). Physiologically based pharmacokinetic (PBPK) modeling of pharmaceutical nanoparticles. AAPS J.19, 2642. 10.1208/s12248-016-0010-3

  • 134

    LiX.WangR.LiuY.LiuY.ZhengH.FengY.et al. (2017). Pharmacokinetic study of single- and multiple-dosing with metolazone tablets in healthy Chinese population. BMC Pharmacol. Toxicol.18:73. 10.1186/s40360-017-0178-x

  • 135

    LiY. F.GuoH. Y.YangF.ZhouL. G.HuangX. H.DingH. Z.et al. (2017). Comparative pharmacokinetics of diaveridine in pigs and chickens following single intravenous and oral administration. J. Vet. Pharmacol. Ther.40, 500504. 10.1111/jvp.12384

  • 136

    LinJ.SahakianD. C.de MoraisS. M.XuJ. J.PolzerR. J.WinterS. M. (2003). The role of absorption, distribution, metabolism, excretion and toxicity in drug discovery. Curr. Top. Med. Chem.3, 11251154. 10.2174/1568026033452096

  • 137

    LindstromA.AnantpadmaM.BakerL.RaghavendraN. M.DaveyR.DavissonV. J. (2018). Phenotypic prioritization of diphyllin derivatives that block filoviral cell entry by vacuolar (H+)-ATPase inhibition. Chemmedchem13, 26642676. 10.1002/cmdc.201800587

  • 138

    LuoJ.LaiT.GuoT.ChenF.ZhangL.DingW.et al. (2018). Synthesis and acaricidal activities of scopoletin phenolic ether derivatives: QSAR, molecular docking study and in Silico ADME predictions. Molecules23:995. 10.3390/molecules23050995

  • 139

    MaY.WeiH. Y.ZhangY. Z.JinW. Y.LiH. L.ZhouH.et al. (2017). Synthesis, bioactivity, 3D-QSAR studies of novel dibenzofuran derivatives as PTP-MEG2 inhibitors. Oncotarget8, 3846638481. 10.18632/oncotarget.16595

  • 140

    MacalinoS. J.GosuV.HongS.ChoiS. (2015). Role of computer-aided drug design in modern drug discovery. Arch. Pharm. Res.38, 16861701. 10.1007/s12272-015-0640-5

  • 141

    MalikR.MehtaP.SrivastavaS.ChoudharyB. S.SharmaM. (2016). Pharmacophore modeling, 3D-QSAR, and in silico ADME prediction of N-pyridyl and pyrimidine benzamides as potent antiepileptic agents. J. Recept. Signal. Transduct. Res.37, 259266. 10.1080/10799893.2016.1217883

  • 142

    MaltarolloV. G.GertrudesJ. C.OliveiraP. R.HonorioK. M. (2015). Applying machine learning techniques for ADME-Tox prediction: a review. Expert Opin. Drug Metab. Toxicol.11, 259271. 10.1517/17425255.2015.980814

  • 143

    MaunzA.GutleinM.RautenbergM.VorgrimmlerD.GebeleD.HelmaC. (2013). lazar: a modular predictive toxicology framework. Front. Pharmacol.4:38. 10.3389/fphar.2013.00038

  • 144

    MerlotC. (2010). Computational toxicology—a tool for early safety evaluation. Drug Discov. Today15, 1622. 10.1016/j.drudis.2009.09.010

  • 145

    MertensB.SimonC.Van BossuytM.OnghenaM.VandermarkenT.Van LangenhoveK.et al. (2016). Investigation of the genotoxicity of substances migrating from polycarbonate replacement baby bottles to identify chemicals of high concern. Food Chem. Toxicol.89, 126137. 10.1016/j.fct.2016.01.009

  • 146

    MichielanL.MoroS. (2010). Pharmaceutical perspectives of nonlinear QSAR strategies. J. Chem. Inf. Model.50, 961978. 10.1021/ci100072z

  • 147

    ModiS. (2003). Computational approaches to the understanding of ADMET properties and problems. Drug Discov. Today8, 621623. 10.1016/S1359-6446(03)02739-9

  • 148

    MohanA.KirubakaranR.ParrayJ. A.SivakumarR.MurugeshE.GovarthananM. (2020). Ligand-based pharmacophore filtering, atom based 3D-QSAR, virtual screening and ADME studies for the discovery of potential ck2 inhibitors. J. Mol. Struct.1205:127670. 10.1016/j.molstruc.2019.127670

  • 149

    MojD.BritzH.BurhenneJ.StewartC. F.EgererG.HaefeliW. E.et al. (2017). A physiologically based pharmacokinetic and pharmacodynamic (PBPK/PD) model of the histone deacetylase (HDAC) inhibitor vorinostat for pediatric and adult patients and its application for dose specification. Cancer Chemother. Pharmacol.80, 10131026. 10.1007/s00280-017-3447-x

  • 150

    MolinaI.SalvadorF.Sanchez-MontalvaA.ArtazaM. A.MorenoR.PerinL.et al. (2017). Pharmacokinetics of benznidazole in healthy volunteers and implications in future clinical trials. Antimicrob. Agents Chemother.61, e0191216. 10.1128/AAC.01912-16

  • 151

    MondalS.MandalS. M.MondalT. K.SinhaC. (2017). Spectroscopic characterization, antimicrobial activity, DFT computation and docking studies of sulfonamide schiff bases. J. Mol. Struct.1127, 557567. 10.1016/j.molstruc.2016.08.011

  • 152

    MoroyG.MartinyV. Y.VayerP.VilloutreixB. O.MitevaM. A. (2012). Toward in silico structure-based ADMET prediction in drug discovery. Drug Discov. Today17, 4455. 10.1016/j.drudis.2011.10.023

  • 153

    MusterW.BreidenbachA.FischerH.KirchnerS.MullerL.PahlerA. (2008). Computational toxicology in drug development. Drug Discov. Today13, 303310. 10.1016/j.drudis.2007.12.007

  • 154

    NanJ.HaoH.XieS.PanY.XiC.MaoF.et al. (2017). Pharmacokinetic and pharmacodynamic integration and modeling of acetylkitasamycin in swine for clostridium perfringens. J. Vet. Pharmacol. Ther.40, 641655. 10.1111/jvp.12404

  • 155

    NandekarP. P.KhomaneK.ChaudharyV.RathodV. P.BorkarR. M.BhandiM. M.et al. (2016). Identification of leads for antiproliferative activity on MDA-MB-435 human breast cancer cells through pharmacophore and CYP1A1-mediated metabolism. Eur. J. Med. Chem.115, 8293. 10.1016/j.ejmech.2016.02.061

  • 156

    NasabR. R.MansourianM.HassanzadehF.ShahlaeiM. (2018). Exploring the interaction between epidermal growth factor receptor tyrosine kinase and some of the synthesized inhibitors using combination of in-silico and in-vitro cytotoxicity methods. Res. Pharm. Sci.13, 509522. 10.4103/1735-5362.245963

  • 157

    NestorovI. (2003). Whole body pharmacokinetic models. Clin. Pharmacokinet.42, 883908. 10.2165/00003088-200342100-00002

  • 158

    NiederaltC.KuepferL.SolodenkoJ.EissingT.SiegmundH. U.BlockM.et al. (2017). A generic whole body physiologically based pharmacokinetic model for therapeutic proteins in PK-sim. J. Pharmacokinet. Pharmacodyn.45, 235257. 10.1007/s10928-017-9559-4

  • 159

    NishaC. M.KumarA.VimalA.BaiB. M.PalD.KumarA. (2016). Docking and ADMET prediction of few GSK-3 inhibitors divulges 6-bromoindirubin-3-oxime as a potential inhibitor. J. Mol. Graph Model.65, 100107. 10.1016/j.jmgm.2016.03.001

  • 160

    NiuL.DingL.LuC.ZuoF.YaoK.XuS.et al. (2016). Flavokawain A inhibits Cytochrome P450 in in vitro metabolic and inhibitory investigations. J. Ethnopharmacol.191, 350359. 10.1016/j.jep.2016.06.039

  • 161

    NosratiM.ShakeranZ.ShakeranZ. (2018). Frangulosid as a novel hepatitis B virus DNA polymerase inhibitor: a virtual screening study. In Silico Pharmacol.6:10. 10.1007/s40203-018-0047-3

  • 162

    OlguinH. J.MartinezH. O.PerezC. F.MendiolaB. R.EspinosaL. R.PachecoJ. L. C.et al. (2017). Pharmacokinetics of sildenafil in children with pulmonary arterial hypertension. World J. Pediatr.13, 588592. 10.1007/s12519-017-0043-4

  • 163

    OnghenaM.Van HoeckE.Van LocoJ.IbanezM.ChertaL.PortolesT.et al. (2015). Identification of substances migrating from plastic baby bottles using a combination of low-resolution and high-resolution mass spectrometric analysers coupled to gas and liquid chromatography. J. Mass Spectrom.50, 12341244. 10.1002/jms.3644

  • 164

    OnghenaM.van HoeckE.VervlietP.ScippoM. L.SimonC.van LocoJ.et al. (2014). Development and application of a non-targeted extraction method for the analysis of migrating compounds from plastic baby bottles by GC-MS. Food Addit. Contam. Part A Chem. Anal. Control Expo. Risk Assess. 31, 20902102. 10.1080/19440049.2014.979372

  • 165

    OveringtonJ. (2009). ChEMBL. an interview with john overington, team leader, chemogenomics at the European bioinformatics institute outstation of the European molecular biology laboratory (EMBL-EBI). interview by wendy A. Warr. J. Comput. Aided Mol. Des.23, 195198. 10.1007/s10822-009-9260-9

  • 166

    PadariH.MetsvahtT.GermovsekE.BarkerC. I.KipperK.HerodesK.et al. (2018). Pharmacokinetics of penicillin g in preterm and term neonates. Antimicrob. Agents Chemother.62, 19952000. 10.1128/AAC.02238-17

  • 167

    ParrottN. J.YuL. J.TakanoR.NakamuraM.MorcosP. N. (2016). Physiologically based absorption modeling to explore the impact of food and gastric pH changes on the pharmacokinetics of alectinib. AAPS J.18, 14641474. 10.1208/s12248-016-9957-3

  • 168

    PatelC. N.KumarS. P.RawalR. M.PatelD. P.GonzalezF. J.PandyaH. A. (2020). A multiparametric organ toxicity predictor for drug discovery. Toxicol. Mech. Methods30, 159166. 10.1080/15376516.2019.1681044

  • 169

    PatelD. B.DarjiD. G.PatelK. R.RajaniD. P.RajaniS. D.PatelH. D. (2020). Synthesis of novel quinoline-thiosemicarbazide hybrids and evaluation of their biological activities, molecular docking, molecular dynamics, pharmacophore model studies, and ADME-tox properties. J. Heterocyclic Chem. 57, 11831200. 10.1002/jhet.3855

  • 170

    PatelN.JainS.LinS. (2016). Transdermal iontophoretic delivery of tacrine hydrochloride: correlation between in vitro permeation and in vivo performance in rats. Int. J. Pharm.513, 393403. 10.1016/j.ijpharm.2016.09.038

  • 171

    PathakS. M.SchaeferK. J.JameiM.TurnerD. B. (2019). Biopharmaceutic IVIVE-mechanistic modeling of single- and two-phase in vitro experiments to obtain drug-specific parameters for incorporation into PBPK models. J. Pharm. Sci.108, 16041618. 10.1016/j.xphs.2018.11.034

  • 172

    PatlewiczG.JeliazkovaN.SaffordR. J.WorthA. P.AleksievB. (2008). An evaluation of the implementation of the cramer classification scheme in the toxtree software. SAR QSAR Environ. Res.19, 495524. 10.1080/10629360802083871

  • 173

    PattersonD. E.CramerR. D.FergusonA. M.ClarkR. D.WeinbergerL. E. (1996). Neighborhood behavior: a useful concept for validation of “molecular diversity” descriptors. J. Med. Chem. 39, 30493059. 10.1021/jm960290n

  • 174

    PerdaemsN.BlascoH.VinsonC.ChenelM.WhalleyS.CazadeF.et al. (2010). Predictions of metabolic drug-drug interactions using physiologically based modelling: two cytochrome P450 3A4 substrates coadministered with ketoconazole or verapamil. Clin. Pharmacokinet.49, 239258. 10.2165/11318130-000000000-00000

  • 175

    Perez-NuenoV. I.RitchieD. W. (2011). Using consensus-shape clustering to identify promiscuous ligands and protein targets and to choose the right query for shape-based virtual screening. J. Chem. Inf. Model.51, 12331248. 10.1021/ci100492r

  • 176

    PetrescuA. M.PaunescuV.IliaG. (2019). The antiviral activity and cytotoxicity of 15 natural phenolic compounds with previously demonstrated antifungal activity. J. Environ. Sci. Health B54, 498504. 10.1080/03601234.2019.1574176

  • 177

    PilariS.GaubT.BlockM.GorlitzL. (2017). Development of physiologically based organ models to evaluate the pharmacokinetics of drugs in the testes and the thyroid gland. CPT Pharmacometr. Syst. Pharmacol.6, 532542. 10.1002/psp4.12205

  • 178

    PiresD. E.BlundellT. L.AscherD. B. (2015). pkCSM: predicting small-molecule pharmacokinetic and toxicity properties using graph-based signatures. J. Med. Chem.58, 40664072. 10.1021/acs.jmedchem.5b00104

  • 179

    PoirierD.RoyJ.Cortes-BenitezF.DutourR. (2016). Targeting cytochrome P450 (CYP) 1B1 with steroid derivatives. Bioorg. Med. Chem. Lett.26, 52725276. 10.1016/j.bmcl.2016.09.046

  • 180

    PolepallyA. R.BrundageR. C.RemmelR. P.LeppikI. E.PennellP. B.WhiteJ. R.et al. (2018). Lamotrigine pharmacokinetics following oral and stable-labeled intravenous administration in young and elderly adult epilepsy patients: effect of age. Epilepsia59, 17181726. 10.1111/epi.14519

  • 181

    PrabhuS. V.SinghS. K. (2019). Identification of potential dual negative allosteric modulators of group I mGluR family: a shape based screening, ADME prediction, induced fit docking and molecular dynamics approach against neurodegenerative diseases. Curr. Top. Med. Chem.19, 26872707. 10.2174/1568026619666191105112800

  • 182

    PradeepP.PovinelliR. J.WhiteS.MerrillS. J. (2016). An ensemble model of QSAR tools for regulatory risk assessment. J. Cheminform8:48. 10.1186/s13321-016-0164-0

  • 183

    ProkopJ. W.LazarJ.CrapittoG.SmithD. C.WortheyE. A.JacobH. J. (2017). Molecular modeling in the age of clinical genomics, the enterprise of the next generation. J. Mol. Model.23:75. 10.1007/s00894-017-3258-3

  • 184

    PuttaS.BerozaP. (2007). Shapes of things: computer modeling of molecular shape in drug discovery. Curr. Top. Med. Chem.7, 15141524. 10.2174/156802607782194770

  • 185

    RacharV.CzejkaM.KitzmuellerM.BuchnerP.LichtneckertM.GreilR.et al. (2016). Assessment of pharmacokinetic interaction between capecitabine and cetuximab in metastatic colorectal cancer patients. Anticancer Res.36, 47154723. 10.21873/anticanres.11026

  • 186

    RamatenkiV.DumpatiR.VadijaR.VellankiS.PotlapallyS. R.RondlaR.et al. (2016). Targeting the ubiquitin-conjugating enzyme E2D4 for cancer drug discovery-a structure-based approach. J. Chem. Biol.10, 5167. 10.1007/s12154-016-0164-6

  • 187

    RamatenkiV.DumpatiR.VadijaR.VellankiS.PotlapallyS. R.RondlaR.et al. (2017). Identification of new lead molecules against UBE2NL enzyme for cancer therapy. Appl. Biochem. Biotechnol.182, 14971517. 10.1007/s12010-017-2414-7

  • 188

    RambiritchV.NaidooP.MaharajB.PillaiG. (2016). Population pharmacokinetic modeling of glibenclamide in poorly controlled South African type 2 diabetic subjects. Clin. Pharmacol.8, 8392. 10.2147/CPAA.S102676

  • 189

    RawatR.VermaS. M. (2020). High-throughput virtual screening approach involving pharmacophore mapping, ADME filtering, molecular docking and MM-GBSA to identify new dual target inhibitors of PfDHODH and PfCytbc1 complex to combat drug resistant malaria. J. Biomol. Struct. Dyn.10.1080/07391102.2020.1784288. [Epub ahead of print].

  • 190

    ReddyK. K.SinghS. K.TripathiS. K.SelvarajC.SuryanarayananV. (2013). Shape and pharmacophore-based virtual screening to identify potential cytochrome P450 sterol 14 alpha-demethylase inhibitors. J. Recept. Signal. Transd.33, 234243. 10.3109/10799893.2013.789912

  • 191

    RichardA. M.WilliamsC. R. (2002). Distributed structure-searchable toxicity (DSSTox) public database network: a proposal. Mutat. Res.499, 2752. 10.1016/S0027-5107(01)00289-5

  • 192

    RomanM.RomanD. L.OstafeV.CiorsacA.IsvoranA. (2018). Computational assessment of pharmacokinetics and biological effects of some anabolic and androgen steroids. Pharm. Res.35:41. 10.1007/s11095-018-2353-1

  • 193

    RomanskiM.KasprzykA.KarbownikA.GlowkaF. K. (2016). Ocular disposition of treosulfan and its active epoxy-transformers following intravenous administration in rabbits. Drug Metab. Pharmacokinet.31, 356362. 10.1016/j.dmpk.2016.07.001

  • 194

    RomanskiM.KasprzykA.WalczakM.ZiolkowskaA.GlowkaF. (2017). Disposition of treosulfan and its active monoepoxide in a bone marrow, liver, lungs, brain, and muscle: studies in a rat model with clinical relevance. Eur. J. Pharm. Sci.109, 616623. 10.1016/j.ejps.2017.09.011

  • 195

    Rostami-HodjeganA. (2012). Physiologically based pharmacokinetics joined with in vitro-in vivo extrapolation of ADME: a marriage under the arch of systems pharmacology. Clin. Pharmacol. Ther.92, 5061. 10.1038/clpt.2012.65

  • 196

    RowlandM.BalantL.PeckC. (2004). Physiologically based pharmacokinetics in drug development and regulatory science: a workshop report (Georgetown University, Washington, DC, May 29-30, 2002). AAPS PharmSci.6:E6. 10.1208/ps060106

  • 197

    RowlandM.PeckC.TuckerG. (2011). Physiologically-based pharmacokinetics in drug development and regulatory science. Annu. Rev. Pharmacol. Toxicol.51, 4573. 10.1146/annurev-pharmtox-010510-100540

  • 198

    RoyK.MitraI.KarS.OjhaP. K.DasR. N.KabirH. (2012). Comparative studies on some metrics for external validation of QSPR models. J. Chem. Inf. Model.52, 396408. 10.1021/ci200520g

  • 199

    SabaN.SealA. (2018). Identification of a less toxic vinca alkaloid derivative for use as a chemotherapeutic agent, based on in silico structural insights and metabolic interactions with CYP3A4 and CYP3A5. J. Mol. Model.24:82. 10.1007/s00894-018-3611-1

  • 200

    SangleG. V.PatilM.DeshmukhN. J.ShenguleS. A.KambleS.VuppalavanchuK. K.et al. (2018). Evaluation of pharmacokinetic and pharmacodynamic parameters following single dose of sitagliptin in healthy Indian males. Eur. J. Clin. Pharmacol.74, 561569. 10.1007/s00228-018-2433-5

  • 201

    SanniD. M.FatokiT. H. (2017). Computational evaluation of pharmacokinetics and potential protein targets of ginger (Zingiber officinale). J. Microbiol. Biotechnol. Res. 7:14. 10.24896/jmbr.2017713

  • 202

    SanniD. M.FatokiT. H.KolawoleA. O.AkinmoladunA. C. (2017). Xeronine structure and function: computational comparative mastery of its mystery. In Silico Pharmacol.5:8. 10.1007/s40203-017-0028-y

  • 203

    Santes-PalaciosR.Romo-MancillasA.Camacho-CarranzaR.Espinosa-AguirreJ. J. (2016). Inhibition of human and rat CYP1A1 enzyme by grapefruit juice compounds. Toxicol. Lett.258, 267275. 10.1016/j.toxlet.2016.07.023

  • 204

    SasaharaK.MashimaA.YoshidaT.ChumanH. (2015). Molecular dynamics and density functional studies on the metabolic selectivity of antipsychotic thioridazine by cytochrome P450 2D6: connection with crystallographic and metabolic results. Bioorg. Med. Chem.23, 54595465. 10.1016/j.bmc.2015.07.043

  • 205

    SchlenderJ. F.MeyerM.ThelenK.KraussM.WillmannS.EissingT.et al. (2016). Development of a whole-body physiologically based pharmacokinetic approach to assess the pharmacokinetics of drugs in elderly individuals. Clin. Pharmacokinet.55, 15731589. 10.1007/s40262-016-0422-3

  • 206

    SchmidtU.StruckS.GrueningB.HossbachJ.JaegerI. S.ParolR.et al. (2009). SuperToxic: a comprehensive database of toxic compounds. Nucleic Acids Res.37, D295D299. 10.1093/nar/gkn850

  • 207

    SchmittW.WillmannS. (2004). Physiology-based pharmacokinetic modeling: ready to be used. Drug Discov. Today Technol.1, 449456. 10.1016/j.ddtec.2004.09.006

  • 208

    Schrödinger Suite (2012). Schrödinger, LLC, New York, NY.

  • 209

    SchymanP.LiuR.DesaiV.WallqvistA. (2017). vNN web server for ADMET predictions. Front. Pharmacol.8:889. 10.3389/fphar.2017.00889

  • 210

    SelickH. E.BeresfordA. P.TarbitM. H. (2002). The emerging importance of predictive ADME simulation in drug discovery. Drug Discov. Today7, 109116. 10.1016/S1359-6446(01)02100-6

  • 211

    ShaabanS.Vervandier-FasseurD.AndreolettiP.ZarroukA.RichardP.NegmA.et al. (2018). Cytoprotective and antioxidant properties of organic selenides for the myelin-forming cells, oligodendrocytes. Bioorg. Chem.80, 4356. 10.1016/j.bioorg.2018.05.019

  • 212

    ShaikM.ShaikS.KilariE. K. (2018). Population pharmacokinetics of gliclazide in normal and diabetic rabbits. Biopharm. Drug Dispos.39, 265274. 10.1002/bdd.2132

  • 213

    ShajiD. (2018). Molecular docking studies of human MCT8 protein with soy isoflavones in allan-herndon-dudley syndrome (AHDS). J. Pharm. Anal.8, 318323. 10.1016/j.jpha.2018.07.001

  • 214

    ShakhnovichV.SmithP. B.GuptillJ. T.JamesL. P.CollierD. N.WuH.et al. (2018). Obese children require lower doses of pantoprazole than nonobese peers to achieve equal systemic drug exposures. J. Pediatr.193, 102108.e101. 10.1016/j.jpeds.2017.10.011

  • 215

    ShanQ.FanJ.WangJ.ZhuX.YinY.ZhengG. (2017). Pharmacokinetics of enrofloxacin after oral, intramuscular and bath administration in crucian carp (Carassius auratus gibelio). J. Vet. Pharmacol. Ther.41, 159162. 10.1111/jvp.12428

  • 216

    ShanQ.WangJ.YangF.MaL.YinY.LiuS.et al. (2018). Pharmacokinetics of cefquinome in crucian carp (Carassius auratus gibelio) after oral, intramuscular, intraperitoneal, and bath administration. J. Vet. Pharmacol. Ther.41, 734738. 10.1111/jvp.12510

  • 217

    ShangJ.SunH.LiuH.ChenF.TianS.PanP.et al. (2017). Comparative analyses of structural features and scaffold diversity for purchasable compound libraries. J. Cheminform.9:25. 10.1186/s13321-017-0212-4

  • 218

    SheridanR. P.FeustonB. P.MaiorovV. N.KearsleyS. K. (2004). Similarity to molecules in the training set is a good discriminator for prediction accuracy in QSAR. J. Chem. Inf. Comput. Sci.44, 19121928. 10.1021/ci049782w

  • 219

    ShuklaR.SinghA. P.SonarP. K.MishraM.SarafS. K. (2016). Schiff bases of benzothiazol-2-ylamine and thiazolo[5,4-b] pyridin-2-ylamine as anticonvulsants: synthesis, characterization and toxicity profiling. Cent. Nerv. Syst. Agents Med. Chem.16, 240248. 10.2174/1871524916666160428110728

  • 220

    SilvaG.ZuravskiL.DuarteM.MachadoM.OliveiraL. (2019). Fluconazole induces genotoxicity in cultured human peripheral blood mononuclear cells via immunomodulation of TNF-α, IL-6, and IL-10: new challenges for safe therapeutic regimens. Immunopharmacol. Immunotoxicol.41:123129. 10.1080/08923973.2019.1566357

  • 221

    Silva-JuniorE. F.AquinoT. M.Araujo-JuniorJ. X. (2017). Quantum mechanical (QM) calculations applied to ADMET drug prediction: a review. Curr. Drug Metab.18, 511526. 10.2174/1389200218666170316094514

  • 222

    Simhadri VsdnaN.MuniappanM.KannanI.ViswanathanS. (2017). Phytochemical analysis and docking study of compounds present in a polyherbal preparation used in the treatment of dermatophytosis. Curr. Med. Mycol.3, 614. 10.29252/cmm.3.4.6

  • 223

    SimoneauC.Van den EedeL.ValzacchiS. (2012). Identification and quantification of the migration of chemicals from plastic baby bottles used as substitutes for polycarbonate. Food Addit. Contam. Part A Chem. Anal. Control Expo. Risk Assess.29, 469480. 10.1080/19440049.2011.644588

  • 224

    SongJ. C.GaoH.QiuH. B.ChenQ. B.CaiM. H.ZhangM. Z.et al. (2018). The pharmacokinetics of dexmedetomidine in patients with obstructive jaundice: a clinical trial. PLoS ONE13:e0207427. 10.1371/journal.pone.0207427

  • 225

    SongL.YuZ.XuY.LiX.LiuX.LiuD.et al. (2020). Preliminary physiologically based pharmacokinetic modeling of renally cleared drugs in chinese pregnant women. Biopharm. Drug Dispos. 41:248267. 10.1002/bdd.2243

  • 226

    StillhartC.ParrottN. J.LindenbergM.ChalusP.BentleyD.SzepesA. (2017). Characterising drug release from immediate-release formulations of a poorly soluble compound, basmisanil, through absorption modelling and dissolution testing. AAPS J.19, 827836. 10.1208/s12248-017-0060-1

  • 227

    SunH.ScottD. O. (2010). Structure-based drug metabolism predictions for drug design. Chem. Biol. Drug Des.75, 317. 10.1111/j.1747-0285.2009.00899.x

  • 228

    SushkoI.NovotarskyiS.KornerR.PandeyA. K.CherkasovA.LiJ.et al. (2010). Applicability domains for classification problems: benchmarking of distance to models for ames mutagenicity set. J. Chem. Inf. Model.50, 20942111. 10.1021/ci100253r

  • 229

    TabeshpourJ.SahebkarA.ZirakM. R.ZeinaliM.HashemzaeiM.RakhshaniS.et al. (2018). Computer-aided drug design and drug pharmacokinetic prediction: a mini-review. Curr. Pharm. Des.24, 30143019. 10.2174/1381612824666180903123423

  • 230

    TakacM. J.-M.Casimiro MaginaJ. D.TakacT. (2019). Evaluation of phenylethylamine type entactogens and their metabolites relevant to ecotoxicology - a QSAR study. Acta Pharm.69, 563584. 10.2478/acph-2019-0038

  • 231

    TanB. H.AhemadN.PanY.PalanisamyU. D.OthmanI.YiapB. C.et al. (2018). Cytochrome P450 2C9-natural antiarthritic interactions: evaluation of inhibition magnitude and prediction from in vitro data. Biopharm. Drug Dispos.39, 205217. 10.1002/bdd.2127

  • 232

    TangY.ZhuW.ChenK.JiangH. (2006). New technologies in computer-aided drug design: toward target identification and new chemical entity discovery. Drug Discov. Today Technol.3, 307313. 10.1016/j.ddtec.2006.09.004

  • 233

    TaoL.ZhangP.QinC.ChenS. Y.ZhangC.ChenZ.et al. (2015). Recent progresses in the exploration of machine learning methods as in-silico ADME prediction tools. Adv Drug Deliv. Rev.86, 83100. 10.1016/j.addr.2015.03.014

  • 234

    TaroncherM.TolosaJ.ProsperiniA.RuizM. J. (2018). In silico and in vitro prediction of the toxicological effects of individual and combined mycotoxins. Food Chem. Toxicol.122, 194202. 10.1016/j.fct.2018.09.055

  • 235

    TaxakN.PrasadK. C.BharatamP. V. (2013). Mechanistic insights into the bioactivation of phenacetin to reactive metabolites: a DFT study. Comput. Theor. Chem.1007, 4856. 10.1016/j.comptc.2012.11.018

  • 236

    TeorellT. (1935). Studies on the “Diffusion Effect” upon ionic distribution. Some theoretical considerations. Proc. Natl. Acad. Sci. U.S.A.21, 152161. 10.1073/pnas.21.3.152

  • 237

    ThielC.CordesH.CondeI.CastellJ. V.BlankL. M.KuepferL. (2016). Model-based contextualization of in vitro toxicity data quantitatively predicts in vivo drug response in patients. Arch. Toxicol.91, 865883. 10.1007/s00204-016-1723-x

  • 238

    TianS.LiY.WangJ.ZhangJ.HouT. (2011). ADME evaluation in drug discovery. 9. Prediction of oral bioavailability in humans based on molecular properties and structural fingerprints. Mol. Pharm.8, 841851. 10.1021/mp100444g

  • 239

    T'JollynH.VermeulenA.Van BocxlaerJ. (2018). PBPK and its virtual populations: the impact of physiology on pediatric pharmacokinetic predictions of tramadol. AAPS J.21:8. 10.1208/s12248-018-0277-7

  • 240

    TripathiA. C.UpadhyayS.PaliwalS.SarafS. K. (2018). N1-benzenesulfonyl-2-pyrazoline hybrids in neurological disorders: syntheses, biological screening and computational studies. EXCLI J.17, 126148. 10.17179/excli2017-871

  • 241

    TrivediA.DaveS.LuhanaK.TrivediS. (2020). Design and prediction of absorption, distribution, metabolism, excretion and toxicity properties of drugs for H1N1 Flu (Swine Flu) using in-silico approaches. Biosci. Biotech. Res. Commun.13, 326334. 10.21786/bbrc/13.1/52

  • 242

    TrivediA.DaveS. M.PatelJ. S. (2019). Design and prediction of absorption, distribution, metabolism, excretion and toxicity properties of drugs for dengue fever using in-silico approaches. Biosci. Biotech. Res. Commun.12, 12011209. 10.21786/bbrc/12.4/45

  • 243

    UbaA. I.YelekciK. (2017). Identification of potential isoform-selective histone deacetylase inhibitors for cancer therapy: a combined approach of structure-based virtual screening, ADMET prediction and molecular dynamics simulation assay. J. Biomol. Struct. Dyn.36, 32313245. 10.1080/07391102.2017.1384402

  • 244

    VaidhyanathanS.WangX.CrisonJ.VariaS.GaoJ. Z. H.SaxenaA.et al. (2018). Bioequivalence comparison of pediatric dasatinib formulations and elucidation of absorption mechanisms through integrated PBPK modeling. J. Pharm. Sci.108, 741749. 10.1016/j.xphs.2018.11.005

  • 245

    Van BossuytM.Van HoeckE.RaitanoG.VanhaeckeT.BenfenatiE.MertensB.et al. (2018). Performance of in silico models for mutagenicity prediction of food contact materials. Toxicol. Sci.163, 632638. 10.1093/toxsci/kfy057

  • 246

    van de WaterbeemdH.GiffordE. (2003). ADMET in silico modelling: towards prediction paradise?Nat. Rev. Drug Discov.2, 192204. 10.1038/nrd1032

  • 247

    VazR. J.LiY.MunsonM.ElliotM.ThurairatnamS. (2018). Amelioration of mechanism-based inactivation of CYP3A4 by a H-PGDS inhibitor. Bioorg. Med. Chem. Lett.28, 30463049. 10.1016/j.bmcl.2018.07.049

  • 248

    VerheyenG. R.BraekenE.Van DeunK.Van MiertS. (2017). Evaluation of existing (Q)SAR models for skin and eye irritation and corrosion to use for REACH registration. Toxicol. Lett.265, 4752. 10.1016/j.toxlet.2016.11.007

  • 249

    VermaS. K.TharejaS. (2017). Structure based comprehensive modelling, spatial fingerprints mapping and ADME screening of curcumin analogues as novel ALR2 inhibitors. PLoS ONE12:e0175318. 10.1371/journal.pone.0175318

  • 250

    VoelknerN. M. F.VoelknerA.CostaJ.SyS. K. B.HermesJ.WeitzelJ.et al. (2017). Dermal pharmacokinetics of pyrazinamide determined by microdialysis sampling in rats. Int. J. Antimicrob. Agents51, 190196. 10.1016/j.ijantimicag.2017.10.001

  • 251

    VuuI.ColesL. D.MaglalangP.LeppikI. E.WorrellG.CrepeauD.et al. (2016). Intravenous topiramate: pharmacokinetics in dogs with naturally occurring epilepsy. Front. Vet. Sci.3:107. 10.3389/fvets.2016.00107

  • 252

    WangJ.HouT. (2009). Recent advances on in silico ADME modeling. Annu. Rep. Comput. Chem.5, 101127. 10.1016/S1574-1400(09)00505-2

  • 253

    WangJ.PengH.KongJ.ZhaoT.ZhangS.CaoX. (2017). Pharmacokinetic profile of ceftiofur hydrochloride injection in lactating holstein dairy cows. J. Vet. Pharmacol. Ther.41, 301306. 10.1111/jvp.12469

  • 254

    WangJ.XueJ.KongJ.LiJ.ZhangS.CaoX. (2018a). Evaluation of pharmacokinetic properties of vitacoxib in fasted and fed horses. J. Vet. Pharmacol. Ther.41, 843847. 10.1111/jvp.12703

  • 255

    WangJ.ZhaoT.SunX.LiuY.ZhuJ.ZhangS.et al. (2018b). Pharmacokinetics of tildipirosin in beagle dogs. J. Vet. Pharmacol. Ther.41, E49E52. 10.1111/jvp.12453

  • 256

    WangJ. H.EndsleyA. N.GreenC. E.MatinA. C. (2016). Utilizing native fluorescence imaging, modeling and simulation to examine pharmacokinetics and therapeutic regimen of a novel anticancer prodrug. BMC Cancer16:524. 10.1186/s12885-016-2508-6

  • 257

    WangJ. Y.ChenH.WangY. Y.WangX. Q.ChenH. Y.ZhangM.et al. (2017). Network pharmacological mechanisms of Vernonia anthelmintica (L.) in the treatment of vitiligo: isorhamnetin induction of melanogenesis via up-regulation of melanin-biosynthetic genes. BMC Syst. Biol.11:103. 10.1186/s12918-017-0486-1

  • 258

    WangL.HaiY.AnL.ChenJ.LiangR.HeX. (2017). Rapid screening the potential mechanism-based inhibitors of CYP3A4 from Tripterygium wilfordi based on computer approaches combined with in vitro bioassay. Bioorg. Med. Chem.25, 26892700. 10.1016/j.bmc.2017.03.037

  • 259

    WangL.YueH.HuangN.XueG.LiuW.XinH. (2018). Human cytochrome P450 enzyme inhibition profile of three flavonoids isolated from Psoralea corylifolia: in silico predictions and experimental validation. N J. Chem.42, 1092210934. 10.1039/C7NJ00884H

  • 260

    WangS.SunH.LiuH.LiD.LiY.HouT. (2016). ADMET evaluation in drug discovery. 16. predicting hERG blockers by combining multiple pharmacophores and machine learning approaches. Mol. Pharm.13, 285566. 10.1021/acs.molpharmaceut.6b00471

  • 261

    WangY.LinW.WuN.HeX.WangJ.FengZ.et al. (2018). An insight into paracetamol and its metabolites using molecular docking and molecular dynamics simulation. J. Mol. Model.24:243. 10.1007/s00894-018-3790-9

  • 262

    WangY.XingJ.XuY.ZhouN.PengJ.XiongZ.et al. (2015). In silico ADME/T modelling for rational drug design. Q. Rev. Biophys.48, 488515. 10.1017/S0033583515000190

  • 263

    WaseemD.ButtA. F.HaqI. U.BhattiM. H.KhanG. M. (2017). Daru25:8. 10.1186/s40199-017-0174-0

  • 264

    WatanabeY.FukuyoshiS.KatoK.HiratsukaM.YamaotsuN.HironoS.et al. (2017). Investigation of substrate recognition for cytochrome P450 1A2 mediated by water molecules using docking and molecular dynamics simulations. J. Mol. Graph Model.74, 326336. 10.1016/j.jmgm.2017.04.006

  • 265

    WatfordS.Ly LyP.WignallJ.ShinR.MartinM. T.FriedmanK. P. (2019). ToxRefDB version 2.0: improved utility for predictive and retrospective toxicology analyses. Reprod. Toxicol.89, 145158. 10.1016/j.reprotox.2019.07.012

  • 266

    WhiteJ.WrzesinskiC.GreenM.JohnsonG. T.McCluskeyJ. D.AbritisA.et al. (2016). A novel method for deriving thresholds of toxicological concern for vaccine constituents. Toxicol. Mech. Methods26, 270275. 10.3109/15376516.2016.1170250

  • 267

    WilliamsA. J.GrulkeC. M.EdwardsJ.McEachranA. D.MansouriK.BakerN. C.et al. (2017). The comptox chemistry dashboard: a community data resource for environmental chemistry. J. Cheminform.9:61. 10.1186/s13321-017-0247-6

  • 268

    WishartD. S. (2007). Improving early drug discovery through ADME modelling: an overview. Drugs R D8, 349362. 10.2165/00126839-200708060-00003

  • 269

    WishartD. S.FeunangY. D.GuoA. C.LoE. J.MarcuA.GrantJ. R.et al. (2018). DrugBank 5.0: a major update to the drugbank database for 2018. Nucleic Acids Res.46, D1074D1082. 10.1093/nar/gkx1037

  • 270

    WishartD. S.KnoxC.GuoA. C.ChengD.ShrivastavaS.TzurD.et al. (2008). DrugBank: a knowledgebase for drugs, drug actions and drug targets. Nucleic Acids Res.36, D901D906. 10.1093/nar/gkm958

  • 271

    WishartD. S.KnoxC.GuoA. C.ShrivastavaS.HassanaliM.StothardP.et al. (2006). DrugBank: a comprehensive resource for in silico drug discovery and exploration. Nucleic Acids Res.34, D668D672. 10.1093/nar/gkj067

  • 272

    WuZ.LeiT.ShenC.WangZ.CaoD.HouT. (2019). ADMET evaluation in drug discovery. 19. reliable prediction of human cytochrome P450 inhibition using artificial intelligence approaches. J. Chem. Inf. Model.59, 45874601. 10.1021/acs.jcim.9b00801

  • 273

    XiaoH.PengH.ZhaoT.KongJ.XueJ.WangJ.et al. (2018). The pharmacokinetics of moxidectin following intravenous and topical administration to swine. J. Vet. Pharmacol. Ther.42, 111115. 10.1111/jvp.12693

  • 274

    XuY.PeiJ.LaiL. (2017). Deep learning based regression and multiclass models for acute oral toxicity prediction with automatic chemical feature extraction. J. Chem. Inf. Model.57, 26722685. 10.1021/acs.jcim.7b00244

  • 275

    YadavS.PandeyS. K.SinghV. K.GoelY.KumarA.SinghS. M. (2017). Molecular docking studies of 3-bromopyruvate and its derivatives to metabolic regulatory enzymes: implication in designing of novel anticancer therapeutic strategies. PLoS ONE12:e0176403. 10.1371/journal.pone.0176403

  • 276

    YamashitaF.HashidaM. (2004). In silico approaches for predicting ADME properties of drugs. Drug Metab. Pharmacokinet.19, 327338. 10.2133/dmpk.19.327

  • 277

    YangF.WangB.LiuZ.XiaX.WangW.YinD.et al. (2017). Prediction of a therapeutic dose for buagafuran, a potent anxiolytic agent by physiologically based pharmacokinetic/pharmacodynamic modeling starting from pharmacokinetics in rats and human. Front. Pharmacol.8:683. 10.3389/fphar.2017.00683

  • 278

    YangH.LouC.SunL.LiJ.CaiY.WangZ.et al. (2018a). admetSAR 2.0: web-service for prediction and optimization of chemical ADMET properties. Bioinformatics35, 10671069. 10.1093/bioinformatics/bty707

  • 279

    YangH.SunL.LiW.LiuG.TangY. (2018b). In silico prediction of chemical toxicity for drug design using machine learning methods and structural alerts. Front. Chem.6:30. 10.3389/fchem.2018.00129

  • 280

    YangH.SunL.WangZ.LiW.LiuG.TangY. (2018c). ADMETopt: a web server for ADMET optimization in drug design via scaffold hopping. J. Chem. Inf. Model.58, 20512056. 10.1021/acs.jcim.8b00532

  • 281

    YangZ. Y.YangZ. J.DongJ.WangL. L.ZhangL. X.DingJ. J.et al. (2019). Structural analysis and identification of colloidal aggregators in drug discovery. J. Chem. Inf. Model.59, 37143726. 10.1021/acs.jcim.9b00541

  • 282

    YeN.MonkS. A.DagaP.BenderD. M.RosenL. B.MullenJ.et al. (2018). Clinical bioavailability of the novel BACE1 inhibitor lanabecestat (AZD3293): assessment of tablet formulations versus an oral solution and the impact of gastric pH on pharmacokinetics. Clin. Pharmacol. Drug Dev.7, 233243. 10.1002/cpdd.422

  • 283

    YongyeA. B.Medina-FrancoJ. L. (2013). Systematic characterization of structure-activity relationships and ADMET compliance: a case study. Drug Discov. Today18, 732739. 10.1016/j.drudis.2013.04.002

  • 284

    YuH.AdedoyinA. (2003). ADME-Tox in drug discovery: integration of experimental and computational technologies. Drug Discov. Today8, 852861. 10.1016/S1359-6446(03)02828-9

  • 285

    YuY.ZhouY. F.SunJ.ShiW.LiaoX. P.LiuY. H. (2017). Pharmacokinetic and pharmacodynamic modeling of sarafloxacin against avian pathogenic Escherichia coli in muscovy ducks. BMC Vet. Res.13:47. 10.1186/s12917-017-0964-0

  • 286

    YugandharP.KumarK. K.NeerajaP.SavithrammaN. (2017). Isolation, characterization and in silico docking studies of synergistic estrogen receptor a anticancer polyphenols from Syzygium alternifolium (Wt.) Walp. J. Intercult. Ethnopharmacol.6, 296310. 10.5455/jice.20170709031835

  • 287

    ZhangC.DingS.FangY.ZhangL.HuW.LuJ.et al. (2017). Iohexol clearance for determination of glomerular filtration rate in cynomolgus monkeys (Macaca fascicularis). J Am Assoc Lab Anim Sci.56, 330333.

  • 288

    ZhangL.ZhangH.AiH. X.HuH.LiS. M.ZhaoJ.et al. (2018). Applications of machine learning methods in drug toxicity prediction. Curr. Top. Med. Chem.18, 987997. 10.2174/1568026618666180727152557

  • 289

    ZhangL.ZhaoL.LiuY.LiuJ.LiX. (2017). Pharmacokinetics of tilmicosin in healthy pigs and in pigs experimentally infected with haemophilus parasuis. J. Vet. Sci.18, 431437. 10.4142/jvs.2017.18.4.431

  • 290

    ZhangN.WuY.HuangZ.YaoL.ZhangL.CaiQ.et al. (2017). The PK-PD relationship and resistance development of danofloxacin against mycoplasma gallisepticum in an in vivo infection model. Front. Microbiol.8:926. 10.3389/fmicb.2017.00926

  • 291

    ZhangP.HaoH.LiJ.AhmadI.ChengG.ChenD.et al. (2016). The epidemiologic and pharmacodynamic cutoff values of tilmicosin against haemophilus parasuis. Front. Microbiol.7:385. 10.3389/fmicb.2016.00385

  • 292

    ZhaoP.ZhangL.GrilloJ. A.LiuQ.BullockJ. M.MoonY. J.et al. (2011). Applications of physiologically based pharmacokinetic (PBPK) modeling and simulation during regulatory review. Clin. Pharmacol. Ther.89, 259267. 10.1038/clpt.2010.298

  • 293

    ZhouD.AfzeliusL.GrimmS. W.AnderssonT. B.ZauharR. J.ZamoraI. (2006). Comparison of methods for the prediction of the metabolic sites for CYP3A4-mediated metabolic reactions. Drug Metab. Dispos.34, 976983. 10.1124/dmd.105.008631

  • 294

    ZhouQ.ZhangG.WangQ.LiuW.HuangY.YuP.et al. (2017). Pharmacokinetic/pharmacodynamic modeling of tulathromycin against pasteurella multocida in a porcine tissue cage model. Front. Pharmacol.8:392. 10.3389/fphar.2017.00392

  • 295

    ZhuJ.WangJ.YuH.LiY.HouT. (2011). Recent developments of in silico predictions of oral bioavailability. Comb. Chem. High Throughput Screen.14, 362374. 10.2174/138620711795508368

  • 296

    ZhuL.ZhangY.YangJ.WangY.ZhangJ.ZhaoY.et al. (2016). Prediction of the pharmacokinetics and tissue distribution of levofloxacin in humans based on an extrapolated PBPK model. Eur. J. Drug Metab. Pharmacokinet.41, 395402. 10.1007/s13318-015-0271-8

  • 297

    ZhuangJ.XingX.WangD.DuZ.WangJ.DongY.et al. (2018). Toxicity assessment of the extractables from multi-layer coextrusion poly ethylene bags exposed to pH=5 solution containing 4% benzyl alcohol and 0.1M sodium acetate. Regul. Toxicol. Pharmacol.94, 4756. 10.1016/j.yrtph.2018.01.015

  • 298

    ZhuangS.ZhangL.ZhanT.LuL.ZhaoL.WangH.et al. (2018). Binding specificity determines the cytochrome P450 3A4 mediated enantioselective metabolism of metconazole. J. Phys. Chem. B122, 11761184. 10.1021/acs.jpcb.7b11170

  • 299

    ZhuangX.LuC. (2016). PBPK modeling and simulation in drug research and development. Acta Pharm. Sin. B6, 430440. 10.1016/j.apsb.2016.04.004

Summary

Keywords

drug discovery, pre-clinical studies, ADMET, pharmacokinetics, PBPK modeling

Citation

Wu F, Zhou Y, Li L, Shen X, Chen G, Wang X, Liang X, Tan M and Huang Z (2020) Computational Approaches in Preclinical Studies on Drug Discovery and Development. Front. Chem. 8:726. doi: 10.3389/fchem.2020.00726

Received

29 March 2020

Accepted

14 July 2020

Published

11 September 2020

Volume

8 - 2020

Edited by

Simone Brogi, University of Pisa, Italy

Reviewed by

Hongbin Yang, University of Cambridge, United Kingdom; Bruno Villoutreix, Institut National de la Santé et de la Recherche Médicale (INSERM), France

Updates

Copyright

*Correspondence: Zunnan Huang

This article was submitted to Medicinal and Pharmaceutical Chemistry, a section of the journal Frontiers in Chemistry

†These authors have contributed equally to this work

Disclaimer

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics