SSE Backbone Hydrogen Bonding

Determine backbone C=O to N-H hydrogen bonding in secondary structure elements.

sse.bridges(sse, type="helix", hbond=TRUE, energy.cut=-1.0)

Arguments

sse

an sse object as obtained with dssp.
type

character string specifying ‘helix’ or ‘sheet’.
hbond

use hbond records in the dssp output.
energy.cut

cutoff for the dssp hbond energy.

Details

Simple functionality to parse the ‘BP’ and ‘hbond’ records of the DSSP output.

Requires input from function dssp with arguments resno=FALSE and full=TRUE.

Value

Returns a numeric matrix of two columns containing the residue ids of the paired residues.

References

Grant, B.J. et al. (2006) Bioinformatics 22, 2695--2696.

Author

Lars Skjaerven

See also

read.pdb, dssp

Examples

if (FALSE) {
# Read a PDB file
pdb <- read.pdb( system.file("examples/1hel.pdb", package="bio3d") )
sse <- dssp(pdb, resno=FALSE, full=TRUE)

sse.bridges(sse, type="helix")
}