Previous Thread
Next Thread
Print Thread
Page 4 of 6 1 2 3 4 5 6
lennart #273 06/02/08 04:31 PM
Joined: Jul 2004
Posts: 165
B
Forum Member
Offline
Forum Member
B
Joined: Jul 2004
Posts: 165
perhaps it is a naive question.
Quote:

It is typically not necessary to define a cutoff for the angle.




Since there are many difference criteria have been used in the lipid simulation, is this criteria is also reasonable for the lipid system. I am not so sure becuse we are here mainly focused on the peptide/protein systems.
thank you very much

beginner #274 06/02/08 05:07 PM
Joined: Sep 2003
Posts: 8,659
Likes: 26
rmv Online Content
Forum Member
Online Content
Forum Member
Joined: Sep 2003
Posts: 8,659
Likes: 26
The forces between atoms are computed in CHARMM w/o regard to the molecule type; for H-bonds, the partial charges and VDW radii play a major role, esp. in the close approach distance.

If you are comparing to another simulation result, you should probably try to use the same distance and angle criteria.

Otherwise, I suggest trying the COOR HBOND defaults; we've done that in a publication to evaluate lipid:water and lipid:sugar H-bonds.

Making sure you've properly accounted for image interactions may be a more important issue, esp. if the peptide reaches the edge of the simulation cell during the simulation.

rmv #275 10/09/08 12:28 AM
Joined: Mar 2007
Posts: 23
D
Forum Member
Offline
Forum Member
D
Joined: Mar 2007
Posts: 23
Hi,
I used the code listed here for estimating the number of hydrogen bonds and lifetime. I am wondering is there any way to get the average distance between the donors and acceptors as well ?

Thanks a lot
Diwakar

diwakar #276 10/09/08 06:07 AM
Joined: Sep 2003
Posts: 4,883
Likes: 12
lennart Offline OP
Forum Member
OP Offline
Forum Member
Joined: Sep 2003
Posts: 4,883
Likes: 12
As of CHARMM version 32 the COOR HBONd command can also give the distance distribution, from which you should easily be able to compute the average. From corman.doc:
"Distance and lifetime histograms can be computed for all (putative) hydrogen
bonds encountered in the analysis"


Lennart Nilsson
Karolinska Institutet
Stockholm, Sweden
lennart #277 10/14/08 02:00 AM
Joined: Mar 2007
Posts: 23
D
Forum Member
Offline
Forum Member
D
Joined: Mar 2007
Posts: 23
Thanks a lot for the help.
I have a couple more questions about calculating the number and lifetime of hydrogen bonds.
1) Is it possible to exclude the intra-molecular hbonds?
I calculated the hbonds by selecting each molecule separately and then counting the hbonds of this molecule with other molecule. This process excludes the intra-molecular bonds but I was wondering if there is a better way of doing it.

2) How can I get the number of hbonds as a function of time?

Regards,
Diwakar

diwakar #278 10/14/08 05:12 AM
Joined: Sep 2003
Posts: 4,883
Likes: 12
lennart Offline OP
Forum Member
OP Offline
Forum Member
Joined: Sep 2003
Posts: 4,883
Likes: 12
1/ You did the right thing, there is no other way to do it. Using the two selections is not all that difficult, is it?
2/ Use a COOR HBOND (without keyword VERBose) in a CHARMM loop and write the number of hydrogen bonds (available as a CHARMM substitution variable, see subst.doc) to a file.


Lennart Nilsson
Karolinska Institutet
Stockholm, Sweden
lennart #279 11/14/08 11:41 PM
Joined: Mar 2007
Posts: 23
D
Forum Member
Offline
Forum Member
D
Joined: Mar 2007
Posts: 23
I am trying to get water water hydrogen bonds. I have added the donor list to TIP3 water model as described in the hbonds.inp. I get the same psf without any extra donors on regenerating the psf using this code:

read rtf card name top_all27.rtf
read para card name par_all27.prm

read psf card name @pdbfile@qual.psf
read coor card name @pdbfile@qual.crd

write psf card name @pdbfile@qual-water-donor.psf
write coor card name @pdbfile@qual-water-donor.crd

I want to know if anything more needs to be done to regenrate the psf?

Thanks
Diwakar

diwakar #280 11/15/08 12:20 AM
Joined: Sep 2003
Posts: 8,659
Likes: 26
rmv Online Content
Forum Member
Online Content
Forum Member
Joined: Sep 2003
Posts: 8,659
Likes: 26
Yes, you must use the GENERATE command instead of READ PSF

An alternative is the DONOR and ACCEPTOR commands; they seem to work as long as there only one of each in a given CHARMM input.

rmv #281 11/15/08 08:56 AM
Joined: Sep 2003
Posts: 4,883
Likes: 12
lennart Offline OP
Forum Member
OP Offline
Forum Member
Joined: Sep 2003
Posts: 4,883
Likes: 12
DONOr/ACCEptor ADD have some issues, and do not work properly in older versions of CHARMM. It is really very easy to regenerate the PSF so I strongly recommend that route; if your water segment is the last one all it takes is:
read rtf card name ... ! the new RTF w/ DONO, ACCE for TIP3
read param card name ..
read psf card name old.psf
delet atom sele segid wat
read seque tip3 345
gene wat noangle nodihe
write psf card name new.psf
* PSF for ... with DONO/ACCE for TIP3
*

If water residue numbering in the original PSF was not sequential one may change it to sequential with the JOIN command (struct.doc).


Lennart Nilsson
Karolinska Institutet
Stockholm, Sweden
lennart #282 11/15/08 04:53 PM
Joined: Sep 2003
Posts: 8,659
Likes: 26
rmv Online Content
Forum Member
Online Content
Forum Member
Joined: Sep 2003
Posts: 8,659
Likes: 26
A good point; you don't have to start over from the beginning of your model building. What's most important is that the number and ordering of atom types in the PSF does not change. Reading the old PSF and deleting the water segment alone won't work if additional segments were added after the water, such as ions. In that case, you'd need to also delete the segments that follow the water segment in the PSF, and then add them back via GENERATE, in the same order as was done initially.

Page 4 of 6 1 2 3 4 5 6

Moderated by  lennart, rmv 

Link Copied to Clipboard
Powered by UBB.threads™ PHP Forum Software 7.7.5
(Release build 20201027)
Responsive Width:

PHP: 7.3.31-1~deb10u5 Page Time: 0.014s Queries: 34 (0.011s) Memory: 0.7851 MB (Peak: 0.8628 MB) Data Comp: Off Server Time: 2023-12-03 01:00:13 UTC
Valid HTML 5 and Valid CSS