Storage space for a contact map
ContactMap
(id)[source]¶Bases: conkit.core.Entity.Entity
A contact map object representing a single prediction
The ContactMap
class represents a data structure to hold a single
contact map prediction in one place. It contains functions to store,
manipulate and organise Contact
instances.
Examples
>>> from conkit.core import Contact, ContactMap
>>> contact_map = ContactMap("example")
>>> contact_map.add(Contact(1, 10, 0.333))
>>> contact_map.add(Contact(5, 30, 0.667))
>>> print(contact_map)
ContactMap(id="example" ncontacts=2)
Attributes
coverage |
The sequence coverage score |
id |
The ID of the selected entity |
ncontacts |
The number of Contact instances in the ContactMap |
precision |
The precision (Positive Predictive Value) score |
repr_sequence |
The representative Sequence associated with the ContactMap |
repr_sequence_altloc |
The representative altloc Sequence associated with the ContactMap |
sequence |
The Sequence associated with the ContactMap |
top_contact |
The first Contact entry in ContactMap |
Methods
add (entity) |
Add a child to the Entity |
assign_sequence_register ([altloc]) |
Assign the amino acids from Sequence to all Contact instances |
calculate_jaccard_index (other) |
Calculate the Jaccard index between two ContactMap instances |
calculate_scalar_score () |
Calculate a scaled score for the ContactMap |
copy () |
Create a shallow copy of Entity |
deepcopy () |
Create a deep copy of Entity |
find (indexes[, altloc]) |
Find all contacts associated with index |
match (other[, remove_unmatched, renumber, ...]) |
Modify both hierarchies so residue numbers match one another. |
plot_map (*args, **kwargs) |
Produce a 2D contact map plot |
remove (id) |
Remove a child |
remove_neighbors ([min_distance, inplace]) |
Remove contacts between neighboring residues |
rescale ([inplace]) |
Rescale the raw scores in ContactMap |
sort (kword[, reverse, inplace]) |
Sort the ContactMap |
assign_sequence_register
(altloc=False)[source]¶Assign the amino acids from Sequence
to all Contact
instances
Parameters: | altloc : bool
|
---|
calculate_jaccard_index
(other)[source]¶Calculate the Jaccard index between two ContactMap
instances
This score analyzes the difference of the predicted contacts from two maps,
where \(x\) and \(y\) are the sets of predicted contacts from two different predictors, \(\left|x \cap y\right|\) is the number of elements in the intersection of \(x\) and \(y\), and the \(\left|x \cup y\right|\) represents the number of elements in the union of \(x\) and \(y\).
The J-score has values in the range of \([0, 1]\), with a value of \(1\) corresponding to identical contact maps and \(0\) to dissimilar ones.
Parameters: | other :
|
---|---|
Returns: | float
|
Warning
The Jaccard distance ranges from \([0, 1]\), where \(1\) means the maps contain identical contacts pairs.
Notes
The Jaccard index is different from the Jaccard distance mentioned in [1]. The Jaccard distance corresponds to \(1-Jaccard_{index}\).
[1] | Q. Wuyun, W. Zheng, Z. Peng, J. Yang (2016). A large-scale comparative assessment of methods for residue-residue contact prediction. Briefings in Bioinformatics, [doi: 10.1093/bib/bbw106]. |
calculate_scalar_score
()[source]¶Calculate a scaled score for the ContactMap
This score is a scaled score for all raw scores in a contact map. It is defined by the formula
where \(x\) corresponds to the raw score of each predicted contact and \(\overline{d}\) to the mean of all raw scores.
The score is saved in a separate Contact
attribute called
scalar_score
This score is described in more detail in [2].
[2] | S. Ovchinnikov, L. Kinch, H. Park, Y. Liao, J. Pei, D.E. Kim, H. Kamisetty, N.V. Grishin, D. Baker (2015). Large-scale determination of previously unsolved protein structures using evolutionary information. Elife 4, e09248. |
coverage
¶The sequence coverage score
The coverage score is calculated by analysing the number of residues covered by the predicted contact pairs.
The coverage score is calculated by dividing the number of contacts \(x_{cov}\) by the number of residues in the sequence \(L\).
Returns: | cov : float
|
---|
See also
find
(indexes, altloc=False)[source]¶Find all contacts associated with index
Parameters: | index : list, tuple
altloc : bool
|
---|---|
Returns: |
|
match
(other, remove_unmatched=False, renumber=False, inplace=False)[source]¶Modify both hierarchies so residue numbers match one another.
This function is key when plotting contact maps or visualising contact maps in 3-dimensional space. In particular, when residue numbers in the structure do not start at count 0 or when peptide chain breaks are present.
Parameters: | other :
remove_unmatched : bool, optional
renumber : bool, optional
inplace : bool, optional
|
---|---|
Returns: | hierarchy_mod
|
Raises: | ValueError
|
ncontacts
¶The number of Contact
instances in the ContactMap
Returns: | ncontacts : int
|
---|
plot_map
(*args, **kwargs)[source]¶Produce a 2D contact map plot
Warning
This function has been deprecated. Please use conkit.plot.ContactMapFigure
instead.
precision
¶The precision (Positive Predictive Value) score
The precision value is calculated by analysing the true and false postive contacts.
The status of each contact, i.e true or false positive status, can be
determined by running the match()
function providing a reference
structure.
Returns: | ppv : float
|
---|
See also
remove_neighbors
(min_distance=5, inplace=False)[source]¶Remove contacts between neighboring residues
Parameters: | min_distance : int, optional
inplace : bool, optional
|
---|---|
Returns: | contact_map :
|
repr_sequence
¶The representative Sequence
associated with the ContactMap
The peptide sequence constructed from the available contacts using the normal res_seq positions
Returns: | sequence : |
---|---|
Raises: | TypeError
|
See also
repr_sequence_altloc
¶The representative altloc Sequence
associated with the ContactMap
The peptide sequence constructed from the available contacts using the altloc res_seq positions
Returns: | sequence : |
---|---|
Raises: | ValueError
|
See also
rescale
(inplace=False)[source]¶Rescale the raw scores in ContactMap
Rescaling of the data is done to normalize the raw scores to be in the range [0, 1]. The formula to rescale the data is:
\(x\) is the original value and \(d\) are all values to be rescaled.
Parameters: | inplace : bool, optional
|
---|---|
Returns: | contact_map :
|
sequence
¶The Sequence
associated with the ContactMap
Returns: | Sequence |
---|
See also
sort
(kword, reverse=False, inplace=False)[source]¶Sort the ContactMap
Parameters: | kword : str
reverse : bool, optional
inplace : bool, optional
|
---|---|
Returns: | contact_map :
|
Raises: | ValueError
|
top_contact
¶The first Contact
entry in ContactMap
Returns: | top_contact :
|
---|