Molecular Dynamics Simulations Of Biological Macromolecu

计算生物物理学部分使用几种理论技术研究生物学意义的问题：分子动力学，分子力学，建模，小分子结构的从头计算分析和分子图形。这些技术适用于各种各样的大分子系统。应用于生物医学感兴趣的分子的特定项目使用分子动力学模拟来预测肽和蛋白质的功能或结构。这些项目包括： - 天然和突变型vnd/NK-2同源结构域-DNA复合物的分子动力学 - C3 a过敏毒素和抗体结合位点 - 人鼻病毒的蛋白质结构稳定性和活性 - 基于QM/MM方法的腺苷激酶催化机理模拟 - 用QM/MM方法研究D5-3-酮类异构酶的催化机理 - 用QM/MM方法研究N-乙酰基转移酶的催化机理 - 用QM/MM方法研究Chorismate Mutase的催化机理 - β-内酰胺水解催化途径的追踪 - 确定人类基因NFAT 5的编码区和启动子区。 - 探索确定一组蛋白质序列的相互关联性的方法。 - 辛基葡萄糖苷胶束的模拟 - 将辛基葡糖苷模拟扩展到肽Mastoparan X，一种黄蜂毒液 - 确定果蝇中13，500个基因中的哪一个编码信号肽 基础研究正在进行中，以更好地了解大分子系统。这些项目包括研究： - 伴侣蛋白介导的蛋白质折叠 - NMR屏蔽张量计算 - 脂质双层凝胶相模拟 - 研究核酸结构的环境依赖性 - 亮氨酸拉链的建模：卷曲螺旋平行与反平行取向的起源 - 肌红蛋白和溶菌酶晶体和溶液的模拟 - CI 2分子动力学模拟 利用计算机模拟研究了伴侣蛋白分子介导的蛋白质折叠。我们的重点是大肠杆菌的GroEL-GroES伴侣蛋白复合物，其相关结构是已知的。使用高性能的科学计算方法，例如结合最先进的CHARMM模拟程序对蛋白质进行粗粒度和全原子描述。这些研究将有助于阐明分子伴侣对蛋白质结构的影响，分子伴侣系统的机制，以及分子伴侣循环的时间尺度。 对1H和15 N核的NMR屏蔽张量进行了从头计算，并与实验结果进行了比较。酰胺1H的计算值和实验值之间几乎没有相关性，并研究了基组和同位素效应等因素的作用作为差异的潜在原因。目前正在进行一项15 N研究，以1H结果为起点。 在生理条件和连续介质模型下用分子动力学方法计算酵母转录因子GCN 4的亮氨酸拉链结构域（GCN 4-p1）的自由能。亮氨酸拉链基序是由两个α-螺旋组成的平行左手超螺旋。据估计，本机（平行）的对齐是积极更稳定的非本地反平行对齐，其中静电能有助于显着的整体充满活力的图片的两个方向，以及平行与反平行方向的偏好。 腺苷被ATP磷酸化是由腺苷激酶催化的。目前正在利用量子力学/分子力学技术研究该机制的不同途径。用我们的带高斯模糊的双键原子方法计算了腺苷5'醇的酸度和天冬氨酸的质子亲合能。这些都是整个机制中的关键步骤，并提供了信心，我们的QM/MM方法将提供合理的答案，同时研究整个系统。 采用全原子模型和显式水相结合的SGMD方法研究蛋白质折叠机理。我们已经成功地观察到了两态蛋白质折叠现象。

The Computational Biophysics Section studies problems of biological significance using several theoretical techniques: molecular dynamics, molecular mechanics, modeling, ab initio analysis of small molecule structure, and molecular graphics. These techniques are applied to a wide variety of macromolecular systems. Specific projects applied to molecules of biomedical interest uses molecular dynamics simulations to predict function or structures of peptides and proteins. Such projects include: - Molecular dynamics of native and mutant vnd/NK-2 homeodomain--DNA complexes - C3a anaphylatoxin and antibody binding sites - Protein structure stabilization and activity in human rhinovirus - Modeling the catalytic mechanism of adenosine kinase with QM/MM methods - The study of the catalytic mechanism of D5-3-Ketosteroid Isomerase using QM/MM methods - The study of the catalytic mechanism of N-acetyltransferase using QM/MM methods - The study of the catalytic mechanism of Chorismate Mutase using QM/MM methods - Tracing the catalytic pathway of b-lactam hydrolysis - Determine the coding and promoter regions for the human gene NFAT5. - exploring methods for determining the inter-relatedness of a group of protein sequences. - Simulations of octyl glucoside micelles - Extend octyl glucoside simulations to the peptide Mastoparan X, a wasp venom - Determining which of the 13,500 genes in the Drosophila code for signal peptides Basic research is underway to provide a better understanding of macromolecular systems. The projects include studies of: - Chaperonin-mediated protein folding - NMR Shielding Tensor calculations - Lipid bilayer gel phase simulations - Investigating the environmental dependence of nucleic acid structure - Modeling leucine zippers: Origins of parallel vs. antiparallel orientation of coiledcoils - Simulations of myoglobin and lysozyme crystals and solutions - Molecular dynamics simulations of CI2 Protein folding mediated by chaperonin molecules is studied using computer simulations. Our focus is on the GroEL-GroES chaperonin complex of the Escherichia coli, for which the associated structures are known. High-performance scientific computing methods are used, such as coarse-grained and all-atom descriptions of proteins in conjunction with the state-of-the-art CHARMM simulation program. These studies will shed light on the effect of chaperones on the protein structure, the mechanism of the chaperonin system, and the timescales in the chaperonin cycle. Ab initio calculations of NMR shielding tensors were performed for comparison with experimental studies for 1H and 15N nuclei. There was little correlation between calculated and experimental values for amide 1H, and role of factors such as basis set, and isotope effect were investigated as potential causes for the discrepancy. An 15N study is currently underway, using the 1H results as a starting point. Free energy calculations on the leucine zipper domain (GCN4-p1) of the yeast transcription factor GCN4 using Molecular Dynamics (MD) under physiological conditions and continuum models. The leucine zipper motif is a parallel left-handed supercoil composed of two a-helices. It is estimated that the native (parallel) alignment is energetically more stable than the non-native antiparallel alignment where electrostatic energies contribute significantly in the overall energetic picture of both orientation as well as to the preference of the parallel vs the antiparallel orientation. The phosphorylation of adenosine by ATP is catalyzed by Adenosine Kinase. The different pathways in which the mechanism can proceed are being worked out by using quantum mechanical/molecular mechanical techniques. Using our double link atom method with gaussian blur, we have calculated the acidity of the 5' alcohol of adenosine and the proton affinity of aspartate. These are crucial steps in the overall mechanism and provides confidence that our QM/MM method will provide reasonable answers while studying the entire system. Studying protein folding mechanism through SGMD simulation with all-atom model and explicit water. We have successfully observed two-state protein folding phenomena.