基于机器学习的非编码RNA-蛋白质相互作用关系预测研究：方法、实践与展望

基于机器学习的非编码RNA-蛋白质相互作用关系预测研究：方法、实践与展望一、引言1.1研究背景在生命科学领域，非编码RNA（ncRNA）与蛋白质的相互作用是生物过程的关键环节，对维持细胞正常功能和生物机体的稳态起着不可或缺的作用。随着后基因组时代的到来，人们逐渐认识到非编码RNA并非基因组中的“暗物质”，而是在基因表达调控、细胞分化、发育进程以及疾病发生发展等诸多方面发挥着重要的调控功能。非编码RNA是指不编码蛋白质的RNA分子，虽然它们不直接参与蛋白质的合成，但通过与蛋白质、DNA及其他RNA分子相互作用，形成了复杂的调控网络。这些相互作用广泛存在于转录前、转录及转录后等各个水平，对基因表达进行精细调控。例如，微小RNA（miRNA）能够通过与靶mRNA的互补配对，抑制mRNA的翻译过程或促使其降解，从而在转录后水平调控基因表达；长链非编码RNA（lncRNA）可以通过与蛋白质形成复合物，招募或阻断转录因子与DNA的结合，进而影响基因转录；环状RNA（circRNA）则能够作为miRNA的“海绵”，竞争性结合miRNA，解除miRNA对其靶基因的抑制作用，间接调控基因表达。这些复杂的调控机制在细胞的增殖、分化、凋亡以及代谢等基本生命活动中均发挥着重要作用，确保了细胞功能的正常行使和生物个体的生长发育。蛋白质作为生命活动的主要执行者，参与了细胞内几乎所有的生理过程，包括信号转导、物质代谢、免疫防御等。非编码RNA与蛋白质的相互作用赋予了蛋白质更为多样化的功能和调控机制。一方面，非编码RNA可以引导蛋白质定位到特定的细胞区域或作用靶点，使其能够准确地发挥生物学功能。例如，某些lncRNA能够与特定的蛋白质结合，将其招募到细胞核内的特定染色体区域，参与基因的转录调控。另一方面，非编码RNA还可以通过改变蛋白质的构象或活性，影响其功能的发挥。例如，一些miRNA与蛋白质结合后，能够调节蛋白质的磷酸化水平或蛋白质-蛋白质相互作用，从而改变蛋白质的生物学活性。在疾病发生发展过程中，非编码RNA与蛋白质相互作用的失调往往起着关键作用。研究表明，许多人类疾病，如癌症、心血管疾病、神经系统疾病等，都与非编码RNA-蛋白质相互作用网络的异常密切相关。在癌症中，某些miRNA与癌基因或抑癌基因编码的蛋白质相互作用异常，导致细胞增殖、凋亡和转移等过程的紊乱；在心血管疾病中，lncRNA与心血管相关蛋白质的相互作用失调，可能影响心肌细胞的功能和心血管系统的稳态；在神经系统疾病中，circRNA与神经退行性疾病相关蛋白质的异常相互作用，可能参与神经细胞的损伤和死亡过程。深入研究这些异常相互作用的机制，不仅有助于揭示疾病的发病机制，还为疾病的早期诊断、预后评估和治疗提供了新的靶点和策略。传统上，研究非编码RNA-蛋白质相互作用主要依赖于实验方法，如免疫共沉淀（Co-IP）、RNA免疫沉淀（RIP）、交联免疫沉淀（CLIP）等。这些实验方法虽然能够直接检测到相互作用的存在，但存在成本高、耗时长、通量低等局限性，难以满足大规模研究的需求。此外，实验方法还受到样本来源、实验条件等因素的限制，对于一些低丰度或瞬时相互作用的检测存在一定困难。因此，开发高效、准确的计算方法来预测非编码RNA-蛋白质相互作用关系具有重要的理论和实际意义。机器学习作为人工智能领域的重要分支，近年来在生物信息学中得到了广泛应用。通过构建合适的机器学习模型，可以从大量的生物数据中挖掘出潜在的模式和规律，实现对非编码RNA-蛋白质相互作用关系的有效预测。机器学习方法具有高效、快速、可扩展性强等优点，能够弥补实验方法的不足，为非编码RNA-蛋白质相互作用的研究提供了新的手段。利用机器学习算法对非编码RNA和蛋白质的序列、结构、表达等特征进行分析和建模，可以预测它们之间是否存在相互作用以及相互作用的强度和特异性。这种计算预测方法不仅能够为实验研究提供有价值的线索和指导，还能够加速对非编码RNA-蛋白质相互作用网络的解析，推动生命科学的发展。1.2研究目的与意义本研究旨在运用机器学习方法，构建高效准确的预测模型，以识别非编码RNA与蛋白质之间的相互作用关系。通过整合多种生物信息学数据和特征，深入挖掘非编码RNA与蛋白质相互作用的潜在规律，为生命科学领域的相关研究提供新的理论支持和技术手段。从基础研究角度来看，深入了解非编码RNA-蛋白质相互作用关系有助于完善我们对生命基本过程的认知。非编码RNA与蛋白质相互作用参与了众多关键生物过程，如基因转录调控、RNA剪接加工、蛋白质翻译等。揭示这些相互作用的机制，能够帮助我们更好地理解细胞内遗传信息的传递和调控网络，进一步阐明生命活动的本质。以转录调控为例，许多转录因子与非编码RNA相互作用，共同调节基因的转录起始、延伸和终止过程。准确识别这些相互作用关系，有助于我们绘制更加精确的基因转录调控图谱，为研究细胞分化、发育以及疾病发生发展过程中的基因表达变化提供重要线索。在疾病研究方面，非编码RNA-蛋白质相互作用的异常与多种人类疾病的发生发展密切相关，包括癌症、心血管疾病、神经退行性疾病等。通过预测这些相互作用关系，可以为疾病的早期诊断、预后评估和治疗提供新的靶点和生物标志物。在癌症研究中，某些miRNA与癌基因或抑癌基因编码的蛋白质相互作用异常，可能导致肿瘤细胞的增殖、侵袭和转移。通过机器学习预测这些异常相互作用，有助于发现新的癌症诊断标志物和治疗靶点，为癌症的精准诊断和个性化治疗提供依据。对于心血管疾病，研究非编码RNA-蛋白质相互作用网络的变化，能够帮助我们揭示疾病的发病机制，寻找潜在的治疗干预点，为开发新型心血管疾病治疗药物提供理论基础。药物研发是另一个重要应用领域。目前，大多数药物研发主要针对蛋白质靶点，但非编码RNA-蛋白质相互作用作为一类新的药物作用靶点，具有巨大的潜力。预测非编码RNA-蛋白质相互作用关系，可以为药物研发提供新的思路和方向，加速新型药物的开发进程。通过筛选能够干扰异常非编码RNA-蛋白质相互作用的小分子化合物或生物制剂，有望开发出针对特定疾病的创新治疗药物。此外，基于机器学习预测结果，可以对潜在药物靶点进行优先级排序，提高药物研发的效率和成功率，降低研发成本。机器学习方法在预测非编码RNA-蛋白质相互作用关系中具有独特的优势。它能够处理大规模、高维度的生物数据，快速准确地挖掘数据中的潜在模式和规律。与传统实验方法相比，机器学习方法具有高效、低成本、高通量等特点，能够在短时间内对大量非编码RNA和蛋白质对进行预测，为实验研究提供有价值的候选对象，从而加速非编码RNA-蛋白质相互作用的研究进程。本研究通过机器学习方法预测非编码RNA-蛋白质相互作用关系，不仅具有重要的理论意义，能够推动生命科学基础研究的发展，还具有广泛的应用前景，将为生物医学研究和药物研发带来新的机遇和突破，有望为人类健康事业做出积极贡献。1.3国内外研究现状近年来，随着机器学习技术的飞速发展，其在预测非编码RNA-蛋白质相互作用关系领域的应用日益广泛，吸引了国内外众多科研团队的关注，取得了一系列具有重要意义的研究成果。在国外，一些研究团队致力于开发基于机器学习的预测模型，以揭示非编码RNA与蛋白质之间复杂的相互作用关系。2022年，早稻田大学的研究人员[具体文献]提出了BERT-RBP模型，通过对人类参考基因组上预训练的BERT架构进行调整，来预测RNA序列与RNA结合蛋白（RBPs）之间的相互作用。该模型在使用154个RBPs的eCLIP-seq数据时，展现出优于当时最先进预测模型的性能。其创新之处在于应用在大型语料库上预训练的BERT模型，提高了模型的可解释性，仅从RNA序列信息就能识别转录本区域类型和RNA二级结构，为预测RNA-RBP相互作用提供了新的思路和方法。2025年，华中师范大学物理研究所与美国华盛顿大学和达特茅斯学院合作[具体文献]，推出了ZHMolGraph模型。该模型集成了图神经网络和无监督大型语言模型来预测RNA-蛋白质相互作用，在两个基准数据集上验证时超越了当前的最佳方法。对于完全未知的RNA和蛋白质的数据集，ZHMolGraph在实现79.8%的高AUROC和82.0%的AUPRC方面有所提高，成为全基因组RNA-蛋白质预测的可靠选择，在建模和设计RNA-蛋白质复合物方面具有广泛潜力。国内学者在该领域也开展了深入研究，并取得了显著进展。上海第二工业大学的研究团队[具体文献]提出一种基于伪特定位点记分矩阵（PsePSSM）和几何保留投影（GPP）降维算法的方法。首先采用PsePSSM算法表达蛋白质序列，编码后的蛋白质特征向量维数很高，接着采用GPP流形学习方法对其进行维数约简，约简后的特征向量输入支持向量机（SVM）分类器训练，用于预测未知的蛋白质与RNA是否相互作用。实验结果表明，该方法能有效解决传统实验方法的问题，为蛋白质与RNA相互作用的研究提供了新途径。辽宁科技大学的赵琪教授团队[具体文献]应用整合策略方法研究长链非编码RNA与蛋白质相互作用。该策略将多种机器学习算法基于多种特征组合方案训练的模型进行整合，相比现有模型具有更广泛的适用性，能够更加全面地揭示潜在的长链非编码RNA与蛋白质相互作用关系。同时，利用随机配对方法构建负数据集，进一步降低了模型的假阳性率，提高了模型预测的准确性和可靠性。尽管国内外在利用机器学习预测非编码RNA-蛋白质相互作用关系方面取得了诸多成果，但当前研究仍存在一些不足和空白。一方面，现有的机器学习模型大多依赖于特定的数据集和特征选择，模型的泛化能力和跨物种适用性有待提高。不同物种的非编码RNA和蛋白质序列、结构及相互作用机制存在差异，如何开发能够有效处理不同物种数据的通用模型，仍是亟待解决的问题。另一方面，对于非编码RNA-蛋白质相互作用的动态变化和环境因素的影响，目前的研究还相对较少。在生物体内，非编码RNA-蛋白质相互作用受到多种因素的调控，如细胞状态、生理病理条件和外界环境刺激等。深入研究这些动态变化和影响因素，对于全面理解非编码RNA-蛋白质相互作用的生物学功能至关重要，但相关研究尚处于起步阶段。此外，虽然机器学习方法在预测相互作用关系方面取得了一定成效，但对于预测结果的生物学验证和功能解析还不够深入。如何将机器学习预测结果与实验验证相结合，进一步阐明非编码RNA-蛋白质相互作用在生物过程中的具体作用机制，也是未来研究需要重点关注的方向。同时，目前的研究主要集中在常见的非编码RNA类型，如miRNA、lncRNA等，对于一些新型非编码RNA与蛋白质相互作用的研究还相对匮乏，存在较大的探索空间。二、非编码RNA-蛋白质相互作用关系概述2.1非编码RNA的分类与功能非编码RNA（ncRNA）是一类不编码蛋白质的RNA分子，在生物体内广泛存在且发挥着至关重要的调控作用。随着研究的不断深入，越来越多的非编码RNA被发现并鉴定，根据其长度、结构和功能的差异，主要可分为微小RNA（miRNA）、长链非编码RNA（lncRNA）、环状RNA（circRNA）以及其他一些小分子非编码RNA等，它们在基因表达调控、细胞分化、发育以及疾病发生发展等过程中扮演着不可或缺的角色。微小RNA（miRNA）是一类长度约为21-23个核苷酸的内源性非编码单链RNA分子。其生物合成过程较为复杂，首先在细胞核内由RNA聚合酶II转录生成初级miRNA（pri-miRNA），pri-miRNA经过核酸酶Drosha及其辅助因子DGCR8的加工，形成长度约为70-100个核苷酸的发夹结构前体miRNA（pre-miRNA）。随后，pre-miRNA被转运蛋白Exportin-5转运至细胞质中，再由核酸酶Dicer切割生成成熟的miRNA。成熟的miRNA会与AGO蛋白等组装形成RNA诱导沉默复合体（RISC），通过与靶mRNA的3'非翻译区（3'-UTR）特异性互补配对，抑制mRNA的翻译过程，或者促使靶mRNA降解，从而在转录后水平对基因表达进行精细调控。例如，在细胞增殖过程中，miR-21通过抑制其靶基因PTEN的表达，激活PI3K/AKT信号通路，促进细胞的增殖和存活；在细胞分化过程中，miR-124则通过靶向调控一系列神经发育相关基因的表达，促进神经干细胞向神经元的分化。miRNA还广泛参与了生物体的代谢、免疫、凋亡等生理过程，并且在多种疾病，如癌症、心血管疾病、神经系统疾病等的发生发展中发挥着重要作用，其表达水平的异常往往与疾病的进程和预后密切相关。长链非编码RNA（lncRNA）是长度大于200个核苷酸的非编码RNA分子。lncRNA的转录过程与mRNA类似，主要由RNA聚合酶II催化转录产生，但与mRNA不同的是，lncRNA通常不具有完整的开放阅读框，不编码蛋白质。lncRNA的结构较为复杂，具有多种二级和三级结构，并且其在基因组上的分布广泛，可位于基因间区、内含子区域、外显子区域或与编码基因重叠。lncRNA通过与DNA、RNA和蛋白质相互作用，在多个层面上参与基因表达调控。在染色质水平，lncRNA可以招募染色质修饰复合物，如多梳蛋白复合物（PRC）等，对染色质进行修饰，从而影响基因的转录活性。例如，XISTlncRNA在X染色体失活过程中发挥关键作用，它通过与PRC2复合物结合，使X染色体上的基因发生组蛋白H3赖氨酸27三甲基化修饰（H3K27me3），导致基因沉默，实现X染色体的失活。在转录水平，lncRNA可以与转录因子或RNA聚合酶II相互作用，促进或抑制基因的转录。在转录后水平，lncRNA可以与mRNA结合，影响mRNA的稳定性、剪接和转运等过程。此外，lncRNA还参与了细胞周期调控、细胞分化、胚胎发育以及肿瘤发生发展等多种生物学过程。在肿瘤研究中发现，许多lncRNA在肿瘤组织中呈现特异性表达，如HOTAIRlncRNA在乳腺癌、结直肠癌等多种肿瘤中高表达，它通过调控相关基因的表达，促进肿瘤细胞的增殖、侵袭和转移。环状RNA（circRNA）是一类具有共价闭合环状结构的非编码RNA分子，其形成过程主要通过反向剪接机制实现，即上游外显子的供体端与下游外显子的受体端反向连接，形成闭环结构。circRNA由于其特殊的环状结构，不具有5'端帽子和3'端poly（A）尾巴，因此对核酸酶具有较高的抗性，在细胞内相对稳定，半衰期较长。circRNA主要定位于细胞质中，少数也可存在于细胞核内。circRNA的功能具有多样性，其中较为突出的是作为miRNA的“海绵”，通过竞争性结合miRNA，解除miRNA对其靶基因的抑制作用，从而间接调控基因表达。例如，ciRS-7（也称为CDR1as）含有大量与miR-7互补的结合位点，能够吸附miR-7，上调miR-7靶基因的表达，在神经系统发育和肿瘤发生发展中发挥重要作用。此外，circRNA还可以与蛋白质相互作用，影响蛋白质的功能和定位；部分circRNA还能够参与转录调控，通过与DNA结合，调节基因的转录起始和延伸过程。近年来的研究表明，circRNA在多种疾病，尤其是神经系统疾病、心血管疾病和癌症中呈现异常表达，有望成为疾病诊断和治疗的新型生物标志物和潜在靶点。除了上述几种常见的非编码RNA外，还有一些小分子非编码RNA也在生物过程中发挥着重要作用。小核RNA（snRNA）主要存在于细胞核内，长度一般在100-300个核苷酸之间，它与蛋白质结合形成小核核糖蛋白颗粒（snRNPs），参与mRNA前体的剪接过程，确保正确的mRNA序列被转录和翻译，对于维持基因表达的准确性和细胞稳态具有重要意义。小核仁RNA（snoRNA）主要位于核仁中，长度约为60-300nt，分为C/D盒snoRNA和H/ACA盒snoRNA两类，主要负责rRNA的转录后修饰，如甲基化修饰和假尿苷化修饰，这些修饰对于rRNA的成熟和核糖体的组装至关重要。转运RNA（tRNA）是一类长度约为70-90个核苷酸的小分子RNA，其结构呈三叶草型，具有携带氨基酸的功能，在蛋白质合成过程中，tRNA通过反密码子与mRNA上的密码子互补配对，将相应的氨基酸转运至核糖体，参与多肽链的合成。核糖体RNA（rRNA）是核糖体的主要组成部分，直接参与核糖体中蛋白质的合成过程，起到识别、选择tRNA以及催化肽键形成等多种作用。此外，还有小干扰RNA（siRNA）、Piwi相互作用RNA（piRNA）等，它们在基因沉默、转座子调控、维持基因组稳定性等方面发挥着重要作用。非编码RNA的种类丰富多样，各自具有独特的结构和功能，它们通过复杂的调控网络，在生物体内的基因表达调控、细胞生理过程以及疾病发生发展等方面发挥着关键作用，深入研究非编码RNA的分类与功能，对于揭示生命奥秘和攻克人类疾病具有重要的理论和实际意义。2.2蛋白质的结构与功能蛋白质是由氨基酸通过肽键连接而成的生物大分子，其结构复杂多样，具有多个层次，从一级结构到四级结构，每一个层次的结构都对蛋白质的功能起着至关重要的作用。这些复杂的结构赋予了蛋白质多种多样的生物学功能，使其成为生命活动的主要执行者。蛋白质的一级结构是指多肽链中氨基酸的排列顺序，这是蛋白质最基本的结构层次，由基因中的核苷酸序列决定。不同的氨基酸通过肽键依次相连，形成线性的多肽链。例如，血红蛋白的α-链由141个氨基酸组成，β-链由146个氨基酸组成，它们特定的氨基酸排列顺序决定了血红蛋白独特的功能。一级结构中的氨基酸序列包含了蛋白质折叠和形成高级结构所需的全部信息，氨基酸序列的改变可能会导致蛋白质功能的异常，如镰状细胞贫血就是由于血红蛋白β-链上的一个氨基酸由谷氨酸变为缬氨酸，导致血红蛋白的结构和功能发生改变，红细胞变形能力下降，易破裂，从而引发贫血症状。蛋白质的二级结构是指多肽链主链原子的局部空间排列，不涉及侧链的构象。常见的二级结构有α-螺旋、β-折叠、β-转角和无规卷曲等。α-螺旋呈右手螺旋状，每3.6个氨基酸残基上升一圈，螺距为0.54nm，肽链中的酰胺氢和羰基氧之间形成氢键，以维持螺旋结构的稳定。许多纤维状蛋白质，如角蛋白，富含α-螺旋结构，赋予毛发、指甲等坚韧的特性。β-折叠是由若干条肽链或肽段平行排列，通过链间的氢键维系而成的片层结构，分为平行式和反平行式两种。蚕丝中的丝心蛋白主要由β-折叠结构组成，使其具有柔软、光滑的特性。β-转角通常由4个氨基酸残基组成，常出现在多肽链的回折处，起到连接不同二级结构的作用。无规卷曲则是指多肽链中没有规律的松散部分，其结构较为灵活，可参与蛋白质与其他分子的相互作用。在二级结构的基础上，多肽链进一步折叠卷曲，形成更为复杂的三维空间结构，即蛋白质的三级结构。三级结构主要由非共价键（如氢键、离子键、疏水作用、范德华力等）以及二硫键来维持。蛋白质的三级结构决定了其整体形状和表面特征，使其能够特异性地结合底物、配体或其他蛋白质分子。例如，肌红蛋白是一种球状蛋白质，其三级结构包含一个疏水核心，内部包裹着血红素辅基，能够可逆地结合氧气，为肌肉组织储存和运输氧气。酶蛋白的三级结构具有高度特异性的活性中心，能够与底物分子精确匹配，催化特定的化学反应。对于由两条或两条以上具有独立三级结构的多肽链组成的蛋白质，这些多肽链之间通过非共价键相互作用，形成特定的空间排布，这种结构被称为蛋白质的四级结构。组成四级结构的每条多肽链称为亚基，亚基单独存在时一般不具有生物学活性，只有通过特定的方式组装成完整的四级结构后，蛋白质才具有完整的生物学功能。例如，血红蛋白由4个亚基（2个α-亚基和2个β-亚基）组成，4个亚基之间通过离子键、氢键等相互作用形成稳定的四级结构。在与氧气结合的过程中，4个亚基之间存在着协同效应，当第一个亚基与氧气结合后，会引起整个血红蛋白分子构象的改变，使得其他亚基对氧气的亲和力增强，从而更有利于血红蛋白在肺部高效地结合氧气，并在组织中释放氧气，满足机体对氧的需求。蛋白质在生命活动中承担着多种多样的功能，是维持生物体正常生理活动的关键物质。在催化功能方面，绝大多数酶都是蛋白质，它们能够显著降低化学反应的活化能，加速生物体内的各种化学反应。例如，淀粉酶能够催化淀粉水解为葡萄糖，蛋白酶可以将蛋白质分解为氨基酸，这些酶促反应对于生物体的消化、吸收和代谢过程至关重要。在物质运输功能方面，一些蛋白质能够结合并运输特定的物质，如血红蛋白负责运输氧气，将氧气从肺部输送到全身各个组织细胞；血清白蛋白可以结合和运输脂肪酸、胆红素等多种小分子物质，维持体内物质的平衡和正常代谢。在结构支持方面，许多蛋白质是构成细胞和组织的重要结构成分，如胶原蛋白是结缔组织（如皮肤、骨骼、肌腱等）的主要成分，赋予组织韧性和强度；角蛋白构成毛发、指甲、羽毛等，起到保护和支撑的作用。在免疫防御功能方面，抗体是一类重要的免疫球蛋白，能够特异性地识别和结合外来病原体（如细菌、病毒等），通过免疫反应清除病原体，保护机体免受感染。在信号传导功能方面，细胞表面的受体蛋白能够识别并结合细胞外的信号分子（如激素、神经递质等），将信号传递到细胞内，引发一系列的细胞内信号转导事件，调节细胞的生理活动，如胰岛素受体与胰岛素结合后，启动细胞对葡萄糖的摄取和利用，调节血糖水平。在基因表达调控方面，一些蛋白质作为转录因子，能够结合到DNA的特定区域，调节基因的转录过程，控制蛋白质的合成，从而影响细胞的分化、发育以及对环境变化的响应。蛋白质的结构和功能密切相关，其复杂的结构层次赋予了蛋白质丰富多样的生物学功能，使其在生命活动的各个方面都发挥着不可或缺的作用。深入研究蛋白质的结构与功能，对于理解生命过程的本质、揭示疾病的发病机制以及开发新型治疗药物具有重要的意义。2.3相互作用关系的生物学意义非编码RNA与蛋白质之间的相互作用在生物体内具有广泛而深远的生物学意义，它们参与了基因表达调控、细胞周期调控、疾病发生发展等多个关键生物学过程，对维持细胞的正常功能和生物机体的稳态起着不可或缺的作用。在基因表达调控方面，非编码RNA-蛋白质相互作用构成了一个复杂而精细的调控网络，在转录前、转录及转录后等各个水平对基因表达进行精准调控。在转录前水平，一些非编码RNA可以与蛋白质结合形成复合物，招募或阻碍转录因子与DNA的结合，从而影响基因转录的起始。例如，增强子RNA（eRNA）可以与转录激活因子结合，增强转录因子与增强子区域的结合能力，促进基因的转录起始；而某些长链非编码RNA（lncRNA）则可以通过与抑制性蛋白质结合，形成转录抑制复合物，阻止转录因子与启动子区域的结合，抑制基因的转录。在转录水平，非编码RNA与蛋白质的相互作用可以影响RNA聚合酶的活性和转录延伸的速率。例如，一些lncRNA能够与RNA聚合酶II结合，调节其在DNA模板上的移动速度，从而影响基因转录的效率。在转录后水平，非编码RNA-蛋白质相互作用更是发挥着重要作用。微小RNA（miRNA）通过与AGO蛋白等组装形成RNA诱导沉默复合体（RISC），与靶mRNA的3'非翻译区（3'-UTR）特异性互补配对，抑制mRNA的翻译过程或促使其降解，实现对基因表达的转录后调控。此外，一些lncRNA和环状RNA（circRNA）也可以通过与mRNA结合，影响mRNA的稳定性、剪接和转运等过程，进而调控基因表达。例如，某些lncRNA可以与mRNA形成双链结构，保护mRNA免受核酸酶的降解，延长其半衰期；circRNA则可以作为miRNA的“海绵”，竞争性结合miRNA，解除miRNA对其靶基因的抑制作用，间接调控基因表达。细胞周期调控是细胞生命活动的重要过程，非编码RNA-蛋白质相互作用在其中发挥着关键的调节作用。细胞周期的正常进行依赖于一系列细胞周期蛋白（Cyclin）和细胞周期蛋白依赖性激酶（CDK）的有序激活和失活，而这些过程受到多种非编码RNA-蛋白质复合物的调控。例如，一些miRNA可以通过靶向调控Cyclin和CDK的表达，影响细胞周期的进程。miR-15和miR-16可以通过抑制CyclinD1的表达，使细胞周期阻滞在G1期，从而抑制细胞增殖；而miR-21则可以通过抑制PTEN的表达，激活PI3K/AKT信号通路，促进细胞周期从G1期向S期的转换，加速细胞增殖。此外，lncRNA也参与了细胞周期调控。例如，HOTAIRlncRNA可以通过与PRC2复合物结合，调控相关基因的表达，影响细胞周期的进程；UCA1lncRNA则可以通过与EZH2蛋白结合，抑制p21基因的表达，促进细胞周期的进展。这些非编码RNA与蛋白质的相互作用，共同维持了细胞周期的正常节律，确保细胞的有序增殖和分化。非编码RNA-蛋白质相互作用的失调与多种疾病的发生发展密切相关，在癌症、心血管疾病、神经系统疾病等复杂疾病的发病机制中扮演着重要角色。在癌症中，非编码RNA-蛋白质相互作用网络的异常往往导致细胞增殖、凋亡、侵袭和转移等过程的紊乱。许多癌基因和抑癌基因编码的蛋白质与非编码RNA之间存在着复杂的相互作用关系。例如，在乳腺癌中，miR-125b可以通过与HER2蛋白结合，抑制HER2信号通路的激活，从而抑制乳腺癌细胞的增殖和侵袭；而在肺癌中，lncRNAMALAT1可以通过与YBX1蛋白结合，促进肺癌细胞的增殖、迁移和侵袭。此外，circRNA也在癌症中发挥着重要作用。例如，circRNA-0001649可以通过吸附miR-125a-3p，上调其靶基因E2F3的表达，促进结直肠癌细胞的增殖和转移。在心血管疾病中，非编码RNA-蛋白质相互作用的异常也参与了疾病的发生发展过程。例如，在心肌梗死中，miR-1和miR-133可以通过与相关蛋白质相互作用，调节心肌细胞的凋亡、增殖和分化，影响心肌梗死后的心脏修复和重构；而lncRNAMIAT则可以通过与多种蛋白质结合，调控心肌细胞的电生理特性和心脏功能，与心律失常等心血管疾病的发生密切相关。在神经系统疾病中，非编码RNA-蛋白质相互作用的失调同样起着关键作用。例如，在阿尔茨海默病中，miR-128可以通过与APP蛋白结合，影响APP的代谢和Aβ的生成，参与阿尔茨海默病的发病过程；而lncRNABACE1-AS则可以通过与BACE1mRNA形成双链结构，稳定BACE1mRNA，促进Aβ的生成，加重阿尔茨海默病的病情。此外，circRNA在神经系统疾病中也具有重要作用。例如，circRNA-0001649可以通过吸附miR-7，上调其靶基因EGFR的表达，促进神经胶质瘤细胞的增殖和侵袭。非编码RNA与蛋白质之间的相互作用在基因表达调控、细胞周期调控、疾病发生发展等方面具有重要的生物学意义。深入研究这些相互作用关系，不仅有助于我们揭示生命活动的本质和疾病的发病机制，还为疾病的诊断、治疗和预防提供了新的靶点和策略，具有广阔的应用前景。2.4现有研究方法的局限性在非编码RNA-蛋白质相互作用关系的研究历程中，传统实验方法作为早期探索这一领域的主要手段，发挥了重要作用，但随着研究的深入，其固有的局限性逐渐凸显，为该领域的进一步发展带来了阻碍。传统实验方法中，免疫共沉淀（Co-IP）、RNA免疫沉淀（RIP）和交联免疫沉淀（CLIP）等技术是检测非编码RNA-蛋白质相互作用的经典手段。Co-IP技术通过抗体特异性识别并沉淀目标蛋白质，进而分离出与之相互作用的非编码RNA，但该方法存在一定的局限性。在实验过程中，由于细胞裂解后，非编码RNA与蛋白质之间的相互作用可能发生改变或解离，导致检测到的相互作用结果不能准确反映细胞内的真实情况。而且，该方法对于低丰度的非编码RNA-蛋白质相互作用的检测灵敏度较低，容易遗漏一些重要的相互作用信息。例如，某些非编码RNA在细胞内的表达水平较低，与蛋白质的结合亲和力也较弱，使用Co-IP技术很难有效地捕获到这些相互作用对。RIP技术利用针对目标蛋白质的抗体将RNA-蛋白复合物沉淀下来，从而对复合物中的RNA进行分析。然而，RIP技术也面临诸多挑战。在细胞裂解过程中，RNA与蛋白质的结合可能受到破坏，导致部分相互作用信息丢失。同时，RIP技术依赖于高质量的抗体，抗体的特异性和亲和力直接影响实验结果的准确性。如果抗体的特异性不佳，可能会沉淀出非特异性结合的RNA-蛋白复合物，增加实验结果的假阳性率。此外，RIP技术通常需要大量的细胞样本，对于一些难以获取大量样本的研究对象，如珍稀的临床样本或特定发育阶段的细胞，其应用受到了很大限制。CLIP技术通过紫外照射使RNA与蛋白质发生共价交联，增强了RNA与蛋白质的结合能力，在一定程度上提高了检测的准确性。但CLIP技术操作复杂，实验步骤繁琐，需要专业的技术人员和昂贵的实验设备。交联过程中可能会产生非特异性交联，导致背景信号增加，干扰对真实相互作用的判断。而且，CLIP技术对实验条件的要求非常严格，不同的实验条件可能会导致结果的差异较大，重复性较差。例如，紫外照射的强度和时间、交联剂的浓度等因素都会对实验结果产生显著影响，使得不同实验室之间的实验结果难以进行比较和验证。除了上述技术自身存在的问题外，传统实验方法还普遍存在成本高、耗时长的缺点。这些实验通常需要使用大量的试剂和耗材，如高质量的抗体、交联剂、核酸提取试剂等，导致实验成本居高不下。而且，从样本准备、实验操作到结果分析，整个实验周期较长，一般需要数周甚至数月的时间。这对于大规模研究非编码RNA-蛋白质相互作用关系来说，效率过低，难以满足快速发展的研究需求。传统实验方法的通量较低，一次实验只能检测少数几个非编码RNA与蛋白质的相互作用，无法同时对大量的非编码RNA和蛋白质对进行系统性研究。在生物体内，非编码RNA-蛋白质相互作用构成了一个庞大而复杂的网络，需要从全局角度进行分析和研究。传统实验方法的低通量特性限制了对这一复杂网络的全面解析，难以揭示非编码RNA-蛋白质相互作用的全貌和内在规律。综上所述，传统实验方法在研究非编码RNA-蛋白质相互作用时存在诸多局限性，迫切需要开发新的研究方法来弥补这些不足。机器学习方法作为一种高效、快速、可扩展性强的计算方法，为解决这些问题提供了新的途径，有望在非编码RNA-蛋白质相互作用关系的研究中发挥重要作用，推动该领域的深入发展。三、机器学习方法原理及应用3.1机器学习基础概念机器学习作为人工智能领域的核心技术之一，近年来在各个学科领域得到了广泛应用，为解决复杂问题提供了全新的思路和方法。其基本概念涵盖了从数据处理到模型构建、训练以及预测的一系列过程，通过让计算机自动从大量数据中学习模式和规律，实现对未知数据的准确预测和决策。根据数据的标注情况和学习目标的不同，机器学习主要分为监督学习、无监督学习和半监督学习三大类，每一类都具有独特的特点和适用场景，在生物信息学领域中也展现出了巨大的应用潜力。监督学习是机器学习中最为常见的一种类型，其核心特点是在训练过程中使用带有标签（标注）的数据。这些标签明确指示了输入数据所对应的输出结果，例如在疾病诊断中，输入的患者临床特征数据所对应的疾病类别标签。监督学习的目标是通过对这些有标签数据的学习，构建一个能够准确预测未知数据标签的模型。在预测非编码RNA-蛋白质相互作用关系时，可以将已知相互作用的非编码RNA和蛋白质对作为正样本，已知不相互作用的对作为负样本，每个样本都带有明确的“相互作用”或“不相互作用”标签。通过这些有标签的样本数据训练模型，如支持向量机（SVM）、决策树、随机森林等，使模型学习到非编码RNA和蛋白质特征与相互作用关系之间的映射规律。当输入新的非编码RNA和蛋白质特征时，模型能够根据学习到的规律预测它们是否相互作用。监督学习的优点在于其预测结果具有较高的准确性和可解释性，因为模型是基于明确的标签信息进行学习的。然而，它也存在一些局限性，例如需要大量的有标签数据进行训练，而获取高质量的标注数据往往需要耗费大量的时间和人力成本。在生物信息学中，实验测定非编码RNA-蛋白质相互作用关系需要复杂的实验技术和高昂的实验费用，这限制了有标签数据的规模。此外，如果训练数据的标注存在错误或偏差，可能会导致模型的预测性能下降。无监督学习则是在没有标签数据的情况下进行学习。其主要目的是发现数据中的潜在结构、模式或分组，例如聚类分析、主成分分析（PCA）等。在生物信息学中，无监督学习可以用于分析大量的非编码RNA和蛋白质数据，挖掘它们之间潜在的关联模式。可以对不同组织或细胞类型中的非编码RNA和蛋白质表达数据进行聚类分析，将具有相似表达模式的非编码RNA和蛋白质聚为一类，从而发现可能存在相互作用的非编码RNA-蛋白质对。无监督学习的优势在于它能够处理大量的未标注数据，发现数据中隐藏的信息和规律，不需要依赖于先验知识和标注信息。但是，由于没有明确的标签指导，无监督学习的结果解释相对困难，其发现的模式和结构可能需要进一步的实验验证和生物学解释。而且，不同的无监督学习算法对数据的假设和处理方式不同，可能会得到不同的结果，需要根据具体问题选择合适的算法和参数。半监督学习结合了监督学习和无监督学习的特点，使用少量的有标签数据和大量的无标签数据进行训练。在生物信息学研究中，获取大量有标签的非编码RNA-蛋白质相互作用数据较为困难，但未标注的非编码RNA和蛋白质数据却相对容易获得。半监督学习算法可以利用这些未标注数据中的信息，辅助少量有标签数据进行模型训练，从而提高模型的泛化能力和预测性能。半监督支持向量机（Semi-supervisedSupportVectorMachines）可以在有限的有标签样本和大量无标签样本上进行训练，通过对无标签样本的学习，进一步优化模型的决策边界，使其能够更好地适应未知数据。半监督学习在一定程度上缓解了监督学习对大量有标签数据的依赖，同时又避免了无监督学习结果难以解释的问题。然而，半监督学习算法的设计和实现较为复杂，需要平衡有标签数据和无标签数据的利用，并且对数据的分布和特征有一定的假设要求，如果假设不成立，可能会影响模型的性能。这三种机器学习类型在生物信息学研究非编码RNA-蛋白质相互作用关系中都具有重要的应用价值。监督学习能够利用已知的相互作用数据进行准确预测，但受限于标注数据的获取；无监督学习可以挖掘潜在的相互作用模式，但结果解释困难；半监督学习则在两者之间寻求平衡，利用少量标注数据和大量未标注数据提高模型性能。在实际应用中，需要根据具体的研究问题和数据特点，合理选择和运用不同的机器学习方法，以实现对非编码RNA-蛋白质相互作用关系的有效预测和深入理解。3.2适用于生物关系预测的机器学习算法3.2.1支持向量机（SVM）支持向量机（SupportVectorMachine，SVM）是一种基于统计学习理论的监督学习模型，在机器学习领域中占据着重要地位，尤其在处理小样本、非线性及高维模式识别问题时表现出色，因此在非编码RNA-蛋白质相互作用预测中得到了广泛应用。SVM的基本原理是在特征空间中寻找一个最优分类超平面，使得位于该超平面两侧的数据点尽可能远，这个距离被称为间隔（margin）。在二维平面中，超平面是一条直线；在三维空间中，超平面是一个平面；而在高维空间中，超平面则是一个维度比样本空间低一维的子空间。对于线性可分的数据集，假设存在一个超平面w^Tx+b=0（其中w是超平面的法向量，决定了超平面的方向；b是偏置项，决定了超平面的位置；x是样本特征向量），能够将两类样本正确分开。为了找到这个最优超平面，SVM通过最大化两类数据中离超平面最近的样本点（即支持向量）到超平面的距离之和，来确定超平面的参数w和b。这个过程可以转化为一个凸二次规划问题，通过求解该问题，可以得到最优的超平面参数，从而实现对数据的分类。在实际应用中，许多数据集往往是线性不可分的，即无法找到一个线性超平面将不同类别的数据完全分开。为了解决这个问题，SVM引入了核函数（KernelFunction）的概念。核函数的作用是将低维空间中的非线性问题映射到高维空间中，使得在高维空间中数据变得线性可分。常用的核函数有线性核（K(x,x')=x^Tx'）、多项式核（K(x,x')=(x^Tx'+c)^d，其中c是常数，d是多项式的次数）、径向基核（RBF，K(x,x')=\exp(-\gamma\|x-x'\|^2)，其中\gamma是核函数的参数）等。通过选择合适的核函数，SVM可以有效地处理非线性分类问题。在非编码RNA-蛋白质相互作用预测中，SVM被广泛应用于构建预测模型。研究人员会将非编码RNA和蛋白质的各种特征，如序列特征、结构特征、表达特征等，作为输入数据，将已知的相互作用关系作为标签，使用SVM进行模型训练。上海第二工业大学的研究团队提出一种基于伪特定位点记分矩阵（PsePSSM）和几何保留投影（GPP）降维算法的方法。首先采用PsePSSM算法表达蛋白质序列，编码后的蛋白质特征向量维数很高，接着采用GPP流形学习方法对其进行维数约简，约简后的特征向量输入支持向量机（SVM）分类器训练，用于预测未知的蛋白质与RNA是否相互作用。实验结果表明，该方法能有效解决传统实验方法的问题，为蛋白质与RNA相互作用的研究提供了新途径。通过SVM模型的学习，能够挖掘出这些特征与相互作用关系之间的潜在模式和规律，从而对未知的非编码RNA-蛋白质对进行相互作用预测。SVM在非编码RNA-蛋白质相互作用预测中具有较高的准确性和泛化能力，能够为生物学家提供有价值的参考信息，有助于加速对非编码RNA-蛋白质相互作用网络的解析。3.2.2随机森林（RandomForest）随机森林（RandomForest）是一种基于决策树的集成学习算法，通过构建多个决策树并将它们的结果进行整合，从而提高模型的预测性能和泛化能力。在生物信息学领域，尤其是在预测非编码RNA-蛋白质相互作用关系方面，随机森林展现出了独特的优势，为研究人员提供了一种强大的分析工具。随机森林的构建过程主要包括以下几个关键步骤：从原始数据集中使用自助采样法（bootstrapsampling）抽取多个样本子集。自助采样法是一种有放回的抽样方法，每次从原始数据集中随机抽取一个样本，重复抽取n次（n为原始数据集的样本数量），得到一个与原始数据集大小相同的样本子集。由于是有放回抽样，每个样本子集可能包含原始数据集中的重复样本，也可能有部分样本未被抽到。这样通过多次自助采样，可以得到多个不同的样本子集，增加了数据的多样性。对每个样本子集，构建一个决策树。在构建决策树的过程中，每个节点在进行分裂时，不是考虑所有的特征，而是随机选择一部分特征，然后从这部分特征中选择一个最优的特征进行分裂。通常，随机选取的特征数量等于总特征数量的平方根或对数值。这种特征随机选择的方式引入了随机性，使得每个决策树都有所不同，增加了模型的多样性，避免了过拟合的问题。重复以上步骤，直到生成指定数量的决策树，这些决策树共同构成了随机森林。在进行预测时，对于分类问题，随机森林通过投票的方式来确定最终的预测结果，即每个决策树对样本进行分类预测，然后统计所有决策树预测结果中出现次数最多的类别作为随机森林的预测类别；对于回归问题，则通过对所有决策树的预测结果进行平均来得到最终的预测值。随机森林之所以能够提高预测性能，主要是因为它通过集成多个决策树，有效地降低了模型的方差。每个决策树可能会对某些数据点做出错误的预测，但由于不同决策树之间具有一定的独立性和多样性，它们的错误预测往往不会集中在同一数据点上。当将多个决策树的预测结果进行整合时，这些错误预测会相互抵消，从而提高了整体的预测准确性。在生物关系预测中，随机森林在处理高维数据和复杂非线性关系方面表现出色。在预测非编码RNA-蛋白质相互作用时，生物数据通常包含大量的特征，如非编码RNA和蛋白质的序列特征、结构特征、表达特征等，这些特征之间可能存在复杂的非线性关系。随机森林能够自动处理这些高维特征，无需进行复杂的特征选择和降维操作，同时能够有效地捕捉到特征之间的非线性关系，从而提高预测的准确性。有研究利用随机森林算法对基因表达数据进行分析，以预测与特定疾病相关的基因。通过将基因表达数据作为特征，疾病状态作为标签，训练随机森林模型，成功地识别出了与疾病密切相关的基因。在非编码RNA-蛋白质相互作用预测中，也有研究采用随机森林算法，结合多种生物信息学特征，构建预测模型。通过对大量已知相互作用和非相互作用的非编码RNA-蛋白质对进行训练，随机森林模型能够学习到这些特征与相互作用关系之间的模式，从而对未知的非编码RNA-蛋白质对进行准确的相互作用预测。随机森林还可以衡量每个特征的重要性，帮助研究人员识别出对预测结果影响较大的关键特征，为进一步研究非编码RNA-蛋白质相互作用的机制提供了有价值的线索。3.2.3神经网络（NeuralNetwork）神经网络（NeuralNetwork），尤其是深度学习中的人工神经网络，近年来在各个领域取得了显著的进展，在生物信息学中处理复杂生物数据时展现出独特的优势，为预测非编码RNA-蛋白质相互作用关系提供了强有力的工具。神经网络的基本结构由输入层、隐藏层和输出层组成，各层之间通过神经元相互连接，信息在神经元之间传递并进行处理，从而实现对数据的特征提取和模式识别。输入层是神经网络与外部数据的接口，负责接收输入数据。在预测非编码RNA-蛋白质相互作用关系时，输入数据可以是经过编码的非编码RNA和蛋白质的序列信息、结构信息、表达信息等。这些信息被转化为数字特征向量后输入到神经网络中。例如，对于非编码RNA和蛋白质的序列信息，可以采用独热编码（One-HotEncoding）、k-mer编码等方式将其转化为数值向量，以便神经网络进行处理。隐藏层是神经网络的核心部分，位于输入层和输出层之间，可以包含一层或多层。隐藏层中的神经元通过权重与输入层和其他隐藏层的神经元相连，每个神经元接收来自前一层神经元的输入，并根据权重对输入进行加权求和，然后通过激活函数（如ReLU、Sigmoid、Tanh等）进行非线性变换，输出处理后的结果。激活函数的作用是为神经网络引入非线性特性，使其能够学习到复杂的非线性关系。在预测非编码RNA-蛋白质相互作用关系时，隐藏层通过对输入数据的层层处理和特征提取，能够自动学习到非编码RNA和蛋白质特征与相互作用关系之间的复杂模式和规律。随着隐藏层数量的增加，神经网络能够学习到更高级、更抽象的特征表示，从而提高对复杂数据的处理能力。例如，在深层神经网络中，前几个隐藏层可能学习到一些基础的序列特征和结构特征，而后面的隐藏层则能够学习到这些特征之间的组合关系和相互作用模式，进而推断出非编码RNA-蛋白质之间是否存在相互作用。输出层负责输出神经网络的预测结果。对于非编码RNA-蛋白质相互作用预测问题，输出层通常采用二分类的方式，即输出一个表示相互作用可能性的概率值。如果概率值大于某个阈值（如0.5），则预测为相互作用；否则，预测为不相互作用。在实际应用中，还可以根据具体需求对输出层进行调整，如采用多分类输出方式，以预测不同类型的非编码RNA-蛋白质相互作用。神经网络在处理复杂生物数据时具有多方面的优势。它能够自动学习数据中的复杂模式和特征表示，无需人工进行复杂的特征工程。生物数据通常具有高度的复杂性和非线性，传统的机器学习方法往往需要依赖人工设计和提取特征，而神经网络可以通过自身的学习机制，从原始数据中自动提取出有效的特征，从而提高模型的性能和泛化能力。神经网络具有强大的非线性建模能力，能够很好地捕捉非编码RNA和蛋白质特征与相互作用关系之间的非线性关系。在生物体内，非编码RNA-蛋白质相互作用受到多种因素的影响，这些因素之间的关系往往是非线性的，神经网络能够有效地处理这种非线性关系，从而实现更准确的预测。神经网络还具有较好的扩展性和适应性，可以通过调整网络结构和参数来适应不同的数据集和任务需求。通过增加隐藏层的数量、调整神经元的数量以及采用不同的训练算法和优化策略，可以优化神经网络的性能，使其更好地应用于非编码RNA-蛋白质相互作用预测等生物信息学任务。在非编码RNA-蛋白质相互作用预测的研究中，神经网络已被广泛应用并取得了一系列成果。2022年，早稻田大学的研究人员提出了BERT-RBP模型，通过对人类参考基因组上预训练的BERT架构进行调整，来预测RNA序列与RNA结合蛋白（RBPs）之间的相互作用。该模型在使用154个RBPs的eCLIP-seq数据时，展现出优于当时最先进预测模型的性能。2025年，华中师范大学物理研究所与美国华盛顿大学和达特茅斯学院合作，推出了ZHMolGraph模型。该模型集成了图神经网络和无监督大型语言模型来预测RNA-蛋白质相互作用，在两个基准数据集上验证时超越了当前的最佳方法。这些研究成果表明，神经网络在预测非编码RNA-蛋白质相互作用关系方面具有巨大的潜力，能够为生物学家提供更准确、更深入的预测结果，有助于推动非编码RNA-蛋白质相互作用机制的研究。3.2.4其他算法除了上述几种常用的机器学习算法外，还有一些其他算法也可用于生物关系预测，它们各自具有独特的特点和优势，在不同的应用场景中发挥着重要作用。朴素贝叶斯（NaiveBayes）算法是基于贝叶斯定理与特征条件独立假设的分类方法。其核心思想是在已知类别标签的条件下，假设各个特征之间相互独立，通过计算每个类别在给定特征下的条件概率，选择概率最大的类别作为预测结果。在预测非编码RNA-蛋白质相互作用时，朴素贝叶斯算法可以利用非编码RNA和蛋白质的特征信息，如序列特征、结构特征等，根据训练数据中不同特征与相互作用关系之间的概率分布，对未知样本进行分类预测。朴素贝叶斯算法的优点是算法简单、计算效率高，对小规模数据集表现出较好的性能，并且对缺失数据不太敏感。然而，由于其假设特征之间相互独立，在实际应用中，当特征之间存在较强的相关性时，朴素贝叶斯算法的性能可能会受到影响。逻辑回归（LogisticRegression）虽然名字中带有“回归”，但它实际上是一种用于解决二分类问题的线性分类模型。它通过构建一个逻辑函数（sigmoid函数），将线性回归的输出映射到0到1之间的概率值，从而实现对样本的分类。在非编码RNA-蛋白质相互作用预测中，逻辑回归可以将非编码RNA和蛋白质的特征作为输入变量，通过训练模型来学习特征与相互作用关系之间的线性关系，预测样本属于相互作用或不相互作用类别的概率。逻辑回归模型简单易懂，计算复杂度较低，训练速度快，并且具有较好的可解释性，可以通过系数来分析各个特征对预测结果的影响程度。然而，逻辑回归假设数据具有线性可分性，对于复杂的非线性关系的建模能力有限，在处理高维数据时可能需要进行特征选择或降维操作，以避免过拟合问题。决策树（DecisionTree）是一种基于树状结构的分类和回归模型。它通过对数据进行递归地划分，将数据分为不同的子集，直到满足一定的停止条件。在每个节点上，决策树根据某个特征的取值来决定数据的划分方向，每个分支代表一个特征值，每个叶子节点代表一个类别或预测值。在预测非编码RNA-蛋白质相互作用时，决策树可以根据非编码RNA和蛋白质的各种特征，如序列特征、结构特征、表达特征等，构建决策树模型，对未知样本进行分类预测。决策树的优点是模型直观、易于理解和解释，能够处理离散型和连续型数据，并且不需要对数据进行归一化等预处理操作。但是，决策树容易出现过拟合问题，对噪声数据比较敏感，泛化能力相对较弱。为了提高决策树的性能，可以采用剪枝等方法来防止过拟合，或者将多个决策树集成起来，如随机森林算法。这些其他机器学习算法在生物关系预测中各有优劣，研究人员可以根据具体的研究问题、数据特点和需求，选择合适的算法或算法组合来构建预测模型，以实现对非编码RNA-蛋白质相互作用关系的有效预测和分析。3.3机器学习在生物分子相互作用预测中的成功案例机器学习在生物分子相互作用预测领域取得了众多令人瞩目的成功案例，为深入理解生物分子间的复杂关系提供了有力支持。这些成功案例不仅展示了机器学习方法的有效性和潜力，也为预测非编码RNA-蛋白质相互作用关系提供了宝贵的经验和借鉴。在蛋白质-蛋白质相互作用预测方面，机器学习发挥了重要作用。2022年，BrunoE.Correia教授课题组在预印本平台bioRxiv上发表文章，提出利用蛋白质表面互作指纹（MaSIF，MolecularSurfaceInteractionFingerprinting）的方法，将蛋白质表面模式转化为数字阵列，用于识别倾向于形成蛋白-蛋白相互作用的区域，以及与特定靶点互补的表面。该方法可指导蛋白质结合位点的预测和蛋白结合剂的从头设计。随后，该团队又在Nature发文，表示利用此方法为SARS-CoV-2病毒的刺突蛋白和癌症免疫治疗的三个重要靶点（PD-1、PD-L1和CTLA-4）创造了新的蛋白质结合剂，这些结合剂都具有接近传统抗体的结合亲和力。这项研究的成功得益于对蛋白质表面特征的有效提取和机器学习算法的精确建模，通过挖掘蛋白质表面模式与相互作用之间的潜在关系，实现了对蛋白质-蛋白质相互作用的精准预测和新型结合剂的设计。其成功经验在于深入挖掘生物分子的结构特征，并将其转化为适合机器学习算法处理的数字特征，同时选择合适的机器学习模型进行训练和预测。这启示我们在预测非编码RNA-蛋白质相互作用时，也应注重对非编码RNA和蛋白质的结构、序列等特征的深入分析和有效提取，以提高预测模型的准确性。在DNA-蛋白质相互作用预测领域，湖南大学信息科学与工程学院彭绍亮教授课题组取得了显著进展。他们提出基于大型通用蛋白质语言模型和领域自适应预训练的DNA结合蛋白语言模型ESM-DBP，系统地研究了从蛋白质初级序列出发有效地预测DNA结合蛋白质和残基的问题。在多个基准测试集上，ESM-DBP与现有预测方法相比展现出优异的预测精度，还在只有少量相似同源序列的蛋白质序列上表现出较好的预测性能，远远超过通用蛋白质语言模型和其他预测方法。团队通过对预测模型的可解释性分析发现神经网络对DNA结合域的高度关注，从而导致ESM-DBP在DNA结合蛋白质预测任务上的高准确率，大大提升了蛋白质语言模型黑箱的可解释性。该研究的成功关键在于利用大规模蛋白质语言模型进行领域自适应预训练，充分挖掘蛋白质序列中的信息，并通过可解释性分析进一步优化模型。这为预测非编码RNA-蛋白质相互作用关系提供了重要的借鉴，即在构建预测模型时，可以借助大规模的生物分子数据和先进的机器学习模型，同时注重模型的可解释性研究，以更好地理解模型的决策过程和预测结果，提高模型的可靠性和实用性。这些成功案例表明，机器学习在生物分子相互作用预测中具有强大的能力和广阔的应用前景。通过深入挖掘生物分子的特征信息，选择合适的机器学习算法和模型，并注重模型的可解释性和性能优化，能够实现对生物分子相互作用的准确预测，为生命科学研究提供有力的支持。在预测非编码RNA-蛋白质相互作用关系时，可以借鉴这些成功经验，从多个方面入手，不断优化预测模型，提高预测的准确性和可靠性，推动非编码RNA-蛋白质相互作用领域的研究取得更大的进展。四、基于机器学习预测非编码RNA-蛋白质相互作用关系的方法构建4.1数据收集与预处理数据收集与预处理是基于机器学习预测非编码RNA-蛋白质相互作用关系的基础步骤，其质量直接影响后续模型的训练和预测效果。本研究通过多种途径收集了非编码RNA和蛋白质的序列数据、结构数据以及表达数据，并运用一系列严格的数据清洗和预处理方法，确保数据的准确性、完整性和可用性。在序列数据收集方面，非编码RNA序列主要来源于权威的数据库，如miRBase、LncRNAdb、circBase等。这些数据库收录了大量经过实验验证和注释的非编码RNA序列信息，涵盖了多种物种和组织类型。对于微小RNA（miRNA），从miRBase数据库中获取其成熟序列和前体序列，该数据库定期更新，保证了数据的时效性和可靠性。蛋白质序列则主要从UniProt数据库中获取，UniProt是全球最大的蛋白质序列和功能信息数据库之一，包含了丰富的蛋白质序列注释信息，如蛋白质的功能描述、结构域信息、翻译后修饰位点等，为后续的特征提取和分析提供了重要依据。在收集过程中，充分考虑了数据的多样性和代表性，确保涵盖不同物种（如人类、小鼠、大鼠等）、不同组织（如肝脏、心脏、大脑等）以及不同生理病理状态下的非编码RNA和蛋白质序列。对于疾病相关的研究，特别收集了疾病样本和正常对照样本中的非编码RNA和蛋白质序列，以便分析疾病状态下非编码RNA-蛋白质相互作用关系的变化。结构数据对于理解非编码RNA和蛋白质的功能及相互作用机制至关重要。非编码RNA的二级结构数据通过RNAfold等软件进行预测获得。RNAfold基于最小自由能原理，利用动态规划算法预测RNA的二级结构，能够准确地预测出RNA分子中的茎环结构、发卡结构等。蛋白质的三级结构数据主要从蛋白质数据库（PDB）中获取，PDB收录了大量通过X射线晶体学、核磁共振等实验技术解析得到的蛋白质三维结构信息。对于一些尚未有实验结构数据的蛋白质，利用同源建模方法，如SWISS-MODEL、MODELLER等，基于已知结构的同源蛋白质构建其三维结构模型。在获取结构数据时，对数据的质量进行了严格筛选。对于RNA二级结构预测结果，评估其预测的可靠性，如通过计算预测结构的自由能、与已知结构的相似性等指标进行判断。对于蛋白质结构数据，优先选择分辨率高、结构完整性好的结构模型，以确保后续分析的准确性。表达数据能够反映非编码RNA和蛋白质在不同组织、不同发育阶段以及不同生理病理条件下的表达水平变化，为研究非编码RNA-蛋白质相互作用关系提供了重要线索。表达数据主要通过高通量测序技术（如RNA-Seq、ChIP-Seq等）和微阵列技术获得。对于非编码RNA的表达数据，利用RNA-Seq技术对不同样本进行测序，然后通过数据处理和分析，计算出每个非编码RNA在不同样本中的表达量，常用的分析工具包括TopHat、Cufflinks等。蛋白质的表达数据则可以通过蛋白质组学技术，如液相色谱-质谱联用（LC-MS/MS）、抗体芯片等进行检测。通过这些技术，可以获得蛋白质在不同样本中的相对表达水平或绝对表达量。在数据收集过程中，为了保证数据的一致性和可比性，严格控制实验条件和样本处理流程。对于RNA-Seq实验，确保样本的采集、RNA提取、文库构建和测序等步骤按照标准化的操作规程进行；对于蛋白质组学实验，同样严格控制样本的处理、分离和检测过程，以减少实验误差对数据质量的影响。数据清洗是数据预处理的关键环节，主要目的是去除数据中的噪声、错误和异常值，提高数据的质量。对于序列数据，检查序列的完整性，去除含有缺失碱基或氨基酸的序列。对于结构数据，检查结构的合理性，如蛋白质结构中的键长、键角是否符合化学原理，RNA二级结构中的碱基配对是否合理等，去除结构异常的数据。对于表达数据，通过数据分布分析，识别并去除表达量异常高或异常低的样本，这些异常值可能是由于实验误差或样本污染等原因导致的。为了使不同类型的数据具有可比性，对数据进行标准化处理。对于表达数据，采用Z-score标准化方法，将每个样本的表达量转化为标准正态分布，使得数据的均值为0，标准差为1。对于数值型的特征数据，如蛋白质的理化性质参数等，也采用类似的标准化方法进行处理。对于分类数据，如非编码RNA的类型（miRNA、lncRNA、circRNA等）和蛋白质的功能分类等，采用独热编码（One-HotEncoding）方法进行编码，将其转化为数值型数据，以便机器学习模型能够处理。数据收集与预处理是构建基于机器学习的非编码RNA-蛋白质相互作用关系预测模型的重要前提。通过广泛收集高质量的数据，并运用科学的清洗和预处理方法，为后续的模型训练和预测提供了可靠的数据基础，有助于提高模型的准确性和可靠性。4.2特征提取与选择4.2.1序列特征提取序列特征提取是预测非编码RNA-蛋白质相互作用关系的关键环节，通过从非编码RNA和蛋白质序列中挖掘有效的特征信息，能够为后续的机器学习模型提供丰富的数据支持，从而提高模型的预测准确性。常见的序列特征提取方法包括k-mer频率、序列相似性等，这些方法从不同角度对序列信息进行量化和分析，为揭示非编码RNA与蛋白质之间的相互作用机制提供了重要线索。k-mer频率是一种常用的序列特征提取方法，它将序列划分为固定长度为k的子序列（k-mer），并统计每个k-mer在序列中出现的频率。对于非编码RNA序列，k通常取值在2-6之间；对于蛋白质序列，k的取值范围一般为2-3。以非编码RNA序列“AGCUAGCU”为例，当k=2时，对应的2-mer包括“AG”“GC”“CU”“UA”“AG”“GC”“CU”，统计这些2-mer的出现频率，即可得到该非编码RNA序列的一个特征向量。在蛋白质序列中，2-mer也被称为双肽，如蛋白质序列“MALWMRLLPLLALLALWGPDPAAAFVNQHLCGSHLVEALYLVCGERGFFYTPKT”，其中的双肽“MA”“AL”“LW”等的出现频率构成了蛋白质序列的k-mer特征。k-mer频率能够反映序列中局部片段的组成和分布情况，不同的非编码RNA和蛋白质可能具有独特的k-mer频率模式，这些模式与它们的结构和功能密切相关，从而有助于预测它们之间的相互作用关系。序列相似性也是一种重要的序列特征。通过计算非编码RNA或蛋白质序列之间的相似性，可以判断它们在进化上的亲缘关系以及功能上的相关性。常用的序列相似性计算方法有BLAST（BasicLocalAlignmentSearchTool）和Smith-Waterman算法。BLAST是一种基于启发式算法的快速序列比对工具，它能够在数据库中快速搜索与查询序列相似的序列，并给出相似性得分和比对结果。在预测非编码RNA-蛋白质相互作用时，可以利用BLAST将待预测的非编码RNA或蛋白质序列与已知相互作用的序列数据库进行比对，通过分析比对结果中的相似性得分、比对长度、一致性百分比等指标，来推断待预测序列与已知相互作用序列之间的相似程度，进而预测它们之间是否存在相互作用。Smith-Waterman算法则是一种基于动态规划的全局序列比对算法，它能够找出两条序列之间的最优局部比对，得到的比对结果更加准确，但计算复杂度较高，适用于对准确性要求较高且序列长度较短的情况。例如，在研究某种新发现的非编码RNA与蛋白质的相互作用时，通过Smith-Waterman算法将其与已知相互作用的非编码RNA序列进行比对，若发现具有较高的相似性，则提示该新非编码RNA可能与相应的蛋白质存在相互作用。除了k-mer频率和序列相似性外，还可以提取其他一些序列特征。可以计算序列的碱基组成或氨基酸组成，即统计序列中各种碱基（A、T、C、G）或氨基酸的比例，这些组成信息能够反映序列的基本特征，不同类型的非编码RNA和蛋白质往往具有不同的碱基或氨基酸组成模式，对预测相互作用关系具有一定的参考价值。还可以考虑序列的保守性，通过多序列比对分析，确定序列中保守区域和变异区域，保守区域可能在非编码RNA-蛋白质相互作用中发挥重要作用，因为它们往往与功能密切相关，在进化过程中受到选择压力的影响而保持相对稳定。在预测过程中，可以将保守性特征与其他特征相结合，提高预测模型的性能。通过k-mer频率、序列相似性等方法提取非编码RNA和蛋白质的序列特征，能够从多个维度挖掘序列信息，为预测非编码RNA-蛋白质相互作用关系提供丰富的数据基础。这些序列特征不仅能够反映序列的基本组成和结构特点，还能在一定程度上揭示它们在进化和功能上的联系，对于深入理解非编码RNA与蛋白质之间的相互作用机制具有重要意义。在实际应用中，应根据具体的研究需求和数据特点，合理选择和组合序列特征提取方法，以提高预测模型的准确性和可靠性。4.2.2结构特征提取结构特征提取在预测非编码RNA-蛋白质相互作用关系中起着至关重要的作用，因为非编码RNA和蛋白质的结构与其功能密切相关，直接影响它们之间的相互作用方式和亲和力。通过获取和分析非编码RNA和蛋白质的二级、三级结构特征，如二级结构元件、结构域等，可以深入了解它们的结构特点和相互作用机制，为机器学习模型提供更丰富、更准确的信息，从而提高预测的准确性和可靠性。非编码RNA的二级结构由碱基之间的互补配对形成，主要包括茎环结构、发卡结构、凸环结构和内部环结构等。这些二级结构元件在非编码RNA与蛋白质的相互作用中发挥着重要作用。茎环结构中的环区通常是与蛋白质结合的关键部位，因为环区的核苷酸序列相对灵活，更容易与蛋白质表面的氨基酸残基形成特异性的相互作用。发卡结构的稳定性和形状也会影响非编码RNA与蛋白质的结合能力。为了提取非编码RNA的二级结构特征，可以使用RNAfold等软件进行预测。RNAfold基于最小自由能原理，通过动态规划算法计算RNA序列形成不同二级结构的自由能，从而预测出最稳定的二级结构。在得到二级结构预测结果后，可以采用多种方法进行特征提取。可以将二级结构表示为点括号表示法，如“((((....))))”，其中括号表示配对的碱基，点表示未配对的碱基，然后通过统计不同类型括号和点的数量、位置以及它们之间的组合关系，来提取二级结构特征。还可以将二级结构转化为图形表示，如将茎环结构表示为节点和边的图，通过图论方法计算图的特征，如节点度、最短路径等，作为非编码RNA的二级结构特征。蛋白质的二级结构主要包括α-螺旋、β-折叠、β-转角和无规卷曲等。这些二级结构元件通过氢键等相互作用维持其稳定性，并构成了蛋白质的基本结构框架。α-螺旋和β-折叠是蛋白质中最常见的二级结构，它们的分布和排列方式对蛋白质的三维结构和功能具有重要影响。在蛋白质与非编码RNA相互作用时，α-螺旋和β-折叠可以通过其表面的氨基酸残基与非编码RNA形成相互作用。为了提取蛋白质的二级结构特征，可以利用DSSP（DefineSecondaryStructureofProteins）等软件进行分