Yes, you need the activation barrier to estimate your rate.
This can be obtained by computing a PMF, either via Umbrella Sampling,
SMD, metadynamics, ...
You should (must) obtain similar PMF with these methods, provided that
you choose the same reaction coordinate approximation (a distance for
example). The PMF is a physical quantity and, in principle, does not
depend on the way you compute it.
We are actually making comparison between US and SMD to study gas
diffusions within a protein and we obtain the same PMFs.
Jerome
On Thu, 2010-05-20 at 17:29 +0800, Catein Catherine wrote:
> Hi Jerome and Thomas,
>
>
>
> Thanks for your comments. I am interested to study the rate of drug binding and drug removal. To my understanding, it is related to the activation energy of the processes, which is pathway-dependent. Am I correct?
>
>
>
> Can I use PMF to calculate the activation energy?
>
>
>
> If the answer is yes, should I expect to see the activation energies are different when I use these methods?
>
>
>
> If the answer is no, are there any reliable method to calculate activation energy of the drug binding and drug removal process.
>
>
>
> Best regards and thanks a lot,
>
>
>
> Cat
>
> > Subject: Re: [AMBER] umbrella sampling method or steered MD should be used?
> > From: jerome.golebiowski.unice.fr
> > To: amber.ambermd.org
> > Date: Thu, 20 May 2010 10:24:54 +0200
> >
> > As Thomas said, you should see the same PMF for both the forward and the
> > backward directions.
> > Note that there is a (generally small) difference between the binding
> > free energy and the difference between two points of the PMF curve.
> > Some authors (including me) consider that the delta(PMF) equals the
> > delta(G), but in principle, it's not totally true.
> > Best
> > Jerome
> >
> >
> > On Thu, 2010-05-20 at 14:27 +0800, Catein Catherine wrote:
> > > Dear Sir/Madam,
> > >
> > >
> > >
> > > I would like to determine the potential energy changes for drug binding the receptor and for the drug remove from the receptor. In practice, the PMF of binding should be lower than that of removing, otherwise no drug-receptor complex could be formed.
> > >
> > >
> > >
> > > I tried to use umbrella sampling method to determine the PMF as a function of distance between the drug and receptor. As the drug removal and drug binding process use the same path (defined as the distance between the drug and the receptor). I can only see PMF for binding and removal are the same.
> > >
> > >
> > >
> > > I noted that I can define steered MD by start and end point. Should I expect to see the PMF obtained by binding and removal of drugs are different using these two different SMD restraints?
> > >
> > >
> > >
> > > (1)
> > >
> > > # Change distance between atoms 485 and 134 from 15 A to 20 A
> > > &rst iat=485,134, r2=15., rk2 = 5000., r2a=20. /
> > >
> > >
> > > (2)
> > >
> > > # Change distance between atoms 485 and 134 from 20 A to 15 A
> > > &rst iat=485,134, r2=20., rk2 = 5000., r2a=15. /
> > >
> > > Please kindly help.
> > >
> > > Cat
> > >
> > >
> > >
> > > _________________________________________________________________
> > > Hotmail: Trusted email with Microsoft’s powerful SPAM protection.
> > > https://signup.live.com/signup.aspx?id=60969
> > > _______________________________________________
> > > AMBER mailing list
> > > AMBER.ambermd.org
> > > http://lists.ambermd.org/mailman/listinfo/amber
> > --
> > Jérôme Golebiowski, PhD (Hab. Dir. Rech.)
> > LCMBA, UMR CNRS-UNS 6001
> > Molecular Modeling team
> > University of Nice Sophia Antipolis, Parc Valrose
> > 06108 Nice cedex 2
> > http://www.unice.fr/lcmba/golebiowski
> > tel +33 (0)4 92 07 61 03
> > fax +33 (0)4 92 07 61 25
> >
> >
> > _______________________________________________
> > AMBER mailing list
> > AMBER.ambermd.org
> > http://lists.ambermd.org/mailman/listinfo/amber
>
> _________________________________________________________________
> Your E-mail and More On-the-Go. Get Windows Live Hotmail Free.
> https://signup.live.com/signup.aspx?id=60969_______________________________________________
> AMBER mailing list
> AMBER.ambermd.org
> http://lists.ambermd.org/mailman/listinfo/amber
>
--
Jérôme Golebiowski, PhD (Hab. Dir. Rech.)
LCMBA, UMR CNRS-UNS 6001
Molecular Modeling team
University of Nice Sophia Antipolis, Parc Valrose
06108 Nice cedex 2
http://www.unice.fr/lcmba/golebiowski
tel +33 (0)4 92 07 61 03
fax +33 (0)4 92 07 61 25
_______________________________________________
AMBER mailing list
AMBER.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber
Received on Thu May 20 2010 - 03:30:03 PDT
