从 pdb2pqr 生成的 PQR 文件与 GROMACS 不同

计算科学 分子动力学 计算生物学 电磁学 数据库
2021-12-17 00:56:53

我有一个 pdb 文件,我想做一个 pqr 文件。在 GROMACS 中,我editconf -mead用来生成 pqr 文件,topol.top 取自 pdb2gmx 并使用它来生成 pdb 文件,然后使用 pdb2pqr。哪个更好?

由 GROMACS 制成的前 10 行 PQR

ATOM  1 N    PRO A   1     0.63000  39.63000  17.25000 -0.2020  1.6250      
ATOM  2 H1   PRO A   1     1.37000  39.93000  17.89000  0.3120  0.5345      
ATOM  3 H2   PRO A   1     0.24000  40.48000  16.88000  0.3120  0.5345      
ATOM  4 CD   PRO A   1     1.21000  38.79000  16.18000 -0.0120  1.6998      
ATOM  5 HD1  PRO A   1     2.30000  38.76000  16.29000  0.1000  0.9800      
ATOM  6 HD2  PRO A   1     0.96000  39.21000  15.20000  0.1000  0.9800      
ATOM  7 CG   PRO A   1     0.63000  37.39000  16.30000 -0.1210  1.6998      
ATOM  8 HG1  PRO A   1     1.28000  36.80000  16.95000  0.1000  1.3248      
ATOM  9 HG2  PRO A   1     0.52000  36.90000  15.33000  0.1000  1.3248      
ATOM  10 CB  PRO A   1    -0.72000  37.64000  16.99000 -0.1150  1.6998      

由 pqd2pqr 生成的前 10 行 PQR:

ATOM      1  N   PRO     1       0.630  39.630  17.250 -0.2020 1.8240
ATOM      2  H   PRO     1       1.370  39.930  17.890  0.3120 0.6000
ATOM      3  H2  PRO     1       0.240  40.480  16.880  0.3120 0.6000
ATOM      4  CD  PRO     1       1.210  38.790  16.180 -0.0120 1.9080
ATOM      5  HD3 PRO     1       2.300  38.760  16.290  0.1000 1.1000
ATOM      6  HD2 PRO     1       0.960  39.210  15.200  0.1000 1.1000
ATOM      7  CG  PRO     1       0.630  37.390  16.300 -0.1210 1.9080
ATOM      8  HG3 PRO     1       1.280  36.800  16.950  0.1000 1.4870
ATOM      9  HG2 PRO     1       0.520  36.900  15.330  0.1000 1.4870
ATOM     10  CB  PRO     1      -0.720  37.640  16.990 -0.1150 1.9080

这些来自 3B7V,一种 HIV-1 蛋白酶。它已被最小化的能量。

1个回答

Van der Walls 半径(输出中的最后一列)是根据力场计算得出的。可能 pdb2pqr 和 editconf 使用不同的力场,因此不同的半径。我不使用 pdb2pqr,但似乎 (1) AMBER99 是默认力场,尽管支持 CHARMM、PARSE 和 TYL06。Gromacs 的 editconf 从 pdb2gmx (2) 生成的拓扑文件中读取力场。在我的系统中 pdb2gmx 默认为 AMBER03。尝试

pdb2gmx -ff select

哪个更好取决于您的问题。您是否正在为 APBS 生成输入?如果是这样,我会使用 pdb2pqr。

参考:

(1) http://www.poissonboltzmann.org/pdb2pqr/user-guide/using-pdb2pqr

(2) http://manual.gromacs.org/online/pdb2gmx.html