Where does the alignment score distribution shape come from?
Ontology highlight
ABSTRACT:
SUBMITTER: Ortet P
PROVIDER: S-EPMC3023300 | biostudies-literature | 2010 Dec
REPOSITORIES: biostudies-literature
Publications
Where does the alignment score distribution shape come from?
Ortet Philippe P Bastien Olivier O
Evolutionary bioinformatics online 20101212
Alignment algorithms are powerful tools for searching for homologous proteins in databases, providing a score for each sequence present in the database. It has been well known for 20 years that the shape of the score distribution looks like an extreme value distribution. The extremely large number of times biologists face this class of distributions raises the question of the evolutionary origin of this probability law.WE INVESTIGATED THE POSSIBILITY OF DERIVING THE MAIN PROPERTIES OF SEQUENCE A ...[more]
Similar Datasets
Project description:The extracellular matrix (ECM) and integrins collaborate to regulate gene expression associated with cell growth, differentiation and survival. Biochemical and molecular analyses of integrin signalling pathways have uncovered several critical cytoplasmic proteins that link the ECM and integrins to intracellular pathways that may contribute to anchorage-dependent growth. A large body of evidence now indicates that the non-receptor protein kinases focal adhesion kinase (FAK) and specific members of the mitogen-activated protein kinases (MAPKs), including the extracellular-signal-regulated kinases (ERKs), mediate these ECM- and integrin-derived signalling events. However, little is known about how FAK and MAPKs contribute to biological processes other than cell proliferation or migration. In addition, remarkably little is known concerning the signalling events that occur in cells that adhere to complex multivalent extracellular matrices via multiple integrin receptors. Given the stringent requirement for attaining a proper morphology in ECM/integrin-directed cell behaviour, it is still not clear how cell shape and tissue architecture impact upon intracellular signalling programmes involving FAK and MAPKs. However, the recent discovery that members of the Rho family of small GTPases are able to regulate ECM/integrin pathways that modulate both cell shape and intracellular signalling provides new insights into how cell morphology and signal transduction become integrated, especially within three-dimensional differentiated tissues.
| S-EPMC1220180 | biostudies-other
Project description:Cardiac organogenesis and pathogenesis are both characterized by changes in myocyte shape, cytoskeletal architecture, and the extracellular matrix (ECM). However, the mechanisms by which the ECM influences myocyte shape and myofibrillar patterning are unknown. We hypothesized that geometric cues in the ECM align sarcomeres by directing the actin network orientation. To test our hypothesis, we cultured neonatal rat ventricular myocytes on islands of micro-patterned ECM to measure how they remodeled their cytoskeleton in response to extracellular cues. Myocytes spread and assumed the shape of circular and rectangular islands and reorganized their cytoskeletons and myofibrillar arrays with respect to the ECM boundary conditions. Circular myocytes did not assemble predictable actin networks nor organized sarcomere arrays. In contrast, myocytes cultured on rectangular ECM patterns with aspect ratios ranging from 1:1 to 7:1 aligned their sarcomeres in predictable and repeatable patterns based on highly localized focal adhesion complexes. Examination of averaged alpha-actinin images revealed invariant sarcomeric registration irrespective of myocyte aspect ratio. Since the sarcomere sub-units possess a fixed length, this observation indicates that cytoskeleton configuration is length-limited by the extracellular boundary conditions. These results indicate that modification of the extracellular microenvironment induces dynamic reconfiguring of the myocyte shape and intracellular architecture. Furthermore, geometric boundaries such as corners induce localized myofibrillar anisotropy that becomes global as the myocyte aspect ratio increases.
| S-EPMC4492320 | biostudies-literature
Project description:Global change is increasingly challenging the sustainable provisioning of ecosystem services to society. Addressing future uncertainty and risk has therefore become a central problem of ecosystem management. With risk management and resilience-based stewardship, two contrasting approaches have been proposed to address this issue. Whereas one is concentrated on anticipating and mitigating risks, the other is focused on fostering the ability to absorb perturbations and maintain desired properties. While they have hitherto been discussed largely separately in the literature, I here propose a unifying framework of anticipating risks and fostering resilience in ecosystem management. Anticipatory action is advocated when the predictability of risk is high and sufficient knowledge to address it is available. Conversely, in situations in which predictability and knowledge are limited, resilience-based measures are paramount. I conclude that, by adopting a purposeful combination of insights from risk and resilience research, we can make ecosystem services provisioning more robust to future uncertainty and change.
| S-EPMC4340566 | biostudies-literature
Project description:It is well established that the risk of developing multiple sclerosis is substantially increased in the relatives of affected individuals and that most of this increase is genetically determined. The observed pattern of familial recurrence risk has long suggested that multiple variants are involved, but it has proven difficult to identify individual risk variants and little has been established about the genetic architecture underlying susceptibility. By using data from two independent genome-wide association studies (GWAS), we demonstrate that a substantial proportion of the thousands of variants that individually fail to show statistically significant evidence of association have allele frequencies in cases that are skewed away from the null distribution through the effects of multiple as-yet-unidentified risk loci. The collective effect of 12,627 SNPs with Cochran-Mantel-Haenszel test (p < 0.2) in our discovery GWAS set optimally explains approximately 3% of the variance in MS risk in our independent target GWAS set, estimated by Nagelkerke's pseudo-R(2). This model has a highly significant fit (p = 9.90E-19). These results statistically demonstrate a polygenic component to MS susceptibility and suggest that the risk alleles identified to date represent just the tip of an iceberg of risk variants likely to include hundreds of modest effects and possibly thousands of very small effects.
| S-EPMC2850422 | biostudies-literature
Project description:OBJECTIVE:To develop a quantitative method for characterizing gestational sac shape. METHODS:Twenty first-trimester gestational sacs in normal pregnancies were studied with three-dimensional (3D) ultrasonography. The 3D coordinates of surface-point sets were obtained for each sac using 30-, 15- and six-slice sampling. Cubic spline interpolation was used with the 15- and six-slice surface-point samples to generate coordinates for those 30-slice surface points not measured. Interpolated and measured values, the latter from the 30-slice sample, were compared and the percent error calculated. Cubic spline interpolation was used to determine the coordinates of a standard surface-point sample (3660) for each sac in each slice sample. These coordinate data were used to give each sac a standard configuration by moving its center of gravity to the origin, aligning its inertial axes along the coordinate axes and converting its volume to 1.0 mL. In this form, a volume shape descriptor could be generated for each sac that was then transformed into a vector containing only shape information. The 20 shape vectors of each slice sample were subjected to principal components analysis, and principal component scores (PCSs) calculated. The first four PCSs were used to define a gestational sac shape score (GSSS-30, GSSS-15 or GSSS-6) for each sac in a given slice sample. The characteristics of each set of GSSSs were determined and those for the GSSS-15 and GSSS-6 were compared with the GSSS-30 characteristics. RESULTS:Cubic spline interpolations were very accurate in most cases, with means close to 0%, and approximately 95% of the errors being less than 10%. GSSS-30 accounted for 67.6% of the shape variance, had a mean of zero and an SD of 1.1, was normally distributed and was not related to menstrual age (R=-0.16, P=0.51). GSSS-15 and GSSS-6 had essentially the same characteristics. No significant differences between individual GSSS-30 values and those for GSSS-15 or GSSS-6 were found, indicating the absence of a slice sample effect. CONCLUSION:Using sophisticated mathematical methods, the gestational sac shape, initially represented by the 3D coordinates of 3660 surface points, was converted to a single number, the GSSS. This score had the appropriate properties for quantitatively characterizing normal, first-trimester gestational sac shapes. As it can be obtained from as few as six slices, it should be useful in many clinical situations. This novel approach has the potential for providing quantitative shape information about a variety of biological shapes and how they change over time.
| S-EPMC3516405 | biostudies-literature
Project description:Due to dramatic advances in DNA technology, quantitative measures of annotation data can now be obtained in continuous coordinates across the entire genome, allowing various heterogeneous 'genomic landscapes' to emerge. Although much effort has been devoted to comparing DNA sequences, not much attention has been given to comparing these large quantities of data comprehensively. In this article, we introduce a method for rapidly detecting local regions that show high correlations between genomic landscapes. We overcame the size problem for genome-wide data by converting the data into series of symbols and then carrying out sequence alignment. We also decomposed the oscillation of the landscape data into different frequency bands before analysis, since the real genomic landscape is a mixture of embedded and confounded biological processes working at different scales in the cell nucleus. To verify the usefulness and generality of our method, we applied our approach to well investigated landscapes from the human genome, including several histone modifications. Furthermore, by applying our method to over 20 genomic landscapes in human and 12 in mouse, we found that DNA replication timing and the density of Alu insertions are highly correlated genome-wide in both species, even though the Alu elements have amplified independently in the two genomes. To our knowledge, this is the first method to align genomic landscapes at multiple scales according to their shape.
| S-EPMC3413149 | biostudies-literature
Project description:Although membrane shape varies greatly throughout the cell, the contribution of membrane curvature to transmembrane protein targeting is unknown because of the numerous sorting mechanisms that take place concurrently in cells. To isolate the effect of membrane shape, we used cell-sized giant unilamellar vesicles (GUVs) containing either the potassium channel KvAP or the water channel AQP0 to form membrane nanotubes with controlled radii. Whereas the AQP0 concentrations in flat and curved membranes were indistinguishable, KvAP was enriched in the tubes, with greater enrichment in more highly curved membranes. Fluorescence recovery after photobleaching measurements showed that both proteins could freely diffuse through the neck between the tube and GUV, and the effect of each protein on membrane shape and stiffness was characterized using a thermodynamic sorting model. This study establishes the importance of membrane shape for targeting transmembrane proteins and provides a method for determining the effective shape and flexibility of membrane proteins.
| S-EPMC3954780 | biostudies-literature
Project description:MOTIVATION: Based on the next generation genome sequencing technologies, a variety of biological applications are developed, while alignment is the first step once the sequencing reads are obtained. In recent years, many software tools have been developed to efficiently and accurately align short reads to the reference genome. However, there are still many reads that can't be mapped to the reference genome, due to the exceeding of allowable mismatches. Moreover, besides the unmapped reads, the reads with low mapping qualities are also excluded from the downstream analysis, such as variance calling. If we can take advantages of the confident segments of these reads, not only can the alignment rates be improved, but also more information will be provided for the downstream analysis. RESULTS: This paper proposes a method, called RAUR (Re-align the Unmapped Reads), to re-align the reads that can not be mapped by alignment tools. Firstly, it takes advantages of the base quality scores (reported by the sequencer) to figure out the most confident and informative segments of the unmapped reads by controlling the number of possible mismatches in the alignment. Then, combined with an alignment tool, RAUR re-align these segments of the reads. We run RAUR on both simulated data and real data with different read lengths. The results show that many reads which fail to be aligned by the most popular alignment tools (BWA and Bowtie2) can be correctly re-aligned by RAUR, with a similar Precision. Even compared with the BWA-MEM and the local mode of Bowtie2, which perform local alignment for long reads to improve the alignment rate, RAUR also shows advantages on the Alignment rate and Precision in some cases. Therefore, the trimming strategy used in RAUR is useful to improve the Alignment rate of alignment tools for the next-generation genome sequencing. AVAILABILITY: All source code are available at http://netlab.csu.edu.cn/bioinformatics/RAUR.html.
| S-EPMC4402702 | biostudies-literature
Project description:BACKGROUND: The chemical property and biological function of a protein is a direct consequence of its primary structure. Several algorithms have been developed which determine alignment and similarity of primary protein sequences. However, character based similarity cannot provide insight into the structural aspects of a protein. We present a method based on spectral similarity to compare subsequences of amino acids that behave similarly but are not aligned well by considering amino acids as mere characters. This approach finds a similarity score between sequences based on any given attribute, like hydrophobicity of amino acids, on the basis of spectral information after partial conversion to the frequency domain. RESULTS: Distance matrices of various branches of the human kinome, that is the full complement of human kinases, were developed that matched the phylogenetic tree of the human kinome establishing the efficacy of the global alignment of the algorithm. PKCd and PKCe kinases share close biological properties and structural similarities but do not give high scores with character based alignments. Detailed comparison established close similarities between subsequences that do not have any significant character identity. We compared their known 3D structures to establish that the algorithm is able to pick subsequences that are not considered similar by character based matching algorithms but share structural similarities. Similarly many subsequences with low character identity were picked between xyna-theau and xyna-clotm F/10 xylanases. Comparison of 3D structures of the subsequences confirmed the claim of similarity in structure. CONCLUSION: An algorithm is developed which is inspired by successful application of spectral similarity applied to music sequences. The method captures subsequences that do not align by traditional character based alignment tools but give rise to similar secondary and tertiary structures. The Spectral Similarity Score (SSS) is an extension to the conventional similarity methods and results indicate that it holds a strong potential for analysis of various biological sequences and structural variations in proteins.
| S-EPMC1131888 | biostudies-literature