GROMACS Tutorial

Step One: Prepare the Topology

The protein we will be working with is the KALP model peptide, denoted KALP15, which has a sequence of: Ac-GKK(LA)4LKKA-NH2. The protocol described here is based on a system built by Kandasamy and Larson in a study of hydrophobic mismatch. The original reference can be found here.

The peptide was prepared in-house using the xLeap module of AmberTools, using ideal backbone geometry of an α-helix (φ = -60°, ψ = -40°). The .pdb file was oriented along the z-axis using editconf -princ, followed by a rotation about the y axis. Note that in GROMACS-3.3.x, the -princ option oriented the long axis of the structure (in this case, the helix axis) along the z-axis by default, but this option has changed as of GROMACS-4.0.4, which orients the long axis along the x-axis. If you want to skip the construction of this peptide, the properly oriented structure can be found here.

Execute pdb2gmx by issuing the following command:

gmx pdb2gmx -f KALP-15_princ.pdb -o KALP-15_processed.gro -ignh -ter -water spc

When prompted, choose the GROMOS96 53A6 parameter set. Choose "None" for the termini; since we have added acetyl and amide capping groups to the N- and C-termini, respectively, we do not want pdb2gmx to build the normal amine and carboxyl groups. Instead, we want pdb2gmx to add connectivity to our capping groups. The -ignh flag tells pdb2gmx to ignore the H atoms in the input. By default, xLeap gave us an all-atom structure (since the AMBER force fields use explicit hydrogen representation). Due to AMBER naming conventions, these H atoms may not have the same nomenclature as those of the GROMOS96 force field. If we tell pdb2gmx to ignore all input H atoms, it will add back only those that it needs.

Now we will need to make some alterations to the topology.

Back: Introduction

Next: The Topology

Site design and content copyright Justin Lemkul

Problems with the site? Send them to the Webmaster