Re: alchemical fep amber topology

From: Chris Harrison (char_at_ks.uiuc.edu)
Date: Sat Mar 14 2009 - 08:33:48 CDT

As far as I'm aware, unless dual topology has been recently implemented, the AMBER MD code uses single topology and cannot utilize dual topology approaches. So naturally, as pointed out by Jerome and Floris, there is no corresponding way to define dual topologies in the Amber file formats and therefore if you are using Amber format files of the topology and initial coordinates then NAMD cannot do dual topology FEPs (for that format).

However, there are available converted topology and parameter files for
the Amber forcefield written in the CHARMM format, so with work it is in
principle possible to do dual topology FEPs with the Amber forcefield,
just not with the Amber file formats.

C.

--
Chris Harrison, Ph.D.
Theoretical and Computational Biophysics Group
NIH Resource for Macromolecular Modeling and Bioinformatics
Beckman Institute for Advanced Science and Technology
University of Illinois, 405 N. Mathews Ave., Urbana, IL 61801
char_at_ks.uiuc.edu                            Voice: 217-244-1733
http://www.ks.uiuc.edu/~char               Fax: 217-244-6078
Floris Buelens <floris_buelens_at_yahoo.com> writes:
> Date: Sat, 14 Mar 2009 03:53:37 -0700 (PDT)
> From: Floris Buelens <floris_buelens_at_yahoo.com>
> To: George Madalin Giambasu <giambasu_at_gmail.com>, namd-l_at_ks.uiuc.edu
> Subject: Re: namd-l: alchemical fep amber topology
> Return-Path: char_at_halifax.ks.uiuc.edu
> Message-ID: <98263.42498.qm_at_web110803.mail.gq1.yahoo.com>
> X-Spam-Status: No, score=-2.6 required=5.0 tests=AWL,BAYES_00 
> 	autolearn=unavailable version=3.1.7-0+tcb1
> 
> 
> Hi George,
> 
> In principle this is absolutely fine in both released and CVS versions of NAMD - just do everything according to the tutorial, but with "amber on" and "parmfile" instead of a PSF file. However, there is as yet no equivalent to 'alchemify' for parm7 files. This is only an issue if you're working with dual topology alchemical groups that are bonded to the rest of the system (e.g. a dual topology residue definition) - for groups that aren't bonded to the rest of the system everything should work as you would expect without modifying the parm7.
> 
> Floris
> 
> 
> 
> 
> ----- Original Message ----
> From: George Madalin Giambasu <giambasu_at_gmail.com>
> To: namd-l_at_ks.uiuc.edu
> Sent: Friday, 13 March, 2009 18:55:22
> Subject: namd-l: alchemical fep amber topology
> 
> Hi,
> 
> I was wondering if there is a way to run alchemical fep simulations
> using Amber topologies (*.parm7 files) using the cvs version of namd?
> 
> I found a post by Jerome Henin 
> (http://www.ks.uiuc.edu/Research/namd/mailing_list/namd-l/1348.html) on
> the list in which it was stated that the fep setup "does not need CHARMM
> anymore", however the Alchemical FEP tutorial states that we have to use
> psf files in conjunction the alchemify program.
> 
> Thanks in advance for any  suggestions or updates!
> 
> 
> George
> 
> 
>       
> 

This archive was generated by hypermail 2.1.6 : Wed Feb 29 2012 - 15:52:29 CST