核酸–小分子相互作用的打分函数和分子对接算法研究

来源 :华中科技大学 | 被引量 : 0次 | 上传用户:zjh73
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
核酸是构成细胞的基本物质之一,其基本单元是核苷酸。核酸作为遗传信息的载体,在信息的存储和传递过程中行使着关键功能,对生物的生成、发育和繁殖等活动有着至关重要的作用。对核酸的研究是当前对生命科学研究的一个重要分支,其中,有一类核酸可以作为药物靶标,例如核糖开关,与配体结合,从而调节基因表达,激活/失活生物学功能。那么就有必要研究核酸–配体相互作用,从而加深对核酸功能的理解。结构决定功能,但由于实验方法昂贵的经济成本和技术难度,经实验确定核酸–配体复合物结构的数量十分有限。分子对接等计算方法因其便捷和较少的花费而受到了广泛关注,在核酸–配体复合物结构建模方面发挥着互补作用。在分子对接中,打分函数通过评估配体分子与核酸结合的紧密性和特异性来指导复合物构象的生成和筛选,打分函数精确与否会直接影响建模结果的好坏。基于玻尔兹曼统计,我们训练了一组基于知识的打分函数,命名为ITScore-NL,用于计算核酸–小分子相互作用能量,并在打分函数中,显性的考虑两个分子之间的堆叠和静电相互作用。堆叠是核酸结构的主要特征之一,也是稳定核酸结构的一种关键作用。在一定距离范围内,当核酸碱基与配体芳香基团处于近似平行的两个平面时,就会存在着堆叠作用。此外,由于核酸的高电荷性,静电相互作用对核酸–小分子配体的结合尤为重要。因此,也在构建打分函数时,引入静电项,以减少一些原子对在某些距离上,因频率不足而造成径向分布函数中的统计误差。从结果中可以看到,堆叠势和静电相互作用对提升打分函数性能的贡献是不同的。加入堆叠项主要提高了结合亲和力的预测精度,相关系数得到提升。而加入静电项主要提高了ITScore在结构筛选中的能力,在宽松和严格的标准下,复合物结构的筛选结果均取得了改善。受启发于我们课题组新训练的打分函数(ITScore-NL)在指导复合物结构建模中的出色表现,进而开发了一种新的RNA/DNA–小分子对接算法,称为NLDock。新的对接方案是对蛋白质–小分子配体对接算法MDock的重大修改。NLDock通过对接多个小分子配体构象来近似模拟配体的柔性。测试结果表明,在局部小分子配体刚性对接和柔性对接方面,NLDock对接成功率均优于其它三种对接方案。NLDock的优异表现可能与ITScore-NL打分函数在RNA/DNA–小分子结合亲和力预测方面的优势有关。同时,在全局对接方面,即结合位点未知情况下,采取基于几何标准,对可能的小分子结合区域进行筛选,在预测的结合区域内再进行局部对接,从而实现全局对接目的。从对接的结果来看,NLDock依然要好于其它对接软件。此外,对接效率也是NLDock的一大优势,平均的对接时间比其它对接软件要快将近一个量级。在未知结合位点的情况下进行分子对接可以称之为盲对接,此时,需要基于个人经验或者实验提供的信息去定义候选配体可能的结合区域,这对于初次接触对接领域的人来说是比较困难的。盲对接一直以来是一项比较有挑战性的工作,也是分子对接过程中所面临最为直接的问题。为了准确的筛选核酸–配体结合位点,我们开发了一种位点预测算法,称之为NLsite。NLsite是一种新颖的,整合几何和能量预测RNA–配体结合位点的方法。首先,RNA受体被放置在一个三维立体网格中,位于RNA结构凹陷区域的格点被筛选出来;其次,在剩余每个格点上放置已给定原子类型的球体,然后,由所训练的原子对势(ITScore-NL),计算不同类型球体位于每个格点的能量。对格点进行能量累加并再次筛选。最后,基于距离映射的方式,将格点能量投影到每个核苷酸上,如此,每个核苷酸都有了其相应的分数,设置一定的分数阈值,定义分数阈值以内的核苷酸是位于结合位点区域的核苷酸。结果表明,NLsite是一个有效的位点预测方法,经过在五个测试集上测试并与其它三种核酸–配体位点预测工具进行比较,可以发现NLsite在四个指标上预测的准确性均优于现有的传统方法。
其他文献
在生物医学和传染病学的研究中,对某些事件的失效时间及其影响因素的研究是普遍存在的.然而,由于经费限制或者过高的技术难度要求常造成对失效时间有主要影响的某些协变量的观测数据有缺失.但在很多情况下,可能有一些与缺失的主要协变量相关的变量却由于测量简单便宜而可以得到几乎完整地搜集.这些容易获得的变量可以为主要变量的情况提供辅助,而在主要变量完整时,又不能对失效时间规律的研究提供额外的信息,因此被称为辅助
学位
第一部分HNF1A在食管癌组织中的表达及其与临床病理特征和预后的关系目的:对于不可手术切除的食管鳞癌患者,根治性放化疗是重要的治疗方式。然而,部分患者因放射抵抗导致放射治疗后预后生存较差。本研究分析HNF1A在食管癌组织中的表达情况并探讨其与临床病理特征和预后的关系,为提高放射敏感性和靶向治疗提供新思路。方法:基于生物信息学综合分析HNF1A基因在食管癌组织标本中的表达水平,探讨其与临床病理特征及
学位
髌骨骨折是在骨科门诊中常见的外伤,通常由于突然撞击髌骨而发生,例如跌倒或交通事故。髌骨位于膝关节前方皮下,对膝关节伸肌机制至关重要,在膝关节伸展过程中,髌骨是一个支点,可以提高伸肌肌腱的力矩力臂。当髌骨骨折发生时,力传递功能丧失,膝关节无法正常活动。手术固定髌骨骨折,可以重建伸肌机制的完整性,恢复膝关节屈伸活动,进行早期康复锻炼,是髌骨骨折最主要的治疗手段。目前,固定髌骨骨折的主流技术是张力带技术
学位
成纤维细胞生长因子受体(Fibroblast growth factor receptor,FGFR)现已成为肿瘤精准诊疗领域的重要靶点,在多种肿瘤中呈现较高水平的表达。FGFR主要分为FGFR1-4四种亚型,成纤维细胞生长因子(Fibroblast growth factor,FGF)与之结合后,形成FGF/FGFR信号通路,在肿瘤发生、生存、侵袭、血管生成、转移、复发、上皮间质转化(Epith
学位
核磁共振能够无损地精准获取物质组成和结构信息,是一种重要的物质分析技术。基于感应线圈的核磁共振谱仪的信噪比正比于γ5/2B03/2,工作磁场强度越高,则核磁共振信号的信噪比越高。因此,传统的核磁共振谱仪通常采用超导磁体产生强的工作磁场。然而,基于超导磁体的高场核磁共振波谱仪存在体积大、成本高、匀场难等问题,特别是匀场难,限制了核磁共振谱的频谱测量精度的提高。超低场核磁共振的工作磁场在超低场(<10
学位
静电加速度计具有极高的分辨率,并且可实现六自由度加速度同步测量,在重力卫星,空间引力波探测等高精度空间实验中发挥重要作用。在地面对静电加速度计的性能进行直接测试与验证,对于提升其可靠性具有重要的意义。基于高压悬浮的地面测试方案能实现加速度计检验质量的六自由度同时控制,对待测样机改动相对较小,适合对工程样机整机进行检验。然而在该方案中,水平方向的地面振动噪声会直接耦合进入加速度计,影响测试结果。基于
学位
晶体与电磁场之间相互作用的研究是固体物理领域极具吸引力的课题之一。强度不同的电磁场会在晶体内部诱导产生不同的物理现象。在过去几十年中,人们尽管对于弱激光场与晶体的相互作用进行了大量的研究,并取得了显著的成果,但是对于强激光场与晶体相互作用的研究仍存在很多亟需解决的问题。近年来,得益于飞秒激光技术的发展与成熟,非线性光学应用等得到了快速发展,认识和理解晶体尤其是半导体内的强场过程则变得尤为重要。半导
学位
硅基MEMS加速度计在消费电子、惯性导航和地球资源勘探中都有很广泛的应用。国内外各大研究机构都对其展开了深入的研究,相继提出了不同类型的MEMS加速度计,其中电容检测式MEMS加速度计发展最为迅速。在加速度计的设计、加工、可靠性、测试以及应用等技术中,加速度计的封装技术是其重要的技术环节。无论是从消费电子需求的低成本和批量化,还是惯性导航和资源勘探需求的加速度计的一致性,圆片级封装都是必须要解决的
学位
第一部分基于转录组测序技术的综合生物信息学分析揭示髌骨不稳相关的关键基因和通路目的:髌骨不稳是一种青少年常见的膝关节损伤,然而到目前为止,其关键的生物标志物和分子机制仍不清楚。本研究目的为通过转录组基因测序技术比较髌骨不稳和对照组小鼠之间的全基因表达差异,发现参与髌骨不稳发生发展的显著差异基因,并通过综合生物信息学分析确定了髌骨不稳中的分子机制和关键基因,以揭示髌骨不稳的确切机制并开发一种新的非手
学位
血脂紊乱作为一种常见的代谢性疾病,与动脉粥样硬化性心血管疾病(arteriosclerotic cardiovascular disease,ASCVD)的发生发展密切相关。在临床工作中,即使对于心血管疾病的危险因素进行积极控制,仍会有部分患者出现大血管或微血管疾病,这就是所谓的心血管剩留风险。这种心血管剩留风险主要与高TG血症有关,通过有效的手段降低TG水平则可以极大的减少心血管疾病的患病风险。
学位