Wiggle-predicting functionally flexible regions from primary sequence.
Ontology highlight
ABSTRACT:
SUBMITTER: Gu J
PROVIDER: S-EPMC1500818 | biostudies-literature |
REPOSITORIES: biostudies-literature
Similar Datasets
Project description:Flexible loop regions play a critical role in the biological function of many proteins and have been shown to be involved in ligand binding. In the context of structure-based drug design, using or predicting an incorrect loop configuration can be detrimental to the study if the loop is capable of interacting with the ligand. Three protein systems, each with at least one flexible loop region in close proximity to the known binding site, were selected for loop prediction using the CorLps program; a six residue loop region from phosphoribosylglycinamide formyltransferase (GART), two nine residue loop regions from cytochrome P450 (CYP) 119, and an 11 residue loop region from enolase were selected for loop prediction. The results of this study indicate that the statistically based DFIRE scoring function implemented in the CorLps program did not accurately rank native-like predicted loop configurations in any protein system. In an attempt to improve the ranking of the native-like predicted loop configurations, the MM/GBSA and the optimized MM/GBSA-dsr scoring functions were used to re-rank the predicted loops with and without bound ligand. In general, single snapshot MM/GBSA scoring provided the best ranking of native-like loop configurations. Based on the scoring function analyses presented, the optimal ranking of native-like loop configurations is still a difficult challenge and the choice of the "best" scoring function appears to be system dependent.
| S-EPMC3243024 | biostudies-literature
Project description:Nucleosomes are the fundamental repeating unit of chromatin and comprise the structural building blocks of the living eukaryotic genome. Micrococcal nuclease (MNase) has long been used to delineate nucleosomal organization. Microarray-based nucleosome mapping experiments in yeast chromatin have revealed regularly-spaced translational phasing of nucleosomes. These data have been used to train computational models of sequence-directed nuclesosome positioning, which have identified ubiquitous strong intrinsic nucleosome positioning signals. Here, we successfully apply this approach to nucleosome positioning experiments from human chromatin. The predictions made by the human-trained and yeast-trained models are strongly correlated, suggesting a shared mechanism for sequence-based determination of nucleosome occupancy. In addition, we observed striking complementarity between classifiers trained on experimental data from weakly versus heavily digested MNase samples. In the former case, the resulting model accurately identifies nucleosome-forming sequences; in the latter, the classifier excels at identifying nucleosome-free regions. Using this model we are able to identify several characteristics of nucleosome-forming and nucleosome-disfavoring sequences. First, by combining results from each classifier applied de novo across the human ENCODE regions, the classifier reveals distinct sequence composition and periodicity features of nucleosome-forming and nucleosome-disfavoring sequences. Short runs of dinucleotide repeat appear as a hallmark of nucleosome-disfavoring sequences, while nucleosome-forming sequences contain short periodic runs of GC base pairs. Second, we show that nucleosome phasing is most frequently predicted flanking nucleosome-free regions. The results suggest that the major mechanism of nucleosome positioning in vivo is boundary-event-driven and affirm the classical statistical positioning theory of nucleosome organization.
| S-EPMC2515632 | biostudies-literature
Project description:Secondary structure elements in the coding regions of mRNAs play an important role in gene expression and regulation, but distinguishing functional from non-functional structures remains challenging. Here we investigate the dependence of sequence-structure relationships in the coding regions on temperature based on the recent PARTE data by Wan et al. Our main finding is that the regions with high and low thermostability (high Tm and low Tm regions) are under evolutionary pressure to preserve RNA secondary structure and primary sequence, respectively. Sequences of low Tm regions display a higher degree of evolutionary conservation compared to high Tm regions. Low Tm regions are under strong synonymous constraint, while high Tm regions are not. These findings imply that high Tm regions contain thermo-stable functionally important RNA structures, which impose relaxed evolutionary constraint on sequence as long as the base-pairing patterns remain intact. By contrast, low thermostability regions contain single-stranded functionally important conserved RNA sequence elements accessible for binding by other molecules. We also find that theoretically predicted structures of paralogous mRNA pairs become more similar with growing temperature, while experimentally measured structures tend to diverge, which implies that the melting pathways of RNA structures cannot be fully captured by current computational approaches.
| S-EPMC5449622 | biostudies-literature
Project description:With a large number of disordered proteins and their important functions discovered, it is highly desired to develop effective methods to computationally predict protein disordered regions. In this study, based on Random Forest (RF), Maximum Relevancy Minimum Redundancy (mRMR), and Incremental Feature Selection (IFS), we developed a new method to predict disordered regions in proteins. The mRMR criterion was used to rank the importance of all candidate features. Finally, top 128 features were selected from the ranked feature list to build the optimal model, including 92 Position Specific Scoring Matrix (PSSM) conservation score features and 36 secondary structure features. As a result, Matthews correlation coefficient (MCC) of 0.3895 was achieved on the training set by 10-fold cross-validation. On the basis of predicting results for each query sequence by using the method, we used the scanning and modification strategy to improve the performance. The accuracy (ACC) and MCC were increased by 4% and almost 0.2%, respectively, compared with other three popular predictors: DISOPRED, DISOclust, and OnD-CRF. The selected features may shed some light on the understanding of the formation mechanism of disordered structures, providing guidelines for experimental validation.
| S-EPMC3654632 | biostudies-literature
Project description:At the earliest break of ancient hominins from their primate relatives in vocal communication, we propose a selection pressure on vocal fitness signaling by hominin infants. Exploratory vocalizations, not tied to expression of distress or immediate need, could have helped persuade parents of the wellness and viability of the infants who produced them. We hypothesize that hominin parents invested more in infants who produced such signals of fitness plentifully, neglecting or abandoning them less often than infants who produced the sounds less frequently. Selection for such exploratory vocalization provided a critically important inclination and capability relevant to language, we reason, because the system that encouraged spontaneous vocalization also made vocalization functionally flexible to an extent that has not been observed in any other animal. Although this vocal flexibility did not by itself create language, it provided an essential foundation upon which language would evolve through a variety of additional steps. In evaluating this speculation, we consider presumable barriers to evolving language that are thought to be implications of Darwinian Theory. It has been claimed that communication always involves sender self-interest and that self-interest leads to deceit, which is countered through clever detection by receivers. The constant battle of senders and receivers has been thought to pose an insuperable challenge to honest communication, which has been viewed as a requirement of language. To make communication honest, it has been proposed that stable signaling requires costly handicaps for the sender, and since language cannot entail high cost, the reasoning has suggested an insurmountable obstacle to the evolution of language. We think this presumed honesty barrier is an illusion that can be revealed by recognition of the fact that language is not inherently honest and in light of the distinction between illocutionary force and semantics. Our paper also considers barriers to the evolution of language (not having to do with honesty) that we think may have actually played important roles in preventing species other than humans from evolving language.
| S-EPMC7870697 | biostudies-literature
Project description:DNA methylation is an important mechanism by which gene transcription and hence cellular function are regulated. Here, we provide detailed functional genome-wide methylome maps of 5 primary peripheral blood leukocyte subsets including T cells, B cells, monocytes/macrophages, and neutrophils obtained from healthy individuals. A comparison of these methylomes revealed highly specific cell-lineage and cell-subset methylation profiles. DNA hypomethylation is known to be permissive for gene expression and we identified cell-subset-specific hypomethylated regions (HMRs) that strongly correlate with gene transcription levels suggesting these HMRs may regulate corresponding cell functions. Single-nucleotide polymorphisms associated with immune-mediated disease in genome-wide association studies preferentially localized to these cell-specific regulatory HMRs, offering insight into pathogenesis by highlighting cell subsets in which specific epigenetic changes may drive disease. Our data provide a valuable reference tool for researchers aiming to investigate the role of DNA methylation in regulating primary leukocyte function in health and immune-mediated disease.
| S-EPMC3862273 | biostudies-literature
Project description:BackgroundRNA-binding proteins (RBPs) play diverse roles in eukaryotic RNA processing. Despite their pervasive functions in coding and noncoding RNA biogenesis and regulation, elucidating the sequence specificities that define protein-RNA interactions remains a major challenge. Recently, CLIP-seq (Cross-linking immunoprecipitation followed by high-throughput sequencing) has been successfully implemented to study the transcriptome-wide binding patterns of SRSF1, PTBP1, NOVA and fox2 proteins. These studies either adopted traditional methods like Multiple EM for Motif Elicitation (MEME) to discover the sequence consensus of RBP's binding sites or used Z-score statistics to search for the overrepresented nucleotides of a certain size. We argue that most of these methods are not well-suited for RNA motif identification, as they are unable to incorporate the RNA structural context of protein-RNA interactions, which may affect to binding specificity. Here, we describe a novel model-based approach--RNAMotifModeler to identify the consensus of protein-RNA binding regions by integrating sequence features and RNA secondary structures.ResultsAs an example, we implemented RNAMotifModeler on SRSF1 (SF2/ASF) CLIP-seq data. The sequence-structural consensus we identified is a purine-rich octamer 'AGAAGAAG' in a highly single-stranded RNA context. The unpaired probabilities, the probabilities of not forming pairs, are significantly higher than negative controls and the flanking sequence surrounding the binding site, indicating that SRSF1 proteins tend to bind on single-stranded RNA. Further statistical evaluations revealed that the second and fifth bases of SRSF1octamer motif have much stronger sequence specificities, but weaker single-strandedness, while the third, fourth, sixth and seventh bases are far more likely to be single-stranded, but have more degenerate sequence specificities. Therefore, we hypothesize that nucleotide specificity and secondary structure play complementary roles during binding site recognition by SRSF1.ConclusionIn this study, we presented a computational model to predict the sequence consensus and optimal RNA secondary structure for protein-RNA binding regions. The successful implementation on SRSF1 CLIP-seq data demonstrates great potential to improve our understanding on the binding specificity of RNA binding proteins.
| S-EPMC3287504 | biostudies-literature
Project description:We analyzed functionality and relative distribution of genetic variants across the complete Oryza sativa genome, using the 40 million single nucleotide polymorphisms (SNPs) dataset from the 3,000 Rice Genomes Project (http://snp-seek.irri.org), the largest and highest density SNP collection for any higher plant. We have shown that the DNA-binding transcription factors (TFs) are the most conserved group of genes, whereas kinases and membrane-localized transporters are the most variable ones. TFs may be conserved because they belong to some of the most connected regulatory hubs that modulate transcription of vast downstream gene networks, whereas signaling kinases and transporters need to adapt rapidly to changing environmental conditions. In general, the observed profound patterns of nucleotide variability reveal functionally important genomic regions. As expected, nucleotide diversity is much higher in intergenic regions than within gene bodies (regions spanning gene models), and protein-coding sequences are more conserved than untranslated gene regions. We have observed a sharp decline in nucleotide diversity that begins at about 250 nucleotides upstream of the transcription start and reaches minimal diversity exactly at the transcription start. We found the transcription termination sites to have remarkably symmetrical patterns of SNP density, implying presence of functional sites near transcription termination. Also, nucleotide diversity was significantly lower near 3' UTRs, the area rich with regulatory regions.
| S-EPMC5075931 | biostudies-literature
Project description:Adult brains are functionally flexible, a unique characteristic that is thought to contribute to cognitive flexibility. While tools to assess cognitive flexibility during early infancy are lacking, we aimed to assess the spatiotemporal developmental features of "neural flexibility" during the first 2 y of life. Fifty-two typically developing children 0 to 2 y old were longitudinally imaged up to seven times during natural sleep using resting-state functional MRI. Using a sliding window approach, MR-derived neural flexibility, a quantitative measure of the frequency at which brain regions change their allegiance from one functional module to another during a given time period, was used to evaluate the temporal emergence of neural flexibility during early infancy. Results showed that neural flexibility of whole brain, motor, and high-order brain functional networks/regions increased significantly with age, while visual regions exhibited a temporally stable pattern, suggesting spatially and temporally nonuniform developmental features of neural flexibility. Additionally, the neural flexibility of the primary visual network at 3 mo of age was significantly and negatively associated with cognitive ability evaluated at 5/6 y of age. The "flexible club," comprising brain regions with neural flexibility significantly higher than whole-brain neural flexibility, were consistent with brain regions known to govern cognitive flexibility in adults and exhibited unique characteristics when compared to the functional hub and diverse club regions. Thus, MR-derived neural flexibility has the potential to reveal the underlying neural substrates for developing a cognitively flexible brain during early infancy.
| S-EPMC7519318 | biostudies-literature
Project description:Members of the nuclear receptor superfamily differentiate in terms of specificity for DNA recognition and binding, oligomeric state, and ligand binding. The wide range of specificities are impressive given the high degree of sequence conservation in the DNA binding domain (DBD) and moderate sequence conservation with high structural similarity within the ligand binding domains (LBDs). Determining sequence positions that are conserved within nuclear receptor subfamilies can provide important indicators into the structural dynamics that translate to oligomeric state of the active receptor, DNA binding specificity and ligand affinity and selectivity. Here we present a method to analyze sequence data from all nuclear receptors that facilitates detection of co-evolving pairs using Mutual Information (MI). Using this method we demonstrate that MI can reveal functionally important sequence positions within the superfamily and the approach identified three sequence positions that have conserved sequence patterns across all nuclear receptors and subfamilies. Interestingly, two of the sequence positions identified are located within the DBD CII and the third was within Helix c of the DBD. These sequences are located within the heterodimer interface of PPAR? (CII) and RXR? (Helix c) based on PDB:3DZU. Helix c of PPAR?, which is not involved in the DBD dimer interface, binds the minor groove in the 5' flanking region in a consensus PPAR? response element (PPRE) and the corresponding RXR? (CII) is found in the 3' flanking region of RXRE (3DZU). As these three sequence positions represent unique identifiers for all nuclear receptors and they are located within the dimer interface of PPAR?-RXR? DBD (3DZU) interfacing with the flanking regions of the NRRE, we conclude they are critical sequence positions perhaps dictating nuclear receptor (NR) DNA binding specificity.
| S-EPMC3049237 | biostudies-literature