《Journal of Chemical Theory and Computation》：Environment-Conditioned Design of α-Helical Peptides
编辑推荐：
生物界面上的肽常常通过优先兼容特定化学环境，而不是通过结合单一分子靶标来编码功能。研究人员报道了一个物理动机明确、基于知识的哈密顿量（Hamiltonian，体系能量函数），用于在α-螺旋肽中编程这类环境偏好。该模型将氨基酸序列映射到一个归一化能量景观，其中组
生物界面上的肽常常通过优先兼容特定化学环境，而不是通过结合单一分子靶标来编码功能。研究人员报道了一个物理动机明确、基于知识的哈密顿量（Hamiltonian，体系能量函数），用于在α-螺旋肽中编程这类环境偏好。该模型将氨基酸序列映射到一个归一化能量景观，其中组合了螺旋倾向性（helix propensity）、局部序列环境（local sequence context）、残基接触（residue contacts）、溶剂暴露（solvent exposure）、膜电荷耦合（membrane charge coupling）和静电屏蔽（electrostatic screening）。随机参照Z分数（Z-score）将不同长度、残基组合和环境置于同一尺度上。在均匀介质基准中，该哈密顿量能够将实验表征的富含丙氨酸的Glu/Lys水溶性螺旋与WALP/GWALP跨膜模型螺旋区分开。在界面基准中，它识别出经典的膜活性抗菌肽（antimicrobial peptides，AMPs）在阴离子界面处比在中性、阳离子、水性和非极性替代环境中更被优先稳定，同时拒绝非抗菌的螺旋对照。同一评分可生成用于环境特异性兼容性、极性–非极性稳健性和一对多（one-vs-rest）阴离子选择性的序列库。低能集合展现了两亲性分离、电荷调控以及序列间隔为3和4个残基处的疏水性周期性，将设计序列与可解释的化学组织联系起来。由于目标函数是离散哈密顿量，它可以通过精确搜索、经典启发式算法以及二进制编码的量子兼容工作流进行探索，而无需改变评分函数。该框架提供了一条透明的途径，用于设计其环境偏好直接在分子目标中指定的α-螺旋肽。
在物理、化学与计算交叉领域，设计在不同物理环境中保持稳定且功能兼容的生物分子序列是一项基础挑战。序列稳定性并非仅由内在折叠偏好决定，还由与溶剂极性、静电和界面几何等周围环境的相互作用所决定。膜相互作用α-螺旋肽，尤其是抗菌肽（antimicrobial peptides，AMPs），其功能依赖于在非均质界面上选择性稳定两亲构象的能力；它们通常不通过结合特定蛋白靶标发挥作用，而是通过识别并破坏特定膜环境。因此，这类肽的活性取决于其对阴离子细菌或病变细胞膜相对健康哺乳动物中性膜的优先兼容性。现有深度学习方法，如AlphaFold 3等，虽然强大，但多数依赖统计关联、同源信息或大规模训练数据，难以直接揭示新设计场景中溶剂极性、电荷与膜几何之间的竞争；而机制显式打分方法又面临序列–构象空间穷举困难。为此，研究人员引入了一个物理动机明确、基于知识的哈密顿量框架，用于膜平行α-螺旋支架上的环境条件化序列设计。该工作将环境偏好直接写入分子目标，发表于《Journal of Chemical Theory and Computation》。
研究人员使用的核心方法包括：(1) 构建基于知识的哈密顿量（Hamiltonian，体系能量函数），将螺旋倾向性、局部序列环境、成对残基接触、溶剂极性、膜电荷耦合和静电屏蔽组合为归一化能量景观；(2) 用随机参照Z分数（Z-score）标准化，使不同长度、字母表和环境可比较；(3) 用二进制编码将同一评分函数映射为经典精确搜索、经典启发式或量子近似优化算法（Quantum Approximate Optimization Algorithm，QAOA）可执行的代价函数；(4) 用横截疏水矩作为几何对齐和景观分析描述符；(5) 以实验表征的水溶性/跨膜肽和膜活性抗菌肽基准及非AMP阴性对照进行验证。研究未涉及样本队列来源。
### 研究结果
**2.1 水溶性与跨膜模型肽的识别（Recognition of Water-Soluble and Transmembrane Model Peptides）**
研究人员将实验表征的富含丙氨酸Glu/Lys水溶性螺旋肽与WALP/GWALP跨膜模型肽分别在均匀极性和均匀非极性哈密顿量中评分。结果表明，水溶性参考肽在水相具有有利Z分数（Z-score）、在非极性相具有不利Z分数，ΔZ＝Z
apolar?Z
water介于3.16～4.14；跨膜参考肽则相反，ΔZ介于?3.46～?5.46。仅使用ΔZ的符号即可将所有参考序列归入预期环境类别，支持均匀极性与非极性打分项的物理解释。
**2.2 经典膜活性肽的识别（Recognition of Canonical Membrane-Active Peptides）**
研究人员将LL-37、Melittin、Magainin 2、Cecropin A、PGLa、Piscidin 1和BP100七种膜活性肽，与(EAAAK)3和Poly-Ala 15两种非AMP螺旋对照，在阴离子、中性和阳离子界面扫描（5%～95%穿透）中评分。所有AMPs在阴离子界面的最优Z分数（?1.71～?5.56σ）均优于均匀水相或均匀非极性相，也优于中性或阳离子界面最优值；两个非AMP对照则更倾向于水相。这表明该哈密顿量识别的不是通用螺旋性，而是阴离子膜活性两亲组织。
**2.3 环境条件化序列设计（Environment-Conditioned Sequence Design）**
在L=20的联合均匀介质目标下，序列空间退火产生从强极性偏向序列EKEKFEKEKEKEKFEKEKEK（Z
polar=?9.26，Z
apolar=?1.29）到随λ
gap增加而变为混合芳香/带电的两环境兼容序列；λ
gap=4时归一化极性–非极性能量差降至0.02。在一对多（one-vs-rest）阴离子选择性设计中，λ=0.25、0.5、0.75在固定穿透深度下恢复同一顶级候选，50%穿透为最佳折中；硬选择性极限λ=0则产生浅穿透分离更强的不同序列家族。这些序列是供后续验证的模型预测。
**2.4 序列-景观结构：漏斗、周期性与中性网络（Sequence-Landscape Structure: Funnels, Periodicity, and Neutral Networks）**
对50000条随机L=12序列的景观分析显示，低能序列富集较大的横截疏水矩，净电荷受静电项约束。对1000条随机诱饵
