1
##############################################################################
3
# Open Babel file: phmodel.txt #
5
# Copyright (c) 1998-2001 by OpenEye Scientific Software, Inc. #
6
# Some portions Copyright (c) 2001-2003 by Geoffrey R. Hutchison #
7
# Part of the Open Babel package, under the GNU General Public License (GPL)#
9
# pH model data (used by phmodel.cpp:OBPhModel) #
11
# TRANSFORM: chemical transforms can be used to modify formal charges, bond #
12
# orders, and to delete atoms (i.e hydrogens). Changes are applied#
13
# to vector bound atoms (use the [:#] SMARTS notation) and bonds #
14
# between vector bound atoms. #
15
# SEEDCHARGE: used to seed partial charges. Seed partial charges #
16
# are used as initial values in Gasteiger charge calculation #
18
##############################################################################
21
TRANSFORM O=CO[#1:1] >> O=CO
22
TRANSFORM O=C[OD1-0:1] >> O=C[O-:1]
25
TRANSFORM [N^3;!$(N~[!#6;!#1]):1] >> [N+:1]
27
#imidazole: note pKa=7.0, histidine pKa=6.0
28
#if you uncomment this, also uncomment the seedcharge statement below
29
#TRANSFORM [nD2:1]1c[nH]cc1 >> [n+:1]1c[nH]cc1
32
TRANSFORM [ND3+0:1]=[#6] >> [ND3+:1]=[#6]
35
TRANSFORM [nD2:1]([#1:2])1[nD2-0][nD2-0][nD2-0]c1 >> [n-:1]1nnnc1
36
TRANSFORM [nD2-0]1[nD2:1]([#1:2])[nD2-0][nD2-0]c1 >> n1[n-:1]nnc1
37
TRANSFORM [nD2-0:1]1[nD2-0][nD2-0][nD2-0]c1 >> [n-:1]1nnnc1
40
TRANSFORM [ND2:1]=[ND2:2]=A >> [N:1]=[N+:2]=A
41
TRANSFORM [ND1:1]=[ND2:2]=A >> [N-:1]=[N+:2]=A
44
TRANSFORM O=CN[OD1-0:1][#1:2] >> O=CN[O-:1]
45
TRANSFORM O=CN[OD1-0:1] >> O=CN[O-:1]
48
TRANSFORM [SD3](=O)[OD1:1] >> [SD3](=O)[O-:1]
49
TRANSFORM [SD3](=O)[O:1][#1:2] >> [SD3](=O)[O-:1]
52
TRANSFORM [SD4]([!D1])(=O)(=O)[OD1:1] >> [SD4]([!D1])(=O)(=O)[O-:1]
53
TRANSFORM [SD4]([!D1])(=O)(=O)[O:1][#1:2] >> [SD4]([!D1])(=O)(=O)[O-:1]
56
TRANSFORM [SD4]([!D1])(=O)(=O)[OD1:1] >> [SD4]([!D1])(=O)(=O)[O-:1]
57
TRANSFORM [SD4]([!D1])(=O)(=O)[O:1][#1:2] >> [SD4]([!D1])(=O)(=O)[O-:1]
60
TRANSFORM [#6^2+0:1](~[N^2])(~[N^2])* >> [#6+:1](~N)(~N)*
63
TRANSFORM [nD2:1]1c[nD2]cc1 >> [n+:1]1c[nD2]cc1
64
# uncomment for tryptophan
65
# TRANSFORM [nD2:1]1cccc1 >> [n+:1]1cccc1
68
TRANSFORM [PD4](=O)([OD2])([OD2])[OD1:1] >> [PD4](=O)([OD2])([OD2])[O-:1]
69
TRANSFORM [PD4](=O)([OD2])([OD2])[OD1:1][#1:2] >> [PD4](=O)([OD2])([OD2])[O-:1]
72
TRANSFORM O=P([!D1])([O:1][#1:2])[O:3][#1:4] >> O=P([*D2,*D3])([O:1])[O:3]
73
TRANSFORM O=P([!D1])([O:1][#1:2])[OD1] >> O=P([!D1])([O:1])O
74
TRANSFORM O=P([*D2,*D3])([OD1:1])[OD1:2] >> O=P([*D2,*D3])([O-:1])[O-:2]
79
# Seeding partial charges for gasteiger calculation
89
SEEDCHARGE [#15-] -1.0
91
SEEDCHARGE [#16-] -1.0
93
#charges spread over multiple atoms
95
SEEDCHARGE C(=O)[O-] 0.0 -0.5 -0.5
100
SEEDCHARGE [nD2]1[nD2][nD2][nD2]c1 -0.2 -0.2 -0.2 -0.2 -0.2
103
SEEDCHARGE [SD3](=O)[O-] 0.0 -0.5 -0.5
105
SEEDCHARGE [SD4](=O)(=O)([O-])[OH1] 0.0 -0.33 -0.33 -0.33 0.0
107
SEEDCHARGE [SD4]([D2])(~[OD1])(~[OD1])~[OD1] 0.0 0.0 -0.33 -0.33 -0.33
110
SEEDCHARGE [#7^2]~[C^2](~[N^2])~[N^2] 0.33 0.0 0.33 0.33
113
SEEDCHARGE [#6]~[C^2](~[N^2])~[N^2] 0.0 0.0 0.5 0.5
116
SEEDCHARGE [PD4](=O)([OD2])([OD2])[OD1] 0.0 -0.5 0.0 0.0 -0.5
118
SEEDCHARGE O=P([!D1])([O-])[O-] -0.66 0.0 0.0 -0.66 -0.66
122
SEEDCHARGE P(=O)(=O)([O-])[O-] 0.0 0.0 0.0 -0.5 -0.5
124
SEEDCHARGE [#6]P(~[OD1])(~[OD1])~[OD1] 0.0 0.0 -0.33 -0.33 -0.33
127
SEEDCHARGE O=C[N;!$(N(C=O)C=O)][OD1] -0.5 0.0 0.0 -0.5
128
SEEDCHARGE O=CN([OD1])NC=O -0.33 0.0 0.0 -0.33 0.0 -0.33
130
#imidazole: note pKa=7.0, histidine pKa=6.0
131
#if you uncomment this, also uncomment the seedcharge statement below
132
#SEEDCHARGE [n+H]1c[nH]cc1 0.5 0.0 0.5 0.0 0.0