# How to read the dssp.dat file gmx 2023.1

**URL:** <https://gromacs.bioexcel.eu/t/how-to-read-the-dssp-dat-file-gmx-2023-1/9063>\
**Category:** User discussions\
**Created:** [May 5, 2024, 12:21pm UTC](https://gromacs.bioexcel.eu/t/how-to-read-the-dssp-dat-file-gmx-2023-1/9063 "2024-05-05T12:21:45Z")\
**Posts on this page:** 20\
**Page:** 1

<div class="post-metadata">

**Author:** ![chitlurikiran](https://avatars.discourse-cdn.com/v4/letter/c/7bcc69/32.png) [@chitlurikiran](https://gromacs.bioexcel.eu/u/chitlurikiran)\
**Post date:** [May 5, 2024, 12:21pm UTC](https://gromacs.bioexcel.eu/t/how-to-read-the-dssp-dat-file-gmx-2023-1/9063/1 "2024-05-05T12:21:45Z")

</div>

GROMACS version:  
GROMACS modification: Yes/No  
Here post your question

gmx dssp -s md\_0\_1.tpr -f md\_0\_1\_noPBC.xtc -o dssp\_1.dat -xvg xmgrace

unable to plot dssp.dat file

---

<div class="post-metadata">

**Author:** ![MagnusL](https://dub1.discourse-cdn.com/flex017/user_avatar/gromacs.bioexcel.eu/magnusl/32/2380_2.png) [@MagnusL](https://gromacs.bioexcel.eu/u/MagnusL)\
**Post date:** [May 6, 2024, 8:23am UTC](https://gromacs.bioexcel.eu/t/how-to-read-the-dssp-dat-file-gmx-2023-1/9063/2 "2024-05-06T08:23:33Z")

</div>

What is it you want to plot? Have you tried using the `-num` output option?

---

<div class="post-metadata">

**Author:** ![chitlurikiran](https://avatars.discourse-cdn.com/v4/letter/c/7bcc69/32.png) [@chitlurikiran](https://gromacs.bioexcel.eu/u/chitlurikiran)\
**Post date:** [May 6, 2024, 9:43am UTC](https://gromacs.bioexcel.eu/t/how-to-read-the-dssp-dat-file-gmx-2023-1/9063/3 "2024-05-06T09:43:39Z")

</div>

i want to do secondary structure analysis, but unfortunately the output is in .dat format and i am unable to use the .dat to plot

---

<div class="post-metadata">

**Author:** ![chitlurikiran](https://avatars.discourse-cdn.com/v4/letter/c/7bcc69/32.png) [@chitlurikiran](https://gromacs.bioexcel.eu/u/chitlurikiran)\
**Post date:** [May 6, 2024, 9:47am UTC](https://gromacs.bioexcel.eu/t/how-to-read-the-dssp-dat-file-gmx-2023-1/9063/4 "2024-05-06T09:47:55Z")

</div>

could you elaborate the command?

---

<div class="post-metadata">

**Author:** ![MagnusL](https://dub1.discourse-cdn.com/flex017/user_avatar/gromacs.bioexcel.eu/magnusl/32/2380_2.png) [@MagnusL](https://gromacs.bioexcel.eu/u/MagnusL)\
**Post date:** [May 6, 2024, 9:51am UTC](https://gromacs.bioexcel.eu/t/how-to-read-the-dssp-dat-file-gmx-2023-1/9063/5 "2024-05-06T09:51:16Z")

</div>

You can open the .dat file in any text editor. But the information is not meant for plotting.

If you are interested in statistics over time, which could be plotted, the `-num` option will give that, e.g.: `gmx dssp -s md_0_1.tpr -f md_0_1_noPBC.xtc -o dssp_1.dat -num num.xvg`.

---

<div class="post-metadata">

**Author:** ![chitlurikiran](https://avatars.discourse-cdn.com/v4/letter/c/7bcc69/32.png) [@chitlurikiran](https://gromacs.bioexcel.eu/u/chitlurikiran)\
**Post date:** [May 6, 2024, 9:57am UTC](https://gromacs.bioexcel.eu/t/how-to-read-the-dssp-dat-file-gmx-2023-1/9063/7 "2024-05-06T09:57:00Z")

</div>

gmx dssp -s md\_0\_1.tpr -f md\_0\_1\_noPBC.xtc -o dssp\_1.dat -num num.xvg

Error in user input:  
Invalid command-line options  
Unknown command-line option -num

it gave an error

---

<div class="post-metadata">

**Author:** ![chitlurikiran](https://avatars.discourse-cdn.com/v4/letter/c/7bcc69/32.png) [@chitlurikiran](https://gromacs.bioexcel.eu/u/chitlurikiran)\
**Post date:** [May 6, 2024, 10:00am UTC](https://gromacs.bioexcel.eu/t/how-to-read-the-dssp-dat-file-gmx-2023-1/9063/8 "2024-05-06T10:00:30Z")

</div>

is there any way to perform secondary structure analysis using 2023.1

---

<div class="post-metadata">

**Author:** ![MagnusL](https://dub1.discourse-cdn.com/flex017/user_avatar/gromacs.bioexcel.eu/magnusl/32/2380_2.png) [@MagnusL](https://gromacs.bioexcel.eu/u/MagnusL)\
**Post date:** [May 6, 2024, 10:02am UTC](https://gromacs.bioexcel.eu/t/how-to-read-the-dssp-dat-file-gmx-2023-1/9063/9 "2024-05-06T10:02:19Z")

</div>

Sorry, I didn’t notice you are using GROMACS 2023. The `dssp` tool was changed in GROMACS 2024. You cannot get time based statistics (the `num` option) in GROMACS 2023.

---

<div class="post-metadata">

**Author:** ![MagnusL](https://dub1.discourse-cdn.com/flex017/user_avatar/gromacs.bioexcel.eu/magnusl/32/2380_2.png) [@MagnusL](https://gromacs.bioexcel.eu/u/MagnusL)\
**Post date:** [May 6, 2024, 10:04am UTC](https://gromacs.bioexcel.eu/t/how-to-read-the-dssp-dat-file-gmx-2023-1/9063/10 "2024-05-06T10:04:18Z")

</div>

You can still do secondary structure analysis. But then I’m not sure what you expect to plot.

---

<div class="post-metadata">

**Author:** ![chitlurikiran](https://avatars.discourse-cdn.com/v4/letter/c/7bcc69/32.png) [@chitlurikiran](https://gromacs.bioexcel.eu/u/chitlurikiran)\
**Post date:** [May 6, 2024, 10:05am UTC](https://gromacs.bioexcel.eu/t/how-to-read-the-dssp-dat-file-gmx-2023-1/9063/11 "2024-05-06T10:05:49Z")

</div>

![scount](https://europe1.discourse-cdn.com/flex017/uploads/bioexcel1/original/2X/b/bc64247fdf4037785132d272485b0c9a674a07b7.png)

i want to plot an image as attached… Can i use the present version .xtc and .tpr file in gromacs 2024 to plot?

---

<div class="post-metadata">

**Author:** ![MagnusL](https://dub1.discourse-cdn.com/flex017/user_avatar/gromacs.bioexcel.eu/magnusl/32/2380_2.png) [@MagnusL](https://gromacs.bioexcel.eu/u/MagnusL)\
**Post date:** [May 6, 2024, 10:06am UTC](https://gromacs.bioexcel.eu/t/how-to-read-the-dssp-dat-file-gmx-2023-1/9063/12 "2024-05-06T10:06:57Z")

</div>

That kind of time-based analysis should be possible in gromacs 2024, yes. And you can use the input from your current simulations.

---

<div class="post-metadata">

**Author:** ![chitlurikiran](https://avatars.discourse-cdn.com/v4/letter/c/7bcc69/32.png) [@chitlurikiran](https://gromacs.bioexcel.eu/u/chitlurikiran)\
**Post date:** [May 6, 2024, 1:07pm UTC](https://gromacs.bioexcel.eu/t/how-to-read-the-dssp-dat-file-gmx-2023-1/9063/13 "2024-05-06T13:07:40Z")

</div>

WIll try with gromacs 2024, thanks @MagnusL for your prompt reply…

---

<div class="post-metadata">

**Author:** ![chitlurikiran](https://avatars.discourse-cdn.com/v4/letter/c/7bcc69/32.png) [@chitlurikiran](https://gromacs.bioexcel.eu/u/chitlurikiran)\
**Post date:** [May 7, 2024, 3:42pm UTC](https://gromacs.bioexcel.eu/t/how-to-read-the-dssp-dat-file-gmx-2023-1/9063/14 "2024-05-07T15:42:33Z")

</div>

> [@chitlurikiran](#):
>
> gmx dssp -s md\_0\_1.tpr -f md\_0\_1\_noPBC.xtc -o dssp\_1.dat -num num.xvg

i had run the command using 2024.1 but generated only LOOPS

 ![Screenshot from 2024-05-08 02-25-33](https://europe1.discourse-cdn.com/flex017/uploads/bioexcel1/original/2X/d/d2ae3de169ceb9566eff93ad7e8e9fd853c12314.png)

how to generate time based analysis as previously shared image

---

<div class="post-metadata">

**Author:** ![MagnusL](https://dub1.discourse-cdn.com/flex017/user_avatar/gromacs.bioexcel.eu/magnusl/32/2380_2.png) [@MagnusL](https://gromacs.bioexcel.eu/u/MagnusL)\
**Post date:** [May 7, 2024, 4:10pm UTC](https://gromacs.bioexcel.eu/t/how-to-read-the-dssp-dat-file-gmx-2023-1/9063/15 "2024-05-07T16:10:03Z")

</div>

Have you checked the file contents in a text editor? Have you tried `xmgrace -nxy num.xvg` (if you are using `xmgrace`)? I think it only plots the first set by default.

---

<div class="post-metadata">

**Author:** ![chitlurikiran](https://avatars.discourse-cdn.com/v4/letter/c/7bcc69/32.png) [@chitlurikiran](https://gromacs.bioexcel.eu/u/chitlurikiran)\
**Post date:** [May 8, 2024, 2:57am UTC](https://gromacs.bioexcel.eu/t/how-to-read-the-dssp-dat-file-gmx-2023-1/9063/16 "2024-05-08T02:57:17Z")

</div>

It worked !!! thanks a lot @MagnusL you saved my day…

---

<div class="post-metadata">

**Author:** ![MagnusL](https://dub1.discourse-cdn.com/flex017/user_avatar/gromacs.bioexcel.eu/magnusl/32/2380_2.png) [@MagnusL](https://gromacs.bioexcel.eu/u/MagnusL)\
**Post date:** [May 8, 2024, 6:55am UTC](https://gromacs.bioexcel.eu/t/how-to-read-the-dssp-dat-file-gmx-2023-1/9063/17 "2024-05-08T06:55:57Z")

</div>

Great!

---

<div class="post-metadata">

**Author:** ![Dewald21](https://avatars.discourse-cdn.com/v4/letter/d/a5b964/32.png) [@Dewald21](https://gromacs.bioexcel.eu/u/Dewald21)\
**Post date:** [May 21, 2024, 9:39am UTC](https://gromacs.bioexcel.eu/t/how-to-read-the-dssp-dat-file-gmx-2023-1/9063/18 "2024-05-21T09:39:33Z")

</div>

Hi @chitlurikiran and @MagnusL, I just wanted to say thank you SO much for this post and the responses! I also ran my MD simulations with gmx 2023.2, and thanks to this post, I’ve figured out how to plot my dssp.dat files in xmgrace.  
Thanks so much!

---

<div class="post-metadata">

**Author:** ![Sounak\_Biswas](https://avatars.discourse-cdn.com/v4/letter/s/3ec8ea/32.png) [@Sounak\_Biswas](https://gromacs.bioexcel.eu/u/Sounak_Biswas)\
**Post date:** [May 21, 2024, 11:44am UTC](https://gromacs.bioexcel.eu/t/how-to-read-the-dssp-dat-file-gmx-2023-1/9063/19 "2024-05-21T11:44:26Z")

</div>

@Dewald21 i didn’t figure out this.Could you please help with this. I am using gromacs 2023.4.If you kindly share what command you have just run.

---

<div class="post-metadata">

**Author:** ![Dewald21](https://avatars.discourse-cdn.com/v4/letter/d/a5b964/32.png) [@Dewald21](https://gromacs.bioexcel.eu/u/Dewald21)\
**Post date:** [May 21, 2024, 12:33pm UTC](https://gromacs.bioexcel.eu/t/how-to-read-the-dssp-dat-file-gmx-2023-1/9063/20 "2024-05-21T12:33:55Z")

</div>

@Sounak_Biswas, so here’s what I did and what worked for me.  
I ran my MD simulations in gromacs 2023.2 (compiled on our HPC cluster without mpi) and then did the dssp with gromacs 2024.1 (compiled on our HPC cluster WITH mpi). But despite the version differences, my dssp analysis still seemed to work and generated the plot that I was looking for.

I ran the gmx dssp on our cluster with gromacs 2024.1 using the command:  
`mpirun gmx_mpi dssp -f md_fit_300_500ns.xtc -s md.tpr -o dssp_peptide_300_500ns.dat -num dssp_peptide_300_500ns.xvg -tu ns -n index.ndx -sel '17'`

- `mpirun gmx_mpi dssp`: the gmx dssp command for our cluster’s mpi-compiled gromacs 2024.1
- `-num dssp_peptide_300_500ns.xvg`: ensures that your results is output to an xvg file that can be viewed in xmgrace (that’s what I used to view my results as a plot over the 500ns of my simulation)
- `-tu ns`: works with gromacs 2024.1 (not with 2023.2) to show you simulation in ns rather than ps, for example
- `-n index.ndx -sel '17'`: specifies to include the index file, because I specifically wanted to view the secondary structure of the peptide in my complex, which is group 17 in my index file.

Then I copied the dssp\_peptide\_300\_500ns.xvg file back to my local machine and used `xmgrace -nxy dssp_peptide_300_500ns.xvg` to view the plot in xmgrace.  
I hope that helps and that you manage to get your working!

Unfortunately, my plot doesn’t look exactly like the one that @chitlurikiran showed in this thread, but I’m still tweaking it to get what I’m looking for. This is my first time doing secondary structure analysis (and using gmx dssp), so I’m still figuring out how it works and how I can use the different options to get what I’m looking for.

---

<div class="post-metadata">

**Author:** ![Sounak\_Biswas](https://avatars.discourse-cdn.com/v4/letter/s/3ec8ea/32.png) [@Sounak\_Biswas](https://gromacs.bioexcel.eu/u/Sounak_Biswas)\
**Post date:** [June 9, 2024, 10:30am UTC](https://gromacs.bioexcel.eu/t/how-to-read-the-dssp-dat-file-gmx-2023-1/9063/21 "2024-06-09T10:30:17Z")

</div>

@Dewald21 thank you so much, could you kindly share the images of your generated graphs.

[Next page](https://gromacs.bioexcel.eu/t/how-to-read-the-dssp-dat-file-gmx-2023-1/9063.md?page=2)
