Re: Need help preventing hydration of the lipid-protein interface in my system

From: Chris Harrison (charris5_at_gmail.com)
Date: Thu May 03 2012 - 18:31:50 CDT

If I remember correctly, these are waters against which no forces will be
applied by the script. Be cautioned though, it's been sometime since I've
looked at that full script. Use cases for when you wouldn't want forces
applied would include crystallographic waters or other structural waters
that were of importance in the model and their positions from the original
crystal/setup/etc hopefully retained during equilibration.

Best,
Chris

On Thu, May 3, 2012 at 5:48 AM, Syed Kashif Zafar <skashifzafar_at_yahoo.com>wrote:

> Dear Users
>
> Hope this mail will find you well. Dear user i need your favor. I am
> reseearch fellow in HEJ RIC ICCBS pakistan.
> i am working on bilayer project. i am facing same problem when i am trying
> to equilibrate bilayer system with water as same faced by maggie in
> 2010 which are mentioned below
>
>
> When I do this, water molecules move into the hydrophobic part of the
> membrane. The tutorial uses the Tcl script named keep_water_out.tcl to
> prevent this hydration of the membrane-protein interface. I would like to
> modify this script for my protein, but am unsure where to begin. I have
> looked into the “User-defined Forces in NAMD” tutorial but I am still
> stuck, being I am not familiar with Tcl scripts.
> The specific section of the Tcl script that I believe I need to modify is
> under the “user definitions” section, where it lists a number of water
> molecules after “set watexcept.” How can I know which water molecules in my
> system need to be defined here?
> Below is an excerpt of the Tcl script.
> Thanks in advance for your help!
>
>
>
> #!/usr/local/bin/namd2
> # keep water from entering gaps between protein and lipids
> # need membrane be in the XY plane!
> # (C) A. Aksimentiev (alek_at_ks.uiuc.edu<alek_at_ks.uiuc.edu?Subject=Re:%20Need%20help%20preventing%20hydration%20of%20the%20lipid-protein%20interface%20in%20my%20system>)
>
> # This script will prevent water molecules in your namd simulations
> # from entering into the hydrophobic part of the membrane. This can be
> very
> # useful for initial equilibration of a protein/lipid/water system.
> # It works with recent NAMD binaries (version 2.4 and higher), which
> # has clearconfig command.
> # Brief description:
> # At initialization, the arrays of namd atom indices of all water oxygens
> # and of all lipid c21 carbons are created. They will be used later to
> # determine which water molecules to push out of the membrane and to
> compute
> # membrane dimensions. At every simulation step the pushing forces are
> applied
> # to a minor number of atoms, which doesn't affect NAMD performance. The
> # membrane dimensions are recalculated every N steps (regulated by the
> global
> # parameter $lipidCheckFreq) as well as the indices of the
> # water molecules that will be pushed during the next N steps (regulated
> by
> # the global parameter $waterCheckFreq). The membrane should be in XY
> plane.
> # To use this script one has to add the following lines to the namd
> # configuration script (without "#" symbol):
> #
> #
> # set waterCheckFreq 100
> # set lipidCheckFreq 100
> # set allatompdb a1b2c5_init.pdb
> # tclForces on
> # tclForcesScript keep_water_out.tcl
> #
> # # lipidCheckFreq must be an integer multiple of waterCheckFreq .
> ###################################################################
> # user definitions begin
> ###################################################################
> print "Starting Tcl forces"
> # define force constant per group (Kcal/(mol*A))
> set fconst 0.1
> set pressure 1.0
> set fconstUp [expr $fconst*$pressure]
> set fconstDown [expr -$fconst*$pressure]
> # define water and lipid names
> set watResName "TIP3"
> set watAtomName "OH2"
> set lipResName "POPC"
> set lipAtomName "C1 "
> # define exceptions
> set watexcept {30621 30624 30627 30630 30633 30642 30648 30654 30669 30678
> 30687 30690 30711 30717 30720 30744 30747 30756 30759 30762 30768 30801
> 30807 30816 30828 30843 30846 30852 30864 30888 30906 30927 30933 30936
> 30957 30960 30972 32546 32549 32552 32555 32558 32567 32579 32594 32603
> 32612 32615 32636 32642 32645 32669 32672 32681 32684 32687 32693 32726
> 32732 32741 32753 32768 32771 32777 32789 32813 32831 32852 32858 32861
> 32882 32885 32897 34471 34474 34477 34480 34483 34492 34498 34504 34519
> 34528 34537 34540 34561 34567 34570 34594 34597 34606 34609 34612 34618
> 34651 34657 34666 34678 34684 34693 34696 34702 34714 34738 34756 34777
> 34783 34786 34798 34807 34810 34822 36396 36399 36402 36405 36408 36417
> 36429 36444 36453 36462 36465 36486 36492 36495 36519 36522 36531 36534
> 36537 36543 36576 36582 36591 36603 36618 36621 36627 36639 36663 36681
> 36702 36708 36711 36732 36735 36747 36754 36757 36760}
> ###################################################################
> # user definitions stop here!
>
>
>
>
> Any support will be highly appreciable for me.
> advance in thanks
>
> regards
> kashif
>

-- 
Chris Harrison, Ph.D.
NIH Center for Macromolecular Modeling and Bioinformatics
Theoretical and Computational Biophysics Group
Beckman Institute for Advanced Science and Technology
University of Illinois, 405 N. Mathews Ave., Urbana, IL 61801
http://www.ks.uiuc.edu/Research/namd       Voice: 773-570-6078
http://www.ks.uiuc.edu/~char                          Fax: 217-244-6078

This archive was generated by hypermail 2.1.6 : Mon Dec 31 2012 - 23:21:29 CST