Force-fields CHARMM

The AMBER and CHARMM force fields are best suited for protein and nucleic acid studies. 

ChemSketch has some special-purpose building functions. The peptide builder creates a line structure from the protein sequence defined with the typical three-letter abbreviations. The carbohydrate builder creates a structure from a text string description of the molecule. The nucleic acid builder creates a structure from the typical one-letter abbreviations. There is a function to clean up the shape of the structure (i.e., make bond lengths equivalent). There is also a three-dimensional optimization routine, which uses a proprietary modification of the CHARMM force field. It is possible to set the molecule line drawing mode to obey the conventions of several different publishers. 

HyperChem uses the improper dihedral angle formed by central atom - neighbor 1 - neighbor 2 - neighbor 3, where the order of neighbors is how they appear in a HIN file. Not all planar atoms customarily have associated improper torsions. The order of atoms is arbitrary but has been consistently chosen by the original authors of the CHARMM force field. The templates contain equivalent CHARMM definitions of improper torsions for amino acids. Improper dihedral angles cannot be defined that do not have a central atom, as is sometimes done in CHARMM calculations. 

The BIOh- force field option in HyperChem has no hydrogen bonding term. This is consistent with evolution and common use of the CHARMM force field (even the 1983 paper did not use a hydrogen bonding term in its example calculations and mentioned that the functional form used then was unsatisfactory and under review). 

To date, a number of simulation studies have been performed on nucleic acids and proteins using both AMBER and CHARMM. A direct comparison of crystal simulations of bovine pancreatic trypsin inliibitor show that the two force fields behave similarly, although differences in solvent-protein interactions are evident [24]. Side-by-side tests have also been performed on a DNA duplex, showing both force fields to be in reasonable agreement with experiment although significant, and different, problems were evident in both cases [25]. It should be noted that as of the writing of this chapter revised versions of both the AMBER and CHARMM nucleic acid force fields had become available. Several simulations of membranes have been performed with the CHARMM force field for both saturated [26] and unsaturated [27] lipids. The availability of both protein and nucleic acid parameters in AMBER and CHARMM allows for protein-nucleic acid complexes to be studied with both force fields (see Chapter 20), whereas protein-lipid (see Chapter 21) and DNA-lipid simulations can also be performed with CHARMM. 

It should also be noted that a force field for a wide variety of small molecules, CHARMm (note the small m, indicating the commercial version of the program and parameters), is available [39] and has been applied to protein simulations with limited success. Efforts are currently under way to extend the CHARMm small molecule force field to make the nonbonded parameters consistent with those of the CHARMM force fields, thereby allowing for a variety of small molecules to be included in computational smdies of biological systems. 

Finally, the parametrization of the van der Waals part of the QM-MM interaction must be considered. This applies to all QM-MM implementations irrespective of the quantum method being employed. From Eq. (9) it can be seen that each quantum atom needs to have two Lennard-Jones parameters associated with it in order to have a van der Walls interaction with classical atoms. Generally, there are two approaches to this problem. The first is to derive a set of parameters, e, and G, for each common atom type and then to use this standard set for any study that requires a QM-MM study. This is the most common aproach, and the derived Lennard-Jones parameters for the quantum atoms are simply the parameters found in the MM force field for the analogous atom types. For example, a study that employed a QM-MM method implemented in the program CHARMM [48] would use the appropriate Lennard-Jones parameters of the CHARMM force field [52] for the atoms in the quantum region. 

CHARMM force field https //rxsecure.umaryland.edu/research/amackere/re-... 

In an attempt to aid interpretation of the IR spectrum of MbCO we decided to model the full protein by use of a hybrid quantum mechanics/molecular mechanics approach (QM/MM), to evaluate changes in the CO stretching frequency for different protein conformations. The QM/MM method used [44] combines a first-principles description of the active center with a force-field treatment (using the CHARMM force field) of the rest of the protein. The QM-MM boundary is modeled by use of link atoms (four in the heme vinyl and propionate substituents and one on the His64 residue). Our QM region will include the CO ligand, the porphyrin, and the axial imidazole (Fig. 3.13). The vinyl and propionate porphyrin substituents were not included, because we had previously found they did not affect the properties of the Fe-ligand bonds (Section 3.3.1). It was, on the other hand, crucial to include the imidazole of the proximal His (directly bonded to the... 

Classical simulations of MbCO using the CHARMM force field were performed for different tautomerization states of the distal histidine residue (His64) [33], These simulations showed that when His64 is protonated at N,5 (denoted the tautomer) it often rotates such that it exposes either the N,>—H bond or the un-protonated N atom to the CO, as depicted in Scheme 3.4. We... 

Hypercube, Inc. at http //www.hyper.com offers molecular modeling packages under the HyperChem name. HyperChem s newest version, Hyper-Chem Release 7.5, is a full 32-bit application, developed for the Windows 95, 98, NT, ME, 2000, and XP operating systems. Density Functional Theory (DFT) has been added as a basic computational engine to complement Molecular Mechanics, Semiempirical Quantum Mechanics and ab initio Quantum Mechanics. The DFT engine includes four combination or hybrid functions, such as the popular B3-LYP or Becke-97 methods. The Bio+ force field in HyperChem represents a version of the Chemistry at HARvard using Molecular Mechanics (CHARMM) force field. Release 7.5 of HyperChem updates... 

For the construction of molecular structures, a 2D formula editor is provided in combination with 3D conversion. Standard potential energy minimization is performed using the modified parameter set of the CHARMm force field [68] the conformational models are built using Monte Carlo conformational analysis together with poling as described in the next section. 

Another set of early studies came from the work of Judson and coworkers [35, 36], which emphasized using GAs for search problems on small molecules and peptides, especially cyclic peptides. A dihedral angle representation was used for the peptides with values encoded as binary strings, and the energy function used the standard CHARMM force field. Mutations were implemented as bit flips and crossovers were introduced by a cut-and-paste of the strings. The small size of the system enabled a detailed investigation of the various parameters and policies chosen. In Ref. [37], a comparison between a GA and a direct search minimization was performed and showed the advantages and weaknesses of each method. As many concepts are shared between search problems on small peptides and complete proteins, these studies have contributed to subsequent attempts on full proteins. 

The first computational examination of COMT was the MD study by Lau and Bruice. They employed the CHARMM force field to the entire system embedded within a bath of TIP3P waters and examined the coordinates every 200 fs for a duration of 1 ns. They looked for the formation of NACs, defined here as separation of less than 3.2 A between the nucleophilic oxyanion and methyl carbon of 5, along with an 0-C-S angle greater than 165°. Formation of the NACs occurred about 7.6 percent of the time, and Tyr68 aids in forcing the cofactor into a proper position for attack. Lau and Bruice conclude that desolvation and the positioning of the substrate and cofactor into proximity are the major roles played by the enzyme in... 

We conclude this section by giving a topical example of the utility of conditional averages in considering molecularly complex systems (Ashbaugh et al, 2004). We considered the RPLC system discussed above (p. 5), but without methanol n-Ci8 alkyl chains, tethered to a planar support, with water as the mobile phase. The backside of the liquid water phase contacts a dilute water vapor truncated by a repulsive wall see Fig. 1.2, p. 7. Thus, it is appropriate to characterize the system as consistent with aqueous liquid-vapor coexistence at low pressure. A standard CHARMM force-field model (MacKerell Jr. et al, 1998) is used, as are standard molecular dynamics procedures - including periodic boimdary conditions - to acquire the data considered here. Our interest is in the interface between the stationary alkyl and the mobile liquid water phases at 300 K. 

---