Developing Improved Methods for Modeling and Simulating Protein Structures

开发蛋白质结构建模和模拟的改进方法

• 批准号: 7733457
• 负责人: HOMER ROBERT GUY
• 金额: $ 16.27万
• 依托单位: DIVISION OF BASIC SCIENCES - NCI
• 依托单位国家: 美国
• 资助国家: 美国
• 起止时间: 至
• 项目状态: 未结题
• 来源: https://reporter.nih.gov/project-details/7733457
• 关键词: Algorithms; Alzheimer' s Disease; Amyloid beta-Protein; Arts; Boston; Computer Systems; Computing Methodologies; Drug Binding Site; Drug Design; Goals; Homology Modeling; Ion Channel; Ions; Lead; Lipid Bilayers; Mathematics; Membrane Proteins; Methods; Modeling; Molecular; Molecular Conformation; Names; Numbers; Peptides; Pharmacology; Potassium Channel; Process; Property; Proteins; Purpose; Side; Simulate; Structural Models; Structure; Testing; Time; Universities; Water; Work; base; improved; method development; models and simulation; molecular dynamics; molecular modeling; nanosecond; neurotoxicity; programs; protein structure; restraint; simulation; success

Our first computational methods subproject involves using symmetry constraints in molecular modeling and simulations. For those cases in which proteins or protein assemblies are composed of identical subunits, the subunits usually have identical conformations and identical interactions with neighboring subunits. We have incorporated this constraint into our models of ion channels for over two decades, with considerable success. For example our four-fold symmetric model of the ion selective region of potassium channels was virtually identical to those of subsequently determined potassium channel crystal structures. However, until recently, we have not incorporated symmetry constraints into an automated algorithm, nor demonstrated that doing so actually improved the quality of structural models. One of our collaborators, Andriy Anishkin, has developed a program that uses harmonic restraints to restore symmetry during molecular dynamic simulations. We have tested the utility of this program for homology modeling of ion channels by analyzing crystal structures of three distantly related ion channels named KcsA, KirBac, NaK. The first step of this project was to develop two homology models of each channel using the other two channel structures as templates. Next, molecular dynamics simulations of these six models, along with the three crystal structures, were performed with the proteins embedded in an explicit lipid bilayer with water and ions on each side and within the pore of the channels. After eight nanoseconds of unrestrained simulations, symmetry restraints were imposed. This process was repeated two more times for each model and crystal structure. The models developed without molecular dynamic simulations, with unrestrained simulations, and with symmetry-restrained simulations were then compared to the original structures to determine which was better. The major finding was that although unrestrained molecular dynamic simulations did not improve the homology models, imposition of symmetry restraints did lead to substantial improvement. The greatest improvement occurred for the pore lining segments, where interaction among adjacent subunits is extensive. Our principal purpose for developing homology models of ion channels is to understand their structure and pharmacology well enough to utilize the models in structure-based drug design. Thus, it is noteworthy that the symmetry restraints substantially improve models of the pore region that forms the principal drug binding sites. Our second computational methods subproject involves simplifying molecular representations of proteins so that properties of proteins and peptides can be simulated for much longer times. A major limitation of conventional molecular dynamic simulations is that they can be performed for only short periods of time, typically less than a tenth of a microsecond. This is substantially shorter than the time required for most conformational changes, and many orders of magnitude shorter than the time required for assembling amyloid beta structures. To overcome this limit, Dr. Sijung Yun of our group has been using discrete molecular dynamics that is about seven orders of magnitude more efficient than conventional molecular dynamics. He helped develop this program during his doctoral work at the University of Boston. Three simplifications are introduced in discrete molecular dynamics in order to increase efficiency while still preserving accuracy. First, interaction potentials are simplified into discrete steps. This greatly simplifies the mathematics and number of calculations required during the simulations. Second, interactions due to water molecules are replaced by effective potentials. This greatly reduces the number of atoms simulated by eliminating water molecules. Third, we used four-bead representation of a residue in a protein instead of representing all the atoms of the residue. This also reduces the number of simulated atoms. While here, Dr. Yun has worked to improve the parameters and test how well known protein crystal structures are maintained when subjected to these simulations. He has also been using this approach to test models of amyloid-beta hexamers developed by our group. Amyloid-beta hexamers are involved in neurotoxicity of Alzheimers disease. Our third computational methods subproject involves developing models of protein structures on regular lattices. We are using lattices that have the types of symmetry that occur in our models of amyloid beta assemblies. The lattice approach should allow us to examine and evaluate all possible ways of constructing peptide assemblies that have this specific type symmetry. Such exhaustive analysis is not computationally feasible even with simplified discrete molecular dynamics. This project has just been started and we are still evaluating the feasibility of the approach.

我们的第一个计算方法子项目涉及在分子建模和模拟中使用对称约束。对于由相同亚基组成的蛋白质或蛋白质组合，这些亚基通常具有相同的构象，并且与邻近亚基具有相同的相互作用。二十多年来，我们已经将这一限制纳入到离子通道模型中，并取得了相当大的成功。例如，我们的钾通道离子选择区域的四重对称模型几乎与随后确定的钾通道晶体结构相同。然而，直到最近，我们还没有将对称约束纳入自动化算法，也没有证明这样做实际上提高了结构模型的质量。我们的一个合作者，Andriy Anishkin，开发了一个程序，使用谐波约束来恢复分子动力学模拟中的对称性。我们通过分析KcsA， KirBac， NaK三个远亲离子通道的晶体结构，测试了该程序对离子通道同源性建模的实用性。该项目的第一步是使用其他两个通道结构作为模板，开发每个通道的两个同源模型。接下来，对这六种模型以及三种晶体结构进行分子动力学模拟，将蛋白质嵌入显性脂质双分子层中，每侧和通道孔内都有水和离子。在8纳秒的无约束模拟后，对称约束被施加。这个过程对每个模型和晶体结构重复两次以上。然后将没有分子动力学模拟、无约束模拟和对称约束模拟的模型与原始结构进行比较，以确定哪种模型更好。主要的发现是，尽管不受约束的分子动力学模拟并没有改进同源性模型，但施加对称约束确实导致了实质性的改进。最大的改善发生在孔隙衬里段，其中相邻亚基之间的相互作用是广泛的。我们开发离子通道同源模型的主要目的是了解它们的结构和药理学，以便在基于结构的药物设计中利用这些模型。因此，值得注意的是，对称约束大大改善了形成主要药物结合位点的孔区域模型。我们的第二个计算方法子项目涉及简化蛋白质的分子表示，以便蛋白质和肽的性质可以模拟更长的时间。传统分子动力学模拟的一个主要限制是，它们只能在很短的时间内进行，通常不到十分之一微秒。这比大多数构象变化所需的时间要短得多，比组装β淀粉样蛋白结构所需的时间要短许多个数量级。为了克服这一限制，我们小组的Yun Sijung博士一直在使用比传统分子动力学效率高7个数量级的离散分子动力学。他在波士顿大学攻读博士学位期间帮助开发了这个项目。为了在保持精度的同时提高效率，在离散分子动力学中引入了三个简化。首先，将相互作用势简化为离散步骤。这大大简化了模拟过程中所需的数学和计算数量。其次，由于水分子的相互作用被有效势所取代。通过消除水分子，这大大减少了模拟的原子数量。第三，我们使用四头表示法来表示蛋白质中的残基，而不是表示残基的所有原子。这也减少了模拟原子的数量。在这里，Yun博士一直致力于改进参数，并测试在这些模拟中如何保持已知的蛋白质晶体结构。他也一直在使用这种方法来测试我们小组开发的淀粉样蛋白- β六聚体模型。淀粉样蛋白-六聚体参与阿尔茨海默病的神经毒性。我们的第三个计算方法子项目涉及在规则晶格上开发蛋白质结构模型。我们使用的晶格具有在淀粉样蛋白组装模型中出现的对称类型。晶格方法应该允许我们检查和评估所有可能的方法来构建具有这种特定类型对称性的肽集。这种详尽的分析在计算上是不可行的，即使是简化的离散分子动力学。这个项目刚刚开始，我们还在评估这个方法的可行性。