Search

Patterns observed by examining the evolutionary relationships among proteins of common origin can reveal the structural and functional importance of specific residue positions. In particular, amino acids that are highly conserved (i.e., their positions evolve at a slower rate than other positions) are particularly likely to be of biological importance, for example, for ligand binding. ConSurf is a bioinformatics tool for accurately estimating the evolutionary rate of each position in a protein family. Here we introduce a new release of ConSurf‐DB, a database of precalculated ConSurf evolutionary conservation profiles for proteins of known structure. ConSurf‐DB provides high‐accuracy estimates of the evolutionary rates of the amino acids in each protein. A reliable estimate of a query protein's evolutionary rates depends on having a sufficiently large number of effective homologues (i.e., nonredundant yet sufficiently similar). With current sequence data, ConSurf‐DB covers 82% of the PDB proteins. It will be updated on a regular basis to ensure that coverage remains high—and that it might even increase. Much effort was dedicated to improving the user experience. The repository is available at https://consurfdb.tau.ac.il/. Broader audience By comparing a protein to other proteins of similar origin, it is possible to determine the extent to which each amino acid position in the protein evolved slowly or rapidly. A protein's evolutionary profile can provide valuable insights: For example, amino acid positions that are highly conserved (i.e., evolved slowly) are particularly likely to be of structural and/or functional importance, for example, for ligand binding and catalysis. We introduce here a new and improved version of ConSurf‐DB, a continually updated database that provides precalculated evolutionary profiles of proteins with known structure.

In mammals, cellular identity is defined through strict regulation of chromatin modifications and DNA methylation that control gene expression. Methylation of cytosines at CpG sites in the genome is mainly associated with suppression; however, the reason for enhancer-specific methylation is not fully understood. We used sequential ChIP-bisulfite-sequencing for H3K4me1 and H3K27ac histone marks. By collecting data from the same genomic region, we identified enhancers differentially methylated between these two marks. We observed a global gain of CpG methylation primarily in H3K4me1-marked nucleosomes during mouse embryonic stem cell differentiation. This gain occurred largely in enhancer regions that regulate genes critical for differentiation. The higher levels of DNA methylation in H3K4me1- versus H3K27ac-marked enhancers, despite it being the same genomic region, indicates cellular heterogeneity of enhancer states. Analysis of single-cell RNA-seq profiles demonstrated that this heterogeneity correlates with gene expression during differentiation. Furthermore, heterogeneity of enhancer methylation correlates with transcription start site methylation. Our results provide insights into enhancer-based functional variation in complex biological systems., Author summary Cellular dynamics are underlined by numerous regulatory layers. The regulatory mechanism of interest in this work are enhancers. Enhancers are regulatory regions responsible, mainly, for increasing the possibility of transcription of a certain gene. Enhancers are marked by two distinct chemical groups-H3K4me1 and H3K27ac on the tail of histones. Histones are the proteins responsible for DNA packaging into condensed chromatin structure. In contrast, DNA methylation is a chemical modification often found on enhancers, and is traditionally associated with repression. A long-debated question revolves around the functional relevance of DNA methylation in the context of enhancers. Here, we combined the two regulatory layers, histone marks and DNA methylation, to a single measurement that can highlight DNA methylation separately on each histone mark but at the same genomic region. When isolated with H3K4me1, enhancers showed higher levels of methylation compared to H3K27ac. As we measured the same genomic locations, we show that differences of DNA methylation between these marks can only be explained by cellular heterogeneity. We also demonstrated that these enhancers tend to play roles in stem cell differentiation and expression levels of the genes they control correlate with cell-to-cell variation.