AMBER Archive (2007)

Subject: Re: AMBER: Radii for water hydrogens are different

From: David A. Case (case_at_scripps.edu)
Date: Fri Jan 19 2007 - 12:13:41 CST


On Tue, Jan 16, 2007, Steve Spronk wrote:

> I stumbled upon what I think is an interesting bug in Amber 8. I used tLEaP
> to set up a system consisting of a protein solvated in TIP3P water, using
> ff03. The PDB that I used still had its crystallographic waters, and I used
> the solvateoct command to add a waterbox.
>
> In the prmtop file that tLEaP wrote out, the %FLAG RADII section showed that
> for the crystallographic waters, the radii for the hydrogens were 0.8 and
> 1.2. They had different radii.

Thanks for the report. We will try to see how to fix this. Note, however,
that these radii are not used for explicit solvent simulations, so that in
your case, the problem is innocuous.

Even for GB simulations (where the RADII are important), we have never
parameterized values for water itself. Still, the program should not be
putting out such weird values.

...regards...dac

-----------------------------------------------------------------------
The AMBER Mail Reflector
To post, send mail to amber_at_scripps.edu
To unsubscribe, send "unsubscribe amber" to majordomo_at_scripps.edu