COMPUTATIONAL STUDIES OF RNA RECOGNITION AND CATALYSIS

这个子项目是利用资源的许多研究子项目之一。 由NIH/NCRR资助的中心拨款提供。对子项目的主要支持 子项目的首席调查员可能是由其他来源提供的， 包括美国国立卫生研究院的其他来源。为子项目列出的总成本可能 表示该子项目使用的中心基础设施的估计数量， 不是由NCRR赠款提供给次级项目或次级项目工作人员的直接资金。 蛋白质、辅因子和其他RNA分子等其他生物分子对RNA的准确识别对许多细胞功能至关重要。利用各种计算化学工具，如分子动力学模拟、量子计算和混合量子力学/分子力学方法，我们的研究考察了三个主要领域，信使RNA-转移RNA(信使RNA-转移RNA)识别，蛋白质翻译的关键步骤，人类免疫缺陷病毒(HIV)中的蛋白质-RNA相互作用和催化RNA分子中的PKA计算。在第一个区域，研究了自然发生的、转录后修饰的碱基在影响tRNA-mRNA识别中的作用。在人类tRNALys，3中，我们发现37位的修饰碱基是维持反密码子碱基(34-36)中典型的阶梯状构象所必需的。使用M05-2X官能团的自然键轨道分析的从头算研究正在进行中，以确定潜在的稳定力和第37位修饰碱基在所有tRNA中保持阶梯状构象的作用。需要在M05-2X/6-31+G(d，p)理论水平上对四核苷酸和三核苷酸(二聚体为~1400个基函数)进行氢位优化，在我们的本地机器上，这可能需要45天以上的计算时间。要在这个项目上取得进展，需要更快的计算资源。在第二个研究领域，我们正在研究水和静电在RNA-肽识别中的作用。在晚期，REV-RRE识别介导部分和未剪接的HIV mRNA的核质输出。从Frankel和他的同事进行的体外筛选研究中，一种被称为RSG-1.2的合成肽被发现比天然的REV肽具有更高的亲和力和特异性结合RRE。我们用Amber模拟了REV和RSG-1.2多肽在显性水中与RRE RNA的络合作用，发现了多肽-RNA络合物中的水结构与结合亲和力之间的关系。需要更多的模拟来证实早先的发现。系统大约有35,000个原子，使用并行的琥珀代码可以更有效地收集数据。最后，与达林·约克合作，我们正在计算被称为核酶的催化RNA分子中的pKas。热力学积分方法需要平衡的启动系统。目前的体系是在显式溶剂(TIP4Pew)中进行的，包括中和RNA以外的150 mM NaC l缓冲溶液，大约有75,000个原子。这些系统需要多次模拟退火轮来平衡离子大气，然后在执行TI计算之前，必须随后在缓冲液存在的情况下平衡RNA。请求这种分配是为了利用并行计算设施，同时还为未来的分配请求探索最佳的TeraGrid平台。

This subproject is one of many research subprojects utilizing the resources provided by a Center grant funded by NIH/NCRR. Primary support for the subproject and the subproject's principal investigator may have been provided by other sources, including other NIH sources. The Total Cost listed for the subproject likely represents the estimated amount of Center infrastructure utilized by the subproject, not direct funding provided by the NCRR grant to the subproject or subproject staff. Accurate RNA recognition by other biomolecules such as proteins, cofactors and other RNA molecules are critical to many cellular functions. Employing a variety of computational chemistry tools such as molecular dynamics simulations, quantum calculations and hybrid quantum mechanical/molecular mechanics methods, our research examines three primary areas, messenger RNA - transfer RNA (mRNA-tRNA) recognition, a key step in the translation of proteins, protein-RNA interactions in human immunodeficiency virus (HIV) and pKa calculations in catalytic RNA molecules. In the first area, the role of naturally occurring, posttranscriptionally modified bases in affecting tRNA-mRNA recognition is examined. In human tRNALys,3, we have found that a modified base at position 37 are required for maintenance of a canonical stair-stepped conformation in the anticodon bases (34-36). Ab initio studies employing natural bond orbital analysis with the M05-2X functional are underway to determine the underlying stabilizing forces and the role of modified bases at the 37th position in retaining a stair-stepped conformation in all tRNAs. Optimization of hydrogen positions at the M05-2X/6-31+G(d,p) theory level needs to be carried out for tetranucleotides and trinucleotides (dimers are ~1400 basis functions), which on our local machines can take greater than 45 days/calculation. Faster computing resources are required to make progress on this project. In the second area of research, we are examining the role of water and electrostatics in RNA-peptide recognition. In late phase Rev-RRE recognition mediates nucleocytoplasmic export of partially and unspliced HIV mRNA. From in vitro selection studies performed by Frankel and coworkers, a synthetic peptide known as RSG-1.2 has been found to bind RRE with greater affinity and specificity than the native Rev peptide. We have simulated both Rev and RSG-1.2 peptides complexed with the RRE RNA in explicit water using AMBER and have found a correlation between water structure in the peptide-RNA complexes and binding affinity. More simulations to corroborate earlier findings are required. Systems are roughly 35,000 atoms and data could be collected more efficiently employing parallel AMBER code. Lastly, in collaboration with Darrin York, we are calculating pKas in catalytic RNA molecules known as ribozymes. The thermodynamic integration methods require equilibrated starting systems. Current systems are carried out in explicit solvent (TIP4Pew), include 150 mM NaCl buffer solution beyond the neutralized RNA and are about 75,000 atoms. These systems require a number of simulated annealing rounds to equilibrate the ion atmosphere and then the RNA must be subsequently equilibrated in the presence of the buffer before TI calculations can be performed. This allocation is requested to take advantage of parallel computing facilities while also exploring optimum Teragrid platforms for future allocation requests.