Re: [AMBER] MM-GBSA calculation in temperature other than 298.15 K

From: Daniel Roe <daniel.r.roe.gmail.com>
Date: Thu, 17 Nov 2016 16:11:13 +0000

What version of AmberTools are you using? I'm pretty certain the 298 K
temperature restriction was removed sone time ago.

-Dan

On Thu, Nov 17, 2016 at 6:45 AM Carlos Simmerling <
carlos.simmerling.gmail.com> wrote:

> then you need to read Jason's email and follow his advice.
>
> On Thu, Nov 17, 2016 at 8:22 AM, Atila Petrosian <
> atila.petrosian.gmail.com>
> wrote:
>
> > Excuse me. Considered temperature is 330 K in my project.
> >
> > On Thu, Nov 17, 2016 at 4:42 PM, Atila Petrosian <
> > atila.petrosian.gmail.com>
> > wrote:
> >
> > > Dear Carlos,
> > >
> > > Thanks for your answer.
> > >
> > > I did MD simulation of my system (protein-ligand complex) in 330 K.
> Thus,
> > > I have the trajectory data at the temperature (300 K) I want. Now, I
> want
> > > to do Binding free energy and have entropy in 300 K. How to do that?
> > >
> > > Any help will highly appreciated to help me to consider both of issues.
> > >
> > > Best,
> > > Atila
> > >
> > > On Thu, Nov 17, 2016 at 4:05 PM, Carlos Simmerling <
> > > carlos.simmerling.gmail.com> wrote:
> > >
> > >> the answers seem to address different things - Adrian mentions needing
> > to
> > >> have trajectory data at the temperature you want. Jason is talking
> about
> > >> the actual MM-PBSA program and how it needs to be changed if you want
> > >> entropies at a different T.
> > >> I think both of these are things you should consider.
> > >>
> > >>
> > >> On Thu, Nov 17, 2016 at 1:46 AM, Atila Petrosian <
> > >> atila.petrosian.gmail.com>
> > >> wrote:
> > >>
> > >> > Dear Carlos,
> > >> >
> > >> > I checked archive. Jason Swail answer is as follows:
> > >> >
> > >> > Not easily. You need to modify source codes to change the
> temperature.
> > >> The
> > >> > vibrational and rotational entropies are computed via statistical
> > >> > mechanical formulae that assume a particular temperature (and for
> > >> > quasi-harmonic entropies, you really need to use the temperature
> that
> > >> the
> > >> > simulation was run at). And last time I checked, this temperature
> was
> > >> > hard-coded directly into the program, so you would need to change
> the
> > >> > original source code and recompile, as well as
> > >> > modify the temperature variable in MMPBSA.py, in order to try a
> > >> different
> > >> > temperature.
> > >> >
> > >> > But I found same question from another person (Shahab Shariati):
> > MM-PBSA
> > >> > calculation in different temprature:
> > >> >
> > >> > Adrian Roitberg answer is as follows:
> > >> >
> > >> > Technically yes, you have to RERUN the MD at the new temperature and
> > >> rerun
> > >> > mmpbsa. You should not reuse the MD at 298 and simply change the
> > >> > temperature in mmbpsa. However, it is extremely unlikely that you
> will
> > >> see
> > >> > statistically meaningful differences between 298 and 310 with
> mmbpsa.
> > >> >
> > >> > I encontered with 2 different answers. Finally, is it possible to do
> > >> > MM-GBSA in temperature other than 298.15 K?
> > >> >
> > >> > Best,
> > >> > Atila
> > >> >
> > >> > On Wed, Nov 16, 2016 at 2:42 PM, Carlos Simmerling <
> > >> > carlos.simmerling.gmail.com> wrote:
> > >> >
> > >> > > Always check the archives before sending email. This has been
> > >> discussed
> > >> > > extensively before, and you can read the answers by typing
> "Mm-gbsa
> > >> > > temperature" in the search box on the amber web site.
> > >> > >
> > >> > > On Nov 16, 2016 12:42 AM, "Atila Petrosian" <
> > >> atila.petrosian.gmail.com>
> > >> > > wrote:
> > >> > >
> > >> > > > Dear Amber users.
> > >> > > >
> > >> > > > I did MD simulation of my system (protein-ligand complex) in 330
> > K.
> > >> > > >
> > >> > > > I used obtained trajectory for MM-GBSA calculation. But in
> output
> > >> file
> > >> > > > (*.dat), there is
> > >> > > >
> > >> > > > |All units are reported in kcal/mole.
> > >> > > > |All entropy results have units kcal/mol (Temperature is 298.15
> > K).
> > >> > > >
> > >> > > > I want to have Binding free energy in 330 K and not 298.15 K.
> > >> > > >
> > >> > > > How to resolve that.
> > >> > > >
> > >> > > > Best,
> > >> > > > Atila
> > >> > > > _______________________________________________
> > >> > > > 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
> > >> > >
> > >> > _______________________________________________
> > >> > 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
> > >>
> > >
> > >
> > _______________________________________________
> > 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
>
_______________________________________________
AMBER mailing list
AMBER.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber
Received on Thu Nov 17 2016 - 08:30:02 PST
Custom Search