一个使用gromacs进行蛋白质模拟的入门教程.docx

一个使用gromacs进行蛋白质模拟的入门教程.docx

《一个使用gromacs进行蛋白质模拟的入门教程.docx》由会员分享，可在线阅读，更多相关《一个使用gromacs进行蛋白质模拟的入门教程.docx（16页珍藏版）》请在冰豆网上搜索。

一个使用gromacs进行蛋白质模拟的入门教程

GROMACSTutorial

LysozymeinWater

JustinLemkul

DepartmentofBiochemistry,VirginiaTech

Thisexamplewillguideanewuserthroughtheprocessofsettingupasimulationsystemcontainingaprotein（lysozyme）inaboxofwater,withions.Eachstepwillcontainanexplanationofinputandoutput,usingtypicalsettingsforgeneraluse.

ThistutorialassumesyouareusingaGROMACSversioninthe4.5.xseries.

GROMACSTutorial

StepOne:

PreparetheTopology

Wemustdownloadtheproteinstructurefilewewillbeworkingwith.Forthistutorial,wewillutilizeheneggwhitelysozyme（PDBcode1AKI）.GototheRCSBwebsiteanddownloadthePDBtextforthecrystalstructure.

Onceyouhavedownloadedthestructure,youcanvisualizethestructureusingaviewingprogramsuchasVMD,Chimera,PyMOL,etc.Onceyou'vehadalookatthemolecule,youaregoingtowanttostripoutthecrystalwaters.Useaplaintexteditorlikevi,emacs（Linux/Mac）,orNotepad（Windows）.Donotusewordprocessingsoftware!

Deletethelinescorrespondingtothesemolecules（residue"HOH"inthePDBfile）.Notethatsuchaprocedureisnotuniversallyappropriate（i.e.,thecaseofaboundactivesitewatermolecule）.Forourintentionshere,wedonotneedcrystalwater.

Alwayscheckyour.pdbfileforentrieslistedunderthecommentMISSING,astheseentriesindicateeitheratomsorwholeresiduesthatarenotpresentinthecrystalstructure.Terminalregionsmaybeabsent,andmaynotpresentaproblemfordynamics.Incompleteinternalsequencesoranyaminoacidresiduesthathavemissingatomswillcausepdb2gmxtofail.Thesemissingatoms/residuesmustbemodeledinusingothersoftwarepackages.Alsonotethatpdb2gmxisnotmagic.Itcannotgeneratetopologiesforarbitrarymolecules,justtheresiduesdefinedbytheforcefield（inthe*.rtpfiles-generallyproteins,nucleicacids,andaveryfiniteamountofcofactors,likeNAD（H）andATP）.

Nowthatthecrystalwatersaregoneandwehaveverifiedthatallthenecessaryatomsarepresent,thePDBfileshouldcontainonlyproteinatoms,andisreadytobeinputintothefirstGROMACStool,pdb2gmx.Thepurposeofpdb2gmxistogeneratethreefiles:

1.Thetopologyforthemolecule.

2.Apositionrestraintfile.

3.Apost-processedstructurefile.

Thetopology（topol.topbydefault）containsalltheinformationnecessarytodefinethemoleculewithinasimulation.Thisinformationincludesnonbondedparameters（atomtypesandcharges）aswellasbondedparameters（bonds,angles,anddihedrals）.Wewilltakeamoredetailedlookatthetopologyonceithasbeengenerated.

Executepdb2gmxbyissuingthefollowingcommand:

pdb2gmx-f1AKI.pdb-o1AKI_processed.gro-waterspce

Thestructurewillbeprocessedbypdb2gmx,andyouwillbepromptedtochooseaforcefield:

SelecttheForceField:

From'/usr/local/gromacs/share/gromacs/top':

1:

AMBER03forcefield（Duanetal.,J.Comp.Chem.24,1999-2012,2003）

2:

AMBER94forcefield（Cornelletal.,JACS117,5179-5197,1995）

3:

AMBER96forcefield（Kollmanetal.,Acc.Chem.Res.29,461-469,1996）

4:

AMBER99forcefield（Wangetal.,J.Comp.Chem.21,1049-1074,2000）

5:

AMBER99SBforcefield（Hornaketal.,Proteins65,712-725,2006）

6:

AMBER99SB-ILDNforcefield（Lindorff-Larsenetal.,Proteins78,1950-58,2010）

7:

AMBERGSforcefield（Garcia&Sanbonmatsu,PNAS99,2782-2787,2002）

8:

CHARMM27all-atomforcefield（withCMAP）-version2.0

9:

GROMOS9643a1forcefield

10:

GROMOS9643a2forcefield（improvedalkanedihedrals）

11:

GROMOS9645a3forcefield（SchulerJCC2001221205）

12:

GROMOS9653a5forcefield（JCC2004vol25pag1656）

13:

GROMOS9653a6forcefield（JCC2004vol25pag1656）

14:

OPLS-AA/Lall-atomforcefield（2001aminoaciddihedrals）

15:

[DEPRECATED]Encadall-atomforcefield,usingfullsolventcharges

16:

[DEPRECATED]Encadall-atomforcefield,usingscaled-downvacuumcharges

17:

[DEPRECATED]Gromacsforcefield（seemanual）

18:

[DEPRECATED]GromacsforcefieldwithhydrogensforNMR

Theforcefieldwillcontaintheinformationthatwillbewrittentothetopology.Thisisaveryimportantchoice!

Youshouldalwaysreadthoroughlyabouteachforcefieldanddecidewhichismostapplicabletoyoursituation.Forthistutorial,wewillusetheall-atomOPLSforcefield,sotype14atthecommandprompt,followedby'Enter'.

Therearemanyotheroptionsthatcanbepassedtopdb2gmx.Somearelistedhere:

-ignh:

IgnoreHatomsinthePDBfile;especiallyusefulforNMRstructures.Otherwise,ifHatomsarepresent,theymustbeinthecorrectorderandnamedexactlyhowGROMACSexpectsthemtobe.

-ter:

InteractivelyassignchargestatesforN-andC-termini.

-inter:

InteractivelyassignchargestatesforGlu,Asp,Lys,Arg,andHis;assigndisulfidestoCys.

Youhavenowgeneratedthreenewfiles:

1AKI_pr