lxhpdx 发表于 2024-6-12 18:53

问题已解决

1. # 选择与名为"MOL"的残基距离在5埃之内的蛋白质原子
   
2. set sel [atomselect top "same residue as {protein within 5 of resname MOL}"]
   
3. 
   
4. # 循环处理200帧的数据
   
5. for {set i 0} {$i < 200} {incr i} {
   
6.     # 设置当前帧
   
7.     $sel frame $i
   
8.     # 更新选择
   
9.     $sel update
   
10.    
    
11.     # 获取当前帧中选中的残基ID
    
12.     set residue_ids [$sel get resid]
    
13.    
    
14.     # 对残基ID进行去重，得到当前帧中的唯一残基ID
    
15.     set unique_residue_ids [lsort -unique $residue_ids]
    
16.    
    
17.     # 遍历唯一的残基ID
    
18.     foreach resid $unique_residue_ids {
    
19.         # 如果该残基ID已经存在于residue_frames数组中，则添加当前帧号
    
20.         if {[info exists residue_frames($resid)]} {
    
21.             lappend residue_frames($resid) $i
    
22.         } else {
    
23.             # 否则，初始化该残基ID的帧集合并添加当前帧号
    
24.             set residue_frames($resid) [list $i]
    
25.         }
    
26.     }
    
27. }
    
28. 
    
29. # 输出每个残基ID以及其出现的帧数
    
30. foreach resid [array names residue_frames] {
    
31.     set frame_count [llength [lsort -unique $residue_frames($resid)]]
    
32.     puts "Residue ID: $resid, Total Frames: $frame_count"
    
33. }

复制代码