Low complexity regions (LCRs) in protein sequences, also defined in some contexts as compositionally biased regions (CBRs), are regions in protein sequences that differ from the composition and complexity of most proteins that is normally associated with globular structure. [1] [2] LCRs have different properties from normal regions regarding structure, function and evolution.
LCRs were originally thought to be unstructured and flexible linkers that served to separate the structured (and functional) domains of complex proteins, [3] but they are also capable of forming secondary structures, like helices (more often) and even sheets. [4] They may play a structural role in proteins such as collagens, myosin, keratins, silk, cell wall proteins. [5] Tandem repeats of short oligopeptides that are rich in glycine, proline, serine or threonine are capable of forming flexible structures that bind ligands under certain pH and temperature conditions. [6] Proline is a well-known alpha-helix breaker, however, amino acid repeats composed of proline may form poly-proline helices. [7]
LCRs were originally thought as ‘junk’ regions or as neutral linkers between domains; however, experimental and computational evidence increasingly indicates that they may play important adaptive and conserved roles, relevant to biotechnology, heterologous protein expression, medicine, as well as to our understanding of protein evolution. [8]
LCRs of eukaryotic proteins have been involved in human diseases, [9] [10] especially neurodegenerative ones, where they tend to form amyloids in humans and other eukaryotes. [11]
They have been reported to have adhesive roles, [12] function in excreted sticky proteins used for prey capture, [13] or have roles as transducers of molecular movement, e.g. in the prokaryotic TonB/TolA systems. [14]
LCRs may form surfaces for interaction with phospholipid bilayers, [15] or as positive charge clusters for DNA binding, [8] [16] [17] or as negative or even histidine-acidic charge clusters for coordinating calcium, magnesium or zinc ions. [8] [16]
They may also play important roles in protein translation, as tRNA ‘sponges’, slowing down translation in order to allow time for the correct folding of the nascent polypeptide chain. [18] They may even function as frame-shift checkpoints, by shifting to an unusual amino acid content that makes the protein highly unstable or insoluble, which in turn triggers fast recycling, before any further cellular damage. [19] [20]
Analyses on model and non-model eukaryotic proteomes have revealed that LCRs are frequently found in proteins involved in binding of nucleic acids (DNA or RNA), in transcription, receptor activity, development, reproduction and immunity whereas metabolic proteins are depleted of LCRs. [3] [21] [22] [23] A bioinformatics study of the Uniprot annotation of LCR containing proteins observed that 44% (9751/22259) of Bacterial and 44% (662/1521) of Archaeal LCRs are detected in proteins of unknown function, however, a significant number of proteins of known function (from many different species), especially those involved in translation and the ribosome, nucleic acid binding, metal-ion binding, and protein folding were also found to contain LCRs. [8]
LCRs are more abundant in eukaryotes, but they also have a significant presence in many prokaryotes. [8] On average, 0.05 and 0.07% of the bacterial and archaeal proteomes (total amino acids of LCRs in a given proteome/total amino acids of that proteome) form LCRs whereas for five model eukaryotic proteomes (human, fruitfly, yeast, fission yeast, Arabidopsis) this coverage was significantly higher (on average, 0.4%; between 2 and 23 times higher than prokaryotes). [8]
Eukaryotic LCRs tend to be longer than prokaryotic LCRs. [8] The average size of a eukaryotic LCR is 42 amino acids long, whereas bacterial, archaeal and phage LCRs are 38, 36 and 33 amino acids long, respectively. [8]
In the Archaea, the halobacterium Natrialba magadii has the highest number of LCRs and the highest enrichment for LCRs. [8] In Bacteria, Enhygromyxa salina, a delta proteobacterium that belongs to myxobacteria has the highest number of LCRs and the highest enrichment for LCRs. [8] Intriguingly, four of the top five bacteria with the highest enrichment for LCRs are also myxobacteria. [8]
The three most enriched amino acids within LCRs of Bacteria are proline, glycine and alanine, whereas in Archaea they are threonine, aspartate and proline. [8] In Phages, they are alanine, glycine and proline. [8] Glycine and proline emerge as very enriched amino acids in all three evolutionary lineages, whereas alanine is highly enriched in Bacteria and Phages but not enriched in Archaea. On the other hand, hydrophobic (M, I, L, V) and aromatic amino acids (F, Y, W) as well as cysteine, arginine and asparagine are heavily under-represented in LCRs. [8] Very similar trends for amino acids with a high (G, A, P, S, Q) and low (M, V, L, I, W, F, R, C) occurrence within LCRs have been observed in eukaryotes as well. [24] [21] This observed pattern of certain amino acids being over-represented (enriched for) or under-represented in LCRs could be partially explained by the energy cost for synthesis or metabolism of each of the amino acids. [8] Another possible explanation, which does not exclude the previous explanation of energy cost could be the reactivity of certain amino acids. [8] For example, Cysteine is a very reactive amino acid that would not be tolerated in high numbers within a small region of a protein. [25] Similarly, extremely hydrophobic regions can form non-specific protein–protein interactions among themselves and with other moderately hydrophobic regions [26] [27] in mammalian cells. Thus, their presence may disturb the balance of protein-protein interaction networks within the cell, especially if the carrier proteins are highly expressed. [8] A third explanation may be based on micro-evolutionary forces and, more specifically, on the bias of DNA polymerase slippage for certain di- tri- or tetra-nucleotides . [8]
A bioinformatics analysis of prokaryotic LCRs identified 5 types of amino acid enrichment, for certain functional categories of LCRs: [8]
Based on the above observations and analyses, a Neural Network webserver named LCR-hound has been developed to predict LCRs and their function. [8]
LCRs are very interesting from a micro and macro evolutionary perspective. [8] They may be generated by DNA slippage, recombination and repair. [28] Thus, they are linked to recombination hotspots and may even possibly facilitate cross-over. [29] [30] By originating from genetic instability, they may cause, at the DNA level, a certain region of the protein to expand or contract and even cause frame-shifts (phase-variants) that affect microbial pathogenicity or provide raw material for evolution. [31] Most intriguingly, they may provide a window into the very early evolution of life. [8] [32] During early evolution, when only few amino acids were available and the primary genetic code was still expanding its repertoire, the first proteins were assumed to be short, repetitive and therefore, of low complexity. [33] [34] Thus, modern LCRs could represent primordial aspects of the evolution towards the protein world and may provide clues about the functions of the early proto-peptides. [8]
Most studies have focused on the evolution, functional and structural role of eukaryotic LCRs. [8] However, a comprehensive study of prokaryotic LCRs from many diverse prokaryotic lineages provides a unique opportunity to understanding the origin, evolution and nature of these regions. Due to the high effective population size and short generation times of prokaryotes, the de novo emergence of a mildly or moderately deleterious amino acid repeat or LCR should quickly be filtered out by strong selective forces. [8] This must be especially the case for LCRs found in highly expressed proteins, since they should also have a great impact on the energy burden of protein translation. [35] [36] Thus, any prokaryotic LCRs that constitute evolutionary accidents with no functional significance should not be fixed by genetic drift and consequently should not demonstrate any levels of conservation among moderately distant evolutionary relatives. [8] On the contrary, any LCR found among homologs of several moderately distant prokaryotic species should very probably reserve a functional role. [8]
The amino acids with the highest frequency in LCRs are glycine and alanine, with their respective codons GGC and GCC being the most frequent, as well as complementary. [8] In eukaryotes and more specifically in chordates (such as human, mouse, chicken, zebrafish and sea squirt), alanine- and glycine-rich LCRs are over-represented in recently formed LCRs and probably are better tolerated by the cell. [37] Intriguingly, it has also been suggested that they represent the very first two amino acids [38] and codons [34] [39] [40] of the early genetic code. Thus, these two codons and their respective amino acids must have been constituents of the earliest oligopeptides, with a length of 10–55 amino acids [41] and very low complexity. Based on several different criteria and sources of data, Higgs and Pudritz [38] suggest G, A, D, E, V, S, P, I, L, T as the early amino acids of the genetic code. Trifonov's work largely agrees with this categorization and proposes that the early amino acids in chronological order are G, A, D, V, S, P, E, L, T, R. An evolutionary analysis observed that many of the amino acids of the suggested very early genetic code (with the exception of the hydrophobic ones) are significantly enriched in bacterial LCRs. [8] Most of the later additions to the genetic code are significantly under-represented in bacterial LCRs. [8] They thus hypothesize and propose that, in a cell-free environment, the early genetic code may have also produced low complexity oligo-peptides from valine and leucine. [8] However, later on, within a more complex cellular environment, these highly hydrophobic LCRs became inappropriate or even toxic from a protein interaction perspective and have been selected against ever since. [8] In addition, they further hypothesize that the very early protopeptides did not have a nucleic acid binding role, [8] because DNA and RNA-binding LCRs are highly enriched in glycine, arginine and lysine, however, arginine and lysine are not among the amino acids of the proposed early genetic code.
Low complexity regions in proteins can be computationally detected from sequence using various methods and definitions, as reviewed in. [2] Among the most popular methodologies to identify LCRs is by measuring their Shannon entropy. [1] The lower the value of the calculated entropy, the more homogeneous the region is in terms of amino acid content. In addition, a Neural Network webserver, LCR-hound has been developed to predict the function of an LCR, based on its amino acid or di-amino acid content. [8] Compression-based tools have also been used to perform such analysis providing higher sensitivity while mitigating the risk of overestimation inherent in other methods. [42]
Amino acids are organic compounds that contain both amino and carboxylic acid functional groups. Although over 500 amino acids exist in nature, by far the most important are the 22 α-amino acids incorporated into proteins. Only these 22 appear in the genetic code of life.
The genetic code is the set of rules used by living cells to translate information encoded within genetic material into proteins. Translation is accomplished by the ribosome, which links proteinogenic amino acids in an order specified by messenger RNA (mRNA), using transfer RNA (tRNA) molecules to carry amino acids and to read the mRNA three nucleotides at a time. The genetic code is highly similar among all organisms and can be expressed in a simple table with 64 entries.
In molecular biology, messenger ribonucleic acid (mRNA) is a single-stranded molecule of RNA that corresponds to the genetic sequence of a gene, and is read by a ribosome in the process of synthesizing a protein.
The SRC Homology 3 Domain is a small protein domain of about 60 amino acid residues. Initially, SH3 was described as a conserved sequence in the viral adaptor protein v-Crk. This domain is also present in the molecules of phospholipase and several cytoplasmic tyrosine kinases such as Abl and Src. It has also been identified in several other protein families such as: PI3 Kinase, Ras GTPase-activating protein, CDC24 and cdc25. SH3 domains are found in proteins of signaling pathways regulating the cytoskeleton, the Ras protein, and the Src kinase and many others. The SH3 proteins interact with adaptor proteins and tyrosine kinases. Interacting with tyrosine kinases, SH3 proteins usually bind far away from the active site. Approximately 300 SH3 domains are found in proteins encoded in the human genome. In addition to that, the SH3 domain was responsible for controlling protein-protein interactions in the signal transduction pathways and regulating the interactions of proteins involved in the cytoplasmic signaling.
Non-coding DNA (ncDNA) sequences are components of an organism's DNA that do not encode protein sequences. Some non-coding DNA is transcribed into functional non-coding RNA molecules. Other functional regions of the non-coding DNA fraction include regulatory sequences that control gene expression; scaffold attachment regions; origins of DNA replication; centromeres; and telomeres. Some non-coding regions appear to be mostly nonfunctional, such as introns, pseudogenes, intergenic DNA, and fragments of transposons and viruses. Regions that are completely nonfunctional are called junk DNA.
The coding region of a gene, also known as the coding sequence (CDS), is the portion of a gene's DNA or RNA that codes for a protein. Studying the length, composition, regulation, splicing, structures, and functions of coding regions compared to non-coding regions over different species and time periods can provide a significant amount of important information regarding gene organization and evolution of prokaryotes and eukaryotes. This can further assist in mapping the human genome and developing gene therapy.
In biology, translation is the process in living cells in which proteins are produced using RNA molecules as templates. The generated protein is a sequence of amino acids. This sequence is determined by the sequence of nucleotides in the RNA. The nucleotides are considered three at a time. Each such triple results in addition of one specific amino acid to the protein being generated. The matching from nucleotide triple to amino acid is called the genetic code. The translation is performed by a large complex of functional RNA and proteins called ribosomes. The entire process is called gene expression.
In genetics and bioinformatics, a single-nucleotide polymorphism is a germline substitution of a single nucleotide at a specific position in the genome. Although certain definitions require the substitution to be present in a sufficiently large fraction of the population, many publications do not apply such a frequency threshold.
A point mutation is a genetic mutation where a single nucleotide base is changed, inserted or deleted from a DNA or RNA sequence of an organism's genome. Point mutations have a variety of effects on the downstream protein product—consequences that are moderately predictable based upon the specifics of the mutation. These consequences can range from no effect to deleterious effects, with regard to protein production, composition, and function.
A leucine zipper is a common three-dimensional structural motif in proteins. They were first described by Landschulz and collaborators in 1988 when they found that an enhancer binding protein had a very characteristic 30-amino acid segment and the display of these amino acid sequences on an idealized alpha helix revealed a periodic repetition of leucine residues at every seventh position over a distance covering eight helical turns. The polypeptide segments containing these periodic arrays of leucine residues were proposed to exist in an alpha-helical conformation and the leucine side chains from one alpha helix interdigitate with those from the alpha helix of a second polypeptide, facilitating dimerization.
A DNA-binding domain (DBD) is an independently folded protein domain that contains at least one structural motif that recognizes double- or single-stranded DNA. A DBD can recognize a specific DNA sequence or have a general affinity to DNA. Some DNA-binding domains may also include nucleic acids in their folded structure.
In biology, the word gene has two meanings. The Mendelian gene is a basic unit of heredity. The molecular gene is a sequence of nucleotides in DNA that is transcribed to produce a functional RNA. There are two types of molecular genes: protein-coding genes and non-coding genes.
Gene structure is the organisation of specialised sequence elements within a gene. Genes contain most of the information necessary for living cells to survive and reproduce. In most organisms, genes are made of DNA, where the particular DNA sequence determines the function of the gene. A gene is transcribed (copied) from DNA into RNA, which can either be non-coding (ncRNA) with a direct function, or an intermediate messenger (mRNA) that is then translated into protein. Each of these steps is controlled by specific sequence elements, or regions, within the gene. Every gene, therefore, requires multiple sequence elements to be functional. This includes the sequence that actually encodes the functional protein or ncRNA, as well as multiple regulatory sequence regions. These regions may be as short as a few base pairs, up to many thousands of base pairs long.
In molecular biology and genetics, DNA annotation or genome annotation is the process of describing the structure and function of the components of a genome, by analyzing and interpreting them in order to extract their biological significance and understand the biological processes in which they participate. Among other things, it identifies the locations of genes and all the coding regions in a genome and determines what those genes do.
Genome evolution is the process by which a genome changes in structure (sequence) or size over time. The study of genome evolution involves multiple fields such as structural analysis of the genome, the study of genomic parasites, gene and ancient genome duplications, polyploidy, and comparative genomics. Genome evolution is a constantly changing and evolving field due to the steadily growing number of sequenced genomes, both prokaryotic and eukaryotic, available to the scientific community and the public at large.
Edward Nikolayevich Trifonov is a Russian-born Israeli molecular biophysicist and a founder of Israeli bioinformatics. In his research, he specializes in the recognition of weak signal patterns in biological sequences and is known for his unorthodox scientific methods.
Proline-rich protein 21 (PRR21) is a protein of the family of proline-rich proteins. It is encoded by the PRR21 gene, which is found on human chromosome 2, band 2q37.3. The gene exists in several species, both vertebrates and invertebrates, including humans. However, the protein have few conserved regions among species.
An array of protein tandem repeats is defined as several adjacent copies having the same or similar sequence motifs. These periodic sequences are generated by internal duplications in both coding and non-coding genomic sequences. Repetitive units of protein tandem repeats are considerably diverse, ranging from the repetition of a single amino acid to domains of 100 or more residues.
CIMAP1C is a gene in humans that encodes the CIMAP1C protein. It is also often referred to as ODF3L1. CIMAP1C is expressed in low levels throughout the body with high expression levels in the testes. It is highly conserved in mammals and reptiles but not present in birds or amphibians, indicating it arose around 300 million years ago.