Skip to content

Commit 23f5d29

Browse files
Made a note about HIS renaming in the docs
1 parent fdf3356 commit 23f5d29

File tree

1 file changed

+3
-2
lines changed

1 file changed

+3
-2
lines changed

docs/examples.rst

Lines changed: 3 additions & 2 deletions
Original file line numberDiff line numberDiff line change
@@ -200,8 +200,9 @@ works the same as for deca-alanine above, though ``plumed.dat`` and ``plumed_mon
200200
ADK contains three histidine residues. By default, these would be assigned by gmx pdb2gmx, however this is almost never a good idea
201201
because different histidine prototation states between umbrella sampling windows will crash replica exchange (because the atom ordering
202202
and connectivity are not identical). It is therefore recommended to first run pdb2gmx on one of the end states and use the resulting
203-
histidine protonation states for all intermediates, as was done in this example.
204-
203+
histidine protonation states for all intermediates, as was done in this example. In order to avoid problems, with pre-assigned histidine
204+
protonation states in your input coordinate files (which come with residue names like HISE, depending on the forcefield), it is often safest to
205+
rename all histidines to 'HIS' before passing to PyMEMENTO and reassigning them later.
205206

206207

207208
Handling lipids

0 commit comments

Comments
 (0)