
Citation: Cibely Silva Martin, Mateus Dassie Maximino, Matheus Santos Pereira, Clarissa de Almeida Olivati, Priscila Alessio. The role of film composition and nanostructuration on the polyphenol sensor performance[J]. AIMS Materials Science, 2017, 4(1): 27-42. doi: 10.3934/matersci.2017.1.27
[1] | Shusen Yan, Weilin Yu . Planar vortices in a bounded domain with a hole. Electronic Research Archive, 2021, 29(6): 4229-4241. doi: 10.3934/era.2021081 |
[2] | Yiyuan Qian, Haiming Song, Xiaoshen Wang, Kai Zhang . Primal-dual active-set method for solving the unilateral pricing problem of American better-of options on two assets. Electronic Research Archive, 2022, 30(1): 90-115. doi: 10.3934/era.2022005 |
[3] | Abdeljabbar Ghanmi, Hadeel Z. Alzumi, Noureddine Zeddini . A sub-super solution method to continuous weak solutions for a semilinear elliptic boundary value problems on bounded and unbounded domains. Electronic Research Archive, 2024, 32(6): 3742-3757. doi: 10.3934/era.2024170 |
[4] | Yijun Chen, Yaning Xie . A kernel-free boundary integral method for reaction-diffusion equations. Electronic Research Archive, 2025, 33(2): 556-581. doi: 10.3934/era.2025026 |
[5] | Margarida Camarinha . A natural 4th-order generalization of the geodesic problem. Electronic Research Archive, 2024, 32(5): 3396-3412. doi: 10.3934/era.2024157 |
[6] | Matthew Gardner, Adam Larios, Leo G. Rebholz, Duygu Vargun, Camille Zerfas . Continuous data assimilation applied to a velocity-vorticity formulation of the 2D Navier-Stokes equations. Electronic Research Archive, 2021, 29(3): 2223-2247. doi: 10.3934/era.2020113 |
[7] | Qingming Hao, Wei Chen, Zhigang Pan, Chao Zhu, Yanhua Wang . Steady-state bifurcation and regularity of nonlinear Burgers equation with mean value constraint. Electronic Research Archive, 2025, 33(5): 2972-2988. doi: 10.3934/era.2025130 |
[8] | Mingyou Zhang, Qingsong Zhao, Yu Liu, Wenke Li . Finite time blow-up and global existence of solutions for semilinear parabolic equations with nonlinear dynamical boundary condition. Electronic Research Archive, 2020, 28(1): 369-381. doi: 10.3934/era.2020021 |
[9] | Massimo Grossi . On the number of critical points of solutions of semilinear elliptic equations. Electronic Research Archive, 2021, 29(6): 4215-4228. doi: 10.3934/era.2021080 |
[10] | Yuwei Hu, Jun Zheng . Local porosity of the free boundary in a minimum problem. Electronic Research Archive, 2023, 31(9): 5457-5465. doi: 10.3934/era.2023277 |
G-quadruplex or G-tetrad (G4), is a thermodynamically stable structural element that is formed between clusters/stretches/tracts of Guanine (G) residues (|x|≥3) and is intra- or inter-molecular [1,2,3]. The intervening loops whence applicable are composed of one or more nucleotide(s) (N∈{A, U, T, G, C}) (Figure 1). G4 is found in DNA (telomeres, double-strand break sites, transcription start sites) and in the untranslated region(s) (5'-, 3'-UTR, introns) of mRNA [4,5]. In vivo, G4 may function to preserve the telomeric ends of chromosomes, repress or promote transcription and regulate translation [4,5]. The generic representation of an intra-strand G4 may be described as follows:
(((Gt,k)t≥3(Nh,k)h≥1)k=3((Gt,k)t≥3)k=1)m=1 | (Def. 1) |
𝑡 ≔ 𝑁𝑢𝑚𝑏𝑒𝑟 𝑜𝑓 𝐺𝑢𝑎𝑛𝑖𝑛𝑒𝑠 𝑝𝑒𝑟 𝐺 − 𝑟𝑖𝑐ℎ 𝑐𝑙𝑢𝑠𝑡𝑒𝑟
ℎ ≔ 𝑁𝑢𝑚𝑏𝑒𝑟 𝑜𝑓 𝑙𝑜𝑜𝑝 − 𝑓𝑜𝑟𝑚𝑖𝑛𝑔 𝑔𝑒𝑛𝑒𝑟𝑖𝑐 𝑖𝑛𝑡𝑒𝑟𝑣𝑒𝑛𝑖𝑛𝑔 𝑛𝑢𝑐𝑙𝑒𝑜𝑡𝑖𝑑𝑒𝑠
𝑘 ≔ 𝐶𝑙𝑢𝑠𝑡𝑒𝑟 𝑖𝑛𝑑𝑒𝑥
𝑚 ≔ 𝑁𝑢𝑚𝑏𝑒𝑟 𝑜𝑓 𝑠𝑡𝑟𝑎𝑛𝑑𝑠
𝐺 ≔ 𝐺𝑢𝑎𝑛𝑖𝑛𝑒
𝐴 ≔ 𝐴𝑑𝑒𝑛𝑖𝑛𝑒
𝑇 ≔ 𝑇ℎ𝑦𝑚𝑖𝑛𝑒
𝐶 ≔ 𝐶𝑦𝑡𝑜𝑠𝑖𝑛𝑒
𝑁 ≔ 𝐴𝑛𝑦 𝑛𝑢𝑐𝑙𝑒𝑜𝑡𝑖𝑑𝑒
The high melting temperature (Tm~600C) of G4 implies that the mature quadruplex is stable and refractory to unfolding. This is partly due to stabilizing Hoogsteen (N7gu1-N2gu2; O6gu1-N1gu2) and reverse Hoogsteen (N7gu1-N1gu2; O6gu1-N2gu2) hydrogen bonding as well as π-orbital stacking between the purine rings of non-contiguous guanine pairs (gu1, gu2) (Figure 1) [6,7]. Additionally, the presence of Adenine residues in the intervening loops, variable loop length (h~1-30 Mer) and permutation have all been shown to contribute to the stability and thence persistence of the mature quadruplex [8,9,10,11].
Tm∝(#Adenine/h)=τ.(#Adenine/h)Tm:=Melting temperatureτ:=constant of proportionalityh:=Length of intervening loops | (1) |
Despite the wide range of methods available that can predict G4 formation in DNA/RNA, there is poor agreement between sequence-based motif locators and empirically derived biophysical data [12,13,14,15]. Motif-independent methods such as those that directly measure the GC-content or the GC-/AT-skew of a query sequence and utilize this data to train machine learning algorithms may address some of these discrepancies [16,17,18].
Investigations into transcribed RNA suggests that secondary and tertiary forms (5'- and 3'-UTRs) may not only coexist with stretches of unfolded ribonucleotides, but can also be read by the ribosomal machinery. Non-canonical translation is described as: a) translation from atypical start sites AUG→{CUG, GUG} or b) peptides (≤100 aa) of short open reading frames (sORF)-encoded polypeptides (SEPs) and upstream open reading frames (uORFs) [19,20,21,22]. The latter are rarely silent and can function as modulators of metabolism (S-Adenosylmethionine decarboxylase, AMD1) or transcription (activating transcription factor, ATF4, H19; yeast AP-1 like, YAP1) and as generic transcription factors (general control protein, GCN4) [19]. G4 has also been observed in one or more exons of the prion protein (PRNP, exon 2), zinc finger protein (ZNF669, exon 1), β-amyloid secretase (BACE1, exon 3) and the estrogen receptor 1 (ESR1, exon 4) among several others [16,23,24,25,26,27,28,29].
Whilst the presence of segments of folded mRNA may have a significant influence on the yield of the protein product(s), the effect on sequence whence part of the protein coding segment (PCS) is largely unknown [4,5,30,31,32,33]. Proteopathies are diseases that result directly from agammaegates of truncated and misfolded proteins. These may occur secondary to a faulty translation machinery such as a ribosome that has stalled on encountering a secondary or tertiay folded mRNA sub segment. Recent data suggests ~45% of the human genome may code for proteins that are either intrinsically disordered (IDPs) or comprise one or more sub-segments that are disordered (IDRs) [34]. The absence of delineable structural features notwithstanding, disordered regions are characterized by short linear motifs (SLiMS) and/or molecular recognition features (MoRFs) [34,35]. The improper folding and heightened degradation rates could lead to perturbed proteostasis and thence contribute to the pathogenesis of proteopathies [34,35]. Primary proteopathies are likely to result directly from mutations (point, chromosomal translocations) in the PCS of a gene. These include sickle cell disease (βE6→V6-mediated defective polymerization), amylin-based type Ⅱ Diabetes Mellitus, Cystic Fibrosis (cystic fibrosis transmembrane conductance regulator), Alzheimer's disease (Amyloid β-peptide) and Parkinson's disease (α-synuclein) [36,37]. Secondary proteopathies, in contrast, result from motif or molecular mimicry of a host protein(s) by a pathogen. These are further classified into acute and chronic variants depending on the onset, genesis and/or resolution of the resultant infection or infestation [34,35].
G4 is known to stall the ribosome during translation and the resultant protein is truncated and/or degraded at an accelerated rate. The manuscript subsumes ribosomal read-through of mRNA with a G-quadruplex and assesses influence of the translated product to proteostasis. Here, I present a mathematical model of a short G4 (20–60 Mer) in the PCS, i.e., translatable G-quadruplex (TG4), in the mRNA of a hypothetical gene. The mapping uses several novel indices to annotate, classify and select suitable Guanine-containing codons (α) and amino acids (β). A generic algorithm then computes and validates, as proof-of-principle, possible peptides (pTG4ij) that correspond to the modeled TG4 (pTG4ij∈PTG4~TG4). Co-occurrence, homology and the distribution of overlapping/shared amino acids between PTG4 and the disorder promoting SLiMS are used to infer probable mechanisms of TG4~PTG4 facilitated misfolding. Standard bioinformatics indices (accuracy, precision, recall, p-value) are used to arrive at these conclusions.
The objective of this investigation is to model a short G4 in an arbitrary PCS (TG4) which when translated will result in a set of peptides (PTG4) with an average length that is less than 100 amino acids. The hypothesis explored in this manuscript is that in the event of a ribosomal read-through, the translated mRNA, with its G4 will result in a modified protein product. This protein will then exhibit considerable propensity to misfold on account of the presence of one or more members of the PTG4.
2.1.1 Model of a translatable G-quadruplex (TG4)
SEPs-derived peptides with the lowest molecular weight (~2.5 KDa) and with lengths varying from ~7–20 aa were identified and used to define the boundaries of the peptides that comprise PTG4 [20,21]. The TG4 (m = 1) is therefore, modeled as an intra-strand sub sequence of the mRNA of a hypothetical gene and has a length of ~20–60 Mer. This is represented (with symbols and variable names as explained in Def. 1) as follows:
TG4:=(((Gt,k)3≤t≤9(Nh,k)2≤h≤7)k=3((Gt,k)3≤t≤9)k=1)m=1 | (Def. 2) |
Since the Guanine-rich clusters and loops are contiguous, the aforementioned model (Def.2) of the TG4 may be approximated with a sequence of codons and is as under:
TG4:=(CODq)Lq∈N|COD∈COD | (Def. 3) |
The algorithm to compute L, which is the number of codons needed to model TG4 is presented and is as follows:
1:N←{u∈[20,60)}2:r←N mod 33:e←N−((N mod 3)/3)4:If e<(⌊e⌋+⌈e⌉)/2 then5:L=⌈e/3⌉6:else If e≥(⌊e⌋+⌈e⌉)/2 then7:L=⌊e/3⌋8:end If |
𝑁 ≔ 𝑁𝑢𝑚𝑏𝑒𝑟 𝑜𝑓 𝑟𝑖𝑏𝑜𝑛𝑢𝑐𝑙𝑒𝑜𝑡𝑖𝑑𝑒𝑠 𝑟𝑒𝑞𝑢𝑖𝑟𝑒𝑑 𝑡𝑜 𝑚𝑜𝑑𝑒𝑙 𝑇𝐺4
𝐿 : = 𝑇𝑜𝑡𝑎𝑙 𝑛𝑢𝑚𝑏𝑒𝑟 𝑜𝑓 𝑐𝑜𝑑𝑜𝑛𝑠 𝑛𝑒𝑒𝑑𝑒𝑑 𝑡𝑜 𝑚𝑜𝑑𝑒𝑙 𝑇𝐺4 (7 ≤ 𝐿 < 21)
𝑞 ≔ 𝑞𝑡ℎ 𝑐𝑜𝑑𝑜𝑛
𝑟 ≔ 𝑅𝑒𝑚𝑎𝑖𝑛𝑑𝑒𝑟 = {0, 1, 2}
𝑒, 𝑢 : = 𝐺𝑒𝑛𝑒𝑟𝑖𝑐 𝑣𝑎𝑟𝑖𝑎𝑏𝑙𝑒𝑠
𝑪𝑶𝑫 ≔ 𝑆𝑒𝑡 𝑜𝑓 𝑣𝑒𝑟𝑡𝑒𝑏𝑟𝑎𝑡𝑒 𝑐𝑜𝑑𝑜𝑛𝑠
The codons selected for modelling TG4 (Def. 3) comprised suitably scored Guanine-containing vertebrate codons (gCOD+n⊂COD) for the Guanine-rich clusters/stretches/tracks (3≤t≤9;Defs. 1 and 2) and generic/no-stop codons for the intervening loops (Figures 2 and 3). Briefly, a Guanine-containing codon (gCODn) is scored by considering its association with two similar flanking codons, i.e., gCODn-1, gCODn, gCODn+1 such that there is at least one occurrence of 'GGGG' (δ≥1.0) (Figures 2 and 3). This non-trivial case (4≤t≤9) is chosen since its trivial equivalent t = 3, is already subsumed (Defs. 1 and 2). Numerically,
αaminocodon=γ.θ.δ+Ω | (2) |
𝛾 : = 𝑃𝑟𝑜𝑏𝑎𝑏𝑖𝑙𝑖𝑡𝑦 𝑜𝑓 𝑐𝑜𝑑𝑜𝑛 𝑜𝑐𝑐𝑢𝑟𝑟𝑒𝑛𝑐𝑒 (𝛾 = 1/64 ≈ 0.02)
𝜃 ≔ 𝑃𝑟𝑜𝑏𝑎𝑏𝑖𝑙𝑖𝑡𝑦 𝑜𝑓 𝑐𝑜𝑑𝑜𝑛 𝑜𝑐𝑐𝑢𝑟𝑟𝑒𝑛𝑐𝑒 𝑤𝑖𝑡ℎ𝑖𝑛 𝑎 𝑔𝑟𝑜𝑢𝑝 (𝜃 = {0.04, 0.11, 0.33, 1})
𝛿 ≔ 𝐷𝑖𝑠𝑡𝑖𝑛𝑐𝑡 𝑜𝑐𝑐𝑢𝑟𝑟𝑒𝑛𝑐𝑒𝑠 𝑜𝑓 ′𝐺𝐺𝐺𝐺′(𝛿 = {0, 1, 2, 6})
Ω ≔ 𝑁𝑢𝑚𝑏𝑒𝑟 𝑜𝑓 𝑎𝑑𝑗𝑎𝑐𝑒𝑛𝑡 𝑐𝑜𝑑𝑜𝑛𝑠 𝑤𝑖𝑡ℎ 𝛿 (Ω = {0, 1, 2})
Since the genetic code is degenerate, amino acids mapped from the selected codons are further scored and grouped (g1, g2, g3) (Figures 3 and 4).
βamino=|gCOD+amino|/|CODamino| | (6) |
gCOD+amino:= set of optimal codons for each amino acid (αaminocodon>0.0000)
CODamino:= Set of codons for each amino acid
Whilst amino acids from groups 1 and 2 (β > 0.00) (3) can represent the modeled G-rich clusters (y∈g1∪g2 = Y), no constraint was imposed on the amino acids (z) used to model the loops (z∈g1∪g2∪g3 = Z) (Figures 3 and 4). The peptidome (PTG4) evaluated by this study is a combinatorial association of peptides such that the molecular weight is ~0.8–2.3 KDa and length of any arbitrary member is ~7–20 aa (Figure 4). This may be represented as follows:
pTG4ij=((((yi,k)1≤i≤3(zi,k)1≤i≤2)k=3(yi)1≤i≤3)(zi)1≤i≤2)j | (Def. 4) |
PTG4=⋃i=20i=7⋃j=Jj=1|pTG4ij| | (Def. 5) |
𝑷𝑻𝑮𝟒 : = 𝑃𝑒𝑝𝑡𝑖𝑑𝑜𝑚𝑒 𝑐𝑜𝑟𝑟𝑒𝑠𝑝𝑜𝑛𝑑𝑖𝑛𝑔 𝑡𝑜 𝑇𝐺4
𝑝𝑇𝐺4𝑖𝑗 : = 𝑗𝑡ℎ 𝑐𝑎𝑛𝑜𝑛𝑖𝑐𝑎𝑙 𝑎𝑚𝑖𝑛𝑜 𝑎𝑐𝑖𝑑 𝑓𝑜𝑟𝑚 𝑜𝑓 𝑃𝑇𝐺4 𝑤𝑖𝑡ℎ "i" 𝑎𝑚𝑖𝑛𝑜 𝑎𝑐𝑖𝑑𝑠
𝑖 : = 𝑁𝑢𝑚𝑏𝑒𝑟 𝑜𝑓 𝑎𝑚𝑖𝑛𝑜 𝑎𝑐𝑖𝑑𝑠 𝑡ℎ𝑎𝑡 𝑐𝑜𝑚𝑝𝑟𝑖𝑠𝑒 𝑡ℎ𝑒 𝑚𝑜𝑑𝑒𝑙𝑙𝑒𝑑 𝑷𝑻𝑮𝟒
𝐽 : = 𝑀𝑎𝑥𝑖𝑚𝑢𝑚 𝑛𝑢𝑚𝑏𝑒𝑟 𝑜𝑓 𝑐𝑎𝑛𝑜𝑛𝑖𝑐𝑎𝑙 𝑝𝑇𝐺4 𝑓𝑜𝑟 "𝑖" 𝑎𝑚𝑖𝑛𝑜 𝑎𝑐𝑖𝑑𝑠
A dataset that comprises experimentally validated G4-forming mRNA segments of several genes (n = 99) was downloaded (http://scottgroup.med.usherbrooke.ca/G4RNA/) and used to investigate the distribution of G4 [16]. Genes which possess non-redundant RNA (R) sub sequences in the PCS are translated in 6 reading frames using an online tool (http://web.expasy.org/translate). The peptides generated are classified as those: ⅰ) with one or more uninterrupted stretch of N-terminal amino acids of length ≥7 aa (~A), ⅱ) with an in-frame termination signal designated as 'STOP' (~B) and ⅲ) without any termination signal, i.e., absence of a 'STOP' in their sequence (~C). The translated peptides are classified as "VALID" ((B∩A)∪(C∩A)) and then queried for matches with pTG4ij(7≤i≤20,j∈N). The PERL scripts that are required to parse and process the resulting data files have been developed in house and the pseudocode for the same is presented as additional information (Pseudocode, PS1: Supplementary Text 1).
This is done by examining the occurrence of PTG4 in amino acid/protein sequences of disordered regions (IDRs) and full-length proteins with disordered regions (IDPs). DisProt 7.0 (http://disprot.org), is a database of experimentally validated and non-redundant sequences of IDRs and IDPs [38]. The sequences (|IDR| = 1445;|IDP| = 800) that comprise these are queried for occurrences of pTG4ij(7≤i≤20,j∈N) (Supplementary Texts 2 and 3). A preliminary partitioning schema divides these datasets into two distinct subsets, i.e., #pTG4ij≥1 (PT+≡PPOS⊂{IDR, IDP}; (Def.6)) and #pTG4ij = 0 (PT-≡PNEG⊂{IDR, IDP}; (Def. 7)). The extent of co-occurrence of one or more SLiMSw≡SL (w = {1, 2, 3}) with pTG4ij (SL±∈{PPOS, PNEG}) (Defs.8 and 9) is then evaluated to infer relevance of PTG4 to misfolding induced proteostasis. The distribution of overlapping/shared sequences of amino acids ((zn)n≥2∈(pTG4ij∩SLiMSw); zn∈Z; ) (Def.10), is examined in protein sequences from taxonomically diverse organisms with ScanProsite (https://prosite.expasy.org/scanprosite). The proof behind this rationale is presented:
(zn)n≥2 ∈(pTG4ij∩SLiMSw)=((zn)n≥2∈pTG4ij)∩((zn)n≥2∈SLiMSw)Let zn=z′n and zn=z′′n.Rewriting =((z′n)n≥2∈pTG4ij)∩((z′′n)n≥2∈SLiMSw)=((z′n)n≥2,(z′′n)n≥2)=pTG4ij×SLiMSw |
𝒁 : = 𝑆𝑒𝑡 𝑜𝑓 𝑎𝑚𝑖𝑛𝑜 𝑎𝑐𝑖𝑑𝑠 (𝑧𝑛 ∈ 𝒁)
𝑝𝑇𝐺4𝑖𝑗 : = 𝐶𝑎𝑛𝑜𝑛𝑖𝑐𝑎𝑙 𝑎𝑚𝑖𝑛𝑜 𝑎𝑐𝑖𝑑 𝑓𝑜𝑟𝑚 𝑜𝑓 𝑷𝑻𝑮𝟒
𝑺𝑳𝒊𝑴𝑺 : = 𝑆𝑒𝑡 𝑜𝑓 𝑠ℎ𝑜𝑟𝑡 𝑙𝑖𝑛𝑒𝑎𝑟 𝑚𝑜𝑡𝑖𝑓𝑠 (𝑆𝐿𝑖𝑀𝑆𝑤 ∈ 𝑺𝑳𝒊𝑴𝑺)
𝑖, 𝑗, 𝑛, 𝑤 : = 𝐼𝑛𝑑𝑖𝑐𝑒𝑠 𝑜𝑓 𝑚𝑒𝑚𝑏𝑒𝑟𝑠 𝑜𝑓 𝒁, 𝑷𝑻𝑮𝟒, 𝑺𝑳𝒊𝑴𝑺
The indices utilized by this study to establish relevance of matched instances of various motifs/co-motifs in the peptide/protein sequences of interest include the accuracy (A), precision (P), recall (R) and the p-value. A 2X2 table which represents the categorized data (2.1.4) is constructed and used to compute various bioinformatics indices. This is outlined as under:
![]() |
𝑇𝑁 ≔ 𝑇𝑟𝑢𝑒 𝑛𝑒𝑔𝑎𝑡𝑖𝑣𝑒 (|𝒔𝑵𝑬𝑮 ∩ 𝑷𝑵𝑬𝑮| = |𝒔𝑵𝑬𝑮|) ≡ 𝑆𝐿−𝑃𝑇− (𝐷𝑒𝑓. 11)
𝐹𝑃 ≔ 𝐹𝑎𝑙𝑠𝑒 𝑝𝑜𝑠𝑖𝑡𝑖𝑣𝑒 (|𝑺𝑵𝑬𝑮 ∩ 𝑷𝑷𝑶𝑺| = |𝑺𝑵𝑬𝑮|) ≡ 𝑆𝐿−𝑃𝑇+ (𝐷𝑒𝑓. 12)
𝐹𝑁 ≔ 𝐹𝑎𝑙𝑠𝑒 𝑛𝑒𝑔𝑎𝑡𝑖𝑣𝑒 (|𝒔𝑷𝑶𝑺 ∩ 𝑷𝑵𝑬𝑮| = |𝒔𝑷𝑶𝑺|) ≡ 𝑆𝐿+𝑃𝑇− (𝐷𝑒𝑓. 13)
𝑇𝑃 ≔ 𝑇𝑟𝑢𝑒 𝑝𝑜𝑠𝑖𝑡𝑖𝑣𝑒 (|𝑺𝑷𝑶𝑺 ∩ 𝑷𝑷𝑶𝑺| = |𝑺𝑷𝑶𝑺|) ≡ 𝑆𝐿+𝑃𝑇+ (𝐷𝑒𝑓. 14)
The equations may then be written as:
(A)=(TN+TP/TN+FP+FN+TP)X100 | (4) |
(P)=(TP/FP+TP)X100 | (5) |
(R)=(TP/FN+TP)X100 | (6) |
The p-values for these analyses are computed by comparing the frequency of occurrence of all pTG4ij in a test sequence (ϕpTG4ij) with the same in randomly-generated (v∈V) sequences of similar lengths (ϕpTG4vij), i.e., 7-50 aa (1≤v≤10000) and > 50 aa (1≤v≤100000) (Pseudocode, PS2: Supplementary Text 1):
p−value=ϕpTG4vij/ϕpTG4ij=(v=|V|∑v=1i=21∑i=7j=J∑j=1pTG4vij)/(i=21∑i=7j=J∑j=1pTG4ij)=(v=|V|∑v=1i=21∑i=7j=J∑j=1pTG4vij/i=21∑i=7j=J∑j=1pTG4ij) | (7) |
The frequency of occurrence of overlapping sequences of amino acids ((zn)n≥2∈(pTG4ij∩SLiMSw); zn∈Z) in pre-compiled and curated protein sequences (ϕ(zn)) across taxa is compared with randomly chosen sequences of comparable lengths (ϕ(vzn); n = 5000). These are used to estimate statistical significance, i.e., p-value = ϕ(vzn)/ϕ(zn) (8).
The data presented discusses implementation of a model of short intra-strand TG4 for various values of α and β, populates PTG4 and establishes the equivalence TG4~PTG4. Co-occurrence and homology studies between PTG4 and the SLiMS in IDRs/IDPs and generic protein sequences across taxa are used to infer probable mechanisms of TG4~PTG4 facilitated misfolding-induced proteostasis.
An association-competent codon not only takes into account the presence of a Guanine residue, but also gives weightage to its position (Figures 1–3, Table 1). This schema partitions standard vertebrate codons into those with a high- (Ranks 1-4;α > 0.0000) or low- (Rank 5;α = 0.0000) propensity to form a contiguous cluster of Guanine residues (Figures 3 and 4, Table 1). Whilst, 'GGG' (Rank 1;α = 2.12) can associate with ({GGG, GxG, xGG, GGx, Gxx, xxG}) bilaterally (δ = 6;Ω = 2), 'GxG' (Rank 2; α = 2.0066) can do so only with 'GGG' (δ = 1;Ω = 2). On the other hand, the codon subsets 'GGx' and 'xGG' (Rank 3;α = 1.0132) can form two clusters of contiguous Guanine residues with 'GGG' and 'xGG'/ 'GGx' unilaterally (δ = 2;Ω = 1). Similarly, the subsets 'xxG' or 'Gxx' (Rank 4;α = 1.0022), can form contiguous Guanines with a single occurrence of 'GGG' (δ = 1;Ω = 1) (Figures 3 and 4, Table 1). Conversely, codons with either a single occurrence of a central Guanine residue 'xGx' or no Guanine residues 'xxx' (Rank 5;α = 0.0000) are unable to form the 'GGGG' and are excluded from this study (Figures 3 and 4, Table 1).
Rank | Codon set, Cardinality | Codon | γ | θ | δ | Ω | α=γ.θ.δ+Ω | aa |
1 | GGG, 1 | GGG | 0.02 | 1.00 | 6 | 2 | 2.1200 | Gly |
2 | GxG, 3 | GUG | 0.02 | 0.33 | 1 | 2 | 2.0066 | Val |
GCG | 0.02 | 0.33 | 1 | 2 | 2.0066 | Ala | ||
GAG | 0.02 | 0.33 | 1 | 2 | 2.0066 | Glu | ||
3 | xGG, 3 | UGG | 0.02 | 0.33 | 2 | 1 | 1.0132 | Trp |
CGG | 0.02 | 0.33 | 2 | 1 | 1.0132 | Arg | ||
AGG | 0.02 | 0.33 | 2 | 1 | 1.0132 | Arg | ||
3 | GGx, 3 | GGU | 0.02 | 0.33 | 2 | 1 | 1.0132 | Gly |
GGC | 0.02 | 0.33 | 2 | 1 | 1.0132 | Gly | ||
GGA | 0.02 | 0.33 | 2 | 1 | 1.0132 | Gly | ||
4 | xxG, 9 | UUG | 0.02 | 0.11 | 1 | 1 | 1.0022 | Leu |
UCG | 0.02 | 0.11 | 1 | 1 | 1.0022 | Ser | ||
UAG | 0.02 | 0.11 | 1 | 1 | 1.0022 | Ter | ||
CUG | 0.02 | 0.11 | 1 | 1 | 1.0022 | Leu | ||
CCG | 0.02 | 0.11 | 1 | 1 | 1.0022 | Pro | ||
CAG | 0.02 | 0.11 | 1 | 1 | 1.0022 | Gln | ||
AUG | 0.02 | 0.11 | 1 | 1 | 1.0022 | Met | ||
ACG | 0.02 | 0.11 | 1 | 1 | 1.0022 | Thr | ||
AAG | 0.02 | 0.11 | 1 | 1 | 1.0022 | Lys | ||
4 | Gxx, 9 | GUU | 0.02 | 0.11 | 1 | 1 | 1.0022 | Val |
GCU | 0.02 | 0.11 | 1 | 1 | 1.0022 | Ala | ||
GAU | 0.02 | 0.11 | 1 | 1 | 1.0022 | Asp | ||
GUC | 0.02 | 0.11 | 1 | 1 | 1.0022 | Val | ||
GCC | 0.02 | 0.11 | 1 | 1 | 1.0022 | Ala | ||
GAC | 0.02 | 0.11 | 1 | 1 | 1.0022 | Asp | ||
GUA | 0.02 | 0.11 | 1 | 1 | 1.0022 | Val | ||
GCA | 0.02 | 0.11 | 1 | 1 | 1.0022 | Ala | ||
GAA | 0.02 | 0.11 | 1 | 1 | 1.0022 | Glu | ||
5 | xGx, 9 | UGU | 0.02 | 0.11 | 0 | 0 | 0.0000 | Cys |
UGC | 0.02 | 0.11 | 0 | 0 | 0.0000 | Cys | ||
UGA | 0.02 | 0.11 | 0 | 0 | 0.0000 | Ter | ||
CGU | 0.02 | 0.11 | 0 | 0 | 0.0000 | Arg | ||
CGC | 0.02 | 0.11 | 0 | 0 | 0.0000 | Arg | ||
CGA | 0.02 | 0.11 | 0 | 0 | 0.0000 | Arg | ||
AGU | 0.02 | 0.11 | 0 | 0 | 0.0000 | Ser | ||
AGC | 0.02 | 0.11 | 0 | 0 | 0.0000 | Ser | ||
AGA | 0.02 | 0.11 | 0 | 0 | 0.0000 | Arg | ||
5 | xxx, 27 | UUU | 0.02 | 0.04 | 0 | 0 | 0.0000 | Phe |
UCU | 0.02 | 0.04 | 0 | 0 | 0.0000 | Ser | ||
UAU | 0.02 | 0.04 | 0 | 0 | 0.0000 | Tyr | ||
UUC | 0.02 | 0.04 | 0 | 0 | 0.0000 | Phe | ||
UCC | 0.02 | 0.04 | 0 | 0 | 0.0000 | Ser | ||
UUA | 0.02 | 0.04 | 0 | 0 | 0.0000 | Leu | ||
UCA | 0.02 | 0.04 | 0 | 0 | 0.0000 | Ser | ||
UAA | 0.02 | 0.04 | 0 | 0 | 0.0000 | Ter | ||
CUU | 0.02 | 0.04 | 0 | 0 | 0.0000 | Leu | ||
CCU | 0.02 | 0.04 | 0 | 0 | 0.0000 | Pro | ||
CAU | 0.02 | 0.04 | 0 | 0 | 0.0000 | His | ||
CUC | 0.02 | 0.04 | 0 | 0 | 0.0000 | Leu | ||
CCC | 0.02 | 0.04 | 0 | 0 | 0.0000 | Pro | ||
CAC | 0.02 | 0.04 | 0 | 0 | 0.0000 | His | ||
CUA | 0.02 | 0.04 | 0 | 0 | 0.0000 | Leu | ||
CCA | 0.02 | 0.04 | 0 | 0 | 0.0000 | Pro | ||
CAA | 0.02 | 0.04 | 0 | 0 | 0.0000 | Gln | ||
AUU | 0.02 | 0.04 | 0 | 0 | 0.0000 | Ile | ||
ACU | 0.02 | 0.04 | 0 | 0 | 0.0000 | Thr | ||
AAU | 0.02 | 0.04 | 0 | 0 | 0.0000 | Asn | ||
AUC | 0.02 | 0.04 | 0 | 0 | 0.0000 | Ile | ||
ACC | 0.02 | 0.04 | 0 | 0 | 0.0000 | Thr | ||
AAC | 0.02 | 0.04 | 0 | 0 | 0.0000 | Asn | ||
AUA | 0.02 | 0.04 | 0 | 0 | 0.0000 | Ile | ||
ACA | 0.02 | 0.04 | 0 | 0 | 0.0000 | Thr | ||
AAA | 0.02 | 0.04 | 0 | 0 | 0.0000 | Lys |
Abbreviations
𝛾: General probability of a codon (𝛾 = 1/64 ≅ 0.02)
𝜃: Probability of codon within a group (𝜃 = {0.04, 0.33, 0.11, 1.00})
𝛿: Number of distinct codon sets that could complete 'GGGG' (𝛿 = {0, 1, 2, 6})
Ω: Number of adjacent positions that contain 𝛿 (Ω = {0, 1, 2})
𝛼: Threshold for selecting codons that may favour G-quadruplex formation
x: Codon specific generic ribonucleotide {𝐴, 𝐺, 𝑈, 𝐶}
aa: Amino acid
Ter: Stop codons {𝑈𝐴𝐺, 𝑈𝐺𝐴, 𝑈𝐴𝐴}
An estimate of the possible combinations of the simplest peptide (∑i=7∑j=Jj=1pTG4ij = 8.00E+03;GlyzGlyzGlyzGly; J = (20)3; i = lengthp(TG4ij) = 7 aa; z∈Z) (Figures 3 and 4, Table 2). This justifies usage of PTG4 (pTG4ij∈PTG4) as a generic representation of the putative peptidome encoded by the TG4 (PTG4). Approximately ~12% (n = 11) of in silico translated amino acid sequences from exon-derived TG4 possesses one of more "STOP" signals and include ESR1, longer RNA variants of PRNP (85 nt) and BCL2 (29 n t, 33 nt, 34 nt) (Table 3; Table 1, Supplementary Text 2). With the exceptions of KCNH2/ZNF669 and the shorter variants of PRNP (14 nt, 15 nt, 20 nt, 24 nt), "VALID" sub sequences are found for BACE1, BCL2, ESR1, PRNP (long) and TERF2 (Table 3; Tables 1A and 1C). Interestingly, all the genes considered possessed at least one occurrence of PTG4 (P = 100%, n = 6) (Table 3; Table 1B). This finding, despite the small sample size is proof-of-principal that the TG4 can be mapped to definite peptide sequences, i.e., TG4~PTG4. Since this can occur only after a ribosomal read through of the G4 containing mRNA, it raises the intriguing possibility that PTG4 whence part of a larger protein may increase its propensity to undergo misfolding. This notion is investigated in non-redundant sequences of IDRs (PTG4~10%, n = 145;0.00≤p-value≤0.20) and IDPs (PTG4~34%, n = 269;0.00≤p-value < 0.5) (Table 4; Tables 2 and 3).
aa | CODamino | gCOD+amino | β | |
Group 1 (n=7) | Ala | 4 | 4 | 1.00 |
Val | 4 | 4 | 1.00 | |
Asp | 2 | 2 | 1.00 | |
Glu | 2 | 2 | 1.00 | |
Trp | 1 | 1 | 1.00 | |
Met | 1 | 1 | 1.00 | |
Gly | 4 | 4 | 1.00 | |
Group 2 (n=7) | Leu | 6 | 2 | 0.3333 |
Gln | 2 | 1 | 0.5 | |
Arg | 6 | 2 | 0.3333 | |
Lys | 2 | 1 | 0.5 | |
Ser | 6 | 1 | 0.1667 | |
Thr | 4 | 1 | 0.25 | |
Pro | 4 | 1 | 0.25 | |
Group 3 (n=6) | Cys | 2 | 0 | 0.00 |
Asn | 2 | 0 | 0.00 | |
Ile | 3 | 0 | 0.00 | |
His | 2 | 0 | 0.00 | |
Phe | 2 | 0 | 0.00 | |
Tyr | 2 | 0 | 0.00 |
Abbreviations
gCOD+amino : Guanine-containing optimal codons excluding STOP (UAG) (𝛼 > 0.000)
COD−amino : Non-optimal codon excluding STOP (UGA, UAA) (𝛼 = 0.000)
𝑪𝑶𝑫𝒂𝒎𝒊𝒏𝒐 = gCOD+amino + COD−amino : All codons for an amino acid
GENE | NAME | G4 (nt) | Ex | STOP(n=11) | VALID(n=59) | |PTG4| |
BACE1 | Beta-secretase 1 | 33 | 3 | n=0 | n=6 | n=2 |
BCL2 | B-cell lymphoma 2 | 33 | 2 | n=1 | n=6 | n=1 |
23 | n=0 | n=6 | ||||
28 | n=0 | n=6 | ||||
29 | n=1 | n=5 | ||||
34 | n=1 | n=5 | ||||
33 | 3 | n=2 | n=5 | |||
ESR1 | Estrogen receptor alpha (ERα) | 36 | 4 | n=1 | n=5 | n=2 |
KCNH2 | Potassium Voltage-Gated Channel sub family H | 18 | 12 | n=0 | n=0 | NA |
ZNF669 | Member 2 Zinc Finger Protein 669 |
1 | ||||
PRNP | Prion protein | 14 | 2 | n=0 | n=0 | n=1 |
15 | n=0 | n=0 | ||||
20 | n=0 | n=0 | ||||
24 | n=0 | n=6 | ||||
85 | n=6 | n=3 | ||||
TERF2 | Telomeric repeat-binding factor 2 | 55 | 1 | n=0 | n=6 | n=1 |
Disordered regions (IDRs; n=1445;0.00≤p-value < 0.05) | |||||||||||||||||
SL-PT- | SL-PT+ | SL+PT- | SL+PT+ | R1T | R2T | C1T | C2T | A (%) | P (%) | R (%) | |||||||
SLiMS1 | 1078 | 64 | 58 | 9 | 1142 | 67 | 1136 | 73 | 89.90 | 12.32 | 13.43 | ||||||
SLiMS2 | 749 | 18 | 121 | 29 | 767 | 150 | 870 | 47 | 84.84 | 61.70 | 19.33 | ||||||
SLiMS3 | 1212 | 108 | 34 | 9 | 1320 | 43 | 1246 | 117 | 89.58 | 7.69 | 20.93 | ||||||
Proteins with disordered segments (IDPs; n=800;0.00≤p-value < 0.05) | |||||||||||||||||
SL-PT- | SL-PT+ | SL+PT- | SL+PT+ | R1T | R2T | C1T | C2T | A (%) | P (%) | R (%) | |||||||
SLiMS1 | 86 | 12 | 28 | 18 | 98 | 46 | 114 | 30 | 72.22 | 60.00 | 39.10 | ||||||
SLiMS2 | 1 | 1 | 96 | 66 | 2 | 162 | 97 | 67 | 40.85 | 98.50 | 40.74 | ||||||
SLiMS3 | 250 | 57 | 26 | 25 | 307 | 51 | 276 | 82 | 76.81 | 30.48 | 49.01 |
Abbreviations
𝐼𝐷𝑅𝑠: Intrinsically disordered regions
𝐼𝐷𝑃s: Intrinsically disordered proteins
𝑧: Any amino acid
𝑆𝐿𝑖𝑀𝑆1: [𝑆𝑇]𝑃𝑧𝑅
𝑆𝐿𝑖𝑀𝑆2: [𝐸𝐷]𝑧𝑧[𝐷𝐸][𝐴𝐺𝑆]
𝑆𝐿𝑖𝑀𝑆3: [𝐾𝑅]𝑧𝑃𝑧𝑧𝑃
𝑆𝐿−𝑃𝑇−: |𝑺𝑵𝑬𝑮 ∩ 𝑷𝑵𝑬𝑮|
𝑆𝐿−𝑃𝑇+: |𝑺𝑵𝑬𝑮 ∩ 𝑷𝑷𝑶𝑺|
𝑆𝐿+𝑃𝑇−: |𝑺𝑷𝑶𝑺 ∩ 𝑷𝑵𝑬𝑮|
𝑆𝐿+𝑃𝑇+: |𝑺𝑷𝑶𝑺 ∩ 𝑷𝑷𝑶𝑺|
𝑅1𝑇: 𝑆𝐿−𝑃𝑇− + 𝑆𝐿−𝑃𝑇+
𝑅2𝑇: 𝑆𝐿+𝑃𝑇− + 𝑆𝐿+𝑃𝑇+
𝐶1𝑇: 𝑆𝐿−𝑃𝑇−+ 𝑆𝐿+𝑃𝑇−
𝐶2𝑇: 𝑆𝐿−𝑃𝑇+ + 𝑆𝐿+𝑃𝑇+
𝐴: Accuracy
𝑃: Precision
𝑅: Recall
The amino acids that comprise the peptide members of PTG4 and the short linear motifs (g1, g2 vs SLiMS) are well conserved. The co-occurrence of PTG4 with SLiMS in the IDRs (A~85-89%;0.00 < p-value≤0.05) suggests that this association is non-trivial and may favor all purported mechanisms of misfolding (hyperphosphorylation, proteolytic cleavage, complex formation) (Table 4; Tables 2 and 3). However, the higher precision of PTG4 with the proteolytic-SLiMS suggests that this may predominate (Table 4; Tables 2 and 3). The data with the IDPs suggests a similar predilection for proteolytic cleavage (A~40-77%;P~99%;0.00 < p-value < 0.05, although hyperphosphorylation (P~60%;0.00 < p-value < 0.05) and complex-promotion (P~30%;0.00 < p-value < 0.05) may constitute viable alternatives to the genesis of misfolding (Table 4; Tables 2 and 3). The presence of overlapping sequences of amino acids between PTG4 and the SLiMS when examined in protein sequences from taxonomically diverse organisms is degenerate for SLiMS1 (number of matches = 6251) and SLiMS3 (number of matches = 1480) (Table 5; Table 4). In contrast, the corresponding data for SLiMS2 (number ofmatches = 3759;0.00 < p-value < 0.05) is statistically significant (Table 5). The taxonomic spread includes archaea (n = 150), bacteria (n = 1735), viruses (n = 84), green land plants (n = 199), fungi (n = 182), eukaryotic invertebrates (n = 43) and vertebrates (n = 700) (Table 4).
SLiMS | Sample | (zn)n≥2∈pTG4ij∩SLiMSw (p-value) | |
SLiMS1=[ST]PzR | Pz | PG (n=1) | PG (Degenerate) |
SLiMS2=[ED]zzD[AGS] | z[DE] | G[DE] (n=2) | [WGRVAELMKQSTP][AG][DE]z2EG[VADE](p-value=0.00069) |
[DE]z | [DE]G (n=1) | ||
zz[DE] | [LMKQSTP]G[DE] (n=14) [VAE]G[DE] (n=6) [WGR][AG][DE] (n=6) |
||
[DE]zz[DE] | [VAE]G[DE]zzEG[VADE] (n=28) [WGR][AG][DE]zzEG[VADE] (n=24) [WGR][AG][DE]zzEG (n=6) GEzzEG[VADE] (n=4) GEzzEG (n=1) |
||
SLiMS3=[KR]zPzzP | Pzz | PG[VADE] (n=4) | PGV (Degenerate) PGA (Degenerate) PGD (Degenerate) PGE (Degenerate) |
Abbreviations
𝑝𝑇𝐺4𝑖𝑗: Members of putative peptidome (𝑝𝑇𝐺4𝑖𝑗 ∈ 𝑷𝑻𝑮𝟒)
𝑆𝐿𝑖𝑀𝑆𝑤: Short linear motifs (𝑆𝐿𝑖𝑀𝑆𝑤 ∈ 𝑺𝑳𝒊𝑴𝑺)
𝑧𝑛: Shared sequence(s) of amino acids between 𝑷𝑻𝑮𝟒 and 𝑺𝑳𝒊𝑴𝑺
𝑖, 𝑗, 𝑤, 𝑛: Indices to characterize members of 𝑷𝑻𝑮𝟒, 𝑺𝑳𝒊𝑴𝑺, 𝒁
The significant association and homology between PTG4 and the SLiMS along with the equivalence data (PTG4~TG4) suggest that TG4 may influence proteostasis in a multitude of ways (Tables 1–5; Tables 1–4, Supplementary Text 2–4).
The short TG4 modeled in this study has an average loop length (h~2 Mer) which may contribute to thermodynamic stability by restricting the mobility of the participating strands (1) [8,9,10,11]. The physical presence of TG4 will result in a stalled ribosome and translation which is prolonged, inefficient and incomplete [31,32,33]. Interestingly, this analysis also includes UAG (Amber; α > 0.0000), which when present in-frame will prematurely terminate translation and result in a truncated protein (Table 1) [39]. Whilst nonsense-mediated mRNA decay may be triggered if the stop codon is within ±50 Mer of the exon-junction complex (EJC), a read-through may occur nonetheless. The resulting protein sequences may be modified which in tandem with one or more occurrences of PTG4 and/or SLiMS would predispose the same to agammaegate and result in a proteopathy [39,40].
Whilst the preponderance of Glycine (Gly) might impart heightened flexibility and limit the formation of stabilizing secondary structural elements in the hypothetical protein, Proline (Pro) confers rigidity and may retard proper folding. There is also remarkable conservation between the amino acids that comprise PTG4 and the SLiMS. These include the complex-promoting hydrophobic (Ala, Val, Met, Trp) and ionic (Asp, Glu, Lys, Arg) residues, along with nucleophile-favoring Serine and Threonine (Figures 3 and 4, Tables 2–5). Whilst, the former may favor agammaegation by non-covalent interactions, the latter may promote phosphorylation-mediated charge imbalance and thence misfolding. Interestingly, the loops of G4 whence modeled by Adenine-containing codons (Axx) are translated to Lysine (K), Arginine (R), Serine (S), Threonine (T) and Isoleucine (I); all of which may also promote misfolding (Figures 3 and 4, Tables 2-5) [8,9,10,11,34,35]. The distribution of PTG4 amongst physiologically relevant proteins further suggests that the peptide-mediated misfolding may influence/regulate signal transduction, cytoskeleton organization, metabolism, synaptic transmission and transcription/translation (Table 6; Table 5).
Cellular function | Disordered regions of proteins | |
1. | Signal transduction | DP00274, DP00224, DP00141, DP00332, DP01063, DP00506, DP00418, DP00341, DP00435, DP00613, DP00463, DP00954, DP00959, DP01104, DP00611, DP00519, DP00086, DP00707, DP00712 |
2. | Endocytosis | DP01073, DP01065, DP01066, DP00225 |
3. | Calcium-calmodulin | DP00092, DP00132, DP00561, DP00118, DP00253 |
4. | Myofibril assembly | DP01090 |
5. | Cytoskeleton | DP01056 DP00240, DP01022, DP00169, DP00716, DP00717, DP01100, DP00122 |
6. | Nuclear pore | DP01075, DP01077, DP01079 |
7. | Phototransduction | DP00768, DP00347 |
8. | Targeting | DP00893, DP00609, DP00610, DP01058 |
9. | Transcription | DP00062, DP00177, DP00633, DP00348, DP00786, DP00049, DP00231, DP00873, DP00720, DP00217, DP00081 |
10. | Translation | DP00082, DP00164, DP00229 DP00949, DP00134 |
11. | Synaptic transmission | DP00943 |
12. | Supercoiling | DP00076 |
13. | Binding | DP00539, DP00854, DP01052, DP00659, DP00656 |
14. | Peptide bond formation | DP00944 |
15. | Enzymes | DP00557, DP00032, DP00095, DP00337, DP00379, DP00787, DP00427, DP00429 |
16. | Bacterial/parasitic virulence | |
Secreted toxins | DP00345, DP00591 | |
Cytoadherence | DP00025, DP00065, DP01096 | |
17. | Viral infectivity | |
Cyclophilin interaction | DP00615, DP01031 | |
Chaperones | DP00699, DP00700, DP00674 | |
Capsid assembly | DP00133, DP00876 | |
Membrane fusion | DP01043 | |
Latency | DP01060 | |
18. | Unknown | DP00119 |
Note: DP≔DisProt ID
The distribution of overlapping/shared amino acids in protein sequences of non-vertebrates suggests that PTG4 is either completely degenerate with the SLiMS or present in proportions that is statistically significant (Tables 5 and 6; Tables 4 and 5). These data imply that motif-mimicry too, might constitute a probable cause (tropism, oncogenic potential, virulence) of infection/infestation-mediated acute/chronic proteopathies [34,35,41,42]. The contribution(s) of misfolding to the pathogenesis of secondary proteopathies is however, debatable. Whilst, there is evidence that mislocalization of proteins can precipitate misfolding, mimicry itself may result in exonuclease-mediated proteolytic cleavage and thence trigger an infective proteopathy [43,44]. Additionally, the presence of sequences of amino acids such as Proline and Threonine in viral or fungal proteins may be responsible for creating and/or maintaining a milieu conducive to the genesis of infective/transmissible proteopathies, viz., a high charge density and imbalance of electrostatic interactions [43,44].
The coexistence of potentially translatable G-quadruplexes (TG4) with unfolded ribonucleotides in the PCS of an mRNA transcript may have important consequences for protein homeostasis. Here, I have investigated the contribution of a short intra-strand translatable G-quadruplex and its associated peptidome (TG4~PTG4) to the genesis of misfolding-induced proteostasis. The co-occurrence, homology and distribution of overlapping/shared amino acids of PTG4 with the SLiMS suggests that this may occur by truncation, complex formation, increased charge density and/or accelerated degradation. An additional mechanism that is also supported is motif-mimicry by pathogens which may trigger the development of infective proteopathies. The putative peptidome (~7–20 aa) that corresponds to the short translatable G-quadruplex delineated by this investigation may be utilized as novel markers of both the primary and secondary proteopathies.
SK outlined and designed the study, designed and conceptualized the algorithm(s) and formulae for prediction, wrote mathematical proofs to establish rigor, collated the data, constructed the models, formulated the filters, carried out the computational analysis, wrote all necessary code and the manuscript.
The author declares no conflict of interest.
[1] |
Engelkamp H, Middelbeek S, Nolte JMR, et al. (1999) Self-Assembly of Disk-Shaped Molecules to Coiled-Coil Aggregates with Tunable Helicity. Science 284: 785–788. doi: 10.1126/science.284.5415.785
![]() |
[2] |
Mondal T, Basak D, Al Ouahabi A, et al. (2015) Extended supramolecular organization of [small pi]-systems using yet unexplored simultaneous intra- and inter-molecular H-bonding motifs of 1,3-dihydroxy derivatives. Chem Commun 51: 5040–5043. doi: 10.1039/C4CC10335A
![]() |
[3] |
Lu H, Kobayashi N (2016) Optically Active Porphyrin and Phthalocyanine Systems. Chem Rev 116: 6184–6261. doi: 10.1021/acs.chemrev.5b00588
![]() |
[4] |
Volpati D, Alessio P, Zanfolim AA, et al. (2008) Exploiting Distinct Molecular Architectures of Ultrathin Films Made with Iron Phthalocyanine for Sensing. J Phys Chem B 112: 15275–15282. doi: 10.1021/jp804159h
![]() |
[5] |
Dey S, Pal AJ (2011) Layer-by-Layer Electrostatic Assembly with a Control over Orientation of Molecules: Anisotropy of Electrical Conductivity and Dielectric Properties. Langmuir 27: 8687–8693. doi: 10.1021/la201471f
![]() |
[6] |
Eccher J, Zajaczkowski W, Faria GC, et al. (2015) Thermal Evaporation versus Spin-Coating: Electrical Performance in Columnar Liquid Crystal OLEDs. ACS Appl Mater Interface 7: 16374–16381. doi: 10.1021/acsami.5b03496
![]() |
[7] | Cea P, Ballesteros Luz M, Martín S (2014) Nanofabrication techniques of highly organized monolayers sandwiched between two electrodes for molecular electronics. Nanofabrication. |
[8] |
Camacho SA, Aoki PHB, Assis FFd, et al. (2014) Supramolecular arrangements of an organometallic forming nanostructured films. Mater Res 17: 1375–1383. doi: 10.1590/1516-1439.279014
![]() |
[9] |
Yang Y, Zhang Y, Wei Z (2013) Supramolecular Helices: Chirality Transfer from Conjugated Molecules to Structures. Adv Mater 25: 6039–6049. doi: 10.1002/adma.201302448
![]() |
[10] |
Li W-S, Aida T (2009) Dendrimer Porphyrins and Phthalocyanines. Chem Rev 109: 6047–6076. doi: 10.1021/cr900186c
![]() |
[11] | Chen Z, Lohr A, Saha-Moller CR, et al. (2009) Self-assembled [small pi]-stacks of functional dyes in solution: structural and thermodynamic features. Chem Soc Rev 38: 564–584. |
[12] | Wurthner F (2004) Perylene bisimide dyes as versatile building blocks for functional supramolecular architectures. Chem Commun 1564–1579. |
[13] | Würthner F, Thalacker C, Diele S, et al. (2001) Fluorescent J-type Aggregates and Thermotropic Columnar Mesophases of Perylene Bisimide Dyes. Chemistry–A Eur J 7: 2245–2253. |
[14] |
Voitechovič E, Bratov A, Abramova N, et al. (2015) Development of label-free impedimetric platform based on new conductive polyaniline polymer and three-dimensional interdigitated electrode array for biosensor applications. Electrochim Acta 173: 59–66. doi: 10.1016/j.electacta.2015.05.011
![]() |
[15] |
Bisquert J, Garcia-Belmonte G, Bueno P, et al. (1998) Impedance of constant phase element (CPE)-blocked diffusion in film electrodes. J Electroanal Chem 452: 229–234. doi: 10.1016/S0022-0728(98)00115-6
![]() |
[16] |
Yang L (2008) Electrical impedance spectroscopy for detection of bacterial cells in suspensions using interdigitated microelectrodes. Talanta 74: 1621–1629. doi: 10.1016/j.talanta.2007.10.018
![]() |
[17] |
Bertoncello P, Peruffo M (2008) An investigation on the self-aggregation properties of sulfonated copper(II) phthalocyanine (CuTsPc) thin films. Colloid Surface A 321: 106–112. doi: 10.1016/j.colsurfa.2008.01.054
![]() |
[18] |
Lee SK, Zu Y, Herrmann A, et al. (1999) Electrochemistry, Spectroscopy and Electrogenerated Chemiluminescence of Perylene, Terrylene, and Quaterrylene Diimides in Aprotic Solution. J Am Chem Soc 121: 3513–3520. doi: 10.1021/ja984188m
![]() |
[19] |
Arıcı M, Arıcan D, Uğur AL, et al. (2013) Electrochemical and spectroelectrochemical characterization of newly synthesized manganese, cobalt, iron and copper phthalocyanines. Electrochim Acta 87: 554–566. doi: 10.1016/j.electacta.2012.09.045
![]() |
[20] |
Apetrei C, Alessio P, Constantino CJL, et al. (2011) Biomimetic biosensor based on lipidic layers containing tyrosinase and lutetium bisphthalocyanine for the detection of antioxidants. Biosens Bioelectron 26: 2513–2519. doi: 10.1016/j.bios.2010.10.047
![]() |
[21] |
Alessio P, Martin CS, de Saja JA, et al. (2016) Mimetic biosensors composed by layer-by-layer films of phospholipid, phthalocyanine and silver nanoparticles to polyphenol detection. Sensor Actuat B-Chem 233: 654–666. doi: 10.1016/j.snb.2016.04.139
![]() |
[22] |
Song Y, Yang T, Zhou X, et al. (2016) A microsensor for hydroquinone and catechol based on a poly(3,4-ethylenedioxythiophene) modified carbon fiber electrode. Anal Method 8: 886–892. doi: 10.1039/C5AY02532J
![]() |
[23] | García-Hernández C, García-Cabezón C, Medina-Plaza C, et al. (2015) Electrochemical behavior of polypyrrol/AuNP composites deposited by different electrochemical methods: sensing properties towards catechol. Beilstein J Nanotechno 6: 2052–2061. |
[24] |
Liu W, Wu L, Zhang X, et al. (2014) Highly-selective electrochemical determination of catechol based on 3-aminophenylboronic acid-3,4,9,10-perylene tetracarboxylic acid functionalized carbon nanotubes modified electrode. Anal Method 6: 718–724. doi: 10.1039/C3AY41633J
![]() |
[25] |
Janeiro P, Oliveira Brett AM (2004) Catechin electrochemical oxidation mechanisms. Anal Chim Acta 518: 109–115. doi: 10.1016/j.aca.2004.05.038
![]() |
1. | Zhihua Huang, Jianfu Yang, Weilin Yu, Boundary plasmas for a confined plasma problem in dimensional two, 2023, 62, 0944-2669, 10.1007/s00526-022-02421-2 |
Rank | Codon set, Cardinality | Codon | γ | θ | δ | Ω | α=γ.θ.δ+Ω | aa |
1 | GGG, 1 | GGG | 0.02 | 1.00 | 6 | 2 | 2.1200 | Gly |
2 | GxG, 3 | GUG | 0.02 | 0.33 | 1 | 2 | 2.0066 | Val |
GCG | 0.02 | 0.33 | 1 | 2 | 2.0066 | Ala | ||
GAG | 0.02 | 0.33 | 1 | 2 | 2.0066 | Glu | ||
3 | xGG, 3 | UGG | 0.02 | 0.33 | 2 | 1 | 1.0132 | Trp |
CGG | 0.02 | 0.33 | 2 | 1 | 1.0132 | Arg | ||
AGG | 0.02 | 0.33 | 2 | 1 | 1.0132 | Arg | ||
3 | GGx, 3 | GGU | 0.02 | 0.33 | 2 | 1 | 1.0132 | Gly |
GGC | 0.02 | 0.33 | 2 | 1 | 1.0132 | Gly | ||
GGA | 0.02 | 0.33 | 2 | 1 | 1.0132 | Gly | ||
4 | xxG, 9 | UUG | 0.02 | 0.11 | 1 | 1 | 1.0022 | Leu |
UCG | 0.02 | 0.11 | 1 | 1 | 1.0022 | Ser | ||
UAG | 0.02 | 0.11 | 1 | 1 | 1.0022 | Ter | ||
CUG | 0.02 | 0.11 | 1 | 1 | 1.0022 | Leu | ||
CCG | 0.02 | 0.11 | 1 | 1 | 1.0022 | Pro | ||
CAG | 0.02 | 0.11 | 1 | 1 | 1.0022 | Gln | ||
AUG | 0.02 | 0.11 | 1 | 1 | 1.0022 | Met | ||
ACG | 0.02 | 0.11 | 1 | 1 | 1.0022 | Thr | ||
AAG | 0.02 | 0.11 | 1 | 1 | 1.0022 | Lys | ||
4 | Gxx, 9 | GUU | 0.02 | 0.11 | 1 | 1 | 1.0022 | Val |
GCU | 0.02 | 0.11 | 1 | 1 | 1.0022 | Ala | ||
GAU | 0.02 | 0.11 | 1 | 1 | 1.0022 | Asp | ||
GUC | 0.02 | 0.11 | 1 | 1 | 1.0022 | Val | ||
GCC | 0.02 | 0.11 | 1 | 1 | 1.0022 | Ala | ||
GAC | 0.02 | 0.11 | 1 | 1 | 1.0022 | Asp | ||
GUA | 0.02 | 0.11 | 1 | 1 | 1.0022 | Val | ||
GCA | 0.02 | 0.11 | 1 | 1 | 1.0022 | Ala | ||
GAA | 0.02 | 0.11 | 1 | 1 | 1.0022 | Glu | ||
5 | xGx, 9 | UGU | 0.02 | 0.11 | 0 | 0 | 0.0000 | Cys |
UGC | 0.02 | 0.11 | 0 | 0 | 0.0000 | Cys | ||
UGA | 0.02 | 0.11 | 0 | 0 | 0.0000 | Ter | ||
CGU | 0.02 | 0.11 | 0 | 0 | 0.0000 | Arg | ||
CGC | 0.02 | 0.11 | 0 | 0 | 0.0000 | Arg | ||
CGA | 0.02 | 0.11 | 0 | 0 | 0.0000 | Arg | ||
AGU | 0.02 | 0.11 | 0 | 0 | 0.0000 | Ser | ||
AGC | 0.02 | 0.11 | 0 | 0 | 0.0000 | Ser | ||
AGA | 0.02 | 0.11 | 0 | 0 | 0.0000 | Arg | ||
5 | xxx, 27 | UUU | 0.02 | 0.04 | 0 | 0 | 0.0000 | Phe |
UCU | 0.02 | 0.04 | 0 | 0 | 0.0000 | Ser | ||
UAU | 0.02 | 0.04 | 0 | 0 | 0.0000 | Tyr | ||
UUC | 0.02 | 0.04 | 0 | 0 | 0.0000 | Phe | ||
UCC | 0.02 | 0.04 | 0 | 0 | 0.0000 | Ser | ||
UUA | 0.02 | 0.04 | 0 | 0 | 0.0000 | Leu | ||
UCA | 0.02 | 0.04 | 0 | 0 | 0.0000 | Ser | ||
UAA | 0.02 | 0.04 | 0 | 0 | 0.0000 | Ter | ||
CUU | 0.02 | 0.04 | 0 | 0 | 0.0000 | Leu | ||
CCU | 0.02 | 0.04 | 0 | 0 | 0.0000 | Pro | ||
CAU | 0.02 | 0.04 | 0 | 0 | 0.0000 | His | ||
CUC | 0.02 | 0.04 | 0 | 0 | 0.0000 | Leu | ||
CCC | 0.02 | 0.04 | 0 | 0 | 0.0000 | Pro | ||
CAC | 0.02 | 0.04 | 0 | 0 | 0.0000 | His | ||
CUA | 0.02 | 0.04 | 0 | 0 | 0.0000 | Leu | ||
CCA | 0.02 | 0.04 | 0 | 0 | 0.0000 | Pro | ||
CAA | 0.02 | 0.04 | 0 | 0 | 0.0000 | Gln | ||
AUU | 0.02 | 0.04 | 0 | 0 | 0.0000 | Ile | ||
ACU | 0.02 | 0.04 | 0 | 0 | 0.0000 | Thr | ||
AAU | 0.02 | 0.04 | 0 | 0 | 0.0000 | Asn | ||
AUC | 0.02 | 0.04 | 0 | 0 | 0.0000 | Ile | ||
ACC | 0.02 | 0.04 | 0 | 0 | 0.0000 | Thr | ||
AAC | 0.02 | 0.04 | 0 | 0 | 0.0000 | Asn | ||
AUA | 0.02 | 0.04 | 0 | 0 | 0.0000 | Ile | ||
ACA | 0.02 | 0.04 | 0 | 0 | 0.0000 | Thr | ||
AAA | 0.02 | 0.04 | 0 | 0 | 0.0000 | Lys |
aa | CODamino | gCOD+amino | β | |
Group 1 (n=7) | Ala | 4 | 4 | 1.00 |
Val | 4 | 4 | 1.00 | |
Asp | 2 | 2 | 1.00 | |
Glu | 2 | 2 | 1.00 | |
Trp | 1 | 1 | 1.00 | |
Met | 1 | 1 | 1.00 | |
Gly | 4 | 4 | 1.00 | |
Group 2 (n=7) | Leu | 6 | 2 | 0.3333 |
Gln | 2 | 1 | 0.5 | |
Arg | 6 | 2 | 0.3333 | |
Lys | 2 | 1 | 0.5 | |
Ser | 6 | 1 | 0.1667 | |
Thr | 4 | 1 | 0.25 | |
Pro | 4 | 1 | 0.25 | |
Group 3 (n=6) | Cys | 2 | 0 | 0.00 |
Asn | 2 | 0 | 0.00 | |
Ile | 3 | 0 | 0.00 | |
His | 2 | 0 | 0.00 | |
Phe | 2 | 0 | 0.00 | |
Tyr | 2 | 0 | 0.00 |
GENE | NAME | G4 (nt) | Ex | STOP(n=11) | VALID(n=59) | |PTG4| |
BACE1 | Beta-secretase 1 | 33 | 3 | n=0 | n=6 | n=2 |
BCL2 | B-cell lymphoma 2 | 33 | 2 | n=1 | n=6 | n=1 |
23 | n=0 | n=6 | ||||
28 | n=0 | n=6 | ||||
29 | n=1 | n=5 | ||||
34 | n=1 | n=5 | ||||
33 | 3 | n=2 | n=5 | |||
ESR1 | Estrogen receptor alpha (ERα) | 36 | 4 | n=1 | n=5 | n=2 |
KCNH2 | Potassium Voltage-Gated Channel sub family H | 18 | 12 | n=0 | n=0 | NA |
ZNF669 | Member 2 Zinc Finger Protein 669 |
1 | ||||
PRNP | Prion protein | 14 | 2 | n=0 | n=0 | n=1 |
15 | n=0 | n=0 | ||||
20 | n=0 | n=0 | ||||
24 | n=0 | n=6 | ||||
85 | n=6 | n=3 | ||||
TERF2 | Telomeric repeat-binding factor 2 | 55 | 1 | n=0 | n=6 | n=1 |
Disordered regions (IDRs; n=1445;0.00≤p-value < 0.05) | |||||||||||||||||
SL-PT- | SL-PT+ | SL+PT- | SL+PT+ | R1T | R2T | C1T | C2T | A (%) | P (%) | R (%) | |||||||
SLiMS1 | 1078 | 64 | 58 | 9 | 1142 | 67 | 1136 | 73 | 89.90 | 12.32 | 13.43 | ||||||
SLiMS2 | 749 | 18 | 121 | 29 | 767 | 150 | 870 | 47 | 84.84 | 61.70 | 19.33 | ||||||
SLiMS3 | 1212 | 108 | 34 | 9 | 1320 | 43 | 1246 | 117 | 89.58 | 7.69 | 20.93 | ||||||
Proteins with disordered segments (IDPs; n=800;0.00≤p-value < 0.05) | |||||||||||||||||
SL-PT- | SL-PT+ | SL+PT- | SL+PT+ | R1T | R2T | C1T | C2T | A (%) | P (%) | R (%) | |||||||
SLiMS1 | 86 | 12 | 28 | 18 | 98 | 46 | 114 | 30 | 72.22 | 60.00 | 39.10 | ||||||
SLiMS2 | 1 | 1 | 96 | 66 | 2 | 162 | 97 | 67 | 40.85 | 98.50 | 40.74 | ||||||
SLiMS3 | 250 | 57 | 26 | 25 | 307 | 51 | 276 | 82 | 76.81 | 30.48 | 49.01 |
SLiMS | Sample | (zn)n≥2∈pTG4ij∩SLiMSw (p-value) | |
SLiMS1=[ST]PzR | Pz | PG (n=1) | PG (Degenerate) |
SLiMS2=[ED]zzD[AGS] | z[DE] | G[DE] (n=2) | [WGRVAELMKQSTP][AG][DE]z2EG[VADE](p-value=0.00069) |
[DE]z | [DE]G (n=1) | ||
zz[DE] | [LMKQSTP]G[DE] (n=14) [VAE]G[DE] (n=6) [WGR][AG][DE] (n=6) |
||
[DE]zz[DE] | [VAE]G[DE]zzEG[VADE] (n=28) [WGR][AG][DE]zzEG[VADE] (n=24) [WGR][AG][DE]zzEG (n=6) GEzzEG[VADE] (n=4) GEzzEG (n=1) |
||
SLiMS3=[KR]zPzzP | Pzz | PG[VADE] (n=4) | PGV (Degenerate) PGA (Degenerate) PGD (Degenerate) PGE (Degenerate) |
Cellular function | Disordered regions of proteins | |
1. | Signal transduction | DP00274, DP00224, DP00141, DP00332, DP01063, DP00506, DP00418, DP00341, DP00435, DP00613, DP00463, DP00954, DP00959, DP01104, DP00611, DP00519, DP00086, DP00707, DP00712 |
2. | Endocytosis | DP01073, DP01065, DP01066, DP00225 |
3. | Calcium-calmodulin | DP00092, DP00132, DP00561, DP00118, DP00253 |
4. | Myofibril assembly | DP01090 |
5. | Cytoskeleton | DP01056 DP00240, DP01022, DP00169, DP00716, DP00717, DP01100, DP00122 |
6. | Nuclear pore | DP01075, DP01077, DP01079 |
7. | Phototransduction | DP00768, DP00347 |
8. | Targeting | DP00893, DP00609, DP00610, DP01058 |
9. | Transcription | DP00062, DP00177, DP00633, DP00348, DP00786, DP00049, DP00231, DP00873, DP00720, DP00217, DP00081 |
10. | Translation | DP00082, DP00164, DP00229 DP00949, DP00134 |
11. | Synaptic transmission | DP00943 |
12. | Supercoiling | DP00076 |
13. | Binding | DP00539, DP00854, DP01052, DP00659, DP00656 |
14. | Peptide bond formation | DP00944 |
15. | Enzymes | DP00557, DP00032, DP00095, DP00337, DP00379, DP00787, DP00427, DP00429 |
16. | Bacterial/parasitic virulence | |
Secreted toxins | DP00345, DP00591 | |
Cytoadherence | DP00025, DP00065, DP01096 | |
17. | Viral infectivity | |
Cyclophilin interaction | DP00615, DP01031 | |
Chaperones | DP00699, DP00700, DP00674 | |
Capsid assembly | DP00133, DP00876 | |
Membrane fusion | DP01043 | |
Latency | DP01060 | |
18. | Unknown | DP00119 |
Rank | Codon set, Cardinality | Codon | γ | θ | δ | Ω | α=γ.θ.δ+Ω | aa |
1 | GGG, 1 | GGG | 0.02 | 1.00 | 6 | 2 | 2.1200 | Gly |
2 | GxG, 3 | GUG | 0.02 | 0.33 | 1 | 2 | 2.0066 | Val |
GCG | 0.02 | 0.33 | 1 | 2 | 2.0066 | Ala | ||
GAG | 0.02 | 0.33 | 1 | 2 | 2.0066 | Glu | ||
3 | xGG, 3 | UGG | 0.02 | 0.33 | 2 | 1 | 1.0132 | Trp |
CGG | 0.02 | 0.33 | 2 | 1 | 1.0132 | Arg | ||
AGG | 0.02 | 0.33 | 2 | 1 | 1.0132 | Arg | ||
3 | GGx, 3 | GGU | 0.02 | 0.33 | 2 | 1 | 1.0132 | Gly |
GGC | 0.02 | 0.33 | 2 | 1 | 1.0132 | Gly | ||
GGA | 0.02 | 0.33 | 2 | 1 | 1.0132 | Gly | ||
4 | xxG, 9 | UUG | 0.02 | 0.11 | 1 | 1 | 1.0022 | Leu |
UCG | 0.02 | 0.11 | 1 | 1 | 1.0022 | Ser | ||
UAG | 0.02 | 0.11 | 1 | 1 | 1.0022 | Ter | ||
CUG | 0.02 | 0.11 | 1 | 1 | 1.0022 | Leu | ||
CCG | 0.02 | 0.11 | 1 | 1 | 1.0022 | Pro | ||
CAG | 0.02 | 0.11 | 1 | 1 | 1.0022 | Gln | ||
AUG | 0.02 | 0.11 | 1 | 1 | 1.0022 | Met | ||
ACG | 0.02 | 0.11 | 1 | 1 | 1.0022 | Thr | ||
AAG | 0.02 | 0.11 | 1 | 1 | 1.0022 | Lys | ||
4 | Gxx, 9 | GUU | 0.02 | 0.11 | 1 | 1 | 1.0022 | Val |
GCU | 0.02 | 0.11 | 1 | 1 | 1.0022 | Ala | ||
GAU | 0.02 | 0.11 | 1 | 1 | 1.0022 | Asp | ||
GUC | 0.02 | 0.11 | 1 | 1 | 1.0022 | Val | ||
GCC | 0.02 | 0.11 | 1 | 1 | 1.0022 | Ala | ||
GAC | 0.02 | 0.11 | 1 | 1 | 1.0022 | Asp | ||
GUA | 0.02 | 0.11 | 1 | 1 | 1.0022 | Val | ||
GCA | 0.02 | 0.11 | 1 | 1 | 1.0022 | Ala | ||
GAA | 0.02 | 0.11 | 1 | 1 | 1.0022 | Glu | ||
5 | xGx, 9 | UGU | 0.02 | 0.11 | 0 | 0 | 0.0000 | Cys |
UGC | 0.02 | 0.11 | 0 | 0 | 0.0000 | Cys | ||
UGA | 0.02 | 0.11 | 0 | 0 | 0.0000 | Ter | ||
CGU | 0.02 | 0.11 | 0 | 0 | 0.0000 | Arg | ||
CGC | 0.02 | 0.11 | 0 | 0 | 0.0000 | Arg | ||
CGA | 0.02 | 0.11 | 0 | 0 | 0.0000 | Arg | ||
AGU | 0.02 | 0.11 | 0 | 0 | 0.0000 | Ser | ||
AGC | 0.02 | 0.11 | 0 | 0 | 0.0000 | Ser | ||
AGA | 0.02 | 0.11 | 0 | 0 | 0.0000 | Arg | ||
5 | xxx, 27 | UUU | 0.02 | 0.04 | 0 | 0 | 0.0000 | Phe |
UCU | 0.02 | 0.04 | 0 | 0 | 0.0000 | Ser | ||
UAU | 0.02 | 0.04 | 0 | 0 | 0.0000 | Tyr | ||
UUC | 0.02 | 0.04 | 0 | 0 | 0.0000 | Phe | ||
UCC | 0.02 | 0.04 | 0 | 0 | 0.0000 | Ser | ||
UUA | 0.02 | 0.04 | 0 | 0 | 0.0000 | Leu | ||
UCA | 0.02 | 0.04 | 0 | 0 | 0.0000 | Ser | ||
UAA | 0.02 | 0.04 | 0 | 0 | 0.0000 | Ter | ||
CUU | 0.02 | 0.04 | 0 | 0 | 0.0000 | Leu | ||
CCU | 0.02 | 0.04 | 0 | 0 | 0.0000 | Pro | ||
CAU | 0.02 | 0.04 | 0 | 0 | 0.0000 | His | ||
CUC | 0.02 | 0.04 | 0 | 0 | 0.0000 | Leu | ||
CCC | 0.02 | 0.04 | 0 | 0 | 0.0000 | Pro | ||
CAC | 0.02 | 0.04 | 0 | 0 | 0.0000 | His | ||
CUA | 0.02 | 0.04 | 0 | 0 | 0.0000 | Leu | ||
CCA | 0.02 | 0.04 | 0 | 0 | 0.0000 | Pro | ||
CAA | 0.02 | 0.04 | 0 | 0 | 0.0000 | Gln | ||
AUU | 0.02 | 0.04 | 0 | 0 | 0.0000 | Ile | ||
ACU | 0.02 | 0.04 | 0 | 0 | 0.0000 | Thr | ||
AAU | 0.02 | 0.04 | 0 | 0 | 0.0000 | Asn | ||
AUC | 0.02 | 0.04 | 0 | 0 | 0.0000 | Ile | ||
ACC | 0.02 | 0.04 | 0 | 0 | 0.0000 | Thr | ||
AAC | 0.02 | 0.04 | 0 | 0 | 0.0000 | Asn | ||
AUA | 0.02 | 0.04 | 0 | 0 | 0.0000 | Ile | ||
ACA | 0.02 | 0.04 | 0 | 0 | 0.0000 | Thr | ||
AAA | 0.02 | 0.04 | 0 | 0 | 0.0000 | Lys |
aa | CODamino | gCOD+amino | β | |
Group 1 (n=7) | Ala | 4 | 4 | 1.00 |
Val | 4 | 4 | 1.00 | |
Asp | 2 | 2 | 1.00 | |
Glu | 2 | 2 | 1.00 | |
Trp | 1 | 1 | 1.00 | |
Met | 1 | 1 | 1.00 | |
Gly | 4 | 4 | 1.00 | |
Group 2 (n=7) | Leu | 6 | 2 | 0.3333 |
Gln | 2 | 1 | 0.5 | |
Arg | 6 | 2 | 0.3333 | |
Lys | 2 | 1 | 0.5 | |
Ser | 6 | 1 | 0.1667 | |
Thr | 4 | 1 | 0.25 | |
Pro | 4 | 1 | 0.25 | |
Group 3 (n=6) | Cys | 2 | 0 | 0.00 |
Asn | 2 | 0 | 0.00 | |
Ile | 3 | 0 | 0.00 | |
His | 2 | 0 | 0.00 | |
Phe | 2 | 0 | 0.00 | |
Tyr | 2 | 0 | 0.00 |
GENE | NAME | G4 (nt) | Ex | STOP(n=11) | VALID(n=59) | |PTG4| |
BACE1 | Beta-secretase 1 | 33 | 3 | n=0 | n=6 | n=2 |
BCL2 | B-cell lymphoma 2 | 33 | 2 | n=1 | n=6 | n=1 |
23 | n=0 | n=6 | ||||
28 | n=0 | n=6 | ||||
29 | n=1 | n=5 | ||||
34 | n=1 | n=5 | ||||
33 | 3 | n=2 | n=5 | |||
ESR1 | Estrogen receptor alpha (ERα) | 36 | 4 | n=1 | n=5 | n=2 |
KCNH2 | Potassium Voltage-Gated Channel sub family H | 18 | 12 | n=0 | n=0 | NA |
ZNF669 | Member 2 Zinc Finger Protein 669 |
1 | ||||
PRNP | Prion protein | 14 | 2 | n=0 | n=0 | n=1 |
15 | n=0 | n=0 | ||||
20 | n=0 | n=0 | ||||
24 | n=0 | n=6 | ||||
85 | n=6 | n=3 | ||||
TERF2 | Telomeric repeat-binding factor 2 | 55 | 1 | n=0 | n=6 | n=1 |
Disordered regions (IDRs; n=1445;0.00≤p-value < 0.05) | |||||||||||||||||
SL-PT- | SL-PT+ | SL+PT- | SL+PT+ | R1T | R2T | C1T | C2T | A (%) | P (%) | R (%) | |||||||
SLiMS1 | 1078 | 64 | 58 | 9 | 1142 | 67 | 1136 | 73 | 89.90 | 12.32 | 13.43 | ||||||
SLiMS2 | 749 | 18 | 121 | 29 | 767 | 150 | 870 | 47 | 84.84 | 61.70 | 19.33 | ||||||
SLiMS3 | 1212 | 108 | 34 | 9 | 1320 | 43 | 1246 | 117 | 89.58 | 7.69 | 20.93 | ||||||
Proteins with disordered segments (IDPs; n=800;0.00≤p-value < 0.05) | |||||||||||||||||
SL-PT- | SL-PT+ | SL+PT- | SL+PT+ | R1T | R2T | C1T | C2T | A (%) | P (%) | R (%) | |||||||
SLiMS1 | 86 | 12 | 28 | 18 | 98 | 46 | 114 | 30 | 72.22 | 60.00 | 39.10 | ||||||
SLiMS2 | 1 | 1 | 96 | 66 | 2 | 162 | 97 | 67 | 40.85 | 98.50 | 40.74 | ||||||
SLiMS3 | 250 | 57 | 26 | 25 | 307 | 51 | 276 | 82 | 76.81 | 30.48 | 49.01 |
SLiMS | Sample | (zn)n≥2∈pTG4ij∩SLiMSw (p-value) | |
SLiMS1=[ST]PzR | Pz | PG (n=1) | PG (Degenerate) |
SLiMS2=[ED]zzD[AGS] | z[DE] | G[DE] (n=2) | [WGRVAELMKQSTP][AG][DE]z2EG[VADE](p-value=0.00069) |
[DE]z | [DE]G (n=1) | ||
zz[DE] | [LMKQSTP]G[DE] (n=14) [VAE]G[DE] (n=6) [WGR][AG][DE] (n=6) |
||
[DE]zz[DE] | [VAE]G[DE]zzEG[VADE] (n=28) [WGR][AG][DE]zzEG[VADE] (n=24) [WGR][AG][DE]zzEG (n=6) GEzzEG[VADE] (n=4) GEzzEG (n=1) |
||
SLiMS3=[KR]zPzzP | Pzz | PG[VADE] (n=4) | PGV (Degenerate) PGA (Degenerate) PGD (Degenerate) PGE (Degenerate) |
Cellular function | Disordered regions of proteins | |
1. | Signal transduction | DP00274, DP00224, DP00141, DP00332, DP01063, DP00506, DP00418, DP00341, DP00435, DP00613, DP00463, DP00954, DP00959, DP01104, DP00611, DP00519, DP00086, DP00707, DP00712 |
2. | Endocytosis | DP01073, DP01065, DP01066, DP00225 |
3. | Calcium-calmodulin | DP00092, DP00132, DP00561, DP00118, DP00253 |
4. | Myofibril assembly | DP01090 |
5. | Cytoskeleton | DP01056 DP00240, DP01022, DP00169, DP00716, DP00717, DP01100, DP00122 |
6. | Nuclear pore | DP01075, DP01077, DP01079 |
7. | Phototransduction | DP00768, DP00347 |
8. | Targeting | DP00893, DP00609, DP00610, DP01058 |
9. | Transcription | DP00062, DP00177, DP00633, DP00348, DP00786, DP00049, DP00231, DP00873, DP00720, DP00217, DP00081 |
10. | Translation | DP00082, DP00164, DP00229 DP00949, DP00134 |
11. | Synaptic transmission | DP00943 |
12. | Supercoiling | DP00076 |
13. | Binding | DP00539, DP00854, DP01052, DP00659, DP00656 |
14. | Peptide bond formation | DP00944 |
15. | Enzymes | DP00557, DP00032, DP00095, DP00337, DP00379, DP00787, DP00427, DP00429 |
16. | Bacterial/parasitic virulence | |
Secreted toxins | DP00345, DP00591 | |
Cytoadherence | DP00025, DP00065, DP01096 | |
17. | Viral infectivity | |
Cyclophilin interaction | DP00615, DP01031 | |
Chaperones | DP00699, DP00700, DP00674 | |
Capsid assembly | DP00133, DP00876 | |
Membrane fusion | DP01043 | |
Latency | DP01060 | |
18. | Unknown | DP00119 |