Identification of ligands for RNA targets via structure-based virtual screening: HIV-1 TAR

Binding of the Tat protein to TAR RNA is necessary for viral replication of HIV-1. We screened the Available Chemicals Directory (ACD) to identify ligands to bind to a TAR RNA structure using a four-step docking procedure: rigid docking first, followed by three steps of flexible docking using a pseudobrownian Monte Carlo minimization in torsion angle space with progressively more detailed conformational sampling on a progressively smaller list of top-ranking compounds. To validate the procedure, we successfully docked ligands for five RNA complexes of known structure. For ranking ligands according to binding avidity, an empirical binding free energy function was developed which accounts, in particular, for solvation, isomerization free energy, and changes in conformational entropy. System-specific parameters for the function were derived on a training set of RNA/ligand complexes with known structure and affinity. To validate the free energy function, we screened the entire ACD for ligands for an RNA aptamer which binds l-arginine tightly. The native ligand ranked 17 out of ca. 153,000 compounds screened, i.e., the procedure is able to filter out >99.98% of the database and still retain the native ligand. Screening of the ACD for TAR ligands yielded a high rank for all known TAR ligands contained in the ACD and suggested several other potential TAR ligands. Eight of the highest ranking compounds not previously known to be ligands were assayed for inhibition of the Tat-TAR interaction, and two exhibited a CD₅₀ of ca. 1 M.     

Exploring hierarchical refinement techniques for induced fit docking with protein and ligand flexibility

We present a series of molecular‐mechanics‐based protein refinement methods, including two novel ones, applied as part of an induced fit docking procedure. The methods used include minimization; protein and ligand sidechain prediction; a hierarchical ligand placement procedure similar to a‐priori protein loop predictions; and a minimized Monte Carlo approach using normal mode analysis as a move step. The results clearly indicate the importance of a proper opening of the active site backbone, which might not be accomplished when the ligand degrees of freedom are prioritized. The most accurate method consisted of the minimized Monte Carlo procedure designed to open the active site followed by a hierarchical optimization of the sidechain packing around a mobile flexible ligand. The methods have been used on a series of 88 protein‐ligand complexes including both cross‐docking and apo‐docking members resulting in complex conformations determined to within 2.0 Å heavy‐atom RMSD in 75% of cases where the protein backbone rearrangement upon binding is less than 1.0 Å α‐carbon RMSD. We also demonstrate that physics‐based all‐atom potentials can be more accurate than docking‐style potentials when complexes are sufficiently refined. © 2009 Wiley Periodicals, Inc. J Comput Chem, 2010     

FDS: flexible ligand and receptor docking with a continuum solvent model and soft-core energy function

The docking of flexible small molecule ligands to large flexible protein targets is addressed in this article using a two-stage simulation-based method. The methodology presented is a hybrid approach where the first component is a dock of the ligand to the protein binding site, based on deriving sets of simultaneously satisfied intermolecular hydrogen bonds using graph theory and a recursive distance geometry algorithm. The output structures are reduced in number by cluster analysis based on distance similarities. These structures are submitted to a modified Monte Carlo algorithm using the AMBER-AA molecular mechanics force field with the Generalized Born/Surface Area (GB/SA) continuum model. This solvent model is not only less expensive than an explicit representation, but also yields increased sampling. Sampling is also increased using a rotamer library to direct some of the protein side-chain movements along with large dihedral moves. Finally, a softening function for the nonbonded force field terms is used, enabling the potential energy function to be slowly turned on throughout the course of the simulation. The docking procedure is optimized, and the results are presented for a single complex of the arabinose binding protein. It was found that for a rigid receptor model, the X-ray binding geometry was reproduced and uniquely identified based on the associated potential energy. However, when side-chain flexibility was included, although the X-ray structure was identified, it was one of three possible binding geometries that were energetically indistinguishable. These results suggest that on relaxing the constraint on receptor flexibility, the docking energy hypersurface changes from being funnel-like to rugged. A further 14 complexes were then examined using the optimized protocol. For each complex the docking methodology was tested for a fully flexible ligand, both with and without protein side-chain flexibility. For the rigid protein docking, 13 out of the 15 test cases were able to find the experimental binding mode; this number was reduced to 11 for the flexible protein docking. However, of these 11, in the majority of cases the experimental binding mode was not uniquely identified, but was present in a cluster of low energy structures that were energetically indistinguishable. These results not only support the presence of a rugged docking energy hypersurface, but also suggest that it may be necessary to consider the possibility of more than one binding conformation during ligand optimization.     

Fully automated molecular mechanics based induced fit protein-ligand docking method

We describe a method for docking a ligand into a protein receptor while allowing flexibility of the protein binding site. The method employs a multistep procedure that begins with the generation of protein and ligand conformations. An initial placement of the ligand is then performed by computing binding site hotspots. This initial placement is followed by a protein side-chain refinement stage that models protein flexibility. The final step of the process is an energy minimization of the ligand pose in the presence of the rigid receptor. Thus the algorithm models flexibility of the protein at two stages, before and after ligand placement. We validated this method by performing docking and cross docking studies of eight protein systems for which crystal structures were available for at least two bound ligands. The resulting rmsd values of the 21 docked protein-ligand complexes showed values of 2 A or less for all but one of the systems examined. The method has two critical benefits for high throughput virtual screening studies. First, no user intervention is required in the docking once the initial binding site selection has been made in the protein. Second, the initial protein conformation generation needs to be performed only once for a given binding region. Also, the method may be customized in various ways depending on the particular scenario in which dockings are being performed. Each of the individual steps of the method is fully independent making it straightforward to explore different variants of the high level workflow to further improve accuracy and performance.     

Automated docking using a Lamarckian genetic algorithm and an empirical binding free energy function

A novel and robust automated docking method that predicts the bound conformations of flexible ligands to macromolecular targets has been developed and tested, in combination with a new scoring function that estimates the free energy change upon binding. Interestingly, this method applies a Lamarckian model of genetics, in which environmental adaptations of an individual's phenotype are reverse transcribed into its genotype and become heritable traits (sic). We consider three search methods, Monte Carlo simulated annealing, a traditional genetic algorithm, and the Lamarckian genetic algorithm, and compare their performance in dockings of seven protein–ligand test systems having known three-dimensional structure. We show that both the traditional and Lamarckian genetic algorithms can handle ligands with more degrees of freedom than the simulated annealing method used in earlier versions of AUTO DOCK , and that the Lamarckian genetic algorithm is the most efficient, reliable, and successful of the three. The empirical free energy function was calibrated using a set of 30 structurally known protein–ligand complexes with experimentally determined binding constants. Linear regression analysis of the observed binding constants in terms of a wide variety of structure-derived molecular properties was performed. The final model had a residual standard error of 9.11 kJ mol⁻¹ (2.177 kcal mol⁻¹) and was chosen as the new energy function. The new search methods and empirical free energy function are available in AUTO DOCK , version 3.0. © 1998 John Wiley & Sons, Inc. J Comput Chem 19: 1639–1662, 1998     

A QXP-based multistep docking procedure for accurate prediction of protein-ligand complexes

The two great challenges of the docking process are the prediction of ligand poses in a protein binding site and the scoring of the docked poses. Ligands that are composed of extended chains in their molecular structure display the most difficulties, predominantly because of the torsional flexibility. On the basis of the molecular docking program QXP-Flo+0802, we have developed a procedure particularly for ligands with a high degree of rotational freedom that allows the accurate prediction of the orientation and conformation of ligands in protein binding sites. Starting from an initial full Monte Carlo docking experiment, this was achieved by performing a series of successive multistep docking runs using a local Monte Carlo search with a restricted rotational angle, by which the conformational search space is limited. The method was established by using a highly flexible acetylcholinesterase inhibitor and has been applied to a number of challenging protein-ligand complexes known from the literature.     

Automated docking of ligands to an artificial active site: augmenting crystallographic analysis with computer modeling

The W191G cavity of cytochrome c peroxidase is useful as a model system for introducing small molecule oxidation in an artificially created cavity. A set of small, cyclic, organic cations was previously shown to bind in the buried, solvent-filled pocket created by the W191G mutation. We docked these ligands and a set of non-binders in the W191G cavity using AutoDock 3.0. For the ligands, we compared docking predictions with experimentally determined binding energies and X-ray crystal structure complexes. For the ligands, predicted binding energies differed from measured values by +/- 0.8 kcal/mol. For most ligands, the docking simulation clearly predicted a single binding mode that matched the crystallographic binding mode within 1.0 A RMSD. For 2 ligands, where the docking procedure yielded an ambiguous result, solutions matching the crystallographic result could be obtained by including an additional crystallographically observed water molecule in the protein model. For the remaining 2 ligands, docking indicated multiple binding modes, consistent with the original electron density, suggesting disordered binding of these ligands. Visual inspection of the atomic affinity grid maps used in docking calculations revealed two patches of high affinity for hydrogen bond donating groups. Multiple solutions are predicted as these two sites compete for polar hydrogens in the ligand during the docking simulation. Ligands could be distinguished, to some extent, from non-binders using a combination of two trends: predicted binding energy and level of clustering. In summary, AutoDock 3.0 appears to be useful in predicting key structural and energetic features of ligand binding in the W191G cavity.     

Multiple ligand simultaneous docking: Orchestrated dancing of ligands in binding sites of protein

Present docking methodologies simulate only one single ligand at a time during docking process. In reality, the molecular recognition process always involves multiple molecular species. Typical protein–ligand interactions are, for example, substrate and cofactor in catalytic cycle; metal ion coordination together with ligand(s); and ligand binding with water molecules. To simulate the real molecular binding processes, we propose a novel multiple ligand simultaneous docking (MLSD) strategy, which can deal with all the above processes, vastly improving docking sampling and binding free energy scoring. The work also compares two search strategies: Lamarckian genetic algorithm and particle swarm optimization, which have respective advantages depending on the specific systems. The methodology proves robust through systematic testing against several diverse model systems: E. coli purine nucleoside phosphorylase (PNP) complex with two substrates, SHP2NSH2 complex with two peptides and Bcl‐xL complex with ABT‐737 fragments. In all cases, the final correct docking poses and relative binding free energies were obtained. In PNP case, the simulations also capture the binding intermediates and reveal the binding dynamics during the recognition processes, which are consistent with the proposed enzymatic mechanism. In the other two cases, conventional single‐ligand docking fails due to energetic and dynamic coupling among ligands, whereas MLSD results in the correct binding modes. These three cases also represent potential applications in the areas of exploring enzymatic mechanism, interpreting noisy X‐ray crystallographic maps, and aiding fragment‐based drug design, respectively. © 2010 Wiley Periodicals, Inc. J Comput Chem, 2010     

Hammerhead: fast,fully automated docking of flexible ligands to protein binding sites

Background: Molecular docking seeks to predict the geometry and affinity of the binding of a small molecule to a given protein of known structure. Rigid docking has long been used to screen databases of small molecules, because docking techniques that account for ligand flexibility have either been too slow or have required significant human intervention. Here we describe a docking algorithm, Hammerhead, which is a fast, automated tool to screen for the binding of flexible molecules to protein binding sites.Results: We used Hammerhead to successfully dock a variety of positive control ligands into their cognate proteins. The empirically tuned scoring function of the algorithm predicted binding affinities within 1.3 log units of the known affinities for these ligands. Conformations and alignments close to those determined crystallographically received the highest scores. We screened 80 000 compounds for binding to streptavidin, and biotin was predicted as the top-scoring ligand, with other known ligands included among the highest-scoring dockings. The screen ran in a few days on commonly available hardware.Conclusions: Hammerhead is suitable for screening large databases of flexible molecules for binding to a protein of known structure. It correctly docks a variety of known flexible ligands, and it spends an average of only a few seconds on each compound during a screen. The approach is completely automated, from the elucidation of protein binding sites, through the docking of molecules, to the final selection of compounds for assay.     

BALLDock/SLICK: a new method for protein-carbohydrate docking

Protein-ligand docking is an essential technique in computer-aided drug design. While generally available docking programs work well for most drug classes, carbohydrates and carbohydrate-like compounds are often problematic for docking. We present a new docking method specifically designed to handle docking of carbohydrate-like compounds. BALLDock/SLICK combines an evolutionary docking algorithm for flexible ligands and flexible receptor side chains with carbohydrate-specific scoring and energy functions. The scoring function has been designed to identify accurate ligand poses, while the energy function yields accurate estimates of the binding free energies of these poses. On a test set of known protein-sugar complexes we demonstrate the ability of the approach to generate correct poses for almost all of the structures and achieve very low mean errors for the predicted binding free energies.     

A molecular docking study of potential inhibitors and repurposed drugs against SARS-CoV-2 main protease enzyme

COVID-19 has affected millions of people. Although many drugs are in use to combat disease, there is not any sufficient treatment yet. Having critical role in propagation of the novel coronavirus (SARS-CoV-2) works Main Protease up into a significant drug target. We have performed a molecular docking study to define possible inhibitor candidates against SARS-CoV-2 Main Protease enzyme. Besides docking Remdesivir, Ribavirin, Chloroquine and 28 other antiviral inhibitors (totally 31 inhibitors) to Main Protease enzyme, we have also performed a molecular docking study of 2177 ligands, which are used against Main Protease for the first time by using molecular docking program Autodock4. All ligands were successfully docked into Main Protease enzyme binding site. Among all ligands, EY16 coded ligand which previously used as EBNA1-DNA binding blocker candidate showed the best score for Main Protease with a binding free energy of −10.83 ​kcal/mol which was also lower than re-docking score of N3 ligand (−10.72 ​kcal/mol) contained in crystal structure of Main Protease. After analyzing the docking modes and docking scores we have found that our ligands have better binding free energy values than the inhibitors in use of treatment. We believe that further studies such as molecular dynamics or Molecular Mechanic Poisson Boltzmann Surface Area studies can make contribution that is more exhaustive to the docking results.     

Molecular modelling prediction of ligand binding site flexibility

We have investigated the efficacy of generating multiple sidechain conformations using a rotamer library in order to find the experimentally observed ligand binding site conformation of a protein in the presence of a bound ligand. We made use of a recently published algorithm that performs an exhaustive conformational search using a rotamer library to enumerate all possible sidechain conformations in a binding site. This approach was applied to a dataset of proteins whose structures were determined by X-ray and NMR methods. All chosen proteins had two or more structures, generally involving different bound ligands. By taking one of these structures as a reference, we were able in most cases to successfully reproduce the experimentally determined conformations of the other structures, as well as to suggest alternative low-energy conformations of the binding site. In those few cases where this procedure failed, we observed that the bound ligand had induced a high-energy conformation of the binding site. These results suggest that for most proteins that exhibit limited backbone motion, ligands tend to bind to low energy conformations of their binding sites. Our results also reveal that it is possible in most cases to use a rotamer search-based approach to predict alternative low-energy protein binding site conformations that can be used by different ligands. This opens the possibility of incorporating alternative binding site conformations to improve the efficacy of docking and structure-based drug design algorithms.     

Computational chemistry in drug lead discovery and design

The main contributions of our group during the last 15 years developing and using biomolecular simulation tools in drug lead discovery and design, in close collaboration with experimental researchers, are presented. Special emphasis has been given to methodological improvements in the following areas: (1) target homology modeling incorporating knowledge about known ligands to accurately characterize the binding site; (2) designing alternative strategies to account for protein flexibility in high-throughput docking; (3) development of stochastic- and normal-mode-based methods to de novo design structurally diverse protein conformers; (4) development and validation of quantum mechanical semi-empirical linear-scaling calculations to correctly estimate ligand binding free energy. Several successful cases of computer-aided drug discovery are also presented, especially our recent work on viral targets.     

QUASI: a novel method for simultaneous superposition of multiple flexible ligands and virtual screening using partial similarity

The structure of many receptors is unknown, and only information about diverse ligands binding to them is available. A new method is presented for the superposition of such ligands, derivation of putative receptor site models and utilization of the models for screening of compound databases. In order to generate a receptor model, the similarity of all ligands is optimized simultaneously taking into account conformational flexibility and also the possibility that the ligands can bind to different regions of the site and only partially overlap. Ligand similarity is defined with respect to a receptor site model serving as a common reference frame. The receptor model is dynamic and coevolves with the ligand alignment until an optimal self-consistent superposition is achieved. When ligand conformational flexibility is permitted, different superposition models are possible and consistent with the data. Clustering of the superposition solutions is used to obtain diverse models. When the models are used to screen a database of compounds, high enrichments are obtained, comparable to those obtained in docking studies.     

Protein:Ligand binding free energies: A stringent test for computational protein design

A computational protein design method is extended to allow Monte Carlo simulations where two ligands are titrated into a protein binding pocket, yielding binding free energy differences. These provide a stringent test of the physical model, including the energy surface and sidechain rotamer definition. As a test, we consider tyrosyl‐tRNA synthetase (TyrRS), which has been extensively redesigned experimentally. We consider its specificity for its substrate l ‐tyrosine (l ‐Tyr), compared to the analogs d ‐Tyr, p‐acetyl‐, and p‐azido‐phenylalanine (ac‐Phe, az‐Phe). We simulate l ‐ and d ‐Tyr binding to TyrRS and six mutants, and compare the structures and binding free energies to a more rigorous “MD/GBSA” procedure: molecular dynamics with explicit solvent for structures and a Generalized Born + Surface Area model for binding free energies. Next, we consider l ‐Tyr, ac‐ and az‐Phe binding to six other TyrRS variants. The titration results are sensitive to the precise rotamer definition, which involves a short energy minimization for each sidechain pair to help relax bad contacts induced by the discrete rotamer set. However, when designed mutant structures are rescored with a standard GBSA energy model, results agree well with the more rigorous MD/GBSA. As a third test, we redesign three amino acid positions in the substrate coordination sphere, with either l ‐Tyr or d ‐Tyr as the ligand. For two, we obtain good agreement with experiment, recovering the wildtype residue when l ‐Tyr is the ligand and a d ‐Tyr specific mutant when d ‐Tyr is the ligand. For the third, we recover His with either ligand, instead of wildtype Gln. © 2015 Wiley Periodicals, Inc.     

Assessment of programs for ligand binding affinity prediction

The prediction of the binding free energy between a ligand and a protein is an important component in the virtual screening and lead optimization of ligands for drug discovery. To determine the quality of current binding free energy estimation programs, we examined FlexX, X-Score, AutoDock, and BLEEP for their performance in binding free energy prediction in various situations including cocrystallized complex structures, cross docking of ligands to their non-cocrystallized receptors, docking of thermally unfolded receptor decoys to their ligands, and complex structures with "randomized" ligand decoys. In no case was there a satisfactory correlation between the experimental and estimated binding free energies over all the datasets tested. Meanwhile, a strong correlation between ligand molecular weight-binding affinity correlation and experimental predicted binding affinity correlation was found. Sometimes the programs also correctly ranked ligands' binding affinities even though native interactions between the ligands and their receptors were essentially lost because of receptor deformation or ligand randomization, and the programs could not decisively discriminate randomized ligand decoys from their native ligands; this suggested that the tested programs miss important components for the accurate capture of specific ligand binding interactions.     

Conformational selection of protein kinase A revealed by flexible-ligand flexible-protein docking

Protein kinases have high structural plasticity: their structure can change significantly, depending on what ligands are bound to them. Rigid-protein docking methods are not capable of describing such effects. Here, we present a new flexible-ligand flexible-protein docking model in which the protein can adopt conformations between two extremes observed experimentally. The model utilized a molecular dynamics-based simulated annealing cycling protocol and a distance-dependent dielectric model to perform docking. By testing this model on docking four diverse ligands to protein kinase A, we found that the ligands were able to dock successfully to the protein with the proper conformations of the protein induced. By imposing relatively soft conformational restraints to the protein during docking, this model reduced computational costs yet permitted essential conformational changes that were essential for these inhibitors to dock properly to the protein. For example, without adequate movement of the glycine-rich loop, it was difficult for the ligands to move from the surface of the protein to the binding site. In addition, these simulations called for better ways to compare simulation results with experiment other than using the popular root-mean-square deviation between the structure of a ligand in a docking pose and that in experiment because the structure of the protein also changed. In this work, we also calculated the correlation coefficient between protein-ligand/protein-protein distances in the docking structure and those in the crystal structure to check how well a ligand docked into the binding site of the protein and whether the proper conformation of the protein was induced.     

ASEDock-docking based on alpha spheres and excluded volumes

ASEDock is a novel docking program based on a shape similarity assessment between a concave portion (i.e., concavity) on a protein and the ligand. We have introduced two novel concepts into ASEDock. One is an ASE model, which is defined by the combination of alpha spheres generated at a concavity in a protein and the excluded volumes around the concavity. The other is an ASE score, which evaluates the shape similarity between the ligand and the ASE model. The ASE score selects and refines the initial pose by maximizing the overlap between the alpha spheres and the ligand, and minimizing the overlap between the excluded volume and the ligand. Because the ASE score makes good use of the Gaussian-type function for evaluating and optimizing the overlap between the ligand and the site model, it can pose a ligand onto the docking site relatively faster and more effectively than using potential energy functions. The posing stage through the use of the ASE score is followed by full atomistic energy minimization. Because the posing algorithm of ASEDock is free from any bias except for shape, it is a very robust docking method. A validation study using 59 high-quality X-ray structures of the complexes between drug-like molecules and the target proteins has demonstrated that ASEDock can faithfully reproduce experimentally determined docking modes of various druglike molecules in their target proteins. Almost 80% of the structures were reconstructed within the estimated experimental error. The success rate of approximately 98% was attained based on the docking criterion of the root-mean-square deviation (RMSD) of non-hydrogen atoms (< or = 2.0 A). The markedly high success of ASEDock in redocking experiments clearly indicates that the most important factor governing the docking process is shape complementarity.     

Flexible matching of test ligands to a 3D pharmacophore using a molecular superposition force field: Comparison of predicted and experimental conformations of inhibitors of three enzymes

Summary A computer procedure TFIT, which uses a molecular superposition force field to flexibly match test compounds to a 3D pharmacophore, was evaluated to find out whether it could reliably predict the bioactive conformations of flexible ligands. The program superposition force field optimizes the overlap of those atoms of the test ligand and template that are of similar chemical type, by applying an attractive force between atoms of the test ligand and template which are close together and of similar type (hydrogen bonding, charge, hydrophobicity). A procedure involving Monte Carlo torsion perturbations, followed by torsional energy minimization, is used to find conformations of the test ligand which cominimize the internal energy of the ligand and the superposition energy of ligand and template. The procedure was tested by applying it to a series of flexible ligands for which the bioactive conformation was known experimentally. The 15 molecules tested were inhibitors of thermolysin, HIV-1 protease or endothiapepsin for which X-ray structures of the bioactive conformation were available. For each enzyme, one of the molecules served as a template and the others, after being conformationally randomized, were fitted. The fitted conformation was then compared to the known binding geometry. The matching procedure was successful in predicting the bioactive conformations of many of the structures tested. Significant deviation from experimental results was found only for parts of molecules where it was readily apparent that the template did not contain sufficient information to accurately determine the bioactive conformation.     

<Emphasis Type="Italic">In-silico</Emphasis> Screening using Flexible Ligand Binding Pockets: A Molecular Dynamics-based Approach

Summary In-silico screening of flexible ligands against flexible ligand binding pockets (LBP) is an emerging approach in structure-based drug discovery. Here, we describe a molecular dynamics (MD) based docking approach to investigate the influence on the high-throughput in-silico screening of small molecules against flexible ligand binding pockets. In our approach, an ensemble of 51 energetically favorable structures of the LBP of human estrogen receptor α (hERα) were collected from 3 ns MD simulations. In-silico screening of 3500 endocrine disrupting compounds against these flexible ligand binding pockets resulted in thousands of ER–ligand complexes of which 582 compounds were unique. Detailed analysis of MD generated structures showed that only 17 of the LBP residues significantly contribute to the overall binding pocket flexibility. Using the flexible LBP conformations generated, we have identified 32 compounds that bind better to the flexible ligand-binding pockets compared to the crystal structure. These compounds, though chemically divergent, are structurally similar to the natural hormone. Our MD-based approach in conjunction with grid–based distributed computing could be applied routinely for in-silico screening of large databases against any given target.