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.
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.
Molecular evolution is the process of change in the sequence composition of cellular molecules such as DNA, RNA, and proteins across generations. The field of molecular evolution uses principles of evolutionary biology and population genetics to explain patterns in these changes. Major topics in molecular evolution concern the rates and impacts of single nucleotide changes, neutral evolution vs. natural selection, origins of new genes, the genetic nature of complex traits, the genetic basis of speciation, the evolution of development, and ways that evolutionary forces influence genomic and phenotypic changes.
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.
Proteinogenic amino acids are amino acids that are incorporated biosynthetically into proteins during translation. The word "proteinogenic" means "protein creating". Throughout known life, there are 22 genetically encoded (proteinogenic) amino acids, 20 in the standard genetic code and an additional 2 that can be incorporated by special translation mechanisms.
In genetics and bioinformatics, a single-nucleotide polymorphism is a germline substitution of a single nucleotide at a specific position in the genome that is present in a sufficiently large fraction of considered population.
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.
In evolutionary biology, conserved sequences are identical or similar sequences in nucleic acids or proteins across species, or within a genome, or between donor and receptor taxa. Conservation indicates that a sequence has been maintained by natural selection.
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.
Directed evolution (DE) is a method used in protein engineering that mimics the process of natural selection to steer proteins or nucleic acids toward a user-defined goal. It consists of subjecting a gene to iterative rounds of mutagenesis, selection and amplification. It can be performed in vivo, or in vitro. Directed evolution is used both for protein engineering as an alternative to rationally designing modified proteins, as well as for experimental evolution studies of fundamental evolutionary principles in a controlled, laboratory environment.
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.
In molecular biology, the cold-shock domain (CSD) is a protein domain of about 70 amino acids which has been found in prokaryotic and eukaryotic DNA-binding proteins. Part of this domain is highly similar to the RNP-1 RNA-binding motif.
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.