Improving virtual screening results with MM/GBSA and MM/PBSA rescoring

Journal of Computer-Aided Molecular Design - Virtual screening (VS) based on molecular docking is one of the most useful methods in computer-aided drug design. By allowing to identify...     

Efficient treatment of out-of-plane bend and improper torsion interactions in MM2, MM3, and MM4 molecular mechanics calculations

Simple and very efficient formulas are presented for four-body out-of-plane bend (used in MM2 and MM3 force fields) and improper torsion (used in the MM4 force field) internal coordinates and their first and second derivatives. The use of a small set of bend and stretch intermediates allows for order of magnitude decreases in calculation time for potential energies and their first and second derivatives, which are required in molecular mechanics calculations. The formulas are eminently suitable for use in molecular simulations of systems with complicated bond networks. © 1997 John Wiley & Sons, Inc. J Comput Chem 18 : 1804–1811, 1997     

MM and QM/MM Modeling of Threonyl-tRNA Synthetase: Model Testing and Simulations

Aminoacyl-tRNA synthetases are centrally important enzymes in protein synthesis. We have investigated threonyl-tRNA synthetase from E. coli, complexed with reactants, using molecular mechanics and combined quantum mechanical/molecular mechanical (QM/MM) techniques. These modeling methods have the potential to provide molecular level understanding of enzyme catalytic processes. Modeling of this enzyme presents a number of challenges. The procedure of system preparation and testing is described in detail. For example, the number of metal ions at the active site, and their positions, were investigated. Molecular dynamics simulations suggest that the system is most stable when it contains only one magnesium ion, and the zinc ion is removed. Two different QM/MM methods were tested in models based on the findings of MM molecular dynamics simulations. AM1/CHARMM calculations resulted in unrealistic structures for the phosphates in this system. This is apparently due to an error of AM1. PM3/CHARMM calculations proved to be more suitable for this enzyme system. These results will provide a useful basis for future modeling investigations of the enzyme mechanism and dynamics.     

Molecular mechanics calculations (MM2 and MM3) on enamines and aniline derivatives

The MM2 and MM3 force fields have been extended to cover this class of compounds. Structures, vibrational spectra, and other data for 13 compounds were examined and can be reproduced satisfactorily by MM3. Except for the spectra, the other data can be reproduced somewhat less well by MM2. © 1994 by John Wiley & Sons, Inc.     

QM/MM vibrational mode tracking

Vibrational spectroscopy is a powerful tool to investigate the structure and dynamics of biomolecules. When small subsystems of large molecules such as active centers of enzymes are studied, quantum chemical calculations based on quantum mechanics/molecular mechanics (QM/MM) coupling schemes are a valuable means to interpret the spectra. The goal of this work is a methodological pilot study on how to selectively and thus efficiently extract certain vibrational information for extended molecular systems described by QM/MM methods. This is achieved by an extension of the mode tracking algorithm and a comparison with the partial Hessian diagonalization approach. After validating the methodology for the CO stretching vibration of 2-butanone and a delocalized CO stretch in acetylacetone, the stretching and bending modes of the CO ligand in CO myoglobin are tracked. Such systems represent an ideal application for mode tracking, because only a few strongly localized vibrations are sought for, while the large remainder of the molecule is of interest only as far as it affects these local vibrations. This influence is treated exactly by mode tracking.     

Density embedded VB/MM: a hybrid ab initio VB/MM with electrostatic embedding

A hybrid QM/MM method that combines ab initio valence-bond (VB) with molecular mechanics (MM) is presented. The method utilizes the ab initio VB approach to describe the reactive fragments and MM to describe the environment thus allows VB calculations of reactions in large biological systems. The method, termed density embedded VB/MM (DE-VB/MM), is an extension of the recently developed VB/MM method. It involves calculation of the electrostatic interaction between the reactive fragments and their environment using the electrostatic embedding scheme. Namely, the electrostatic interactions are represented as one-electron integrals in the ab initio VB Hamiltonian, hence taking into account the wave function polarization of the reactive fragments due to the environment. Moreover, the assumptions that were utilized in an earlier version of the method, VB/MM, to formulate the electrostatic interactions effect on the off-diagonal matrix elements are no longer required in the DE-VB/MM methodology. Using DE-VB/MM, one can calculate, in addition to the adiabatic ground state reaction profile, the energy of the diabatic VB configurations as well as the VB state correlation diagram for the reaction. The abilities of the method are exemplified on the identity SN2 reaction of a chloride anion with methyl chloride in aqueous solution. Both the VB configurations diagram and the state correlation diagram are presented. The results are shown to be in very good agreement with both experimental and other computational data, suggesting that DE-VB/MM is a proper method for application to different reactivity problems in biological systems.     

How Photoisomerization Drives Peptide Folding and Unfolding: Insights from QM/MM and MM Dynamics Simulations

Photoswitchable azobenzene cross‐linkers can control the folding and unfolding of peptides by photoisomerization and can thus regulate peptide affinities and enzyme activities. Using quantum mechanics/molecular mechanics (QM/MM) methods and classical MM force fields, we report the first molecular dynamics simulations of the photoinduced folding and unfolding processes in the azobenzene cross‐linked FK‐11 peptide. We find that the interactions between the peptide and the azobenzene cross‐linker are crucial for controlling the evolution of the secondary structure of the peptide and responsible for accelerating the folding and unfolding events. They also modify the photoisomerization mechanism of the azobenzene cross‐linker compared with the situation in vacuo or in solution.     

Development of an internal searching algorithm for parameterization of the MM2/MM3 force fields

Application of Allinger's MM2/MM3 force fields to molecules of real interest is frequently hindered by the lack of parameters for various heterocyclic systems and for poly-functionalized molecules. A common approach to this problem is to manually choose missing parameters “by analogy” with those that are part of the force field's internal parameter set. Naturally, this is generally attempted only by those possessing extensive experience with force fields. In order to use the MM2/MM3 force fields to study herbicides, an algorithm has been developed to automate this process for the non MM2 specialist. Using a set of “relative cost” criteria for atom type replacement, the algorithm searches the force field parameter set and selects the most appropriate parameters for a given molecule whose MM2 output file contains “missing parameter” errors. The program selects parameter error messages from a standard MM2 output file, finds analogous parameters, asks the user to verify their appropriateness and creates a standard MM2 parameter deck for the molecule of interest.     

QM/MM Studies on Cyclodextrin-Alcohol Interaction

The main objective of this study is to provide an insight into the interactions involved during adsorption of the alcohols on β-CD composite nanostructured membrane. Interactions between β-cyclodextrin (β-CD) and alcohols (methanol, ethanol and butanol) are studied using the QM/MM method. Magnitude of interaction energies show that the alcohols are adsorbed on the membrane. In addition, the thermochemical analysis suggests that the formation of these host-guest complexes is enthalpy driven.     

QM/MM Methods for Biomolecular Systems

Two are better than one : Quantum mechanics/molecular mechanics (QM/MM) methods are the state‐of‐the‐art computational technique for treating reactive and other “electronic” processes in biomolecular systems. This Review presents the general methodological aspects of the QM/MM approach, its use within optimization and simulation techniques, and its areas of application, always with a biomolecular focus.

     

MM3 MODELING OF ALDOPENTOSE PYRANOSE RINGS

MM3 (version 1992, ?=3.0) was used to study the ring conformations of d-xylopyranose, d-lyxopyranose and d-arabinopyranose. The energy surfaces exhibit low-energy regions corresponding to chair and skew forms with high-energy barriers between these regions corresponding to envelope and half-chair forms. The lowest energy conformer is ⁴ C ₁ for α- and β-xylopyranose and α- and β-lyxopyranose, and the lowest energy conformer is ¹ C ₄ for α- and β-arabinopyranose. Only α-lyxopyranose exhibits a secondary low-energy region (¹ C ₄) within 1 kcal/mol of its global minimum. Overall, the results are in good agreement with NMR and crystallographic results. For many of these molecules, skew conformations are found with relatively low energies (2.5 to 4 kcal/mol above lowest energy chair form). The ² S _O and ¹ C ₄ conformers of crystalline benzoyl derivatives of xylopyranose are in secondary low-energy regions on the β-xylopyranose surface, within 3.8 kcal/mol of the global ⁴ C ₁ minimum.     

Binding affinities by alchemical perturbation using QM/MM with a large QM system and polarizable MM model

The most general way to improve the accuracy of binding‐affinity calculations for protein–ligand systems is to use quantum‐mechanical (QM) methods together with rigorous alchemical‐perturbation (AP) methods. We explore this approach by calculating the relative binding free energy of two synthetic disaccharides binding to galectin‐3 at a reasonably high QM level (dispersion‐corrected density functional theory with a triple‐zeta basis set) and with a sufficiently large QM system to include all short‐range interactions with the ligand (744–748 atoms). The rest of the protein is treated as a collection of atomic multipoles (up to quadrupoles) and polarizabilities. Several methods for evaluating the binding free energy from the 3600 QM calculations are investigated in terms of stability and accuracy. In particular, methods using QM calculations only at the endpoints of the transformation are compared with the recently proposed non‐Boltzmann Bennett acceptance ratio (NBB) method that uses QM calculations at several stages of the transformation. Unfortunately, none of the rigorous approaches give sufficient statistical precision. However, a novel approximate method, involving the direct use of QM energies in the Bennett acceptance ratio method, gives similar results as NBB but with better precision, ～3 kJ/mol. The statistical error can be further reduced by performing a greater number of QM calculations. © 2015 Wiley Periodicals, Inc.     

Comparative performance of MM3(92) and two TINKER MM3 versions for the modeling of carbohydrates

The 1992 version of MM3 was largely used for modeling mono-, di-, and trisaccharides. In later versions of MM3 improvements were made in some parameters that may be important for carbohydrates. This corrected MM3 force field is part of the Tinker package, freely available (as its 4.1 version), and included in the Chem 3D Ultra 8.0 package (as the 3.7 version). The latter version lacks the corrections to the standard bond lengths produced by electronegativity and anomeric effects, whereas the Tinker 4.1 version only lacks the latter correction. The present work compares the performance of the three MM3 versions (and in some cases, DFT and/or HF/ab initio procedures) on several carbohydrate model problems as the chair and rotamer equilibria in 2-hydroxy- and 2-methoxytetrahydropyran, hydrogen bonding in cis-2,3-dihydroxytetrahydropyran, and the potential energy surfaces around the glycosidic bonds of two sulfated disaccharides and two trisaccharides. Tinker MM3 can be used accurately to estimate carbohydrate energies and geometries, and-with the help of some programming-to pursue studies on the potential energy surfaces of di- and trisaccharides. In most cases results obtained using the three MM3 versions are similar, although large energy differences are obtained when comparing a rotameric distribution around a O-C-O-H dihedral, which is almost forced to the exo-anomeric position by the Tinker versions. In other systems smaller energy differences are found, but they can nevertheless lead to a different global minimum when comparing conformers of similar energy. MM3(92) establishes better the differences between the bond lengths in both anomers, as an expected expression of the anomeric correction.     

Chaperoned alchemical free energy simulations: a general method for QM, MM, and QM/MM potentials

A general method for alchemical free energy simulations using QM, MM, and QM/MM potential is developed by introducing "chaperones" to restrain the structures, particularly near the end points. A calculation of the free energy difference between two triazole tautomers in aqueous solution is used to illustrate the method.     

Ab initio and molecular mechanics (MM2 and MM3) calculations of nonconjugated positively charged nitrogen-containing compounds

High-level ab initio calculations have been performed on N-methyl-N-methyleneammonium and related compounds to obtain accurate rotational barriers, structures, and vibrational frequencies. The 6-31G** basis set has been utilized at the Hartree-Fock level of theory for these calculations because little experimental data are available. The MM2(91) and MM3(94) force fields have been parameterized to include these nonconjugated charged nitrogen-containing compounds. Molecular mechanics geometries and vibrational frequencies compare well with the ab initio results. © 1995 John Wiley & Sons, Inc.     

Ab initio and molecular mechanics (MM2 and MM3) calculations of positively charged conjugated nitrogen-containing compounds

Ab initio calculations have been carried out on s-trans-N-vinylmethyleneammonium, pyridinium, and related compounds to obtain rotational barriers, structures, and vibrational frequencies. The restricted Hartree-Fock (RHF) level of theory with 6-31G** basis set was used for these calculations. In addition, the MM2(91) and MM3(94) force fields have been parameterized to calculate these positively charged nitrogen-containing compounds. A bond order term was incorporated in the force field to reproduce accurately the rotational barriers of s-trans-N-vinylmethyleneammonium and related compounds. Molecular mechanics geometries and vibrational frequencies compare well with those calculated by ab initio methods. © 1996 by John Wiley & Sons, Inc.     

Electronically-coupled MM quadruply-bonded complexes of molybdenum and tungsten

The reaction of M2(O2CBu(t))4 (M = Mo, W) with a dicarboxylic acid in toluene yields compounds of general formula [M2]-O2C-X-CO2-[M2] ([M2] = M2(O2CBu(t))3; X = conjugated spacer). The M2 units are electronically coupled via interactions between the M2 delta and dicarboxylate pi* orbitals, and the magnitude of this coupling is revealed by electronic structure calculations and spectroscopic data. These compounds show intense metal to ligand charge transfer (MLCT) absorptions in the visible region of the electronic spectrum that are temperature and solvent dependent. Evidence of electronic coupling is seen in their cyclic voltammograms, which show two successive one-electron oxidations. The extent of electronic coupling in the mixed valence radical cations [M2]-O2C-X-CO2-[M2]+, generated by oxidation with one equivalent of AgPF6 or FeCp2PF6, is evaluated by EPR and UV-vis-NIR spectroscopic data, and delocalized behavior is observed in compounds with W2 units separated by up to 13.6 angstroms. The simplicity of the frontier M2 orbital interactions with the bridge pi orbitals provides a convenient system with which to study electron transfer in mixed valence systems, as compared to the extensively studied, but more complicated, dinuclear t(2g)6/t(2g)5 mixed valence compounds. Oligomeric and polymeric compounds incorporating M2 units have also been synthesized, having general formula [M2(O2CR)2(O2C-Thio-CO2)]n (Thio = n-hexyl substituted ter- and quinque-thiophenes). They can be deposited as thin films by spin coating, and show photoluminescence and electroluminescence. These metallo-polythiophenes show potential for application in electronic materials. (