Advances in Drug Discovery 2020
Advances in Drug Discovery 2020
and Development
New Advances
Contributors
Marta Vázquez, Natalia Guevara, Marianela Lorier, Pietro Fagiolino, Iris Feria-Romero, Sandra
Orozco-Suarez, Grace Mugumbate, Lorraine Yamurai Bishi, Dr Sundeep Chaitanya Vedithi, Tom L.
Blundell, Donato Bonifazi, Angelica Intini, Giovanni Migliaccio, Carsten Wrenger, Arne Krüger,
Thales Kronenberger, Vinicius Gonçalves Maltarollo, Thuluz Meza-Menchaca, Rossana C Zepeda, Claudia
Juárez-Portilla, Mallika Alvala, Aaftaab Sethi, Khusbhoo Joshi, Sasikala K, Mohammed Chyad Al-
Noaemi, Vishwanath Gaitonde
Individual chapters of this publication are distributed under the terms of the Creative Commons
Attribution 3.0 Unported License which permits commercial use, distribution and reproduction of
the individual chapters, provided the original author(s) and source publication are appropriately
acknowledged. If so indicated, certain images may not be included under the Creative Commons
license. In such cases users will need to obtain permission from the license holder to reproduce
the material. More details and guidelines concerning content reuse and adaptation can be found at
[Link]
Notice
Statements and opinions expressed in the chapters are these of the individual contributors and not
necessarily those of the editors or publisher. No responsibility is accepted for the accuracy of
information contained in the published chapters. The publisher assumes no responsibility for any
damage or injury to persons or property arising out of the use of any materials, instructions, methods
or ideas contained in the book.
151
Countries delivered to
Top 1%
most cited scientists
12.2%
Contributors from top 500 universities
Preface III
Section 1
Introduction 1
Chapter 1 3
Introductory Chapter: The Modern-Day Drug Discovery
by Partha Karmakar, Ashit Trivedi and Vishwanath Gaitonde
Section 2
Computational Chemistry Advances in Modern Day Drug Discovery 11
Chapter 2 13
Past, Present, and Future of Molecular Docking
by Thuluz Meza Menchaca, Claudia Juárez-Portilla
and Rossana C. Zepeda
Chapter 3 27
Molecular Docking in Modern Drug Discovery: Principles
and Recent Applications
by Aaftaab Sethi, Khusbhoo Joshi, K. Sasikala and Mallika Alvala
Chapter 4 49
Computational Deorphaning of Mycobacterium tuberculosis Targets
by Lorraine Yamurai Bishi, Sundeep Chaitanya Vedithi, Tom L. Blundell
and Grace Chitima Mugumbate
Section 3
Current View towards Pharmacokinetics in Drug Discovery 69
Chapter 5 71
Revisiting Pharmacokinetics and Pharmacogenetics of Methadone
in Healthy Volunteers
by Natalia Guevara, Marianela Lorier, Marta Vázquez, Pietro Fagiolino,
Iris Feria-Romero and Sandra Orozco-Suarez
Chapter 6 85
ADME Profiling in Drug Discovery and a New Path Paved on Silica
by Arne Krüger, Vinicius Gonçalves Maltarollo, Carsten Wrenger
and Thales Kronenberger
Section 4
Recent Case Studies: Advances in Drug Discovery Research 117
Chapter 7 119
Successive Drug Therapy for a Very Rare Autosomal Diseases
by Mohammed Chyad Al-Noaemi and Hassan Ali Daghriri
Chapter 8 135
Challenges and New Frontiers in the Paediatric Drug
Discovery and Development
by Angelica Intini, Donato Bonifazi and Giovanni Migliaccio
XII
Preface
Prevention and cure cases from various diseases reach a billion individuals every
year. This is possible due to the organized efforts in the drug discovery program
across the world led by pharmaceutical companies and academic institutions.
Nevertheless, the majority of deaths globally are attributed to heart disease, stroke,
pulmonary disease, lower respiratory infections, Alzheimer’s, cancer, diabetes,
and tuberculosis. One common aspect of these disease conditions is that they are
preventable to some extent and a comparatively better quality of life can be pro-
vided to the patients with proper and timely treatment. To surmount this challenge,
the drug discovery programs have advanced significantly in the past decade. This
book, “Drug Discovery and Development – New Advances”, brings together some
of the recent progress in the field through the complex process, starting with
computational modeling right through launching a drug in the market. The book
aims to provide expert opinions and experiences from researchers encompassing
diverse key topics critical in drug discovery and development. Case studies of some
selected disease conditions will showcase the challenges and burdens faced by the
pharmaceutical industries. The book covers the latest advances in molecular dock-
ing, pharmacokinetics, and new therapeutic options of some rare diseases along
with pediatric drug discovery and development. Overall, the book represents a view
of the complex fabric of the modern-day drug discovery program.
Partha Karmakar
Research Scientist,
Washington University School of Medicine,
USA
Ashit Trivedi
Clinical Pharmacology Modeling and Simulations Group at Amgen,
USA
Section 1
Introduction
1
Chapter 1
The history of drug discovery and development is as old as some of the oldest
human civilizations. The practice of Ayurveda in India and traditional Chinese
medicines in China are over 5000-year-old therapeutic traditions that are still
in practice at large. Papyrus Ebers is evidence of medicinal practice in Egypt
about 3000 years ago [1–6]. The Greek and Roman medicines became popular
in Europe and western Asia between ~700 BC and 200 BC [7]. The ancient Arab
medicines were in practice to a great extent until 1500 AD and are still in use in the
Mediterranean gulf [8, 9]. The beginning of modern era medicine can be consid-
ered from the time when Edward Jenner discovered immunization for smallpox.
The development in the field was gradual until Sir Alexander Fleming discovered
Penicillin in 1928; since then, the field of medicinal chemistry and drug discovery
has flourished, and by the end of the twentieth century, it became a complex inter-
disciplinary platform primarily based on synthetic organic chemistry expanding
into various biological specificities [10–13]. As a result, the global pharmaceutical
market strengthened to nearly 400 billion US dollars by the year 2001 [14, 15].
3
Drug Discovery and Development - New Advances
Figure 1.
Classical components of drug discovery.
4
Introductory Chapter: The Modern-Day Drug Discovery
DOI: [Link]
Figure 2.
Modern-day symbiotic fabric of drug discovery.
reduced the library sizes for HTS and helped to understand complex multiomic data
[37, 38]. The rapid progress of different ML methods will have considerable impact
on future therapies [39].
The historical prototype for clinical drug development was to conduct a few
Phase 1 studies followed by a couple of Phase 2 studies consequently leading to
multiple expensive Phase 3 trials to demonstrate the efficacy of the drug candidate.
With the changing landscape and regulatory requirements, the number of clinical
studies to elucidate multiple questions related to drug properties such as the
mechanism of actions, pharmacokinetics (PK), pharmacodynamics (PD), and drug
metabolism increased overwhelmingly prior to Phase 3 studies. The increase in the
number of clinical trials has made drug development more lengthy and exorbitant.
To overcome this limitation and reach patients promptly, it is imperative to utilize
advanced technologies and approaches. One such approach is the PK/PD guided
drug development. PK/PD modeling has been extensively employed to generate
first-in-human dose predictions and selecting optimal doses for Phase 2 and Phase
3 trials. PK/PD modeling also plays an instrumental role in identifying if any dose
adjustments are needed in special populations such as pediatrics and geriatrics and
patients with hepatic or renal impairments [40, 41]. Additionally, PK/PD model-
informed drug development (MIDD) has gained increasing momentum in recent
years and is extensively used across pharmaceutical industries globally.
MIDD has become a crucial tool after receiving formal recognition in
Prescription Drug User Fee Act (PDUFA) VI, thus paving a path forward to opti-
mize drug dosing prior to approval and post-marketing and in special populations
5
Drug Discovery and Development - New Advances
in the absence of dedicated clinical trials. Dose optimization and clinical trial design
have been most established domains of MIDD; new technologies such as artificial
intelligence, ML, and real-world data (RWD), wearables along data science, have
the potential to transform MIDD.
ML approaches provide a set of tools that can improve decision-making for
well-specified questions with abundant, high-quality data. While using ML in the
early stages of drug design, target selection, and high-throughput screening is
almost standard today, the potential of ML during drug development has not been
recognized. The observed data/evidence obtained during the developmental phase
does not necessarily answer all the questions; thus the scope of MIDD is largely
expanded with analysis of RWD to generate real-world evidence (RWE) to resolve
these unanswered questions. Although RWD is obtained under less-controlled
settings requiring proper interpretation of the findings, it should be considered as
an attractive tool appealing for MIDD [42].
The emerging new techniques, such as portable devices, wearables, and applica-
tions (apps), may improve the dosing accuracy for patients and the quality of the
collected medical information in real-world medical practice. These tools may
improve the quality of electronic health records, making real-world data a reliable
source for drug development and dose optimization or individualization. All these
tools will make real-world data/real-world evidence a more appealing source for
MIDD [43].
Along with the power of data analytics, advances in computational chemis-
try, and new diagnostic techniques, PK/PD modeling tools have also influenced
the drug discovery research and development. These advances assist to build a
comprehensive protein-receptor database, thereby enabling a defined library size
for designing and optimization of a lead molecule. Along with the classical small
molecule drug discovery and development, many protein and antibody-based
pharmaceuticals have appeared as blockbuster drugs.
Author details
3 Chemical Research and Development, Cambrex High Point, Inc., High Point, NC,
USA
© 2020 The Author(s). Licensee IntechOpen. This chapter is distributed under the terms
of the Creative Commons Attribution License ([Link]
by/3.0), which permits unrestricted use, distribution, and reproduction in any medium,
provided the original work is properly cited.
6
Introductory Chapter: The Modern-Day Drug Discovery
DOI: [Link]
References
[1] Sharma T, Rawal G. Role of [11] Duggleby HJ, Tolley SP, Hill CP,
ayurveda in tumorigenesis: A brief Dodson EJ, Dodson G, Moody PC.
review. International Journal of Green Penicillin acylase has a single-amino-
Pharmacy (IJGP). 2012;6(2):93-101 acid catalytic centre. Nature.
1995;373(6511):264
[2] Mamtani R, Mamtani R.
Ayurveda and yoga in cardiovascular [12] Heath JR, Ribas A, Mischel PS.
diseases. Cardiology in Review. Single-cell analysis tools for drug
2005;13(3):155-162 discovery and development. Nature
Reviews Drug Discovery. 2016;15(3):204
[3] Khan S, Balick MJ. Therapeutic plants
of Ayurveda: A review of selected clinical [13] Kelm JM, Lal-Nag M, Sittampalam
and other studies for 166 species. The GS, Ferrer M. Translational in vitro
Journal of Alternative & Complementary research: Integrating 3D drug discovery
Medicine. 2001;7(5):405-515 and development processes into the
drug development pipeline. Drug
[4] Glazier A. A landmark in the Discovery Today. 2019;24(1):26-30
history of Ayurveda. The Lancet.
2000;356(9235):1119 [14] Walker N. Impacts of Glocalization
on the Pharmaceutical Industry.
[5] Farquhar J. Knowing Practice: The 2019. Available from: [Link]
Clinical Encounter of Chinese Medicine. [Link]/articles/
London: Routledge; 2018 impacts-of-glocalization-on-the-
pharmaceutical-industry
[6] Law BYK, Wu AG, Wang MJ,
Zhu YZ. Chinese medicine: A hope [15] Pharmaceutical market: Worldwide
for neurodegenerative diseases? revenue 2001-2018. 2019. Available
Journal of Alzheimer’s Disease. from: [Link]
2017;60(s1):S151-S160 statistics/263102/pharmaceutical-
market-worldwide-revenue-since-2001
[7] Ackerknecht EH, Haushofer L. A
Short History of Medicine. Baltimore: [16] Yaari Z, Da Silva D, Zinger A,
Johns Hopkins University Press; 2016 Goldman E, Kajal A, Tshuva R, et al.
Theranostic barcoded nanoparticles for
[8] Saad B, Azaizeh H, Said O. Tradition personalized cancer medicine. Nature
and perspectives of Arab herbal Communications. 2016;7:13325
medicine: A review. Evidence-based
Complementary and Alternative [17] Lee H, Lee Y, Song C, Cho HR,
Medicine. 2005;2(4):475-479 Ghaffari R, Choi TK, et al. An endoscope
with integrated transparent
[9] Bynum W, Porter R. Arab-Islamic bioelectronics and theranostic
medicine. In: Companion Encyclopedia nanoparticles for colon cancer treatment.
of the History of Medicine. London: Nature Communications. 2015;6:10059
Routledge; 2013. pp. 702-753
[18] Karmakar P, Lee K, Sarkar S,
[10] Cates J, Christie R. Subacute Wall KA, Sucheck SJ. Synthesis of a
bacterial endocarditis. A review of 442 liposomal MUC1 glycopeptide-based
patients treated in 14 Centres appointed immunotherapeutic and evaluation of
by the penicillin trials Committee of the the effect of L-rhamnose targeting on
Medical Research Council. Quarterly cellular immune responses. Bioconjugate
Journal of Medicine. 1951;20(78):93-130 Chemistry. 2015;27(1):110-120
7
Drug Discovery and Development - New Advances
[19] Karmakar P, Gaitonde V. Promising During, and After the Great Recession.
recent strategies with potential 2015. Available from: [Link]
clinical translational value to combat [Link]/cgi/[Link]?article=1
antibacterial resistant surge. Medicine. 218&context=education_etd; https://
2019;6(1):21 [Link]/education_etd/217
[20] Catalano JG, Gaitonde V, Beesu M, [27] Joshi RR, Sonavane U, Jani V,
Leivers AL, Shotwell JB. Phenoxide Saxena A, Koulgi S, Uppuladinne M,
leaving group SNAr strategy for the et al. Turbo Analytics: Applications of
facile preparation of 7-amino-3-aryl Big Data and HPC in Drug Discovery. In
pyrazolo [1, 5-a] pyrimidines from a Structural Bioinformatics: Applications
3-bromo-7-phenoxypyrazolo [1, 5-a] in Preclinical Drug Discovery
pyrimidine intermediate. Tetrahedron Process. Switzerland: Springer; 2019.
Letters. 2015;56(44):6077-6079 pp. 347-374
8
Introductory Chapter: The Modern-Day Drug Discovery
DOI: [Link]
9
Section 2
Computational Chemistry
Advances in Modern Day
Drug Discovery
11
Chapter 2
Abstract
1. Introduction
13
Drug Discovery and Development - New Advances
14
Past, Present, and Future of Molecular Docking
DOI: [Link]
As one of the most commonly used approaches since the 1980s, the experimen-
tal data obtained through molecular docking techniques have grown at an increas-
ing rate since the approach was first established. Programs configured through
different algorithms for molecular docking analysis have been developed on an
almost yearly basis, significantly improving pharmaceutical research [6]. The first
algorithms were designed for protein-protein interactions. Along with the scoring
function, which is used to determine the best binding poses, algorithms designed to
calculate the best geometrically complementary shapes as rigid bodies are necessary
to identify the most favorable orientations and conformational bindings with the
potential to confer a putative drug candidate.
15
Drug Discovery and Development - New Advances
The gradual achievement of more powerful and complex algorithms with the
addition of further parameters has paralleled computational technological advances
over the last few decades. In order to achieve optimum flexibility, in silico methods
use different tools with different approaches. Docking software depends on the
algorithms employed, which comprise three different kinds: systematic, stochastic,
or deterministic.
In the beginning, calculation algorithms that consider docking complexes to
be rigid structures were used. In rigid docking, the objective is to match the ligand
to the protein receptor, with the main aim being the generation of as many poses
as possible in order to achieve the optimum of all poses. Through this process, all
possibilities are considered heuristically to identify a group of complementary
matches that present the most favorable van der Waals forces between the ligand
and the macromolecule receptor. Intermolecular interaction calculations avoid
any flexibility but nevertheless have a level of freedom dependent on a 3x3 matrix
plus the vector rotation. This means that three rotational and three translational
degrees of freedom cover all possible moves in three-dimensional space within the
active site. However, no binding is permitted, as the macromolecular structures are
simplistically represented as solid structures located under a center of mass and
longitude [8].
The earliest work was performed using structural shape contacts, in which the
fitting of outlines enables the best possible complementary configuration between
two proteins to be identified [9]. A little later, a shape matching strategy algorithm
was used by Kuntz and collaborators in UCSF8 to continue searching for possible
configurations using the geometric distance between the ligand atoms and the
macromolecule or receptor spheres (Figure 1).
In this method, the ideal intersection or match between the ligand and recep-
tor is viewed as a “negative image” that represents the active site. The image is
produced by covering the receptor surface region and overlapping spheres with a
solvent, in which a part of the overlapping spheres comprises the actual binding
site. This constitutes the fundamentals of the DOCK search algorithm [10]. A few
years later, Kuntz also developed a more advanced approach by conferring flex-
ibility to the ligand; however, this variant is still categorized as “flexible docking.”
Figure 1.
Top left, binding site; top right, ligand. Down below conjugate with geometrical fitness functional group related
proposed by the earliest docking algorithm model.
16
Past, Present, and Future of Molecular Docking
DOI: [Link]
Subsequently, the investigation of HIV-1 protease using this approach was notable
for leading to the technique’s exponential use in drug discovery [11].
Following the pioneering work from Kuntz, a different approach was taken a
decade later in order to develop an improved new geometric recognition method,
which was developed through an algorithm called Fourier transformation [12]. For
the first time, the molecules could be described by a digital model, allowing their
interior and exterior parts to be distinguished. This novel method allows faster
calculation by determining the surface of contact, overlap, and approximation
using the six degrees of freedom. In this method, molecules are considered rigid
bodies, and the changes in structure have the degrees of freedom. This technique
makes it possible to process atomic coordinates, and Zdock represents an example
of this approach. Nevertheless, rigid-body algorithms are very erratic and ineffec-
tive in terms of any structural and conformational change arising due to the inter-
face between the ligand and the receptor. In this context, new alternatives to enable
torsions and angle movement became a matter of interest. In the same period, a new
semiflexible docking innovation was achieved using the HADDOCK protocol [13],
which involves rigid-body docking complemented by semiflexible optimization
in order to describe possible torsion angles in the main backbone and side chains.
Unlike the previous Fourier transformation method [12], which uses a grid, this
method adopts a Cartesian approach with particular coordinates, in which one of
the two molecules is flexible and the solvent can be selected. One of the two mole-
cules therefore needs to be small in order to be computationally possible in terms of
the number of conformational variations. Other methods also attempt to describe
flexible bodies undergoing rotational conformational, rotational, and translational
changes, mimicking the nature of biological molecules. In this category, both the
ligand and the receptor that are modeled by simulating protocols are flexible.
However, the flexibility needs to be lowered to make computational configuration
possible. In the end, flexible docking approaches offer a more precise technique
capable of imitating in vivo behavior of the possible structural conformations.
In flexible docking, there are two different logarithmic approaches, determinis-
tic incremental construction and stochastic. Systematic incremental construction
algorithms are most commonly used, which gradually develop binding predictions
on the basis of all possible ligand-binding poses covering all specified areas, e.g.,
DOCK [14], Glide [15], LUDI [16] FlexX [17], Hammerhead [18], and Surflex
[19], in which on-the-fly incremental ligand construction is implemented. In this
method, the number of analyses grows in line with increases in the degrees of
freedom as part of anchor-and-grow methods. In a different example, in eHiTS, the
ligand is fragmented, and each piece is tested for rigid docking, commonly based
on library screening for the best conformations to religate the fragments and test
their flexibility.
A different approach randomizes probabilistic or stochastic algorithms to
selectively reject or accept configurations through the criteria spectrum, in which
computational efforts are optimized, e.g., AutoDock [20], DARWIN [21], Monte
Carlo [22], and GOLD [23]. By the middle of the 1990s, this technique was the point
of origin of a diverse set of methods that are most commonly present in the genetic
algorithm, named after Darwin’s theory of evolution, in which the ligand is inter-
preted as a chromosome and its fragments are considered genes [24]. Every gene
exhibits conformational behavior due to its torsional/translational nature. During
computational analyses, the information is transmitted and altered through sto-
chastic crossover and mutational events evolving through specific parameters. The
changes improve the conformational binding pose from the ligand and the receptor,
e.g., Lamarckian (AutoDock). In the case of the Monte Carlo stochastic variant that
produces randomized translational conformations, the most thermodynamically
17
Drug Discovery and Development - New Advances
stable potential bindings are explored by focusing on the local minimum energy
using a decision criteria parameter that is based on a temperature reaction, called
Metropolis. The flexibility also alternates with rigid rotation, displaying several
parameters at once. A more recent development is the deterministic method, which
has been used for Newton equation simulations and also employs Monte Carlo
methods that can measure trajectories, using Amber, Charm, and GROMACS;
however, this scope forms the focus of the present work, and wide reviews have
been provided by other researchers [25–27].
The drug discovery informatics market had an estimated value of 713.4 million
USD in 2016 [28]. The presence of in silico tools that can allow the computation
of data flowing from diverse methodology pathways in parsimony with medical
chemistry can be synergistic in terms of upgrading the market and are well-known
in the scientific literature. In this manner, molecular docking has been consolidated
as a useful technique among sequence analysis platforms, molecular modeling, and
clinical training management. The use of molecular docking in each of these fields
is enhancing drug discovery in the pharmaceutical and biotechnology sector. As
it comprises several stages and workflows, the discovery of new drugs relies on in
silico tools and molecular docking in particular to simplify the overall process.
A crucial factor is the steadily rising number of structures stored in the Protein
Data Bank (PDB). The PDB is the most robust, currently storing over 151,000
structures and counting. The 3D structure information bank includes a large set
of proteins, lipids, carbohydrates, and nucleic acids, in both single structures and
complexes [29]. On the other hand, nearly a hundred different forms of molecular
docking software are available, which offer analogous implementations with vari-
ous implementation options. There has been rapid progress in developing faster
architecture based on graphics processing unit clusters, more adequate algorithms
for optimized computational analysis, and the tracking of ligand-receptor binding
expressed in scoring functions.
Although there is a need to maintain computational equipment, the associated
expenses are certainly lower than the costs of “wet lab” experiments, and molecular
docking is therefore an affordable technique. One of the most challenging tasks
in bioinformatics sciences is undoubtedly the development of new and effective
drugs, which is currently an almost mandatory step before wet lab experiments. In
structure-based drug modeling, obtaining the most accurate and efficient model of
ligand-receptor binding is a crucial step and is a suitable starting point for further
evaluation to test new compounds or drug candidates, but also and no less impor-
tantly, to discard the improbable candidates. Molecular-ligand docking is a signifi-
cant tool in pharmacology at present and an important area of drug discovery that
has comprised a central node of important achievements over the current century.
As an interdisciplinary process of multiple joint efforts mainly from the pharma-
ceutical sector, biotechnological companies, and academic researchers, as well as
many other fields, the process is highly complex and requires the most accurate and
precise tools and methodologies. This has been enhanced by an increasing number
of protein coordinates and the high number of available software programs that are
constantly evolving with more sophisticated levels and a wider field of applications,
in combination with more numerous candidates. In order to discover new drugs, as
well as improve the existing ones, it is necessary to understand the targets as well
as the nature of the possible drug candidates. In silico bioinformatics approaches
have attracted increased interest due to the results of post-genomic era sequencing.
18
Past, Present, and Future of Molecular Docking
DOI: [Link]
Figure 2.
Chart bar displaying paper publications per year (1982–2020) (NCBI, accessed on January 12, 2019).
Due to the limited set of protein-coding genes, the complexity is much higher due
to posttranscriptional modifications, prosthetic groups, multimeric complexes,
and other various phenomena, clearly demonstrating the need to better understand
their nature to fulfill biomedical objectives. Interestingly this year’s (2019) publica-
tions account for the first time a pause in the upper trend of docking publication
number (Figure 2). This may be symptomatic on how the future holds already
crucial challenges.
The drug discovery informatics market is estimated to grow from 1.5 billion
in 2016 to 2.84 billion by 2022 and may continue expanding. Accordingly, there is
currently a rising demand for the discovery and implementation of novel infor-
matics solutions. The major factors driving the expansion of the global market
include the transition from pure research to clinical treatment. More skilled
professionals, interdisciplinary backgrounds, and the high pricing of informatics
software may have a crucial impact on the growing market. At present, a number of
well-established applications have been made available for free or as paid software
or services. However, many challenges remain to be addressed to enable the full
potential of this powerful technique to be realized.
Nevertheless, in the case of pharmacology, the synergistic aspect is an important
chemical phenomenon in which two different biomolecules with different origins
can have an exponential effect in combination that is greater than their separate
effects. If it is determined that a particular structure is more favorable [30] in terms
of the docking score and it may be correlated with synergism, this can be second-
ary, due to the fact that a molecular docking procedure has not been developed to
examine it in a particular scoring function. A linear/quadratic formula could be
developed to measure synergy by discriminating between synergistic, additive, or
antagonistic effects, which can be expressed both qualitatively and quantitatively.
In this sense, further work is needed to investigate how the chemosensitivity
between a macromolecule and ligand could be detected once more than one ligand
19
Drug Discovery and Development - New Advances
20
Past, Present, and Future of Molecular Docking
DOI: [Link]
been well-established. In this sense, we can state that the sensitivity of the original
conformation of the ligands remains unanswered. Furthermore, in the case of mul-
tidomain proteins, proteins are frequently composed of more than a single effector
domain, and this should be taken into consideration.
With regard to a different aspect, how water is placed around the binding site
is not a straightforward problem to solve, although recent studies have proposed
the use of this parameter as functionally valid in specific contexts [34] within and
around the conglomerate binding site. X-ray crystallography is the most extensively
used tool for predicting 3D conformational structure; however, the actual output is
only partially informative, due to the fact that the density limits are out of resolu-
tion and, on occasion, the electron density can be of insufficient quality. Future
efforts need to endorse novel alternatives to increase the capacity and parameters
that can be used in every aspect of a given analysis, not only in terms of water but
also the physiological solutes found in nature and even protonation, in addition to
the pH potency spectra.
An understanding of the biological functions and roles of a protein in a particu-
lar cell or tissue is highly relevant in determining the role of a protein’s structure,
including all of its functional domains. Genome-wide studies have demonstrated
that multidomains are present in over 70% of eukaryotic proteins. Nevertheless,
protein-folding studies usually consider only single domains and are therefore not
focused on the mechanisms in multidomains that can even influence the folding
structure [34]. Very crucial obstacles are involved in multidomain docking analyses.
In some examples, the understanding of intermolecular movement can be restricted
by rigid docking methodologies that lack the ability to consider the effect of mul-
tiple domains in a single macromolecule. A given protein is not always present in a
static and simplistic single conformational shape but can be present in a collection
of scaffolds, stages, and intersections of conformational shapes. As a consequence,
the free energy landscape can be profoundly affected, distinctively changing the
scoring function’s output. This continues to present a major issue [35].
To improve modeling, the role played by multiple molecules in the context of
a certain reaction is an indispensable step that must be considered. At the current
stage of technology, this does not fall under the current scope of molecular docking,
due to the fact that the processes are far too complex and it is difficult to manage
all of the interactions that occur during a molecular binding and reaction. In order
to mimic how chemistry works in nature, the inclusion of more than two factors
(ligand/macromolecule) where methodologically possible would be a priority to
enable the possible interactions in a molecular group to be predicted. Although a
few software packages use this approach, in the future, it needs to become more
common in other methods to address the binding modes of ligands in assessments
with higher stoichiometry using multiple ligand complexes against the molecular
target. Additionally, as stated earlier in this work, it would be of great interest to
evaluate the synergy of ligand combination conjugates.
5. Conclusion
Over the last four decades, molecular docking has improved quite remarkably,
contributing to the enhancement and improvement of pharmacology in addition
to many different areas of applied and molecular biology. After the first complete
draft of the Human Genome Project was announced in 2003, the scientific com-
munity concluded that there are far fewer protein-coding genes than expected and
it has therefore been swift to study how molecules interact by investigating more
possible target bindings of a given molecule. The increasing demand for molecular
21
Drug Discovery and Development - New Advances
Author details
© 2020 The Author(s). Licensee IntechOpen. This chapter is distributed under the terms
of the Creative Commons Attribution License ([Link]
by/3.0), which permits unrestricted use, distribution, and reproduction in any medium,
provided the original work is properly cited.
22
Past, Present, and Future of Molecular Docking
DOI: [Link]
References
23
Drug Discovery and Development - New Advances
and assessment of docking accuracy. [24] Oshiro CM, Kuntz ID, Dixon JS.
Journal of Medicinal Chemistry. Flexible ligand docking using a genetic
2004;47(7):1739-1749 algorithm. Journal of Computer-Aided
Molecular Design. 1995;9(2):113-130
[16] Bohm HJ. The computer program
LUDI: A new method for the de novo [25] Salomon-Ferrer R, Case DA,
design of enzyme inhibitors. Journal Walker RC. An overview of the Amber
of Computer-Aided Molecular Design. biomolecular simulation package.
1992;6(1):61-78 Wiley Interdisciplinary Reviews:
Computational Molecular Science.
[17] Kramer B, Rarey M, Lengauer T. 2013;3(2):198-210
Evaluation of the FLEXX incremental
construction algorithm for [26] Brooks BR, Brooks CL 3rd,
protein–ligand docking. Proteins: Mackerell AD Jr, Nilsson L, Petrella RJ,
Structure, function. Bioinformatics. Roux B, et al. CHARMM: The
1999;37(2):228-241 biomolecular simulation program.
Journal of Computational Chemistry.
[18] Welch W, Ruppert J, Jain AN. 2009;30(10):1545-1614
Hammerhead: Fast, fully automated
[27] Berendsen HJC, van der Spoel D,
docking of flexible ligands to protein
binding sites. Chemistry & Biology. van Drunen R. GROMACS - a
1996;3(6):449-462 message-passing parallel molecular-
dynamics implementation.
[19] Spitzer R, Jain AN. Surflex-
Computer Physics Communications.
dock: Docking benchmarks and 1995;91(1-3):14
real-world application. Journal of
[28] Drug Discovery Informatics
Computer-Aided Molecular Design.
Market Size, S. T. A. R. Drug Discovery
2012;26(6):687-699
Informatics Market Size, Share &
Trends Analysis Report by Workflow
[20] Morris GM, Huey R, Lindstrom W,
(Discovery, Development), by Mode
Sanner MF, Belew RK, Goodsell DS,
(In-House, Outsourced), by Services,
et al. AutoDock4 and AutoDockTools4:
by Region, Vendor Landscape, and
Automated docking with selective
Segment Forecasts 2018-2025. 2016.
receptor flexibility. Journal
978-1- 68038-749-0
of Computational Chemistry.
2009;30(16):2785-2791 [29] Berman HM, Bhat TN, Bourne PE,
Feng Z, Gilliland G, Weissig H, et al.
[21] Taylor JS, Burnett RM. DARWIN: A The protein data bank and the challenge
program for docking flexible molecules. of structural genomics. Nature
Proteins. 2000;41(2):173-191 Structural Biology. 2000;7:957-959
24
Past, Present, and Future of Molecular Docking
DOI: [Link]
25
Chapter 3
Abstract
The process of hunt of a lead molecule is a long and a tedious process and one is
often demoralized by the endless possibilities one has to search through. Fortu-
nately, computational tools have come to the rescue and have undoubtedly played a
pivotal role in rationalizing the path to drug discovery. Of all techniques, molecular
docking has played a crucial role in computer aided drug design and has swiftly
gained ranks to secure a valuable position in the modern scenario of structure-based
drug design. In this chapter, the principle, sampling algorithms, scoring functions
and diverse available software’s for molecular docking have been summarized. We
demonstrate the interplay of docking, classical techniques of structure-based design
and X-ray crystallography in the process of drug discovery. In addition, we dwell
upon some of the limitations faced in docking studies. Finally, several success
stories of molecular docking approaches in drug discovery have been highlighted,
concluding with remarks on molecular docking for the future.
1. Introduction
The path to drug discovery is a long, challenging & arduous task not to mention
the overburdening finances it demands. As of 2014, the average cost of developing a
new drug from scratch was found to be an estimated $2.5 billion, an increase of
145% from the previous study done by the same organization in 2003. The major
reasons for this drastic increase in the cost is mainly attributed to high failure rate of
drugs among others [1]. Understanding of the drug discovery process is important
to handle the challenges faced by the pharma companies in terms of cost and
innovation.
The process of identifying a target, synthesizing an active compound with
suitable characteristics like minimal toxicity, high bioavailability, cost-effective
synthesis, etc., and finally developing it to introduce in the market is a time-
consuming, extremely complex and risky endeavor [2]. Initially, a target is identi-
fied which plays a key role in progress of the disease. Once a link between the target
and the disease has been established, the next step is to identify potential candidates
which can stop or reverse the progress of the disease [3]. This process starts with the
27
Drug Discovery and Development - New Advances
discovery of molecules that show efficacy in a simple screen, called “hits.” Screen-
ing is a process in which normally a large number of compounds from natural
products and online databases are examined for biological activity in high-
throughput assays. This step in the drug discovery process is very crucial and
demands maintaining huge molecular libraries and carrying out thousands or mil-
lions of assays, which leaves the academicians and small pharmaceutical companies
at a disadvantage and also shoots up the cost for larger industries. Next, the “hits”
found are chemically modified to give improved pharmaceutical properties, such
compounds are often called “leads.” But, it is quite apparent that the method stated
above for discovery of a drug has a number of pitfalls. From an academic point of
view, carrying out high throughput screens (HTS) is costly, time-consuming and
not feasible; while, from an industrial perspective, it does nothing to improve the
eminent danger of market saturation.
Truly innovative and blockbuster drugs are what drive the pharmaceutical
industry forward but, over the past few years introduction of new molecular enti-
ties (NMEs) has vastly reduced. For example, in 2007 only 19 NMEs were approved
by the US Food and Drug Administration (FDA), the least since 1983 [4]. Cur-
rently, and even in the future it is expected that only slight modifications of the
existing blockbuster drugs would be carried out which would further aggravate this
problem [5]. HTS would not help in either curbing the rising costs of discovering
hits or the problem of finding truly innovative and blockbuster NMEs, the two
major hurdles that the pharmaceutical industry faces now-a-days.
To overcome these challenges, molecular docking is an exemplary tool. During
the first step to find hits from existing chemicals for a drug discovery and develop-
ment project, virtual screening (VS) is a perfectly viable and an alternative or
complementary approach to HTS for fulfilling the screening of thousands or mil-
lions of compounds within a few days. In addition, the speed of VS helps in kick-
starting projects for newer targets for which no leads are available [6]. Molecular
docking is one of the most applied virtual screening methods and has become
increasingly useful overtime on account of immense growth in 3D X-ray and NMR
structures and their improved resolution (physics and knowledge based docking
algorithms depend on it) reported in the Protein Data Bank (PDB). As an example,
in total 46,541 X-ray structures were reported at the end of 2008 in PDB, but by the
end of 2018 it had grown to a staggering figure of 131,993 [7]. In addition, it is a
resource saving technique which provides accessibility of screening to academia and
small industries which were earlier limited to large pharmaceutical giants.
In this chapter, we will discuss a particular class of molecular design, i.e.,
“Docking” along with the various algorithms, techniques, success stories and limi-
tations related to it. In the end, we will conclude with its scope in the near future.
2. Molecular docking
Two molecules can interact in a number of ways let alone the interaction of a
protein and protein or a protein and small molecule. Molecular docking helps us in
predicting the intermolecular framework formed between a protein and a small
molecule or a protein and protein and suggest the binding modes responsible for
inhibition of the protein. To accurately carry out docking studies one requires the
high-resolution X-ray, NMR or homology-modeled structure with known/predicted
binding site in the biomolecule. To date, 148,827 are available in the database (PDB)
[3]. Docking methods fit a ligand into a binding site by combining and optimizing
variables like steric, hydrophobic and electrostatic complementarity and also
estimating the free energy of binding (scoring) [8].
28
Molecular Docking in Modern Drug Discovery: Principles and Recent Applications
DOI: [Link]
There are two basic components which distinguish the variety of docking soft-
wares available to choose from—One is sampling algorithm and the other is scoring
function, these are discussed in detail.
As pointed out earlier, there are a huge number of modes of binding between
two molecules and even with advances in parallel computing and higher clock
speeds of modern computers it would be expensive and time-consuming to generate
all the possible modes. Therefore, algorithms were needed which could fish out
valuable conformations from the fruitless ones.
Various algorithms were developed in this regard and can be classified by the
number of degrees of freedom they ignore. The simplest of the algorithms intro-
duced treated the molecules as two rigid bodies thereby reducing the degree of
freedom to just six (three translational and three rotational). A very well cited
example of a program using this algorithm is DOCK [9]. This program was designed
to find molecules which had a huge extent of shape similarity to the pockets/
grooves or binding sites. It derives an image of suspected binding sites present on
the surface of the protein. This image consists of several overlapping spheres of
varying radii which touch the molecular surface of the macromolecule at just two
points. The ligand molecule is also considered as a set of spheres which approxi-
mately fill the space occupied by the ligand. Once the respective representations of
the protein surface and the ligand in terms of sphere are complete, the pairing rule
is applied. The pairing rule is based on the principle that ligand sphere can be paired
with a protein sphere if the internal distances of all the spheres in the ligand set
match all the internal distances within the protein set, allowing some user specified
tolerance. Thus, it allows the program to identify geometrically similar cluster of
spheres on the protein site and the ligand. Many other programs were developed
later which make use of such matching algorithm (MA) which include LibDock [8],
LIDAEUS [10], PhDOCK [11], Ph4DOCK [12], Q-fit [13], SANDOCK [14], etc.
All these programs based on MA have the advantage of speed but have several
limitations such as prior need for detailed receptor geometry and lack of molecular
flexibility which does not accurately define many aspects of ligand-protein
interactions.
The second algorithm is that of incremental construction (IC), wherein the
ligand is fragmented from rotatable bonds into various segments. One of the seg-
ments is anchored to the receptor surface. The anchor is generally considered to be
the fragment which shows maximum interactions with the receptor surface, has
minimum number of alternate conformations and fairly rigid such as the ring
system. Once the base/anchor has been established, the next step is to add each of
the fragments step by step. Ideally, those fragments are added first which have a
greater chance of showing interactions like hydrogen bonding since they are direc-
tional in nature and are responsible for specificity of the ligand. In addition, hydro-
gen bonds lead to more accurate prediction of geometry. Once a particular fragment
is added, the poses with the least energies are considered for the next iteration,
making the algorithm extremely fast and robust [15]. IC has been used in programs
like DOCK 4.0 [16], FlexX [15], Hammerhead [17], SLIDE [18] and eHiTS [19],
SKELGEN [20], ProPose [21], PatchDock [22], MacDock [23], FLOG [24], etc. One
major limitation of this program is that it is restricted to medium sized ligands and
is not feasible for large size ligands where the number of fragments generated pose a
problem.
Another useful algorithm is the Monte Carlo (MC) technique. In this approach, a
ligand is modified gradually using bond rotation and translation or rotation of the
29
Drug Discovery and Development - New Advances
entire ligand. More than one parameter can also be changed at a time to get a
particular conformation. That conformation is then evaluated at the binding site
based on energy calculation using molecular mechanics and is then rejected or
accepted for the next iteration based on Boltzmann’s probability constant. Accep-
tance or rejection of the conformation is a function of the change in energy with
respect to a parameter T, which can be physically interpreted as temperature (sim-
ulated annealing). This criterion of acceptance or rejection makes this method
strikingly different than the others. Whereas the other algorithm favor decrease in
energy, in MC method increases are also possible. For higher values of T increases
are likely. If one starts at a high value of T, then small energy barriers can be jumped
and the configuration can move beyond local minima and is therefore particularly
useful in situations where a global minimum is sought among many local minima
[25]. An interesting spin-off of the MC approach is the Tabu search, which main-
tains a record of the search space of the binding site which has already been visited
and thus ensures that the binding site is explored to the maximum [26]. MC
approach has been made use of in programs like DockVision 1.0.3 [25], FDS [27],
GlamDock [28], ICM [29], MCDOCK [30], PRODOCK [31], QXP [32],
ROSETTALIGAND [33], RiboDock [34], Yucca [35], AutoDock [36], etc. One of
the major concerns with MC approach is the uncertainty of convergence, which can
be improved by performing multiple independent runs.
Genetic algorithm (GA) is quite similar to MC method and is basically used to
find the global minima [37]. These are much inspired by the Darwin’s Theory of
Evolution [38]. GA maintains a population of ligands with an associated fitness
determined by the scoring function. Each ligand represents a potential hit. The GA
alters the ligands of the population by mutation or crossover. In the first stage, a
new population is created by accessing and then selecting the more fit ligands from
the previous step. The members of the populations are then transformed in the
alteration step. The mutation operator creates new ligands from a single ligand by
randomly changing a fragment in its representation while the crossover operator
exchanges information between two (occasionally more) members of the popula-
tion [39–41]. GA has been incorporated in programs like Autodock 4.0 [42], DAR-
WIN [43], DIVALI [39], FITTED [44], FLIPDock [45], GAMBLER [46], GAsDock
[47], GOLD 3.1 [48], PSI-DOCK [49]. GA has a similar limitation like that of MC
method wherein the uncertainty of convergence is a major drawback.
Another approach is the hierarchical method. In this approach, the low energy
conformations of the ligand are pre-computed and aligned. The populations of the
pre-generated ligand conformations are merged into a hierarchy such that similar
conformations are positioned adjacent to each other within the hierarchy. After-
wards, on carrying out rotation or translation of the ligand, the docking program
will make use of this hierarchical data structure and thus minimize the outcomes.
Let us understand with a simple example—if an atom near the rigid center of the
ligand is found to clash with the protein in a given rotation/translation, then this
approach can reject all of the conformations lying below in the hierarchy to that of
the conformation under scrutiny, because the descendants must contain the same
clash as well [50]. GLIDE software makes use of the hierarchical method [51, 52].
30
Molecular Docking in Modern Drug Discovery: Principles and Recent Applications
DOI: [Link]
Here, rij stands for the distance between protein atom i and ligand atom j, Aij and
Bij are the van der Waal parameters, qi and qj are the atomic charges and ε(rij) is the
distance dependent dielectric constant.
One common example of a FF scoring function is that of the program DOCK
[56] represented in Eq. (1), where, the effect of solvent is indirectly considered by
the distance dependent dielectric constant e(rij) seen in the Coulombic potential.
One major drawback of this function is that it does not consider an important
solvent effect that charged groups favor aqueous environments whereas non-polar
groups tend to stay in non-aqueous environments, commonly referred to as the
desolvation effect [57]. Ignorance could lead to biased results as the function would
now be totally dependent on the coulombic interactions and would thus favor
highly charged ligands. In other words, it only takes into account the interaction of
protein and ligand, which is inadequate. To build a more robust function one needs
to also evaluate how both interact with water before the formation of the complex
and how water mediates this process.
Later the Shoichet group [58] improved upon the existing function by adding the
effects of the solvent on protein-ligand interactions using implicit solvent models.
They employed the Poisson-Boltzmann approach to model the electrostatic potential
of the protein. The van der Waals interactions were calculated using the Lennard-
Jones potential; the electrostatic interaction between the ligand and the protein was
estimated using a precomputed receptor potential determined with DelPhi [59].
Ligand desolvation penalties were calculated with HYDREN [60]. The solvent-
corrected scores were found to be closer to experimental binding free energies than
the DOCK program scores, but were still too favorable. The overestimation of
complex stability could be due to the neglect of solute entropic terms [58].
There a few scoring functions which be classified in this category such
as DockScore [56], GoldScore [61], HADDOCK Score [62], ICM SF [29],
QXP SF [32], etc.
31
Drug Discovery and Development - New Advances
Empirical scoring functions: the basis of this scoring function is that the binding
energies of a complex can be approximated by a sum of individual uncorrelated
terms. The coefficients of the various terms involved in calculation of binding
energy are obtained from regression analysis using experimentally determined
binding energies or potentially from X-ray structural information. Empirical func-
tions have simpler energy terms to evaluate when compared to force field scoring
functions and thus are much faster in binding score calculations.
The first empirical scoring function developed to predict binding free energies
was implemented in LUDI, credited to the pioneering work of Bohm [63]. The
energy was derived using experimental binding free energies and protein-ligand
crystal structures for 45 complexes.
ΔGbind ¼ ΔGO þ ΔGhb ∑ f ðΔR; ΔαÞ þ ΔGionic ∑ f ðΔR; ΔαÞ þ ΔGlipo Alipo
h�bonds ionic int:
þ ΔGrot NROT:
(2)
32
Molecular Docking in Modern Drug Discovery: Principles and Recent Applications
DOI: [Link]
density in the reference state where interatomic interactions are zero and g(r) is
pair distribution function.
Popular knowledge based functions include DrugScore [79], PMF [72, 80],
MScore [81], SMoG [71], BLEEP [74], ITScore/SE [75], etc. The computational
simplicity of such functions is a major advantage especially when one has large
databases at hand however, the accuracy of predicting the reference state and
underrepresentation of interactions with halogens and metals are the major hurdles.
Each of the above classified have their inherent drawbacks, in this regard,
combination of more than one scoring functions has given improved results. This
approach has been widely regarded as “Consensus Scoring” [46].
Another set of scoring functions which have recently started to attract attention
are based on machine learning. One of the programs based on functions incorpo-
rating machine learning was able to achieve an astounding hit rate of 88.6% [82].
The nexus of machine learning and scoring functions is promising but the develop-
ment of such a tool is slow owing to its complexity.
In order to compare the variety of scoring functions that have been developed
up until now, comparative assessment of scoring functions (CASF) is an incredible
platform to begin with [83].
There is another set of classification proposed for the scoring functions namely
physics-based methods, empirical scoring functions, knowledge-based potentials,
and descriptor-based scoring functions but there is still no clear consensus on which
classification of scoring functions would be appropriate [84].
3. Applications
Molecular docking has been developed and improving for many years, but its
ability to generate a viable drug is still generally questioned. In the section below,
you will find examples where docking approach lead to recognition of active hits for
a variety of different receptors/targets.
HIV 1 Integrase—a new binding site for drugs treating AIDS was discovered
by Schames et al. using docking while considering the flexibility of the receptor
through molecular dynamics. The group used AutoDock in conjunction with
the relaxed-complex method to discover novel mode of inhibition of HIV
integrase [85].
α1A Adrenergic receptor—Evers et al. generated a model of the receptor using
homology modeling based on the X-ray crystallographic structure of bovine rho-
dopsin. Hierarchical virtual screening method was performed by them on the
Aventis in-house compound repository in a stepwise manner. 22,950 filtered com-
pounds were then docked into the α1A receptor homology model with the program
GOLD and scored with PMF. The top scoring compounds were finally clustered
according to their unity fingerprint similarity, and a diverse set of 80 compounds
was tested in a radio ligand displacement assay. Thirty-seven compounds displayed
a Ki < 10 μM with the most active having Ki = 1.4 nM [86].
Type I TGF-beta receptor kinase—A striking example and a proof of the benefit
of in silico approach over classical high-throughput screening involves the discov-
ery of novel Type I TGF-beta receptor kinase inhibitor. The same molecule (HTS-
466284); Figure 1, a 27 nM inhibitor, was discovered independently using virtual
screening [87] and also by traditional enzyme and cell-based high-throughput
screening in the same year [88]. The compound discovered experimentally required
in vitro screening of a large library of compounds in a TGF-β-dependent cell-based
assay which required more time, proved to be costlier and required usage of a
variety of chemicals when compared to its computational counterpart.
33
Drug Discovery and Development - New Advances
Figure 1.
Structure of HTS-46628, type I TGF-beta receptor kinase inhibitor.
Figure 2.
Structures for Aurora Kinase A inhibitor with IC50 12 and 43 pM respectively.
Aurora Kinase A—A major improvement was seen in the inhibitory activity of
Aurora Kinase A inhibitors which were designed using in silico techniques by Park
et al. [89]. This research group made use of a genetic algorithm to carry out the
sampling while the scoring function involved the energy terms from the AutoDock
program with a slight modification of the dehydration energy term. The design
strategy and tools used to carry out the study proved to be immensely successful
with some inhibitors revealing exceptionally high potency at low picomolar levels;
Figure 2 [89].
Dopamine D3 receptor—The 3D structure of the Dopamine 3 (D3) subtype
receptor was modeled by Varady et al. from the X-ray crystallographic structure of
rhodopsin and validated using experimental data. A D3 pharmacophore model was
devised by them from 10 selective and potent known D3 receptor ligands. Using
their model, 250,251 compound were screened from the National Cancer Institute
(NCI) 3D database. The hit list of 2478 potential ligands was then filtered for known
chemotypes. After removal of all compounds that were structurally similar to
known D3 receptor ligands, 1314 candidates remained. At the end, 20 compounds
supplied by NCI to the group were tested, out of which eight had Ki values below
500 nM, among which one of the compounds had Ki = 11 nM; Figure 3 [90].
Serotonin receptor (5HT1A)—Due to lack of structural information available for
the receptor, Becker et al. made use of PREDICT, to develop a unique non-
homology model for building a virtual 3D structure of the receptor. Using the
model, 40,000 compounds from Predix’s compound library were screened for
molecular docking and 78 virtual hits were discovered and then purchased by them
from respective vendors. The in vitro 5-HT1A binding assays elucidated that 16 of
the 78 compounds tested by the group were found to be hits with Ki < 5 μM,
reflecting a 21% hit rate, 9 of which had a Ki < 1 μM. The most potent molecule had
Ki = 1 nM (Figure 4) and was selected as a lead molecule for further optimization.
One significant feature of the study which highlights the utility of docking was that
34
Molecular Docking in Modern Drug Discovery: Principles and Recent Applications
DOI: [Link]
Figure 3.
Structure of dopamine D3 receptor inhibitor with Ki = 11 nM.
Figure 4.
Structure of serotonin receptor inhibitor with Ki = 1 nM.
the complete discovery process, i.e., from in silico screening through lead optimi-
zation, preclinical, and into clinical studies, was very rapid, requiring less than a
couple of years from program initiation to Phase I clinical trial [91].
Crystal structure prediction challenge—The International Blind Test is a chal-
lenge organized by the Cambridge Crystallographic Data Center wherein a previ-
ously determined crystal structure is only revealed once all the participants submit
their respective structures. In the Fifth International Blind Test, the challenge was
toughened by including flexible molecules with 50–60 atoms. The successful pre-
diction by two participants of the crystal structure of molecule XX in the blind test
indicated that search methods and models for lattice energy are capable of provid-
ing worthwhile results, both in terms of the range of structures considered in the
search and relative energies of the structures and thus can act as efficient ranking
systems [92].
Muscarinic M3 receptor—A pharmacophore model was constructed by Marriot
et al. from the known molecules showing significant M3 potency [93]. The research
group utilized the program DISCO, which generated five models. Three models
were rejected based on structural overlay. 3D screening was performed by Unity 3D
of the Astra compound database. The first model developed by them gave 176 hits
while the second model gave 173 hits; 172 compounds were common to the two sets
and were tested for their M3-antagonistic potency. Several compounds with
35
Drug Discovery and Development - New Advances
Figure 5.
Structure of muscarinic M3 receptor antagonist.
Figure 6.
Structures of checkpoint kinase 1 inhibitor with IC50 450 nM and 4 μM respectively.
36
Molecular Docking in Modern Drug Discovery: Principles and Recent Applications
DOI: [Link]
variety of functional groups was reduced to a data set of just 343 test compounds.
Molecular docking was performed by them and the top scoring poses of the
GoldScore ranking list were taken into account for the manual selection of the
virtual hits based on visual inspection of the appropriate fit of the molecule in the
active site. A data set of 44 compounds including the five low scoring compounds
were finally selected for experimental evaluation. The activity of 21 out of the
selected 39 in silico hits was experimentally confirmed and four out of the five
structures predicted as inactive showed no activity on cathepsin K. This study
demonstrated to a huge extent the ability of docking to generate positive outcomes
(Figure 7) [95].
Human aldose reductase (ALR2)—ALR2 catalyzes a key reaction in the polyol
pathway of glucose metabolism, a process implicated in the long-term complica-
tions of diabetes. Its inhibitors were designed by Wang et al. using molecular
dynamic (MD) simulations and virtual screening [96]. A major challenge encoun-
tered by them in the in silico studies was that the binding site of the enzyme
underwent large conformational changes and adopted distinct configurations upon
binding different classes of ligands. To address this issue, the group sampled poten-
tially accessible binding site conformations by MD simulations based on the avail-
able crystallographic structures of ALR2. After this procedure, three average
conformations were selected for the docking. FlexX was utilized to carry out
docking of 7200 compounds of which 128 compounds were selected by them for
further screening. Out of these 72 molecules were selected which had RMSD < 3.00
A for experimental assay, of which 15 novel ALR2 inhibitors hits were discovered.
The most potent inhibitor had an IC50 = 0.24 μM; Figure 8 [96].
Cyclooxygenase-2 (COX-2) and β-amyloid aggregation inhibitors—Dadashpour
et al. made use of AutoDock4.2 to carry out docking studies of designed molecules
Figure 7.
Respective structures for active and inactive covalent binders of human cathepsin K.
Figure 8.
Structure of human aldose reductase inhibitor with IC50 = 0.24 μM.
37
Drug Discovery and Development - New Advances
Figure 9.
Structure of cyclooxygenase-2 inhibitor with IC50 = 10.1 μM.
4. Limitations
The major limitation of molecular docking is due to the lack of confidence on the
ability of scoring functions to give accurate binding energies. This stems from the
fact that some intermolecular interaction terms are hardly predicted accurately,
such as solvation effect and entropy change [98]. In addition, some intermolecular
interactions are rarely considered in scoring functions which have been proven to
be of significance. For instance, halogen bonding is verified to make a contribution
to protein-ligand binding affinity [99] and so do guanidine-arginine interactions
[100], but are not considered.
Transthyretin-thyroxine complex—One critical example wherein energy func-
tions failed is that of transthyretin-thyroxine complex. The docking simulations
with energy functions resulted in generation of two binding modes, one similar to
the native binding mode of thyroxine and the other belonging to an alternate
binding domain with a root mean square deviation (RMSD) of 8.97 Å from native
binding state. The energy simulation was carried out and the lower energy solution
picked by the docking program was the one with higher RMSD. Thus, in this case
molecular docking failed to make the correct prediction of binding mode. Thereby,
it would be fair to conclude that we might get many false negatives during the
process of VS. [101].
It is still an unsolved problem to accurately deal with the water molecules in
binding pocket during docking process, which is tough task and needs a lot of
attention in the near future due to two reasons. Firstly, the x-ray crystal structures
lack the coordinate information of hydrogen, due to inefficient scattering by
smaller atoms. Not knowing the exact position of hydrogen leads to inaccuracies in
identifying water molecules which might be acting as a bridging molecule between
the ligand and the receptor. Secondly, no reliable theoretical approach is available
to accurately predict how water molecules are affected by ligands and how strong
the effect is. On top of that, it impossible with our current knowledge to predict
how many water molecules in the binding pocket would be replaced by potential
ligands and how the hydrogen bonding network would be disturbed by ligand
binding [102].
38
Molecular Docking in Modern Drug Discovery: Principles and Recent Applications
DOI: [Link]
One of the major challenges faced in the field of docking is that of rigid receptor.
A protein can adopt many different conformations depending upon the ligand to
which it binds. As a result, docking performed using a rigid receptor will corre-
spond to a single receptor conformation, which leads to false negatives in many
cases where later the ligand was found to be active. This happens because a protein
can exist in constant motion between different conformational states having similar
energies, which is usually neglected in docking [58].
Finally, the spectrum of activity against off-target proteins is something rarely
seen even in computational screens and is only dealt by animal and human trials.
5. Conclusion
Thus, it is quite evident from the case studies highlighted above and many more
success stories that one can find in literature related to computer aided drug design,
that in silico approaches in combination with biophysical data, experimental high
throughput screening and biology/toxicology/clinical studies are an indispensable
tool in the process of drug discovery. It assists in decision making, conceptualizing
new ideas and exploring them in a rapid manner to test a hypothesis, bringing
solutions to problems that cannot be assessed experimentally either because the
experiments is too difficult to design or because it would cost too much.
Undoubtedly, many challenges still remain to be addressed such as role of water
molecules, solvent effects, entropic effects, and receptor flexibility.
There is more than sufficient information now that proves the utility of compu-
tational tools in drug design and there is no scope for any debate regarding the
effectiveness and advantage of computational tools in the process of drug discovery.
Author details
© 2019 The Author(s). Licensee IntechOpen. This chapter is distributed under the terms
of the Creative Commons Attribution License ([Link]
by/3.0), which permits unrestricted use, distribution, and reproduction in any medium,
provided the original work is properly cited.
39
Drug Discovery and Development - New Advances
References
[1] Tufts Center for the Study of Drug [10] Wu SY, McNae I, Kontopidis G,
Development. Cost to develop and win McClue SJ, McInnes C, Stewart KJ, et al.
marketing approval for a new drug is Discovery of a novel family of CDK
$2.6 billion; 2014 inhibitors with the program LIDAEUS:
Structural basis for ligand-induced
[2] Muntha P. Drug discovery & disordering of the activation loop.
development. Journal of Pharmacy and Structure. 2003;11(4):399-410. DOI:
Pharmaceutical Sciences. March, 2016;5 10.1016/S0969-2126(03)00060-1
(1):135-142
[11] Joseph-McCarthy D, Thomas BE IV,
[3] Available from: [Link]
Belmarsh M, Moustakas D, Alvarez JC.
ov/forpatients/approvals/drugs/uc
Pharmacophore-based molecular
[Link]
docking to account for ligand flexibility.
[4] Paul SM, Mytelka DS, Dunwiddie Proteins: Structure, Function, and
CT, Persinger CC, Munos BH, Lindborg Bioinformatics. 2003;51(2):172-188.
SR, et al. How to improve R&D DOI: 10.1002/prot.10266
productivity: The pharmaceutical
industry’s grand challenge. Nature [12] Goto J, Kataoka R, Hirayama N.
Reviews Drug Discovery. 2010;9(3):203 Ph4Dock: Pharmacophore-based
protein�ligand docking. Journal of
[5] Karamehic J, Ridic O, Ridic G, Jukic Medicinal Chemistry. 2004;47(27):
T, Coric J, Subasic D, et al. Financial 6804-6811. DOI: 10.1021/jm0493818
aspects and the future of the
pharmaceutical industry in the United [13] Jackson RM. Q-fit: A probabilistic
States of America. Materia Socio method for docking molecular
Medica. 2013;25(4):286. DOI: 10.5455/ fragments by sampling low energy
msm.2013.25.286-290 conformational space. Journal of
Computer-Aided Molecular Design.
[6] Alvarez JC. High-throughput 2002;16(1):43-57. DOI: 10.1023/A:
docking as a source of novel drug leads. 1016307520660
Current Opinion in Chemical Biology.
2004;8(4):365-370. DOI: 10.1016/j.
[14] Burkhard P, Taylor P, Walkinshaw
cbpa.2004.05.001
M. An example of a protein ligand found
[7] Available from: [Link]
by database mining: Description of the
org/stats/growth/xray\ docking method and its verification by a
2.3 Å X-ray structure of a thrombin-
[8] Diller DJ, Merz KM Jr. High ligand complex1. Journal of Molecular
throughput docking for library design Biology. 1998;277(2):449-466. DOI:
and library prioritization. Proteins: 10.1006/jmbi.1997.1608
Structure, Function, and Bioinformatics.
2001;43(2):113-124. DOI: 10.1002/ [15] Rarey M, Kramer B, Lengauer T,
1097-0134(20010501)43:2<113::AID- Klebe G. A fast flexible docking method
PROT1023>[Link];2-T using an incremental construction
algorithm. Journal of Molecular Biology.
[9] Kuntz ID, Blaney JM, Oatley SJ, 1996;261(3):470-489. DOI: 10.1006/
Langridge R, Ferrin TE. A geometric jmbi.1996.0477
approach to macromolecule-ligand
interactions. Journal of Molecular [16] Ewing TJ, Makino S, Skillman AG,
Biology. 1982;161(2):269-288. DOI: Kuntz ID. DOCK 4.0: Search strategies
10.1016/0022-2836(82)90153-X for automated molecular docking of
40
Molecular Docking in Modern Drug Discovery: Principles and Recent Applications
DOI: [Link]
[18] Schnecke V, Kuhn LA. Virtual [25] Hart TN, Read RJ. A multiple-start
screening with solvation and ligand- Monte Carlo docking method. Proteins:
induced complementarity. In: Virtual Structure, Function, and Bioinformatics.
Screening: An Alternative or 1992;13(3):206-222. DOI: 10.1002/
Complement to High Throughput prot.340130304
Screening? Netherlands: Springer; 2000.
pp. 171-190. DOI: 10.1007/0-306- [26] Baxter CA, Murray CW, Clark DE,
46883-2_10 Westhead DR, Eldridge MD. Flexible
docking using Tabu search and an
[19] Zsoldos Z, Reid D, Simon A, Sadjad empirical estimate of binding affinity.
BS, Peter Johnson A. eHiTS: An Proteins: Structure, Function, and
innovative approach to the docking and Bioinformatics. 1998;33(3):367-382.
scoring function problems. Current DOI: 10.1002/(SICI)1097-0134
Protein and Peptide Science. 2006;7(5): (19981115)33:3<367::AID-PROT6>3.0.
421-435. DOI: 10.2174/13892030677 CO;2-W
8559412
[27] Taylor RD, Jewsbury PJ, Essex JW.
[20] Alberts IL, Todorov NP, Dean PM. FDS: Flexible ligand and receptor
Receptor flexibility in de novo ligand docking with a continuum solvent
design and docking. Journal of model and soft-core energy function.
Medicinal Chemistry. 2005;48(21): Journal of Computational Chemistry.
6585-6596. DOI: 10.1021/jm050196j 2003;24(13):1637-1656. DOI: 10.1002/
jcc.10295
[21] Seifert MH. ProPose: Steered virtual
screening by simultaneous protein- [28] Tietze S, Apostolakis J. GlamDock:
ligand docking and ligand-ligand Development and validation of a new
alignment. Journal of Chemical docking tool on several thousand
Information and Modeling. 2005;45(2): protein-ligand complexes. Journal of
449-460. DOI: 10.1021/ci0496393 Chemical Information and Modeling.
2007;47(4):1657-1672. DOI: 10.1021/
[22] Schneidman-Duhovny D, Inbar Y, ci7001236
Nussinov R, Wolfson HJ. PatchDock
and SymmDock: Servers for rigid and [29] Totrov M, Abagyan R. Flexible
symmetric docking. Nucleic Acids protein-ligand docking by global energy
Research. 2005;33(suppl_2):W363- optimization in internal coordinates.
W367. DOI: 10.1093/nar/gki481 Proteins: Structure, Function, and
Bioinformatics. 1997;29(S1):215-220.
[23] Fradera X, Kaur J, Mestres J. DOI: 10.1002/(SICI)1097-0134
Unsupervised guided docking of
covalently bound ligands. Journal of [30] Liu M, Wang S. MCDOCK: A Monte
Computer-Aided Molecular Design. Carlo simulation approach to the
41
Drug Discovery and Development - New Advances
[31] Trosset JY, Scheraga HA. [38] Yang JM, Kao CY. Flexible ligand
PRODOCK: Software package for docking using a robust evolutionary
protein modeling and docking. Journal algorithm. Journal of Computational
of Computational Chemistry. 1999; Chemistry. 2000;21(11):988-998. DOI:
20(4):412-427. DOI: 10.1002/(SICI) 10.1002/1096-987X(200008)21:11<
1096-987X(199903)20:4<412::AID- 988::AID-JCC8>[Link];2-H
JCC3>[Link];2-N
[39] Clark KP. Flexible ligand docking
[32] McMartin C, Bohacek RS. QXP: without parameter adjustment across
Powerful, rapid computer algorithms four ligand-receptor complexes. Journal
for structure-based drug design. Journal of Computational Chemistry. 1995;
of Computer-Aided Molecular Design. 16(10):1210-1226. DOI: 10.1002/
1997;11(4):333-344. DOI: 10.1023/A: jcc.540161004
1007907728892
[40] Oshiro CM, Kuntz ID, Dixon JS.
[33] Meiler J, Baker D. Flexible ligand docking using a genetic
ROSETTALIGAND: Protein-small algorithm. Journal of Computer-Aided
molecule docking with full side-chain Molecular Design. 1995;9(2):113-130.
flexibility. Proteins: Structure, Function, DOI: 10.1007/BF00124402
and Bioinformatics. 2006;65(3):
538-548. DOI: 10.1002/prot.21086 [41] Jones G, Willett P, Glen RC.
Molecular recognition of receptor sites
[34] Morley SD, Afshar M. Validation of using a genetic algorithm with a
an empirical RNA-ligand scoring description of desolvation. Journal of
function for fast flexible docking using Molecular Biology. 1995;245(1):43-53.
RiboDock®. Journal of Computer-Aided DOI: 10.1016/S0022-2836(95)80037-9
Molecular Design. 2004;18(3):189-208.
DOI: 10.1023/B:JCAM.0000035199. [42] Österberg F, Morris GM, Sanner
48747.1e MF, Olson AJ, Goodsell DS. Automated
docking to multiple target structures:
[35] Choi V. Yucca: An efficient Incorporation of protein mobility and
algorithm for small-molecule docking. structural water heterogeneity in
Chemistry & Biodiversity. 2005;2(11): AutoDock. Proteins: Structure,
1517-1524. DOI: 10.1002/cbdv. Function, and Bioinformatics. 2002;
200590123 46(1):34-40. DOI: 10.1002/prot.10028
[36] Goodsell DS, Olson AJ. Automated [43] Taylor JS, Burnett RM. DARWIN: A
docking of substrates to proteins by program for docking flexible molecules.
simulated annealing. Proteins: Proteins: Structure, Function, and
Structure, Function, and Bioinformatics. Bioinformatics. 2000;41(2):173-191.
1990;8(3):195-202. DOI: 10.1002/ DOI: 10.1002/1097-0134(20001101)41:
prot.340080302 2<173::AID-PROT30>[Link];2-3
[37] Morris GM, Goodsell DS, Halliday [44] Corbeil CR, Englebienne P,
RS, Huey R, Hart WE, Belew RK, et al. Moitessier N. Docking ligands into
Automated docking using a Lamarckian flexible and solvated macromolecules. 1.
genetic algorithm and an empirical Development and validation of FITTED
binding free energy function. Journal of 1.0. Journal of Chemical Information
42
Molecular Docking in Modern Drug Discovery: Principles and Recent Applications
DOI: [Link]
43
Drug Discovery and Development - New Advances
44
Molecular Docking in Modern Drug Discovery: Principles and Recent Applications
DOI: [Link]
[73] Feher M, Deretey E, Roy S. BHB: A [81] Yang C-Y, Wang R, Wang S. M-
simple knowledge-based scoring score: A knowledge-based potential
function to improve the efficiency of scoring function accounting for protein
database screening. Journal of Chemical atom mobility. Journal of Medicinal
Information and Computer Sciences. Chemistry. 2006;49(20):5903-5911.
2003;43(4):1316-1327. DOI: 10.1021/ DOI: 10.1021/jm050043w
ci030006i
[82] Wójcikowski M, Ballester PJ,
[74] Mitchell JB, Laskowski RA, Alex A, Siedlecki P. Performance of machine-
Thornton JM. BLEEP-potential of mean learning scoring functions in structure-
force describing protein-ligand based virtual screening. Scientific
interactions: I. Generating potential. Reports. 2017;7:46710. DOI: 10.1038/
Journal of Computational Chemistry. srep46710
1999;20(11):1165-1176. DOI: 10.1002/
(SICI)1096-987X(199908)20:11 < [83] Su M, Yang Q, Du Y, Feng G, Liu Z,
1165::AID-JCC7>[Link];2-A Li Y, et al. Comparative assessment of
scoring functions: The CASF-2016
[75] Huang S-Y, Zou X. Inclusion of
update. Journal of Chemical
solvation and entropy in the knowledge-
Information and Modeling. 25 Feb 2019;
based scoring function for protein-
59(2):895-913. DOI: 10.1021/acs.
ligand interactions. Journal of Chemical
jcim.8b00545. Epub 2018 Dec 11
Information and Modeling. 2010;50(2):
262-273. DOI: 10.1021/ci9002987
[84] Liu J, Wang R. Classification of
[76] Kumar A, Goyal R, Jain S. Docking current scoring functions. Journal of
Methodologies and Recent Advances, Chemical Information and Modeling.
Oncology: Breakthroughs in Research and 2015;55(3):475-482. DOI: 10.1021/
Practice. Hershey, PA: IGI Global; 2017. ci500731a
pp. 804-828. DOI: 10.4018/978-1-
5225-0549-5.ch031 [85] Schames JR, Henchman RH, Siegel
JS, Sotriffer CA, Ni H, McCammon JA.
[77] Brás N, Cerqueira N, Sousa S, Discovery of a novel binding trench in
Fernandes P, Ramos M. Protein ligand HIV integrase. Journal of Medicinal
docking docking in drug discovery drug Chemistry. 2004;47(8):1879-1881. DOI:
discovery. In: Protein Modelling. 10.1021/jm0341913
London: Springer; 2014. pp. 249-286.
DOI: 10.1007/978-3-319-09976-7_11 [86] Evers A, Klabunde T. Structure-
based drug discovery using GPCR
[78] Náray-Szabó G. Protein Modelling. homology modeling: Successful virtual
London: Springer; 2014. DOI: 10.1007/ screening for antagonists of the alpha1A
978-3-319-09976-7 adrenergic receptor. Journal of
Medicinal Chemistry. 2005;48(4):
[79] Gohlke H, Hendlich M, Klebe G. 1088-1097. DOI: 10.1021/jm0491804
Knowledge-based scoring function to
predict protein-ligand interactions1. [87] Singh J, Chuaqui CE, Boriack-Sjodin
Journal of Molecular Biology. 2000; PA, Lee W-C, Pontz T, Corbley MJ,
295(2):337-356. DOI: 10.1006/ et al. Successful shape-based virtual
jmbi.1999.3371 screening: The discovery of a potent
inhibitor of the type I TGFβ receptor
[80] Muegge I. PMF scoring revisited. kinase (TβRI). Bioorganic & Medicinal
Journal of Medicinal Chemistry. 2006; Chemistry Letters. 2003;13(24):
45
Drug Discovery and Development - New Advances
46
Molecular Docking in Modern Drug Discovery: Principles and Recent Applications
DOI: [Link]
47
Chapter 4
Computational Deorphaning of
Mycobacterium tuberculosis Targets
Lorraine Yamurai Bishi, Sundeep Chaitanya Vedithi,
Tom L. Blundell and Grace Chitima Mugumbate
Abstract
1. Introduction
49
Drug Discovery and Development - New Advances
2. Method
Cole and co-workers [14] in 1998 reported the complete sequence of Mtb, which
comprises of 4,411,529 base pairs. The genome has an evenly distributed guanine-
cysteine content of 65.6% and represents the second-largest bacterial genome
sequence currently available. Additionally, the genome is rich in repetitive DNA,
particularly insertion sequences, and in new multi-gene families and duplicated
housekeeping genes, providing evidence for horizontally-transferred pathogenicity
islands of a particular base composition [14].
The genome of Mtb has some exceptional features, for example there are over
200 genes that encode enzymes for the metabolism of fatty acids, comprising 6% of
the total (Table 1). Among these, about 100 are predicted to function in the oxida-
tion of fatty acids. This large number of Mtb enzymes that putatively have fatty
acids as substrates may be linked to the ability of this pathogen to grow in the tissues
of the infected host, where fatty acids maybe the major carbon source. Another
50
Computational Deorphaning of Mycobacterium tuberculosis Targets
DOI: [Link]
Table 1.
General classification of Mtb genes. Adopted from [15].
unusual feature of the Mtb genome is the presence of the unrelated Pro-Glu (PE)
and Pro-Pro-Glu (PPE) families of proteins that have conserved N-terminal
domains of 100 and 180 amino acids respectively. The antigenicity of these proteins
has led to the assumption that at least some of these proteins may be involved in
antigenic variation of Mtb during infection [15].
51
Drug Discovery and Development - New Advances
their TB treatment before beginning HIV treatment [20]. Patients who contract
MDR-TB with HIV have a very poor prognosis due to the duration of treatment;
these individuals frequently succumb within a few months. Therefore, there is an
urgent need to develop continually new active agents to combat MDR-TB which has
been compounded by the emergence of XDR-TB. Furthermore, cases of TDR-TB
have been noted in China, India, Africa, and Eastern Europe. In TDR-TB, the
Mycobacterium are resistant to all available therapeutics [19]. To address this, in
2012 the U.S. Food and Drug Agency (FDA) approved bedaquiline for MDR-TB [21]
and later delamanid was approved as a compassionate care option for XDR-TB and
TDR-TB infections, nonetheless the EMA approved both agents for MDR-TB [22].
The biggest challenge is that these drugs have reported human ether-a-go-go related
gene (hERG) toxicity, as well as multiple absorption, distribution, metabolism and
excretion (ADME) issues due to their high lipophilicity [21]. This leads to an urgent
need for development of new agents that have successful therapeutic effects.
Figure 1.
Circular diagram of the genome of Mtb genes, essential proteins and the number of proteins that are drug
targets.
52
Computational Deorphaning of Mycobacterium tuberculosis Targets
DOI: [Link]
Table 2.
Mtb drug targets and the current used drugs [24–26].
53
Drug Discovery and Development - New Advances
54
Computational Deorphaning of Mycobacterium tuberculosis Targets
DOI: [Link]
55
Drug Discovery and Development - New Advances
function. In addition, the conformations are ranked based on their calculated bind-
ing energy. Highly ranked compounds are selected as ligands for the target. On the
other hand, reverse or inverse docking is used for identifying targets of drug pheno-
typic hits from a sea of targets. In this way, structure-based screening helps to iden-
tify and explain polypharmacology, molecular mechanism of action of substances,
facilitate drug repurposing, detect adverse drug reactions and hence toxicity.
Figure 2.
(A) Binding modes of two fragment-like molecules inside the long cylindrical allosteric binding pocket of
EthR defined by five helices. Yellow sticks depict the molecule occupying the upper binding site close to the
entrance of the pocket and cyan sticks represent a molecule occupying the inner binding site close to the HTH
domain. (B) EthR-ligand interactions involving Trp103 (yellow) at the entrance of the binding pocket of the
protein. Ligand atoms and bonds are in pink, grey rings are hydrophobic interactions, red rings show hydrogen
bonds. (C) EthR-ligand (pink) interactions involving Phe110 located at the center of the binding pocket of
EthR.
56
Computational Deorphaning of Mycobacterium tuberculosis Targets
DOI: [Link]
The Structural Genomics Consortium (SGC) [39] which has focused on proteins
of interest to medicine, has impressive achievements, in 2011 defining ~40% of
the structures of proteins from human parasites deposited in the PDB [40]. The
Tuberculosis Structural Genomics Consortium (TBSGC), an international col-
laboration involving 53 countries, has focused on 3D structures of Mtb [40]. This
activity and others working on Mtb proteomes have deposited 2274 structures in
the PDB, but still representing less than 583 gene products, only 13.97% of genome.
Although this is a small percentage, it compares impressively with knowledge
of protein structures of two other mycobacterial pathogens where there is great
clinical interest: for M. leprae causing leprosy there are experimentally-defined 3D
structures for 15 gene products and for M. abscessus, a free living Mycobacterium,
which is a growing challenge for cystic fibrosis patients, there are 53 experimen-
tally-defined 3D structures in the PDB.
57
Drug Discovery and Development - New Advances
Rapid progress in this and other related software coupled with increasing com-
puting power has enabled genome scale prediction of protein structures, as a viable
alternative to experimental determination. In order to construct computational models
of all gene products, which we here refer to as the structural proteome, we identify
templates by a sequence-structure homology search using Fugue [47], which uses
local-structural-environment-specific substitution tables to predict the likelihood of a
common 3D structure. We have incorporated Fugue into a pipeline (Vivace), in which
templates are selected from TOCCATA (Ochoa Montaño and Blundell, unpublished),
a database of consensus profiles built from CATH 3.5 [48] and SCOP 1.75A [49] based
classification of proteins structures (PDB files). PDBs within each profile are clustered
based on sequence similarity using CD-HIT [50] and structures are aligned using
BATON, a modified version of COMPARER [45]. After further optimization of the
clusters by discarding templates with more than 20% difference in sequence identity to
the maximum hit, remaining templates are classified into states based on ligand binding
and oligomerization. Five different states, known as “liganded-monomeric,” “liganded-
complexed,” “apo-monomeric,” “apo-complexed” and “any,” are generated in each
profile hit. Models are built in each of these states using Modeller 9.10 [46] and refined.
Later NDOPE, GA341 [51] Molprobity [52] and SSAG [53] are used to determine the
quality of the models.
The first application of this approach was to construct the Chopin Database
([Link] a database of protein structures for
H37Rv strain of Mtb. This has provided structures that are reasonably certain for around
65% of gene products. These have proved reliable indicators of the overall structures but
may have some uncertainties especially in loop regions and domain-domain relation-
ships. A further ~19% probably have correct folds while the remaining would unlikely
to be correct. Nevertheless, compared to those structures defined experimentally by
X-ray analysis, this represents a 6-fold increase of structural information available that
might be useful in assessing druggability and the impacts of mutations.
Similar models of the structural proteome for M. abscessus (Skwark et al.,
unpublished) and M. leprae (Vedithi et al., unpublished) have been developed in
the group. In M. leprae, of the 1615 gene products, templates were identified for
1429 gene products and we were able to model 1161 proteins with high confidence.
A total of 36,408 models were built in different ligand bound and oligomeric states
for the 1161 proteins. The distribution of Fugue Z score across models indicates
that only 4% of the proteome has no hits and 15% has poor scores. ~80% of the
proteome has acceptable and good hits, and the corresponding Z scores. Around
47% of the protein queries identified templates with identity and coverage greater
than 40 and 67% of the models in the proteome are of best quality as estimated by
NDOPE, GA341, Molprobity and Secondary Structure Agreement (SSAG).
58
Computational Deorphaning of Mycobacterium tuberculosis Targets
DOI: [Link]
We have also spent time over 2 decades analysing the impacts of mutations
evident in the increasing wealth of available genome sequences for pathogenic myco-
bacteria and cancers. We originally developed SDM [56] in 1997, a method depend-
ing on statistical analysis of environment-dependent amino-acid substitution tables
[57, 58]. In 2013 machine learning was introduced with the arrival of Douglas Pires in
Cambridge, developing first mCSM for stability [59] followed by several “flavours”
including mCSM-PPI for impacts on protein-protein interactions, mCSM-NA [60]
for nucleic acid interactions and mCSM-lig for impacts on small-molecule ligand
interactions useful for understanding drug resistance [61]. A critical part of using
machine learning is to have an extensive database of experimentally-defined impacts
of mutations on stability and interactions, such as Platinum by David Ascher when
in Cambridge [62], a database of experimentally measured effects of mutations on
structurally defined protein-ligand complexes that was developed for mCSM-lig.
These two structural approaches to predicting the impacts of mutations (SDM &
mCSM) have proved complementary and more reliable than most sequence-only
Figure 3.
(A) Protomeric model of PyrG (CTP-Synthase) of M. leprae modelled with a quality of 4.25 (best).
(B) Homo-8-mer of PyrG of M. leprae modelled with a quality of 4.25 (best).
59
Drug Discovery and Development - New Advances
Figure 4.
Indicates the maximum destabilising effect a mutation can induce on the stability of RNA-polymerase
β-subunit of M. leprae (target for rifampin) measured by mCSM-stability.
60
Computational Deorphaning of Mycobacterium tuberculosis Targets
DOI: [Link]
indicating donor, acceptor and lipophilic interactions correlate well with experimental
binding sites of fragments that can be elaborated in fragment-based discovery. For the
ligand bound structures, lower contouring can provide “warm spots” for the binding
sites, indicating possibilities for elaborating the fragment in the binding pocket.
The models of individual molecules of the modelled proteome can be individually
decorated with the hotspot maps. They give a good indication of the known func-
tional sites on experimentally defined structures of proteins, often demonstrating
that a functional site comprises several hotspots involved in binding substrates and
cofactors. They also provide a good indication of the location of allosteric sites [70].
7. Conclusion
Acknowledgements
LYB and GCM are grateful to Chinhoyi University of Technology for their
support in introducing computational drug discovery and development research
work at the University and all our collaborators. TLB and SCV thank the Gates
Foundation, the Cystic Fibrosis Trust and the American Leprosy Mission for their
funding of computational and experimental work on approaches to combating
disease from mycobacterial infections. They also thank colleagues in Cambridge
and elsewhere who have contributed over the years to our efforts to develop new
approaches to structural biology, computational bioinformatics and drug discovery.
61
Drug Discovery and Development - New Advances
Author details
© 2019 The Author(s). Licensee IntechOpen. This chapter is distributed under the terms
of the Creative Commons Attribution License ([Link]
by/3.0), which permits unrestricted use, distribution, and reproduction in any medium,
provided the original work is properly cited.
62
Computational Deorphaning of Mycobacterium tuberculosis Targets
DOI: [Link]
References
63
Drug Discovery and Development - New Advances
[16] Hoagland DT, Liu J, Lee RB, Lee [25] Kaneko T, Cooper C, Mdluli K.
RE. New agents for the treatment Challenges and opportunities in
of drug-resistant Mycobacterium developing novel drugs for TB. Future
tuberculosis ☆. Advanced Drug Delivery Medicinal Chemistry. 2011;3:1373-1400.
Reviews. 2016;102:55-72. DOI: 10.1016/j. DOI: 10.4155/fmc.11.115
addr.2016.04.026
[26] Campaniço A, Moreira R, Lopes F.
[17] World Health Organisation. Global Drug discovery in tuberculosis. New
tuberculosis Report 2017. Geneva; 2017 drug targets and antimycobacterial
agents. European Journal of Medicinal
[18] Zhang Y, Post-Martens K, Denkin S. Chemistry. 2018;150:525-545. DOI:
New drug candidates and therapeutic 10.1016/[Link].2018.03.020
targets for tuberculosis therapy. Drug
Discovery Today. 2006;11:21-27. DOI: [27] Chiliza TE, Pillay M, Pillay B.
10.1016/S1359-6446(05)03626-3 Identification of unique essential proteins
from a Mycobacterium tuberculosis F15/
[19] Pawlowski A, Jansson M, Sköld M, LAM4/KZN phage secretome library.
Rottenberg ME, Källenius G. Tuberculosis Pathogens and Disease. 2017;75:1-10. DOI:
and HIV co-infection. PLoS Pathogens. 10.1093/femspd/ftx001
2012;8(2):e1002464. [Link]
org/10.1371/[Link].1002464 [28] Kapetanovic IM. Computer-aided
drug discovery and development
[20] Metcalfe JZ, Porco TC, Westenhouse (CADDD): In silico-chemico-biological
J, Damesyn M, Facer M, Hill J, et al. approach. Chemico-Biological
Tuberculosis and HIV co-infection, Interactions. 2008;171:165-176. DOI:
California, USA, 1993-2008. Emerging 10.1016/[Link].2006.12.006
Infectious Diseases. 2013;19:400-406.
DOI: 10.3201/eid1903.121521 [29] Machado E, Cerdeira C, de Miranda
AB, Catanho M. Web resources on
[21] Worley MV, Estrada SJ. Bedaquiline: tuberculosis: Information, research,
A novel antitubercular agent for the and data analysis. In: Mycobacterium-
treatment of multidrug-resistant research and development. IntechOpen
tuberculosis. The Journal of Human
Pharmacology and Drug Therapy. [30] Gaulton A, Bellis LJ, Bento AP,
2014;34(11):1187-1197 Chambers J, Davies M, Hersey A, et al.
{ChEMBL}: A large-scale bioactivity
[22] Gawad J, Bonde C. Current affairs, database for drug discovery. Nucleic
future perspectives of tuberculosis Acids Research. 2012;40:D1100-D1107.
and antitubercular agents. The Indian DOI: 10.1093/nar/gkr777
Journal of Tuberculosis. 2018;65:15-22.
DOI: 10.1016/[Link].2017.08.011 [31] Malhotra S, Mugumbate G,
Blundell TL, Higueruelo AP. TIBLE: A
[23] Lamichhane G. Novel targets web-based, freely accessible resource
in M. tuberculosis: Search for new for small-molecule binding data for
drugs. Trends in Molecular Medicine. mycobacterial species. Database
2011;17:25-33. DOI: 10.1016/j. (Oxford). 2017;2017:1-7. DOI: 10.1093/
molmed.2010.10.004 database/bax041
[24] Lewis K. Platforms for antibiotic [32] Mugumbate G, Newton AS, Rosenthal
discovery. Nature Reviews. Drug PJ, Gut J, Moreira R, Chibale K, et al.
64
Computational Deorphaning of Mycobacterium tuberculosis Targets
DOI: [Link]
65
Drug Discovery and Development - New Advances
[46] Šali A, Blundell TL. Comparative [54] Dejesus MA, Gerrick ER, Xu W,
protein modelling by satisfaction of Park SW, Long JE, Boutte CC, et al.
spatial restraints. Journal of Molecular Comprehensive essentiality analysis
Biology. 1993;234:779-815. DOI: of the Mycobacterium tuberculosis
10.1006/jmbi.1993.1626 genome via saturating transposon
mutagenesis. MBio Journal.
[47] Shi J, Blundell TL, Mizuguchi K. 2017;8:e02133-e02116. DOI: 10.1128/
FUGUE: Sequence-structure homology mBio.02133-16
recognition using environment-specific
substitution tables and structure- [55] Mori G, Chiarelli LR, Esposito M,
dependent gap penalties. Journal of Makarov V, Bellinzoni M, Hartkoorn
Molecular Biology. 2001;310:243-257. RC, et al. Thiophenecarboxamide
DOI: 10.1006/jmbi.2001.4762 derivatives activated by EthA kill
Mycobacterium tuberculosis by inhibiting
[48] Orengo C, Michie A, Jones S, Jones the CTP synthetase PyrG. Chemistry &
D, Swindells M, Thornton J. Biology. 2015;22:917-927. DOI: 10.1016/j.
CATH—A hierarchic classification of chembiol.2015.05.016
protein domain structures. Structure.
1997;5:1093-1109. DOI: 10.1016/ [56] Topham CM, Srinivasan N, Blundell
S0969-2126(97)00260-8 TL. Prediction of the stability of
protein mutants based on structural
[49] Rost B, Brenner SE, Chothia C, environment-dependent amino acid
Hubbard TJP, Murzin AG, Li WW, et al. substitution and propensity tables.
CKAAPs DB: A conserved key amino Protein Engineering. 1997;10:7-21. DOI:
acid positions database. Methods in 10.1093/protein/10.1.7
Enzymology. 2002;28:409-411. DOI:
10.1016/S0076-6879(96)66039-X [57] Worth CL, Preissner R, Blundell TL.
SDM—A server for predicting effects
[50] Fu L, Niu B, Zhu Z, Wu S, Li W. of mutations on protein stability and
CD-HIT: Accelerated for clustering malfunction. Nucleic Acids Research.
the next-generation sequencing data. 2011;39:W215-W222. DOI: 10.1093/nar/
Bioinformatics. 2012;28:3150-3152. DOI: gkr363
10.1093/bioinformatics/bts565
[58] Pandurangan AP, Ochoa-Montaño B,
[51] Melo F, Sali A. Fold assessment for Ascher DB, Blundell TL. SDM: A server
comparative protein structure modeling. for predicting effects of mutations
Protein Science. 2007;16:2412-2426. on protein stability. Nucleic Acids
DOI: 10.1110/ps.072895107 Research. 2017;45:W229-W235. DOI:
10.1093/nar/gkx439
[52] Davis IW, Murray LW, Richardson
JS, Richardson DC. MolProbity: [59] Pires DEV, Ascher DB, Blundell
Structure validation and all-atom TL. MCSM: Predicting the effects of
contact analysis for nucleic acids mutations in proteins using graph-based
and their complexes. Nucleic Acids signatures. Bioinformatics. 2014;30:
Research. 2004;32:W615-W619. DOI: 335-342. DOI: 10.1093/bioinformatics/
10.1093/nar/gkh398 btt691
[53] Eramian D, Shen M, Devos D, [60] Pires DEV, Ascher DB. MCSM-NA:
Melo F, Sali A, Marti-Renom MA. A Predicting the effects of mutations
composite score for predicting errors on protein-nucleic acids interactions.
in protein structure models. Protein Nucleic Acids Research. 2017;45:
Science. 2006;15:1653-1666. DOI: W241-W246. DOI: 10.1093/nar/
10.1110/ps.062095806 gkx236
66
Computational Deorphaning of Mycobacterium tuberculosis Targets
DOI: [Link]
[61] Pires DEV, Blundell TL, Ascher [68] Verdonk ML, Cole JC, Taylor R.
DB. MCSM-lig: Quantifying the SuperStar: A knowledge-based approach
effects of mutations on protein-small for identifying interaction sites in
molecule affinity in genetic disease and proteins. Journal of Molecular Biology.
emergence of drug resistance. Scientific 1999;289:1093-1108. DOI: 10.1006/
Reports. 2016;6:29575. DOI: 10.1038/ jmbi.1999.2809
srep29575
[69] Radoux CJ, Olsson TSG, Pitt WR,
[62] Pires DEV, Blundell TL, Ascher Groom CR, Blundell TL. Identifying
DB. Platinum: A database of interactions that determine fragment
experimentally measured effects of binding at protein hotspots.
mutations on structurally defined Journal of Medicinal Chemistry.
protein-ligand complexes. Nucleic Acids 2016;59:4314-4325. DOI: 10.1021/acs.
Research. 2015;43:D387-D391. DOI: jmedchem.5b01980
10.1093/nar/gku966
[70] Thomas SE, Mendes V, Kim
[63] Pires DEV, Chen J, Blundell TL, SY, Malhotra S, Ochoa-Montaño B,
Ascher DB. In silico functional dissection Blaszczyk M, et al. Structural biology
of saturation mutagenesis: Interpreting and the design of new therapeutics:
the relationship between phenotypes and From HIV and cancer to mycobacterial
changes in protein stability, interactions infections: A paper dedicated to John
and activity. Scientific Reports. Kendrew. Journal of Molecular Biology.
2016;6:19848. DOI: 10.1038/srep19848 2017;429:2677-2693. DOI: 10.1016/j.
jmb.2017.06.014
[64] Desaphy J, Azdimousa K,
Kellenberger E, Rognan D. Comparison [71] Jones G, Willett P, Glen RC, Leach
and druggability prediction of AR, Taylor R. Development and
protein-ligand binding sites from validation of a genetic algorithm for
pharmacophore-annotated cavity flexible docking. Journal of Molecular
shapes. Journal of Chemical Information Biology. 1997;267:727-748. DOI:
and Modeling. 2012;52:2287-2299. DOI: 10.1006/jmbi.1996.0897
10.1021/ci300184x
67
Section 3
69
Chapter 5
Abstract
1. Introduction
71
Drug Discovery and Development - New Advances
2. Pharmacokinetic study
72
Revisiting Pharmacokinetics and Pharmacogenetics of Methadone in Healthy Volunteers
DOI: [Link]
arterial blood to tissues such as MTD, such increased tissue/venous plasma ratio
would explain the toxicity of MTD in certain tissues and the lack of correlation
between venous MTD concentrations and adverse effects [15]. When elimination
predominates, the opposite is observed [30]. However, if MTD recycling is operat-
ing at the monoexponential decay of levels, increased arterial/venous plasma drug
concentration ratios will also be observed due to drug reabsorptions.
So its storage in body tissues and the slow release to plasma as well as its
recycling process could be responsible for its prolonged elimination half-life.
This last fact is exploited in preventing withdrawal symptoms. However, the long
half-life does not seem to correlate with the observed shorter duration of analgesia
(6–12 hours) after steady state is reached [31].
As measuring drug levels in arteries is an uncommon practice, our group has
been working for a long time [32–35] using saliva in order to surrogate arterial free
plasma drug concentrations as this biological fluid highly correlates with arte-
rial plasma due to the fact that it is produced by ultrafiltration of the latter [36].
Salivary peaks during the elimination phase would be indicating reentry processes
as it was observed in a study carried out with patients [35].
It is important to study the stereoselectivity of MTD metabolism once the blood-
gastrointestinal tract cycling is operating, and to investigate whether the intestinal
metabolism of MTD could be assessed as relevant in relation with the hepatic one.
For this purpose, our research group has carried out an in vivo study.
Blood samples were withdrawn from the antecubital vein through cannula-
tion and saliva samples were collected in Salivette® tubes at the following times:
0–0.5–1–2–3–4–6–8–10–12–16–24–36–48–72 and 96 hours post dose. Urine was col-
lected at 0 (before dose intake) and at the end of the following intervals: 0–2, 2–4,
4–7, 7–8.5, 8.5–10, 10–11.5, 11.5–13, 13–14.5, 14.5–16, and 16–24 hours after dosing
and sample volumes were recorded. Aliquots of urine samples were kept in order to
73
Drug Discovery and Development - New Advances
measure the analyte content. Immediately after sampling, pH was measured using
a portable pH meter for urine samples. All samples were kept in a freezer at −25°C
until the time of analysis.
When the pre-dose blood sample was taken, another blood sample was taken
to obtain genomic DNA in order to determine the genotype of the CYP2B6 and
CY2C19 isoenzymes of the subjects.
MTD enantiomers in plasma, saliva, and urine were quantified. EDDP enantiomer
quantification was performed in urine. MTD and EDDP were extracted with a mixture
of hexane and isoamyl alcohol from 2.0 mL of plasma or 1.0 mL of urine or saliva
samples that were previously alkalinized. Then, the organic phase was evaporated
under a stream of nitrogen, and the residue was reconstituted with the mobile phase.
Imipramine (10.00 μg/mL) was used as the internal standard and 50 mL was added
to plasma or urine or saliva. MTD (in all the three fluids) and EDDP (only in urine)
quantification was performed using a validated HPLC-UV chiral method, which was
an adaptation of a previously published methodology [38]. The mobile phases con-
sisted of phosphate buffer 20 mM pH 6.0 + 2 mM diisopropylamine: acetonitrile (92:8)
for urine analysis and phosphate buffer 20 mM pH 7.0 + 2 mM diisopropylamine:
acetonitrile (82:18) for plasma and saliva analysis. The flow rate was 0.7 mL/min. The
separation of the compounds was performed on a CHIRALPACK AGP™ (100 × 4 mm;
5 μm) column with a silica guard column. Detection was performed at a wavelength of
215 nm. The analysis was carried out at 25°C and the injection volume was 80 μL.
The HPLC method was linear for MTD between 4.0 and 160 ng/mL and
between 19.0 and 3280 ng/mL for plasma or saliva and urine samples, respectively.
The linearity for EDDP in urine was proven from 52.0 to 4200 ng/mL. Inter- and
intra-day precision and accuracy were below 14% for both compounds.
The following pharmacokinetic parameters were obtained from the MTD plasma
and saliva concentration versus time curves for both enantiomers of MTD:
Experimental Cmax and Tmax were computed and the AUC was estimated by
the trapezoid method up to 96 hours, or until the last quantifiable concentration
time. As for most of the subjects, the concentrations were not quantifiable for times
longer than 24 h and AUC was determined up to 24 h. The R/S concentration ratio
was computed as an indicator of possible stereoselective metabolic changes because
of drug recycling.
From the urinary concentrations of MTD and EDDP and the volumes of urine
recorded, the amounts excreted in the time interval between two consecutive
micturitions were calculated. Excretion rates versus time were plotted and the R/S
ratios of MTD and EDDP were calculated for this parameter.
Statistical significances between means were assessed by a nonpaired (between
sexes) and a paired (between enantiomers) t-student test.
74
Revisiting Pharmacokinetics and Pharmacogenetics of Methadone in Healthy Volunteers
DOI: [Link]
Figure 1.
Mean R- and S-MTD plasma concentration-time profiles for treatments A and B in women and men.
R S R S R S R S
b b b b
Plasma Women 27.3 ± 12.6 22.1 ± 7.8 190 ± 61 302 ± 89 18.6 ± 10.3 32.2 ± 8.4 3.5 2.0
(2.0–8.0) (1.0–8.0)
Men 25.0 ± 4.7 24.6 ± 6.8 192b ± 101 304b ± 118 17.5b ± 7.9 31.0b ± 6.3 3.0 2.0
(1.0–10) (1.0–8.0)
a
median (range).
b
p < 0.01, paired t-student test between R and S.
Table 1.
Mean (± standard deviation) pharmacokinetic parameters of MTD obtained in women and in men.
75
Drug Discovery and Development - New Advances
Figure 2.
Mean urinary excretion rates of R- and S-MTD versus time in men and women.
Figure 3.
Mean urinary excretion rates of R- and S-EDDP versus time in men and women.
76
Revisiting Pharmacokinetics and Pharmacogenetics of Methadone in Healthy Volunteers
DOI: [Link]
Figure 4.
Mean R- and S-MTD saliva concentration-time profiles for treatments A and B in women and men.
Figure 5.
Mean (± 95%CI) R-to-S MTD and EDDP urinary excretion rate ratios after oral administration of MTD.
had finished. During the absorption and rapid disposition phase, this ratio is
increasing. However, R/S ratios of EDDP were constant from the beginning, except
after food intake (mainly between 3 and 7 hours post lunch intake when MTD recir-
culation is taking place) when the ratio decreased and this might explain differences
in EDDP systemic formation.
The molecules of MTD present in the systemic circulation undergo both intestinal
and hepatic stereoselective metabolism by CYP2C19 and CYP2B6 enzymes. CYP2C19
is stereoselective towards the R-isomer while CYP2B6 towards the S-isomer. After food
intake, when a process of drug reentry is operating, the molecules of MTD that had
been secreted into the gastric juice can be reabsorbed in the intestine. Consequently,
a greater number of molecules enter the enterocyte. The change observed in the R/S
ratio of EDDP after the ingestion of meals evidences a different stereoselectivity
between intestinal and hepatic metabolism, possibly due to a relative differential
content of CYP3A4 and CYP2C19 in enterocytes and hepatocytes, being the relative
presence of CYP3A4 greater at the intestine. In the case of MTD, the metabolism of the
S-enantiomer is favored after the passage of MTD through the intestine compared to
its passage through the liver. Although during food intake there is an increased blood
flow to the splanchnic area, and the liver and the other organs in this region receive a
greater number of molecules from the blood coming from areas that do not belong to
the splanchnic region, for drugs secreted in the gastric juice, the fraction of molecules
that the intestine receives is even greater because there is a supplementary quantity of
molecules that enter the intestine coming from the gastric juice. If no secretion was
taking place, the molecules would be transferred from the stomach directly to the liver
through the portal bloodstream without passing through the enterocytes.
Therefore, by favoring recirculation rather than bypassing hepatic metabolism,
the intestinal metabolism would be increasing. Our research reveals an important
role of the intestine in the systemic (and pre-systemic) metabolism of MTD,
presenting a greater stereoselectivity towards the S-isomer. Although this isomer
has little or no activity as an opioid agonist, it is able to inhibit the reuptake of
77
Drug Discovery and Development - New Advances
3. Pharmacogenetic study
CYP2B6 and not CYP3A4 is the principle determinant of clinical MTD elimina-
tion and is one of the most polymorphic cytochrome P450 (P450) genes in humans
and, currently, it has 30 defined alleles with over 100 described polymorphisms
[39]. According to Kharasch et al. [40], CYP2B6*6 allele carriers showed higher
MTD concentrations and slower elimination, whereas CYP2B6*4 carriers had lower
concentrations and faster elimination.
CYP2C19 plays an important role in MTD metabolism and CYP2C19 gene is highly
polymorphic as well. Loss of enzyme activity results from the CYP2C19*2 allele and
the CYP2C19*17 allele is associated with increased enzymatic activity [41, 42].
3.1 Methodology
Once the genomic DNA was obtained from the leukocyte fraction, the individu-
als were genotyped for the CYP2B6 and CYP2C19 genes by massive sequencing,
which was carried out at the Institute of Genomic Medicine (INMEGEN) in Mexico.
In order to be processed by massive sequencing, genomic DNA samples should
have a concentration higher than 10 ng/μL, and the ratio of absorbances 260/280
and 260/230 should be approximately 2 to be able to consider that the DNA
obtained was of good quality. In cases in which the sample did not meet these
requirements, purification was performed using the Mag Jet Genomic DNA Kit
(Thermo Scientific) which includes incubation with proteinase and RNAse and
purification with magnetic beads.
As a result of this processing, the genotype of the 12 volunteers was obtained for
CYP2B6 and CYP2C19 enzymes. Considering the polymorphisms found and based
on the literature, we determined the phenotype that would be expected, that is,
increased, normal, or decreased enzyme activity.
Regarding the polymorphisms in the gene that encodes CYP2C19, 5 of the volun-
teers in our study presented the allelic variant * 2 (rs4244285), which is associated with
a decrease in the activity of the enzyme, whereas 2 volunteers presented the allelic
variant * 17 (rs3758581), which is associated with an increase in the activity. Regarding
the polymorphisms in the gene that encodes CYP2B6, 6 volunteers presented the allelic
variant * 4 (rs2279343), which determines an increased enzymatic activity.
S/R ratios for MTD in plasma and urine and the S/R ratios for EDDP in urine
were calculated. The individuals were grouped into two. Group 1 included those
volunteers in whom the activity of CYP2B6 was increased and CYP2C19 activity
was normal or decreased as well as those volunteers in whom CYP2B6 activ-
ity was normal but CYP2C19 activity was diminished. Group 2 included those
individuals with normal activity of both enzymes as well as those in which the
activity of CYP2B6 was normal but that of CYP2C19 was increased and a case
78
Revisiting Pharmacokinetics and Pharmacogenetics of Methadone in Healthy Volunteers
DOI: [Link]
Table 2.
S/R ratios for MTD in plasma and urine and for EDDP in urine obtained in treatment A and the activity of
CYP2B6 and CYP2C19.
in which the activity of both enzymes was increased. This classification allowed
grouping those individuals, in whom a preferential biotransformation was
expected on the S isomer, considering the activity of the enzyme together with
its stereoselectivity. The averages of the S/R ratios for each group were calculated,
both for treatment A and for treatment B, and the results are shown in Tables 2
and 3, respectively.
The three average S/R ratios were compared by a t-student test, and no signifi-
cant differences were obtained in any of the cases. However, the S/R ratios of MTD
either in plasma or in urine are lower in Group 1 compared to Group 2, which is
in agreement with the stereoselectivity of CYP2B6 towards the S-MTD since the
metabolism of the S-isomer is greater compared to the R-isomer when the activity
of CYP2B6 is increased and the activity of CYP2C19 decreased. The results obtained
for the S/R ratios of EDDP are different, probably because the biotransformation of
MTD mediated by these enzymes also leads to the formation of other metabolites.
Moreover, there is a lack of information in the literature about the stereoselectivity
of EDDP clearance.
Genetic variation of CYP2C19 mainly affects MTD metabolism, and it has a
minor effect on the metabolite, maybe because it contributes very little to EDDP
formation (1/10 compared to CYP2B6 contribution).
79
Drug Discovery and Development - New Advances
Table 3.
S/R ratios for MTD in plasma and urine and for EDDP in urine obtained in treatment B and activity of
CYP2B6 and CYP2C19.
4. Conclusions
Our results confirm MTD recirculation via gastric secretion and subsequent
intestinal reabsorption. MTD is extensively metabolized in the liver but intestinal
metabolism is not only relevant but also stereoselective.
Although the opioid effect of MTD is mainly due to the R-isomer, the S-isomer
also has an analgesic action by inhibiting the reuptake of serotonin and nor-
adrenaline and by exhibiting a noncompetitive antagonism of the NMDA receptor.
The latter action is also responsible for preventing or attenuating tolerance and
withdrawal syndrome. Therefore, in those patients who have an increased activity
of the CYP2B6 enzyme or a normal activity of this enzyme in combination with a
decreased activity of CYP2C19, (situations that favor the S-isomer metabolism), the
analgesic effect could be diminished and the development of tolerance as well as the
withdrawal symptoms could be exacerbated.
Despite the fact that blood-gastrointestinal-blood recycling extends the resi-
dence of a drug in the body, in this case, the elimination of the S-isomer is increased
with each passage through the enterocyte. Consequently, the recycling process of
MTD would not be favoring an increased analgesic effect as it would be expected.
This is in agreement with the shorter duration of analgesia observed in the clinical
setting after steady state is reached.
80
Revisiting Pharmacokinetics and Pharmacogenetics of Methadone in Healthy Volunteers
DOI: [Link]
The occurrence of frequent adverse effects such as nausea was observed only
in women, even after receiving two doses of metoclopramide prior to the dose of
MTD. Although tolerance to nausea and vomits develop with chronic use, the physi-
cian should consider a lower starting dose of 5 mg/day for women. Apparently, an
initial dose of 10 mg/day for men could be appropriate.
Acknowledgements
The genetic part of this work was supported by Consejo Nacional de Ciencia
y Tecnología (Conacyt, México) Research [grant number 248513]. Natalia
Guevara received fundings to perform this work from the Comisión Sectorial
de Investigación Científica (CSIC), Programa de Iniciación a la Investigación,
Universidad de la República (Uruguay).
Conflict of interest
Author details
© 2018 The Author(s). Licensee IntechOpen. This chapter is distributed under the terms
of the Creative Commons Attribution License ([Link]
by/3.0), which permits unrestricted use, distribution, and reproduction in any medium,
provided the original work is properly cited.
81
Drug Discovery and Development - New Advances
References
[9] Haumann J, Geurts JW, van Kuijk [17] Crettol S, Digon P, Golay KP,
SM, Kremer B, Joosten EA, van den Brawand M, Eap CB. In vitro P
82
Revisiting Pharmacokinetics and Pharmacogenetics of Methadone in Healthy Volunteers
DOI: [Link]
[19] Routledge PA. The plasma protein [27] Vinson RK. Pharmacokinetics of
binding of basic drugs. British a new immediate-release methadone
Journal of Clinical Pharmacology. tablet formulation with decreased
1986;22:499-506 in vitro solubility. Clinical Drug
Investigation. 2012;32:487-495
[20] Belpaire FM, De Rick A, Dello
C, Fraeyman N, Bogaert MG. Alpha [28] Vázquez M, Fagiolino P.
1-acid glycoprotein and serum binding Pharmacotherapy of chronic pain.
of drugs in healthy and diseased dogs. In: Maldonado C, editor. Pain Relief.
Journal of Veterinary Pharmacology and Rijeka, Croatia: IntechOpen; 2017.
Therapeutics. 1987;10(1):43-48 DOI: 10.5772/66444. Available
from: [Link]
[21] Foster DJ, Somogyi AA, Dyer com/books/pain-relief-from-
KR, White JM, Bochner F. Steady- analgesics-to-alternative-therapies/
state pharmacokinetics of (R)- and pharmacotherapy-of-chronic-pain
(S)-methadone in methadone
maintenance patients. British [29] Lynn RK, Olsen GD, Leger RM,
Journal of Clinical Pharmacology. Gordon WP, Smith RG, Gerber N.
2000;50(5):427-440 The secretion of methadone and
its major metabolite in the gastric
[22] Inturrisi CE, Verebely K. Disposition juice of humans: Comparison with
of methadone in man after a single blood and salivary concentrations.
oral dose. Clinical Pharmacology and Drug Metabolism and Disposition.
Therapeutics. 1972;13:923-930 1976;4(5):504-509
[23] Totah RA, Sheffels P, Roberts T, [30] Lam G, Chiou WL. Determination
Whittington D, Thummel K, Kharasch of the steady-state volume of
ED. Role of CYP2B6 in stereoselective distribution using arterial and venous
human methadone metabolism. plasma data from constant infusion
Anesthesiology. 2008;108(3):363-374 studies with procainamide. The Journal
of Pharmacy and Pharmacology.
[24] Gerber JG, Rhodes RJ, Gal J. 1982;34:132-134
Stereoselective metabolism of
methadone N-demethylation by [31] Sunilkumar MM, Lockman
cytochrome P4502B6 and 2C19. K. Practical pharmacology of
Chirality. 2004;16(1):36-44 methadone: A long-acting opioid.
Indian Journal of Palliative Care.
[25] Foster DJ, Somogyi AA, Bochner 2018;24:S10-S14
F. Methadone N-demethylation in
human liver microsomes: Lack of [32] Maldonado C, Fagiolino P, Vázquez
stereoselectivity and involvement of M, et al. Therapeutic carbamazepine
83
Drug Discovery and Development - New Advances
84
Chapter 6
Abstract
The drug discovery and development pipeline have more and more relied on
in vitro testing and in silico predictions to reduce investments and optimize lead
compounds. A comprehensive set of in vitro assays is available to determine key
parameters of absorption, distribution, metabolism, and excretion, for example,
lipophilicity, solubility, and plasma stability. Such test systems aid the evaluation
of the pharmacological properties of a compound and serve as surrogates before
entering in vivo testing and clinical trials. Nowadays, computer-aided techniques
are employed not just in the discovery of new lead compounds but embedded as
part of the entire drug development process where the ADME profiling and big data
analyses add a new layer of complexity to those systems. Herein, we give a short
overview of the history of the drug development pipeline presenting state-of-the-
art ADME in vitro assays as established in academia and industry. We will further
introduce the underlying good practices and give an example of the compound
development pipeline. In the next step, recent advances at in silico techniques will
be highlighted with special emphasis on how pharmacogenomics and in silico PK
profiling can enhance drug monitoring and individualization of drug therapy.
1. Introduction
85
Drug Discovery and Development - New Advances
The current state of pharmaceutical DDD estimates that only up to ten com-
pounds out of thousand screened hits would result in optimized leads and enter
preclinical testing, with a chance of 9.6% to pass the clinical testing phase [1, 2].
Additionally, the drug approval process is estimated to last in average 15 years, with
major expenses in phases II and III of clinical trials, which highlights the drawback
a failure in (pre-) clinical testing causes [3–6], where the overall DDD cost for each
drug can reach as high as ~$2.56 billion preapproval rising to $2.87 billion including
postapproval investments [6–8]. From the initial small molecule screened as hit to
the optimized lead, a variety of in vitro tests are performed to guarantee efficacy
and safety, but also to find structure-activity relationships (SAR), which can then
be connected to specific physicochemical properties of the compound and further
aid in the lead optimization phase [8–10].
The drug development phase starts with preclinical testing followed by the clini-
cal stage comprising phase I–III human trials. Each of the phases aims to answer
a specific question. Initially, preclinical trials are conducted in animals and can
provide information about whether a drug is toxic or not. Compounds that show no
toxicity in animals then advance to phase I trials, which will study whether the drug
is also safe in healthy humans and provide an initial idea for appropriate dosage.
In phase II, the efficacy of the drug is examined in parallel to potential side effects
to answer the question if it principally meets the expected performance. Phase II
presents the biggest hurdle with a transition success rate as low as 30%. Ultimately,
drug candidates enter clinical phase III in which the preliminary results found so far
need to be proofed and any adverse reactions monitored to make sure that the drug
really helps treating the disease [2, 11].
Starting from the generation of a lead compound assessment, and optimiza-
tion of pharmacokinetic properties and correlation to pharmacodynamic effects
increases in importance as one of the three major attrition causes among toxicity
and efficacy [8, 12]. In this sense, it is not surprising that the period between lead
and the clinical candidate is sometimes referred to as “valley of death” due to the
often occurring failures and dead ends during this time of the DDD process, which
results in high costs and missing deadlines [13].
In a long ongoing effort, more and more in silico techniques are being integrated
into several points of DDD with different purposes. In silico techniques can ease the
process of SAR assessment as well as the generation of compound series by guiding
combinatorial chemistry since they allow fast and easy evaluation of compounds
prior to synthesis from big libraries. For instance, combinatorial chemistry offered
an option to readily produce a broad range of potentially pharmaceutical active
small molecules in a short time, while SAR data in combination with complex math-
ematical algorithms, such as regression analyses based or machine-learning–based
approaches, allow to determine the potential effects of the analogues and deriva-
tive’s structures a priori [8].
Latter approach can save time and resources by eliminating in early stages
molecules that have predicted low efficacy against the target or to suggest the
next round of chemical modifications [14, 15]. Still, lead generation and/or
optimization will eventually also include in vivo testing after no toxic side effect
was shown in vitro. In vivo efficacy testing will be carried out as proof of concept
followed by PK assessment and ultimately animal models of human disease to
find correlations between preliminary data and potential performance later on in
humans [12].
86
ADME Profiling in Drug Discovery and a New Path Paved on Silica
DOI: [Link]
In silico ADME prediction aims to generate tools and models based on experi-
mental data to calculate in vivo behavior of compounds by finding quantitative
structure-property relationships (QSPRs), which connect structural information
to physical and chemical characteristics or even biological behavior (quantita-
tive structure-activity relationship; QSAR). Gained empirical data are then
related to descriptors/properties thereby supporting the process of hit-to-lead
optimization [10, 16].
When using in silico methods for prediction, it is important to keep in mind that
algorithms and tools applied are only models thus being only as good as the data
and idea they are based on. That implies a continuous experimental validation and
improvement as a basic principle that is supported by an interdisciplinary team. In
this sense, frequently used models include QSPR predictors, matched molecular
pair (MMP), and data trend analysis since they allow comparably easy application
and are based on a high amount of (end) point data. For instance, some experi-
ments offer highly convenient data but do not contribute much to model design,
whereas others show high variability but lead to impactful models. Considering
the nature of data, it is important to know which type can be used as input from
different sources (low variability biological and activity data or homogeneously cal-
culated chemical descriptors) in contrast to data that should only be used from one
source (Caco-2, MDCK). A sophisticated approach to generate reliable data or to
determine differences between individual experiments is to use assays with control
compounds [10]. The target property must be obtained under the same experimen-
tal condition and, in the best scenario, obtained from the same laboratory aiming to
avoid interlaboratory and interpersonal data noise [17].
Furthermore, the choice of the number and type of molecular descriptors has a
high impact, since it influences the accuracy and interpretability of the model. One
would expect that using the maximum number of descriptors would be beneficial,
but in reality, the risk of overfitting the data or losing the interpretability is a trade-
off. This leads to the point that it is fundamental for a “good” model to find the per-
fect compromise between quality and quantity. Nevertheless, it is crucial to test and
train a model and to evaluate its predictability by different means, such as statistical
measures and internal and external validation as recommended by organizations as
OECD [18], and also includes outlier analysis to reduce the noise in the model. An
extensive review of different adequate validation methods is discussed in [19].
As a result of newly achieved advances in computational capability, more com-
plex models and algorithms can now be applied. Despite this, it is still a challenge
to create a model for the pharmacokinetic and pharmacodynamic phenomena and
interactions within an organism as complex as a mammal, let alone humans [10].
Finally, notwithstanding the apparent linearity, the development of a new chemical
entity into a drug is an iterative process, even more, where modeling is concerned,
with data from failed attempts being integrated into the new predictions [13].
Pharmacology is a major part of the DDD process and describes the interaction
of an organism and the drug. It can be divided into two main branches: while phar-
macodynamics (PD) describes what the drug does to the body, pharmacokinetics
(PK) is interested in what the body does to the drug [20]. The main processes of PK
are absorption, distribution, metabolism, and excretion (ADME), finally comple-
mented by toxicity (ADMET). While ADME tries to maximize the pharmacological
performance of a small molecule, toxicology aims to ensure that it causes no harm
in any kind of side effect [21].
87
Drug Discovery and Development - New Advances
Table 1.
Tools for ADME evaluation.
88
ADME Profiling in Drug Discovery and a New Path Paved on Silica
DOI: [Link]
For instance, although the perfect approach of PK profiling would also reflect
the kinetics of drug administration and concentration at the site of action, most in
vivo systems rely on plasma sampling as a medium of drug equilibrium since it is
easily accessible. As a consequence, results are highly influenced by intrinsic and
extrinsic factors such as interpersonal variances as already stated above [20].
Each compound possesses individual physicochemical properties, such as solu-
bility or lipophilicity, which are influenced by biochemical properties of the body
as the different pH of tissues. Although they can be similar, each compound will
behave differently, and it is futile to address in vivo behavior without any prelimi-
nary knowledge of the basic PK parameters in vitro [26].
Furthermore, every PK assessment varies depending on the route of admin-
istration and requires different models and assays. While some routes depend on
absorption mechanisms like oral and transdermal administration, others (i.e.,
intravenous) directly target the bloodstream and the bioavailability is essentially
equal to 100%. Hereafter, we will discuss oral administration parameters of small
molecules as the most common form due to many advantages like reliability, safety,
price, their experimental approaches, and most common prediction modes [27, 28].
Passive transport across membranes is defined as permeability, which is depen-
dent on lipophilicity, since biological membranes are virtually lipid bilayers, and is
by far the most important transport for small molecules, especially in oral absorp-
tion [8, 24, 29]. Nonlipophilic compounds normally do not traverse membranes
passively, while highly lipophilic molecules run the risk to get stuck within the
membranes [30].
Properties utilized for measuring lipophilicity are the logarithm of the partition
coefficient (logP) and the distribution coefficient (logD) with the first not dif-
ferentiating between ionized and nonionized species. Both are normally applied for
n-octanol/water representing an organic and aqueous phase, respectively [21, 26].
Ionizability and lipophilicity provide a strong indication if a compound is likely
to be orally absorbed or not [21].
Ultimately, also the molecular size of the compound is involved in successful
absorption due to the aforementioned effects on permeability and solubility [31].
Usually, increasing molecular weight by adding new chemical moieties leads to
decreased solubility in aqueous solutions [32] and while big lipophilic compounds
partition passively along membranes (transcellular), small charged molecules can
also cross membranes via tight junctions (paracellular) [26]. For oral absorption in
terms of permeability, Lipinski and collaborators already proposed in 1997 [33, 34]
that orally active compounds should fit at least three of observed four parameters:
molecular weight < 500 g mol−1, logP < 5; number of hydrogen bond acceptors <10;
number of hydrogen bond donors <10; the well-known Lipinski’s rule of 5 (Ro5).
In other words, Ro5 stated a physicochemical space in which molecules outside its
domain has a low probability to become orally active. Other rules, as Veber rules
[35], Daina and Zoete [36], Egan and collaborators [37], Lovering et al. [38], and
Ritchie and colleagues’ [39] works for example, also included other properties as the
sum of hydrogen bond acceptor and donors, rotatable bonds count, polar surface
area, number of aromatic rings, and fraction of sp3 carbon atoms.
Despite the criticism and overinterpretation of Lipinski and derived rules, the
influence of physicochemical parameters on oral bioavailability and related param-
eters (as logP and aqueous solubility) is notable. Moreover, these rules are still being
employed nowadays in virtual screening campaigns aiming to reduce the number
of compounds from massively large available libraries (e.g., ZINC, which contains
more than 750 millions of compounds) [40–42]. Furthermore, those initial steps
instigate the generation of more complex models to predict not just oral bioavail-
ability but other PK-related parameters as Caco-2 permeability, aqueous solubility,
89
Drug Discovery and Development - New Advances
90
ADME Profiling in Drug Discovery and a New Path Paved on Silica
DOI: [Link]
Modifications in terms of the GSE have been proposed, for instance with the
SCRATCH model, which replaces the melting point by molar aqueous activity
coefficient, with comparable accuracy (R2 = 0.956, RMSE = 0.859 in a dataset of
883 compounds) [54]. Ali and collaborators suggested replacing the melting point
descriptor of the GSE with TPSA, aiming to overcome the issues with compounds
with high melting points and also to explicitly take into account the effect of polar
and polarizable atoms on the aqueous solubility [55].
The argument that real drugs are actually more soluble than drug-like molecules,
filtered by Lipinski’s rule of five [56], pointed out the studies in the direction of more
complex models. Indeed, nowadays, the quantitative structure-property relationship
(QSPR) models correlating the aqueous solubility with various molecular descriptors
are often employed. As an example, Chevillard et al. reported the use of a random
forest protocol to select the most accurate model among several available, both in
commercial or free software packages, for each compound [57]. They report that the
multimodel approach can enlarge the applicability domain given that more accurate
results for solubility prediction were obtained in comparison to using individual
models. This approach agrees with other reports that consensus of local QSAR
models can generate predictive workflows, especially for datasets with large struc-
tural diversity [58, 59]. It is worth noting that Lipinski himself recently revisited his
own rules [60], in vision of new potential classes of drugs, such as natural products,
peptide-like, and fragments, which, despite the validated effect, would defy the
original Ro5 limits.
91
Drug Discovery and Development - New Advances
92
ADME Profiling in Drug Discovery and a New Path Paved on Silica
DOI: [Link]
Finally, beyond the usual ADME parameters of interest in DDD, there are
several other unusual ones that also can be predicted; as examples, we here point
the permeability of the models for skin permeability, which evolved from simple
diffusion models based on molecular weight and n-octanol/water partition coef-
ficient [83, 84], until more sophisticated models, such as (non)-linear QSPR models
and even molecular dynamics simulation (as extensively reviewed by [85, 86]).
93
Drug Discovery and Development - New Advances
4.1 Absorption
Oral bioavailability is defined as the amount of drug that reaches the site of
action after oral administration and is influenced by factors like drug solubility and
dissolution, chemical and enzymatic stability in the gastric and intestinal lumen,
interacting luminal contents (food), gastrointestinal transit time, enterocyte per-
meability, and intestinal and hepatic metabolism [24]. Recently, bioavailability has
been also described as the rate and speed of the drug to reach systemic bloodstream,
considering the initial formulation as the starting point.
Oral administration includes a pharmaceutical phase—prior to PK and PD
phases—that comprises disintegration and dissolution of the dosage form. When
using oral dosage forms, the shape and chemical composition (e.g., tablets) play
an important role since they contribute to the time needed for disintegration and
dissolution.
Following the pharmaceutical phase, absorption is the first step in the phar-
macokinetic phase and is defined as the movement of the drug from the site of
administration to the bloodstream. The main properties determining the rate of oral
absorption for small molecules are permeability and solubility [87].
As such, the rate of dissolution and ionization, which are described by the
Noyes-Whitney and Henderson-Hasselbalch equation, respectively, is the key fac-
tors in lead optimization for oral administration and is complemented by lipophilic-
ity as an additional factor influencing membrane permeation and solubility of the
compound [31].
Dissolution can be expressed by a function of the aqueous solubility of a com-
pound, the surface area of the administered tablet (or the particles in other solid
formulation), and a specific dissolution rate constant. Altering any of these param-
eters directly affects the dissolution profile [26]. While solubility is an endpoint
value indicating the amount of a compound that is soluble in a solvent, dissolution
describes the kinetic process of a compound being solved in a solvent [88].
On the other hand, ionization reflects if a compound is present in the charged
or uncharged state and is at least influenced by two major parameters. The physico-
chemical property responsible for ionization is the pKa and describes the ionization
state of that entity at a given pH. It is also referred to as aqueous ionization constant
[30]. Thereby, it is directly influenced by the pH of the environment, the second
parameter, which drastically changes on the way through the GI tract, from about
pH 1 to 8 in the stomach and ileum, respectively.
The determination of the ionization state of a compound in the gastrointestinal
system (stomach, jejunum, ileum, and colon) is crucial for absorption since it
not only influences the solubility of a compound but also the lipophilicity and
permeability [26, 30, 89]. About 60–70% of all drugs (effective 1999) are ioniz-
able, which underlines the role that ionization plays in ADME assessment [30, 90].
While charged molecules easily dissolve in aqueous systems (GI tract), they do
not permeate membranes via passive diffusion and are reliant on active transport.
The contrary is true for uncharged molecules, which pass biological membranes
passively but show low solubility in aqueous solutions. Mechanisms of drug absorp-
tion include passive diffusion, active transport, and receptor-mediated endocy-
tosis, which are influenced by different factors and can themselves influence the
bioavailability.
Similar to model and prediction, the absorption of a drug is a complex process,
which is influenced not only by the physicochemical properties of drugs themselves
but also by the physiological state of the tissue in question. As such, there are a large
94
ADME Profiling in Drug Discovery and a New Path Paved on Silica
DOI: [Link]
number of prediction models available, which were generated based on the physico-
chemical properties involved in the absorption process, such as membrane perme-
ability and drug solubility. These models can help formulation scientists to optimize
drugs with poor absorption due to low aqueous solubility.
Initial absorption models can be separated into dispersion and compartmental
models [91]. While dispersion models treat the GI as a continuous system, with
variable pH and surface area, compartmental models take into account physi-
ological factors such as transporters. The compartmental absorption transit (CAT)
was one of the first models to regard distinct physiological properties, such as the
minimal absorption in the stomach and colon, while assuming some mathemati-
cal simplifications, such as the instant dissolution of the drug and linear kinetics
[92]. CAT was further modified as advanced CAT (ACAT), by including nonlinear
absorption kinetics and the effects of the first-pass metabolism. ACAT also consid-
ers the gastrointestinal tract as nine subsections, each with unique physicochemical
properties, such as pH, allowed solubility, particle size, and permeability [93].
Novel developments have included other absorption routes other than the GI,
which have been recently included in commercially available software, such as oral
absorption for the development of sublingual zolpidem tablets [94]. The absorption
constant (Ka, expressed in terms of h−1 min−1), or also called first-order absorp-
tion rate constant (to not be confounded with pKa), is employed in most of the
aforementioned models and is determined as a result from the changes in mass of
absorbable drug over time at the site of administration. Ka can be derived from the
decrease in the drug amount of absorbable present at the site of administration over
time; however, it is often indirectly determined by the drug amounts measured in
the blood and/or urine.
Along physicochemical models, which have a global application, machine-learning
techniques were extensively employed to model absorption (as comprehensively
reviewed by Kumar et al. [95]) and are inclined to be local models, since they are
mostly based on a small, homogeneous dataset that influences their applicability
domain.
After being absorbed and entering the circulatory system, the drug moves
reversibly between different compartments within the body, which is described
as distribution and influenced by several physicochemical properties of the
drug and biological factors of the body. One of the most important properties is
lipophilicity, and as such logP/logD, since it reflects the ability of the compound
to pass biological membranes to reach other sites, tissues, and organs within the
body [25]. Additional factors include phospholipid and (plasma) protein binding,
which reduces the free drug concentration within the body, can prevent the migra-
tion to the receptor side/side of action, and causes drug-drug interactions [25,
96]. Interestingly, binding to plasma proteins can also prolong the drug action by
releasing the drug over a longer period of time. It is also important to note that the
influence of lipophilicity on plasma protein binding is hypothesized to be higher for
acidic compounds than for bases, meaning that negative charges contribute highly
to plasma protein binding and prevent tissue binding, which leads to diminished
volumes of distribution (Vd, Eq. (3)). The Vd is the amount of drug that is freely
available in the blood, thus not bound to plasma proteins or other components [25,
97, 98]. Vd is an apparent volume that increases proportionally to the extravascular
drug binding and not an anatomically defined volume. Consequently, exten-
sive drug binding outside the bloodstream leads to increasing values of volume
distribution.
95
Drug Discovery and Development - New Advances
Volume distribution (Vd) is defined by the ratio between the amount of drug in
the body (A) and the drug concentration in plasma (C, comprising both free drug
and protein-bound drug):
A
Vd = __ (3)
C
The parameter describing protein binding is the plasma protein affinity con-
stant Ki. Many efforts to determine distribution led to chromatography-based meth-
ods, such as (RP-)HPLC to mimic n-octanol/water logP or lipophilicity to measure
distribution. In general, chromatographic methods are believed to resemble biologi-
cal partition processes more than octanol/water partition [25]. In the beginning,
stationary phases in (RP-)HPLC were either silica-based or polymer-based but both
had difficulties to reproduce logP and logD values despite several additives in the
mobile phases [99]. The introduction of biomimetic (stationary) phases coated with
human serum albumin (HAS), α1 acid glycoprotein (AGP), or immobilized artificial
membranes (IAM) revolutionized the methodology since they allowed a better
approximation of the biological system [25, 100].
A method to address plasma protein binding is the use of HSA and other plasma
proteins (e.g., α1 acid glycoprotein) coupled with RP-HPLC [25, 101]. On the other
hand, HPLC combined with IAMs is a popularly accepted technique for phospho-
lipid interaction and partition and several IAM columns are commercially available
for DDD projects. Both techniques represent good assay systems to model in vivo Vd
in high-throughput scale [98]. Problems with HPLC techniques, which are also true
for biomimetic phases, include the lack of a gold standard that is needed to calibrate
and later standardize results to make a comparison possible [25].
In vitro standard methods for unbound plasma fraction calculation include equi-
librium dialysis and ultrafiltration among several others as the two most commonly
used methods and are considered the gold standard for binding assessment [26].
To calculate the Vd “a priori”/nonexperimentally, plasma protein binding,
experimental logD and pKa are necessary. Then again, based on the Vd, the half-life
(t1/2) of a compound can be calculated [102]. Apart from protein binding, tissue
binding is also involved in the distribution of the compound. Generally, “tissue”
here comprises several components of the human body such as lipids, DNA, or RNA
and is also referred to as nonspecific binding [26].
In silico models to predict the Vd are often based on lipophilicity and solubility
descriptors, which correlate with the fractions of the drugs that are either bound to
plasma proteins or freely available. The work of del Amo et al. not just accurately
predicts Vd and unbound drug fraction but also compares the model’s performance
against the commercially available software VolSurf+ with comparable accuracy
(R2 = 0.70 and 0.71, respectively) [103].
Expanding these studies, the work of Lombardo and Jing generated a set of
models to predict the Vd in the steady state (Vss), using a dataset of 1096 diverse
compounds [104]. They compared models generated by linear (PLS) with nonlinear
(Random Forest) models, recommending the latter, with 33 descriptors, as the
optimal method for Vss prediction.
The Vd of drugs is greatly influenced by binding to plasma proteins with
several machine-learning–based models generated to predict this interaction.
Protein-protein interaction (PPI) information derived from molecular docking
was employed to derive a PPI-QSAR model for a small dataset of antibiotics (65
unique compounds), which resulted in an accurate model (R2 = 0.86 for the test set)
[105]. Additionally, global quantitative models using an array of classification and
regression models using physicochemical and molecular descriptors derived from
96
ADME Profiling in Drug Discovery and a New Path Paved on Silica
DOI: [Link]
a dataset of 794 compounds were shown to correctly classify the binding status of
the test set compounds and could be used as a prescreening [106]. Another recent
QSAR study using an extensively curated training set of 967 diverse pharmaceuti-
cals aimed to predict plasma protein-bound fractions (fb) using models generated
by six machine-learning algorithms with 26 molecular descriptors [107]. This study
is particularly interesting where the applicability domain is concerned allowing to
differentiate whether the classification derives from (un-)favorable regions.
del Amo et al. recently reported one of the first QSPR models to predict intra-
vitreal volume of distribution and clearance of small molecules [108]; the model
relies on the LogD and hydrogen bond capacity to understand phenomena such as
intraocular pressure and guide drug discovery. Complementarily, the prediction of
the drug passage through the blood-ocular barrier was described to be an important
factor to evaluate volume distribution in this organ [109].
Recently, as a novel approach bridging the animal experiments with human
results, it was shown that in PXB mice, a chimeric mice linage with a humanized
liver, plasma concentration-time profiles could be used to infer human’s compound
half-life [110].
Volume of distribution is also closely related to half-life and clearance param-
eters. As the Vd is a relative measurement of the free concentration of drug in the
blood, this same amount could be excreted by kidneys in the glomerular filtration
(clearance). Consecutively, the rate of clearance (discussed below in Excretion sec-
tion) directly influences the amount of available drug. Naturally, the concentration
of free drug that can bind its molecular target is related to the therapeutic dosage
and the half-life of the administered drug (as seen in Eq. (4)).
Half-life definition. Half-life is calculated by a ratio between the Napierian
logarithm multiplied by the volume of distribution (Vd) and renal clearance (CL):
[Link]
t1 / 2 = _______ (4)
CL
4.3 Metabolism
97
Drug Discovery and Development - New Advances
4.4 Excretion
Excretion is guided by one of the major clearance organs, and the assessment
of clearance behavior sometimes involves isolated organs or tissues [24]. Humans
rely on the kidney clearance as a major route for xenobiotic excretion, despite other
available routes such as feces, bile, sweat, and breath. The excretion pathways
directly impact the concentration of available drugs and are often measured in
terms of half-life and the initial administered dose.
The renal clearance of a drug is another important parameter, which is usually
employed to predict drug excretion. Experimentally, clearance is defined by the
drug concentration drug along a defined time of renal excretion by a linear equa-
tion (Eq. (5)).
Equation for renal clearance. m is the substance’s mass generation rate, K is the
clearance and C is the concentration at the time, and V is the volume where the drug
is distributed, or for systemic approaches the total body water.
dC = − K.C + m
V.___ (5)
dt
98
ADME Profiling in Drug Discovery and a New Path Paved on Silica
DOI: [Link]
Gombar et al. developed SVM- and MLR-based QSAR models to predict both
systemic clearance and apparent volume of distribution from intravenous data
[121] using as input structural fingerprints and electro-topological states (so-called
E-states), respectively. The model performed with high accuracy, despite the highly
diverse initial dataset employed for its generation, which points the importance of
those models in early steps of the drug-discovery pipeline.
Also, the work of Kusama et al. established a chemoinformatic-based classifica-
tion model to predict the major clearance pathways of 141 approved drugs based
on four physicochemical parameters: charge, molecular weight, lipophilicity, and
protein unbound fraction in plasma, resulting in a final model with an accuracy
of 88% [122]. This model approach was further refined by using support vector
machine and increasing the number of relevant descriptors [123]. In order to better
model the biotransformation processes, often the major triggers of excretion,
the work of Berellini et al. used ELASTICO (Enhanced Leave Analog-Structural,
Therapeutic, Ionization Class Out) to provide an appropriate sampling during the
validation process. Their partial least-square models resulted in a highly accurate
model derived from 754 compounds [124].
On another topic, ABCB1, also known as P-glycoprotein (P-gp or MDR1), is
a membrane protein member of the ATP-binding cassette (ABC) transporters
superfamily. Together with the hERG channel and CYP3A4, P-gp is one of the most
widely studied antitarget, where its inhibition could bring consequences for several
processes, such as the absorption, distribution, and excretion of drugs. Classical
studies used chemometric methods to describe bioavailability in terms of P-gp
and CYP enzyme activities, generating QSAR models based on 805 unique drug
molecules with high accuracy (R2 = 0.80 for the test set) [125]. Alternatively, an
approach to predict P-glycoprotein inhibition using molecular interaction fields,
derived from a literature collection of more than 1200 structures, generated a
pharmacophore model for competitive P-gp inhibition [126].
The most recent reported studies involving prediction of drug clearance, both
from human and rat hepatic in vitro systems, were based on microsomes, with a
recent emphasis on the use of hepatocytes. Wood et al. discuss the inherent limita-
tion of using human hepatocyte predictions, due to underprediction when com-
pared to in vivo clearance data, and the comments on the potential causes for those
divergences [127].
As the pinnacle of ADME in silico approaches, the holistic physiologically based
pharmacokinetic (PBPK) modeling was initially conceptualized by Teorell [128],
aiming to enable the prediction of drugs’ pharmacokinetic behavior in humans using
preclinical data. Recent PBPK models benefit from the large amount of available
ADME data not only to aid the drug discovery process and dose regiment selection
but also to guide the risk assessment for regulatory reviews [129]. PBPK models are
compartmentalized representations of the different organs, and each compartment
can be described by a specific tissue volume and blood flow rate, which communi-
cates with the blood (venous and arterial). Each organ/tissue has a unique volume,
permeability, and eliminating anatomical constants and terms, which are deter-
mined independently from the studied drug, while other physiological drug-specific
parameters are later incorporated, such as affinity toward plasma proteins, tissue-to-
plasma distribution rate, and even on target activity (Km, Vmax, or binding kinetics).
One of PBPK models’ important features is the perspective for the mechanistic
and prospective prediction of a drug’s pharmacokinetic profiles. The use of drug-
dependent parameters includes, but is not limited to, physicochemical properties,
solubility and permeability values, and also the role of individual enzymes and
transporters in the metabolism. Those parameters can be determined in vitro
99
Drug Discovery and Development - New Advances
ADME prediction
vNN-ADMET Public web server for ADMET property [Link]
prediction based on 15 nearest neighbor [132]
models.
Swiss-ADME Public web server for ADME property [Link]
prediction. It has a very unique LogP [133]
calculation ([Link]) based on free energy.
pkCSM ADME web server based on chemical [Link]
fragment similarity (the so-called [134]
graph-based signatures).
ADMETlab Web server using similarity-based [Link]
ADME calculator models and drug- [135]
likeness analyses.
Schrodinger— Calculates pKa; LogP; water [Link]
QikProp solubility—Schrodinger also offers QikProp, Schrödinger, LLC, NY, 2019
other tools for property calculation.
DDI-Predictor DDI-Predictor is able to make [Link]
quantitative predictions of drug
exposure even in cases where the
interaction has not been studied yet.
PBPK models and platforms
GastroPlus Comprises 10 different modules [Link]
including PBPK modeling and in
vitro vs. in vivo correlation, can be
parameterized for different disease
states and age groups.
PKSIM PBK modeling tool with integrated [Link] [Link]/products/
database of anatomical and [Link]
physiological parameters for humans,
mouse, rat, dog, and monkey. Can
model different scenarios depending on
the chosen building blocks.
Simcyp Incorporates databases of genetic, [Link]
physiological, and epidemiological
information to enable simulation of
different populations and species,
ultimately is able to predict ADME
parameters.
ADMEWORKS As a differential is able to predict [Link]
DDI Simulator drug-drug interactions using nonlinear materials_life_science/products/
models. ddi_simulator
Table 2.
Tools for ADME prediction and PBPK modeling.
100
ADME Profiling in Drug Discovery and a New Path Paved on Silica
DOI: [Link]
Early PBPK models, such as the work of Varma et al., described another layer of
complexity by including drug-drug interactions (DDI). The dosing time-dependent
model considering the interaction between repaglinide with rifampicin was able to pre-
dict repaglinide plasma concentrations along a day. The model also predicted the drug
interaction with other CYP3A4 and OATP1B1 inhibitors, which could result in further
DDIs. Reports of DDI leading to complications in patients with particular genotype
stimulated studies such as the one performed by Fermier et al. [131], where the effects
of polymorphic cytochromes provided the basis for a more accurate DDI prediction.
During recent years, larger molecules (LM) have gained in significance and
popularity, due to achievements and approvals, as new molecular entities. These
“biologics” are normally biotechnologically synthesized or recombinantly pro-
duced compounds of biological origin such as peptides, antibodies, and nucleic
acids [136]. From a historical perspective, drug discovery and development of
LMs are heavily delayed in comparison to SMs with their first approved entity
happening in the 1980s [137]. At about the same time, two major inventions
allowed huge progress in pharmacokinetics assessment of small molecules,
contributing to smaller drop-out rates in later DDD stages [136]. One of them was
the improved understanding of CYP450 mechanism and the other, the invention
of (HP)LC-MS technology, fueled the assessment of the ADME parameters. LMs’
discovery and development face many challenges, which demand high efforts
to overcome but also offer unique opportunities in comparison to those of small
molecules [138, 139].
The main differences between small and large molecules, despite the molecular
weight, the number of heavy atoms, and torsions, can be found in the physicochem-
ical properties, such as permeability, oral bioavailability, stability, specificity, and
immunogenicity [138, 139]. New parameters, unique for large molecules, are also of
interest, such as the physical particle size and the hydrodynamic radius, which has
a dramatic effect on the absorption. Both parameters are related to the overall shape
and correlate well with MW for globular proteins, but not necessarily for unstruc-
tured or highly modified entities. As a result, biologics are normally administered
parenterally, only targeting extracellular structures; they are also more likely to trig-
ger an immune response; and their production costs are considerably higher [139].
Interestingly, with the exception of the costs, these disadvantages can potentially
be circumvented by appropriate delivery systems, for example, nanoparticle-based
delivery to facilitate membrane permeation.
Other parameters, such as charges, which were previously modeled by pKa
in case of small molecules, are heavily heterogeneous in LMs. The charge can be
represented by the use of isoelectric points (pI), which are calculated from the
available amino-acid sequence, and surface charge, which can use individual pKa’s
and structural information to be inferred. Overall protein charge often influences
the biologic excretion [140], since negatively charged molecules undergo less renal
filtration disregarding size effect [141].
While representing difficulties in the development of new molecular entities,
the aforementioned properties also offer special advantages that small molecules
cannot cover. As such, LMs normally have longer t1/2, slower clearance, and higher
selectivity; are multifunctional; and rarely expose drug interactions [139]. Apart
from those, it was suggested that only 2–5% of the human genes can be targeted by
small molecules, offering a niche for LMs’ application against several diseases [138].
101
Drug Discovery and Development - New Advances
The increasing effort and development of new technologies, driven by the belief
in higher success rates, enabled the latest advances in the field [138]. For instance,
currently, peptide drugs only account for ~2% of the drug market but are in use in
a wide range of diseases such as acromegaly and multiple sclerosis, together with
different cancer types such as prostate and breast cancer.
Several other biologics are currently in use, namely monoclonal antibodies
(mAbs) and bispecific antibodies (bsAbs), as example agents that activate or
enhance the immunologic response. Of special interest in cancer therapy is a sub-
class of bsAbs, so-called bispecific T-cell engager (BiTEs), which can recruit CD3
cells at the tumor site by binding to both cell types thereby directing the immuno-
logical response [142].
Other interesting examples for biologics comprise hormones (e.g., insulin),
cytokines (such as erythropoietin, EPO; IL-1; IL-2; IL-6) [143], nucleic acids such
as siRNA (ONPATTRO) [144], and aptamers (Pegaptanib) [145]. While such a
broad spectrum of molecule classes offers also a wide range of treatments, at the
same time, it exacerbates the need for new developments since every molecule type
exhibits different properties. In the field of predicting the biologics activity against
specific targets, classical modeling tools, such as Monte Carlo sampling, genetic algo-
rithms, docking, and molecular dynamics simulation, were adapted or even devel-
oped anew to accommodate the specifics (as extensively reviewed by [146, 147]).
On the other hand, the absence of standard techniques to assess ADME proper-
ties hampers the PK profiling and thus further development [136]. In fact, the
current knowledge of LM pharmacokinetics is even impaired compared to the basic
knowledge of ADME principles for small molecules in the 1980s [136]. Although
the basic PK principles are similar between SMs and LMs, the specific mechanisms
influencing each step of ADME are different. To begin with, the route of adminis-
tration between them can differ, which leads to different mechanisms of absorption
and first-pass metabolism. Furthermore, LMs are not metabolized by CYPs but can
still trigger the release of pro-inflammatory cytokines leading to heavy side effects
known as cytokine storm [136, 139]. Also, other modifications play a role in biolog-
ics ADME, namely glycosylation, PEGylation, and neonatal Fc receptor (FcRn)
interactions [139, 148]. Unfortunately, up until now, most of the evaluation of those
factors is only addressed on in vivo level systems, which are not suited for HTS, are
expensive and labor intensive, and require longer bioethics evaluation.
In this regard, the development of in vitro and in silico methods to evaluate
ADME should be a high-profile goal. One of the main challenges will be to find
a way to integrate as many of the biologics into the process in order to facilitate
ADME assessment and guide large molecules’ DDD as already implemented for
their smaller counterparts.
6. Conclusions
The main difficulties in PK profiling lie in the high costs and comparable low
throughput of in vivo models. The extensive use of animals in DDD also raises ethi-
cal issues and is further affected since animal models not always translate readily to
the humans, especially in terms of metabolism [149, 150]. Furthermore, the advent
of combinatorial chemistry coupled with HTS for efficacy evaluation leads to an
explosion in the number to an extent that the classical PK assays could not compen-
sate [29, 47]. In vitro PK screens are supposed to offer a solution to the problem by
complementing in vivo assessment to reduce costs while increasing efficiency, but
they also suffer from shortcomings. In general, one must distinguish between two
main forms of in vitro systems: static and dynamic models. Only dynamic models
102
ADME Profiling in Drug Discovery and a New Path Paved on Silica
DOI: [Link]
are suited for PK evaluation because they allow variation of compound concentra-
tions, a key factor in pharmacokinetics. In this sense, diffusion-based dynamic in
vitro models offer a solution but still are quite limited in terms of high throughput
and costs. An alternative was presented by Lockwood and colleagues in the form
of a 3D-printed fluidic device utilizing trans-well technique generating dynamic in
vitro PK profiles also applicable for HTS infrastructure [149].
What distinguishes the DDD “then” and “now” is principally two main changes.
First, in the past, pharmaceutical companies as well as academic laboratories were not
that concerned with ADMET assessment in the early stages of drug discovery (hit and
lead generation) and only addressed PK from preclinical stages on forward. Instead,
HTS/HCS, genomics, and computational chemistry were high-profile areas. Today,
almost all pharmaceutical big-players have shifted pharmacokinetic profiling to dis-
covery phases. However, only the future will tell whether those changes will yield fruit.
Second, CADD became more and more part of the DDD pipeline in dif-
ferent stages facilitating fast screening of compounds in silico and supporting
QSAR. Although bioinformatics techniques already substituted many in vitro tests,
basically all of them require in vitro and/or in vivo validation and standardization
to guarantee trustable predictions. Another important aspect, recently addressed
by the work of Ferreria and Andricopulo [151], is the importance of translating
those models into well-structured and user-friendly (online) platforms that can be
accessed and used by the drug discovery community. Still, the efficacy and reli-
ability of computer simulations increase permanently and drastically, and many see
a future of solely virtual drug discovery. Thankfully, these failures resulted in the
consequence of addressing safety and efficacy concern earlier in the drug discovery
process, for instance, via in vitro screens to assess metabolic stability and absorption
properties and diminish failure rates later on [13].
Acknowledgements
Notes
Conflict of interest
103
Drug Discovery and Development - New Advances
Author details
© 2019 The Author(s). Licensee IntechOpen. This chapter is distributed under the terms
of the Creative Commons Attribution License ([Link]
by/3.0), which permits unrestricted use, distribution, and reproduction in any medium,
provided the original work is properly cited.
104
ADME Profiling in Drug Discovery and a New Path Paved on Silica
DOI: [Link]
References
[2] BIO Industry Analysis Reports. [10] Lombardo F, Desai PV, Arimoto R,
BIO n.d. Available from: [Link] Desino KE, Fischer H, Keefer
[Link]/bio-industry-analysis-reports CE, et al. In silico absorption,
[Accessed February 13, 2019] distribution, metabolism, excretion,
and pharmacokinetics (ADME-PK):
[3] Mullard A. Parsing clinical success Utility and best practices. An industry
rates. Nature Reviews Drug Discovery. perspective from the international
2016;15:447. DOI: 10.1038/nrd.2016.136 consortium for innovation through
quality in pharmaceutical development.
[4] Hughes JP, Rees S, Kalindjian SB, Journal of Medicinal Chemistry.
Philpott KL. Principles of early 2017;60:9097-9113. DOI: 10.1021/acs.
drug discovery. British Journal of jmedchem.7b00487
Pharmacology. 2011;162:1239-1249.
DOI: 10.1111/j.1476-5381.2010.01127.x [11] FDA. The Drug Development
Process. n.d. [Link]
[5] Van Norman GA. Drugs, devices, ForPatients/Approvals/Drugs/
and the FDA: Part 1: An overview [Accessed February 13, 2019]
of approval processes for drugs.
JACC: Basic to Translational Science. [12] Sinha S, Vohora D. Chapter 2—
2016;1:170-179. DOI: 10.1016/j. Drug discovery and development:
jacbts.2016.03.002 An overview. In: Vohora D, Singh G,
editors. Pharmaceutical Medicine and
[6] Paul SM, Mytelka DS, Dunwiddie Translational Clinical Research. Boston:
CT, Persinger CC, Munos BH, Lindborg Academic Press; 2018. pp. 19-32. DOI:
SR, et al. How to improve R&D 10.1016/B978-0-12-802103-3.00002-X
productivity: The pharmaceutical
industry’s grand challenge. Nature [13] Hefti FF. Requirements for
Reviews. Drug Discovery. 2010;9:203- a lead compound to become
214. DOI: 10.1038/nrd3078 a clinical candidate. BMC
Neuroscience. 2008;9:S7. DOI:
[7] DiMasi JA, Grabowski HG, Hansen 10.1186/1471-2202-9-S3-S7
RW. Innovation in the pharmaceutical
industry: New estimates of R&D [14] Andrews KA, Wesche D, McCarthy
costs. Journal of Health Economics. J, Möhrle JJ, Tarning J, Phillips L, et al.
2016;47:20-33. DOI: 10.1016/j. Model-informed drug development for
jhealeco.2016.01.012 malaria therapeutics. Annual Review
of Pharmacology and Toxicology.
[8] Shearer TW, Smith KS, Diaz D, Asher 2018;58:567-582. DOI: 10.1146/
C, Ramirez J. The role of in vitro ADME annurev-pharmtox-010715-103429
105
Drug Discovery and Development - New Advances
[15] Lombardino JG, Lowe JA. The role of Berlin, Heidelberg: Springer Berlin
the medicinal chemist in drug discovery— Heidelberg; 2013. pp. 1125-1138. DOI:
Then and now. Nature Reviews. Drug 10.1007/978-3-642-25240-2_48
Discovery. 2004;3:853-862. DOI: 10.1038/
nrd1523 [23] Vaddady PK, Lee RE, Meibohm B. In
vitro pharmacokinetic/pharmacodynamic
[16] Maltarollo VG, Kronenberger T, models in anti-infective drug
Espinoza GZ, Oliveira PR, Honorio development: Focus on TB. Future
KM. Advances with support Medicinal Chemistry. 2010;2:1355-1369.
vector machines for novel drug DOI: 10.4155/fmc.10.224
discovery. Expert Opinion on Drug
Discovery. 2019;14:23-33. DOI: [24] Dowty ME, Messing DM, Lai Y,
10.1080/17460441.2019.1549033 Kirkovsky L. Adme. ADMET for
Medicinal Chemists. John Wiley &
[17] Huang J, Fan X. Why QSAR fails: An Sons, Ltd; 2010. pp. 145-200. DOI:
empirical evaluation using conventional 10.1002/9780470915110.ch4
computational approach. Molecular
Pharmaceutics. 2011;8:600-608. DOI: [25] Valkó KL. Lipophilicity and
10.1021/mp100423u biomimetic properties measured by
HPLC to support drug discovery.
[18] Validation of (Q )SAR Models— Journal of Pharmaceutical and
OECD. n.d. [Link] Biomedical Analysis. 2016;130:35-54.
chemicalsafety/risk-assessment/ DOI: 10.1016/[Link].2016.04.009
[Link] [Accessed
February 13, 2019] [26] Fan J, de Lannoy IAM.
Pharmacokinetics. Biochemical
[19] Gramatica P, Sangion A. A Pharmacology. 2014;87:93-120. DOI:
historical excursus on the statistical 10.1016/[Link].2013.09.007
validation parameters for QSAR
models: A clarification concerning [27] Amidon GL, Lennernäs H, Shah VP,
metrics and terminology. Journal of Crison JR. A theoretical basis for a
Chemical Information and Modeling. biopharmaceutic drug classification:
2016;56:1127-1131. DOI: 10.1021/acs. The correlation of in vitro drug product
jcim.6b00088 dissolution and in vivo bioavailability.
Pharmaceutical Research.
[20] Abdel-Rahman SM, Kauffman RE. 1995;12:413-420
The integration of pharmacokinetics
and pharmacodynamics: Understanding [28] Lennernäs H. Human intestinal
dose-response. Annual Review permeability. Journal of Pharmaceutical
of Pharmacology and Toxicology. Sciences. 1998;87:403-410. DOI:
2004;44:111-136. DOI: 10.1146/annurev. 10.1021/js970332a
pharmtox.44.101802.121347
[29] Avdeef A. Physicochemical profiling
[21] Hodgson J. ADMET—Turning (solubility, permeability and charge
chemicals into drugs. Nature state). Current Topics in Medicinal
Biotechnology. 2001;19:722-726. DOI: Chemistry. 2001;1:277-351
10.1038/90761
[30] Comer J, Box K. High-throughput
[22] Tietgen H, Walden M. measurement of drug pK a values
Physicochemical properties. In: Vogel for ADME screening. Journal of
HG, Maas J, Hock FJ, Mayer D, editors. the Association for Laboratory
Drug Discovery and Evaluation: Automation. 2003;1:55-59. DOI: 10.1016/
Safety and Pharmacokinetic Assays. S1535-5535(04)00243-6
106
ADME Profiling in Drug Discovery and a New Path Paved on Silica
DOI: [Link]
[38] Lovering F, Bikker J, Humblet C. [44] Kibbey CE, Poole SK, Robinson B,
Escape from flatland: Increasing Jackson JD, Durham D. An
107
Drug Discovery and Development - New Advances
[50] Valkó K. Chapter 12. Measurements [58] Bergström CAS, Wassvik CM,
of physical properties for drug design in Norinder U, Luthman K, Artursson P.
industry. In: Valkó K, editor. Handbook Global and local computational models
of Analytical Separations. Vol. 1. for aqueous solubility prediction
108
ADME Profiling in Drug Discovery and a New Path Paved on Silica
DOI: [Link]
[65] Jensen JH, Swain CJ, Olsen L. [73] Kansy M, Avdeef A, Fischer H.
Prediction of pKa values for druglike Advances in screening for membrane
molecules using semiempirical quantum permeability: High-resolution PAMPA
chemical methods. The Journal of for medicinal chemists. Drug Discovery
109
Drug Discovery and Development - New Advances
110
ADME Profiling in Drug Discovery and a New Path Paved on Silica
DOI: [Link]
111
Drug Discovery and Development - New Advances
112
ADME Profiling in Drug Discovery and a New Path Paved on Silica
DOI: [Link]
113
Drug Discovery and Development - New Advances
[130] Madden JC, Pawar G, Cronin MTD, [137] Ortho Multicenter Transplant
Webb S, Tan Y-M, Paini A. In silico Study Group. A randomized clinical
resources to assist in the development trial of OKT3 monoclonal antibody
and evaluation of physiologically- for acute rejection of cadaveric renal
based kinetic models. Computational transplants. The New England Journal
Toxicology. 2019;11:33-49. DOI: of Medicine. 1985;313:337-342. DOI:
10.1016/[Link].2019.03.001 10.1056/NEJM198508083130601
114
ADME Profiling in Drug Discovery and a New Path Paved on Silica
DOI: [Link]
115
Section 4
117
Chapter 7
Abstract
1. Introduction
1.1 Tyrosine
119
Drug Discovery and Development - New Advances
Figure 1.
The steps of tyrosine degradation.
2.1 History
2.2 Pathophysiology
120
Successive Drug Therapy for a Very Rare Autosomal Diseases
DOI: [Link]
Figure 2.
Steps of tyrosine and its inborn enzymatic error of metabolism causing tyrosinemia type I, II, and III.
121
Drug Discovery and Development - New Advances
2.4 Diagnosis
122
Successive Drug Therapy for a Very Rare Autosomal Diseases
DOI: [Link]
Figure 3.
The structure of nitisinone (NTBC).
Sven Lindstedt recognized the potential value of NTBC for the treatment of
HT-1. By blocking the proximal tyrosinemia pathway, NTBC minimizes the forma-
tion of FAA and maleylacetoacetic. It was this keen insight that led to the original
clinical trial with five patients, which documented the rapid reversal of clinical
symptoms [58].
NTBC dosing should be sufficient to completely suppress plasma and urine SA
detection and normalize liver and renal function. SA, in either plasma or urine,
should be below detectable limits (or within the limits of normal established by
123
Drug Discovery and Development - New Advances
the reference laboratory). The dose of NTBC should be increased if the SA level
increases once patient adherence has been confirmed [43].
The standard recommended dosage of NTBC is 1 mg/kg body weight [43, 59, 60].
The half-life of NTBC has been measured in healthy human subjects and found
to be approximately 54 h [61]. Because of this long half-life, a single daily dose of
NTBC is satisfactory for maintaining inhibition of HPD [62, 63].
In the evaluation of its safety profile, rats and dogs exposed to NTBC developed
elevated plasma tyrosine levels and ocular lesions. The ocular lesions (keratopathy)
were caused by tyrosine crystals within the cornea, which on cessation of the diet
recovered [43, 57].
The FDA approved NTBC in January 2002 [64].
The combined nitisinone and low phenylalanine and tyrosine diet treatment
should be initiated as soon as possible following the diagnosis of HT-1, to maintain
PAA concentrations within the treatment range. Phenylalanine must be restricted
in the diets of affected patients since approximately 75% of dietary phenylalanine
is hydroxylated to form tyrosine [65, 66]. The combined diet restriction and NTBC
treatment resulted in a greater than 90% survival rate, normal growth, improved
liver function, prevention of cirrhosis, correction of renal tubular acidosis, and
improvement in secondary rickets [67–69].
Before the availability of nitisinone for the treatment of tyrosinemia type I, the
only definitive therapy was liver transplantation. The first case of HT-1 treated with
liver transplant was in 1978 performed by Fisch and his colleagues [70]. The patient
died 3 months later, but the biochemical derangements improved. Subsequently,
the use of liver transplants in HRT cases has increased, and the benefits appear to be
confirmed [71].
Liver transplantation should be reserved for those children who (1) have a severe
liver failure at clinical presentation and fail to respond to nitisinone therapy or (2)
have documented evidence of malignant changes in hepatic tissue [72]. Transplant
recipients may also benefit from low-dose (0.1 mg/kg/day) nitisinone therapy to
prevent continued renal tubular and glomerular dysfunction resulting from the
persistence of succinylacetone in the plasma and urine [73].
Diet restriction for the treatment of HT-1 patients was introduced by Halvorsen
and Gjessing in 1964 [19] and for a long time was the only treatment available. It
had a beneficial effect on the renal tubular defects but did not cure the liver disease.
A girl with HT-1, diagnosed at 6 months of age, was treated with a diet restricted in
phenylalanine and tyrosine. At 9½ years of age, she developed an acutely enlarged
liver and spleen, and the diagnosis of hepatocarcinoma was made [70].
124
Successive Drug Therapy for a Very Rare Autosomal Diseases
DOI: [Link]
We are reporting three Saudi siblings who have diagnosed as patients with HT-1.
They are living in Najran city, the southern province of Saudi Arabia.
The first case is 27-year-old male patient. He has been diagnosed at the age of
4 months in Great Ormond Street Hospital (London-UK) and treated by nitisinone
in combination with tyrosine and phenylalanine-free diet. For the next 15 years,
he used to visit the clinic for regular checking. The final report showed that he
has good general health, with normal liver and renal function test and normal
alpha-fetoprotein. Since then and until now, he attends the National Laboratory for
Newborn Screening, Department of Genetics, King Faisal Specialist Hospital and
Research Centre for a routine checkup. Still, he is on nitisinone and diet restriction.
He has graduated from Najran Technical College, and now he is doing very well
in his job as the vice director of staff affairs at Al-Ghad International College for
Applied Medical Sciences in Najran-KSA.
The second case is his brother, who is 20 years old. He has been diagnosed in
the first few days after birth by the National Laboratory for Newborn Screening,
Department of Genetics, King Faisal Specialist Hospital and Research Centre,
Riyadh, Saudi Arabia. Since then he has been treated by nitisinone in combination
with tyrosine and phenylalanine-free diet. He is having good general and mental
conditions. He is now a 2nd-year university student and doing well in his study.
The third case is their younger brother, who is 15 years old. Diagnosed since
birth by the National Laboratory for Newborn Screening, Department of Genetics,
King Faisal Specialist Hospital and Research Centre, Riyadh, Saudi Arabia. Since
then, he has also been put on nitisinone in combination with tyrosine and phenyl-
alanine-free diet. He is also having a good general and mental conditions. He is a
secondary school student, and he is also doing very well in his study.
The three brothers are now on nitisinone in a dose of 1 mg/kg/day. It is marketed
under the brand name Orfadin by the company Swedish Orphan Biovitrum (Sobi).
Also, they are on tyrosine and a phenylalanine-free diet supplemented with HCU
Anamix Junior LQ [Nutricia Advanced Medical Nutrition Company].
They have a regular visit to the National Laboratory for Newborn Screening,
Department of Genetics, King Faisal Specialist Hospital, and Research Centre,
Riyadh, Saudi Arabia for regular checking, and the last visit was a few weeks ago.
Their elder two siblings (sister and brother) died in their early life of unknown
cause. The sister died at the age of 26 months. From the history of her mother, it
seems that the daughter developed abdominal distension, and their doctors told
her that the baby has hepatosplenomegaly. She also developed jaundice and became
reluctant to milk or food. Until the day of death at the age of 26 months, she could
not sit or walk; she was very week as her mother described.
Then after, they have a boy who developed after birth hepatomegaly, and
jaundice, abdominal distension with very thin extremities and died early after birth
at the age of 4 months, as his mother said.
125
Drug Discovery and Development - New Advances
After the birth of the third baby (our first case), they went to Great Ormond
Street Hospital and diagnosed as having tyrosinemia-1. At the same time, the
parents diagnosed as a carrier of the disease.
The clinical history of their previous dead siblings suggested that they did have
undiagnosed HT-1. Furthermore, the diagnosis of the parents that they are a carrier
of the disease is highly supporting that the death of their two children early life was
due to the lack of diagnosis and treatment of HT-1.
3. Discussion
HT-1 is a rare but clinically severe and fatal inborn error that principally affects
the liver, kidney, and peripheral nerve [32]. In general, the most diagnosed patients
of inborn errors of metabolism (IEMs) including HT-1, were born from consan-
guineous married parents. As HT-1 is a rare inherited autosomal recessive disorder,
it explains why it is more common in population with a high rate of consanguineous
marriages, such as in United Arab Emirates, Oman, Kuwait, and Saudi Arabia,
in which the rate of consanguineous marriages reaches up to approximately 60%
[33, 77–80], and even the first reported case of HT-1 in Japan (1957) was a child
from parents of consanguineous marriage [13, 14]. Furthermore, in our study, the
reported three family cases are also born from a consanguineous married parent.
In the sixties and early seventies HT-1 patients were treated by phenylalanine-
and tyrosine-free died; it delayed the mortality and morbidity for a few years but
it did not prevent the development of hepatic failure, renal complications and
hepatocellular carcinoma (HCC) even if begun within the first month of life [43].
Then after, treatment with hepatic transplantation was the only option for survival.
But it has been reported that the development of HCC was observed in 17–37% of
affected children [81, 82]. Furthermore, after transplantation, urine and plasma
SA decreases but is not completely suppressed [43, 63], and even other scientists
reported that plasma succninylacetone is persistently raised after liver transplanta-
tion [73] presumably because of continued production in the kidneys which could
cause damage to the liver and the kidney. Also, some of the patients died from
the complications of hepatic transplant, whether the surgical or the immuno-
suppressive drug complications [43, 63, 69]. Therefore, in 1992, the introduction
of nitisinone in combination with diet restriction was the ideal therapy for HT-1
patients especially if started in their early neonatal days [58].
The three HT-1 patients in our study used nitisinone in combination with
phenylalanine and tyrosine diet restriction from early days of their neonatal life and
till now, which is 27, 20, and 15 years, respectively. They are not only still alive but
also doing very well in their living.
These results justify implementing prevention programs that incorporate
genetic counseling and neonatal diagnostic screening tests, especially in the sus-
pected families of consanguineous marriages to detect the neonatal patients with
HT-1 as early as possible and then to start treatment which will minimize the lethal
consequences of the disease.
All subsequent children of the parents of a child with HT-1 should have urine
and blood succinylacetone analyzed as soon as possible after birth to enable the ear-
liest possible diagnosis and initiation of therapy. Early detection of newborn babies
with HT-1, followed by prompt treatment with nitisinone in combination with a low
phenylalanine and tyrosine diet has improved the survival to over 90% and resulted
in normal growth, improved liver function, prevention of cirrhosis, correction of
kidney disease, and improvement in rickets [41, 42, 83, 84]. In 2012, Larochelle
126
Successive Drug Therapy for a Very Rare Autosomal Diseases
DOI: [Link]
et al. reported that patients who receive nitisinone treatment before 1 month had no
detectable liver disease after more than 5 years [69].
These data suggested that early neonatal diagnosis of HT-1 and treatment with
nitisinone and diet restriction not only keep the survival of the patients but also
keep them in good general, physical, and mental conditions.
4. Conclusion
Nitisinone (NTBC) has been used since 1992 and proves to be an effective and
safe pharmacological treatment for HT-1 in combination with phenylalanine- and
tyrosine-free diet.
In this paper, we are reporting three cases (brothers) treated safely and succes-
sively with NTBC in combination with diet. All of them are in very good condi-
tions. The elder brother is on NTBC since 27 years ago. He is one of the few cases
worldwide treated since 1992 and till now, and he is living with a very good general
health.
HT-1 is not only a rare and fatal autosomal disease, but it is a very rare genetic
disease that can be successfully, effectively, and safely treated by drug therapy,
which is nitisinone (NTBC).
4.1 Recommendation
None.
Ethical consideration
The protocol was submitted and approved by the Research Ethical Committee of
Al-Ghad International College for Applied Medical Sciences-Najran, KSA.
Consent
Written informed consent was obtained from the family for the anonymized
information to be published in this article.
127
Drug Discovery and Development - New Advances
Author details
© 2019 The Author(s). Licensee IntechOpen. This chapter is distributed under the terms
of the Creative Commons Attribution License ([Link]
by/3.0), which permits unrestricted use, distribution, and reproduction in any medium,
provided the original work is properly cited.
128
Successive Drug Therapy for a Very Rare Autosomal Diseases
DOI: [Link]
References
129
Drug Discovery and Development - New Advances
liver, renal tubular defects with rickets From disease to cancer. Acta Biochimica
and impaired tyrosine metabolism. Acta Polonica. 1996;43:209-216
Paediatrica Scandinavica. 1964;53:18-32
[28] Jorquera R, Tanguay RM. The
[19] Halvorsen S, Gjessing LR. Studies on mutagenicity of the tyrosine metabolite,
tyrosinosis. 1. Effect of low-tyrosine and fumarylacetoacetate, is enhanced by
low-phenylalanine diet. British Medical glutathione depletion. Biochemical and
Journal. 1964;2:1171-1173 Biophysical Research Communications.
1997;232:42-48
[20] Gjessing LR, Halvorsen S.
Hypermethioninemia in acute [29] Mitchell G, Bartlett DW, Fraser TE,
tyrosinosis. Lancet. 1965:1132-1133 Hawkes TR, Holt DC, Townson JK, et al.
Mesotrione: A new selective herbicide
[21] Scriver CR, Clow CL, for use in maize. Pest Management
Siverberg M. Hypermethioninemia Science. 2001;57(2):120-128
in acute tyrosinosis. Lancet.
1966;1:153-154 [30] Rashed MS, Al-Ahaidib LY,
Al-Dirbashi OY. Tandem mass
[22] Larochelle J, Mortezai A, spectrometric assay of succinylacetone
Belanger M, Tremblay M, Claveau JC, in urine for diagnosis of
Aubin G. Experience with 37 infants hepatorenaltyrosinemia. Analytical
with tyrosinemia. Canadian Medical Biochemistry. 2005;339:310-317
Association Journal.
1967;97(18):1051-1054 [31] Croffie JM, Gupta SK, Chong SKF,
Fitzgerald JF. Tyrosinemia type 1 should
[23] Nakamura K, Ito M, Shigematsu Y,
be suspected in infants with severe
Endo F. Tyrosinemia type I in Japan: coagulopathy even in the absence of
A report of five cases. Advances in other signs of liver failure. Pediatrics.
Experimental Medicine and Biology. 1999;103(3):675-678
2017;959:133-138
[32] Russo PA, Mitchell GA,
[24] Lindblad B, Lindstedt S,
Tanguay RM. Tyrosinemia: A review.
Steen G. On the enzymic defects in
Pediatric and Developmental Pathology.
hereditary tyrosinemia. Proceedings
2001;4:212-221
of the National Academy of Sciences
of the United States of America.
[33] Al-Shamsi A, Hertecant JL,
1977;74:4641-4645
Al-Hamad S, Souid AK, Al-Jasmi.
[25] De Laet C, Dionisi-Vici C, Mutation spectrum and birth prevalence
Leonard JV, McKiernan P, Mitchell G, of inborn errors of metabolism
Monti L, et al. Recommendations for among Emiratis: A study from Tawam
the management of tyrosinaemia type 1. Hospital Metabolic Center, United Arab
Orphanet Journal of Rare Diseases. Emirates. Sultan Qaboos University
2013;8:8 Medical Journal. 2014;14(1):e42-e49
[26] Berube D, Phaneuf D, Tanguay RM, [34] Udawat P, Gupta RK, Shambhavi,
Gagne R. Assignment of the Sitaraman S. Case Report. Tyrosinemia
fumarylacetoacetate hydrolase gene to type I: A case report and review of the
chromosome 15q23-15q25 (Abstract). literature. Indian Journal of Medical
Cytogenetics and Cell Genetics. Case Reports. 2016;5(1):9-12
1989;51:962
[35] Bliksrud YT, Brodtkorb E, Backe PH,
[27] Tanguay RM, Jorquera R, Poudrier J, Woldseth B, Rootwelt H. Hereditary
St Louis M. Tyrosine and its catabolites: tyrosinaemia type I in Norway:
130
Successive Drug Therapy for a Very Rare Autosomal Diseases
DOI: [Link]
131
Drug Discovery and Development - New Advances
132
Successive Drug Therapy for a Very Rare Autosomal Diseases
DOI: [Link]
133
Drug Discovery and Development - New Advances
134
Chapter 8
Abstract
Drug discovery and development advances in the last decades allowed to find
a treatment for many preventable diseases. However, all too often, children are
excluded from these progresses since most of the new medicines have been discov-
ered and developed for the adult population. Even if paediatricians routinely give
drugs to children ‘off-label’, researchers have demonstrated that children do not
respond to medications in the same way as adults. Furthermore, certain specific
disorders are unique to children or occur in children differently than in adults.
Besides specifically testing medicines in children in proper clinical studies taking
into due account the peculiarity of this population, there is a growing recognition
of the need to develop paediatric medicines having in mind the specificities of this
vulnerable population. In this chapter, we will provide an overview on the drug
discovery and development path for children highlighting challenges and new fron-
tiers of each phase from the discovery to the preclinical and clinical development
as well as we will provide a slightest hint about paediatric biomarkers discovery,
age-appropriate formulation, pregnancy, and perinatal pharmacology and in silico
pharmacology. Finally, pricing and reimbursement policies for medicines and new
and existing research initiatives in the field will be discussed.
1. Introduction
135
Drug Discovery and Development - New Advances
The lack of a regulatory framework that obliged to test medications in the pae-
diatric population taking into account the specificities of children and the ethical
concerns behind resulted in several examples of therapeutic tragedies in paediatric
patients. A new liquid formulation of the antibiotic sulphanilamide was developed
in 1938 to allow oral dosing for paediatric patients who could not swallow the tablet
form. Unluckily, the solvent used to dissolve the active substance was a toxin that
caused many adverse events with a 30% mortality rate [3]. And again, Thalidomide
was marketed in Europe in the late 1950s for the treatment of nausea in pregnant
women causing severe birth defects in thousands of children including severe
shortening of the extremities, malformations of ears, heart, intestines and other
structures, depending on the embryologic stage at the time of exposure [4].
These tragedies are just an example of the high risk to which children have been
exposed for years and have led to the increasing awareness that new medications
for children should be carefully studied before they could be approved, defining
the proper requirements and ethical issues to guarantee efficacious and safer drugs
for children. As a consequence, regulations have been adopted independently in
the most developed countries, but in accordance with unified guidelines suggested
by the International Conference on Harmonisation of Technical Requirements for
Registration of Pharmaceuticals for Human Use (ICH), an organisation working on
the harmonisation of pharmaceutical regulatory requirements within the European
Union (EU), Japan and the United States (US) [5].
The European Paediatric Regulation was adopted in 2006 and entered into force
in 2007 [6] imposing to pharmaceutical companies developing drugs of potential
interest for children to prepare a paediatric investigational plan (PIP) to obtain a mar-
keting authorisation for an indication in adults, unless they were granted a product-
specific waiver by the Paediatric Committee of the European Medicines Agency
(EMA), for example if the indication does not occur in children [7]. The paediatric
regulation has defined rules concerning the development of medicinal products for
paediatric use and introduced rewards and incentives for the development of paediat-
ric drugs (i.e. the paediatric-use marketing authorization—PUMA) [5].
In the US, two main acts have complemented each other ruling the evaluation of
drugs in infants and children and increasing the paediatric clinical studies and drug
labelling for children: the PREA of 2003 [8] and the BPCA of 2002 [9], both amended
in the FDAAA of 2007 [10]. A different approach has been taken in Japan, more
focused on premiums granted to pharmaceutical companies as rewards for developing
paediatric medicines without a regulatory framework specifically addressing paedi-
atric clinical research. As an effect of these premiums, the price of those drugs is not
reduced as normally occurred every 2 years in the Japanese system [11].
2. Paediatric diseases
136
Challenges and New Frontiers in the Paediatric Drug Discovery and Development
DOI: [Link]
Obviously, childhood is the period of life when the physiological and physical
changes are the most important and the fastest. Physiological systems and functions
are immature in neonates at birth with the degree of immaturity depending on
gestational age. These systems develop progressively and changes can be observed,
for example, in gastrointestinal motility and function, body composition and size,
activities of transporters and metabolism enzymes, and renal function. The process
is dynamic and nonlinear with progressive rapid growth and maturation in the first
weeks/months of life, and slower thereafter. These developmental changes affect
drug disposition, as discussed later, with differences among neonates, children,
adolescents and adults [18, 19].
Therefore, defining the paediatric population is a very complex issue since it rep-
resents an extremely heterogeneous population. To address the peculiarity of each
age group and to provide guidance for regulatory and clinical matters, the interna-
tional regulation on paediatric clinical trials [20] has described four subsets: pre-
term and term neonates (0–27 days), infants (1–23 months), children (2–11 years)
and adolescents (12–18 years). In addition, the recently revised EMA guideline
‘Ethical considerations for clinical trials on medicinal products conducted with
minors’ issued on September 2017 has further sub-divided the age group 2–18 years
into pre-schoolers (2–5 years), schoolers (6–9 years) and adolescents (10–18 years)
[5, 21]. It has to be underlined that this definition, although useful to unify the
system of rules and law in this field, does not always reflect the maturity of the
child, which is something that is generally recognised as crucial aspect to be taken
into account during the conduct of paediatric clinical trials [5].
When it comes to talking about paediatric diseases, we cannot exclude the rare
disease field since many rare diseases are diagnosed during childhood. Rare diseases
include a very heterogeneous group of disorders, affecting any body system. A
disease is defined rare if it affects fewer than 1 in 2000 people in Europe and fewer
than 200,000 people in the United State. A high percentage of rare diseases (about
80%) affect children, and in 50% of cases, all rare diseases are characterised by
a childhood-onset with a significant impact on the well-being of the patients and
families [22–24].
Although rare diseases have by definition a low prevalence, with some having a
single identified case worldwide, collectively they affect about 6–8% of the human
population with a number of diseases recognised as rare comprised between 6000
and 8000 diseases [25].
Despite the high impact they have on the worldwide population, few treatments
are available on the market. Drug development for rare diseases poses unique scien-
tific and ethical challenges, most of which in common with the obstacles described
137
Drug Discovery and Development - New Advances
in this chapter for the paediatric population. Since they affect a small population,
heterogeneous and widely dispersed, it is more difficult to enrol enough patients in
clinical studies and pharmaceutical company shows a scarce interest in this field for
the low return they may have.
Moreover, considering the high incidence and prevalence during the childhood, the
ethical issue is predominant in this field. And additional challenges may result from
the frequently progressive, life-limiting or life-threatening nature of these diseases.
As described below, new approaches in all the phases of the drug development
process may offer valuable solutions to overcome these difficulties in the rare
diseases as well as paediatric diseases field.
Drug discovery and development path represents the long process starting with
the identification of new target molecules (discovery phase), going through stud-
ies on microorganisms and animals (preclinical development) and finally testing
the new medicines in the target population (clinical development) to bring them
to the market (authorization and commercialization). Considering the differences
between children and adults above mentioned, a new drug to be used in children
should be specifically tested in children themselves in controlled clinical studies.
At the same extent, medicines for children should be developed having in mind the
specificities of this vulnerable population starting from the very initial phase of
discovery.
138
Challenges and New Frontiers in the Paediatric Drug Discovery and Development
DOI: [Link]
Xu et al. provided an overview on the in vitro models used to study paediatric brain
tumours underlined that in the initial drug screening for new therapies, it is critical
to use cell lines more closely related to the tumour and organism being studied. The
authors listed 60 paediatric brain tumour cell lines reported in the literature, of which
only a small number can be obtained from central repositories such as ATCC [29] or
Children’s Oncology Group (COG) [30], thus rendering more difficult for the research
community to have access to the most adequate cell lines [31].
Considering these findings, novel preclinical models should be evaluated as plat-
form for drug discovery for paediatric diseases, such as induced pluripotent stem
cells (iPSCs) or innovative techniques including organoids and organs-on-a-chip.
Disease-specific iPSCs represent a promising platform to understand pathologi-
cal progression in patient-derived cells presenting many advantages: iPSCs are an
unlimited source of patient-specific cells for drug testing and for the development
of personalised medicine [32]. Advances in human pluripotent stem cell (hPSC) or
tissue-resident adult stem cell (AdSC) research have led to the possibility to mimic
any tissue in the human body through three-dimensional (3D) model including
organoids and organs-on-a-chip that can be used as in vitro screening models [33].
However, to confirm the adherence of these in vitro models with their normal coun-
terparts in vivo, we need a much deeper understanding of the physiology of human
development than what is currently available.
In addition, as regarding the animal models, the number of comprehensive
studies describing the normal development of different physiological systems and
processes in laboratory animals from molecular to system levels is very limited,
and such studies usually do not exist in animal models of paediatric diseases. Thus,
questions of comparability of developmental stages across species continue to
create debate. The need to use juvenile animal models will be better discussed in the
following section.
In addition to the need of developing cellular and animal models more suitable
to study paediatric diseases and the instruments to work with immature animals,
all the new emergent technologies should be timely applied to the paediatric drug
discovery in order to speed up the pharmacological research, including pluripotent
stem cell, 3D cell cultures, target validation, patient-derived cell assays, micro-
fluidics, high-throughput cell image analysis, non-invasive drug delivery systems
and devices to measure drug safety or efficacy non-invasively.
139
Drug Discovery and Development - New Advances
and excretion (ADME) [40]. The rate and extent of the bioavailability of a drug
may vary as a consequence of the development changes that occur in absorptive sur-
faces, especially the gastrointestinal tract. Dissimilarities have also been reported
in drug metabolism, transporters expression, biliary function and renal clearance,
resulting in differences in drug disposition and elimination [41].
Similarly to PK profile, PD profile is also affected by human development and
drug targets may vary under developmental control: their level of expression,
affinity or activity may diverge according to the patient’s age, resulting in variable
drug responses depending on patients’ age group. This is particularly important
in younger infants, more vulnerable to drug toxicity and related adverse events by
modifying drug therapeutic windows [42].
Another aspect to be taken into due account is represented by the effect of the
ontogeny and genetic variation interactions on drug response, known as pharma-
cogenetics [43]. Several pharmacogenetics studies have indeed demonstrated the
differences in response to drugs between children and adults [44].
To take into consideration these aspects, age-appropriated technologies and
models in paediatric drug development should be applied: appropriate cellular
models, juvenile animal model, administration of sub-pharmacologic doses (micro-
dosing) to evaluate PK in a first-in-paediatric study, modelling and simulations and
pharmacogenetics biomarkers.
Juvenile animal models should be used to take into due account the specificities
of the paediatric population as described above and to fill the gap between devel-
opmental and mature toxicity. Indeed, the same drugs can have a different safety
profile in children compared to adults due to many aspects such as body weight,
developmental differences in growth and function of target organs, immune system
maturation and different expression of receptors system. For example, adult models
of epilepsy cannot be simply applied to the study of paediatric epilepsy and key
differences exist in human and rodent brain maturation process [45].
Extrapolation of data from adults or studies using adult animals is not always
adequate to predict these differences in safety profile for paediatric age groups.
For this aim, ‘Guideline on the need for non-clinical testing in juvenile animals of
pharmaceuticals for paediatric indications’ has been adopted in January 2008 by
the EMA. The guidelines recommend the ‘use of juvenile animal models when a
drug safety cannot be appropriately defined in the intended paediatric age group
on the basis of human data or previous animal studies’ and provide recommenda-
tion on the ‘timing and utility of juvenile animal studies in relation to phases of
drug development process’. In particular, the document points out that studies in
juvenile animals should be performed on a case-by-case basis rather than using
standardised study protocols and describes the key aspects to take into consid-
eration in the study’s design: age of the animal and duration of the studies, route
of administration, selection of species, PK and toxicokinetics, dose selection,
endpoint [46]. Juvenile studies are especially recommended when it has been
demonstrated that a medicine causes toxicity in adult at the target organ level and/
or to tissues that undergo significant post-natal development (CNS, immune, or
reproductive systems). As also underlined by Anderson et al., it is important to
conduct the preclinical experiments in the most appropriate species at the most
relevant age on the basis of comparability of the specific organ system develop-
ment in question [47]. And many issues have to be considered in juvenile toxicol-
ogy studies: difficulties in the dose administration due to the small size of the
animals, in blood and tissue sample collection, and in distinguishing direct versus
latent effects [48].
Therefore, proper animal models should be developed. As an example, Lohi
et al. described the zebrafish as a model for paediatric diseases, with particular
140
Challenges and New Frontiers in the Paediatric Drug Discovery and Development
DOI: [Link]
Finally, efficacy and safety of the new medicine should be tested in appropriate
clinical trials. When it comes to the clinical development of a drug, several issues
related to the peculiarity of the paediatric population have to be faced. Conducting
a paediatric clinical trial raises several scientific and operational challenges.
First of all, the low prevalence of many paediatric diseases leads to a limited
number of children affected by each condition. In addition, the ethical issues are
also to be considered to obtain clinical benefits for children assuring the best pos-
sible protection for these vulnerable subjects. Moreover, considering the heteroge-
neous nature of the paediatric population, the population subsets to be included in a
study should be chosen with great attention in order to be sure to consider the most
likely target population for the medicine being tested.
Another issue to be considered in the design of a paediatric clinical trial is the
lack of tools and/or methods for quantitative and qualitative assessment tailored for
the paediatric population and its sub-groups (study endpoints, questionnaires and
scales for the measurement of psychophysical parameters and tools for the assess-
ment of adverse reactions).
The difficulties described above, in testing appropriate drugs in children,
have brought to an increased use of off-label drugs with high risks for adverse
safety events and efficacy failures and to a general knowledge gap in paediatric
research [1].
The US and EU Regulatory agencies foster the drug clinical development
through regulations and incentives and the increasing number of paediatric trials
and specific label changes and dosing recommendations.
Ground-breaking methodologies such as innovative trial design, application of
modelling and simulation and other tools supporting paediatric trials (such as spe-
cific outcomes measures, biomarkers, statistical methods, etc.) can help researchers
to overcome obstacles faced in planning, initiating and conducting a clinical trial
involving children.
141
Drug Discovery and Development - New Advances
For example, to reduce the number of samples required for a study, sparse and
scavenged sampling approach can be used. Sparse sampling uses a lower number of
samples per patient compared with traditional PK sampling methods. Scavenged
sampling consists in the use of residual blood/plasma samples remaining after the
laboratory tests obtained in the course of medical care. These approaches reduce
the risk for the child and eliminate the need for vascular punctures specifically for
the study and, as a consequence, increase the rate of parental consent and the avail-
ability of several samples per infant [56].
Statistical methods, such as the Bayesian design, allow the extrapolation of
results out of fewer children than in the conventional, fixed-number design, also
considering evidences in adults [5].
Modelling and simulation approaches allow to successfully predict the optimal
dosing regimens from the preclinical to the clinical phase [57].
More innovative trial design methods are being developed to overcome the
limits related to small samples and to the acceptability of the trial. These alterna-
tive approaches, limiting the amount of experimentation in children, represent a
promising way of ultimately improving paediatric care [58].
142
Challenges and New Frontiers in the Paediatric Drug Discovery and Development
DOI: [Link]
143
Drug Discovery and Development - New Advances
Most studies on human placental biology have been conducted on tissue obtained
after term delivery, or earlier, often from pathological pregnancies at various stages
of disease, or from ex vivo model system. Less information can be obtained about the
earlier phases of gestation and the normal development and functions of human pla-
centa [67]. Behind these difficulties in obtaining the tissue, the studies on placenta
require high level of expertise. To overcome these limitations, some initiatives have
been undertaken. A 3D in vitro model of human placenta has been developed by a
research group at the University of Vienna. The 3D model shows self-organisation,
self-renewal and constant growth capacity and can be also pharmacologically and
genetically manipulated allowing to study the physiological and pathophysiological
processes of human placenta [68]. Another attempt to develop a model of human
placenta has been carried out by the Huh Lab at the University of Pennsylvania,
which developed the first placenta-on-a-chip to study drug delivery to the placenta
and preterm birth. It consists of a small block of silicone that contains two overlap-
ping layers of microchannels that are lined with trophoblast cells isolated from the
outer surface of the placental barrier and separated by a porous membrane [69].
These advancements will allow a better understanding of the transport processes
through the placenta and a better designing of new obstetrician drug.
144
Challenges and New Frontiers in the Paediatric Drug Discovery and Development
DOI: [Link]
145
Drug Discovery and Development - New Advances
(Innovative Medicines Initiative 2) project [75] and the PedCRIN project [76]. C4C
is a European project aimed to implement an infrastructure of clinical sites organ-
ised as a pan-European network to test medicines through well-organised, moni-
tored and evaluated profit and non-profit paediatric clinical trials. PedCRIN project
is intended to develop tools and actions for paediatric and neonatal trials in order to
better address the real needs and gaps of the paediatric research community.
Behind the initiatives mentioned above, other actions have been taken to foster
the early drug discovery and preclinical development phases. In this field, we can
mention the European Paediatric Translational Research Infrastructures (EPTRI)
project [77], aimed to design a research infrastructure (RI) completely dedicated
to paediatrics to be included in the landscape of the ESFRI RIs. EPTRI aims to be
complementary and fully integrated in the context of the existing RIs providing
services, competences, expertise in the paediatric drug discovery and develop-
ment. EPTRI will provide support to the paediatric research community through
its thematic platform: Human Development and Paediatric Medicines Discovery,
Paediatric Biomarkers and Biosamples, Paediatric Pharmacology, and Paediatric
Medicines Formulations and Medical Devices. Through them, EPTRI will promote
a translational approach from the bedside to the bench side, to make available
more efficacious and safer drugs for children. In the formulation field, it has to be
mentioned that the European Paediatric Formulation Initiative (EuPFI) [78] is a
consortium working in a pre-competitive way on paediatric drug formulations and
aimed to speed up the development of better and safer medicines for children by
identifying issues and challenges in paediatric formulation development. EuPFI has
set up the database Safety and Toxicity of Excipients for Paediatrics (STEP) that
provides updated information on excipients safety and toxicity in children.
To address specifically the rare diseases, the European Joint Programme Rare
Disease (EJP RD) [79], recently founded by the European Commission, brings over
130 institutions from 35 countries to create virtuous circle among research, care and
medical innovation in the rare disease landscape. In particular, the project will improve
the integration, the efficacy, the development and the social impact of research on
rare disease and will implement an efficient model of financial support for all types
of research on RD (fundamental, clinical, epidemiological, social, economic and
health service), providing support to accelerate the exploitation of research results
for the benefit of patients. To more specifically focus on the drug development in rare
diseases, a task force has been created within International Rare Diseases Research
Consortium(IRDiRC) [80], the Orphan Drug Development Guidebook Taskforce,
aimed at providing support to academic and industrial drug developers and describing
the available tools and initiatives specific for rare disease drug development [81].
As described, many initiatives exist as a result of the growing understanding
that children cannot be considered as small adults, but need to be addressed specifi-
cally in the drug development path. But more efforts and the involvement of the
national and international policy bodies are still needed to make the development
of medicines for children a priority.
6. Conclusion
146
Challenges and New Frontiers in the Paediatric Drug Discovery and Development
DOI: [Link]
recognised lack of medicines specifically targeted for them. Moreover, the enor-
mous progresses and advancements reached in the pharmaceutical field have not
been applied to the paediatric population at the same extent of the adults.
The gap in the availability of proper medicines for children can be traced back to
ethical, practical and economic reasons. As discussed in the chapter, the main prac-
tical reasons can be associated with the differences existing in the diseases affecting
children compared to adults, as well as in the different physiology itself of the chil-
dren compared to adults, the low number of patients affected, the need to take into
account different age groups and the need to make available appropriate formula-
tions. Moreover, the ethical concerns make more difficult to obtain the parents’
consent. In addition, the pharmaceutical companies are not interested in this niche
market, since they cannot foresee an adequate economic return. Furthermore, more
challenges have to be faced when considering paediatric rare diseases. Complex
aetiology, small affected population and subsequently small market size, high cost,
and possibly low return on investment led to a large gap between basic research and
patient unmet needs for rare disease drug discovery.
Many initiatives have been taken over the years, also at institutional levels, to
promote a ‘good research’ in the paediatric field, in order to involve children and at
the same time preserve them by unnecessary risks. Only increasing our understand-
ing about human development processes and about how these processes impact
on the onset and progression of diseases will able us to develop specific medicines
targeted for children. The knowledge of these processes will allow us to transfer in
the paediatrics all the advancements and innovative technologies nowadays avail-
able in the adults’ pharmacological research. Thus, more efforts are needed in terms
of capitals, human resources, and technological expertise to speed up both the
preclinical and clinical drug development in children and make available to children
new medicines and appropriate treatments.
Conflict of interest
Author details
© 2019 The Author(s). Licensee IntechOpen. This chapter is distributed under the terms
of the Creative Commons Attribution License ([Link]
by/3.0), which permits unrestricted use, distribution, and reproduction in any medium,
provided the original work is properly cited.
147
Drug Discovery and Development - New Advances
References
[3] Adams DP. The penicillin [12] Kearns GL. Selecting the proper
mystique and the popular press pediatric dose: It is more than size that
(1935-1950). Pharmacy in History. matters. Clinical Pharmacology and
1984;26(3):134-142 Therapeutics. 2015;98(3):238-240
[4] Ward RM et al. The need for [13] Grobner SN et al. The landscape of
pediatric drug development. The genomic alterations across childhood
Journal of Pediatrics. 2018;192:13-21 cancers. Nature. 2018;555(7696):321-327
148
Challenges and New Frontiers in the Paediatric Drug Discovery and Development
DOI: [Link]
149
Drug Discovery and Development - New Advances
150
Challenges and New Frontiers in the Paediatric Drug Discovery and Development
DOI: [Link]
[59] Biomarkers Definitions Working [67] Grigsby PL. Animal models to study
Group. Biomarkers and surrogate placental development and function
endpoints: Preferred definitions throughout normal and dysfunctional
and conceptual framework. Clinical human pregnancy. Seminars in
Pharmacology and Therapeutics. Reproductive Medicine. 2016;34(1):11-16
2001;69(3):89-95
[68] 3d-Model-of-Human-Placenta. 2018.
[60] Bai JP et al. Strategic biomarkers Available from: [Link]
for drug development in treating [Link]/web/en/about-us/news/
rare diseases and diseases in neonates detailsite/2018/news-im-august-2018/3d-
and infants. The AAPS Journal. model-of-human-placenta-developed/
2013;15(2):447-454
[69] Placenta-on-a-Chip. 2018. Available
[61] Goldman J et al. Development of from: [Link]
biomarkers to optimize pediatric patient org/2018/09/placenta-on-a-chip/
management: What makes children
different? Biomarkers in Medicine. [70] Shekhar C. In silico
2011;5(6):781-794 pharmacology: Computer-aided
methods could transform drug
[62] Shores DR, Everett AD. Children as development. Chemistry & Biology.
biomarker orphans: Progress in the field 2008;15(5):413-414
of pediatric biomarkers. The Journal of
Pediatrics. 2018;193:14-20.e31 [71] Johnson TN, Rostami-Hodjegan A,
Tucker GT. Prediction of the clearance
[63] European Medicines Agency. of eleven drugs and associated
Committee for Medicinal Products variability in neonates, infants and
for Human Use (CHMP). Reflection children. Clinical Pharmacokinetics.
Paper: Formulations of Choice for the 2006;45(9):931-956
Paediatric Population. EMEA/CHMP/
PEG/194810/2005. 2005. Available [72] Maharaj AR, Barrett JS, Edginton
from: [Link] AN. A workflow example of PBPK
documents/scientific-guideline/ modeling to support pediatric research
reflection-paper-formulations-choice- and development: Case study with
paediatric-population_en.pdf lorazepam. The AAPS Journal.
2013;15(2):455-464
[64] Tronnes JN, Lupattelli A, Nordeng
H. Safety profile of medication [73] EnprEMA. Available from:
used during pregnancy: Results of [Link]
a multinational European study. partners-networks/networks/european-
Pharmacoepidemiology and Drug network-paediatric-research-european-
Safety. 2017;26(7):802-811 medicines-agency-enpr-ema
151
Drug Discovery and Development - New Advances
152
Edited by Vishwanath Gaitonde,
Partha Karmakar and Ashit Trivedi
ISBN
ISBN 978-1-78985-219-6
978-1-78923-975-1
Published in London, UK
© 2020 IntechOpen
© Rachandewwa / iStock