Re: Asymmetric concentration gradient

From: Gumbart, JC (gumbart_at_physics.gatech.edu)
Date: Tue Aug 24 2021 - 00:11:42 CDT

You can give solvate a minmax to only solvate a region above or below the membrane mid-plane (make sure to give them different segnames for both the water and ions). Then a few lines using psfgen (readpsf ___ pdb ___) will combine them. There’s no script though, but if you play around with the options, you should quickly get a sense of what to do.

Best,
JC

On Aug 24, 2021, at 12:10 AM, Ropon-Palacios G. <biodano.geo_at_gmail.com<mailto:biodano.geo_at_gmail.com>> wrote:

Dear Prof. Gumbart,

My question was on that side, how to do the solvation up and down independently, could you give me an example of the script, to be able to do it please. I find it a bit unfamiliar as this is my first time doing this.

From: "Gumbart, JC" <gumbart_at_physics.gatech.edu<mailto:gumbart_at_physics.gatech.edu>>
Date: Monday, August 23, 2021 at 11:05 PM
To: NAMD list <namd-l_at_ks.uiuc.edu<mailto:namd-l_at_ks.uiuc.edu>>, "Ropon-Palacios G." <biodano.geo_at_gmail.com<mailto:biodano.geo_at_gmail.com>>
Subject: Re: namd-l: Asymmetric concentration gradient

Hi Geo,

It’s been a long time since I ran this. But one way I can think of would be to solvate and ionize above and below separately and then combine them with psfgen. I attach an example tclBC script and NAMD input file for running the actual simulation once you’ve built the system.

Best,
JC



On Aug 23, 2021, at 10:08 PM, Ropon-Palacios G. <biodano.geo_at_gmail.com<mailto:biodano.geo_at_gmail.com>> wrote:

Dear users,

I am trying to simulate a membrane protein system where it has a concentration difference on one side as well as the other, similar to the work of Professor JC Gumbart “Molecular dynamics simulations of membrane proteins under asymmetric ionic concentrations”

My question is, how do I ionize or how do I place the ions on both sides with a difference in concentration? Could you please share a script to do this?.

Best,

Geo.

This archive was generated by hypermail 2.1.6 : Fri Dec 31 2021 - 23:17:11 CST