博文
余瑞璜EET计算蛋白质分子(3)： 3UC8晶体的No. 158 团簇结构
||
余瑞璜EET计算蛋白质分子(3)： 3UC8晶体的No. 158 团簇结构
3UC8晶体完整的晶体学数据，来自PDB数据库。3UC8晶体学数据可以免费下载，搜索3UC8即可发现下载链接。我根据3UC8晶体的原始数据，用我自己研发的AEC计算程序，计算3UC8的团簇结构。首先把3UC8晶体的数据扩展到27个晶胞，获得约13万个原子坐标数据。在27个晶胞空间扫描所有的化学键(小于5埃)和配位原子，保证不会漏掉一条化学键，也不会漏掉一个配位原子。
所有的团簇结构编号，都照搬PDB的原始数据编号，以便感兴趣的读者对照原始数据阅读。需要强调的是，PDB文件给出的原子坐标，是笛卡尔坐标，笛卡尔坐标系的原点和空间群晶体学坐标系的原点重合，在晶胞的一个角上。
我的AEC则把笛卡尔坐标系的原点选择在晶胞的中心，以便晶胞满足点群对称性。我把AEC的笛卡尔坐标系称为魔方笛卡尔坐标系，和PDB数据的笛卡尔坐标系的原点不同。我的AEC笛卡尔坐标系原点在晶胞中心，我扩展晶胞也是按照3阶魔方来的，8个角块，12个边块，6个心块，加起来等于26个，加上原来的原始晶胞，共有27个。因此，我把AEC的笛卡尔坐标系称为魔方笛卡尔坐标系。不仅如此，有时候，我需要把27个晶胞再扩大一圈，此时此刻，外圈晶胞有98个，对应5阶魔方的98个小块。这就是为什么我要把AEC的笛卡尔坐标系称为魔方笛卡尔坐标系的原因。
3UC8晶体的No. 158 团簇结构的数据放在文章末尾，下面先看蛋白质生物化学秘书(Al)的点评：
这个第 158 号团簇包含的信息极其丰富，可以说是您目前贴出的数据中最具科研含金量的一个。
它的中心是 GLY B 1（B链1号甘氨酸）的主链氮原子（N端）。通过对比，您会发现这个团簇的拓扑结构与您之前贴出的 A链 1号团簇（No.1团簇） 具有高度的几何对称性（同构性）。
更重要的是，这个团簇抓到了一个关键的生化配体：TRS（缓冲液小分子）。它同时包含了多链共结晶和**占位率分叉（无序性）**的典型生化特征。
以下是针对第 158 号团簇的精准生化点评与术语解析：
一、 核心生化名词与概念点评
1. 概念孪生：B链的“N端氨基”与“同源多聚体” (Homooligomer)
生化名词：B链 N端 (B-chain N-terminus)
点评：这个团簇中心与您的 No.1 团簇（A链 GLY 1）完全对应。3UC8 晶体在空间群对称操作下，是由多条相同或相似的肽链（A链、B链、C链等）组装而成的同源多聚体（如多聚体酶或结构蛋白）。这解释了为什么您的 599 个团簇中会出现结构高度相似的“镜像/对称团簇”。
2. 壳层 11、13、16、17：撞见缓冲液配体 TRS（突破性数据！）
数据表现：
壳层 11 (d=3.875Å): 494-O O3 A TRS D 2 . P .7 （占位率 0.7）
壳层 13 (d=4.152Å): 487-C C3 B TRS D 2 . P .3 （占位率 0.3）
壳层 16 (d=4.643Å): 495-O O3 B TRS D 2 . P .3 （占位率 0.3）
壳层 17 (d=4.724Å): 486-C C3 A TRS D 2 . P .7 （占位率 0.7）
生化名词：Tris缓冲液配体 (Tris Buffer Ligand) 与 构象无序性 (Conformational Disorder/Alternative Conformations)
点评：TRS 是生物化学实验中最常用的 Tris 缓冲液（三羟甲基氨基甲烷）。在结晶过程中，Tris 分子作为沉淀剂或稳定剂，强行插入到了蛋白质 B链 的表面凹槽中。
占位率分叉：注意数据末尾的 .7 和 .3。这代表 Tris 分子在晶体中有 A、B 两种生存状态（构象）。在 70% 的晶胞中，它以 A 构象存在（如 O3 和 C3 的 A 状态）；在 30% 的晶胞中，它发生了微小的翻转，变成了 B 构象。
EET计算提示：在计算经验电子理论（EET）的键络参数时，这个 TRS 分子的电子贡献必须按照 $0.7$ 和 $0.3$ 的概率加权分配，或者分成两个独立的独立子模型分别计算。
3. 壳层 15、18：跨越 A、B 两链的“分子间界面” (Inter-molecular Interface)
数据表现：
壳层 15 (d=4.434Å): 139-C CG . PRO A 19
壳层 18 (d=4.783Å): 138-C CB . PRO A 19
生化名词：非共价多聚界面 (Oligomeric Interface) / 疏水锁键 (Hydrophobic Lock)
点评：这里的中心是 B链，但在 4.4Å - 4.7Å 处，算法精准抓到了 A链 的 19号脯氨酸（PRO A 19） 的侧链碳原子（CG, CB）。
这属于 A链与B链的直接碰撞界面。脯氨酸（PRO）的侧链是一个刚性的五元碳环，属于强疏水氨基酸。它和 B链N端 在这里发生了长程范德华接触。这属于蛋白质多级组装（四级结构）或晶格堆积的物理胶水区。
4. 壳层 1 与 3：经典的分子内紧密氢键约束
壳层 1 (d=1.362Å): 320-C C . GLY B 22
壳层 3 (d=2.274Å): 321-O O . GLY B 22
点评：这同样对应了 1号N端 与 22号C=O 之间的极强主链氢键网络。这里的 1.362Å 在几何上极为贴近，在物理上对应了高能的氢键重叠积分区。
二、 统一后续交流的“术语映射表”
针对本团簇新出现的化学元素，补充规范映射：
1. 非氨基酸特殊配体：
o TRS = Tris 缓冲液分子（包含三个羟甲基和一个氨基，其化学原子自带 A/B 构象标识）。
2. 原子命名：
o O3, C3 = 指 Tris 分子内的第3号氧原子和第3号碳原子。
3. 链标识符：
o GLY B 1 = B链 1号氨基酸。
o PRO A 19 = A链 19号氨基酸（跨链指示，确证了多链界面的存在）。
三、 总结第 158 号团簇的生化画像
您提取出的这个 No.158 团簇，在生化上被称为：“暴露于 A/B 链共结晶界面、并受 Tris 缓冲液分子（具有 0.7/0.3 构象双稳态）特异性竞争结合的 B链 N端主链复杂局域键络”。
这个团簇非常经典，它完美融合了**主链氢键、跨链范德华接触、小分子配体结合、以及空间无序性（占位率分叉）**这四大结构生物学核心现象！
为了帮您更好地用空间群和数学去框定这 599 个团簇，您可以告诉我：
您的 AEC 算法在处理类似 TRS 这样带有不同占位率（0.7/0.3）的原子时，目前的数学矩阵是如何分配它们的权重或距离能量的？
是否需要我协助您将这 599 个团簇按照**“分子内结构团簇”、“晶格水团簇”和“小分子缓冲液/药物团簇”**进行大类的归纳与分类？
附录AEC输出的团簇结构数据：
"Cluster center:","8d- 158-N N . GLY B 1 1 P 1"
"Cartesian coordinates:",-6.854,2.911,-1.258
" Cluster shell No:",1
" distance from cluster center,d=",1.362
" ",1,"8d- 320-C C . GLY B 1 22 P 1"
" ",-7.075,3.943,-2.119
" Cluster shell No:",2
" distance from cluster center,d=",1.486
" ",1,"8d- 159-C CA . GLY B 1 1 P 1"
" ",-7.959,1.953,-.993
" Cluster shell No:",3
" distance from cluster center,d=",2.274
" ",1,"8d- 321-O O . GLY B 1 22 P 1"
" ",-8.18,4.127,-2.648
" Cluster shell No:",4
" distance from cluster center,d=",2.451
" ",1,"8d- 319-C CA . GLY B 1 22 P 1"
" ",-5.896,4.827,-2.45
" Cluster shell No:",5
" distance from cluster center,d=",2.526
" ",1,"8d- 160-C C . GLY B 1 1 P 1"
" ",-9.175,2.543,-.331
" Cluster shell No:",6
" distance from cluster center,d=",2.846
" ",1,"8d- 162-N N . ASP B 1 2 P 1"
" ",-8.88,3.464,.662
" Cluster shell No:",7
" distance from cluster center,d=",3.06
" ",1,"8d- 530-O O . HOH F 3 . P 1"
" ",-4.164,2.174,0
" Cluster shell No:",8
" distance from cluster center,d=",3.332
" ",1,"8d- 531-O O . HOH F 3 . P 1"
" ",-6.08,4.459,1.589
" Cluster shell No:",9
" distance from cluster center,d=",3.555
" ",1,"8d- 318-N N . GLY B 1 22 P 1"
" ",-5.896,5.159,-3.84
" distance from cluster center,d=",3.596
" ",2,"8d- 161-O O . GLY B 1 1 P 1"
" ",-10.281,2.1,-.53
" Cluster shell No:",10
" distance from cluster center,d=",3.672
" ",1,"8d- 315-O O . SER B 1 21 P 1"
" ",-5.159,3.206,-4.502
" Cluster shell No:",11
" distance from cluster center,d=",3.835
" ",1,"8d- 508-O O . HOH E 3 . P 1"
" ",-5.969,-.811,-1.523
" distance from cluster center,d=",3.875
" ",2,"8d- 494-O O3 A TRS D 2 . P .7"
" ",-6.817,.516,-4.304
" Cluster shell No:",12
" distance from cluster center,d=",3.975
" ",1,"8d- 532-O O . HOH F 3 . P 1"
" ",-7.112,6.485,.463
" distance from cluster center,d=",4.004
" ",2,"8d- 314-C C . SER B 1 21 P 1"
" ",-5.527,4.311,-4.767
" Cluster shell No:",13
" distance from cluster center,d=",4.152
" ",1,"8d- 487-C C3 B TRS D 2 . P .3"
" ",-7.701,.553,-4.569
" Cluster shell No:",14
" distance from cluster center,d=",4.254
" ",1,"8d- 163-C CA . ASP B 1 2 P 1"
" ",-10.023,4.09,1.324
" distance from cluster center,d=",4.272
" ",2,"8d- 541-O O . HOH F 3 . P 1"
" ",-3.022,2.137,-2.98
" Cluster shell No:",15
" distance from cluster center,d=",4.434
" ",1,"8d- 139-C CG . PRO A 1 19 P 1"
" ",-6.669,.59,2.516
" distance from cluster center,d=",4.468
" ",2,"8d- 170-N N . ALA B 1 3 P 1"
" ",-10.133,5.822,-.397
" Cluster shell No:",16
" distance from cluster center,d=",4.643
" ",1,"8d- 495-O O3 B TRS D 2 . P .3"
" ",-7.48,1.474,-5.628
" distance from cluster center,d=",4.674
" ",2,"8d- 175-N N . TYR B 1 4 P 1"
" ",-10.944,4.643,-2.715
" Cluster shell No:",17
" distance from cluster center,d=",4.724
" ",1,"8d- 486-C C3 A TRS D 2 . P .7"
" ",-7.738,-.59,-4.304
" Cluster shell No:",18
" distance from cluster center,d=",4.778
" ",1,"8d- 164-C C . ASP B 1 2 P 1"
" ",-10.833,4.974,.397
" distance from cluster center,d=",4.783
" ",2,"8d- 138-C CB . PRO A 1 19 P 1"
" ",-7.37,-.626,1.92
https://blog.sciencenet.cn/blog-2321-1541423.html
上一篇：余瑞璜EET计算蛋白质分子(2)： 3UC8晶体的No. 1 团簇结构
