Re: ABF Questions

From: Giacomo Fiorin (giacomo.fiorin_at_gmail.com)
Date: Mon Jul 23 2018 - 08:27:06 CDT

Hi Kelly, all of the methods implemented (including ABF) use the "colvars"
keyword to list the variables to use (can be more than one) out of those
defined. For the syntax, please read the manual:
http://colvars.github.io/colvars-refman-namd/colvars-refman-namd.html#sec:colvarbias

Thanks,
Giacomo

On Sat, Jul 21, 2018 at 3:06 PM McGuire, Kelly <mcg05004_at_byui.edu> wrote:

> Thanks Giacomo, that is good information. Could those additional
> collective variables all go in the same .in file and be biased
> simultaneously all in one ABF simulation? Or would these be separate
> simulations?
>
>
> *Kelly L. McGuire*
>
> *PhD Scholar*
>
> *Biophysics*
>
> *Department of Physiology and Developmental Biology*
>
> *Brigham Young University*
>
> *LSB 3050*
>
> *Provo, UT 84602*
>
>
> ------------------------------
> *From:* Giacomo Fiorin <giacomo.fiorin_at_gmail.com>
> *Sent:* Saturday, July 21, 2018 12:44:52 PM
> *To:* NAMD list; McGuire, Kelly
> *Cc:* Jérôme Hénin
> *Subject:* Re: namd-l: ABF Questions
>
> This is a clear sign that you need to include other degrees of freedom
> (i.e. other collective variables), which are limited by kinetic barriers
> and need to be biased as well. For a transmembrane channel pore, these can
> be movements of groups (side chains, or even entire protein domains) that
> need to get out of the way of the ligand before it can pass. This movement
> is a rare event, but likely much lower than 50 kcal/mol: that number is
> corrupted because the orthogonal degrees of freedom do not follow a
> canonical distribution in your simulation, but ABF (and all other free
> energy methods, really) assumes that they do.
>
> The solution is highly specific to the system of study: check the
> literature for this and other similar channels, to see whether what you
> plan can be done, and how did others treat it. You may have to invest
> significant effort or computer time, or restrict the investigation to a
> region of the PMF that can be sampled reliably.
>
> Giacomo
>
> On Sat, Jul 21, 2018 at 10:54 AM McGuire, Kelly <mcg05004_at_byui.edu> wrote:
>
> Thanks for the response! I'm trying to figure out why the free energy
> profile of my ion channel and ligand has about 50 kcal/mol higher energy
> with the ligand exiting into bulk water on the C-terminus side vs exiting
> the N-terminus side into bulk water. They should be about equal energy,
> right? That doesn't seem correct...
>
>
> *Kelly L. McGuire*
>
> *PhD Scholar*
>
> *Biophysics*
>
> *Department of Physiology and Developmental Biology*
>
> *Brigham Young University*
>
> *LSB 3050*
>
> *Provo, UT 84602*
>
>
> ------------------------------
> *From:* Jérôme Hénin <jerome.henin_at_ibpc.fr>
> *Sent:* Saturday, July 21, 2018 5:03:39 AM
> *To:* Namd Mailing List; McGuire, Kelly
> *Subject:* Re: namd-l: ABF Questions
>
>
>
> On Sat, 21 Jul 2018 at 04:17, McGuire, Kelly <mcg05004_at_byui.edu> wrote:
>
> Two questions about ABF:
>
>
> 1) Right now I have the center of mass of my ligand as the main atoms and
> the center of mass of my protein as the ref atoms (i.e. all of the CA C N O
> atom types of the backbone). But if for my ref atoms I just the center of
> mass of one of the side chains near the middle of the protein, could this
> change my PMF results drastically? How critical is the choice of ref atoms
> within the protein for ABF with a ligand in the protein?
>
>
> That is actually a question about your protein, not about ABF.
>
>
> (Is there a limit of atom numbers that can be used for the ref atoms?)
>
> There is no limit. Large numbers will incur a performance penalty, but a
> few dozens or even a couple hundred should be fine - especially with a
> center-of-mass-based coordinate.
>
>
> 2) The lowerboundarywall and upperboundarywall are used essentially to
> keep say a ligand within the specified window, but how much does the choice
> of value itself have an effect on the PMF result? Would 10 kcal/mol vs 100
> kcal/mol for these values have a large effect on my results?
>
>
> Most likely not. The only case that could be a problem is a PMF that is
> steeply decreasing towards the boundary, but it is more typically
> increasing.
>
> Jerome
>
>
>
> --
> Giacomo Fiorin
> Associate Professor of Research, Temple University, Philadelphia, PA
> Contractor, National Institutes of Health, Bethesda, MD
> http://goo.gl/Q3TBQU
> https://github.com/giacomofiorin
>

-- 
Giacomo Fiorin
Associate Professor of Research, Temple University, Philadelphia, PA
Contractor, National Institutes of Health, Bethesda, MD
http://goo.gl/Q3TBQU
https://github.com/giacomofiorin

This archive was generated by hypermail 2.1.6 : Sat Sep 14 2019 - 23:19:43 CDT