外文翻译-- The Morphology Prediction of Lysozyme Crystals Deduced from the BFDH Law and.PDF
TheMorphologyPredictionofLysozymeCrystalsDeducedfromtheBFDHLawandAttachmentEnergyModelBasedontheIntermolecularInteractionZhanzhongWang,PingpingJiangSchoolofAgricultureandBioengineeringTianjinUniversityTianjin300072,PeoplesRepublicofChinawzz7698tju.edu.cnLepingDangSchoolofChemicalEngineeringandTechnologyTianjinUniversityTianjin300072,PeoplesRepublicofChinadlepinghotmail.comAbstractThecrystalmorphologyoforthorhombiclysozymeispredictedusingtheBravais-Friedel-Donnay-Harker(BFDH)andtheattachmentenergy(AE)modelsofmolecularsimulationsoftwareCerius2invacuo.Themorphologypredictedbytwomodelsisapproximatelyconsistent.ThemorphologypredictedbyAEmodelisingoodagreementwiththemorphologyofcrystalsgrownfromsolutionatpH6.5.Themaincrystalfaces011,101and110areobservedinmorphologypredictedbyAEmodel.BycleavingrevealablecrystalfacesinmorphologypredictedbyAEmodel,surfacechemistryvisualizationandtheoreticalanalysisbasedoninteractionofinintra-moleculesorinter-moleculesfortheimportantmorphologicalformsareperformed.TheresultshowsthatsterichindranceandH-bandinteractionplayscriticalrolefortheplate-likemorphologyoforthorhombiclysozyme.Keywords-lysozyme;morphologypredcition;modeling;IntermolecularInteractionINTRODUCTIONTheshapeofacrystalisgovernedbytherelativegrowthratesofeachofthecrystalfacespresent.Themostprominentfaceofacrystalistheslowestgrowing,whilethesmallestfaceisthefastestgrowing.Thus,themorphologicalimportanceofaparticularfaceisinverselyproportionaltothegrowthrateofcrystalface.Identificationofpossiblegrowthdirectionsorfacesisthefirstpartofmorphologyprediction1.EarlyworkincrystalmorphologypredictionwasundertakenbyGibbsin1875inwhichheproposedthattheshapeofacrystalwillbeonetominimizethetotalfreeenergyassociatedwiththesurfaceenergiesofthegrowncrystal.LaterBravais,Friedel,andDannayandHarker1,2proposedmorphologicalsimulationsbasedmerelyoncrystallatticegeometry.TheBFDHlaw,substantiatedbyempiricalobservation,proposedthatthecrystalwillbeformedbyfacesboundingthedirectionsofthegreatestinterplanarspacingofthecrystallattice.Accordingtothislaw,thelargertheinterplanardistancedhklis,thelargerthemorphologicalimportance(MI)ofthecorrespondinghklfaceis.Thus,basedonunitcellparameters,spacegroupandhence,extinctioncondition,andthereciprocaloftheinterplanardistances1/dhkl,usedasthedistancesfromthecentreofthecrystaltotherespectivehklsurfaces,onemayobtainatheoreticalBFDHmorphology.TheBFDHmethodestimatesthemorphologyfromthecrystalsymmetryandthelatticeparameterswithouttakingintoaccountthechemicalnatureandpackingoftheatomsormoleculesthatformthecrystalandneglectingthepossibilityofspecificenergeticinteractionsbetweenthesurfaceatomsinfluencingthecrystalmorphology.Itmaybeusedasaquickscreeningtooltoestimatethehabitofacrystal.In1955,HartmanandPerdokproposedperiodic-bondchain(PBC)theorytoaccountforenergyinteractionsbetweencrystallizingspeciesinthederivationofcrystalmorphology.LaterthistheorywasfurtherdevelopedbyHartmanandHartmanandBennema3.APBCisastrongbondbetweenmoleculesrunningparalleltoacrystallographicdirection.AcrystalismadeofanetworkofPBCshavingdifferentenergiesdependingonthebondenergybetweenthemoleculestheycontain.Consequently,itispossibletoclassifythePBCsasafunctionoftheirenergyanddistinguishstrongPBCsfromweakPBCs4.Growthrateisrelatedtodirectionofbondchainandthefastestdirectionofcrystalgrowthisthatofthestrongestchemicalbond.In1980,basedonPBCmodel,HartmanandPerdokproposedattachmentenergymodel.Buildonincorporatingenergytermsinthecalculationofthetheoreticalcrystalmorphology,advancesinmorphologypredictionwererealizedwhencrystalgrowthwasviewedastheattachmentofslicesorlayersoforderedmoleculestoanexistingcrystalface1.Thisconcept,knownastheattatchmentenergy(AE)method,requiresthedeterminationofthelatticeandsliceenergiesforcalculationoftheattachmentenergyineachimportantcrystallographicdirection.Thisattachmentenergydirectlyinfluencesthemacroscopicshapeofthecrystal4.Thelatticeenergyistheenergycalculatedbysummingpairinteractionsfortheentire,perfectcrystal.Theenergyreleasedontheformationofagrowthsliceofthickness,dhkl,isthesliceenergy,whiletheenergyreleasedontheattachmentofagrowthslicetothecrystalsurfaceistheattachmentenergy.Thetermofinterestformorphologicalpredictionsistheattachmentenergyofeachcrystalface.Subtractingthesliceenergyforeachfacefromthelatticeenergyprovidestheattachmentenergyforeachofthefaces.Therelativegrowthrateforeachfacewastakentobeproportionaltoattachmentenergyandhencefaceswiththelowestattachmentenergieswerepredictedtobetheslowestgrowingsurfacesandhencetohavethehighestmorphologicalimportance5.Thenamecrystallogenesisarosewhenitbecameevidentthatthefieldofcrystallizationofproteinsandotherbiologicalmacromoleculeswasnotrestrictedsimplytocrystalproductionfordiffractionstudies,butitencompassed,infact,allphasesofstructuralbiology,fromproteinexpressionandpurification,torecordingofdiffractiondata.Withregardto978-1-4244-4713-8/10/$25.00©2010IEEEbiologicalmacromolecule,X-raycrystallographyhasplayedafundamentalroleinconnectingthedotsbetweengenomicdataandbiologicalfunctionbyprovidingaccuratestructuralinformationtoresolveseveralsignificantresearchproblems.TheearliestX-raydiffractionstudiesofthestructuresofbiologicalmacromoleculesbeganintheearlypartofthe20thcentury,verylittleprogresshasbeenmadeinourunderstandingofhowtofacilitatetheprocessofcrystallizingsuchmacromoleculesforstructuralanalyses.Asaresult,obtaininghigh-qualitymacromolecularcrystalsremainsdifficult,unpredictable,andfrustratingandbecomespersistentbottlenecktothegreaterapplicationofX-raycrystallographyinstructuralbiology6,7.SolvingproteinstructuresbyX-raycrystallographyiscontingentupontheavailabilityofordered,diffraction-qualitycrystals8.Inthefieldofcrystalgrowth,inordertoobtainhigh-qualitylaboratory-growncrystals,theresearcherwhoisinterestedincrystalmorphologyoftenneedstovisualizeacrystalhabitresultingfromasetofobservedorcalculateddata,suchasgrowthrates,surfaceenergiesandhabitcontrollingfactors.Crystallizationiseffectedviamolecularrecognitionattheinterfacebetweenthegrowingcrystalanditsmotherphase.Sucheffectshavelongbeenobservedandhavebeenrelatedtogrowthsolutionthermodynamics,specificfacestructureandgrowthmechanism,andmolecularrecognitionprocessatspecificplanes9.Increasedunderstandingofcrystalgrowthfromsolutioncanenhancetheperformanceoftheseparationandpurificationprocessesinmanyindustries.Althoughthecrystalgrowthhabitdiffersfromthetheoreticalmorphologyduetotheinfluenceofcrystalgrowthenvironment,thecrystalmorphologyisverymuchconnectedwiththecrystalstructure.Crystalmorphologyhasbeenthefocusofnumerousresearchefforts10-12.Crystallizationofheneggwhitelysozymehasalreadybeenstudiedforsomeyearssinceitisanidealmodelsystemforcrystallizationofproteinsingeneral.Otherkindsoflysozymecrystalshaveseldombeenused,exceptforthetheorthorhombicform13.Themorphologyanalysisoftheorthorhombiclysozymecrystalsmayincreaseourknowledgeofcrystalgrowthoflysozymetoinstructtoobtainhigh-qualitycrystals,whichcanbeappliedstructuralanlysistoprovideaccurateinformationforgenomicdataandbiologicalfunction.Theaimofthepresentresearchiswastoapplyadvancedmolecularmodellingtechniquestosimulatethetheoreticalcrystalmorphologyoforthorhombiclysozyme,andcomparesthetheoreticalmorphologywithactualhabitoflysozymecrystalgrowninanaqueoussolution.Inthecourseofrefiningthemodel,valuableinformationconcerningmolecularinteractionswithinthelysozymecrystalwasextractedfromthesimulation.EXPERIMENTALANDCOMPUTATIONALMOLECULARMODELINGMETHODOLOGYMaterialsHenegg-whitelysozyme,recrystallizedandlyopholizedwaspurchasedfromSigmaandusedwithoutfurtherpurification.Itsmolecularweightwasassumedtobe14.3kgmol-1.Otherchemicalreagentswereanalyticalpurity.DistilleddeionizedwaterofHPLCgradewasused.Crystallizationexperimentandmicroscopymeasurement.Protein-waterstocksolutionswithproteinconcentrationof10-20%wereprepared,andtheirmeasuredpHwas6.5orslightlyhigher.Ifnecessary,pHwasloweredto6.5bysmalladditionsofHClaqueoussolutions.Supersaturatedsolutionswereobtainedbymixingprotein-waterstocksolutionswithconcentratedsalt-waterstocksolutions.Astocksolutionofsodiumacetatebuffercontaining0.02%sodiumazidewasalsoaddedtominimizepossiblepHchanges.Thesupersaturatedsolutionswerethenleftat4°C(3-15daysdependingonsamplesupersaturation)toproducecrystals.Amicroscopicobservationofthecrystalmorphologywascarriedouttakingadvantageofelectronicmicroscopy.DigitalimagesanalysiswasperformedusingaPanasonicLumixDMC-FZ20systemoperatingthePanasonicimageanalysisconnectedtoa3CCDcolorvisioncameramountedonanOlympusBH2opticalmicroscope.Cerius2modelingCrystalstructuresdataoforthorhombiclysozymewasachievedbyProteinDataBanktocalculatetheoreticalmorphology.Experimentalmorphologywasusedtocomparewiththeresultofthesimulatedtheoreticalmorphology.CrystalMorphologyPredictionUsingBFDHandAEmodelClassicmolecularmechanicsanddynamicssimulationswerecarriedoutusingtheCerius2softwarepackage.Thecrystalbuilder,molecularmechanics,andmorphologypredictionmoduleswereemployedtoaccomplishthemodelinggoalsofthisresearch.Universalforcefieldwasused.RESULTSANDDISCUSSIONStructureanalysisThecrystalstructureoflysozymereportedincludedorthorhombic,monoclinic,tetragonalandtriclinicforms,whichhasbeenresolvedusingX-raydiffractionstudies.Inoursimulations,thecrystalstructureoftheorthorhombiclysozymewasobtainedfromtheProteinDataBank14.Thestere-structureoflysozymemoleculeisshowninFigure1.Theunitcellofitscrystalstructureareasfollows:a=30.47,b=59.39,c=68.78,=90°,andspacegroupisP212121,whichisshowninFigure2.Figure1.Thestere-structureoflysozymemolecule:ballandstickstyle(left),stickstyle(right)Figure2.TheunitcelloforthorhombiclysozymeMorphologypredictionbyBFDHmodelStructureanalysisThethreefaces110,002,011,101dominatedthecrystalhabitasrevealedfromtheBFDHmorphologicalprediction.TheresultsaresummarizedinTable.AccordingtoBFDHlaw,thelargertheinterplanardistancedhklis,thelargerthemorphologicalimportance(MI)ofthecorrespondinghklfaceis.FromtheTable1,itcanbeseenthatdhklof011faceisthelargestanditsfacetareais63.49%oftotalfacetarea,whichmeans011faceownsthelargestthemorphologicalimportance(MI).Thedhklof101faceisslightlylargerthan110faceanddhklof002isthesmallest.CrystalmorphologypredictedbyBFDHmodelcanbededucedandwasshowninFigure3(a),revealingallfaceswhichbelongtothecrystallographiczones110,002,011,and101,respectively.FromFigure3(a),thereisawellagreementwithcalculationresultinTable1.The011facewasfoundtobethelargestinarea,subsequently101face,110faceand002face.TABLE.CALCUTIONRESULTSOFBFDHMODELhklMultiplicitydhklTotalfacetareaTotalfacetarea/%011444.9590.9763.49002234.399.0406.309101427.8626.8718.75110427.1116.4111.45MorphologypredictionbyAEmodelTheresultsofattachmentenergycalculationswerereportedinTable.Basedonattachmentenergymodel,theMIsequenceisbasedontheassumptionthatthegrowthratesoffacesareproportionaltoattachmentenergy.Hence,thegreaterattachmentenergy,thefasterthecorrespondingfacegrowsandthesmalleritsmorphologicalimportance.Oncetherelativegrowthratesofthesignificantfacesareknown,amacroscopicimageofthecrystalmorphologycanbepostulated.CrystalmorphologypredictedbyAEmodelcanbededucedandwasshowninFigure3(b).Themorphology(Figure3a)predictedbyBFDHmodelisapproximatelyconsistentwiththatpredictedAEmodel.Themaindifferencesexistindisappearanceof002faceinAEmodel.FromTableandFigure3(b),itcanbeseenthatthereisagoodcorrelationbetweengrowthfaceareaandtheattachmentenergyofthedominantface.Themorphologypredictedisplatelike,suggestingthatasinglefaceisindeeddominatingtheoverallcrystalgrowth.Inthesecases,themostdominantfacehasanattachmentenergywhoseabsolutevalueismuchlowerthanthatfortheotherfaces,andtherelativegrowthrateisnotablysmaller.Hence,onefacedominatesthecrystalgrowth.TABLE.CALCULATIONRESULTSOFAEMODELhklMultiplicitydhklEatt(total)TotalfacetareaTotalfacetarea%/Å/(kcal/mol)/Å2/Å2011444.95-74.1116130076.80101427.8-148.74285020.40110427.11-167.758782.799(a)(b)Figure3.Predictedcrystalmorphology(a)BFDHmodel,(b)AEmodelComparisonofpredictedmorphologywithexperimentalgrowthcrystalThemorphologicalsimulationswereconfrontedwithexperimentaldataobtainedfromsolutiongrowncrystals.Thecrystalmorphologywasassessedbyopticalmicroscopy.MorphologicalsketchwasshowninFigure4.Figure4.ExperimentalcrystalmorphologyThepredictedmorphologybasedupontheAEmodelandBFDHmodelcomparesfavorablywiththeexperimentalmorphology.Themaindifferencebetweentheexperimentalmorphologyandpredictedmorphologyisthat002facedisappearsintheexperimentalmorphology,whichismoreconsistentwiththemorphologypredictedbyAEmodelthanbyBFDHmodel.Despitedifferencebetweenpredictedandexperimentalmorphology,itshouldberememberedthat,overall,thepredictedandsolutiongrownmorphologyarerathercomparable.Byfurtheranalyzingmorphologypredictedbytwomodelsandexperimentalmorphology,itcanbeseenthat011faceisthemostdominantinarea.Thiswasperhapsduetothefactthatthisfacehadsmallergrowthpromotinghydrogen-bondingcomponentinvolvedintheintermolecularinteractionsinvolvedinitsattachmentenergy,incontrasttothoseonthe110,002,101faces.SurfacechemistryvisualizationfortheimportantmorphologicalformsThesurfacechemistryofthethreemainsurfaceswasinvestigatedtoprovidebetterunderstandingtointeractionandrecognitionamonglysozyme:thethreefaces011,101,and110expecteddominatethecrystalhabitasrevealedfromtheAEmorphologicalprediction.(a)(011)(b)(101)(c)(110)Figure5.Cleavagestructureofmainexposurecrystalplane(a)011,(b)101,(c)110TheresultsaresummarizedinFigure5.The011surface(Figure5a)wasfoundtoberoughonthemolecularlevel.Incontrastto101crystalface,chemicalgrouprevealingisbigger,whichcanbringaboutagreatersterichindrancetoholdbackpackingoflysozymemoleculeinsolutionin011face.Asaresult,thegrowthrateof011crystalfaceiscomparativeslow,showinghighmorphologicalimportance.Examinationofthe101surface(Figure5b)revealsasmoothersurfacecomparedwith011faceonthemolecularlevel.Thechemicalgrouprevealingissmall,whichallowslysozymemoleculeeasilytopackin101surface,leadingto101surfacehaverapidgrowthrate.The110surface(Figure5c)wasfoundtobeveryopenandroughonthemolecularlevelrevealingadiagonalpatterndownthroughthesurfacewithalternatingorientationsofthelysozymemolecule.Thepatternproducedwasdiagonalinnaturewithhydrogenbondinganddonororacceptoratomswereactiveonthesurfaceforbindingoncomingmolecules.Inaddition,withthissurface,bothamino-groupandhydroxylcomponentsoflysozymemoleculewerefoundtoberevealabletopossiblebindingwithappropriatemolecules.Asaresult,thegrowthrateof110surfaceisthegreatest,showingthesmallestmorphologicalimportance.CONCLUSIONSThecrystalmorphologyoforthorhombiclysozymeispredictedusingAEmodelinconjunctionwithBFDHmodel.Accuratemodelingandmorphologypredictionofthiscrystalareachieved.Themorphologypredictedbytwomodelsisapproximatelyconsistent.ThemorphologypredictedbyAEmodelismoreconsistentwiththemorphologyofcrystalsgrownfromsolutionthanBFDHmodel.Themaincrystalfaces011,101and110areobservedanddominantinmorphologypredictedbyAEmodel.Bycleavingthesedominantcrystalfaces,surfacechemistryvisualizationandtheoreticalanalysisbasedoninteractioninintra-moleculesorinter-moleculesfortheimportantmorphologicalformsareperformed.TheresultindicatesthatsterichindranceeffectandH-bandinteractionplayscriticalrolefortheplate-likemorphologyoforthorhombiclysozyme,whichprovidesaimportantinstructioninmolecularlevelforpreparationofhigh-qualitylysozymecrystal.ACKNOWLEDGMENTThisworkwasfundedbytheNationalNaturalScienceFoundationofChina(No.20806053),DoctoralFundofMinistryofEducationofChina(No.200800561029)andChinaPostdoctoralScienceFoundation(No.20080440677)REFERENCES1J.C.Givand,R.W.Rousseau,P.J.Ludovice,“CharacterizationofL-isoleucinecrystalmorphologyfrommolecularmodeling,”J.Cryst.Growth,vol.194,pp.228238,1998.2J.Prywer,“ExplanationofsomepeculiaritiesofcrystalmorphologydeducedfromtheBFDHlaw,”J.Cryst.Growth,vol.270,pp.699710,2004.3G.Pfefer,R.Boistelle,“Theoreticalmorphologyofadipicacidcrystals,”J.Cryst.Growth,vol.208,pp.615622,2000.4P.Hartman,“Theattachmentenergyasahabitcontrollingfactor:I.Theoreticalconsiderations,”J.Cryst.Growth,vol.49,pp.145156,1980.5S.David,C.Coombes,A.Richard,“CalculationofAttachmentEnergiesandRelativeVolumeGrowthRatesasanAidtoPolymorphPrediction,”Cryst.GrowthDes,vol.5,pp.879885,2005.6C.M.Li,K.L.Kirkwood,G.D.Brayer,“TheBiologicalCrystallizationResource:FacilitatingKnowledge-BasedProteinCrystallizations,”Cryst.GrowthDes,vol.7,pp.21472152,2007.7S.X.Lin,A.McPherson,R.Giegé,“GoodCrystals,StillaChallengeforStructuralBiology,”J.Cryst.Growth,vol.7,pp.21242125,2007.8A.Warke,C.Momany,“AddressingtheProteinCrystallizationBottleneckByCocrystallization,”J.Cryst.Growth,vol.7,pp.22192225,2007.9M.W.Elaine,J.R.Kevin,J.M.Stephen,“AMolecularDynamicsStudyofSolventandImpurityInteractionontheCrystalHabitSurfacesof-Caprolactam,”Langmuir,vol.14,pp.56205630,1998.10N.Kubota,J.W.Mullin,”Akineticmodelforcrystalgrowthfromaqueoussolutioninthepresenceofimpurity,”J.Cryst.Growth,vol.152,pp.203208,1995.11H.C.Koolman,R.W.Rousseau,“EffectsofisomorphiccompoundsonthepurityandmorphologyofL-isoleucinecrystals,”AIChEJournal,vol.42,pp.147153,1996.12J.Z.Chen,N.F.Zhuang,S.K.Lin,“Theoreticalmorphologyandgrowthhabitofrubidiumhydrogenselenatecrystals,”J.Cryst.Growth,vol.205,pp.584589,1999.13Y.Matsuzukia,T.Kubotab,X.Y.Liua,“AFMobservationofthesurfacemorphologyandimpurityeffectsonorthorhombichenegg-whitelysozymecrystals,”J.Cryst.Growth,vol.242,pp.199208,2002.14http:/www.rcsb.org/pdb/home/home.do