Using physics-based pose predictions and free energy perturbation calculations to predict binding poses and relative binding affinities for FXR ligands in the D3R Grand Challenge 2

De novo drug design through gradient-based regularized search in information-theoretically controlled latent space

Over the last decade, automatic chemical design frameworks for discovering molecules with drug-like properties have significantly progressed. Among them, the variational autoencoder (VAE) is a cutting-edge approach that models the tractable latent space of the molecular space. In particular, the usage of a VAE along with a property estimator has attracted considerable interest because it enables gradient-based optimization of a given molecule. However, although successful results have been achieved experimentally, the theoretical background and prerequisites for the correct operation of this method have not yet been clarified. In view of the above, we theoretically analyze and rigorously reconstruct the entire framework. From the perspective of parameterized distribution and the information theory, we first describe how the previous model overcomes the limitations of the beta VAE in discovering molecules with the desired properties. Furthermore, we describe the prerequisites for training the above model. Next, from the log-likelihood perspective of each term, we reformulate the objectives for exploring latent space to generate drug-like molecules. The distributional constraints are defined in this study, which will break away from the invalid molecular search. We demonstrated that our model could discover a novel chemical compound for targeting BCL-2 family proteins in de novo approach. Through the theoretical analysis and practical implementation, the importance of the aforementioned prerequisites and constraints to operate the model was verified.

Computational Investigation of BMAA and Its Carbamate Adducts as Potential GluR2 Modulators

Beta-N-methylamino-l-alanine (BMAA) is a potential neurotoxic nonprotein amino acid, which can reach the human body through the food chain. When BMAA interacts with bicarbonate in the human body, carbamate adducts are produced, which share a high structural similarity with the neurotransmitter glutamate. It is believed that BMAA and its l-carbamate adducts bind in the glutamate binding site of ionotropic glutamate receptor 2 (GluR2). Chronic exposure to BMAA and its adducts could cause neurological illness such as neurodegenerative diseases. However, the mechanism of BMAA action and its carbamate adducts bound to GluR2 has not yet been elucidated. Here, we investigate the binding modes and the affinity of BMAA and its carbamate adducts to GluR2 in comparison to the natural agonist, glutamate, to understand whether these can act as GluR2 modulators. Initially, we perform molecular dynamics simulations of BMAA and its carbamate adducts bound to GluR2 to examine the stability of the ligands in the S1/S2 ligand-binding core of the receptor. In addition, we utilize alchemical free energy calculations to compute the difference in the free energy of binding of the beta-carbamate adduct of BMAA to GluR2 compared to that of glutamate. Our findings indicate that carbamate adducts of BMAA and glutamate remain stable in the binding site of the GluR2 compared to BMAA. Additionally, alchemical free energy results reveal that glutamate and the beta-carbamate adduct of BMAA have comparable binding affinity to the GluR2. These results provide a rationale that BMAA carbamate adducts may be, in fact, the modulators of GluR2 and not BMAA itself.

Modern Alchemical Free Energy Methods for Drug Discovery Explained

This Perspective provides a contextual explanation of the current state-of-the-art alchemical free energy methods and their role in drug discovery as well as highlights select emerging technologies. The narrative attempts to answer basic questions about what goes on "under the hood" in free energy simulations and provide general guidelines for how to run simulations and analyze the results. It is the hope that this work will provide a valuable introduction to students and scientists in the field.

Assessing the performance of docking, FEP, and MM/GBSA methods on a series of KLK6 inhibitors

Kallikrein 6 (KLK6) is an attractive drug target for the treatment of neurological diseases and for various cancers. Herein, we explore the accuracy and efficiency of different computational methods and protocols to predict the free energy of binding (ΔG_bind) for a series of 49 inhibitors of KLK6. We found that the performance of the methods varied strongly with the tested system. For only one of the three KLK6 datasets, the docking scores obtained with rDock were in good agreement (R² ≥ 0.5) with experimental values of ΔG_bind. A similar result was obtained with MM/GBSA (using the ff14SB force field) calculations based on single minimized structures. Improved binding affinity predictions were obtained with the free energy perturbation (FEP) method, with an overall MUE and RMSE of 0.53 and 0.68 kcal/mol, respectively. Furthermore, in a simulation of a real-world drug discovery project, FEP was able to rank the most potent compounds at the top of the list. These results indicate that FEP can be a promising tool for the structure-based optimization of KLK6 inhibitors.

Accuracy and precision of binding free energy prediction for a tacrine related lead inhibitor of acetylcholinesterase with an arsenal of supercomputerized molecular modelling methods: a comparative study

Nowadays, advanced computational chemistry methods offer various strategies for revealing prospective hit structures in drug development essentially through accurate binding free energy predictions. After the era of molecular docking and quantitative structure-activity relationships, much interest has been lately oriented to perturbed molecular dynamic approaches like replica exchange with solute tempering and free energy perturbation (REST/FEP) and the potential of the mean force with adaptive biasing and accelerated weight histograms (PMF/AWH). Both of these receptor-based techniques can exploit exascale CPU&GPU supercomputers to achieve high throughput performance. In this fundamental study, we have compared the predictive power of a panel of supercomputerized molecular modelling methods to distinguish the major binding modes and the corresponding binding free energies of a promising tacrine related potential antialzheimerics in human acetylcholinesterase. The binding free energies were estimated using flexible molecular docking, molecular mechanics/generalized Born surface area/Poisson-Boltzmann surface area (MM/GBSA/PBSA), transmutation REST/FEP with 12 x 5 ns/λ windows, annihilation FEP with 20 x 5 ns/λ steps, PMF with weight histogram analysis method (WHAM) and 40 x 5 ns samples, and PMF/AWH with 10 x 100 ns replicas. Confrontation of the classical approaches such as canonical molecular dynamics and molecular docking with alchemical calculations and steered molecular dynamics enabled us to show how large errors in ΔG predictions can be expected if these in silico methods are employed in the elucidation of a common case of enzyme inhibition.Communicated by Ramaswamy H. Sarma.

Ligand and Structure-Based Virtual Screening in Combination, to Evaluate Small Organic Molecules as Inhibitors for the XIAP Anti-Apoptotic Protein: The Xanthohumol Hypothesis

Herein, we propose two chalcone molecules, (E)-1-(4-methoxyphenyl)-3-(p-tolyl) prop-2-en-1-one and (E)-3-(4-hydroxyphenyl)-1-(2,4,6-trihydroxyphenyl) prop-2-en-1-one, based on the anticancer bioactive molecule Xanthohumol, which are suitable for further in vitro and in vivo studies. Their ability to create stable complexes with the antiapoptotic X-linked IAP (XIAP) protein makes them promising anticancer agents. The calculations were based on ligand-based and structure-based virtual screening combined with the pharmacophore build. Additionally, the structures passed Lipinski’s rule for drug use, and their reactivity was confirmed using density functional theory studies. ADMET studies were also performed to reveal the pharmacokinetic potential of the compounds. The candidates were chosen from 10,639,400 compounds, and the docking protocols were evaluated using molecular dynamics simulations.

FEPrepare: A Web-Based Tool for Automating the Setup of Relative Binding Free Energy Calculations

Relative binding free energy calculations in drug design are becoming a useful tool in facilitating lead binding affinity optimization in a cost- and time-efficient manner. However, they have been limited by technical challenges such as the manual creation of large numbers of input files to set up, run, and analyze free energy simulations. In this Application Note, we describe FEPrepare, a novel web-based tool, which automates the setup procedure for relative binding FEP calculations for the dual-topology scheme of NAMD, one of the major MD engines, using OPLS-AA force field topology and parameter files. FEPrepare provides the user with all necessary files needed to run a FEP/MD simulation with NAMD. FEPrepare can be accessed and used at https://feprepare.vi-seem.eu/.

Thermodynamics and Kinetics of Drug-Target Binding by Molecular Simulation

Computational studies play an increasingly important role in chemistry and biophysics, mainly thanks to improvements in hardware and algorithms. In drug discovery and development, computational studies can reduce the costs and risks of bringing a new medicine to market. Computational simulations are mainly used to optimize promising new compounds by estimating their binding affinity to proteins. This is challenging due to the complexity of the simulated system. To assess the present and future value of simulation for drug discovery, we review key applications of advanced methods for sampling complex free-energy landscapes at near nonergodicity conditions and for estimating the rate coefficients of very slow processes of pharmacological interest. We outline the statistical mechanics and computational background behind this research, including methods such as steered molecular dynamics and metadynamics. We review recent applications to pharmacology and drug discovery and discuss possible guidelines for the practitioner. Recent trends in machine learning are also briefly discussed. Thanks to the rapid development of methods for characterizing and quantifying rare events, simulation's role in drug discovery is likely to expand, making it a valuable complement to experimental and clinical approaches.

Machine learning classification can reduce false positives in structure-based virtual screening

With the recent explosion in the size of libraries available for screening, virtual screening is positioned to assume a more prominent role in early drug discovery's search for active chemical matter. In typical virtual screens, however, only about 12% of the top-scoring compounds actually show activity when tested in biochemical assays. We argue that most scoring functions used for this task have been developed with insufficient thoughtfulness into the datasets on which they are trained and tested, leading to overly simplistic models and/or overtraining. These problems are compounded in the literature because studies reporting new scoring methods have not validated their models prospectively within the same study. Here, we report a strategy for building a training dataset (D-COID) that aims to generate highly compelling decoy complexes that are individually matched to available active complexes. Using this dataset, we train a general-purpose classifier for virtual screening (vScreenML) that is built on the XGBoost framework. In retrospective benchmarks, our classifier shows outstanding performance relative to other scoring functions. In a prospective context, nearly all candidate inhibitors from a screen against acetylcholinesterase show detectable activity; beyond this, 10 of 23 compounds have IC₅₀ better than 50 μM. Without any medicinal chemistry optimization, the most potent hit has IC₅₀ 280 nM, corresponding to K_i of 173 nM. These results support using the D-COID strategy for training classifiers in other computational biology tasks, and for vScreenML in virtual screening campaigns against other protein targets. Both D-COID and vScreenML are freely distributed to facilitate such efforts.

Rigorous Free Energy Simulations in Virtual Screening

Virtual high throughput screening (vHTS) in drug discovery is a powerful approach to identify hits: when applied successfully, it can be much faster and cheaper than experimental high-throughput screening approaches. However, mainstream vHTS tools have significant limitations: ligand-based methods depend on knowledge of existing chemical matter, while structure-based tools such as docking involve significant approximations that limit their accuracy. Recent advances in scientific methods coupled with dramatic speedups in computational processing with GPUs make this an opportune time to consider the role of more rigorous methods that could improve the predictive power of vHTS workflows. In this Perspective, we assert that alchemical binding free energy methods using all-atom molecular dynamics simulations have matured to the point where they can be applied in virtual screening campaigns as a final scoring stage to prioritize the top molecules for experimental testing. Specifically, we propose that alchemical absolute binding free energy (ABFE) calculations offer the most direct and computationally efficient approach within a rigorous statistical thermodynamic framework for computing binding energies of diverse molecules, as is required for virtual screening. ABFE calculations are particularly attractive for drug discovery at this point in time, where the confluence of large-scale genomics data and insights from chemical biology have unveiled a large number of promising disease targets for which no small molecule binders are known, precluding ligand-based approaches, and where traditional docking approaches have foundered to find progressible chemical matter.

D3R grand challenge 4: blind prediction of protein-ligand poses, affinity rankings, and relative binding free energies

The Drug Design Data Resource (D3R) aims to identify best practice methods for computer aided drug design through blinded ligand pose prediction and affinity challenges. Herein, we report on the results of Grand Challenge 4 (GC4). GC4 focused on proteins beta secretase 1 and Cathepsin S, and was run in an analogous manner to prior challenges. In Stage 1, participant ability to predict the pose and affinity of BACE1 ligands were assessed. Following the completion of Stage 1, all BACE1 co-crystal structures were released, and Stage 2 tested affinity rankings with co-crystal structures. We provide an analysis of the results and discuss insights into determined best practice methods.

Free-energy landscape of molecular interactions between endothelin 1 and human endothelin type B receptor: fly-casting mechanism

The free-energy landscape of interaction between a medium-sized peptide, endothelin 1 (ET1), and its receptor, human endothelin type B receptor (hETB), was computed using multidimensional virtual-system coupled molecular dynamics, which controls the system's motions by introducing multiple reaction coordinates. The hETB embedded in lipid bilayer was immersed in explicit solvent. All molecules were expressed as all-atom models. The resultant free-energy landscape had five ranges with decreasing ET1-hETB distance: completely dissociative, outside-gate, gate, binding pocket, and genuine-bound ranges. In the completely dissociative range, no ET1-hETB interaction appeared. In the outside-gate range, an ET1-hETB attractive interaction was the fly-casting mechanism. In the gate range, the ET1 orientational variety decreased rapidly. In the binding pocket range, ET1 was in a narrow pathway with a steep free-energy slope. In the genuine-bound range, ET1 was in a stable free-energy basin. A G-protein-coupled receptor (GPCR) might capture its ligand from a distant place.

Molecular Mechanism for Attractant Signaling to DHMA by E. coli Tsr

The attractant chemotaxis response of Escherichia coli to norepinephrine requires that it be converted to 3,4-dihydroxymandelic acid (DHMA) by the monoamine oxidase TynA and the aromatic aldehyde dehydrogenase FeaB. DHMA is sensed by the serine chemoreceptor Tsr, and the attractant response requires that at least one subunit of the periplasmic domain of the Tsr homodimer (pTsr) has an intact serine-binding site. DHMA that is generated in vivo by E. coli is expected to be a racemic mixture of the (R) and (S) enantiomers, so it has been unclear whether one or both chiral forms are active. Here, we used a combination of state-of-the-art tools in molecular docking and simulations, including an in-house simulation-based docking protocol, to investigate the binding properties of (R)-DHMA and (S)-DHMA to E. coli pTsr. Our studies computationally predicted that (R)-DHMA should promote a stronger attractant response than (S)-DHMA because of a consistently greater-magnitude piston-like pushdown of the pTsr α-helix 4 toward the membrane upon binding of (R)-DHMA than upon binding of (S)-DHMA. This displacement is caused primarily by interaction of DHMA with Tsr residue Thr156, which has been shown by genetic studies to be critical for the attractant response to L-serine and DHMA. These findings led us to separate the two chiral species and test their effectiveness as chemoattractants. Both the tethered cell and motility migration coefficient assays validated the prediction that (R)-DHMA is a stronger attractant than (S)-DHMA. Our study demonstrates that refined computational docking and simulation studies combined with experiments can be used to investigate situations in which subtle differences between ligands may lead to diverse chemotactic responses.

Effect of set up protocols on the accuracy of alchemical free energy calculation over a set of ACK1 inhibitors

Hit-to-lead virtual screening frequently relies on a cascade of computational methods that starts with rapid calculations applied to a large number of compounds and ends with more expensive computations restricted to a subset of compounds that passed initial filters. This work focuses on set up protocols for alchemical free energy (AFE) scoring in the context of a Docking–MM/PBSA–AFE cascade. A dataset of 15 congeneric inhibitors of the ACK1 protein was used to evaluate the performance of AFE set up protocols that varied in the steps taken to prepare input files (using previously docked and best scored poses, manual selection of poses, manual placement of binding site water molecules). The main finding is that use of knowledge derived from X-ray structures to model binding modes, together with the manual placement of a bridging water molecule, improves the R² from 0.45 ± 0.06 to 0.76 ± 0.02 and decreases the mean unsigned error from 2.11 ± 0.08 to 1.24 ± 0.04 kcal mol^-1. By contrast a brute force automated protocol that increased the sampling time ten-fold lead to little improvements in accuracy. Besides, it is shown that for the present dataset hysteresis can be used to flag poses that need further attention even without prior knowledge of experimental binding affinities.

Inhibiting two cellular mutant epidermal growth factor receptor tyrosine kinases by addressing computationally assessed crystal ligand pockets

Aim: Blocking receptor tyrosine kinases is a useful strategy for inhibiting the overexpression of EGFR. However, the quality of crystal pocket is an essential issue for virtually identifying new leads for surviving resistance cancer cells. Results: With the examinating crystal pocket quality by the self-docking root-mean-square deviation (RMSD) calculation, we used the two best kinase pockets of mutant EGFR kinases, T790M/L858R and G719S, for virtual screening. After sorting all the docking poses of the 57,177 library compounds by consensus scores, three evidently blocked cellular EGFR phosphorylation in the H1975 and SW48 cell lines. Conclusion: The computationally assessed qualities of crystal pockets of crystal EGFR kinases can help identify new cellular active and target-specific ligands rapidly and at low cost.

Modeling the value of predictive affinity scoring in preclinical drug discovery

Drug discovery is widely recognized to be a difficult and costly activity in large part due to the challenge of identifying chemical matter which simultaneously optimizes multiple properties, one of which is affinity for the primary biological target. Further, many of these properties are difficult to predict ahead of expensive and time-consuming compound synthesis and experimental testing. Here we highlight recent work to develop compound affinity prediction models, and extensively investigate the value such models may provide to preclinical drug discovery. We demonstrate that the ability of these models to improve the overall probability of success is crucially dependent on the shape of the error distribution, not just the root-mean-square error. In particular, while scoring more molecule ideas generally improves the probability of project success when the error distribution is Gaussian, fat-tail distributions such as a Cauchy distribution, can lead to a situation where scoring more ideas actually decreases the overall probability of success.

Relative Binding Free Energy Calculations in Drug Discovery: Recent Advances and Practical Considerations

Accurate in silico prediction of protein-ligand binding affinities has been a primary objective of structure-based drug design for decades due to the putative value it would bring to the drug discovery process. However, computational methods have historically failed to deliver value in real-world drug discovery applications due to a variety of scientific, technical, and practical challenges. Recently, a family of approaches commonly referred to as relative binding free energy (RBFE) calculations, which rely on physics-based molecular simulations and statistical mechanics, have shown promise in reliably generating accurate predictions in the context of drug discovery projects. This advance arises from accumulating developments in the underlying scientific methods (decades of research on force fields and sampling algorithms) coupled with vast increases in computational resources (graphics processing units and cloud infrastructures). Mounting evidence from retrospective validation studies, blind challenge predictions, and prospective applications suggests that RBFE simulations can now predict the affinity differences for congeneric ligands with sufficient accuracy and throughput to deliver considerable value in hit-to-lead and lead optimization efforts. Here, we present an overview of current RBFE implementations, highlighting recent advances and remaining challenges, along with examples that emphasize practical considerations for obtaining reliable RBFE results. We focus specifically on relative binding free energies because the calculations are less computationally intensive than absolute binding free energy (ABFE) calculations and map directly onto the hit-to-lead and lead optimization processes, where the prediction of relative binding energies between a reference molecule and new ideas (virtual molecules) can be used to prioritize molecules for synthesis. We describe the critical aspects of running RBFE calculations, from both theoretical and applied perspectives, using a combination of retrospective literature examples and prospective studies from drug discovery projects. This work is intended to provide a contemporary overview of the scientific, technical, and practical issues associated with running relative binding free energy simulations, with a focus on real-world drug discovery applications. We offer guidelines for improving the accuracy of RBFE simulations, especially for challenging cases, and emphasize unresolved issues that could be improved by further research in the field.