OK. Thanks a lot!
Kefa
Thomas Cheatham wrote:
>> Thanks a lot! Dr. Cheatham. If I run constant (isotropic) pressure
>> simulation, will the situation be better? My system is a lipid bilayer
>>
>
> The same.
>
>
>> Another question is, if I don't use the original refc restraint
>> coordinates all the time, but change it every time when restart the
>> simulation, then the structure of the restraint part of the system might
>> change too much if I restart the simulation too many times, right? Will
>>
>
> Yes, but the change is small. If you are restraining a single large
> protein, what is happening as it runs is that the center of mass (assuming
> the default of NSPCAL=1) of the protein is shifting. Using the previous
> restart will have the shift, but the coordinates will not have changed so
> much. If you are restraining multiple molecules, this is a larger issue
> which can iterate offline.
>
>
>>>> GRP1
>>>> 100.0
>>>> ATOM 1 13609
>>>>
>
> Again, drop the force constants to 5-15 kcal/mol. This will still highly
> restrain atoms 1 13609 yet avoid the need for vlimit restrictions.
>
> --tom
>
> _______________________________________________
> AMBER mailing list
> AMBER.ambermd.org
> http://lists.ambermd.org/mailman/listinfo/amber
>
_______________________________________________
AMBER mailing list
AMBER.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber
Received on Wed Aug 19 2009 - 22:42:52 PDT
