Exhaustive docking and solvated interaction energy scoring: lessons learned from the SAMPL4 challenge

We continued prospective assessments of the Wilma–solvated interaction energy (SIE) platform for pose prediction, binding affinity prediction, and virtual screening on the challenging SAMPL4 data sets including the HIV-integrase inhibitor and two host–guest systems. New features of the docking algorithm and scoring function are tested here prospectively for the first time. Wilma–SIE provides good correlations with actual binding affinities over a wide range of binding affinities that includes strong binders as in the case of SAMPL4 host–guest systems. Absolute binding affinities are also reproduced with appropriate training of the scoring function on available data sets or from comparative estimation of the change in target’s vibrational entropy. Even when binding modes are known, SIE predictions lack correlation with experimental affinities within dynamic ranges below 2 kcal/mol as in the case of HIV-integrase ligands, but they correctly signaled the narrowness of the dynamic range. Using a common protein structure for all ligands can reduce the noise, while incorporating a more sophisticated solvation treatment improves absolute predictions. The HIV-integrase virtual screening data set consists of promiscuous weak binders with relatively high flexibility and thus it falls outside of the applicability domain of the Wilma–SIE docking platform. Despite these difficulties, unbiased docking around three known binding sites of the enzyme resulted in over a third of ligands being docked within 2 Å from their actual poses and over half of the ligands docked in the correct site, leading to better-than-random virtual screening results.     

DOCK 4.0: Search strategies for automated molecular docking of flexible molecule databases

In this paper we describe the search strategies developed for docking flexible molecules to macomolecular sites that are incorporated into the widely distributed DOCK software, version 4.0. The search strategies include incremental construction and random conformation search and utilize the existing Coulombic and Lennard-Jones grid-based scoring function. The incremental construction strategy is tested with a panel of 15 crystallographic testcases, created from 12 unique complexes whose ligands vary in size and flexibility. For all testcases, at least one docked position is generated within 2 Å of the crystallographic position. For 7 of 15 testcases, the top scoring position is also within 2 Å of the crystallographic position. The algorithm is fast enough to successfully dock a few testcases within seconds and most within 100 s. The incremental construction and the random search strategy are evaluated as database docking techniques with a database of 51 molecules docked to two of the crystallographic testcases. Incremental construction outperforms random search and is fast enough to reliably rank the database of compounds within 15 s per molecule on an SGI R10000 cpu.     

Development and validation of a modular, extensible docking program: DOCK 5

We report on the development and validation of a new version of DOCK. The algorithm has been rewritten in a modular format, which allows for easy implementation of new scoring functions, sampling methods and analysis tools. We validated the sampling algorithm with a test set of 114 protein-ligand complexes. Using an optimized parameter set, we are able to reproduce the crystal ligand pose to within 2 A of the crystal structure for 79% of the test cases using our rigid ligand docking algorithm with an average run time of 1 min per complex and for 72% of the test cases using our flexible ligand docking algorithm with an average run time of 5 min per complex. Finally, we perform an analysis of the docking failures in the test set and determine that the sampling algorithm is generally sufficient for the binding pose prediction problem for up to 7 rotatable bonds; i.e. 99% of the rigid ligand docking cases and 95% of the flexible ligand docking cases are sampled successfully. We point out that success rates could be improved through more advanced modeling of the receptor prior to docking and through improvement of the force field parameters, particularly for structures containing metal-based cofactors.     

Blinded predictions of standard binding free energies: lessons learned from the SAMPL6 challenge

In the context of the SAMPL6 challenges, series of blinded predictions of standard binding free energies were made with the SOMD software for a dataset of 27 host–guest systems featuring two octa-acids hosts (OA and TEMOA) and a cucurbituril ring (CB8) host. Three different models were used, ModelA computes the free energy of binding based on a double annihilation technique; ModelB additionally takes into account long-range dispersion and standard state corrections; ModelC additionally introduces an empirical correction term derived from a regression analysis of SAMPL5 predictions previously made with SOMD. The performance of each model was evaluated with two different setups; buffer explicitly matches the ionic strength from the binding assays, whereas no-buffer merely neutralizes the host–guest net charge with counter-ions. ModelC/no-buffer shows the lowest mean-unsigned error for the overall dataset (MUE 1.29?<?1.39?<?1.50 kcal mol^?1, 95% CI), while explicit modelling of the buffer improves significantly results for the CB8 host only. Correlation with experimental data ranges from excellent for the host TEMOA (R² 0.91?<?0.94?<?0.96), to poor for CB8 (R² 0.04?<?0.12?<?0.23). Further investigations indicate a pronounced dependence of the binding free energies on the modelled ionic strength, and variable reproducibility of the binding free energies between different simulation packages.     

A Gibbs free energy correlation for automated docking of carbohydrates

Thermodynamic information can be inferred from static atomic configurations. To model the thermodynamics of carbohydrate binding to proteins accurately, a large binding data set has been assembled from the literature. The data set contains information from 262 unique protein-carbohydrate crystal structures for which experimental binding information is known. Hydrogen atoms were added to the structures and training conformations were generated with the automated docking program AutoDock 3.06, resulting in a training set of 225,920 all-atom conformations. In all, 288 formulations of the AutoDock 3.0 free energy model were trained against the data set, testing each of four alternate methods of computing the van der Waals, solvation, and hydrogen-bonding energetic components. The van der Waals parameters from AutoDock 1 produced the lowest errors, and an entropic model derived from statistical mechanics produced the only models with five physically and statistically significant coefficients. Eight models predict the Gibbs free energy of binding with an error of less than 40% of the error of any similar models previously published.     

SAMPL4, a blind challenge for computational solvation free energies: the compounds considered

For the fifth time I have provided a set of solvation energies (1 M gas to 1 M aqueous) for a SAMPL challenge. In this set there are 23 blind compounds and 30 supplementary compounds of related structure to one of the blind sets, but for which the solvation energy is readily available. The best current values of each compound are presented along with complete documentation of the experimental origins of the solvation energies. The calculations needed to go from reported data to solvation energies are presented, with particular attention to aspects which are new to this set. For some compounds the vapor pressures (VP) were reported for the liquid compound, which is solid at room temperature. To correct from VP_{subcooled liquid} to VP_sublimation requires ΔS_fusion, which is only known for mannitol. Estimated values were used for the others, all but one of which were benzene derivatives and expected to have very similar values. The final compound for which ΔS_fusion was estimated was menthol, which melts at 42 °C so that modest errors in ΔS_fusion will have little effect. It was also necessary to look into the effects of including estimated values of ΔCp on this correction. The approximate sizes of the effects of inclusion of ΔCp in the correction from VP_{subcooled liquid} to VP_sublimation were estimated and it was noted that inclusion of ΔCp invariably makes ΔG_S more positive. To extend the set of compounds for which the solvation energy could be calculated we explored the use of boiling point (b.p.) data from Reaxys/Beilstein as a substitute for studies of the VP as a function of temperature. B.p. data are not always reliable so it was necessary to develop a criterion for rejecting outliers. For two compounds (chlorinated guaiacols) it became clear that inclusion represented overreach; for each there were only two independent pressure, temperature points, which is too little for a trustworthy extrapolation. For a number of compounds the extrapolation from lowest temperature at which the VP was reported to 25 °C was long (sometimes over 100°) so that it was necessary to consider whether ΔCp might have significant effects. The problem is that there are no experimental values and possible intramolecular hydrogen bonds make estimation uncertain in some cases. The approximate sizes of the effects of ΔCp were estimated, and it was noted that inclusion of ΔCp in the extrapolation of VP down to room temperature invariably makes ΔGs more negative.     

Predicting hydration free energies with chemical accuracy: the SAMPL4 challenge

An implicit solvent model described by a non-simple dielectric medium is used for the prediction of hydration free energies on the dataset of 47 molecules in the SAMPL4 challenge. The solute is represented by a minimal parameter set model based on a new all atom force-field, named the liquid simulation force-field. The importance of a first solvation shell correction to the hydration free energy prediction is discussed and two different approaches are introduced to address it: either with an empirical correction to a few functional groups (alcohol, ether, ester, amines and aromatic nitrogen), or an ab initio correction based on the formation of a solute/explicit water complex. Both approaches give equally good predictions with an average unsigned error <1 kcal/mol. Chemical accuracy is obtained.     

Specific empirical free energy function for automated docking of carbohydrates to proteins

We present an automated docking protocol specifically optimized to predict the structure and affinity of a protein-carbohydrate complex. A scoring function was developed based on a training set of 30 protein-carbohydrate complexes of known structure and affinity. Combinations of several models for hydrogen bonding, torsional entropy loss, and solvation were tested for their ability to fit the training set data, and the best model was used with AutoDock. The electrostatic empirical coefficient is larger than in a previously obtained model using a training set comprised of various types of protein-ligand complexes, indicating that electrostatic interactions play a more important role in determining the affinity between a carbohydrate and a protein. The differences in the relative weighting of the empirical coefficients in the model yields predicted free energies for the training set with a standard error of 1.403 kcal/mol. The new scoring function was tested on 17 Aspergillus niger glucoamylase inhibitors for which binding energies had been determined experimentally. Free energies of complex formation were predicted with a residual standard error of 1.101 kcal/mol. The new scoring function therefore provides a robust method for predicting free energies of formation and optimal conformations of carbohydrate-protein complexes.     

Critical evaluation of search algorithms for automated molecular docking and database screening

The DOCK program explores possible orientations of a molecule within a macromolecular active site by superimposing atoms onto precomputed site points. Here we compare a number of different search methods, including an exhaustive matching algorithm based on a single docking graph. We evaluate the performance of each method by screening a small database of molecules to a variety of macromolecular targets. By varying the amount of sampling, we can monitor the time convergence of scores and rankings. We not only show that the site point–directed search is tenfold faster than a random search, but that the single graph matching algorithm boosts the speed of database screening up to 60-fold. The new algorithm, in fact, outperforms the bipartite graph matching algorithm currently used in DOCK. The results indicate that a critical issue for rapid database screening is the extent to which a search method biases run time toward the highest-ranking molecules. The single docking graph matching algorithm will be incorporated into DOCK version 4.0. © 1997 John Wiley & Sons, Inc. J Comput Chem 18: 1175–1189     

Manual and automated enrichment procedures for biological samples using lipophilic gels

Aspects of the use of lipophilic gels in manual sample preparation procedures are reviewed. Neutral gels with a controlled hydrophobicity are used for sorbent extraction of non-polar and medium polarity compounds from biological fluids. Acidic amphiphilic compounds can be extracted as ion-pairs with decyltrimethylammonium ions. Solvent or detergent extracts of tissues or faeces can be mixed with hydrophobic gels for transfer of analytes from a solvent to a gel phase, permitting subsequent sample preparation in gel bed systems. Hydrophobic gels, alkyl-bonded silica and polystyrene matrices can be used in series for extraction of compounds with a wide range of polarities. Group fractionations are performed on neutral and ion-exchanging lipophilic gels to yield fractions of neutral, basic and acidic metabolites within selected polarity ranges. Selective isolation of phenolic acids on a strong anion exchanger, of ethynylic steroids on a strong cation exchanger in silver form and of oximes of ketonic steroids on a strong cation exchanger in hydrogen form is possible. A computerized system for automatic sample preparation is also described. It consists of an extraction bed, a cation-exchange column and an anion-exchange column. The pumps and switching valves are arranged so that the columns can operate in series or parallel for isolation of neutral, basic and acidic metabolites of amphiphilic compounds and for regeneration of the column beds. Fractions can be collected, or the effluent from the column beds can be diluted with water to permit sorption on a solid phase. The applicability of the automated method to the analysis of bile acids and metabolites of mono(2-ethylhexyl) phthalate is demonstrated.     

Packing optimization for automated generation of complex system's initial configurations for molecular dynamics and docking

Molecular Dynamics is a powerful methodology for the comprehension at molecular level of many chemical and biochemical systems. The theories and techniques developed for structural and thermodynamic analyses are well established, and many software packages are available. However, designing starting configurations for dynamics can be cumbersome. Easily generated regular lattices can be used when simple liquids or mixtures are studied. However, for complex mixtures, polymer solutions or solid adsorbed liquids (for example) this approach is inefficient, and it turns out to be very hard to obtain an adequate coordinate file. In this article, the problem of obtaining an adequate initial configuration is treated as a "packing" problem and solved by an optimization procedure. The initial configuration is chosen in such a way that the minimum distance between atoms of different molecules is greater than a fixed tolerance. The optimization uses a well-known algorithm for box-constrained minimization. Applications are given for biomolecule solvation, many-component mixtures, and interfaces. This approach can reduce the work of designing starting configurations from days or weeks to few minutes or hours, in an automated fashion. Packing optimization is also shown to be a powerful methodology for space search in docking of small ligands to proteins. This is demonstrated by docking of the thyroid hormone to its nuclear receptor.     

Extended solvent-contact model approach to SAMPL4 blind prediction challenge for hydration free energies

Extended solvent-contact model was applied to the blind prediction of the hydration free energies of 47 organic molecules included in the SAMPL4 data set. To obtain a suitable prediction tool, we constructed a hydration free energy function involving three kinds of atomic parameters. With respect to total 34 atom types introduced to describe all SAMPL4 molecules, 102 atomic parameters were defined and optimized with a standard genetic algorithm in such a way to minimize the difference between the experimental hydration free energies and those calculated with the hydration free energy function. In this parameterization, we used a training set comprising 77 organic molecules with varying sizes and shapes. The estimated hydration free energies for the SAMPL4 molecules compared reasonably well with the experimental results with the associated squared correlation coefficient and root mean square deviation of 0.89 and 1.46 kcal/mol, respectively. Based on the comparative analysis of experimental and computational hydration free energies of the SAMPL4 molecules, the methods for further improvement of the present hydration model are suggested.     

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.     

Distributed automated docking of flexible ligands to proteins: Parallel applications of AutoDock 2.4

Summary AutoDock 2.4 predicts the bound conformations of a small, flexible ligand to a nonflexible macromolecular target of known structure. The technique combines simulated annealing for conformation searching with a rapid grid-based method of energy evaluation based on the AMBER force field. AutoDock has been optimized in performance without sacrificing accuracy; it incorporates many enhancements and additions, including an intuitive interface. We have developed a set of tools for launching and analyzing many independent docking jobs in parallel on a heterogeneous network of UNIX-based workstations. This paper describes the current release, and the results of a suite of diverse test systems. We also present the results of a systematic investigation into the effects of varying simulated-annealing parameters on molecular docking. We show that even for ligands with a large number of degrees of freedom, root-mean-square deviations of less than 1 Å from the crystallographic conformation are obtained for the lowest-energy dockings, although fewer dockings find the crystallographic conformation when there are more degrees of freedom.The AutoDock 2.4 suite is written in ANSI C, and is supplied with Makefiles for the following platforms: Convex, DEC Alpha OSF/1, Hewlett-Packard Precision Architecture, Silicon Graphics, and Sun. The AutoDock suite of programs is freely available to the noncommercial scientific community and to educational establishments. Further information, including additional figures and MPEG animations showing all docked conformations for each test system, can be found at the following URL: http://www.scripps.edu/pub/olson-web/doc/autodock.     

Modern automated microextraction procedures for bioanalytical,environmental, and food analyses

Sample preparation frequently is considered the most critical stage of the analytical workflow. It affects the analytical throughput and costs; moreover, it is the primary source of error and possible sample contamination. To increase efficiency, productivity, and reliability, while minimizing costs and environmental impacts, miniaturization and automation of sample preparation are necessary. Nowadays, several types of liquid-phase and solid-phase microextractions are available, as well as different automatization strategies. Thus, this review summarizes recent developments in automated microextractions coupled with liquid chromatography, from 2016 to 2022. Therefore, outstanding technologies and their main outcomes, as well as miniaturization and automation of sample preparation, are critically analyzed. Focus is given to main microextraction automation strategies, such as flow techniques, robotic systems, and column-switching approaches, reviewing their applications to the determination of small organic molecules in biological, environmental, and food/beverage samples.     

Combining in silico and in cerebro approaches for virtual screening and pose prediction in SAMPL4

The SAMPL challenges provide an ideal opportunity for unbiased evaluation and comparison of different approaches used in computational drug design. During the fourth round of this SAMPL challenge, we participated in the virtual screening and binding pose prediction on inhibitors targeting the HIV-1 integrase enzyme. For virtual screening, we used well known and widely used in silico methods combined with personal in cerebro insights and experience. Regular docking only performed slightly better than random selection, but the performance was significantly improved upon incorporation of additional filters based on pharmacophore queries and electrostatic similarities. The best performance was achieved when logical selection was added. For the pose prediction, we utilized a similar consensus approach that amalgamated the results of the Glide-XP docking with structural knowledge and rescoring. The pose prediction results revealed that docking displayed reasonable performance in predicting the binding poses. However, prediction performance can be improved utilizing scientific experience and rescoring approaches. In both the virtual screening and pose prediction challenges, the top performance was achieved by our approaches. Here we describe the methods and strategies used in our approaches and discuss the rationale of their performances.     

AutoDock4 and AutoDockTools4: Automated docking with selective receptor flexibility

We describe the testing and release of AutoDock4 and the accompanying graphical user interface AutoDockTools. AutoDock4 incorporates limited flexibility in the receptor. Several tests are reported here, including a redocking experiment with 188 diverse ligand‐protein complexes and a cross‐docking experiment using flexible sidechains in 87 HIV protease complexes. We also report its utility in analysis of covalently bound ligands, using both a grid‐based docking method and a modification of the flexible sidechain technique. © 2009 Wiley Periodicals, Inc. J Comput Chem, 2009     

Virtual screening of integrase inhibitors by large scale binding free energy calculations: the SAMPL4 challenge

As part of the SAMPL4 blind challenge, filtered AutoDock Vina ligand docking predictions and large scale binding energy distribution analysis method binding free energy calculations have been applied to the virtual screening of a focused library of candidate binders to the LEDGF site of the HIV integrase protein. The computational protocol leveraged docking and high level atomistic models to improve enrichment. The enrichment factor of our blind predictions ranked best among all of the computational submissions, and second best overall. This work represents to our knowledge the first example of the application of an all-atom physics-based binding free energy model to large scale virtual screening. A total of 285 parallel Hamiltonian replica exchange molecular dynamics absolute protein-ligand binding free energy simulations were conducted starting from docked poses. The setup of the simulations was fully automated, calculations were distributed on multiple computing resources and were completed in a 6-weeks period. The accuracy of the docked poses and the inclusion of intramolecular strain and entropic losses in the binding free energy estimates were the major factors behind the success of the method. Lack of sufficient time and computing resources to investigate additional protonation states of the ligands was a major cause of mispredictions. The experiment demonstrated the applicability of binding free energy modeling to improve hit rates in challenging virtual screening of focused ligand libraries during lead optimization.