The relaxed complex method: Accommodating receptor flexibility for drug design with an improved scoring scheme

Synergizing drug repurposing and target identification for neurodegenerative diseases

Despite dedicated research efforts, the absence of disease-curing remedies for neurodegenerative diseases (NDDs) continues to jeopardize human society and stands as a challenge. Drug repurposing is an attempt to find new functionality of existing drugs and take it as an opportunity to discourse the clinically unmet need to treat neurodegeneration. However, despite applying this approach to rediscover a drug, it can also be used to identify the target on which a drug could work. The primary objective of target identification is to unravel all the possibilities of detecting a new drug or repurposing an existing drug. Lately, scientists and researchers have been focusing on specific genes, a particular site in DNA, a protein, or a molecule that might be involved in the pathogenesis of the disease. However, the new era discusses directing the signaling mechanism involved in the disease progression, where receptors, ion channels, enzymes, and other carrier molecules play a huge role. This review aims to highlight how target identification can expedite the whole process of drug repurposing. Here, we first spot various target-identification methods and drug-repositioning studies, including drug-target and structure-based identification studies. Moreover, we emphasize various drug repurposing approaches in NDDs, namely, experimental-based, mechanism-based, and in silico approaches. Later, we draw attention to validation techniques and stress on drugs that are currently undergoing clinical trials in NDDs. Lastly, we underscore the future perspective of synergizing drug repurposing and target identification in NDDs and present an unresolved question to address the issue.

Molecular Dynamics as a Tool for Virtual Ligand Screening

Rational drug design is essential for new drugs to emerge, especially when the structure of a target protein or nucleic acid is known. To that purpose, high-throughput virtual ligand screening campaigns aim at discovering computationally new binding molecules or fragments to modulate particular biomolecular interactions or biological activities, related to a disease process. The structure-based virtual ligand screening process primarily relies on docking methods which allow predicting the binding of a molecule to a biological target structure with a correct conformation and the best possible affinity. The docking method itself is not sufficient as it suffers from several and crucial limitations (lack of full protein flexibility information, no solvation and ion effects, poor scoring functions, and unreliable molecular affinity estimation).At the interface of computer techniques and drug discovery, molecular dynamics (MD) allows introducing protein flexibility before or after a docking protocol, refining the structure of protein-drug complexes in the presence of water, ions, and even in membrane-like environments, describing more precisely the temporal evolution of the biological complex and ranking these complexes with more accurate binding energy calculations. In this chapter, we describe the up-to-date MD, which plays the role of supporting tools in the virtual ligand screening (VS) process.Without a doubt, using docking in combination with MD is an attractive approach in structure-based drug discovery protocols nowadays. It has proved its efficiency through many examples in the literature and is a powerful method to significantly reduce the amount of required wet experimentations (Tarcsay et al, J Chem Inf Model 53:2990-2999, 2013; Barakat et al, PLoS One 7:e51329, 2012; De Vivo et al, J Med Chem 59:4035-4061, 2016; Durrant, McCammon, BMC Biol 9:71-79, 2011; Galeazzi, Curr Comput Aided Drug Des 5:225-240, 2009; Hospital et al, Adv Appl Bioinforma Chem 8:37-47, 2015; Jiang et al, Molecules 20:12769-12786, 2015; Kundu et al, J Mol Graph Model 61:160-174, 2015; Mirza et al, J Mol Graph Model 66:99-107, 2016; Moroy et al, Future Med Chem 7:2317-2331, 2015; Naresh et al, J Mol Graph Model 61:272-280, 2015; Nichols et al, J Chem Inf Model 51:1439-1446, 2011; Nichols et al, Methods Mol Biol 819:93-103, 2012; Okimoto et al, PLoS Comput Biol 5:e1000528, 2009; Rodriguez-Bussey et al, Biopolymers 105:35-42, 2016; Sliwoski et al, Pharmacol Rev 66:334-395, 2014).

Worth the Weight: Sub-Pocket EXplorer (SubPEx), a Weighted Ensemble Method to Enhance Binding-Pocket Conformational Sampling

Structure-based virtual screening (VS) is an effective method for identifying potential small-molecule ligands, but traditional VS approaches consider only a single binding-pocket conformation. Consequently, they struggle to identify ligands that bind to alternate conformations. Ensemble docking helps address this issue by incorporating multiple conformations into the docking process, but it depends on methods that can thoroughly explore pocket flexibility. We here introduce Sub-Pocket EXplorer (SubPEx), an approach that uses weighted ensemble (WE) path sampling to accelerate binding-pocket sampling. As proof of principle, we apply SubPEx to three proteins relevant to drug discovery: heat shock protein 90, influenza neuraminidase, and yeast hexokinase 2. SubPEx is available free of charge without registration under the terms of the open-source MIT license: http://durrantlab.com/subpex/.

Worth the weight: Sub-Pocket EXplorer (SubPEx), a weighted-ensemble method to enhance binding-pocket conformational sampling

Structure-based virtual screening (VS) is an effective method for identifying potential small-molecule ligands, but traditional VS approaches consider only a single binding-pocket conformation. Consequently, they struggle to identify ligands that bind to alternate conformations. Ensemble docking helps address this issue by incorporating multiple conformations into the docking process, but it depends on methods that can thoroughly explore pocket flexibility. We here introduce Sub-Pocket EXplorer (SubPEx), an approach that uses weighted ensemble (WE) path sampling to accelerate binding-pocket sampling. As proof of principle, we apply SubPEx to three proteins relevant to drug discovery: heat shock protein 90, influenza neuraminidase, and yeast hexokinase 2. SubPEx is available free of charge without registration under the terms of the open-source MIT license: http://durrantlab.com/subpex/.

A Curvilinear-Path Umbrella Sampling Approach to Characterizing the Interactions Between Rapamycin and Three FKBP12 Variants

Rapamycin is an immunosuppressant macrolide that exhibits anti-proliferative properties through inhibiting the mTOR kinase. In fact, the drug first associates with the FKBP12 enzyme before interacting with the FRB domain of its target. Despite the availability of structural and thermodynamic information on the interaction of FKBP12 with rapamycin, the energetic and mechanistic understanding of this process is still incomplete. We recently reported a multiple-walker umbrella sampling simulation approach to characterizing the protein–protein interaction energetics along curvilinear paths. In the present paper, we extend our investigations to a protein-small molecule duo, the FKBP12•rapamycin complex. We estimate the binding free energies of rapamycin with wild-type FKBP12 and two mutants in which a hydrogen bond has been removed, D37V and Y82F. Furthermore, the underlying mechanistic details are analyzed. The calculated standard free energies of binding agree well with the experimental data, and the roles of the hydrogen bonds are shown to be quite different for each of these two mutated residues. On one hand, removing the carboxylate group of D37 strongly destabilizes the association; on the other hand, the hydroxyl group of Y82 is nearly unnecessary for the stability of the complex because some nonconventional, cryptic, indirect interaction mechanisms seem to be at work.

Accommodating receptor flexibility and free energy calculation to reduce false positive binders in the discovery of natural products blockers of SARS-COV-2 spike RBD-ACE2 interface

The emergence of severe acute respiratory syndrome coronavirus 2 (SARS-COV-2), which causes coronavirus disease-19 (COVID-19) has caused more than 2 million deaths around the globe. The high transmissibility rate of the disease is related to the strong interaction between the virus spike receptor-binding domain (RBD) and the human angiotensin-converting enzyme 2 (ACE2) as documented in several reports. In this study, using state-of-the-art computational methods, natural products were screened and their molecular mechanism to disrupt spike RBD-ACE2 recognition was evaluated. There is the sensitivity of results to receptor ensemble docking calculations. Binding free energy and MD simulation are important tools to evaluate the thermodynamics of binding stability and the capacity of top hits to disrupt RBD-ACE2 recognition. The free energy profiles provide a slight decrease in binding affinity of the virus-receptor interaction. Three flavonoids parvisoflavone B (3), alpinumisoflavone (5) and norisojamicin (2) were effective in blocking the viral entry by binding strongly at the spike RBD-ACE2 interface with the inhibition constant of 0.56, 0.78 and 0.93 μM, respectively. The same compounds demonstrated similar effect on free ACE2 protein. Compound (2), also demonstrated ability to bind strongly on free spike RBD. Well-tempered metadynamics established that parvisoflavone B (3) works by binding to three sites namely interface α, β and loop thereby inhibiting viral cell entry. Owing to their desirable pharmacokinetic properties, the presented top hit natural products are suggested for further SARS-COV-2 molecular targets and subsequent in vitro and in vivo evaluations.

Structure-based in silico approaches for drug discovery against Mycobacterium tuberculosis

Mycobacterium tuberculosis is the causative agent of TB and was estimated to cause 1.4 million death in 2019, alongside 10 million new infections. Drug resistance is a growing issue, with multi-drug resistant infections representing 3.3% of all new infections, hence novel antimycobacterial drugs are urgently required to combat this growing health emergency. Alongside this, increased knowledge of gene essentiality in the pathogenic organism and larger compound databases can aid in the discovery of new drug compounds. The number of protein structures, X-ray based and modelled, is increasing and now accounts for greater than > 80% of all predicted M. tuberculosis proteins; allowing novel targets to be investigated. This review will focus on structure-based in silico approaches for drug discovery, covering a range of complexities and computational demands, with associated antimycobacterial examples. This includes molecular docking, molecular dynamic simulations, ensemble docking and free energy calculations. Applications of machine learning onto each of these approaches will be discussed. The need for experimental validation of computational hits is an essential component, which is unfortunately missing from many current studies. The future outlooks of these approaches will also be discussed.

Identification of Novel Cyclin A2 Binding Site and Nanomolar Inhibitors of Cyclin A2-CDK2 Complex

BACKGROUND

Given the diverse roles of cyclin A2 both in cell cycle regulation and in DNA damage response, identifying small molecule regulators of cyclin A2 activity carries significant potential to regulate diverse cellular processes in both ageing/neurodegeneration and in cancer.

OBJECTIVE

Based on cyclin A2's recently discovered role in DNA repair, we hypothesized that small molecule inhibitors that were predicted to bind to both cyclin A2 and CDK2 will be useful as a radiosensitizer of cancer cells. In this study, we used structure-based drug discovery to find inhibitors that target both cyclin A2 and CDK2.

METHODS

Molecular dynamics simulations were used to generate diverse binding pocket conformations for application of the relaxed complex scheme. We then used structure-based virtual screening to find potential dual cyclin A2 and CDK2 inhibitors. Based on a consensus score of docked poses from Glide and AutoDock Vina, we identified about 40 promising hit compounds, where all PAINS scaffolds were removed from consideration. A biochemical luminescence assay of cyclin A2-CDK2 function was used for experimental verification.

RESULTS

Four lead inhibitors of cyclin A2-CDK2 complex have been identified using a relaxed complex scheme virtual screen have been verified in a biochemical luminescence assay of cyclin A2- CDK2 function. Two of the four lead inhibitors had inhibitory concentrations in the nanomolar range.

CONCLUSION

The four cyclin A2-CDK2 complex inhibitors are the first reported inhibitors that were specifically designed not to target the cyclin A2-CDK2 protein-protein interface. Overall, our results highlight the potential of combined advanced computational tools and biochemical verification to discover novel binding scaffolds.

Accelerating AutoDock4 with GPUs and Gradient-Based Local Search

AutoDock4 is a widely used program for docking small molecules to macromolecular targets. It describes ligand-receptor interactions using a physics-inspired scoring function that has been proven useful in a variety of drug discovery projects. However, compared to more modern and recent software, AutoDock4 has longer execution times, limiting its applicability to large scale dockings. To address this problem, we describe an OpenCL implementation of AutoDock4, called AutoDock-GPU, that leverages the highly parallel architecture of GPU hardware to reduce docking runtime by up to 350-fold with respect to a single-threaded process. Moreover, we introduce the gradient-based local search method ADADELTA, as well as an improved version of the Solis-Wets random optimizer from AutoDock4. These efficient local search algorithms significantly reduce the number of calls to the scoring function that are needed to produce good results. The improvements reported here, both in terms of docking throughput and search efficiency, facilitate the use of the AutoDock4 scoring function in large scale virtual screening.

The AutoDock suite at 30

The AutoDock suite provides a comprehensive toolset for computational ligand docking and drug design and development. The suite builds on 30 years of methods development, including empirical free energy force fields, docking engines, methods for site prediction, and interactive tools for visualization and analysis. Specialized tools are available for challenging systems, including covalent inhibitors, peptides, compounds with macrocycles, systems where ordered hydration plays a key role, and systems with substantial receptor flexibility. All methods in the AutoDock suite are freely available for use and reuse, which has engendered the continued growth of a diverse community of primary users and third-party developers.

Discovery of Novel Small-Molecule Calcium Sensitizers for Cardiac Troponin C: A Combined Virtual and Experimental Screening Approach

Heart failure is a leading cause of death throughout the world and is triggered by a disruption of the cardiac contractile machinery. This machinery is regulated in a calcium-dependent manner by the protein complex troponin. Calcium binds to the N-terminal domain of cardiac troponin C (cNTnC) setting into motion the cascade of events leading to muscle contraction. Because of the severity and prevalence of heart failure, there is a strong need to develop small-molecule therapeutics designed to increase the calcium sensitivity of cardiac troponin in order to treat this devastating condition. Molecules that are able to stabilize an open configuration of cNTnC and additionally facilitate the binding of the cardiac troponin I (cTnI) switch peptide have the potential to enable increased calcium sensitization and strengthened cardiac function. Here, we employed a high throughput virtual screening methodology built upon the ability of computational docking to reproduce known experimental results and to accurately recognize cNTnC conformations conducive to small molecule binding using a receiver operator characteristic curve analysis. This approach combined with concurrent stopped-flow kinetic experimental verification led to the identification of a number of sensitizers, which slowed the calcium off-rate. An initial hit, compound 4, was identified with medium affinity (84 ± 30 μM). Through refinement, a calcium sensitizing agent, compound 5, with an apparent affinity of 1.45 ± 0.09 μM was discovered. This molecule is one of the highest affinity calcium sensitizers known to date.

Using machine learning to improve ensemble docking for drug discovery

Ensemble docking has provided an inexpensive method to account for receptor flexibility in molecular docking for virtual screening. Unfortunately, as there is no rigorous theory to connect the docking scores from multiple structures to measured activity, researchers have not yet come up with effective ways to use these scores to classify compounds into actives and inactives. This shortcoming has led to the decrease, rather than an increase in the performance of classifying compounds when more structures are added to the ensemble. Previously, we suggested machine learning, implemented in the form of a naïve Bayesian model could alleviate this problem. However, the naïve Bayesian model assumed that the probabilities of observing the docking scores to different structures to be independent. This approximation might prevent it from achieving even higher performance. In the work presented in this paper, we have relaxed this approximation when using several other machine learning methods-k nearest neighbor, logistic regression, support vector machine, and random forest-to improve ensemble docking. We found significant improvement.

Improving ensemble docking for drug discovery by machine learning

Ensemble docking has provided an inexpensive method to account for receptor flexibility in molecular docking. However, it is still unclear how best to use the docking scores from multiple structures to classify compounds into actives and inactives. Previous studies have also found that the performance of classification could decrease rather than increase with the number of structures included in the ensemble. Machine learning could help to alleviate these problems.

The compromise of virtual screening and its impact on drug discovery

Introduction: Docking and structure-based virtual screening (VS) have been standard approaches in structure-based design for over two decades. However, our understanding of the limitations, potential, and strength of these techniques has enhanced, raising expectations. Areas covered: Based on a survey of reports in the past five years, we assess whether VS: (1) predicts binding poses in agreement with crystallographic data (when available); (2) is a superior screening tool, as often claimed; (3) is successful in identifying chemical scaffolds that can be starting points for subsequent lead optimization cycles. Data shows that knowledge of the target and its chemotypes in postprocessing lead to viable hits in early drug discovery endeavors. Expert opinion: VS is capable of accurate placements in the pocket for the most part, but does not consistently score screening collections accurately. What matters is capitalization on available resources to get closer to a viable lead or optimizable series. Integration of approaches, subjective hit selection guided by knowledge of the receptor or endogenous ligand, libraries driven by experimental guides, validation studies to identify the best docking/scoring that reproduces experimental findings, constraints regarding receptor-ligand interactions, thoroughly designed methodologies, and predefined cutoff scoring criteria strengthen VS's position in pharmaceutical research.

Docking simulation and antibiotic discovery targeting the MlaC protein in Gram-negative bacteria

To maintain the lipid asymmetry of the cell envelope in Gram-negative bacteria, the MlaC protein serves as a lipid transfer factor and delivers phospholipids from the outer to the inner membrane. A strategy of antibiotic discovery is to design a proper compound that can tightly bind to the MlaC protein and inhibit the MlaC function. In this study, we performed virtual screening on multiple MlaC structures obtained from molecular dynamics simulations to identify potential MlaC binders. Our results suggested that clorobiocin is a compound that could bind to the MlaC protein. Through the comparison of the bound geometry between clorobiocin and novobiocin, we pointed out that the methyl-pyrrole group of the noviose sugar in clorobiocin forms hydrophobic interactions with amino acids in the phospholipid binding pocket, which allows the compound to bind deep in the active site. This also explains why clorobiocin shows a tighter binding affinity than novobiocin. Our study highlights a practical path of antibiotic development against Gram-negative bacteria.

Determination of Ligand Binding Epitope Structures Using Polarization Transfer from Hyperpolarized Ligands

Drug discovery processes require the determination of the protein binding site structure, which can be achieved via nuclear magnetic resonance (NMR) spectroscopy. While traditional NMR spectroscopy suffers from low sensitivity, NMR signals can be significantly enhanced through hyperpolarization of nuclear spins. Here, folic acid is hyperpolarized by dissolution dynamic nuclear polarization (D-DNP). Polarization transfer to dihydrofolate reductase is compared to signal evolution predicted for docking-derived structures. The results demonstrate that a scoring function derived from the experimental data improves the ranking of structures. With data from six methyl groups, Spearman's correlation coefficient of the experimental scoring function to the root-mean-square deviation from a reference structure is 0.88 for five individually addressed ligand protons and 0.59 for the entire ligand, while the same correlation coefficient of the energy calculated from docking alone is 0.49. D-DNP NMR-derived ranking, therefore, is capable of determining the ligand structure with a small number of individually addressed source spins.

Molecular modelling and competitive inhibition of a Mycobacterium tuberculosis multidrug-resistance efflux pump

Tuberculosis is a major cause of mortality and morbidity in developing countries, and the emergency of multidrug and extensive drug resistance cases is an utmost issue on the control of the disease. Despite the efforts on the development of new antibiotics, eventually there will be strains resistant to them as well. Efflux plays an important role in the evolution of resistance in Mycobacterium tuberculosis. Tap is an important efflux pump associated with tuberculosis resistant to isoniazid, rifampicine and ofloxacin and with multidrug resistance. The development of efflux inhibitors for Tap could raise the effectiveness of second line drugs and reduce the duration of the current treatment. Therefore the objective of this study is to build a reliable molecular model of Tap efflux pump and test the possible competitive inhibition between efflux inhibitors and antibiotics in the optimized structure. We built twenty five Tap models with molecular modelling to elect the best according to the results of the validation analysis. The elected model went through to a 100 ns molecular dynamics simulation in a lipid bilayer, and the resulting optimized structure was used in docking studies to test if the used efflux inhibitors may act via competitive inhibition on antibiotics. The validation results pointed the model built by Phyre2 as the closest to a possible native Tap structure, and therefore it was the elected model. RSMD analysis revealed the model is stable, where the predicted binding site stabilized between 15 and 20 ns, maintaining the RMSD at around 0.35 Å throughout the molecular dynamics simulation in a lipid bilayer. Therefore this model is reliable and can also be used for further studies. The docking studies showed a possibility of competitive inhibition by NUNL02 on ofloxacin and bedaquiline, and by verapamil on ofloxacin and rifampicin. This presents the possibility that NUNL02 and verapamil are possible inhibitors of Tap efflux and highlights the importance of including efflux inhibitors as adjuvants to the tuberculosis therapy, as it indicates a possible extrusion of ofloxacin, rifampicin and bedaquilin by Tap.

Efficiency of Stratification for Ensemble Docking Using Reduced Ensembles

Molecular docking can account for receptor flexibility by combining the docking score over multiple rigid receptor conformations, such as snapshots from a molecular dynamics simulation. Here, we evaluate a number of common snapshot selection strategies using a quality metric from stratified sampling, the efficiency of stratification, which compares the variance of a selection strategy to simple random sampling. We also extend the metric to estimators of exponential averages (which involve an exponential transformation, averaging, and inverse transformation) and minima. For docking sets of over 500 ligands to four different proteins of varying flexibility, we observe that, for estimating ensemble averages and exponential averages, many clustering algorithms have similar performance trends: for a few snapshots (less than 25), medoids are the most efficient, while, for a larger number, optimal (the allocation that minimizes the variance) and proportional (to the size of each cluster) allocation become more efficient. Proportional allocation appears to be the most consistently efficient for estimating minima.

Bridging Molecular Docking to Molecular Dynamics in Exploring Ligand-Protein Recognition Process: An Overview

Computational techniques have been applied in the drug discovery pipeline since the 1980s. Given the low computational resources of the time, the first molecular modeling strategies relied on a rigid view of the ligand-target binding process. During the years, the evolution of hardware technologies has gradually allowed simulating the dynamic nature of the binding event. In this work, we present an overview of the evolution of structure-based drug discovery techniques in the study of ligand-target recognition phenomenon, going from the static molecular docking toward enhanced molecular dynamics strategies.

Plasticity of the Binding Site of Renin: Optimized Selection of Protein Structures for Ensemble Docking

Protein flexibility poses a major challenge to docking of potential ligands in that the binding site can adopt different shapes. Docking algorithms usually keep the protein rigid and only allow the ligand to be treated as flexible. However, a wrong assessment of the shape of the binding pocket can prevent a ligand from adapting a correct pose. Ensemble docking is a simple yet promising method to solve this problem: Ligands are docked into multiple structures, and the results are subsequently merged. Selection of protein structures is a significant factor for this approach. In this work we perform a comprehensive and comparative study evaluating the impact of structure selection on ensemble docking. We perform ensemble docking with several crystal structures and with structures derived from molecular dynamics simulations of renin, an attractive target for antihypertensive drugs. Here, 500 ns of MD simulations revealed binding site shapes not found in any available crystal structure. We evaluate the importance of structure selection for ensemble docking by comparing binding pose prediction, ability to rank actives above nonactives (screening utility), and scoring accuracy. As a result, for ensemble definition k-means clustering appears to be better suited than hierarchical clustering with average linkage. The best performing ensemble consists of four crystal structures and is able to reproduce the native ligand poses better than any individual crystal structure. Moreover this ensemble outperforms 88% of all individual crystal structures in terms of screening utility as well as scoring accuracy. Similarly, ensembles of MD-derived structures perform on average better than 75% of any individual crystal structure in terms of scoring accuracy at all inspected ensembles sizes.

Molecular Dynamics as a Tool for Virtual Ligand Screening

Rational drug design is essential for new drugs to emerge, especially when the structure of a target protein or catalytic enzyme is known experimentally. To that purpose, high-throughput virtual ligand screening campaigns aim at discovering computationally new binding molecules or fragments to inhibit a particular protein interaction or biological activity. The virtual ligand screening process often relies on docking methods which allow predicting the binding of a molecule into a biological target structure with a correct conformation and the best possible affinity. The docking method itself is not sufficient as it suffers from several and crucial limitations (lack of protein flexibility information, no solvation effects, poor scoring functions, and unreliable molecular affinity estimation).At the interface of computer techniques and drug discovery, molecular dynamics (MD) allows introducing protein flexibility before or after a docking protocol, refining the structure of protein-drug complexes in the presence of water, ions and even in membrane-like environments, and ranking complexes with more accurate binding energy calculations. In this chapter we describe the up-to-date MD protocols that are mandatory supporting tools in the virtual ligand screening (VS) process. Using docking in combination with MD is one of the best computer-aided drug design protocols nowadays. It has proved its efficiency through many examples, described below.

Successful Identification of Cardiac Troponin Calcium Sensitizers Using a Combination of Virtual Screening and ROC Analysis of Known Troponin C Binders

Calcium-dependent cardiac muscle contraction is regulated by the protein complex troponin. Calcium binds to the N-terminal domain of troponin C (cNTnC) which initiates the process of contraction. Heart failure is a consequence of a disruption of this process. With the prevalence of this condition, a strong need exists to find novel compounds to increase the calcium sensitivity of cNTnC. Desirable are small chemical molecules that bind to the interface between cTnC and the cTnI switch peptide and exhibit calcium sensitizing properties by possibly stabilizing cTnC in an open conformation. To identify novel drug candidates, we employed a structure-based drug discovery protocol that incorporated the use of a relaxed complex scheme (RCS). In preparation for the virtual screening, cNTnC conformations were identified based on their ability to correctly predict known cNTnC binders using a receiver operating characteristics analysis. Following a virtual screen of the National Cancer Institute's Developmental Therapeutic Program database, a small number of molecules were experimentally tested using stopped-flow kinetics and steady-state fluorescence titrations. We identified two novel compounds, 3-(4-methoxyphenyl)-6,7-chromanediol (NSC600285) and 3-(4-methylphenyl)-7,8-chromanediol (NSC611817), that show increased calcium sensitivity of cTnC in the presence of the regulatory domain of cTnI. The effects of NSC600285 and NSC611817 on the calcium dissociation rate was stronger than that of the known calcium sensitizer bepridil. Thus, we identified a 3-phenylchromane group as a possible key pharmacophore in the sensitization of cardiac muscle contraction. Building on this finding is of interest to researchers working on development of drugs for calcium sensitization.

Protein structure-based drug design: from docking to molecular dynamics

Recent years have witnessed rapid developments of computer-aided drug design methods, which have reached accuracy that allows their routine practical applications in drug discovery campaigns. Protein structure-based methods are useful for the prediction of binding modes of small molecules and their relative affinity. The high-throughput docking of up to 10⁶ small molecules followed by scoring based on implicit-solvent force field can robustly identify micromolar binders using a rigid protein target. Molecular dynamics with explicit solvent is a low-throughput technique for the characterization of flexible binding sites and accurate evaluation of binding pathways, kinetics, and thermodynamics. In this review we highlight recent advancements in applications of ligand docking tools and molecular dynamics simulations to ligand identification and optimization.

NMR structure-based optimization of Staphylococcus aureus sortase A pyridazinone inhibitors

Staphylococcus aureus is a leading cause of hospital-acquired infections in the USA and is a major health concern as methicillin-resistant S. aureus and other antibiotic-resistant strains are common. Compounds that inhibit the S. aureus sortase (SrtA) cysteine transpeptidase may function as potent anti-infective agents as this enzyme attaches virulence factors to the bacterial cell wall. While a variety of SrtA inhibitors have been discovered, the vast majority of these small molecules have not been optimized using structure-based approaches. Here we have used NMR spectroscopy to determine the molecular basis through which pyridazinone-based small molecules inhibit SrtA. These inhibitors covalently modify the active cysteine thiol and partially mimic the natural substrate of SrtA by inducing the closure of an active site loop. Computational and synthetic chemistry methods led to second-generation analogues that are ~70-fold more potent than the lead molecule. These optimized molecules exhibit broad-spectrum activity against other types of class A sortases, have reduced cytotoxicity, and impair SrtA-mediated protein display on S. aureus cell surface. Our work shows that pyridazinone analogues are attractive candidates for further development into anti-infective agents, and highlights the utility of employing NMR spectroscopy and solubility-optimized small molecules in structure-based drug discovery.

Computational methods in drug discovery

The process for drug discovery and development is challenging, time consuming and expensive. Computer-aided drug discovery (CADD) tools can act as a virtual shortcut, assisting in the expedition of this long process and potentially reducing the cost of research and development. Today CADD has become an effective and indispensable tool in therapeutic development. The human genome project has made available a substantial amount of sequence data that can be used in various drug discovery projects. Additionally, increasing knowledge of biological structures, as well as increasing computer power have made it possible to use computational methods effectively in various phases of the drug discovery and development pipeline. The importance of in silico tools is greater than ever before and has advanced pharmaceutical research. Here we present an overview of computational methods used in different facets of drug discovery and highlight some of the recent successes. In this review, both structure-based and ligand-based drug discovery methods are discussed. Advances in virtual high-throughput screening, protein structure prediction methods, protein-ligand docking, pharmacophore modeling and QSAR techniques are reviewed.

Computational protein-ligand docking and virtual drug screening with the AutoDock suite

Computational docking can be used to predict bound conformations and free energies of binding for small-molecule ligands to macromolecular targets. Docking is widely used for the study of biomolecular interactions and mechanisms, and it is applied to structure-based drug design. The methods are fast enough to allow virtual screening of ligand libraries containing tens of thousands of compounds. This protocol covers the docking and virtual screening methods provided by the AutoDock suite of programs, including a basic docking of a drug molecule with an anticancer target, a virtual screen of this target with a small ligand library, docking with selective receptor flexibility, active site prediction and docking with explicit hydration. The entire protocol will require ∼5 h.

Structure-Based Virtual Ligand Screening on the XRCC4/DNA Ligase IV Interface

The association of DNA Ligase IV (Lig4) with XRCC4 is essential for repair of DNA double-strand breaks (DSBs) by Non-homologous end-joining (NHEJ) in humans. DSBs cytotoxicity is largely exploited in anticancer therapy. Thus, NHEJ is an attractive target for strategies aimed at increasing the sensitivity of tumors to clastogenic anticancer treatments. However the high affinity of the XRCC4/Lig4 interaction and the extended protein-protein interface make drug screening on this target particularly challenging. Here, we conducted a pioneering study aimed at interfering with XRCC4/Lig4 assembly. By Molecular Dynamics simulation using the crystal structure of the complex, we first delineated the Lig4 clamp domain as a limited suitable target. Then, we performed in silico screening of ~95,000 filtered molecules on this Lig4 subdomain. Hits were evaluated by Differential Scanning Fluorimetry, Saturation Transfer Difference-NMR spectroscopy and interaction assays with purified recombinant proteins. In this way we identified the first molecule able to prevent Lig4 binding to XRCC4 in vitro. This compound has a unique tripartite interaction with the Lig4 clamp domain that suggests a starting chemotype for rational design of analogous molecules with improved affinity.

Conformational transition paths harbor structures useful for aiding drug discovery and understanding enzymatic mechanisms in protein kinases

This short article examines the usefulness of fast simulations of conformational transition paths in elucidating enzymatic mechanisms and guiding drug discovery for protein kinases. It applies the transition path method in the MOIL software package to simulate the paths of conformational transitions between six pairs of structures from the Protein Data Bank. The structures along the transition paths were found to resemble experimental structures that mimic transient structures believed to form during enzymatic catalysis or conformational transitions, or structures that have drug candidates bound. These findings suggest that such simulations could provide quick initial insights into the enzymatic mechanisms or pathways of conformational transitions of proteins kinases, or could provide structures useful for aiding structure-based drug design.

Flexible receptor docking for drug discovery

INTRODUCTION

Molecular docking has become a popular method for virtual screening. Docking small molecules to a rigid biological receptor is fast but could produce many false negatives and identify less diverse compounds. Flexible receptor docking has alleviated this problem.

AREAS COVERED

This article focuses on reviewing ensemble docking as an approximate but inexpensive method to incorporate receptor flexibility in molecular docking. It outlines key features and recent advances of this method and points out problem areas that need to be addressed to make it even more useful in drug discovery.

EXPERT OPINION

Among the different methods introduced for flexible receptor docking, ensemble docking represents one of the most popular approaches, especially for high-throughput virtual screening. One can generate structural ensembles by using experimental structures, by structural modeling and by various types of molecular simulations. In building a structural ensemble, a judicious choice of the structures to be included can improve performance. Furthermore, reducing the size of the structural ensemble can cut computational costs, and removing the structures that can bind few ligands well could enrich the number of true actives identified by ensemble docking. The ability of ensemble docking to identify more true positives at the top of a rank-ordered list also depends on the choice of the methods to score and rank compounds, an area that needs further research.

Molecular docking to flexible targets

It is widely accepted that protein receptors exist as an ensemble of conformations in solution. How best to incorporate receptor flexibility into virtual screening protocols used for drug discovery remains a significant challenge. Here, stepwise methodologies are described to generate and select relevant protein conformations for virtual screening in the context of the relaxed complex scheme (RCS), to design small molecule libraries for docking, and to perform statistical analyses on the virtual screening results. Methods include equidistant spacing, RMSD-based clustering, and QR factorization protocols for ensemble generation and ROC analysis for ensemble selection.

Protein Flexibility in Drug Discovery: From Theory to Computation

Nowadays it is widely accepted that the mechanisms of biomolecular recognition are strongly coupled to the intrinsic dynamic of proteins. In past years, this evidence has prompted the development of theoretical models of recognition able to describe ligand binding assisted by protein conformational changes. On a different perspective, the need to take into account protein flexibility in structure-based drug discovery has stimulated the development of several and extremely diversified computational methods. Herein, on the basis of a parallel between the major recognition models and the simulation strategies used to account for protein flexibility in ligand binding, we sort out and describe the most innovative and promising implementations for structure-based drug discovery.

Structure, dynamics, and interaction of Mycobacterium tuberculosis (Mtb) DprE1 and DprE2 examined by molecular modeling, simulation, and electrostatic studies

The enzymes decaprenylphosphoryl-β-D-ribose oxidase (DprE1) and decaprenylphosphoryl-β-D-ribose-2-epimerase (DprE2) catalyze epimerization of decaprenylphosporyl ribose (DPR) todecaprenylphosporyl arabinose (DPA) and are critical for the survival of Mtb. Crystal structures of DprE1 so far reported display significant disordered regions and no structural information is known for DprE2. We used homology modeling, protein threading, molecular docking and dynamics studies to investigate the structural and dynamic features of Mtb DprE1 and DprE2 and DprE1-DprE2 complex. A three-dimensional model for DprE2 was generated using the threading approach coupled with ab initio modeling. A 50 ns simulation of DprE1 and DprE2 revealed the overall stability of the structures. Principal Component Analysis (PCA) demonstrated the convergence of sampling in both DprE1 and DprE2. In DprE1, residues in the 269–330 area showed considerable fluctuation in agreement with the regions of disorder observed in the reported crystal structures. In DprE2, large fluctuations were detected in residues 95–113, 146–157, and 197–226. The study combined docking and MD simulation studies to map and characterize the key residues involved in DprE1-DprE2 interaction. A 60 ns MD simulation for DprE1-DprE2 complex was also performed. Analysis of data revealed that the docked complex is stabilized by H-bonding, hydrophobic and ionic interactions. The key residues of DprE1 involved in DprE1-DprE2 interactions belong to the disordered region. We also examined the docked complex of DprE1-BTZ043 to investigate the binding pocket of DprE1 and its interactions with the inhibitor BTZ043. In summary, we hypothesize that DprE1-DprE2 interaction is crucial for the synthesis of DPA and DprE1-DprE2 complex may be a new therapeutic target amenable to pharmacological validation. The findings have important implications in tuberculosis (TB) drug discovery and will facilitate drug development efforts against TB.

Scoring functions for fragment-based drug discovery

Fragment-based drug design represents a challenge for computational drug design because almost inevitably fragments will be weak binders to the biomolecular targets of a specific disease, and the performances of the scoring functions for weak binders are usually poorer than those for the stronger binders. This protocol describes how to predict the binding modes and binding affinities of fragments towards their binding partner with our refined AutoDock scoring function incorporating a quantum chemical charge model, namely, the restrained electrostatic potential (RESP) model. This scoring function was calibrated by robust regression analysis and has been demonstrated to perform well for general classes of protein-ligand interactions and for weak binders (with root-mean square of error of about 2.1 kcal/mol).

Virtual screening strategies: recent advances in the identification and design of anti-cancer agents

Virtual screening (VS) is a well-established technique, which is now routinely employed in computer aided drug designing process. VS can be broadly classified into two categories, i.e., ligand-based and structure-based approach. In recent years, VS has emerged as a time saving and cost effective technique, capable of screening millions of compounds in a user friendly manner. In the area of cancer drug design, VS methods have been widely used and helped in identifying novel molecules as potential anti-cancer agents. Both ligand-based VS (LBVS) structure-based VS (SBVS) methods have been highly useful in the identification of a number of potential anti-cancer agents exhibiting activities in nanomolar range. In tune with the rapid progress in the enhancement of computational power, VS has witnessed significant change in terms of speed and hit rate and in future it is expected that VS will be a preferential alternative to high throughput screening (HTS). This review, discusses recent trends and contribution of VS in the area of anti-cancer drug discovery.

Drug-induced conformational population shifts in topoisomerase-DNA ternary complexes

Type II topoisomerases (TOP2) are enzymes that resolve the topological problems during DNA replication and transcription by transiently cleaving both strands and forming a cleavage complex with the DNA. Several prominent anti-cancer agents inhibit TOP2 by stabilizing the cleavage complex and engendering permanent DNA breakage. To discriminate drug binding modes in TOP2-α and TOP2-β, we applied our newly developed scoring function, dubbed AutoDock4^RAP, to evaluate the binding modes of VP-16, m-AMSA, and mitoxantrone to the cleavage complexes. Docking reproduced crystallographic binding mode of VP-16 in a ternary complex of TOP2-β with root-mean-square deviation of 0.65 Å. Molecular dynamics simulation of the complex confirmed the crystallographic binding mode of VP-16 and the conformation of the residue R503. Drug-related conformational changes in R503 have been observed in ternary complexes with m-AMSA and mitoxantrone. However, the R503 rotamers in these two simulations deviate from their crystallographic conformations, indicating a relaxation dynamics from the conformations determined with the drug replacement procedure. The binding mode of VP-16 in the cleavage complex of TOP2-α was determined by the conjoint use of docking and molecular dynamics simulations, which fell within a similar binding pocket of TOP2-β cleavage complex. Our findings may facilitate more efficient design efforts targeting TOP2-α specific drugs.

Discovery of Staphylococcus aureus sortase A inhibitors using virtual screening and the relaxed complex scheme

Staphylococcus aureus is the leading cause of hospital-acquired infections in the United States. The emergence of multidrug-resistant strains of S. aureus has created an urgent need for new antibiotics. Staphylococcus aureus uses the sortase A enzyme to display surface virulence factors suggesting that compounds that inhibit its activity will function as potent anti-infective agents. Here, we report the identification of several inhibitors of sortase A using virtual screening methods that employ the relaxed complex scheme, an advanced computer-docking methodology that accounts for protein receptor flexibility. Experimental testing validates that several compounds identified in the screen inhibit the activity of sortase A. A lead compound based on the 2-phenyl-2,3-dihydro-1H-perimidine scaffold is particularly promising, and its binding mechanism was further investigated using molecular dynamics simulations and conducting preliminary structure-activity relationship studies.

Computational delineation of tyrosyl-substrate recognition and catalytic landscapes by the epidermal growth factor receptor tyrosine kinase domain

The epidermal growth factor receptor (EGFR) is a receptor tyrosine kinase (RTK), which catalyzes protein phosphorylation reactions by transferring the γ-phosphoryl group from an ATP molecule to the hydroxyl group of tyrosine residues in protein substrates. EGFR is an important drug target in the treatment of cancers and a better understanding of the receptor function is critical to discern cancer mechanisms. We employ a suite of molecular simulation methods to explore the mechanism of substrate recognition and to delineate the catalytic landscape of the phosphoryl transfer reaction. Based on our results, we propose that a highly conserved region corresponding to Val852-Pro853-Ile854-Lys855-Trp856 in the EGFR tyrosine kinase domain (TKD) is essential for substrate binding. We also provide a possible explanation for the established experimental observation that protein tyrosine kinases (including EGFR) select substrates with a glutamic acid at the P - 1 position and a large hydrophobic amino acid at the P + 1 position. Furthermore, our mixed quantum mechanics/molecular mechanics (QM/MM) simulations show that the EGFR protein kinase favors the dissociative mechanism, although an alternative channel through the formation of an associative transition state is also possible. Our simulations establish some key molecular rules in the operation for substrate-recognition and for phosphoryl transfer in the EGFR TKD.

Discovery of novel inhibitors of HIV-1 reverse transcriptase through virtual screening of experimental and theoretical ensembles

Non-nucleoside reverse transcriptase inhibitors (NNRTIs) are potent anti-HIV chemotherapeutics. Although there are FDA-approved NNRTIs, challenges such as the development of resistance have limited their utility. Here, we describe the identification of novel NNRTIs through a combination of computational and experimental approaches. Based on the known plasticity of the NNRTI binding pocket (NNIBP), we adopted an ensemble-based virtual screening strategy: coupling receptor conformations from 10 X-ray crystal structures with 120 snapshots from a total of 480 ns of molecular dynamics (MD) trajectories. A screening library of 2864 National Cancer Institute (NCI) compounds was built and docked against the ensembles in a hierarchical fashion. Sixteen diverse compounds were tested for their ability to block HIV infection in human tissue cultures using a luciferase-based reporter assay. Three promising compounds were further characterized, using a HIV-1 RT-based polymerase assay, to determine the specific mechanism of inhibition. We found that 2 of the three compounds inhibited the polymerase activity of RT (with potency similar to the positive control, the FDA-approved drug nevirapine). Through a computational approach, we were able to discover two compounds which inhibit HIV replication and block the activity of RT, thus offering the potential for optimization into mature inhibitors.

Drug screening strategy for human membrane proteins: from NMR protein backbone structure to in silica- and NMR-screened hits

About 8000 genes encode membrane proteins in the human genome. The information about their druggability will be very useful to facilitate drug discovery and development. The main problem, however, consists of limited structural and functional information about these proteins because they are difficult to produce biochemically and to study. In this paper we describe the strategy that combines Cell-free protein expression, NMR spectroscopy, and molecular DYnamics simulation (CNDY) techniques. Results of a pilot CNDY experiment provide us with a guiding light towards expedited identification of the hit compounds against a new uncharacterized membrane protein as a potentially druggable target. These hits can then be further characterized and optimized to develop the initial lead compound quicker. We illustrate such "omics" approach for drug discovery with the CNDY strategy applied to two example proteins: hypoxia-induced genes HIGD1A and HIGD1B.

Virtual screening with AutoDock Vina and the common pharmacophore engine of a low diversity library of fragments and hits against the three allosteric sites of HIV integrase: participation in the SAMPL4 protein-ligand binding challenge

To rigorously assess the tools and protocols that can be used to understand and predict macromolecular recognition, and to gain more structural insight into three newly discovered allosteric binding sites on a critical drug target involved in the treatment of HIV infections, the Olson and Levy labs collaborated on the SAMPL4 challenge. This computational blind challenge involved predicting protein-ligand binding against the three allosteric sites of HIV integrase (IN), a viral enzyme for which two drugs (that target the active site) have been approved by the FDA. Positive control cross-docking experiments were utilized to select 13 receptor models out of an initial ensemble of 41 different crystal structures of HIV IN. These 13 models of the targets were selected using our new "Rank Difference Ratio" metric. The first stage of SAMPL4 involved using virtual screens to identify 62 active, allosteric IN inhibitors out of a set of 321 compounds. The second stage involved predicting the binding site(s) and crystallographic binding mode(s) for 57 of these inhibitors. Our team submitted four entries for the first stage that utilized: (1) AutoDock Vina (AD Vina) plus visual inspection; (2) a new common pharmacophore engine; (3) BEDAM replica exchange free energy simulations, and a Consensus approach that combined the predictions of all three strategies. Even with the SAMPL4's very challenging compound library that displayed a significantly lower amount of structural diversity than most libraries that are conventionally employed in prospective virtual screens, these approaches produced hit rates of 24, 25, 34, and 27 %, respectively, on a set with 19 % declared binders. Our only entry for the second stage challenge was based on the results of AD Vina plus visual inspection, and it ranked third place overall according to several different metrics provided by the SAMPL4 organizers. The successful results displayed by these approaches highlight the utility of the computational structure-based drug discovery tools and strategies that are being developed to advance the goals of the newly created, multi-institution, NIH-funded center called the "HIV Interaction and Viral Evolution Center".