Deriving atomic charges and building a force field library for a new molecule are key steps when developing a force field required for conducting structural and energy-based analysis using molecular mechanics. Derivation of popular RESP charges for a set of residues is a complex and error prone procedure, because it depends on numerous input parameters. To overcome these problems, the R.E.D. Tools (RESP and ESP charge Derive, http://q4md-forcefieldtools.org/RED/) have been developed to perform charge derivation in an automatic and straightforward way. The R.E.D. program handles chemical elements up to bromine in the periodic table. It interfaces different quantum mechanical programs employed for geometry optimization and computing molecular electrostatic potential(s), and performs charge fitting using the RESP program. By defining tight optimization criteria and by controlling the molecular orientation of each optimized geometry, charge values are reproduced at any computer platform with an accuracy of 0.0001 e. The charges can be fitted using multiple conformations, making them suitable for molecular dynamics simulations. R.E.D. allows also for defining charge constraints during multiple molecule charge fitting, which are used to derive charges for molecular fragments. Finally, R.E.D. incorporates charges into a force field library, readily usable in molecular dynamics computer packages. For complex cases, such as a set of homologous molecules belonging to a common family, an entire force field topology database is generated. Currently, the atomic charges and force field libraries have been developed for more than fifty model systems and stored in the RESP ESP charge DDataBase. Selected results related to non-polarizable charge models are presented and discussed.

Molecular dynamics simulations describing the solvation process of native and modified cyclodextrins (per-substituted α-, β-, and γ-cyclodextrins, as well as an amino-acid derived β-cyclodextrin) have been performed. A homogeneous force field, namely "q4md-CD", has been built from the development of a new force field topology database and from a combination of the GLYCAM04 and Amber99SB force fields to correctly describe the geometrical, structural, dynamical and hydrogen bonding aspects of heterogeneous cyclodextrin based systems. These include native, organo- and peptidic-linked cyclodextrins. q4md-CD features: (i) geometrical parameters from Amber99SB to describe the protein parts, (ii) geometrical parameters from GLYCAM04 for the carbohydrate and organic parts when available or those of Amber99SB otherwise, (iii) partial atomic charges, embedded in force field libraries for the carbohydrate and organic fragments, were derived using the R.E.D. tools according to the "Amber" strategy and (iv) scaling factors of 1.2 and 2.0 were imposed for the 1-4 electrostatic and 1-4 van der Waals interactions, respectively. Results given by q4md-CD on native cyclodextrins have been compared to those obtained with reference to force fields like GLYCAM04, GLYCAM06 and Amber99SB as well as with experimental data. This work not only gives a global view of the performances of the aforementioned force fields towards a correct description of solvated cyclodextrins, but also extends the capabilities of current force fields by addressing some issues concerning hydrogen bonding and opens new possibilities towards studies of glycoconjugates by molecular dynamics.

The web-based RESP ESP charge DataBase (R.E.DD.B., http://q4md-forcefieldtools.org/REDDB) is a free and new source of RESP and ESP atomic charge values and force field libraries for model systems and/or small molecules. R.E.DD.B. stores highly effective and reproducible charge values and molecular structures in the Tripos mol2 file format, information about the charge derivation procedure, scripts to integrate the charges and molecular topology in the most common molecular dynamics packages. Moreover, R.E.DD.B. allows users to freely store and distribute RESP or ESP charges and force field libraries to the scientific community, via a web interface. The first version of R.E.DD.B., released in January 2006, contains force field libraries for molecules as well as molecular fragments for standard residues and their analogs (amino acids, monosaccharides, nucleotides and ligands), hence covering a vast area of relevant biological applications.

We present the first comparative molecular dynamics investigation for a dodecylphosphocholine (DPC) micelle performed in condensed phase using the CHARMM36, GROMOS53A6, GROMOS54A7, and GROMOS53A6/Berger force fields and a set of parameters developed anew. Our potential consists of newly derived RESP atomic charges, which are associated with the Amber99SB force field developed for proteins. This new potential is expressly designed for simulations of peptides and transmembrane proteins in a micellar environment. To validate this new ensemble, molecular dynamics simulations of a DPC micelle composed of 54 monomers were carried out in explicit water using a "self-assembling" approach. Characteristic micellar properties such as aggregation kinetic, volume, size, shape, surface area, internal structure, surfactant conformation, and hydration were thoroughly examined and compared with experiments. Derived RESP charge values combined with parameters taken from Amber99SB reproduce reasonably well important structural properties and experimental data compared to the other tested force fields. However, the headgroup and alkyl chain conformations or the micelle hydration simulated with the Amber99SB force field display some differences. In particular, we show that Amber99SB slightly overestimates the trans population of the alkyl Csp(3)-Csp(3)-Csp(3)-Csp(3) dihedral angle (i.e., CCCC) and reduces the flexibility of the DPC alkyl chain. In agreement with experiments and previously published studies, the DPC micelle shows a slightly ellipsoidal shape with a radius of gyration of ∼17 Å for the different potentials evaluated. The surface of contact between the DPC headgroup and water molecules represents between 70% and 80% of the total micelle surface independently of the force field considered. Finally, molecular dynamics simulations show that water molecules form various hydrogen-bond patterns with the surfactant headgroup, as noted previously for phospholipids with a phosphatidylcholine headgroup.