Informational and File Sharing ServerModeling of Macromolecular Systems | |
Tutorial 2: Including ligands in protein simulation Basic instructionsYou may want to include ligands in your protein simulation (such as cofactor, substrate or inhibitor). In this case you should load extra ligand.lib (atom types, charges and connectivity) and ligand.frcmod (additional force field parameters) files into tleap. Notice, that atomic names and atomic order in the input protein_lig.pdb file must be the same as in ligand.lib. let's assume that your protein structure contains bound cofactor and substrate. First of all you should inspect Amber parameter database and RESP ESP charge DDataBase for your ligand. Let`s suppose that you`ve found parameters for cofactor but failed with substrate. So you only need to parametrize your substrate (see next section). Once all necessary files are prepared you could load them in tleap to allow inclusion of ligands in your system: tleap -s -f leaprc.ff99SB> loadamberparams cofactor.frcmod > loadamberparams substrate.frcmod > loadoff cofactor.lib > loadoff substrate.lib > mol = loadpdb protein_lig.pdb > solvatebox mol TIP3PBOX 12.0 > addions mol Na+ 0 > saveamberparm mol protein_lig.prmtop protein_lig.inpcrd > quit |