Re: Hi: Re: [gmx-users] Hello: Forcefield and potentials for DNA-zinc finger complex

2009-05-26 Thread Tsjerk Wassenaar
Hi Nehme,

> Also, in our zinc finger models, the zinc plays a structural role and it is
> not implicated in DNA recognition. Furthermore, I will read your paper and
> the references. I looked in the Literature and from NMR studies/X-ray and MD
> done on zinc fingers containing a zinc ion coordinated to 4 cysteines. It
> seems that these cysteines are deprotonated and the remaining negative
> charges are stabilized by nearby positively charged side chains of arginines
> or lysines.

It's hard, although not impossible, to assess the protonation state of
such a center experimentally. Mind that crystal structures and NMR
structures also use force fields and predefined topologies to model
the structure. That means that you will only get out what you put in.
If you set some group to be deprotonated while it should be
protonated, this may be hard or even impossible to track down. You may
still get a reasonable fit to the experimental data. For such issues
Xray/NMR are not the best means. Note that I'm not saying you're
wrong, but you should justify your choices.

> I tried to use CHARMM forcefield, but every time I had the 4 cysteines
> protonated, the coordination bonds removed and the parameters for zinc not
> assigned.
> How can I solve this issue?
>

In Gromacs you can use the file specbond.dat to define special bonds.
Have a look at http://wiki.gromacs.org/index.php/specbond.dat and
search the archives for more information.

Cheers,

Tsjerk


-- 
Tsjerk A. Wassenaar, Ph.D.
Junior UD (post-doc)
Biomolecular NMR, Bijvoet Center
Utrecht University
Padualaan 8
3584 CH Utrecht
The Netherlands
P: +31-30-2539931
F: +31-30-2537623
___
gmx-users mailing listgmx-users@gromacs.org
http://www.gromacs.org/mailman/listinfo/gmx-users
Please search the archive at http://www.gromacs.org/search before posting!
Please don't post (un)subscribe requests to the list. Use the 
www interface or send it to gmx-users-requ...@gromacs.org.
Can't post? Read http://www.gromacs.org/mailing_lists/users.php


Re: Hi: Re: [gmx-users] Hello: Forcefield and potentials for DNA-zinc finger complex

2009-05-20 Thread Tsjerk Wassenaar
Hi Nehme,

I've did simulations on TRAIL, which also contains a zinc-finger
domain, involving three cysteines and a chloride ion. But it's not so
simple. Luckily the exact parameters seemed not that vital, as the
zinc apparently functions to keep the three subunits together and was
in any case not directly related to the properties we were interested
in. Still, I spent some time parameterizing it. For that I did quite a
search of the available literature. You'll find some references in our
2008 Proteins paper. One of the most striking things is the still
unresolved issue of the protonation state of zinc finger domains.
Especially in the acidic environment around DNA, this may become an
issue. It is why I have stayed away from zinc-finger domains lately,
as they'd cost me too much work properly parameterizing them, even
though our lab is doing experimental work on them and it would be
interesting to complement that with MD. That's also why I mentioned to
be interested in such parameters. Not in what you can readily
calculate with QM or what you can get from automatic topology builders
or even from literature, including my own. That doesn't solve the
issue of the protonation state or overall charge of the complex. I
won't trust any report on zinc-finger domain MD, unless it is shown to
properly reproduce experimental properties, including what can be
determined in relation to the protonation state.

Cheers,

Tsjerk


On Wed, May 20, 2009 at 10:40 AM, Nehme El Hachem  wrote:
> Hi Tsjerk,
>
> I thank you for your help. Indeed, I am very new to GROMACS ( I have a Linux
> Red hat 9.0 and I do have some issues installing the software).
>
> If you are interested, please find a link:
> http://www.epress.com/w3jbio/vol2/harris2/harris.htm
> Molecular Dynamics Simulation Of The Estrogen Receptor Protein In Complex
> With A Non-Consensus Estrogen Response Element DNA Sequence In A 10 Angstrom
> Water Layer.
>
> I am trying to do some molecular dynamics on the oestrogen receptor DNA
> binding domain (PDB code: 1HCQ). I have InsightII with Discover3 and CHARMm
> but I am confused with the type of forcefield that I need. Furthermore, I
> have to add counterions to neutralize the system, so it seems complicated.
> I wanted to try GROMACS, because it has many forcfields (ex: ffAmber) for
> modeling this type of system.
>
> I appreciate your help,
>
> Nehme.
>
>



-- 
Tsjerk A. Wassenaar, Ph.D.
Junior UD (post-doc)
Biomolecular NMR, Bijvoet Center
Utrecht University
Padualaan 8
3584 CH Utrecht
The Netherlands
P: +31-30-2539931
F: +31-30-2537623
___
gmx-users mailing listgmx-users@gromacs.org
http://www.gromacs.org/mailman/listinfo/gmx-users
Please search the archive at http://www.gromacs.org/search before posting!
Please don't post (un)subscribe requests to the list. Use the 
www interface or send it to gmx-users-requ...@gromacs.org.
Can't post? Read http://www.gromacs.org/mailing_lists/users.php