# Computing and visually analyzing mutual information in molecular co-evolution

- Sebastian Bremm
^{1}Email author, - Tobias Schreck
^{1}, - Patrick Boba
^{2}, - Stephanie Held
^{2}and - Kay Hamacher
^{2}

**11**:330

https://doi.org/10.1186/1471-2105-11-330

© Bremm et al; licensee BioMed Central Ltd. 2010

**Received: **2 February 2010

**Accepted: **17 June 2010

**Published: **17 June 2010

## Abstract

### Background

Selective pressure in molecular evolution leads to uneven distributions of amino acids and nucleotides. In fact one observes correlations among such constituents due to a large number of biophysical mechanisms (folding properties, electrostatics, ...). To quantify these correlations the mutual information -after proper normalization - has proven most effective. The challenge is to navigate the large amount of data, which in a study for a typical protein cannot simply be plotted.

### Results

To visually analyze mutual information we developed a matrix visualization tool that allows different views on the mutual information matrix: filtering, sorting, and weighting are among them. The user can interactively navigate a huge matrix in real-time and search e.g., for patterns and unusual high or low values. A computation of the mutual information matrix for a sequence alignment in FASTA-format is possible. The respective stand-alone program computes in addition proper normalizations for a null model of neutral evolution and maps the mutual information to *Z*-scores with respect to the null model.

### Conclusions

The new tool allows to compute and visually analyze sequence data for possible co-evolutionary signals. The tool has already been successfully employed in evolutionary studies on HIV1 protease and acetylcholinesterase. The functionality of the tool was defined by users using the tool in real-world research. The software can also be used for visual analysis of other matrix-like data, such as information obtained by DNA microarray experiments. The package is platform-independently implemented in Java and free for academic use under a GPL license.

## Background

### Background of the analysis problem

The understanding of molecular evolution requires a detailed understanding of the dynamics *among* the constituents of a molecule during its evolution [1]. Computational biologists seek evolutionary signals in data sets of sequences of biomolecules and signatures of a correlated dynamic of the evolutionary processes shaping the characteristics of molecules under investigation. Such co-evolution occurs mostly, when either an amino acid or a nucleotide within a biomolecule evolves in concert with another "site" within the same or a partnering molecule. The mutual information (MI) is widely used to detect such correlated evolutionary dynamics. The computation of the MI (see section Information-theoretical measure for details) itself is straightforward. However, proper normalization needs to be carefully taken into account [2–6] for typical, finite-sized data sets. This issue is discussed in more detail below. The problem of analyzing the obtained mutual information values was, however, not tackled until now: as the MI is a quantification between any two sites within a protein, of e.g., N amino acids, for such a molecule we need to compute and to analyze
MI values. Even for modest sized proteins with *N* ~ 100 this means to visualize 5050 real values. Typically one cannot easily grasp structures, scales, etc. in such large amounts of data.

### Visual-Interactive Approach and Matrix Visualization

Recently, visualization has been widely recognized as a promising approach to help analysts and researchers to better understand such large amounts of complex data. The approach suggests to have visual-interactive displays appropriately encode information using visual mappings; and let the user interactively manipulate these displays to navigate, drill-down, and explore [7–9].

According to the structure of the data to be visualized, different visual mappings are appropriate. Matrix visualization [10] is appropriate for large amounts of data elements for which pairwise relationships with quantitative attributes are given. By representing the quantitative value of each relationship by color, matrix visualization is highly scalable, ultimately representing each relationship by a single pixel.

Important problems to address in designing effective matrix visualization systems involve choosing an appropriate color scale [11], data preprocessing steps, and applying suitable sorting on the matrices to be visualized. The latter is specifically important, as is allows to make assessments on the overall structure of the relationships. Matrix sorting usually arranges rows and columns of the matrix by similarity, with an appropriate similarity function defined on the row and column vectors of the matrix.

Matrix visualization as a technique is well-known and to date has found its way into software systems such as R [12] or Matlab [13]. However, many implementations are focused on producing static images, offering only limited support for interactive parameter change and navigation in the matrix display by the user. We therefore developed a fully interactive matrix visualization system in Java. It allows the user to change important parameters and navigate the data on the fly by means of a two-stage zooming mechanism. Furthermore, we support the joint visualization of two matrices, supporting our specific analysis problem.

## Implementation

### Visualization

The implemented application generically supports two data matrices: one data and one weight matrix (see section Normalization & Weighting) which can be inspected individually or jointly. The general approach is to map the normalized matrix values to an appropriate color map and display it as a grid. While using color is typically not the first choice for representing absolute values, it allows comparative analysis of value ranges and provides a highly compact view of the overall data distribution. Specifically in the case of large data matrices and in conjunction with an appropriate matrix sorting mechanism it allows the assessment of the features of the matrix. The application supports detailed visualization of either one of the two matrices. Detailed information is available by interactively *zooming* into parts of the matrix and *restricting* the displayed data to specific value ranges. Thereby, the approach follows Shneiderman's Information Visualization Mantra ("Overview first - zoom, filter, refine - details on demand" [14]). Sorting the matrix by arranging rows and columns by similarity reduces its complexity and allows identification of systematic (similar) relationships between entities of the experiment by homogeneous colors [10].

Our sorting algorithm works by finding a so-called *seed* row according to the maximum of the sum of contained values. This row is made the top row. Then, the sorting algorithm among the remaining rows finds the one that is most similar to the seed row, where the degree of (dis)similarity is measured by the *l*_{1} norm between the respective row vectors. The algorithm places the found row just below the seed row, makes the found row the new seed row, and iteratively continues until all rows have been processed. The same approach is then applied on the columns of the matrix. This sorting algorithm is rather simple, yet provides a useful starting point for the visual analysis. The algorithm has quadratic runtime complexity. For an overview of the application design, please see the system illustration and description provided in section Results and Discussion.

### Mutual information computation

The stand-alone program micato (m utual i nformation ca lculation to ol) reads a sequence file in FASTA format and calculates the MI of the sequence contained therein. By separation from the visualization tool micato can be run on e.g., clusters using job-scheduling systems. This is useful in particular for sampling large instances of null models for normalization.

To this end micato calculates in a first step the sequence entropies of each column of the sequence alignment and stores it. Then the joint entropy of each pairing of two columns is calculated and by equation 1 the MI is calculated and stored in a matrix MI_{
ij
}for a pair of positions (*i*, *j*). Then micato runs a user defined number of independent column shuffles to generate a statistically significant number of instances of the null model (see section Normalization & Weighting for details). The MI matrix is exported as a CSV file, as well as the *Z*-scores of those MI values with respect to the statistics of the null model. The CSV format can be read by the MIMatrixViz program without further conversion.

## Methods

### Information-theoretical measure

where *x* and *y* are outcomes for random variables *X*_{
i
}and *Y*_{
j
}drawn from a symbol set
, taken from a multiple sequence alignment as columns *i* and *j*. The symbols *H*_{
i
}(*x*), *H*_{
j
}(*y*) are the column sequence entropies and *H*_{
ij
}(*x*, *y*) is the entropy computed by the joint probability function *p*_{
ij
}(*x*, *y*). Repeated application of the equation leads to a symmetric MI matrix (MI_{
ij
}) for all pairings (*i*, *j*).

In studies on proteins the symbol set consists of the 20 standard amino acids
, which has to be expanded to include the gap character and an extra character for non-standard amino acids
. We set the probabilities *p*_{
i
}(*x*),*p*_{
j
}(*y*), and *p*_{
ij
}(*x*, *y*) to the observed frequencies of amino acids within the columns of a multiple sequence alignment. This can be done with the supplied routines in the program micato.

### Normalization & Weighting

Although in the post-genomic era [16] we have access to huge databases of sequences, in a typical setting the number of sequences available is still only of the order 10^{2}. We have previously shown that this limitation might lead to substantial finite-size effects in the computed MI values [2]. These effects can be compensated by normalizing the obtained MI values to a null model of evolution [17].

*shuffled*columns: while maintaining the one-column sequence entropies with this procedure we destroy any correlation between any two columns

*i*and

*j*and therefore obtain distributions for MI

_{ ij }with respect to naturally occurring, independent symbol frequencies. From this distribution we can easily compute

*Z*values for any MI

_{ ij }value by

where
is the average of the MI values for the shuffled columns *i* and *j*, and *σ*(*MI*_{
ij
}) the standard deviation of this sample. Note that this has to be done for each column pair (*i*, *j*) independently. This protocol is also implemented in the program micato and it can be loaded into our visualization program as a weighting scheme.

## Results and Discussion

The implementation of all components like user interface, data storage and algorithms is kept modular, so it is easy to adopt the system to upcoming needs. This includes e.g., additional matrix ordering methods, color schemes or data filtering mechanisms. Currently, filtering can be done in two steps. Firstly, by setting a maximum and minimum value of the mutual information. Secondly, by doing the same for the *Z*-score. By this the *Z*-score can be used as filter or additionally as a weight for controlling the color saturation of the corresponding matrix entries.

The visualization system can be used additionally with any matrix valued data set beyond co-evolutionary studies. This is achieved by specification of a simple file format for matrix data. The software assumes the raw matrix and weight data to be contained in plain ASCII files, where each row in the file contains one matrix cell entry of the form <*row column value* >.

### Example Protocol

Codoñer et al. [1], as well as other authors [17] have argued that intramolecular co-evolution typically results from a superposition of various biological and biochemical influences, which depend highly on the system under investigation. In fact, the decomposition is the key analysis task in co-evultionary studies. Explorative analysis can then be used for hypothesis creation [18] on the origin of such influences. We applied the software package to a sequence set of variants of the acetylcholinesterase (AChE) (Held S, Hoffgaard F, Hamacher K: Biophysical Annotation of Molecular Coevolution of Acetylcholinesterase, submitted).

*Z*scores and high MI values as illustrated in Figure 2(c). We show how the software can be used to restrict the display to the supposedly important ranges of MI and

*Z*scores. In part (d) of the same figure we show the residues marked in the molecular structure. These are (for the sequence of

*T. californica*) Y70, V236, N280, F284, F288, G335, S345, V360, Q374. Most of these residues belong to the peripheral anionic site (PAS) site of AChE, which is important for establishing contact with the substrate [19]. Also, various AChE inhibitors bind to this site, suggesting a partial explanation why evolutionary signals occur at these spots. The involved residues are subject to a co-evolutionary pressure, the origin of which one can now start to investigate based on the insight we gained from the application of MIMatrixViz to the AChE sequence data.

## Future Trends & Intended Use

The implemented routines can be used in studies on molecular (co-)evolution by working on provided FASTA-files. The MI computation can be run in batch-mode to allow for compute cluster usage. The output of computed mutual information values and their weighting by *Z*-scores can be opened by the interactive matrix visualization software provided. There, as described above the user can navigate the vast amount of data, applying filters in sequence space and in value space, and using *Z*-scores (or other externally provided weights) to estimate and visualize the statistical significance of the mutual information values.

Typically, the visual analytics approach does not guarantee to reveal all relevant features of empirical data. However, in a generic biological application one does not know beforehand what signals to look for. This renders automatic processing ineffective and one has to resort to visual and interactive inspection. Future work includes extending the functionality of the visualization software with additional functionality. First, additional matrix sorting algorithms with user-settable sorting criteria should be included, allowing the user to take suitable views on the data set. The matrix display should be extended by side views showing the similarity of rows and columns as well as the reordering (confusion) index in case the matrix has been sorted. In the long run, the system should be integrated with additional relevant meta data, and linked with additional viewing components such as 3D molecular viewers. Our software is provided as Java Bytecode. The sourcecode can be made available upon request. We are also open for collaboration aiming at improving the functionality of the software and applying it to new use cases.

## Conclusions

In the MIMatrixViz package we provide routines to compute mutual information of evolutionary dynamics in molecules. The package is capable of normalizing those values and therefore accounts for finite-sized data sets. The visualization part is separated from this to allow batch-usage on servers and clusters for sufficient statistics. The visual approach allows to interactively explore the data, and investigate patterns, structures, and particular interesting spots within the mutual information matrices.

The user can generate graphics and filtered data sets with the package in publication ready quality. To this end, the application allows the user to export matrix images in the lossless PNG file format, and to export selected data subsets as plain ASCII files.

Other matrix-oriented data, as e.g., obtained by DNA microarray experiments, can be visually analyzed with the tool, too. External knowledge can be incorporated by the weight matrix to augment the insight one gains from the expression levels detected at the feature sites. Potential scenarios include phylogenetic likelihoods for particular hits on reporters, gauging bias to cope with potential shortcomings in the production and/or binding processes.

## Availability and Requirements

**Project name:**
MIMatrixViz

**Project home page:**
http://www.gris.informatik.tu-darmstadt.de/projects/vsa/matrixvis/

**Operating system:** Platform independent (Requires a Java Virtual Machine (JVM) on the target system)

**Programming language:**
Java

**Requirements:** Java Runtime Environment ≥ 1.6

**License**: GPL for academic users.

**Any restrictions to use by non-academics:** For commercial applications of MIMatrixViz, please contact the authors.

## Declarations

## Authors’ Affiliations

## References

- Codoñer FM, Fares MA: Why Should We Care About Molecular Coevolution?
*Evolutionary Bioinformatics*2008, 4(2):29–38.Google Scholar - Weil P, Hoffgaard F, Hamacher K: Estimating Sufficient Statistics in Co-Evolutionary Analysis by Mutual Information.
*Computational Biology and Chemistry*2009, 33(6):440–444. 10.1016/j.compbiolchem.2009.10.003View ArticlePubMedGoogle Scholar - Gouveia-Oliveira R, Pedersen A: Finding coevolving amino acid residues using row and column weighting of mutual information and multi-dimensional amino acid representation.
*Algorithms for Molecular Biology*2007, 2: 12. [http://www.almob.org/content/2/1/12] 10.1186/1748-7188-2-12View ArticlePubMedPubMed CentralGoogle Scholar - Hoffgaard F, Weil P, Hamacher K: BioPhysConnectoR: Connecting Sequence Information and Biophysical Models.
*BMC Bioinformatics*2010, 11: 199. [http://www.biomedcentral.com/1471–2105/11/199] 10.1186/1471-2105-11-199View ArticlePubMedPubMed CentralGoogle Scholar - Buslje CM, Santos J, Delfino JM, Nielsen M: Correction for phylogeny, small number of observations and data redundancy improves the identification of coevolving amino acid pairs using mutual information.
*Bioinformatics*2009, 25(9):1125–1131. 10.1093/bioinformatics/btp135View ArticlePubMedPubMed CentralGoogle Scholar - Gloor G, Martin L, Wahl L, Dunn S: Mutual information in protein multiple sequence alignments reveals two classes of coevolving positions.
*Biochemistry*2005, 44(19):7156–7165. 10.1021/bi050293eView ArticlePubMedGoogle Scholar - Hansen CD, Johnson CR (Eds):
*The visualization handbook*Amsterdam [u.a.]: Elsevier; 2005.Google Scholar - Card SC, Mackinlay J, Shneiderman B:
*Readings in Information Visualization: Using Vision to Think.*Morgan Kaufmann Publishers; 1999.Google Scholar - Thomas J, Cook K: Illuminating the Path: The Research and Development Agenda for Visual Analytics. 2005. IEEE Computer SocietyGoogle Scholar
- Wu HM, Tzeng S, Chen CH: Handbook of Data Visualization. 2008, 681–708. Springer. chap. Matrix VisualizationView ArticleGoogle Scholar
- Levkowitz H, Herman GT: Color Scales for Image Data.
*IEEE Comput Graph Appl*1992, 12: 72–80. 10.1109/38.135886View ArticleGoogle Scholar - The R Project for Statistical Computing[http://www.r-project.org/]
- MathWorks T: MATLAB.[http://www.mathworks.com/]
- Shneiderman B: The Eyes Have It: A Task by Data Type Taxonomy for Information Visualizations. In
*VL '96: Proceedings of the 1996 IEEE Symposium on Visual Languages*. Washington, DC, USA: IEEE Computer Society; 1996:336. full_textView ArticleGoogle Scholar - MacKay D:
*Information Theory, Inference, and Learning Algorithms.*2nd edition. Cambridge: Cambridge University Press; 2004.Google Scholar - Lengauer T: Bioinformatics: From the Pre-genomic to the Post-genomic Era.
*ERCIM News*2000, 43: 6–7.Google Scholar - Hamacher K: Relating Sequence Evolution of HIV1-Protease to Its Underlying Molecular Mechanics.
*Gene*2008, 422: 30–36. 10.1016/j.gene.2008.06.007View ArticlePubMedGoogle Scholar - Keim DA: Information Visualization and Visual Data Mining.
*IEEE Transactions on Visualization and Computer Graphics (TVCG)*2002, 8: 1–8. 10.1109/2945.981847View ArticleGoogle Scholar - Bourne Y, Taylor P, Radié Z, Marchot P: Structural insights into ligand interactions at the acetylcholinesterase peripheral anionic site.
*EMBO Journal*2003, 22: 1–12. 10.1093/emboj/cdg005View ArticlePubMedPubMed CentralGoogle Scholar - Humphrey W, Dalke A, Schulten K: VMD - Visual Molecular Dynamics.
*Journal of Molecular Graphics*1996, 14: 33–38. 10.1016/0263-7855(96)00018-5View ArticlePubMedGoogle Scholar

## Copyright

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.