查看: 908  |  回复: 0

Xu_shi_wei

新虫 (小有名气)

[交流] DNA的量子化分析

§7.3 PDB 生物大分子尺度:40 项独立构型

数据来自 PDB 公开晶体结构数据库,涵盖 DNA、RNA、蛋白质、多糖四大类,严格去重后保留 40 项独立构型。

去重标准 [21](Saenger)。 以下情况只保留一个代表:(i) 同一分子的不同 PDB 条目;(ii) 不同方法(X射线、NMR、冷冻电镜)测定的同一结构;(iii) 不同突变体但主链构象相同的结构。去重后的 40 项构型覆盖 DNA(21项,含 B-DNA、Z-DNA [20]、多核苷酸复合物 [22](Arnott & Selsing)及 G-四链体等)、RNA(2项)、蛋白质(12项 [24](Pauling et al.))、多糖(5项),全部为实验测定结构。

基本统计:

策略 命中数 命中率 零假设命中率 超出率 $p$ 值
先验 $L$ 38/40 95.0% 61.5%(条件概率) 33.5 个百分点 $1.21 \times 10^{-6}$
最小 dev 38/40 95.0% 77.0%(联合) 18.0 个百分点 $2.42 \times 10^{-3}$

按结构类型分组:

结构类型 样本量 命中数 命中率
DNA 21 20 95.2%
RNA 2 0 0%
蛋白质 12 11 91.7%
多糖 5 5 100%

按先验 $L$ 策略的 $Q_{\text{nearest}}$ 分布:

$Q_{\text{nearest}}$ 2 3 4 5 6 8
数量 4 8 6 11 10 1

40 个 PDB 结构的 $Q_{\text{nearest}}$ 分布与 $L=1$ 谱系的理论格点 $Q \in \{2,3,4,5,6,8\}$ 完全一致。其中 $Q=5$ 者最多(11 项),全部为 B-DNA 和 G-四链体结构;$Q=3$ 者(8 项)包括 Holliday junction、胶原蛋白 [23](Ramachandran & Kartha)、直链淀粉等;$Q=2$ 者(4 项)为纤维素 II、壳多糖、IHF-DNA bend 等大角度扭转结构;$Q=8$ 者仅 TMV 衣壳(扭转角 $22.5^\circ$)一项。

反例的结构性分析。 两个反例分别为:

· A-RNA(1EHZ,扭转角 $32.7^\circ$,偏离理论谱点 $L=1$、$Q=5$ 对应 $36^\circ$)
· π-螺旋蛋白(1R69,扭转角 $81.8^\circ$,偏差 $10.02\%$,略超容差)

A-RNA 与 A-DNA 共享同一偏离值 $32.7^\circ$,且两者的糖环构象均为 C3'-endo。C2'-endo 构象(B-DNA)的振荡振幅较小,锁相得以保持;C3'-endo 构象的振荡振幅较大,径向呼吸自由度主导了公转-自转锁相,导致扭转角偏离 $\pi/(LQ)$ 谱。两个反例共享同一偏离值,说明偏离根源不在碱基化学,而在糖环构象——这一致性本身即构成对理论的支撑,而非削弱。

在先验 $L$ 策略下,$p = 1.21 \times 10^{-6}$;在最小 dev 策略下,$p = 2.42 \times 10^{-3}$。两种策略均显著拒绝零假设。

https://doi.org/10.5281/zenodo.22688571
https://www.researchgate.net/publication/414173080
发自小木虫手机客户端
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 Xu_shi_wei 的主题更新
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 极限拉扯,如何抉择 +11 Isak0320 2026-10-12 15/750 2026-10-14 21:55 by 猫鱼不是鱼
[基金申请] 申请科研项目纯靠实力能成功吗? 200+5 David1982 2026-10-13 8/400 2026-10-14 21:38 by lhw_GPA4.0
[硕博家园] 中科院海洋所招收海洋材料相关领域博士研究生及研究助理 +3 sujiner 2026-10-08 4/200 2026-10-14 20:40 by qicaiyun
[职场人生] 40+中年人失业后,还能找到好工作吗? +3 xqy2009 2026-10-13 3/150 2026-10-14 20:17 by 43358052
[论文投稿] 二审4个多月了还是Awaiting Reviewer Scores 10+3 纳斯达克ML 2026-10-11 4/200 2026-10-14 17:30 by babero
[教师之家] 课题组会帮你通过考核么? +4 好运来~9527 2026-10-10 7/350 2026-10-13 08:41 by xli1984
[基金申请] 入职后第一次管课题组的专利缴费,对账才发现官费比别人多掏了一倍 +6 13108017953 2026-10-08 6/300 2026-10-12 20:07 by w3w3y3
[公派出国] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 9/450 2026-10-11 22:34 by XQNyDCa1Tk5n
[公派出国] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +5 qR2AdCWTVtkw 2026-10-09 10/500 2026-10-11 22:26 by XQNyDCa1Tk5n
[找工作] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 8/400 2026-10-11 22:26 by XQNyDCa1Tk5n
[考博] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 9/450 2026-10-11 22:15 by XQNyDCa1Tk5n
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 22:11 by XQNyDCa1Tk5n
[教师之家] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 (EPI+-1)(金币-50) +5 qR2AdCWTVtkw 2026-10-09 10/500 2026-10-11 22:06 by XQNyDCa1Tk5n
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 (EPI+-1)(金币-50) +5 qR2AdCWTVtkw 2026-10-09 10/500 2026-10-11 14:25 by VOMexwnDzOFZ
[教师之家] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 (EPI+-1)(金币-50) +4 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 14:05 by VOMexwnDzOFZ
[考研] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 8/400 2026-10-11 13:51 by VOMexwnDzOFZ
[硕博家园] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +5 qR2AdCWTVtkw 2026-10-09 8/400 2026-10-11 13:45 by VOMexwnDzOFZ
[找工作] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 10:54 by VOMexwnDzOFZ
[考博] 招收博士生一名 +3 mandolin 2026-10-10 3/150 2026-10-11 10:21 by 不是山谷!
[有机交流] 各位大佬,路线设计求助 +6 zapen 2026-10-08 8/400 2026-10-10 15:59 by zapen
信息提示
请填处理意见