# Large RMSD Differences of Homologous Proteins
**URL:** https://gromacs.bioexcel.eu/t/large-rmsd-differences-of-homologous-proteins/931
**Category:** User discussions
**Created:** [October 8, 2020, 12:08am UTC](https://gromacs.bioexcel.eu/t/large-rmsd-differences-of-homologous-proteins/931 "2020-10-08T00:08:30Z")
**Posts on this page:** 11
**Page:** 2
### Author:  [@joshhayden](https://gromacs.bioexcel.eu/u/joshhayden)
#### Post date: [October 20, 2020, 1:19pm UTC](https://gromacs.bioexcel.eu/t/large-rmsd-differences-of-homologous-proteins/931/21 "2020-10-20T13:19:51Z")
I made the new .tpr file using gmx grompp, just as I did before running the final md simulation, but using the fixed.gro in the -c option. Was this correct?
Are there any other possible reasons other than the protein starting uncentered to explain the unusual RMSD?
---
### Author:  [@jalemkul](https://gromacs.bioexcel.eu/u/jalemkul)
#### Post date: [October 21, 2020, 1:12pm UTC](https://gromacs.bioexcel.eu/t/large-rmsd-differences-of-homologous-proteins/931/22 "2020-10-21T13:12:49Z")
> [@joshhayden](#):
>
> I made the new .tpr file using gmx grompp, just as I did before running the final md simulation, but using the fixed.gro in the -c option. Was this correct?
Yes.
> [@joshhayden](#):
>
> Are there any other possible reasons other than the protein starting uncentered to explain the unusual RMSD?
Can you provide your input files somewhere, and the first few frames of your trajectory (in any format)? I don’t want to guess wildly about anything and I still have no idea how one would obtain such large values of RMSD immediately at the outset of the simulation.
---
### Author:  [@joshhayden](https://gromacs.bioexcel.eu/u/joshhayden)
#### Post date: [October 23, 2020, 12:53pm UTC](https://gromacs.bioexcel.eu/t/large-rmsd-differences-of-homologous-proteins/931/23 "2020-10-23T12:53:17Z")
Would there be a good way to transfer my files to you?
---
### Author:  [@jalemkul](https://gromacs.bioexcel.eu/u/jalemkul)
#### Post date: [October 23, 2020, 5:42pm UTC](https://gromacs.bioexcel.eu/t/large-rmsd-differences-of-homologous-proteins/931/24 "2020-10-23T17:42:21Z")
If they can’t simply be attached to your post here, upload to Google Drive and provide the link. Make sure sharing options allow for anyone with the link to see the contents.
---
### Author:  [@joshhayden](https://gromacs.bioexcel.eu/u/joshhayden)
#### Post date: [October 24, 2020, 1:50am UTC](https://gromacs.bioexcel.eu/t/large-rmsd-differences-of-homologous-proteins/931/25 "2020-10-24T01:50:39Z")
In a Google Drive folder, I have uploaded all files used to calculate the rmsd, as well as the fixed.gro file made, the previous md\_0\_150.tpr, and an MPG of the filtered trajectory visualized with VMD.
The link to the Google Drive folder is here:
> **[L93K with D-Tartrate - Google Drive](https://drive.google.com/drive/folders/1C74uklkVsfxZEadB0V4DZgXZAi4IO766?usp=sharing)**
---
### Author:  [@jalemkul](https://gromacs.bioexcel.eu/u/jalemkul)
#### Post date: [October 28, 2020, 2:00pm UTC](https://gromacs.bioexcel.eu/t/large-rmsd-differences-of-homologous-proteins/931/26 "2020-10-28T14:00:35Z")
Please reduce the file size of the trajectory by only including the first few frames. I’m pretty sure I know what’s going on but I want to test something. 5-10 frames should be sufficient to diagnose.
---
### Author:  [@joshhayden](https://gromacs.bioexcel.eu/u/joshhayden)
#### Post date: [October 28, 2020, 7:14pm UTC](https://gromacs.bioexcel.eu/t/large-rmsd-differences-of-homologous-proteins/931/27 "2020-10-28T19:14:43Z")
I have uploaded a .xtc file of the first 10 frames to the same Google Drive I shared!
---
### Author:  [@jalemkul](https://gromacs.bioexcel.eu/u/jalemkul)
#### Post date: [October 29, 2020, 6:09pm UTC](https://gromacs.bioexcel.eu/t/large-rmsd-differences-of-homologous-proteins/931/28 "2020-10-29T18:09:45Z")
I can’t reproduce the weird numbers.
`gmx rms -s fixed.gro -f md_0_150_chainA_frames_0_10.xtc`
gives:
```
# This file was created Thu Oct 29 14:08:01 2020
# Created by:
# :-) GROMACS - gmx rms, 2020.3 (-:
#
# Executable: /usr/local/gromacs/current/bin/gmx
# Data prefix: /usr/local/gromacs/current
# Working dir: /Users/jalemkul/Downloads
# Command line:
# gmx rms -s fixed.gro -f md_0_150_chainA_frames_0_10.xtc
# gmx rms is part of G R O M A C S:
#
# Giant Rising Ordinary Mutants for A Clerical Setup
#
@ title "RMSD"
@ xaxis label "Time (ps)"
@ yaxis label "RMSD (nm)"
@TYPE xy
@ subtitle "Backbone after lsq fit to Backbone"
0.0000000 0.0004864
10.0000000 0.0865524
20.0000000 0.0958802
30.0000000 0.0962045
40.0000000 0.0943499
50.0000000 0.0967933
60.0000000 0.0874779
70.0000000 0.0879277
80.0000000 0.0977267
90.0000000 0.1061396
100.0000000 0.1017297
```
These numbers are sensible and don’t look like what you have. Note that I’m selecting the protein backbone for the calculation. Did you do something different?
---
### Author:  [@joshhayden](https://gromacs.bioexcel.eu/u/joshhayden)
#### Post date: [October 30, 2020, 3:54pm UTC](https://gromacs.bioexcel.eu/t/large-rmsd-differences-of-homologous-proteins/931/29 "2020-10-30T15:54:53Z")
Thank you very much for trying to calculate it yourself, Dr. Lemkul. It is greatly appreciated!
I did select the protein backbone for the calculation. The only difference was that I used the centered .tpr file that is also in the Google Drive folder as opposed to the fixed.gro file that you used for the reference structure. I’m assuming that this would make the difference.
Did you manually account for the Mg ion mass in order to use the fixed.gro file?
---
### Author:  [@jalemkul](https://gromacs.bioexcel.eu/u/jalemkul)
#### Post date: [November 2, 2020, 3:02pm UTC](https://gromacs.bioexcel.eu/t/large-rmsd-differences-of-homologous-proteins/931/30 "2020-11-02T15:02:20Z")
I changed MG -\> Mg in fixed.gro since atommass.dat works on a case-sensitive basis.
---
### Author:  [@joshhayden](https://gromacs.bioexcel.eu/u/joshhayden)
#### Post date: [November 4, 2020, 9:24pm UTC](https://gromacs.bioexcel.eu/t/large-rmsd-differences-of-homologous-proteins/931/31 "2020-11-04T21:24:41Z")
That appears to have done the trick when calculating for all the frames as well! The RMSD appears to hover around 0.25 nm.
Thank you so very much for your help!
[Previous page](https://gromacs.bioexcel.eu/t/large-rmsd-differences-of-homologous-proteins/931.md?page=1)