Automatic molecular structure perception for the universal force field

The Universal Force Field (UFF) is a classical force field applicable to almost all atom types of the periodic table. Such a flexibility makes this force field a potential good candidate for simulations involving a large spectrum of systems and, indeed, UFF has been applied to various families of molecules. Unfortunately, initializing UFF, that is, performing molecular structure perception to determine which parameters should be used to compute the UFF energy and forces, appears to be a difficult problem. Although many perception methods exist, they mostly focus on organic molecules, and are thus not well‐adapted to the diversity of systems potentially considered with UFF. In this article, we propose an automatic perception method for initializing UFF that includes the identification of the system's connectivity, the assignment of bond orders as well as UFF atom types. This perception scheme is proposed as a self‐contained UFF implementation integrated in a new module for the SAMSON software platform for computational nanoscience ( http://www.samson-connect.net ). We validate both the automatic perception method and the UFF implementation on a series of benchmarks.     

Development and testing of a general amber force field

We describe here a general Amber force field (GAFF) for organic molecules. GAFF is designed to be compatible with existing Amber force fields for proteins and nucleic acids, and has parameters for most organic and pharmaceutical molecules that are composed of H, C, N, O, S, P, and halogens. It uses a simple functional form and a limited number of atom types, but incorporates both empirical and heuristic models to estimate force constants and partial atomic charges. The performance of GAFF in test cases is encouraging. In test I, 74 crystallographic structures were compared to GAFF minimized structures, with a root-mean-square displacement of 0.26 A, which is comparable to that of the Tripos 5.2 force field (0.25 A) and better than those of MMFF 94 and CHARMm (0.47 and 0.44 A, respectively). In test II, gas phase minimizations were performed on 22 nucleic acid base pairs, and the minimized structures and intermolecular energies were compared to MP2/6-31G* results. The RMS of displacements and relative energies were 0.25 A and 1.2 kcal/mol, respectively. These data are comparable to results from Parm99/RESP (0.16 A and 1.18 kcal/mol, respectively), which were parameterized to these base pairs. Test III looked at the relative energies of 71 conformational pairs that were used in development of the Parm99 force field. The RMS error in relative energies (compared to experiment) is about 0.5 kcal/mol. GAFF can be applied to wide range of molecules in an automatic fashion, making it suitable for rational drug design and database searching.     

Peptide free‐energy profile is strongly dependent on the force field: Comparison of C96 and AMBER95

The C96 and AMBER95 force fields were compared with small model peptides Ac‐(Ala)_n‐NMe (Ac = CH₃CO, NMe = NHCH₃, n=2 and 3) in vacuo and in TIP3P water by computing the free‐energy profiles using multicanonical molecular dynamics method. The C96 force field is a modified version of the AMBER95 force field, which was adjusted to reproduce the energy difference between extended β‐ and constrained α‐helical energies for the alanine tetrapeptide, obtained by the high level ab initio MO method. The slight modification resulted in a large difference in the free energy profiles. The C96 force field prefers relatively extended conformers, whereas the AMBER95 force field favors turn conformations. © 2000 John Wiley & Sons, Inc. J Comput Chem 21: 748–762, 2000     

All-atom and united-atom simulations of guanidinium-based ionic liquids

Ionic liquids (ILs) have been widely used in separation, catalysis, electrochemistry, etc., and one of the most outstanding characteristics is that ILs can be tailored and tuned for specific tasks. In order to design and make better use of ionic liquids, the structures and properties relationship is indispensable. Both molecular dynamics and Monte Carlo simulations have been proved useful to understand the behavior of molecules at the microscale and the properties of the system. However, the quality of such simulations depends on force field parameters describing the interactions between atoms. All-atom (AA) or the united-atom (UA) force fields will be chosen because of the demand for more exact results or the lower computational cost, respectively. In order to make a systematic comparison of the two force fields, molecular simulations for four kinds of acyclic guanidinium-based ionic liquids (cations: (R2N)2C=N+<, anion: nitric or perchloric acid) were performed based on the AA and the UA force fields in this work. AA force field parameters were derived from our previous work (Fluid Phase Equilib., 2008, 272: 1-7), and the UA parameters were proposed in this work. Molecular dynamics simulation results for the AA and UA force fields were compared. Simulation densities are very similar to each other. Center of mass radial distribution functions (RDFs), site to site RDFs and spatial distribution functions (SDFs) were also investigated to depict the microscopic structures of the ILs.     

Parametrization of macrolide antibiotics using the force field toolkit

Macrolides are an important class of antibiotics that target the bacterial ribosome. Computer simulations of macrolides are limited as specific force field parameters have not been previously developed for them. Here, we determine CHARMM‐compatible force field parameters for erythromycin, azithromycin, and telithromycin, using the force field toolkit (ffTK) plugin in VMD. Because of their large size, novel approaches for parametrizing them had to be developed. Two methods for determining partial atomic charges, from interactions with TIP3P water and from the electrostatic potential, as well as several approaches for fitting the dihedral parameters were tested. The performance of the different parameter sets was evaluated by molecular dynamics simulations of the macrolides in ribosome, with a distinct improvement in maintenance of key interactions observed after refinement of the initial parameters. Based on the results of the macrolide tests, recommended procedures for parametrizing very large molecules using ffTK are given. © 2015 Wiley Periodicals, Inc.     

ForceFit: A code to fit classical force fields to quantum mechanical potential energy surfaces

The ForceFit program package has been developed for fitting classical force field parameters based upon a force matching algorithm to quantum mechanical gradients of configurations that span the potential energy surface of the system. The program, which runs under UNIX and is written in C++, is an easy‐to‐use, nonproprietary platform that enables gradient fitting of a wide variety of functional force field forms to quantum mechanical information obtained from an array of common electronic structure codes. All aspects of the fitting process are run from a graphical user interface, from the parsing of quantum mechanical data, assembling of a potential energy surface database, setting the force field, and variables to be optimized, choosing a molecular mechanics code for comparison to the reference data, and finally, the initiation of a least squares minimization algorithm. Furthermore, the code is based on a modular templated code design that enables the facile addition of new functionality to the program. © 2010 Wiley Periodicals, Inc. J Comput Chem, 2010     

The quest for the best nonpolarizable water model from the adaptive force matching method

The recently introduced adaptive force matching (AFM) method is used to develop a significantly improved pair‐wise nonpolarizable potential for water. A rigid version of the potential is also presented to enable larger time steps for biological simulations. In this work, it is demonstrated that the AFM method can be used to systematically assess the importance of each functional term during the construction of a force field. For a water potential, it is established that a single off‐atom charge center (M) in the plane of water outperforms two out‐of‐plane charge sites for reproducing intermolecular forces. The four‐site pair‐wise nonpolarizable force field developed in this work rivals some of the most sophisticated polarizable models in terms of reproducing accurate ab initio forces. The force fields are parameterized to perform best in the temperature range from 0 to 40°C. Equilibrium and dynamical properties calculated with the flexible and rigid force fields are in good agreement with experimental results. For the flexible model, the agreement improves when path integral simulation is performed. These force fields provide high‐quality results at a very low computational cost and are thus well suited to atomistic scale biological simulations. The AFM method provides a mechanism for selecting important terms in force field expressions and is a very promising tool for producing accurate force fields in condensed phases. © 2010 Wiley Periodicals, Inc. J Comput Chem, 2011     

分子力场进展

分子力学（简称ＭＭ）是近年来化学家常用的一种计算方法。与量子力学从头计算和半经验方法相比，用分子力学处理大分子可以大大节省计算时间，而且，在大多数情况下，用分子力学方法计算得到的分子几何构型参数与实验值之间的差值可在实验误差范围之内。所以，分子力学是研究生物化学体系的有效和可行的手段。分子力学的核心是分子力场。本文介绍了分子力场的量子力学背景、分子力场和光谱力场之间的关系。分子力场的一般形式、分力     

Empirische Kraftfeldberechnungen an Cyclosilanen: Perhalogencyclosilane

Empirical force field calculations are used to interpret relative energies and structures of various conformations of perhalocyclotetra-, -penta- and-hexasilanes. The results are in good agreement with experimental data as far as they are known. Comparison is made with some carbon analogs.

     

Force field modeling of conformational energies: Importance of multipole moments and intramolecular polarization

The ability of force fields to reproduce relative conformational energies for seven molecules is probed. The correlation against LMP2/cc‐pVTZ results for standard force field employing fixed partial charges deteriorates as the molecules become more polar. Inclusion of multipole moments and intramolecular polarization can improve the correlation, and both contributions are of similar importance. © 2006 Wiley Periodicals, Inc. Int J Quantum Chem, 2007     

A biomolecular force field based on the free enthalpy of hydration and solvation: the GROMOS force-field parameter sets 53A5 and 53A6

Successive parameterizations of the GROMOS force field have been used successfully to simulate biomolecular systems over a long period of time. The continuing expansion of computational power with time makes it possible to compute ever more properties for an increasing variety of molecular systems with greater precision. This has led to recurrent parameterizations of the GROMOS force field all aimed at achieving better agreement with experimental data. Here we report the results of the latest, extensive reparameterization of the GROMOS force field. In contrast to the parameterization of other biomolecular force fields, this parameterization of the GROMOS force field is based primarily on reproducing the free enthalpies of hydration and apolar solvation for a range of compounds. This approach was chosen because the relative free enthalpy of solvation between polar and apolar environments is a key property in many biomolecular processes of interest, such as protein folding, biomolecular association, membrane formation, and transport over membranes. The newest parameter sets, 53A5 and 53A6, were optimized by first fitting to reproduce the thermodynamic properties of pure liquids of a range of small polar molecules and the solvation free enthalpies of amino acid analogs in cyclohexane (53A5). The partial charges were then adjusted to reproduce the hydration free enthalpies in water (53A6). Both parameter sets are fully documented, and the differences between these and previous parameter sets are discussed.     

Effect of ampere force on electrolyte flow and current passing due to EMF generation in electrochemical cell rotating in magnetic field

An Ampere force acts in the rotating electrochemical cell, which is located in the magnetic field. This force causes an electrolyte flow directed oppositely to the direction of cell rotation. The effect of Ampere force on the distribution of hydrodynamic velocity in the gap between two cylindrical electrodes of rotating cell is theoretically analyzed. Under certain simplifying assumptions, an equation is obtained for calculating the current passing in the cell by action of potential difference, which arises in the cell due to the Lorentz force taking into account the Ampere force.     

Advancing beyond the atom-centered model in additive and nonadditive molecular mechanics

A computational approach to the inclusion of off-center charges in both additive and nonadditive molecular mechanics calculations is presented. The additional sites in the molecular skeleton are placed in the approximate locations of the chemically intuitive electron lone pair, and are treated as formal particles throughout the calculation. The increase in the number of charge sites results in overall improvement in the energy associated with the angular dependence of hydrogen bonds and improved statistical accuracy of the electrostatic potential derived charges. The addition of lone pairs also results in improved accuracy in relative solvation free energy calculation for the pyridine to benzene and methanol to methane mutations. Because the number of atoms that require lone pairs is small, the extra accuracy can be achieved with little computational overhead. © 1997 John Wiley & Sons, Inc. J Comput Chem 18 : 1632–1646, 1997     

A case study in conformation design: learning by doing

Efforts are described to design simple, fully flexible but conformationally preorganised omega-hydroxy-nonanoic acids that could serve as the conformation controlling unit in analogues of the potent protein-kinase C activator aplysiatoxin. Such analogues are macrodilactones incorporating the designed omega-hydroxy-nonanoic acid and 3,4-dihydroxy-pentanoic acid, which contains the pharmacophoric groups. The design process (replacement of CH(2) groups by an oxygen atom, annelation of a six-membered ring and placement of alkyl substituents) of the omega-hydroxy-nonanoic acids was monitored by force-field calculations. In the end of this process simple analogues of aplysiatoxin are proposed in which the proper disposition of the pharmacophoric groups is secured by a conformationally flexible but preorganised template structure as part of the macrodilactone ring.     

The tetracycline: Mg2+ complex: a molecular mechanics force field

Tetracycline (Tc) is an important antibiotic, which binds specifically to the ribosome and several proteins, in the form of a Tc-:Mg2+ complex. To model Tc:protein and Tc:RNA interactions, we have developed a molecular mechanics force field model of Tc, which is consistent with the CHARMM force field for proteins and nucleic acids. We used structures from the Cambridge Crystallographic Data Base to identify the main Tc conformations that are likely to be present in solution and in biomolecular complexes. A conformational search was also done, using the MM3 force field to perform simulated annealing of Tc. Several resulting, low-energy structures were optimized with an ab initio model and used in developing the new Tc force field. Atomic charges and Lennard-Jones parameters were derived from a supermolecule ab initio approach. We considered the ab initio energies and geometries of a probe water molecule interacting with Tc at 36 different positions. We considered both a neutral and a zwitterionic Tc form, with and without bound Mg2+. The final rms deviation between the ab initio and force field energies, averaged over all forms, was just 0.35 kcal/mol. The model also reproduces the ab initio geometry and flexibility of Tc. As further tests, we did simulations of a Tc crystal, of Tc:Mg2+ and Tc:Ca2+ complexes in aqueous solution, and of a solvated complex between Tc:Mg2+ and the Tet repressor protein (TetR). With slight, ad hoc adjustments, the model can reproduce the experimental, relative, Tc binding affinities of Mg2+ and Ca2+. It performs well for the structure and fluctuations of the Tc:Mg2+:TetR complex. The model should therefore be suitable to investigate the interactions of Tc with proteins and RNA. It provides a starting point to parameterize other compounds in the large Tc family.     

Harmonic force field between the (0001) surface of graphite and adsorbed methane

Empirical potentials were used to approximate the interaction between the (0001) surface of graphite and a single methane molecule. The frequency of the vibration perpendicular to the surface was obtained by assuming that the force field under consideration is harmonic, and also by solving the one-dimensional Schrödinger equation of the vibration. It was found that the two approaches yielded nearly identical excitation energies for the fundamental transition.     

Molecular recognition in carnitine acyltransferases

We are designing and synthesizing rigid guests to probe the topography of the carnitine acyltransferases, regulatory enzymes in lipid metabolism. Our designs are based on structural studies of substrates and possible molecular mechanisms of enzymatic activity. Recent X-ray,¹H NMR, and force-field computational studies on carnitine and acetylcarnitine, coupled with the known stereospecificity for activity in carnitine acyltransferases, have led us to propose a molecular mechanism for acyl transfer in these enzymes. The folded conformation of an acylcarnitine is most populated and should be preferred for binding to these enzymes, because, in this conformation, the acyloxy is the most sterically accessible. There are four key recognition sites on the enzymes: I, carboxylate; II, trimethylammonium; III, coenzyme A; IV, acyl. Sites, I, II and III serve as the three loci required to create a chiral environment on the enzymes for carnitine. An addition-elimination reaction involving the formation of a tetrahedral intermediate is suggested as the mechanism for O-to-S acyl transfer. This proposed tetrahedral intermediate is chiral and the enzymes should prefer theR configuration at this center. Based on this proposal, conformationally rigid tetrahedral-intermediate analogues have been designed, synthesized and assayed. Morpholinium and 2-hydroxymorpholinium derivatives inhibit carnitine acetyltransferase and palmitoyltransferase. Because of rigidity at their two chiral centers, these inhibitors serve as probes of molecular topography of recognition sites, I, II, and IV.