Rapid parameterization of small molecules using the force field toolkit

The inability to rapidly generate accurate and robust parameters for novel chemical matter continues to severely limit the application of molecular dynamics simulations to many biological systems of interest, especially in fields such as drug discovery. Although the release of generalized versions of common classical force fields, for example, General Amber Force Field and CHARMM General Force Field, have posited guidelines for parameterization of small molecules, many technical challenges remain that have hampered their wide‐scale extension. The Force Field Toolkit (ffTK), described herein, minimizes common barriers to ligand parameterization through algorithm and method development, automation of tedious and error‐prone tasks, and graphical user interface design. Distributed as a VMD plugin, ffTK facilitates the traversal of a clear and organized workflow resulting in a complete set of CHARMM‐compatible parameters. A variety of tools are provided to generate quantum mechanical target data, setup multidimensional optimization routines, and analyze parameter performance. Parameters developed for a small test set of molecules using ffTK were comparable to existing CGenFF parameters in their ability to reproduce experimentally measured values for pure‐solvent properties (<15% error from experiment) and free energy of solvation (±0.5 kcal/mol from experiment). © 2013 Wiley Periodicals, Inc.     

CHARMM‐GUI ligand reader and modeler for CHARMM force field generation of small molecules

Reading ligand structures into any simulation program is often nontrivial and time consuming, especially when the force field parameters and/or structure files of the corresponding molecules are not available. To address this problem, we have developed Ligand Reader & Modeler in CHARMM‐GUI. Users can upload ligand structure information in various forms (using PDB ID, ligand ID, SMILES, MOL/MOL2/SDF file, or PDB/mmCIF file), and the uploaded structure is displayed on a sketchpad for verification and further modification. Based on the displayed structure, Ligand Reader & Modeler generates the ligand force field parameters and necessary structure files by searching for the ligand in the CHARMM force field library or using the CHARMM general force field (CGenFF). In addition, users can define chemical substitution sites and draw substituents in each site on the sketchpad to generate a set of combinatorial structure files and corresponding force field parameters for throughput or alchemical free energy simulations. Finally, the output from Ligand Reader & Modeler can be used in other CHARMM‐GUI modules to build a protein‐ligand simulation system for all supported simulation programs, such as CHARMM, NAMD, GROMACS, AMBER, GENESIS, LAMMPS, Desmond, OpenMM, and CHARMM/OpenMM. Ligand Reader & Modeler is available as a functional module of CHARMM‐GUI at http://www.charmm-gui.org/input/ligandrm . © 2017 Wiley Periodicals, Inc.     

A CGenFF-based force field for simulations of peptoids with both cis and trans peptide bonds

Peptoids, or poly-n-substituted glycines, are peptide-like polymers composed of a flexible backbone decorated with diverse chemical side chains. Peptoids can form a variety of self-assembling structures based on the type and sequence of the side chains attached to their backbones. All-atom molecular dynamics simulations have been useful in predicting the conformational structures of proteins and will be valuable tools for identifying combinations of peptoid side chains that may form interesting folded structures. However, peptoid models must address a major degree of freedom not common in proteins – the cis/trans isomerization of the peptide bond. This work presents CHARMM general force field (CGenFF) parameters developed to accurately represent peptoid conformational behavior, with an emphasis on a correct representation of both the cis and trans isomers of the peptoid backbone. These parameters are validated against experimental and quantum mechanics data and used to simulate three peptoid side chains in explicitly solvated systems. © 2019 Wiley Periodicals, Inc.