I chose C-alpha as a group for root mean square calculation(rmsf), but why my rmsf calculation in gromacs did still count on atom-based, not residue-based?
The command I used: gmx rmsf -f md_noPBC.xtc -s md.tpr -o rmsf.xvg
Thanks in advanced
You need to use the flag -res, if not the default is measured by atom.
Cheers
I've just found a solution for this:
just need to add flag -res , so the command would be:
gmx rmsf -f md_noPBC.xtc -s md.tpr -o rmsf.xvg -res
then choose C-alpha as a group for rms calculation
You have to create an index:
$ gmx make_ndx -f md_0_10.tpr -o index_ca.ndx
Choose
> 1 & 3
22 Protein_&_C-alpha : 347 atoms
Then, calculate RMSF as follows:
gmx rmsf -s md_0_10.tpr -f md_0_10_center.xtc -n index_ca.ndx -o rmsf_protein.xvg
Thank you! M. Florencia Martini
For atomic level
gmx rmsf -f md_0_1.xtc -s md_0_1.tpr -o rmsf.xvg
and for residue
change .xtc and .tpr file name accordingly.
I used to do this solvation by doing these steps: 1. found the box volume information from gromacs by choosing the distance between the outer side of protein and the box randomly (I set it to 22...
25 August 2020 9,736 0 View
Hi, I did not know really much about xmgrace, so there's an unsolved problem I got: I wanted to open and plot the data within a dat file in xmgrace. I clicked Menu Data --> Import --> ASCII -->...
16 August 2020 2,459 2 View
Hi, I'm kinda new to computational biochemistry. so here are some few questions of mine: I made an amylase system solvated with NaCl 5M via PACKMOL, and the pdb file of solvated amylase was...
09 August 2020 2,606 2 View
We have observed that tube to tube sheet joint leaked in our boiler and needs to overcome same by knowing the root cause.
08 August 2024 3,161 0 View
I have face this problem anyone help me how to solve this issue ?which is below Fatal error: There are inconsistent shifts over periodic boundaries in a molecule type consisting of 78 atoms. The...
07 August 2024 2,598 1 View
How do soil microflora interact with plant roots and influence plant nutrition, health, and productivity?
06 August 2024 9,618 3 View
Hello, I am currently analyzing some phosphoproteomics data, but I have peptides with multiple phosphorylation sites or phosphorylations together with carbamidomethylation or oxidation. How can I...
04 August 2024 8,432 3 View
When we conduct linear regression, there are several assumptions. The assumption of normality is whether the residual errors are normally distributed, not whether a predictor is normal?
31 July 2024 6,164 3 View
How microorganisms are important for maintaining of healthy soil and biodiversity and microorganisms and plant roots contribute to soil formation?
31 July 2024 8,939 5 View
Dear All Gromacs User I would like to restart MD without using checkpint .cpt files? I heard there is a way for acheiveing this. Any input and sugegstions are most welcome. Thanking You Budheswar
29 July 2024 3,084 8 View
I wanted to know whether we can observe the synergistic/antagonistic/additive properties of combinations or mixtures of compounds through docking analysis. But during docking preparation any...
28 July 2024 7,413 6 View
If I want to calculate molecular dynamics (Gromacs), do I first need to optimize the geometry of the molecule? For example, I first draw the molecule in Avogadro, and then what? I guess I can't...
28 July 2024 673 3 View
Dear people,
24 July 2024 1,252 2 View