Description Usage Arguments Details Value Note Author(s) References See Also Examples
View source: R/build.interact.r
Calculation of the interaction strength between all amino acids of a protein.
1 | build.interact(cseq, mj1, mj2 = mj1, d, alpha = 82)
|
cseq |
(coded) amino acid sequence |
mj1 |
matrix for intrachain interactions between amino acids |
mj2 |
matrix for interchain interactions between amino acids |
d |
vector of chain lengths |
alpha |
strength of the peptide bond |
Per default only one matrix with interaction strengths between amino acids is used to compute the interaction map. But it is possible to differentiate between intrachain [Miyazawa and Jernigan (1996)] and interchain [Keskin et al. (1998)] interactions by using an additional interaction matrix. Both matrices are included in the package and can be used directly. The value of d
specifies the lengths of the different chains in the sequence.
Per definition there is no interaction between an amino acid and itself, therefore the diagonal is set to zero. All entries that represent peptide bonds between two amino acides are set to alpha
(only entries in the off-diagonals).
The input of cseq
can be a numeric sequences [0,19] as well as a string vector of standard amino acid sequences in one or three-letter code.
Return value is a symmetric matrix specifying the interactions between the amino acids. If n is the length of the sequence, the matrix dimension results in n \times n.
Choosing one's own matrices for inter- and intrachain interactions requires some thought: this funcion only deals with the standard amino acids thus it requires always symmetric, 20 \times 20 matrices.
Franziska Hoffgaard
Miyazawa and Jernigan (1996) Journal of Molecular Biology 256, 623–644.
Hamacher (2006) Journal of Chemical Theory and Computation 2, 873–878.
Keskin, Bahar, Badretdinov, Ptitsyn and Jernigan (1998) Protein Science 7, 2578–2586.
1 2 3 4 5 6 7 8 9 | #Miyazawa/Jernigan matrix
mat<-as.matrix(read.table(system.file("mj1.txt", package = "BioPhysConnectoR")))
#Keskin matrix
mat2<-as.matrix(read.table(system.file("mj2.txt", package = "BioPhysConnectoR")))
cseq<-rep(1:5, 2)
d<-c(5, 5)
im<-build.interact(cseq = cseq, mj1 = mat, mj2 = mat2, d = d)
|
Loading required package: snow
Loading required package: matrixcalc
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.