From: Victor Kwan (vkwan8_at_uwo.ca)
Date: Thu Feb 27 2020 - 16:26:39 CST
In VMD: Load the last frame of your DCD and the modified pdb. File -> Save Coordinates -> First 0 Last 0, Save as a new pdb.
You'll get a pdb file with modified resid and the coordinate of the system after equilibrium. You might have to fix the O/B column afterward.
On Thu, Feb 27, 2020 at 5:03 PM Francesco Pietra <chiendarret_at_gmail.com<mailto:chiendarret_at_gmail.com>> wrote:
setting up QM-MM needs a psf and a related pdb
On Thu, Feb 27, 2020 at 9:45 PM Victor Kwan <vkwan8_at_uwo.ca<mailto:vkwan8_at_uwo.ca>> wrote:
DCD files does not contain topological information (resid, etc), only the atom number. It does not matter if you load a psf containing resid 0 or resid 999 after you load the dcd.
On Thu, Feb 27, 2020 at 3:28 PM Francesco Pietra <chiendarret_at_gmail.com<mailto:chiendarret_at_gmail.com>> wrote:
Hello Victor, of course, as implied in my "renumbering". 3 ,5 etc are fictitious resid ID, which in reality were not in consecutive order. i am trying to save the MD work that i did with resid 0.
On Thu, Feb 27, 2020, 21:02 Victor Kwan <vkwan8_at_uwo.ca<mailto:vkwan8_at_uwo.ca>> wrote:
ResID does not have to be in numerical order. I have attached a short proof of concept of 216 SWM4-NDP water that starts with resid 999. You can easily change the resid from 0 to 999 in your pdb and psf files.
On Thu, Feb 27, 2020 at 1:40 PM Francesco Pietra <chiendarret_at_gmail.com<mailto:chiendarret_at_gmail.com>> wrote:
My question is whether there is a way to circumvent the inability of <prepare.qm.region.tcl> (vmd text mode) with nucleotides and a small organic ligand, to recognize residue number 0 for the organic ligand
set seltext "(nucleic and resid 0 3 5 7 8 9)"
the generated .QMonly.psf contains all except what relates to resid 0.
I understand that I could renumber the starting pdb and rebuild psf/pdb, but I have already equilibrated the system with resid 0. Segname is not useful here because it does not discriminate the nucleotides.
Thanks for advice
This archive was generated by hypermail 2.1.6 : Fri Dec 31 2021 - 23:17:08 CST