我们描述了一种基于序列多样化的研究方法,用于估算蛋白质-蛋白质相互作用(PPIs)中多特异性结合位点的氨基酸偏好性。该策略通过生成并筛选数千种潜在的肽配体来实现 计算机模拟从而克服了现有实验方法的一些局限性。
需要JoVE订阅才能观看此内容。 请登录或开始免费试用
我们描述了一种基于序列多样化的研究方法,用于估算蛋白质-蛋白质相互作用(PPIs)中多特异性结合位点的氨基酸偏好性。该策略通过生成并筛选数千种潜在的肽配体来实现 计算机模拟从而克服了现有实验方法的一些局限性。
许多蛋白质-蛋白质相互作用涉及短肽段与肽结合结构域的结合。通常,此类相互作用需要识别具有可变保守性的线性基序。在同一配体中,高度保守区域与较可变区域的组合常常导致结合的多特异性,这是酶和细胞信号传导蛋白的常见特性。表征肽结合结构域的氨基酸偏好对于设计蛋白质-蛋白质相互作用(PPIs)的调控分子至关重要。计算方法是对通常成本高昂且操作繁琐的实验技术的一种高效替代方案,能够设计出潜在的调控分子,并可在后续实验中进行验证。本文介绍了一种利用Rosetta分子建模软件包中的Pepspec应用程序来预测肽结合结构域氨基酸偏好的方法。当受体蛋白的结构和肽配体的性质均已知或可推断时,该方法尤为适用。该方法首先从配体中一个特征明确的锚定残基开始,通过随机添加氨基酸残基进行延伸。随后通过柔性骨架肽段对接评估所生成肽段的结合亲和力,以筛选出预测结合评分最优的肽段。这些肽段可用于计算氨基酸偏好,并可选择性地构建位置权重矩阵(PWM),供后续研究使用。为展示该方法的应用,我们研究了人干扰素调节因子5(IRF5)亚基之间的相互作用。该相互作用此前已知具有多特异性,但整体上由一个称为pLxIS的短保守基序引导。所预测的氨基酸偏好与先前关于IRF5结合表面的知识一致。可被磷酸化的丝氨酸残基所在位置表现出天冬氨酸和谷氨酸的高频出现,这很可能是因为它们带负电荷的侧链与磷酸化丝氨酸相似。
两种蛋白质之间的相互作用通常涉及短肽段氨基酸与肽结合结构域的结合,类似于蛋白质-肽界面。参与此类蛋白质-蛋白质相互作用(PPI)的受体蛋白通常能够识别一组重叠但不同的配体序列,这一特性称为多特异性1,2。多特异性识别是许多细胞蛋白的特征,但在酶和细胞信号传导蛋白中尤为显著3。与多特异性结合位点相互作用的蛋白质,其序列中通常包含保守程度较高和较低区域的组合4,5,6。在此情况下,更为保守的序列基序参与严格的分子相互作用;相反,可变性较高的序列则与受体结合位点中具有一定宽容性的表面发生相互作用。通常,这些保守性较低但仍具有功能相关性的片段为缺乏明确二级结构模式的环状区域,或具有更动态构象的区域,如内在无序蛋白所典型的结构7。
鉴定结合位点的潜在肽配体通常是设计能够干扰相应蛋白质相互作用(PPIs)的介导分子的第一步8。然而,在多特异性结合位点的配体中,大多数序列位置上通常难以找到单一出现频率最高的氨基酸残基。相反,这些位点可能根据氨基酸的化学性质,对某一类特定氨基酸具有偏好性,例如带负电荷的酸性氨基酸(如天冬氨酸或谷氨酸)、具有大体积芳香环的氨基酸(如苯丙氨酸),或更疏水的残基(如脂肪族氨基酸丙氨酸、缬氨酸、亮氨酸或异亮氨酸)3。已有多种实验方法可揭示蛋白质结合位点的氨基酸偏好性,包括定向进化9、多密码子扫描突变10和深度突变扫描11。这些方法均基于序列多样化策略,即对原始配体引入突变,并进一步分析其对受体蛋白功能的影响(详见Bratulic和Badran12的综述)。然而,这些方法通常需要检测大规模的序列文库,因而更加繁琐、成本较高且耗时较长。
推断多特异性结合位点氨基酸偏好的计算方法有望克服湿实验方法的局限性。其中,计算机序列多样化方法通过评估配体序列中大量氨基酸替换的能量效应,来表征蛋白质-蛋白质相互作用(PPI)的结构可塑性13。该方法首先以结合于受体结合位点的肽配体结构或模型为基础,随后对配体序列引入突变。接着利用统计学和能量评分函数评估这些突变对稳定性和结合亲和力的影响。评估阶段所得的一组高分配体序列可用于计算氨基酸偏好。该策略能够以高效的方式处理大量配体序列,因此相较于湿实验方法通常所能处理的有限序列数量,可提供更全面且一致的氨基酸偏好推断结果。
Rosetta 分子建模套件中的 Pepspec 应用程序14 是一种在肽设计模式中将序列多样化作为关键步骤的工具。该应用程序需要一个受体蛋白的结构或模型,其中结合有一个长度至少为单个氨基酸残基的肽,该肽将作为后续步骤的锚定点。随后,结合肽的序列会被延伸(如有必要)并进行多样化,以生成大量假定的肽配体。接着通过柔性骨架肽对接来评估这些肽的结合亲和力,从而筛选出预测结合评分最优的肽。尽管该应用程序的主要输出结果是设计阶段末期筛选出的最佳肽候选物,但在该阶段接受的大量肽序列也可用于计算靶标结合位点的氨基酸偏好性。氨基酸偏好性以每个配体序列位置上各氨基酸残基出现的频率来表示,可呈现为位置权重矩阵(PWM)或更直观的序列标志图(sequence logo)。
本文描述了一种用于估算参与蛋白质-蛋白质相互作用(PPI)的受体蛋白结合表面氨基酸偏好的实验方案。该方案适用于已知蛋白配体的线性片段与受体蛋白结合的情形,因此可将该相互作用建模为蛋白-肽界面。在此情形下,配体中的保守基序通常与受体结合位点内的特定口袋相互作用,尽管参与PPI的整个配体片段可能包含保守性较低的区域。协议主要步骤的流程图见图1。该方案从蛋白-蛋白复合物的三维结构出发,将配体蛋白进一步缩减为最有可能发生强相互作用的片段,而保持受体蛋白结构不变。最有可能发生强相互作用的片段通过BUDE Alanine Scan服务器15推断得出,该服务器通过计算丙氨酸扫描突变分析,识别两个相互作用蛋白之间的热点残基。在此方法中,配体中的残基被逐一替换为丙氨酸,随后根据复合物自由能或稳定性的估计变化值(ΔΔG)来推断相应残基在目标PPI中的重要性。确定最有可能发生强相互作用的片段后,将其与受体蛋白形成的复合物作为基础结构提交至Pepspec,以进行序列多样化分析。

图1:本研究提出的实验方案主要步骤概览。 图中标注的数字与方案部分的步骤编号相对应。图中所用示例为文中描述的蛋白质-蛋白质复合物。在此复合物中,被视为受体的蛋白质链以粉红色显示,而被视为配体的链以浅蓝色显示,其中预测的最适相互作用片段以红色突出显示。请点击此处查看该图的放大版本。
本方案的一个局限性在于需要已解析的蛋白质-肽段相互作用界面结构。该方案也可从目标蛋白质-肽段相互作用界面的模型开始,但本文未描述具体的建模步骤。此外,尽管该方案可在任何操作系统的个人计算机上运行,但涉及 Rosetta 应用程序的步骤需要 Linux 环境。由于 Pepspec 通常需要执行大量迭代,因此强烈建议在序列多样化步骤中使用计算机集群。
通过估算人干扰素调节因子(IRF)家族成员IRF5结合表面的氨基酸偏好性,展示本建议方案的应用。我们选择该蛋白作为示例,是因为在其激活过程中,两个亚基结合形成二聚体,其结构已被充分表征16。在IRF二聚体中,结合可被建模为一种蛋白-肽界面,其中一个亚基提供结合表面,另一个亚基则通过一段包含短保守基序pLxIS的区域进行相互作用17,18。此外,与IRF亚基的结合具有多特异性,因此它们可以形成同源二聚体、异源二聚体,以及与其他已知为共激活因子的细胞蛋白形成的复合物18。
访问受限。请登录或开始试用以查看此内容。
1. 蛋白质-肽段相互作用界面的初始准备

图 2:本研究中用作代表性示例的蛋白质数据库(PDB)结构页面。(A)用于输入目标结构 PDB 访问编号的搜索框。(B)以多种格式下载该结构的菜单。(C)当结构以非对称单元形式保存时,下载生物组装体的选项(详见步骤 1.1.2)。请点击此处查看该图的放大版本。

图3:在BUDE丙氨酸扫描服务器中选择受体和配体。(A)蛋白-蛋白复合物的图形表示。(B)用于在服务器中输入任务名称的文本框。(C)用于交互式选择将作为受体和配体的链的面板(详见步骤1.2)。请点击此处查看该图的放大版本。

图4:BUDE丙氨酸扫描服务器的结果页面。 配体序列中潜在的最佳相互作用区段以红色框标出。在左侧面板中,预测能量贡献较高的残基(Leu433)以绿色高亮显示。 请点击此处查看该图的放大版本。
2. 序列多样化
注意:在以下步骤中,rosetta_main 指的是 Rosetta 的主安装目录,通常位于 /opt/rosetta_src__bundle/main/,其中 表示已安装的 Rosetta 版本。此外,假设 Rosetta 应用程序可在系统范围内访问;如果情况并非如此,则必须提供可执行文件的完整路径。当从源代码编译时,这些可执行文件位于 /rosetta_main/source/bin/ 目录中。


3. 氨基酸偏好性的估算

访问受限。请登录或开始试用以查看此内容。
本文介绍了一种用于预测干扰素调节因子5(IRF5)结合表面氨基酸偏好的实验方案。IRF5属于干扰素调节因子家族,该家族是一类被称为人干扰素调节因子的转录因子。这些蛋白质是先天性和适应性免疫应答的调控因子,参与多种免疫细胞的分化与活化过程。IRF亚基具有高度可塑且多特异性的结合表面,能够形成同源二聚体、异源二聚体以及与其他细胞蛋白的复合物17,18。二聚化被认为是这些因子活化的第一步,在大多数家族成员中,该过程由多个丝氨酸/苏氨酸残基的磷酸化所触发18。在二聚化过程中,每个单体通过其序列C端区域附近一个高度保守的基序pLxIS与另一个单体的结合表面相互作用。pLxIS这一缩写部分代表了结合表面的氨基酸偏好,即依次识别一个极性氨基酸(“p”),随后是两个富含亮氨酸(“L”)和异亮氨酸(“I”)的位置,这两个位置之间由任意氨基酸占据的一个位置(“x”)隔开,最后是一个可磷酸化的丝氨酸残基(本例中为Ser436)。包括pLxIS基序在内的多个丝氨酸残基的磷酸化,可促进一个单体C端片...
访问受限。请登录或开始试用以查看此内容。
本文介绍了一种基于计算机序列多样化方法来估算潜在多特异性结合位点氨基酸偏好的实验方案。目前已有少数计算工具被开发用于预测蛋白质-肽相互作用界面的氨基酸偏好14,25,26。这些工具具有预测性质,但在用于预测的计算算法以及为提高准确性所采用的校正方法上有所不同。在本研究中,我们使用了Rosetta分子建模套件中的Pepspec应用程序14。尽管该程序主要面向肽的设计,但它实现了一种序列多样化算法,可用于预测氨基酸偏好。据我们所知,该工具是目前唯一提供内置脚本、可直接从序列多样化结果中计算PWM(位置权重矩阵)的可用工具。需要特别指出的是,本方案聚焦于蛋白质-蛋白质相互作用(PPIs),因此初始结构应为两个蛋白质亚基形成的复合物。在实际进行序列多样化步骤之前,需将被视为配体的蛋白质截取为预期与受体蛋白发生相互作用的片段,并进一步将其作为肽处理。然而,该方案也可应用于蛋白质-肽复合物,在此类情况下,步骤1.1–...
访问受限。请登录或开始试用以查看此内容。
作者无任何利益冲突需要披露。
感谢巴拿马国家研究系统(SNI)(资助编号 SNI-043-2023 和 SNI-170-2021)、巴拿马国家科学技术与创新秘书处(SENACYT)以及人力资源培训与利用研究所(IFARHU)提供的经费支持。作者谨此感谢 Miguel Rodríguez 博士对稿件的仔细审阅。
访问受限。请登录或开始试用以查看此内容。
| 姓名 | 公司 | 目录编号 | 评论 |
|---|---|---|---|
| BUDE丙氨酸扫描服务器 | 爱丁堡大学 | https://pragmaticproteindesign.bio.ed.ac.uk/balas/ | doi: 10.1021/acschembio.9b00560 |
| Rosetta建模软件 | Rosetta Commons | https://www.rosettacommons.org/software | doi: 10.1002/prot.22851 |
| UCSF Chimera | 加州大学旧金山分校 | https://www.cgl.ucsf.edu/chimera/ | doi: 10.1002/jcc.20084 |
访问受限。请登录或开始试用以查看此内容。