If you are asking whether protonation states affect PT/REMD simulations,
then yes.
If you want to manually change the protonation state, edit the PDB residue
name
to be HID or HIP according to whichever you think is the most appropriate
(and delete
any extra hydrogens if necessary). Then load the modified PDB into LEaP.
-Dan
On Wed, Jul 10, 2013 at 2:04 AM, Francesco Pietra <chiendarret.gmail.com>wrote:
> Hello:
> I am preparing files for parallel tempering to hopefully get a
> dominant cluster of conformations for the unstructure (X-ray diffr)
> tails of a protein homotrimer.
>
> A problem: Poissson-Boltzmann calculations of pKa assign diprotonation
> to most histidines of these tails (generated with MODELLER), except
> for one subunit, where one histidine is assigned HIE. The apparent
> reason is that this particular histidine comes closer (according to
> MODELLER) to two other tails, while the other histidines of these
> tails feel the bulk solvent. There is no reason that this represents a
> low-energy situation.
>
> I expect that parallel tempering is biased by the protonation assigned
> to any titrable group for cases as illustrated above.
>
> Grateful for advice.
>
> francesco pietra
>
> _______________________________________________
> AMBER mailing list
> AMBER.ambermd.org
> http://lists.ambermd.org/mailman/listinfo/amber
>
--
Dr. Daniel J. Sindhikara <http://www.dansindhikara.com/Information.html>
Ritsumeikan University <http://www.ritsumei.ac.jp/eng/>
sindhikara.gmail.com <http://www.dansindhikara.com>
--
_______________________________________________
AMBER mailing list
AMBER.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber
Received on Tue Jul 09 2013 - 17:00:03 PDT
