分子力学方法介绍PPT演示课件

1,分子力学简介IntroductiontoMolecularMechanics,2,分子力学是模拟分子行为的一种计算方法。分子力学认为分子体系的势能函数时分子体系中原子位置的函数。分子力学将分子体系作为在势能面上运动的力学体系来处理，求解的是经典力学方程，而不是量子力学的薛定谔方程。所以分子力学方法可以求得分子的平衡结构和热力学性质等，但不能得到分子体系与电子结构有关的其它性质。分子力学的势能函数表达方程很简单，其计算速度很快（约是半经验量子化学计算方法速度的1000倍），能够用于生物大分子体系的计算。对于力场参数成熟的分子力学方法，已经可以达到很高的计算精度。,3,目前，分子力学是模拟蛋白质、核酸等生物大分子结构和性质以及配体和受体相互作用的常用方法，在生命科学领域得到了广泛的应用。随着分子图形学的不断发展，分子力学已经广泛应用于分子模型设计。当今优秀的分子设计程序都将分子力学作为初始模型优化的主要方法，分子力学和分子图形学已经充分地糅合在分子设计中，分子模型的构建也是分子力学为主，分子力学方法是计算机辅助分子设计中常用的方法。特别是在药物设计中，已经离不开分子力学计算和模拟方法。应用分子力学方法，可以迅速得到分子的低能构象，通过构象分析可以获得合理的药效构象和药效基团。如果已知受体的三维结构，可以用分子力学模拟药物与受体的相互作用。在分子的定量结构活性关系研究中，也需要用分子力学方法进行计算。,4,分子力学主要应用分子力场方法计算分子的势能，分子力学使用解析经验势能函数来描述分子的势能面。每个分子都有固定的立场。在分子力学实际计算时，将力场分解成不同的组分，使用理论计算和实验拟合的方法建立力场参数。力场参数要有可移植性，它要适应于同类分子的计算，即同一类分子也要有很高的计算精度。由于分子力学是经验的计算方法，不同的分子力学方法采用不同的势能函数表达式，而且力场参数值也不相同。通常将分子的势函数分解成键伸缩能Estretch、角弯曲能Ebend、二面角扭曲能Etorsion、范德华作用能EVdW和静电作用能Eelectrostatic等项，总能量可以表示为：StericEnergy=Estretch+Ebend+Etorsion+EVdW+Eelectrostatic,分子力学的基本原理,5,其中与化学键有关的有：键伸缩能Estretch、角弯曲能Ebend、二面角扭曲能Etorsion，与化学键无关的有范德华作用能EVdW和氢键以及静电作用能Eelectrostatic。对于更精确的力场，可以在势函数表达式中增加一些交叉项，如键伸缩-键弯曲振动交叉项、二面角扭曲-键伸缩交叉项等，也可以加入氢键函数项等其他修正项。,6,基本概念,7,StericEnergy,Thedifferencebetweentheactualenergyofthesystemandtheenergyofahypotheticalsystemwhereallstructuralparametersareattheirequilibrium(minimumenergy)values.,StericEnergy=Estretch+Ebend+Etorsion+EVdW+Eelectrostatic+Estretch-bend+Etorsion-stretch+,8,AGeneralForceFieldCalculation,9,AtomTypes,Inmolecularmechanicsatomsaregiventypes-thereareoftenseveraltypesforeachelement.Typesdependontheenvironment:,“Transferability”isassumed-forexamplethataC=Obondwillbehavemoreorlessthesameinallmolecules.,10,ASimpleMolecularMechanicsForceField,11,GeneralFeatures,12,键伸缩能,根据经典力学的胡克定律，两个成键原子可以视为有一根弹簧连接的两个小球，其势函数是简谐振动势函数，其计算非常方便，非常实用于生物大分子的计算。,kb键伸缩常数，r0平衡键长,StretchingEnergy,13,StretchingEnergy,14,BondStretching,Harmonicpotential(AMBER),15,BondStretching,Cubic(MM2)andquadratic(MM3)potentials,16,BondStretchingParameters(MM2),17,角弯曲能BendingEnergy,两个相连的化学键（三个成键原子）具有一定的键角，键角弯曲势能函数与键伸缩势函数类似，可用虎克定律描述。键弯曲能常数，实际键角和平均键角。,18,19,AngleBending,AMBER:,MM2:,MM3:,20,AngleBendingParameters(MM2),21,TorsionEnergy（二面角扭曲转动能）,三个相连的化学键（四个成键原子）具有一定的二面角（扭角）。头尾两个成键原子绕中间化学键旋转时，需要一定的能量，存在明显的旋转势垒，在分子力学中，二面角扭曲转动势函数常用傅立叶级数模拟。,22,23,FunctionalForm,Vn=4,n=2,g=180,Vn=2,n=3,g=0,24,FunctionalForm:AMBER,25,FunctionalForm:MM2,MM3,26,Example:Butane,TypeV1V2V3C-C-CC0.2000.2700.093C-C-CH0.0000.0000.267H-C-CH0.0000.0000.237,27,CrossTerms:Stretch-Bend,Couplingbetweeninternalcoordinates.Importantforreproducingstructuresofunusual(e.g.,highlystrained)systemsandvibrationalspectra.,Stretch-Bend:Asabondangleisdecreased,theadjacentbondsstretchtoreducetheinteractionbetweenthe1,3atoms.,28,CrossTerms:Stretch-Torsion,Stretch-Torsion:ForanA-B-C-Dtorsion,thecentralB-CbondelongatesinresponsetoeclipsingoftheA-BandC-Dterminalbonds.,29,CrossTerms,30,Non-BondedEnergy非键相互作用能,非键相互作用能主要包括范德华相互作用能和静电相互作用能。,31,VanDerWaals(VdW)Interactions,范德华相互作用能,范德华相互作用能是分子力学中重要的作用，它可分为进程作用的排斥力和远程的吸引力，一般以Lennard-Jones公式来模拟范德华非键相互作用分子力场。其中r表示原子之间的距离，A和B是非键原子对的范德华常数。R的六次方为原子对的吸引作用，R的12次方位排斥作用。,32,HydrogenBonding,H-bondgeometryindependent:,risthedistancebetweentheH-bonddonorandH-bondacceptor.,H-bondgeometrydependent:,33,PointChargeElectrostaticModels,qi,qjarepointcharges.Charge-chargeinteractionsarelongranged(decayasr-1).Whenenoughpointchargesareusedalltheelectricmomentscanbereproduced.Whenqi,qjarecenteronthenucleitheyarereferredtoaspartialatomiccharges.,PartialatomicchargesFittothermodynamicproperties.AbinitioCalculations(6-31G*,MNDO).Mullikencharge.Reflectmolecularconstitutionratherthanmolecularinteractions.Electrostaticpotential.,34,RapidMethodsforCalculatingAtomicCharges,35,ThermodynamicProperties,36,分子力学总能量并没有严格的物理意义，为了获得精确的计算结果，分子力场的能量函数表达式的复杂程度没有限制。根据适用范围不同，目前已发展了很多的分子力场，分子力学方法根据适用范围可分为小分子和大分子两类。适用于小分子的分子力学的函数形式较复杂和精确。例如MM2/MM3/MM4TINKERUFFMOMECCOSMOS等。适用于大分子的分子力学函数形式较简单，例如：AMBERCHARMMGROMOSOPLSCVFFCFFMMFF等。,37,分子力场的参数化,分子力学计算结果的精确性除了与力场势函数表达式有关外，还与力场参数的数值密切相关。有效的力场势函数和正确的力场函数可使分子力学计算达到很高的精度，有些分子力学方法对分子的性质计算结果非常好，可以达到实验值的精度。由于不同的分子力学方法所用的势函数表达方式不同，其力场参数不能相互代替。由于化合物的原子类型很多，目前分子力场不可能适用于各种的原子类型，满足各种化合物计算的要求。在缺乏力场参数时，就需要进行分子力场参数的优化。,38,分子力学势函数是由一系列的可调参数组成的。对可调参数进行优化，是分子力学的计算值最符合分子的某些性质的实验数值，得到一套力场的优化参数，再使用这套参数去预测相同原子类型的其它分子的结构和性质。传统的分子力场参数化方法是通过拟合实验数据（几何构型、构象能、生成热、光谱数据等）来优化参数。分子的几何结构和力场参数可由实验或量子化学计算获得。键伸缩振动常数可直接由价键力场计算或振动光谱获得。平衡键长、平衡键角和角弯曲常数可由X射线衍射、中子衍射、电子衍射等方法测定，也可由量子化学计算得到。扭转位能参数来自于NMR谱带和弛豫时间。构象能可从光谱和热化学数据得到。非键参数可从晶格参数和液体的物理性质数据得到。,分子力场的参数化,39,优化力场参数时，首先要选择某些代表性化合物，根据其原子和成键化学环境不断修正力场参数，使得力场参数的计算结果与真实分子的结构、能量和分子的振动光谱一致，接着进行同类型化合物分子的计算，验证这套参数的可靠性。力场参数的优化可使用最小二乘拟合法和尝试法。,分子力场的参数化,40,Parameterization,Choosingvaluesfortheparametersinthepotentialfunctionequationstobestreproduceexperimentaldata.,41,Parameters-Quantity,ForNatomtypesrequire:Nnon-bondedparametersN*NstretchparametersN*N*NbendparametersN*N*N*NtorsionparametersExampleMacroModelMM2*:39atomtypes,42,Parameters-Source,Aforcefieldparameterizedaccordingtodatafromonesource(e.g.,experimentalgasphase,experimentalsolidphase,abinitio)willfitonlyqualitativelydatafromothersources.,43,Parameters-Source,44,Parameters-Quality,45,TransferabilityofParameters,Constructgraphicalrepresentationofmoleculetobemodeled(“frontend”)Selectforcefieldmethodandterminationcondition(gradient,#cycles,ortime)PerformgeometryoptimizationExamineoutputgeometry.isitreasonable?Searchforglobalminimum.,StepsinPerformingMolecularMechanicsCalculations,46,Inputisusuallydonegraphically(bysketchingorbuildingstructuresatom-by-atomorbyassemblingcomponentparts).Thisgraphicalmodelisconvertedtoamathematicalmodelbythesoftware.Eachsoftwarepackagehasitsownfiletype,butmosthavesomecommonfeatures.PDBfileismostcommondenominator.,InputFileStructure,47,C1-1.1291.281-0.000C2-2.5581.772-0.000C3-3.5190.606-0.000H4-0.5961.6370.890H5-0.5961.637-0.890H6-2.7332.3920.890H7-2.7332.392-0.890H8-4.5580.9520.000H9-3.359-0.017-0.890H10-3.359-0.0170.890H11-1.1100.183-0.000(thisMAYbethesameasthe.PDBfile,asshownhere,ortheorientationofthemoleculemaybedifferent,makingthenumbersdifferent),Cartesiancoordinate(XYZ)file,48,distanceangledihedralconnectivityN0.000000.000000.00000000H1.020010.000000.00000100H1.02001104.536810.00000120H1.02001104.53681109.579611230(endoffile)(1meansoptimize,0meanskeepconstant,-1meansvaryaccordingtoadesignatedpattern),(sometimescalledZ-matrix),InternalCoordinates(forNH3),49,LocalminimumvsglobalminimumManylocalminima;onlyONEglobalminimumMethods:Newton-Raphson(blockdiagonal),steepestdescent,conjugategradient,others.,globalminimum,localminima,EnergyMinimization,50,“Stericenergy”hasNOphysicalmeaning,anditisdefineddifferentlyindifferentprogramsThereforeitCANNOTbeusedtocomparestructurescalculatedbydifferentprogramsItsuseislimitedtocomparingISOMERICstructureshavingtheSAMEnumberandkindsofbonds(conformers,stereoisomers).,Usesof“StericEnergy”,51,CalculationsareveryfastGeometryoptimizationsofsmalltomedium-sizemoleculescanbeaccomplishedonapcConformationsofmacromolecules(includingbiomacromoleculessuchaspeptidesandpolysaccharides)canbecalculatedusingworkstationsorparallelprocessingcomputers.,SuccessesofMolecularMechanicsCalculations,52,Reasonablegeometriesareusuallyobtained:Bondlengthswithin0.1AngstromofexperimentalvaluesBondangleswithin2ofexperimentalvalues.Calculatedenergiesareusuallyquitegood:Enthalpiesofformationwithin2kcal/mol(8kJ/mol)ofexperimentalvaluesProvidesinputstructureformoreinvolvedcalculations(molecularorbitalmethods).,SuccessesofMolecularMechanicsCalculations,53,Thecalculationsdonotaccountforelectrons!Orbitalinteractionsareignored!Theselectionof“atomtype”iscrucialtothecomputationalresult:e.g.,AMBERhas5typesofOxygen:carbonyl,alcohol,acid,ester/ether,water(seenextslide)NoconsiderationisgiventotheimportanceofdelocalizedpelectronsystemsOnlygroundstatesconsidered.notT.S.or*,LimitationsofMolecularMechanics,54,Obtainingareasonablygoodgeometry(instructureswherepielectronsarenotinvolved.Asastartingpointforfurthercalculations,suchassemi-empirical,abinitio,ordensityfunctional.Searchingtheenergysurfaceforminimumenergyconformations(itisusuallytooexpensivetodothisusingMOmethods).,UsesofMolecularMechanics,55,Optimizedgeometry(minimumenergyconformation)Equilibriumbondlengths,bondangles,anddihedral(torsional)anglesDipolemoment(vectorsumofbonddipoles)EnthalpyofFormation(insomeprograms).,PropertiesCalculated,56,Equalto“stericenergy”plussumofgroupenthalpyvalues(CH2,CH3,C=O,etc.),withafewcorrectiontermsNotcalculatedbyallmolecularmechanicsprograms(e.g.,HyperChemandTitan)Calculatedvaluesaregenerallyquiteclosetoexperimentalvaluesforcommonclassesoforganiccompounds.,EnthalpyofFormation,57,EnthalpyofFormation.,58,EnthalpyofFormation.,59,SybylMM+MM3ExptCH3CH3C-C1.5541.5321.5311.526C-H1.0951.1151.1131.109CH3COCH3C-C1.5181.5171.5161.522C-H1.1071.1141.1111.110C=O1.2231.2101.2111.222,BondLengths,60,SybylMM+MM3CH3CH3H-C-C109.5111.0111.4H-C-H109.4107.9107.5CH3COCH3C-C-C116.9116.6116.1H-C-H109.1108.3107.9C-C-O121.5121.7122.0,BondAngles,61,MM2/MM3(Allinger)best;generalpurposeMMX(Gilbert)addedTSs,otherelements;goodMM+(Ostlund)inHyperChem;general;goodOPLS(Jorgenson)proteinsa