Molecular dynamics simulation of protein crystal with polarized protein-specific force field

来源 :第十二届全国量子化学会议 | 被引量 : 0次 | 上传用户:wenlingqiang6268047
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
  The explicitly polarized force field model is indispensable in the simulation of protein crystal due to its particular electrostatic environment which is different from that in water solution.The polarized protein-specific charge(PPC)has shown important priorities for the study of protein-ligand binding and protein folding.However,the applicability of PPC in the simulation of protein crystal is yet to be examined.Two 250 ns molecular simulations were carried out to study the structure and dynamics of crystal toxin protein Ⅱ from the scorpion Androctonus australis Hector employing PPC as well as the standard AMBER99SB force field to investigate the effect of electrostatic polarization on the crystal stability.Our simulation results showed that under PPC,the monomer in the subunit cell and the lattice in the supercell were more stable with much smaller RMSDs and the lattice atomic fluctuations were also better in line with the crystallographic B-factors.Most of the interactions at interfaces in the X-ray structure were well preserved,underscoring the important effect of polarization on maintaining the crystal stability.However,our results also showed that the hydrogen bond between Asp53 and Gln37 and the cation-л interaction between Arg56 and His64 were not as stable,indicating that further optimization of force field with van der Waals interaction parameters is desired.
其他文献
The hydrogen bond and halogen bond have attracted more and more attention.In this work,12 halogen bonded and 3 hydrogen bonded models are adopted to scan the potential energy surface in the second-ord
会议
运用二级微扰理论探究气相条件下羟基自由基和反式2-氯乙烯基二氯胂(C2H2As Cl3)的反应机理.在MP2/6-311++G(d.p)水平下优化反应物,中间体,过渡态及产物的几何构型,并在同一理论水平下计算了各个驻点的振动频率,零点能.为得到更可靠的势能面信息,基于MP2/6-311++G(d.p)水平下优化的几何构型,采用G3和CCSD(T)的方法进行了单点能的计算.在本研究体系中,所有的能量
会议
The adsorption of UO2(H2O)52+ on the hydroxylated α-SiO2(001)surface was studied by a periodic density functional theory(DFT)and ab initio MD(AIMD)simulations.The effects of pH,CO2 and the aqua soluti
会议
A hydrophobic task-specific ionic liquid designed to extract Cu2+ from water was prepared by appending thiourea-substituted alkyl groups to imidazoles and combining the resulting cationic species with
会议
Based on our previous work,[1,2] the molecular dynamics simulations and first-principles quantum mechanical/molecular mechanical free energy calculations(QM/MM-FE)have been performed to provide the de
会议
Application of the electronic circular dichroism(ECD)has been proved to be practical and reliable,especially when the molecule has high conformational flexibility.This method is superior to optical ro
会议
The catalytic cycle of α-ketoglutarate(α-KG)dependent non-heme Fe(II)halogenase SyrB2 was studied by hybrid QM/MM approach.Our results show that with the native substrate(L-Thr),only the quintet spin
会议
Tetrabromobisphenol A(TBBPA),a widely used brominated flame retardant,is harmful to human and wildlife in various aspects as an endocrine disruptor.The DFT calculations both in gaseous phase and aqueo
会议
本文针对在生命和材料科学中有广泛应用、具有典型的pi-pi 相互作用的苯分子,构建并对比了两种适用于苯分子的粗粒化模型:不带电荷的三点模型和带电荷的四点模型。粗粒化粒子的质量通过转动惯量确定,电荷分布通过分子四极距确定;分子在模拟中用刚性结构,无分子内作用能;分子间作用力用LJ(12-6)和库仑势函数表达,而参数通过拟合径向分布函数、气液相平衡共存曲线和晶体的密度和升华焓确定。三点模型可以准确的描
会议
In order to study the influence of molecular structure on the formation of a monolayer,two molecules have been considered,namely N-stearoyldopamine(DOPA)and 4-stearyl-catechol(ST).The difference betwe
会议