Quantification of Membrane Protein Conformational Free Energy from Mutations and a Single Atom
基于突变和单个原子的膜蛋白构象自由能定量分析
来源：J. Am. Chem. Soc. 2025, 147, 34316−34326
1.摘要
次级主动转运蛋白通过在向内和向外两种构象间切换实现底物跨膜转运，但其构象间的相对自由能差异仍缺乏精确的实验测量。本研究利用¹⁹F核磁共振（NMR）光谱技术，定量测定了大肠杆菌多药转运蛋白EmrE两种构象的自由能差。通过构建仅保留一个色氨酸残基（Trp63）的EmrE突变体，实现了特异性的氟标记，并制备了仅一个单体携带标记或突变的最小异二聚体。研究发现，单个保守氨基酸突变（L51I）可使构象平衡偏移高达1.5 kcal/mol，而单个氢原子被氟原子取代也能导致0.8 kcal/mol的自由能变化。分子动力学模拟表明，氟原子通过差异氢键作用调控构象偏好。这些结果提供了定量证据，证明微小的分子扰动即可显著影响转运蛋白的自由能景观，这种可塑性可能在进化中帮助微调转运机制，为理解膜蛋白的进化驱动力提供了新视角。
2.关键词（中文）
膜蛋白构象自由能、次级主动转运蛋白、EmrE、¹⁹F核磁共振、构象平衡、蛋白质进化
3.研究目的
一是建立基于¹⁹F NMR的膜蛋白构象自由能高精度定量测量方法，填补传统技术难以精确测定构象间微小自由能差的空白；二是定量评估单个保守氨基酸突变和单个原子取代对EmrE内向-外向构象平衡的影响程度；三是从分子层面揭示微小扰动调控膜蛋白构象自由能的作用机制；四是基于构象自由能的可塑性，提出膜转运蛋白从同源寡聚体向单链融合蛋白进化的新假说，解释其进化优势。
4.研究思路
首先构建EmrE的单色氨酸突变体（W31F/W45F/W76F，仅保留底物结合口袋中的Trp63），用于特异性¹⁹F标记，同时验证该突变体保留野生型转运功能；分别用5-氟色氨酸（5FW）和6-氟色氨酸（6FW）标记Trp63，制备不同标记的EmrE蛋白；通过混合标记蛋白与未标记蛋白（或突变体蛋白）形成异二聚体，利用¹⁹F NMR光谱区分不同构象的单体信号；通过光谱去卷积计算两种构象的种群比例，结合热力学公式推导自由能差；引入L51I和E14Q两个保守突变，定量测定突变对构象平衡的影响，并校正氟原子本身的构象效应；结合全原子分子动力学（MD）模拟，分析氟原子与周围残基的相互作用差异，阐明其构象偏好的分子机制；最后整合所有实验结果，提出膜转运蛋白进化的自由能驱动假说。
5.研究亮点
一是实现了膜蛋白构象自由能的超高精度定量测量，检测限低至0.2 kcal/mol，首次捕捉到单个原子取代引起的构象自由能变化；二是发现单个保守突变和单个氟原子取代即可显著改变EmrE的构象平衡，证明膜蛋白的自由能景观具有高度可塑性，挑战了“构象平衡由整体结构决定”的传统认知；三是通过MD模拟阐明了氟原子通过差异氢键作用调控构象偏好的分子机制，5FW在单体A中形成更多分子间氢键，6FW在单体B中形成更多分子内氢键；四是从自由能角度提出了单链融合转运蛋白的进化优势假说：单链蛋白的单点突变可直接调控构象平衡和转运功能，而同源二聚体的双点突变会产生抵消效应，限制了其进化潜力；五是建立的最小异二聚体¹⁹F NMR方法具有普适性，可推广到其他同源寡聚体膜蛋白的构象动力学研究。
6.可延伸方向
一是将该定量方法应用于MFS、ABC等其他家族的次级主动转运蛋白，验证构象自由能可塑性的普遍性；二是系统研究底物结合、质子梯度、膜电位等生理条件对EmrE构象自由能的影响，构建完整的转运循环能量景观模型；三是结合冷冻电镜和单分子荧光技术，在原子水平和单分子层面验证构象自由能的测量结果，实现多尺度的构象动力学表征；四是通过定向进化筛选能精准调控构象平衡的突变体，开发基于转运蛋白的生物传感器和药物递送系统；五是利用祖先序列重建技术，验证单链融合转运蛋白的进化优势假说，追溯膜转运蛋白的进化历程；六是基于构象自由能的定量数据，指导设计能特异性锁定转运蛋白某一构象的新型抗生素和抗癌药物，克服多药耐药性。
7.测量的数据及研究意义
测定了EmrE的三维结构（PDB ID:8UWU）及Trp63在底物结合口袋中的位置，数据来自图1a、图1b。意义：明确了Trp63作为构象探针的合理性，其位于底物结合口袋中心，靠近质子结合关键残基Glu14，能灵敏反映构象变化。
测定了野生型EmrE、EmrE^W63突变体、E14D失活突变体和空载体对照对溴化乙锭、原黄素和甲基紫精的半抑制浓度（IC₅₀），数据来自图1c、图S1。意义：验证了EmrE^W63突变体保留了野生型的多药外排功能，证明其可作为可靠的研究模型，同时建立了功能丧失的基准线。
测定了5FW标记的野生型EmrE和EmrE^W63的一维¹⁹F NMR光谱，数据来自图1d。意义：确认了Trp63的¹⁹F信号能清晰区分EmrE的两种构象，且EmrE^W63的光谱无重叠干扰，适合后续定量分析。
测定了不同比例5FW/6FW标记EmrE与未标记EmrE混合后的一维¹⁹F NMR光谱，数据来自图2a、图2d。意义：首次发现单个氟原子取代即可改变构象平衡，5FW偏向单体A，6FW偏向单体B。
测定了部分标记EmrE的二维¹⁹F/¹⁹F ZZ-交换光谱，数据来自图2b、图2e。意义：直接证明了两种构象之间存在动态交换，验证了光谱峰的归属。
总结了5FW和6FW取代对构象自由能的影响，数据来自图2c、图2f。意义：定量给出了单个氟原子的构象效应，5FW为0.8±0.2 kcal/mol，6FW为-0.2±0.1 kcal/mol。
测定了MD模拟中5FW和6FW在不同单体中的氢键接触频率，数据来自图3a、图3c。意义：从分子层面揭示了氟原子影响构象偏好的机制，5FW在单体A中形成更多分子间氢键，6FW在单体B中形成更多分子内氢键。
展示了MD模拟中5FW和6FW与周围残基的相互作用代表性构象，数据来自图3b、图3d。意义：直观呈现了氢键相互作用的细节，与NMR实验结果相互印证。
测定了L51I突变体异二聚体的一维¹⁹F NMR光谱及去卷积结果，数据来自图4a、图4b。意义：定量测定了单个保守突变L51I的构象效应，表观自由能差达2.3 kcal/mol。
测定了5FW-L51I与6FW混合异二聚体的一维和二维¹⁹F NMR光谱，数据来自图4c、图4d。意义：证明了突变和氟原子取代的效应具有加和性，可通过组合实现更大的构象偏移。
测定了E14Q突变体异二聚体的一维和二维¹⁹F NMR光谱及去卷积结果，数据来自图5a、图5b、图5c、图5d、图5e。意义：定量测定了E14Q突变的构象效应，揭示了Glu14质子化状态对构象平衡的重要性。
总结了所有样品的表观自由能和校正后的突变/氟原子自由能贡献，数据来自表1。意义：系统汇总了所有定量测量结果，清晰展示了不同扰动对构象自由能的影响程度。
展示了单链多肽基因与同源寡聚体基因进化分化的假说示意图，数据来自图6。意义：基于实验结果提出了膜转运蛋白进化的新模型，解释了自然界中单链融合转运蛋白占主导的现象。
8.结论
本研究建立了基于¹⁹F NMR的高精度膜蛋白构象自由能定量方法，成功测定了大肠杆菌多药转运蛋白EmrE内向和外向构象的自由能差。研究发现，膜蛋白的自由能景观具有高度可塑性：单个保守氨基酸突变（L51I）可使构象平衡偏移高达1.5 kcal/mol，而单个氢原子被氟原子取代也能导致0.8 kcal/mol的显著变化。分子动力学模拟表明，氟原子通过与周围残基形成差异氢键作用，选择性稳定某一构象。基于这些发现，提出了单链融合转运蛋白的进化优势假说：单链蛋白的单点突变可直接调控构象平衡和转运功能，而同源二聚体的双点突变会产生抵消效应，限制了其进化潜力。本研究不仅为膜蛋白构象动力学研究提供了新的定量工具，也为理解膜转运蛋白的进化机制和开发靶向药物奠定了重要基础。
9.芬兰Bioscreen仪器测量数据的研究意义
本研究使用芬兰Bioscreen Pro C全自动微生物生长曲线分析仪，在37℃、缓慢振荡条件下，每15分钟自动测定一次600 nm处的光密度（OD₆₀₀），连续监测12-24小时，获得了野生型EmrE、EmrE^W63突变体、E14D失活突变体和空载体对照在不同浓度溴化乙锭、原黄素和甲基紫精存在下的完整生长曲线，数据是图1c和图S1中IC₅₀值的核心基础。需要特别说明的是，本研究中Bioscreen并非用于测定微生物生长曲线本身，而是通过生长曲线间接评估EmrE转运蛋白的功能活性，其研究意义主要体现在以下几个方面：
高通量平行验证突变体功能：Bioscreen Pro C可同时处理多达200个样品，本研究中对4种菌株、3种药物、8个浓度梯度均设置了至少4个生物学重复，在短时间内完成了上百个样品的生长动力学测定，大幅提高了实验效率，避免了传统分光光度计逐个测定的繁琐和人为误差。
精确量化EmrE的转运活性：通过拟合生长曲线的四参数逻辑斯蒂（4PL）方程，计算得到了各菌株对不同药物的半抑制浓度（IC₅₀）。结果显示，EmrE^W63突变体对溴化乙锭的IC₅₀为102.0 μM，与野生型的159.6 μM接近，而E14D失活突变体的IC₅₀与空载体对照无显著差异。这一结果严格证明了EmrE^W63突变体保留了完整的多药外排功能，为后续使用该突变体进行NMR构象研究提供了关键的功能验证，确保了实验结果的生理相关性。
动态反映药物作用的时间进程：仪器的连续监测功能不仅能获得终点OD值，还能记录整个生长过程的动力学变化，反映药物对细菌生长的抑制速率和程度。这对于评估转运蛋白的外排效率至关重要，因为细菌的生长状态直接与细胞内药物浓度相关，而细胞内药物浓度又取决于EmrE的外排活性。
标准化实验条件保证数据可靠性：Bioscreen的全自动化操作（自动孵育、自动振荡、自动读数）消除了人工操作带来的随机误差，统一的温度、振荡速度和通气条件保证了不同批次实验结果的可比性。本研究中所有生长抑制实验均在相同条件下进行，获得的IC₅₀值具有高度的重复性和统计学显著性（P-value < 0.0001）。
建立功能评估的基准体系：通过测定空载体对照和E14D失活突变体的IC₅₀值，建立了EmrE功能完全丧失的基准线；通过测定野生型EmrE的IC₅₀值，建立了正常功能的基准线。这一基准体系为后续评估其他EmrE突变体的功能提供了明确的参照，也为研究底物特异性与构象变化的关系提供了背景信息。