An efficient molecular docking using conformational space annealing

Molecular docking falls into the general category of global optimization problems because its main purpose is to find the most stable complex consisting of a receptor and its ligand. Conformational space annealing (CSA), a powerful global optimization method, is incorporated with the Tinker molecular modeling package to perform molecular docking simulations of six receptor-ligand complexes (3PTB, 1ULB, 2CPP, 1STP, 3CPA, and 1PPH) from the Protein Data Bank. In parallel, Monte Carlo with the minimization (MCM) method is also incorporated into the Tinker package for comparison. The energy function, consisting of electrostatic interactions, van der Waals interactions, and torsional energy terms, is calculated using the AMBER94 all-atom empirical force field. Rigid docking simulations for all six complexes and flexible docking simulations for three complexes (1STP, 3CPA, and 1PPH) are carried out using the CSA and the MCM methods. The simulation results show that the docking procedures using the CSA method generally find the most stable complexes as well as the native-like complexes more efficiently and accurately than those using the MCM, demonstrating that CSA is a promising search method for molecular docking problems.     

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.     

Blind docking method combining search of low‐resolution binding sites with ligand pose refinement by molecular dynamics‐based global optimization

This study describes the development of a new blind hierarchical docking method, bhDock, its implementation, and accuracy assessment. The bhDock method uses two‐step algorithm. First, a comprehensive set of low‐resolution binding sites is determined by analyzing entire protein surface and ranked by a simple score function. Second, ligand position is determined via a molecular dynamics‐based method of global optimization starting from a small set of high ranked low‐resolution binding sites. The refinement of the ligand binding pose starts from uniformly distributed multiple initial ligand orientations and uses simulated annealing molecular dynamics coupled with guided force‐field deformation of protein–ligand interactions to find the global minimum. Assessment of the bhDock method on the set of 37 protein–ligand complexes has shown the success rate of predictions of 78%, which is better than the rate reported for the most cited docking methods, such as AutoDock, DOCK, GOLD, and FlexX, on the same set of complexes. © 2009 Wiley Periodicals, Inc. J Comput Chem 2010     

GalaxyDock2: Protein–ligand docking using beta‐complex and global optimization

In this article, an enhanced version of GalaxyDock protein–ligand docking program is introduced. GalaxyDock performs conformational space annealing (CSA) global optimization to find the optimal binding pose of a ligand both in the rigid‐receptor mode and the flexible‐receptor mode. Binding pose prediction has been improved compared to the earlier version by the efficient generation of high‐quality initial conformations for CSA using a predocking method based on a beta‐complex derived from the Voronoi diagram of receptor atoms. Binding affinity prediction has also been enhanced by using the optimal combination of energy components, while taking into consideration the energy of the unbound ligand state. The new version has been tested in terms of binding mode prediction, binding affinity prediction, and virtual screening on several benchmark sets, showing improved performance over the previous version and AutoDock, on which the GalaxyDock energy function is based. GalaxyDock2 also performs better than or comparable to other state‐of‐the‐art docking programs. GalaxyDock2 is freely available at http://galaxy.seoklab.org/softwares/galaxydock.html . © 2013 Wiley Periodicals, Inc.     

A new method for the gradient‐based optimization of molecular complexes

We present a novel method for the local optimization of molecular complexes. This new approach is especially suited for usage in molecular docking. In molecular modeling, molecules are often described employing a compact representation to reduce the number of degrees of freedom. This compact representation is realized by fixing bond lengths and angles while permitting changes in translation, orientation, and selected dihedral angles. Gradient‐based energy minimization of molecular complexes using this representation suffers from well‐known singularities arising during the optimization process. We suggest an approach new in the field of structure optimization that allows to employ gradient‐based optimization algorithms for such a compact representation. We propose to use exponential mapping to define the molecular orientation which facilitates calculating the orientational gradient. To avoid singularities of this parametrization, the local minimization algorithm is modified to change efficiently the orientational parameters while preserving the molecular orientation, i.e. we perform well‐defined jumps on the objective function. Our approach is applicable to continuous, but not necessarily differentiable objective functions. We evaluated our new method by optimizing several ligands with an increasing number of internal degrees of freedom in the presence of large receptors. In comparison to the method of Solis and Wets in the challenging case of a non‐differentiable scoring function, our proposed method leads to substantially improved results in all test cases, i.e. we obtain better scores in fewer steps for all complexes. © 2008 Wiley Periodicals, Inc. J Comput Chem, 2009     

Fast and accurate grid representations for atom‐based docking with partner flexibility

Macromolecular docking methods can broadly be divided into geometric and atom‐based methods. Geometric methods use fast algorithms that operate on simplified, grid‐like molecular representations, while atom‐based methods are more realistic and flexible, but far less efficient. Here, a hybrid approach of grid‐based and atom‐based docking is presented, combining precalculated grid potentials with neighbor lists for fast and accurate calculation of atom‐based intermolecular energies and forces. The grid representation is compatible with simultaneous multibody docking and can tolerate considerable protein flexibility. When implemented in our docking method ATTRACT, grid‐based docking was found to be ∼35x faster. With the OPLSX forcefield instead of the ATTRACT coarse‐grained forcefield, the average speed improvement was >100x. Grid‐based representations may allow atom‐based docking methods to explore large conformational spaces with many degrees of freedom, such as multiple macromolecules including flexibility. This increases the domain of biological problems to which docking methods can be applied. © 2017 Wiley Periodicals, Inc.     

ELECT++: Faster conformational search method for docking flexible molecules using molecular similarity

We have developed a program, ELECT++ (Effective LEssening of Conformations by Template molecules in C++), to speed up the conformational search for small flexible molecules using the similar property principle. We apply this principle to molecular shape and, importantly, to molecular flexibility. After molecules in a database are clustered according to flexibility and shape (FCLUST++), additional reagents are generated to screen the conformational space of molecules in each cluster (TEMPLATE++). We call these representative reagents of each cluster template reagents. Template reagents and clustered reagents produce, after reaction, template molecules and clustered molecules, respectively (tREACT++). The conformations of a template molecule are searched in the context of a macromolecular target. Acceptable conformational choices are then applied to all molecules in its cluster, thus effectively biasing conformational space to speed up conformational searches (tSEARCH++). In our incremental search method, it is necessary to calculate the root-mean-square deviations (RMSD) matrix of distances between different conformations of the same molecule to reduce the number of conformations. Instead of calculating the RMSD matrix for all molecules in a cluster, the RMSD matrix of a template molecule is chosen as a reference and applied to all the molecules in its cluster. We demonstrate that FCLUST++ clusters the primary amine reagents from the Available Chemicals Directory (ACD) successfully. The program tSEARCH++ was applied to dihydrofolate reductase with virtual molecules generated by tREACT++ using clustered primary amine reagents. The conformational search by the program tSEARCH++ was about 4.8 times faster than by SEARCH++, with an acceptable range of errors. © 1998 John Wiley & Sons, Inc. J Comput Chem 19: 1834–1852, 1998     

The molecular weight distribution for chain polymerization plus side reaction

The kinetics of chain polymerization is investigated for the case of a complicating side reaction. In addition to the polymerization reaction, A_i + M → A_i+1, there is a reversible side reaction, A_i + Q ↔ B_i. Initiation is assumed to be instantaneous. The monomer concentration M, and the concentration of the reacting species Q, are assumed to be constant. The reaction kinetics are solved exactly, yielding the distribution of living and dormant polymer, as well as the molecular weight distribution, as explicit functions of the reaction rate constants and the time. © 1997 John Wiley & Sons, Inc. J Polym Sci B: Polym Phys 35: 1711–1725, 1997     

Ligand-protein docking using a quantum stochastic tunneling optimization method

A novel hybrid optimization method called quantum stochastic tunneling has been recently introduced. Here, we report its implementation within a new docking program called EasyDock and a validation with the CCDC/Astex data set of ligand-protein complexes using the PLP score to represent the ligand-protein potential energy surface and ScreenScore to score the ligand-protein binding energies. When taking the top energy-ranked ligand binding mode pose, we were able to predict the correct crystallographic ligand binding mode in up to 75% of the cases. By using this novel optimization method run times for typical docking simulations are significantly shortened.     

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     

Rapid and efficient screening of adsorbent for oligopeptide using molecular docking and isothermal titration calorimetry

A rapid and efficient method based on molecular docking and isothermal titration calorimetry (ITC) was developed to identify effective adsorbents for the target peptide Ser‐Glu‐Ala‐Asp‐His (SEADH). Preliminary screening of five candidate adsorbents using molecular docking revealed that three suitable structures (A1, A2, and A3) either with or without coordination of Zn²⁺ should be effective. The three selected structures were then prepared and further screened by evaluating their affinities for the peptide SEADH using ITC. The screening results revealed that the adsorbent A2 coordinated with Zn²⁺ was the best adsorbent, and subsequent static adsorption experiments confirmed the reliability of the screening method. Further ITC analysis, combined with molecular docking, was performed to provide the possible adsorption mechanism.     

Oxidation‐assisted structural elucidation of compounds containing a tertiary amine side chain using liquid chromatography mass spectrometry

A novel analytical technique for the structural elucidation of compounds bearing a tertiary amine side chain via “in vial” instantaneous oxidation and liquid chromatography mass spectrometry (LC‐MS) was developed. A series of lidocaine homologs and benzimidazole derivatives with a major/single amine representative base peak in both their EI‐MS and ESI‐MS/MS spectra were subjected to oxidation by a 0.1% solution of hydrogen peroxide (including several ¹⁶O/¹⁸O exchange experiments), followed by LC‐ESI‐MS/MS analysis. The N‐oxide counterparts promoted extensive fragmentation with complete coverage of all parts of the molecule, enabling detailed structural elucidation and unambiguous identification of the unoxidized analytes at low nanogram per milliliter levels.     

A hybrid method of molecular dynamics and harmonic dynamics for docking of flexible ligand to flexible receptor

We have developed a new docking method to consider receptor flexibility, a hybrid method of molecular dynamics and harmonic dynamics. The global motions of the whole receptor were approximately introduced into those of the receptor in the docking simulation as harmonic dynamics. On the other hand, the local flexibility of the side chains was also considered by conventional molecular dynamics. We confirmed that this new method can reproduce the fluctuations of the whole receptor by making a comparison of the directions and amplitudes of the global fluctuations. Then this method was applied to the docking of HIV-1 protease and its ligand. As a result, we observed a docking process where the ligand enters into the binding pocket well, which implies that this method is effective enough to reproduce a molecular complex formation.     

Cluster‐based molecular docking study for in silico identification of novel 6‐fluoroquinolones as potential inhibitors against Mycobacterium tuberculosis

A classical protein sequence alignment and homology modeling strategy were used for building three Mycobacterium tuberculosis‐DNA gyrase protein models using the available topoII‐DNA‐6FQ crystal structure complexes originating from different organisms. The recently determined M. tuberculosis‐DNA gyrase apoprotein structures and topoII‐DNA‐6FQ complexes were used for defining the 6‐fluoroquinolones (6‐FQs) binding pockets. The quality of the generated models was initially validated by docking of the cocrystallized ligands into their binding site, and subsequently by quantitative evaluation of their discriminatory performances (identification of active/inactive 6‐FQs) for a set of 145 6‐FQs with known biological activity values. The M. tuberculosis‐DNA gyrase model with the highest estimated discriminatory power was selected and used afterwards in an additional molecular docking experiment on a mixed combinatorial set of 427 drug‐like 6‐FQ analogs for which the biological activity values were predicted using a prebuilt counter‐propagation artificial neural network model. A novel three‐level Boolean‐based [T/F (true/false)] clustering algorithm was used to assess the generated binding poses: Level 1 (geometry properties assessment), Level 2 (score‐based clustering and selection of the (T)‐signed highly scored Level 1 poses), and Level 3 (activity‐based clustering and selection of the most “active” (T)‐signed Level 2 hits). The frequency analysis of occurrence of the fragments attached at R₁ and R₇ position of the (T)‐signed 6‐FQs selected in Level 3 revealed several novel attractive fragments and confirmed some previous findings. We believe that this methodology could be successfully used in establishing novel possible structure‐activity relationship recommendations in the 6‐FQs optimization, which could be of great importance in the current antimycobacterial hit‐to‐lead processes. © 2012 Wiley Periodicals, Inc.     

A novel scoring function for molecular docking

We present a novel scoring function for docking of small molecules to protein binding sites. The scoring function is based on a combination of two main approaches used in the field, the empirical and knowledge-based approaches. To calibrate the scoring function we used an iterative procedure in which a ligand's position and its score were determined self-consistently at each iteration. The scoring function demonstrated superiority in prediction of ligand positions in docking tests against the commonly used Dock, FlexX and Gold docking programs. It also demonstrated good accuracy of binding affinity prediction for the docked ligands.     

Evaluation of the performance of four molecular docking programs on a diverse set of protein‐ligand complexes

Many molecular docking programs are available nowadays, and thus it is of great practical value to evaluate and compare their performance. We have conducted an extensive evaluation of four popular commercial molecular docking programs, including Glide, GOLD, LigandFit, and Surflex. Our test set consists of 195 protein‐ligand complexes with high‐resolution crystal structures (resolution ≤2.5 Å) and reliable binding data [dissociation constant (K_d) or inhibition constant (K_i)], which are selected from the PDBbind database with an emphasis on diversity. The top‐ranked solutions produced by these programs are compared to the native ligand binding poses observed in crystal structures. Glide and GOLD demonstrate better accuracy than the other two on the entire test set. Their results are also less sensitive to the starting structures for docking. Comparison of the results produced by these programs at three different computation levels reveal that their accuracy are not always proportional to CPU cost as one may expect. The binding scores of the top‐ranked solutions produced by these programs are in low to moderate correlations with experimentally measured binding data. Further analyses on the outcomes of these programs on three suites of subsets of protein‐ligand complexes indicate that these programs are less capable to handle really flexible ligands and relatively flat binding sites, and they have different preferences to hydrophilic/hydrophobic binding sites. Our evaluation can help other researchers to make reasonable choices among available molecular docking programs. It is also valuable for program developers to improve their methods further. © 2010 Wiley Periodicals, Inc. J Comput Chem, 2010     

Theoretical investigation of the molecular structure,vibrational spectra,and molecular docking of tramadol using density functional theory

The optimized molecular structures, harmonic vibrational wavenumbers, and the corresponding vibrational assignments of (1S,2S)-tramadol and (1R,2R)-tramadol are computationally examined using the B3LYP density functional theory method together with the standard 6–311++G(d,p) and def2-TVZP basis sets. The optimized structures show that phenolic rings of both 1R,2R and 1S,2S tramadol adopt planar geometry, which are slightly distorted due to the substitution at the meta-position; and the six-membered cyclohexane adopts a slightly distorted chair conformation. The 1S,2S enantiomer is energetically more favorable than 1R,2R with the energy differences of 1.32 and 1.03 kcal/mol obtained at B3LYP/6–311++G(d,p) and B3LYP/Def2-TVZP levels, respectively. The analysis of the binding pocket in the silico molecular docking with the m-opioid receptor shows that it originated two clusters with the 1S,2S enantiomer and one cluster with the 1R,2R enantiomer of tramadol. The results point to a more stable complex of the m-opioid receptor with the 1R,2R enantiomer of tramadol.     

DNA interaction,antimicrobial and molecular docking studies of biologically interesting Schiff base complexes incorporating 4‐formyl‐N,N‐dimethylaniline and propylenediamine

Four new transition metal complexes incorporating a Schiff base ligand derived from propylenediamine and 4‐formyl‐N ,N ‐dimethylaniline have been synthesized using transition metal salts. The characterization of the newly formed complexes was done from physicochemical parameters and using various techniques like ¹H NMR, ¹³C NMR, IR, UV, electron paramagnetic resonance and mass spectroscopies, powder X‐ray diffraction and magnetic susceptibility. All the complexes were found to be monomeric in nature with square planar geometry. X‐ray powder diffraction illustrates that the complexes have a crystalline nature. The interaction of metal complexes with calf thymus DNA was investigated using UV–visible absorption, viscosity measurements, cyclic voltammetry, emission spectroscopy and docking analysis. The results indicate that the Cu(II), Co(II), Ni(II) and Zn(II) complexes interact with DNA by intercalative binding mode with optimum intrinsic binding constants of 4.3 × 10⁴, 3.9 × 10⁴, 4.7 × 10⁴ and 3.7 × 10⁴ M⁻¹, respectively. These DNA binding results were rationalized using molecular docking in which the docked structures indicate that the metal complexes fit well into the A‐T rich region of target DNA through intercalation. The metal complexes exhibit an effective cleavage with pUC19 DNA by an oxidative cleavage mechanism. The synthesized ligand and the complexes were tested for their in vitro antimicrobial activity. The complexes show enhanced antifungal and antibacterial activities compared to the free ligand.     

Polyesters with semifluorinated side chains: A proposal for the solid‐state structure

A proposal for the solid‐state structure of poly(p‐phenylene isophthalate) with oxydecylperfluorodecyl side chains is presented, which was calculated by the new Rietveld refinement program BGMN®. A triclinic unit cell with a = 0.575 nm, b = 4.06 nm, c = 2.1 nm, α = 91.2°, β = 85.7°, and γ = 66.1° was obtained using space group symmetry P1¯ (No. 2) with Z = 2. All fractional atom coordinates, reflections and structure factors F were determined. The results show that highly occupied netplanes lie in plane and perpendicularly to the side chains. It may be supposed that these net planes will form the low‐energy surface of the polymer. © 2000 John Wiley & Sons, Inc. J Polym Sci B: Polym Phys 38: 1617–1625, 2000