Complementarity between in silico and biophysical screening approaches in fragment-based lead discovery against the A(2A) adenosine receptor

Hit-to-Lead Optimization of Heterocyclic Carbonyloxycarboximidamides as Selective Antagonists at Human Adenosine A3 Receptor

Antagonism of the human adenosine A3 receptor (hA3R) has potential therapeutic application. Alchemical relative binding free energy calculations of K18 and K32 suggested that the combination of a 3-(2,6-dichlorophenyl)-isoxazolyl group with 2-pyridinyl at the ends of a carbonyloxycarboximidamide group should improve hA3R affinity. Of the 25 new analogues synthesized, 37 and 74 showed improved hA3R affinity compared to K18 (and K32). This was further improved through the addition of a bromine group to the 2-pyridinyl at the 5-position, generating compound 39. Alchemical relative binding free energy calculations, mutagenesis studies and MD simulations supported the compounds' binding pattern while suggesting that the bromine of 39 inserts deep into the hA3R orthosteric pocket, so highlighting the importance of rigidification of the carbonyloxycarboximidamide moiety. MD simulations highlighted the importance of rigidification of the carbonyloxycarboximidamide, while suggesting that the bromine of 39 inserts deep into the hA3R orthosteric pocket, which was supported through mutagenesis studies 39 also selectively antagonized endogenously expressed hA3R in nonsmall cell lung carcinoma cells, while pharmacokinetic studies indicated low toxicity enabling in vivo evaluation. We therefore suggest that 39 has potential for further development as a high-affinity hA3R antagonist.

Alchemical Free Energy Calculations on Membrane-Associated Proteins

Membrane proteins have diverse functions within cells and are well-established drug targets. The advances in membrane protein structural biology have revealed drug and lipid binding sites on membrane proteins, while computational methods such as molecular simulations can resolve the thermodynamic basis of these interactions. Particularly, alchemical free energy calculations have shown promise in the calculation of reliable and reproducible binding free energies of protein-ligand and protein-lipid complexes in membrane-associated systems. In this review, we present an overview of representative alchemical free energy studies on G-protein-coupled receptors, ion channels, transporters as well as protein-lipid interactions, with emphasis on best practices and critical aspects of running these simulations. Additionally, we analyze challenges and successes when running alchemical free energy calculations on membrane-associated proteins. Finally, we highlight the value of alchemical free energy calculations calculations in drug discovery and their applicability in the pharmaceutical industry.

Targeting RNA Structure to Inhibit Editing in Trypanosomes

Mitochondrial RNA editing in trypanosomes represents an attractive target for developing safer and more efficient drugs for treating infections with trypanosomes because this RNA editing pathway is not found in humans. Other workers have targeted several enzymes in this editing system, but not the RNA. Here, we target a universal domain of the RNA editing substrate, which is the U-helix formed between the oligo-U tail of the guide RNA and the target mRNA. We selected a part of the U-helix that is rich in G-U wobble base pairs as the target site for the virtual screening of 262,000 compounds. After chemoinformatic filtering of the top 5000 leads, we subjected 50 representative complexes to 50 nanoseconds of molecular dynamics simulations. We identified 15 compounds that retained stable interactions in the deep groove of the U-helix. The microscale thermophoresis binding experiments on these five compounds show low-micromolar to nanomolar binding affinities. The UV melting studies show an increase in the melting temperatures of the U-helix upon binding by each compound. These five compounds can serve as leads for drug development and as research tools to probe the role of the RNA structure in trypanosomal RNA editing.

Evaluating the use of absolute binding free energy in the fragment optimisation process

Key to the fragment optimisation process within drug design is the need to accurately capture the changes in affinity that are associated with a given set of chemical modifications. Due to the weakly binding nature of fragments, this has proven to be a challenging task, despite recent advancements in leveraging experimental and computational methods. In this work, we evaluate the use of Absolute Binding Free Energy (ABFE) calculations in guiding fragment optimisation decisions, retrospectively calculating binding free energies for 59 ligands across 4 fragment elaboration campaigns. We first demonstrate that ABFEs can be used to accurately rank fragment-sized binders with an overall Spearman's r of 0.89 and a Kendall τ of 0.67, although often deviating from experiment in absolute free energy values with an RMSE of 2.75 kcal/mol. We then also show that in several cases, retrospective fragment optimisation decisions can be supported by the ABFE calculations. Comparing against cheaper endpoint methods, namely Nwat-MM/GBSA, we find that ABFEs offer better ranking power and correlation metrics. Our results indicate that ABFE calculations can usefully guide fragment elaborations to maximise affinity.

Adenosine A2A receptor antagonists: from caffeine to selective non-xanthines

A long evolution of knowledge of the psychostimulant caffeine led in the 1960s to another purine natural product, adenosine and its A2A receptor. Adenosine is a short-lived autocrine/paracrine mediator that acts pharmacologically at four different adenosine receptors in a manner opposite to the pan-antagonist caffeine and serves as an endogenous allostatic regulator. Although detrimental in the developing brain, caffeine appears to be cerebroprotective in aging. Moderate caffeine consumption in adults, except in pregnancy, may also provide benefit in pain, diabetes, and kidney and liver disorders. Inhibition of A2A receptors is one of caffeine's principal effects and we now understand this interaction at the atomic level. The A2A receptor has become a prototypical example of utilizing high-resolution structures of GPCRs for the rational design of chemically diverse drug molecules. The previous focus on discovery of selective A2A receptor antagonists for neurodegenerative diseases has expanded to include immunotherapy for cancer, and clinical trials have ensued. LINKED ARTICLES: This article is part of a themed issue on Structure Guided Pharmacology of Membrane Proteins (BJP 75th Anniversary). To view the other articles in this section visit http://onlinelibrary.wiley.com/doi/10.1111/bph.v179.14/issuetoc.

Inverse Mixed-Solvent Molecular Dynamics for Visualization of the Residue Interaction Profile of Molecular Probes

To ensure efficiency in discovery and development, the application of computational technology is essential. Although virtual screening techniques are widely applied in the early stages of drug discovery research, the computational methods used in lead optimization to improve activity and reduce the toxicity of compounds are still evolving. In this study, we propose a method to construct the residue interaction profile of the chemical structure used in the lead optimization by performing “inverse” mixed-solvent molecular dynamics (MSMD) simulation. Contrary to constructing a protein-based, atom interaction profile, we constructed a probe-based, protein residue interaction profile using MSMD trajectories. It provides us the profile of the preferred protein environments of probes without co-crystallized structures. We assessed the method using three probes: benzamidine, catechol, and benzene. As a result, the residue interaction profile of each probe obtained by MSMD was a reasonable physicochemical description of the general non-covalent interaction. Moreover, comparison with the X-ray structure containing each probe as a ligand shows that the map of the interaction profile matches the arrangement of amino acid residues in the X-ray structure.

Structure-Based Virtual Screening for Ligands of G Protein-Coupled Receptors: What Can Molecular Docking Do for You?

G protein-coupled receptors (GPCRs) constitute the largest family of membrane proteins in the human genome and are important therapeutic targets. During the last decade, the number of atomic-resolution structures of GPCRs has increased rapidly, providing insights into drug binding at the molecular level. These breakthroughs have created excitement regarding the potential of using structural information in ligand design and initiated a new era of rational drug discovery for GPCRs. The molecular docking method is now widely applied to model the three-dimensional structures of GPCR-ligand complexes and screen for chemical probes in large compound libraries. In this review article, we first summarize the current structural coverage of the GPCR superfamily and the understanding of receptor-ligand interactions at atomic resolution. We then present the general workflow of structure-based virtual screening and strategies to discover GPCR ligands in chemical libraries. We assess the state of the art of this research field by summarizing prospective applications of virtual screening based on experimental structures. Strategies to identify compounds with specific efficacy and selectivity profiles are discussed, illustrating the opportunities and limitations of the molecular docking method. Our overview shows that structure-based virtual screening can discover novel leads and will be essential in pursuing the next generation of GPCR drugs. SIGNIFICANCE STATEMENT: Extraordinary advances in the structural biology of G protein-coupled receptors have revealed the molecular details of ligand recognition by this large family of therapeutic targets, providing novel avenues for rational drug design. Structure-based docking is an efficient computational approach to identify novel chemical probes from large compound libraries, which has the potential to accelerate the development of drug candidates.

Engineering of Challenging G Protein-Coupled Receptors for Structure Determination and Biophysical Studies

Membrane proteins such as G protein-coupled receptors (GPCRs) exert fundamental biological functions and are involved in a multitude of physiological responses, making these receptors ideal drug targets. Drug discovery programs targeting GPCRs have been greatly facilitated by the emergence of high-resolution structures and the resulting opportunities to identify new chemical entities through structure-based drug design. To enable the determination of high-resolution structures of GPCRs, most receptors have to be engineered to overcome intrinsic hurdles such as their poor stability and low expression levels. In recent years, multiple engineering approaches have been developed to specifically address the technical difficulties of working with GPCRs, which are now beginning to make more challenging receptors accessible to detailed studies. Importantly, successfully engineered GPCRs are not only valuable in X-ray crystallography, but further enable biophysical studies with nuclear magnetic resonance spectroscopy, surface plasmon resonance, native mass spectrometry, and fluorescence anisotropy measurements, all of which are important for the detailed mechanistic understanding, which is the prerequisite for successful drug design. Here, we summarize engineering strategies based on directed evolution to reduce workload and enable biophysical experiments of particularly challenging GPCRs.

Application of In Silico and HTS Approaches to Identify Nuclear Import Inhibitors for Venezuelan Equine Encephalitis Virus Capsid Protein: A Case Study

The development of new drugs is costly and time-consuming, with estimates of over $US1 billion and 15 years for a product to reach the market. As understanding of the molecular basis of disease improves, various approaches have been used to target specific molecular interactions in the search for effective drugs. These include high-throughput screening (HTS) for novel drug identification and computer-aided drug design (CADD) to assess the properties of putative drugs before experimental work begins. We have applied conventional HTS and CADD approaches to the problem of identifying antiviral compounds to limit infection by Venezuelan equine encephalitis virus (VEEV). Nuclear targeting of the VEEV capsid (CP) protein through interaction with the host nuclear import machinery has been shown to be essential for viral pathogenicity, with viruses incapable of this interaction being greatly attenuated. Our previous conventional HTS and in silico structure-based drug design (SBDD) screens were successful in identifying novel inhibitors of CP interaction with the host nuclear import machinery, thus providing a unique opportunity to assess the relative value of the two screening approaches directly. This focused review compares and contrasts the two screening approaches, together with the properties of the inhibitors identified, as a case study for parallel use of the two approaches to identify antivirals. The utility of SBDD screens, especially when used in parallel with traditional HTS, in identifying agents of interest to target the host-pathogen interface is highlighted.

NMR as a "Gold Standard" Method in Drug Design and Discovery

Studying disease models at the molecular level is vital for drug development in order to improve treatment and prevent a wide range of human pathologies. Microbial infections are still a major challenge because pathogens rapidly and continually evolve developing drug resistance. Cancer cells also change genetically, and current therapeutic techniques may be (or may become) ineffective in many cases. The pathology of many neurological diseases remains an enigma, and the exact etiology and underlying mechanisms are still largely unknown. Viral infections spread and develop much more quickly than does the corresponding research needed to prevent and combat these infections; the present and most relevant outbreak of SARS-CoV-2, which originated in Wuhan, China, illustrates the critical and immediate need to improve drug design and development techniques. Modern day drug discovery is a time-consuming, expensive process. Each new drug takes in excess of 10 years to develop and costs on average more than a billion US dollars. This demonstrates the need of a complete redesign or novel strategies. Nuclear Magnetic Resonance (NMR) has played a critical role in drug discovery ever since its introduction several decades ago. In just three decades, NMR has become a "gold standard" platform technology in medical and pharmacology studies. In this review, we present the major applications of NMR spectroscopy in medical drug discovery and development. The basic concepts, theories, and applications of the most commonly used NMR techniques are presented. We also summarize the advantages and limitations of the primary NMR methods in drug development.

Combined Protein- and Ligand-Observed NMR Workflow to Screen Fragment Cocktails against Multiple Proteins: A Case Study Using Bromodomains

As fragment-based drug discovery has become mainstream, there has been an increase in various screening methodologies. Protein-observed ¹⁹F (PrOF) NMR and ¹H CPMG NMR are two fragment screening assays that have complementary advantages. Here, we sought to combine these two NMR-based assays into a new screening workflow. This combination of protein- and ligand-observed experiments allows for a time- and resource-efficient multiplexed screen of mixtures of fragments and proteins. PrOF NMR is first used to screen mixtures against two proteins. Hit mixtures for each protein are identified then deconvoluted using ¹H CPMG NMR. We demonstrate the benefit of this fragment screening method by conducting the first reported fragment screens against the bromodomains of BPTF and Plasmodium falciparum (Pf) GCN5 using 467 3D-enriched fragments. The hit rates were 6%, 5% and 4% for fragments binding BPTF, PfGCN5, and fragments binding both proteins, respectively. Select hits were characterized, revealing a broad range of affinities from low µM to mM dissociation constants. Follow-up experiments supported a low-affinity second binding site on PfGCN5. This approach can be used to bias fragment screens towards more selective hits at the onset of inhibitor development in a resource- and time-efficient manner.

In Silico Drug Design for Purinergic GPCRs: Overview on Molecular Dynamics Applied to Adenosine and P2Y Receptors

Molecular modeling has contributed to drug discovery for purinergic GPCRs, including adenosine receptors (ARs) and P2Y receptors (P2YRs). Experimental structures and homology modeling have proven to be useful in understanding and predicting structure activity relationships (SAR) of agonists and antagonists. This review provides an excursus on molecular dynamics (MD) simulations applied to ARs and P2YRs. The binding modes of newly synthesized A1AR- and A3AR-selective nucleoside derivatives, potentially of use against depression and inflammation, respectively, have been predicted to recapitulate their SAR and the species dependence of A3AR affinity. P2Y12R and P2Y1R crystallographic structures, respectively, have provided a detailed understanding of the recognition of anti-inflammatory P2Y14R antagonists and a large group of allosteric and orthosteric antagonists of P2Y1R, an antithrombotic and neuroprotective target. MD of A2AAR (an anticancer and neuroprotective target), A3AR, and P2Y1R has identified microswitches that are putatively involved in receptor activation. The approach pathways of different ligands toward A2AAR and P2Y1R binding sites have also been explored. A1AR, A2AAR, and A3AR were utilizes to study allosteric phenomena, but locating the binding site of structurally diverse allosteric modulators, such as an A3AR enhancer LUF6000, is challenging. Ligand residence time, a predictor of in vivo efficacy, and the structural role of water were investigated through A2AAR MD simulations. Thus, new MD and other modeling algorithms have contributed to purinergic GPCR drug discovery.

Small Polar Hits against S. aureus: Screening, Initial Hit Optimization, and Metabolomic Studies

The global prevalence of antibacterial resistance requires new antibacterial drugs with novel chemical scaffolds and modes of action. It is also vital to design compounds with optimal physicochemical properties to permeate the bacterial cell envelope. We described an approach of combining and integrating whole cell screening and metabolomics into early antibacterial drug discovery using a library of small polar compounds. Whole cell screening of a diverse library of small polar compounds against Staphylococcus aureus gave compound 2. Hit expansion was carried out to determine structure-activity relationships. A selection of compounds from this series, together with other screened active compounds, was subjected to an initial metabolomics study to provide a metabolic fingerprint of the mode of action. It was found that compound 2 and its analogues have a different mode of action from some of the known antibacterial compounds tested. This early study highlighted the potential of whole cell screening and metabolomics in early antibacterial drug discovery. Future works will require improving potency and performing orthogonal studies to confirm the modes of action.

Binding mode information improves fragment docking

Docking is commonly used in drug discovery to predict how ligand binds to protein target. Best programs are generally able to generate a correct solution, yet often fail to identify it. In the case of drug-like molecules, the correct and incorrect poses can be sorted by similarity to the crystallographic structure of the protein in complex with reference ligands. Fragments are particularly sensitive to scoring problems because they are weak ligands which form few interactions with protein. In the present study, we assessed the utility of binding mode information in fragment pose prediction. We compared three approaches: interaction fingerprints, 3D-matching of interaction patterns and 3D-matching of shapes. We prepared a test set composed of high-quality structures of the Protein Data Bank. We generated and evaluated the docking poses of 586 fragment/protein complexes. We observed that the best approach is twice as accurate as the native scoring function, and that post-processing is less effective for smaller fragments. Interestingly, fragments and drug-like molecules both proved to be useful references. In the discussion, we suggest the best conditions for a successful pose prediction with the three approaches.

Unpacking the Pathogen Box-An Open Source Tool for Fighting Neglected Tropical Disease

The Pathogen Box is a 400-strong collection of drug-like compounds, selected for their potential against several of the world's most important neglected tropical diseases, including trypanosomiasis, leishmaniasis, cryptosporidiosis, toxoplasmosis, filariasis, schistosomiasis, dengue virus and trichuriasis, in addition to malaria and tuberculosis. This library represents an ensemble of numerous successful drug discovery programmes from around the globe, aimed at providing a powerful resource to stimulate open source drug discovery for diseases threatening the most vulnerable communities in the world. This review seeks to provide an in-depth analysis of the literature pertaining to the compounds in the Pathogen Box, including structure-activity relationship highlights, mechanisms of action, related compounds with reported activity against different diseases, and, where appropriate, discussion on the known and putative targets of compounds, thereby providing context and increasing the accessibility of the Pathogen Box to the drug discovery community.

Solvents to Fragments to Drugs: MD Applications in Drug Design

Simulations of molecular dynamics (MD) are playing an increasingly important role in structure-based drug discovery (SBDD). Here we review the use of MD for proteins in aqueous solvation, organic/aqueous mixed solvents (MDmix) and with small ligands, to the classic SBDD problems: Binding mode and binding free energy predictions. The simulation of proteins in their condensed state reveals solvent structures and preferential interaction sites (hot spots) on the protein surface. The information provided by water and its cosolvents can be used very effectively to understand protein ligand recognition and to improve the predictive capability of well-established methods such as molecular docking. The application of MD simulations to the study of the association of proteins with drug-like compounds is currently only possible for specific cases, as it remains computationally very expensive and labor intensive. MDmix simulations on the other hand, can be used systematically to address some of the common tasks in SBDD. With the advent of new tools and faster computers we expect to see an increase in the application of mixed solvent MD simulations to a plethora of protein targets to identify new drug candidates.

Computational insights into factor affecting the potency of diaryl sulfone analogs as Escherichia coli dihydropteroate synthase inhibitors

Dihydropteroate synthase (DHPS) is an alluring target for designing novel drug candidates to prevent infections caused by pathogenic Escherichia coli strains. Diaryl Sulfone (SO) compounds are found to inhibit DHPS competitively with respect to the substrate pABA (p-aminobenzoate). The extra aromatic ring of diaryl sulfone compounds found to stabilize them in highly flexible pABA binding loops. In this present study, a statistically significant 3D-QSAR model was developed using a data set of diaryl sulfone compounds. The favourable and unfavourable contributions of substitutions in sulfone compounds were illustrated by contour plot obtained from the developed 3D-QSAR model. Molecular docking calculations were performed to investigate the putative binding mode of diaryl sulfone compounds at the catalytic pocket. DFT calculations were carried out using SCF approach, B3LYP- 6-31 G (d) basis set to compute the HOMO, LUMO energies and their respective location at pABA binding pocket. Further, the developed model was validated by FEP (Free Energy Perturbation) calculations. The calculated relative free energy of binding between the highly potent and less potent sulfone compound was found to be -3.78 kcal/ mol which is comparable to the experimental value of -5.85 kcal/mol. A 10 ns molecular dynamics simulation of inhibitor and DHPS confirmed its stability at pABA catalytic site. Outcomes of the present work provide deeper insight in designing novel drug candidates for pathogenic Escherichia coli strains.

In Silico Studies in Drug Research Against Neurodegenerative Diseases

Background

Neurodegenerative diseases such as Alzheimer's disease (AD), amyotrophic lateral sclerosis, Parkinson's disease (PD), spinal cerebellar ataxias, and spinal and bulbar muscular atrophy are described by slow and selective degeneration of neurons and axons in the central nervous system (CNS) and constitute one of the major challenges of modern medicine. Computer-aided or in silico drug design methods have matured into powerful tools for reducing the number of ligands that should be screened in experimental assays.

Methods

In the present review, the authors provide a basic background about neurodegenerative diseases and in silico techniques in the drug research. Furthermore, they review the various in silico studies reported against various targets in neurodegenerative diseases, including homology modeling, molecular docking, virtual high-throughput screening, quantitative structure activity relationship (QSAR), hologram quantitative structure activity relationship (HQSAR), 3D pharmacophore mapping, proteochemometrics modeling (PCM), fingerprints, fragment-based drug discovery, Monte Carlo simulation, molecular dynamic (MD) simulation, quantum-mechanical methods for drug design, support vector machines, and machine learning approaches.

Results

Detailed analysis of the recently reported case studies revealed that the majority of them use a sequential combination of ligand and structure-based virtual screening techniques, with particular focus on pharmacophore models and the docking approach.

Conclusion

Neurodegenerative diseases have a multifactorial pathoetiological origin, so scientists have become persuaded that a multi-target therapeutic strategy aimed at the simultaneous targeting of multiple proteins (and therefore etiologies) involved in the development of a disease is recommended in future.

Recent Trends and Applications of Molecular Modeling in GPCR⁻Ligand Recognition and Structure-Based Drug Design

G protein-coupled receptors represent the largest family of human membrane proteins and are modulated by a variety of drugs and endogenous ligands. Molecular modeling techniques, especially enhanced sampling methods, have provided significant insight into the mechanism of GPCR–ligand recognition. Notably, the crucial role of the membrane in the ligand-receptor association process has earned much attention. Additionally, docking, together with more accurate free energy calculation methods, is playing an important role in the design of novel compounds targeting GPCRs. Here, we summarize the recent progress in the computational studies focusing on the above issues. In the future, with continuous improvement in both computational hardware and algorithms, molecular modeling would serve as an indispensable tool in a wider scope of the research concerning GPCR–ligand recognition as well as drug design targeting GPCRs.

Discovery of Novel Adenosine Receptor Antagonists through a Combined Structure- and Ligand-Based Approach Followed by Molecular Dynamics Investigation of Ligand Binding Mode

An intense effort is made by pharmaceutical and academic research laboratories to identify and develop selective antagonists for each adenosine receptor (AR) subtype as potential clinical candidates for "soft" treatment of various diseases. Crystal structures of subtypes A_2A and A₁ARs offer exciting opportunities for structure-based drug design. In the first part of the present work, Maybridge HitFinder library of 14400 compounds was utilized to apply a combination of structure-based against the crystal structure of A_2AAR and ligand-based methodologies. The docking poses were rescored by CHARMM energy minimization and calculation of the desolvation energy using Poisson-Boltzmann equation electrostatics. Out of the eight selected and tested compounds, five were found positive hits (63% success). Although the project was initially focused on targeting A_2AAR, the identified antagonists exhibited low micromolar or micromolar affinity against A_2A/A₃, ARs, or A₃AR, respectively. Based on these results, 19 compounds characterized by novel chemotypes were purchased and tested. Sixteen of them were identified as AR antagonists with affinity toward combinations of the AR family isoforms (A_2A/A₃, A₁/A₃, A₁/A_2A/A₃, and A₃). The second part of this work involves the performance of hundreds of molecular dynamics (MD) simulations of complexes between the ARs and a total of 27 ligands to resolve the binding interactions of the active compounds, which were not achieved by docking calculations alone. This computational work allowed the prediction of stable and unstable complexes which agree with the experimental results of potent and inactive compounds, respectively. Of particular interest is that the 2-amino-thiophene-3-carboxamides, 3-acylamino-5-aryl-thiophene-2-carboxamides, and carbonyloxycarboximidamide derivatives were found to be selective and possess a micromolar to low micromolar affinity for the A₃ receptor.

Exploring G Protein-Coupled Receptors (GPCRs) Ligand Space via Cheminformatics Approaches: Impact on Rational Drug Design

The primary goal of rational drug discovery is the identification of selective ligands which act on single or multiple drug targets to achieve the desired clinical outcome through the exploration of total chemical space. To identify such desired compounds, computational approaches are necessary in predicting their drug-like properties. G Protein-Coupled Receptors (GPCRs) represent one of the largest and most important integral membrane protein families. These receptors serve as increasingly attractive drug targets due to their relevance in the treatment of various diseases, such as inflammatory disorders, metabolic imbalances, cardiac disorders, cancer, monogenic disorders, etc. In the last decade, multitudes of three-dimensional (3D) structures were solved for diverse GPCRs, thus referring to this period as the "golden age for GPCR structural biology." Moreover, accumulation of data about the chemical properties of GPCR ligands has garnered much interest toward the exploration of GPCR chemical space. Due to the steady increase in the structural, ligand, and functional data of GPCRs, several cheminformatics approaches have been implemented in its drug discovery pipeline. In this review, we mainly focus on the cheminformatics-based paradigms in GPCR drug discovery. We provide a comprehensive view on the ligand- and structure-based cheminformatics approaches which are best illustrated via GPCR case studies. Furthermore, an appropriate combination of ligand-based knowledge with structure-based ones, i.e., integrated approach, which is emerging as a promising strategy for cheminformatics-based GPCR drug design is also discussed.

Fragment-Based Discovery and Optimization of Enzyme Inhibitors by Docking of Commercial Chemical Space

Fragment-based lead discovery has emerged as a leading drug development strategy for novel therapeutic targets. Although fragment-based drug discovery benefits immensely from access to atomic-resolution information, structure-based virtual screening has rarely been used to drive fragment discovery and optimization. Here, molecular docking of 0.3 million fragments to a crystal structure of cancer target MTH1 was performed. Twenty-two predicted fragment ligands, for which analogs could be acquired commercially, were experimentally evaluated. Five fragments inhibited MTH1 with IC₅₀ values ranging from 6 to 79 μM. Structure-based optimization guided by predicted binding modes and analogs from commercial chemical libraries yielded nanomolar inhibitors. Subsequently solved crystal structures confirmed binding modes predicted by docking for three scaffolds. Structure-guided exploration of commercial chemical space using molecular docking gives access to fragment libraries that are several orders of magnitude larger than those screened experimentally and can enable efficient optimization of hits to potent leads.

Fragment optimization for GPCRs by molecular dynamics free energy calculations: Probing druggable subpockets of the A 2A adenosine receptor binding site

Fragment-based lead discovery is becoming an increasingly popular strategy for drug discovery. Fragment screening identifies weakly binding compounds that require optimization to become high-affinity leads. As design of leads from fragments is challenging, reliable computational methods to guide optimization would be invaluable. We evaluated using molecular dynamics simulations and the free energy perturbation method (MD/FEP) in fragment optimization for the A_2A adenosine receptor, a pharmaceutically relevant G protein-coupled receptor. Optimization of fragments exploring two binding site subpockets was probed by calculating relative binding affinities for 23 adenine derivatives, resulting in strong agreement with experimental data (R² = 0.78). The predictive power of MD/FEP was significantly better than that of an empirical scoring function. We also demonstrated the potential of the MD/FEP to assess multiple binding modes and to tailor the thermodynamic profile of ligands during optimization. Finally, MD/FEP was applied prospectively to optimize three nonpurine fragments, and predictions for 12 compounds were evaluated experimentally. The direction of the change in binding affinity was correctly predicted in a majority of the cases, and agreement with experiment could be improved with rigorous parameter derivation. The results suggest that MD/FEP will become a powerful tool in structure-driven optimization of fragments to lead candidates.

Ligand Discovery for a Peptide-Binding GPCR by Structure-Based Screening of Fragment- and Lead-Like Chemical Libraries

Peptide-recognizing G protein-coupled receptors (GPCRs) are promising therapeutic targets but often resist drug discovery efforts. Determination of crystal structures for peptide-binding GPCRs has provided opportunities to explore structure-based methods in lead development. Molecular docking screens of two chemical libraries, containing either fragment- or lead-like compounds, against a neurotensin receptor 1 crystal structure allowed for a comparison between different drug development strategies for peptide-binding GPCRs. A total of 2.3 million molecules were screened computationally, and 25 fragments and 27 leads that were top-ranked in each library were selected for experimental evaluation. Of these, eight fragments and five leads were confirmed as ligands by surface plasmon resonance. The hit rate for the fragment screen (32%) was thus higher than for the lead-like library (19%), but the affinities of the fragments were ∼100-fold lower. Both screens returned unique scaffolds and demonstrated that a crystal structure of a stabilized peptide-binding GPCR can guide the discovery of small-molecule agonists. The complementary advantages of exploring fragment- and lead-like chemical space suggest that these strategies should be applied synergistically in structure-based screens against challenging GPCR targets.

Recent Advances and Applications of Molecular Docking to G Protein-Coupled Receptors

The growing number of studies on G protein-coupled receptors (GPCRs) family are a source of noticeable improvement in our understanding of the functioning of these proteins. GPCRs are responsible for a vast part of signaling in vertebrates and, as such, invariably remain in the spotlight of medicinal chemistry. A deeper insight into the underlying mechanisms of interesting phenomena observed in GPCRs, such as biased signaling or allosteric modulation, can be gained with experimental and computational studies. The latter play an important role in this process, since they allow for observations on scales inaccessible for most other methods. One of the key steps in such studies is proper computational reconstruction of actual ligand-receptor or protein-protein interactions, a process called molecular docking. A number of improvements and innovative applications of this method were documented recently. In this review, we focus particularly on innovations in docking to GPCRs.

Structure-Based Screening of Uncharted Chemical Space for Atypical Adenosine Receptor Agonists

Small molecule screening libraries cover only a small fraction of the astronomical number of possible drug-like compounds, limiting the success of ligand discovery efforts. Computational screening of virtual libraries representing unexplored chemical space could potentially bridge this gap. Drug development for adenosine receptors (ARs) as targets for inflammation and cardiovascular diseases has been hampered by the paucity of agonist scaffolds. To identify novel AR agonists, a virtual library of synthetically tractable nucleosides with alternative bases was generated and structure-based virtual screening guided selection of compounds for synthesis. Pharmacological assays were carried out at three AR subtypes for 13 ribosides. Nine compounds displayed significant activity at the ARs, and several of these represented atypical agonist scaffolds. The discovered ligands also provided insights into receptor activation and revealed unknown interactions of endogenous and clinical compounds with the ARs. The results demonstrate that virtual compound databases provide access to bioactive matter from regions of chemical space that are sparsely populated in commercial libraries, an approach transferrable to numerous drug targets.

In search of novel ligands using a structure-based approach: a case study on the adenosine A2A receptor

In this work, we present a case study to explore the challenges associated with finding novel molecules for a receptor that has been studied in depth and has a wealth of chemical information available. Specifically, we apply a previously described protocol that incorporates explicit water molecules in the ligand binding site to prospectively screen over 2.5 million drug-like and lead-like compounds from the commercially available eMolecules database in search of novel binders to the adenosine A_2A receptor (A_2AAR). A total of seventy-one compounds were selected for purchase and biochemical assaying based on high ligand efficiency and high novelty (Tanimoto coefficient ≤0.25 to any A_2AAR tested compound). These molecules were then tested for their affinity to the adenosine A_2A receptor in a radioligand binding assay. We identified two hits that fulfilled the criterion of ~50 % radioligand displacement at a concentration of 10 μM. Next we selected an additional eight novel molecules that were predicted to make a bidentate interaction with Asn253^6.55, a key interacting residue in the binding pocket of the A_2AAR. None of these eight molecules were found to be active. Based on these results we discuss the advantages of structure-based methods and the challenges associated with finding chemically novel molecules for well-explored targets.

Predicting Binding Affinities for GPCR Ligands Using Free-Energy Perturbation

The rapid growth of structural information for G-protein-coupled receptors (GPCRs) has led to a greater understanding of their structure, function, selectivity, and ligand binding. Although novel ligands have been identified using methods such as virtual screening, computationally driven lead optimization has been possible only in isolated cases because of challenges associated with predicting binding free energies for related compounds. Here, we provide a systematic characterization of the performance of free-energy perturbation (FEP) calculations to predict relative binding free energies of congeneric ligands binding to GPCR targets using a consistent protocol and no adjustable parameters. Using the FEP+ package, first we validated the protocol, which includes a full lipid bilayer and explicit solvent, by predicting the binding affinity for a total of 45 different ligands across four different GPCRs (adenosine A_2AAR, β₁ adrenergic, CXCR4 chemokine, and δ opioid receptors). Comparison with experimental binding affinity measurements revealed a highly predictive ranking correlation (average spearman ρ = 0.55) and low root-mean-square error (0.80 kcal/mol). Next, we applied FEP+ in a prospective project, where we predicted the affinity of novel, potent adenosine A_2A receptor (A_2AR) antagonists. Four novel compounds were synthesized and tested in a radioligand displacement assay, yielding affinity values in the nanomolar range. The affinity of two out of the four novel ligands (plus three previously reported compounds) was correctly predicted (within 1 kcal/mol), including one compound with approximately a tenfold increase in affinity compared to the starting compound. Detailed analyses of the simulations underlying the predictions provided insights into the structural basis for the two cases where the affinity was overpredicted. Taken together, these results establish a protocol for systematically applying FEP+ to GPCRs and provide guidelines for identifying potent molecules in drug discovery lead optimization projects.

Identification of Novel Smoothened Ligands Using Structure-Based Docking

The seven transmembrane protein Smoothened is required for Hedgehog signaling during embryonic development and adult tissue homeostasis. Inappropriate activation of the Hedgehog signalling pathway leads to cancers such as basal cell carcinoma and medulloblastoma, and Smoothened inhibitors are now available clinically to treat these diseases. However, resistance to these inhibitors rapidly develops thereby limiting their efficacy. The determination of Smoothened crystal structures enables structure-based discovery of new ligands with new chemotypes that will be critical to combat resistance. In this study, we docked 3.2 million available, lead-like molecules against Smoothened, looking for those with high physical complementarity to its structure; this represents the first such campaign against the class Frizzled G-protein coupled receptor family. Twenty-one high-ranking compounds were selected for experimental testing, and four, representing three different chemotypes, were identified to antagonize Smoothened with IC₅₀ values better than 50 μM. A screen for analogs revealed another six molecules, with IC₅₀ values in the low micromolar range. Importantly, one of the most active of the new antagonists continued to be efficacious at the D473H mutant of Smoothened, which confers clinical resistance to the antagonist vismodegib in cancer treatment.

Fragment Screening of Human Aquaporin 1

Aquaporins (AQPs) are membrane proteins that enable water transport across cellular plasma membranes in response to osmotic gradients. Phenotypic analyses have revealed important physiological roles for AQPs, and the potential for AQP water channel modulators in various disease states has been proposed. For example, AQP1 is overexpressed in tumor microvessels, and this correlates with higher metastatic potential and aggressiveness of the malignancy. Chemical modulators would help in identifying the precise contribution of water channel activity in these disease states. These inhibitors would also be important therapeutically, e.g., in anti-cancer treatment. This perceived importance contrasts with the lack of success of high-throughput screens (HTS) to identify effective and specific inhibitors of aquaporins. In this paper, we have screened a library of 1500 "fragments", i.e., smaller than molecules used in HTS, against human aquaporin (hAQP1) using a thermal shift assay and surface plasmon resonance. Although these fragments may not inhibit their protein target, they bound to and stabilized hAQP1 (sub mM binding affinities (KD), with an temperature of aggregation shift ΔTagg of +4 to +50 °C) in a concentration-dependent fashion. Chemically expanded versions of these fragments should follow the determination of their binding site on the aquaporin surface.

Fragment-Based Discovery of Subtype-Selective Adenosine Receptor Ligands from Homology Models

Fragment-based lead discovery (FBLD) holds great promise for drug discovery, but applications to G protein-coupled receptors (GPCRs) have been limited by a lack of sensitive screening techniques and scarce structural information. If virtual screening against homology models of GPCRs could be used to identify fragment ligands, FBLD could be extended to numerous important drug targets and contribute to efficient lead generation. Access to models of multiple receptors may further enable the discovery of fragments that bind specifically to the desired target. To investigate these questions, we used molecular docking to screen >500 000 fragments against homology models of the A3 and A1 adenosine receptors (ARs) with the goal to discover A3AR-selective ligands. Twenty-one fragments with predicted A3AR-specific binding were evaluated in live-cell fluorescence-based assays; of eight verified ligands, six displayed A3/A1 selectivity, and three of these had high affinities ranging from 0.1 to 1.3 μM. Subsequently, structure-guided fragment-to-lead optimization led to the identification of a >100-fold-selective antagonist with nanomolar affinity from commercial libraries. These results highlight that molecular docking screening can guide fragment-based discovery of selective ligands even if the structures of both the target and antitarget receptors are unknown. The same approach can be readily extended to a large number of pharmaceutically important targets.

A Role for Fragment-Based Drug Design in Developing Novel Lead Compounds for Central Nervous System Targets

Hundreds of millions of U.S. dollars are invested in the research and development of a single drug. Lead compound development is an area ripe for new design strategies. Therapeutic lead candidates have been traditionally found using high-throughput in vitro pharmacological screening, a costly method for assaying thousands of compounds. This approach has recently been augmented by virtual screening (VS), which employs computer models of the target protein to narrow the search for possible leads. A variant of VS is fragment-based drug design (FBDD), an emerging in silico lead discovery method that introduces low-molecular weight fragments, rather than intact compounds, into the binding pocket of the receptor model. These fragments serve as starting points for “growing” the lead candidate. Current efforts in virtual FBDD within central nervous system (CNS) targets are reviewed, as is a recent rule-based optimization strategy in which new molecules are generated within a 3D receptor-binding pocket using the fragment as a scaffold. This process not only places special emphasis on creating synthesizable molecules but also exposes computational questions worth addressing. Fragment-based methods provide a viable, relatively low-cost alternative for therapeutic lead discovery and optimization that can be applied to CNS targets to augment current design strategies.

Discovery of novel inhibitors of human S-adenosylmethionine decarboxylase based on in silico high-throughput screening and a non-radioactive enzymatic assay

Natural polyamines are small polycationic molecules essential for cell growth and development, and elevated level of polyamines is positively correlated with various cancers. As a rate-limiting enzyme of the polyamine biosynthetic pathway, S-adenosylmethionine decarboxylase (AdoMetDC) has been an attractive drug target. In this report, we present the discovery of novel human AdoMetDC (hAdoMetDC) inhibitors by coupling computational and experimental tools. We constructed a reasonable computational structure model of hAdoMetDC that is compatible with general protocols for high-throughput drug screening, and used this model in in silico screening of hAdoMetDC inhibitors against a large compound library using a battery of computational tools. We also established and validated a simple, economic, and non-radioactive enzymatic assay, which can be adapted for experimental high-throughput screening of hAdoMetDC inhibitors. Finally, we obtained an hAdoMetDC inhibitor lead with a novel scaffold. This study provides both new tools and a new lead for the developing of novel hAdoMetDC inhibitors.

Combinatorial Consensus Scoring for Ligand-Based Virtual Fragment Screening: A Comparative Case Study for Serotonin 5-HT(3)A, Histamine H(1), and Histamine H(4) Receptors

In the current study we have evaluated the applicability of ligand-based virtual screening (LBVS) methods for the identification of small fragment-like biologically active molecules using different similarity descriptors and different consensus scoring approaches. For this purpose, we have evaluated the performance of 14 chemical similarity descriptors in retrospective virtual screening studies to discriminate fragment-like ligands of three membrane-bound receptors from fragments that are experimentally determined to have no affinity for these proteins (true inactives). We used a complete fragment affinity data set of experimentally determined ligands and inactives for two G protein-coupled receptors (GPCRs), the histamine H1 receptor (H1R) and the histamine H4 receptor (H4R), and one ligand-gated ion channel (LGIC), the serotonin receptor (5-HT3AR), to validate our retrospective virtual screening studies. We have exhaustively tested consensus scoring strategies that combine the results of multiple actives (group fusion) or combine different similarity descriptors (similarity fusion), and for the first time systematically evaluated different combinations of group fusion and similarity fusion approaches. Our studies show that for these three case study protein targets both consensus scoring approaches can increase virtual screening enrichments compared to single chemical similarity search methods. Our cheminformatics analyses recommend to use a combination of both group fusion and similarity fusion for prospective ligand-based virtual fragment screening.

Improvements, trends, and new ideas in molecular docking: 2012-2013 in review

Molecular docking is a computational method for predicting the placement of ligands in the binding sites of their receptor(s). In this review, we discuss the methodological developments that occurred in the docking field in 2012 and 2013, with a particular focus on the more difficult aspects of this computational discipline. The main challenges and therefore focal points for developments in docking, covered in this review, are receptor flexibility, solvation, scoring, and virtual screening. We specifically deal with such aspects of molecular docking and its applications as selection criteria for constructing receptor ensembles, target dependence of scoring functions, integration of higher-level theory into scoring, implicit and explicit handling of solvation in the binding process, and comparison and evaluation of docking and scoring methods.

In silico design of low molecular weight protein-protein interaction inhibitors: Overall concept and recent advances

Protein-protein interactions (PPIs) are carrying out diverse functions in living systems and are playing a major role in the health and disease states. Low molecular weight (LMW) "drug-like" inhibitors of PPIs would be very valuable not only to enhance our understanding over physiological processes but also for drug discovery endeavors. However, PPIs were deemed intractable by LMW chemicals during many years. But today, with the new experimental and in silico technologies that have been developed, about 50 PPIs have already been inhibited by LMW molecules. Here, we first focus on general concepts about protein-protein interactions, present a consensual view about ligandable pockets at the protein interfaces and the possibilities of using fast and cost effective structure-based virtual screening methods to identify PPI hits. We then discuss the design of compound collections dedicated to PPIs. Recent financial analyses of the field suggest that LMW PPI modulators could be gaining momentum over biologics in the coming years supporting further research in this area.

Molecular docking screening using agonist-bound GPCR structures: probing the A2A adenosine receptor

Crystal structures of G protein-coupled receptors (GPCRs) have recently revealed the molecular basis of ligand binding and activation, which has provided exciting opportunities for structure-based drug design. The A2A adenosine receptor (A2AAR) is a promising therapeutic target for cardiovascular diseases, but progress in this area is limited by the lack of novel agonist scaffolds. We carried out docking screens of 6.7 million commercially available molecules against active-like conformations of the A2AAR to investigate whether these structures could guide the discovery of agonists. Nine out of the 20 predicted agonists were confirmed to be A2AAR ligands, but none of these activated the ARs. The difficulties in discovering AR agonists using structure-based methods originated from limited atomic-level understanding of the activation mechanism and a chemical bias toward antagonists in the screened library. In particular, the composition of the screened library was found to strongly reduce the likelihood of identifying AR agonists, which reflected the high ligand complexity required for receptor activation. Extension of this analysis to other pharmaceutically relevant GPCRs suggested that library screening may not be suitable for targets requiring a complex receptor-ligand interaction network. Our results provide specific directions for the future development of novel A2AAR agonists and general strategies for structure-based drug discovery.

Cell-based and virtual fragment screening for adrenergic α2C receptor agonists

Fragment-based drug discovery has emerged as an alternative to conventional lead identification and optimization strategies generally supported by biophysical detection techniques. Membrane targets like G protein-coupled receptors (GPCRs), however, offer challenges in lack of generic immobilization or stabilization methods for the dynamic, membrane-bound supramolecular complexes. Also modeling of different functional states of GPCRs proved to be a challenging task. Here we report a functional cell-based high concentration screening campaign for the identification of adrenergic α2C receptor agonists compared with the virtual screening of the same ligand set against an active-like homology model of the α2C receptor. The conventional calcium mobilization-based assay identified active fragments with a similar incidence to several other reported fragment screens on GPCRs. 16 out of 3071 screened fragments turned out as specific ligands of α2C, two of which were identified by virtual screening as well and several of the hits possessed surprisingly high affinity and ligand efficiency. Our results indicate that in vitro biological assays can be utilized in the fragment hit identification process for GPCR targets.

Structure-based virtual screening for fragment-like ligands of the G protein-coupled histamine H4 receptor

Structure-based virtual screening using H₁R- and β₂R-based histamine H₄R homology models identified 9 fragments with an affinity ranging from 0.14 to 6.3 μm for H₄R.

Structural and energetic effects of A2A adenosine receptor mutations on agonist and antagonist binding

To predict structural and energetic effects of point mutations on ligand binding is of considerable interest in biochemistry and pharmacology. This is not only useful in connection with site-directed mutagenesis experiments, but could also allow interpretation and prediction of individual responses to drug treatment. For G-protein coupled receptors systematic mutagenesis has provided the major part of functional data as structural information until recently has been very limited. For the pharmacologically important A(2A) adenosine receptor, extensive site-directed mutagenesis data on agonist and antagonist binding is available and crystal structures of both types of complexes have been determined. Here, we employ a computational strategy, based on molecular dynamics free energy simulations, to rationalize and interpret available alanine-scanning experiments for both agonist and antagonist binding to this receptor. These computer simulations show excellent agreement with the experimental data and, most importantly, reveal the molecular details behind the observed effects which are often not immediately evident from the crystal structures. The work further provides a distinct validation of the computational strategy used to assess effects of point-mutations on ligand binding. It also highlights the importance of considering not only protein-ligand interactions but also those mediated by solvent water molecules, in ligand design projects.

Multiple fragment docking and linking in primary and secondary pockets of dopamine receptors

A sequential docking methodology was applied to computationally predict starting points for fragment linking using the human dopamine D3 receptor crystal structure and a human dopamine D2 receptor homology model. Two focused fragment libraries were docked in the primary and secondary binding sites, and best fragment combinations were enumerated. Similar top scoring fragments were found for the primary site, while secondary site fragments were predicted to convey selectivity. Three linked compounds were synthesized that had 9-, 39-, and 55-fold selectivity in favor of D3 and the subtype selectivity of the compounds was assessed on a structural basis.

The receptor concept in 3D: from hypothesis and metaphor to GPCR-ligand structures

The first mentioning of the word "receptor" for the structure with which a bioactive compound should react for obtaining its specific influence on a physiological system goes back to the years around 1900. The receptor concept was adapted from the lock and key theory for the enzyme substrate and blockers interactions. Through the years the concept, in the beginning rather being a metaphor, not a model, was refined and became reality in recent years. Not only the structures of receptors were elucidated, also the receptor machineries were unraveled. Following a brief historical review we will describe how the recent breakthroughs in the experimental determination of G protein-coupled receptor (GPCR) crystal structures can be complemented by computational modeling, medicinal chemistry, biochemical, and molecular pharmacological studies to obtain new insights into the molecular determinants of GPCR-ligand binding and activation. We will furthermore discuss how this information can be used for structure-based discovery of novel GPCR ligands that bind specific (allosteric) binding sites with desired effects on GPCR functional activity.