# compute.dssp()

### Description

Computes and returns a DSSP file

### Usage

compute.dssp(pdb, destfile = './')

### Arguments

pdb is either a PDB id, or the path to a pdb file.

destfile a character string with the path where the DSSP file is going to be saved.

### Value

An online computed dsssp file that is saved at the indicated location (the working directory by default).

### References

parse.dssp(), mkdssp() and acc.dssp()

### Details

The ptm package contains a number of ancillary functions that deal with Protein Data Bank (PDB) files. These functions may be useful when structural 3D data need to be analyzed. The mentioned functions are:

The DSSP (Define Secondary Structure of Proteins) algorithm allows to assign secondary structure to the amino acids of a protein using the atomic coordinates of the protein (a PDB file).

DSSP, based on the identification of intra-backbone hydrogen bonds of the protein can identify eight types of secondary structure that can be grouped in three main categories:

• Helices
G = 3-turn helix $3_{10}$ helix. Min length 3 residues.
H = 4-turn helix α helix. Minimum length 4 residues.
I = 5-turn helix π helix. Minimum length 5 residues.

• Strands
E = extended strand in parallel and/or anti-parallel β-sheet conformation. Min length 2 residues.
B = residue in isolated β-bridge (single pair β-sheet hydrogen bond formation)

• Loops
T = hydrogen bonded turn (3, 4 or 5 turn)
S = bend (the only non-hydrogen-bond based assignment).
C = coil (residues which are not in any of the above conformations).

Given a PDB file (or its 4 letter ID) there are two different ways to get the corresponding dssp file. The one used by the function compute.dssp() consist in taking advantage of the REST API provided by the Centre for Molecular and Biomolecular Informatics, to carry out the computation online (see Facilities that make the PDB data collection more powerful for a recent review).

For instance, let’s obtain the dssp file for the Dynein light chain 2.

compute.dssp(pdb = '2xqq')
## [1] "Work done!. See file at: ./2xqq.dssp"

We can now parse the obtained dssp file:

Dynein <- parse.dssp('./2xqq.dssp')