Correlation analysis for protein evolutionary family based on amino acid position mutations and application in PDZ domain

Background

It has been widely recognized that the mutations at specific directions are caused by the functional constraints in protein family and the directional mutations at certain positions control the evolutionary direction of the protein family. The mutations at different positions, even distantly separated, are mutually coupled and form an evolutionary network. Finding the controlling mutative positions and the mutative network among residues are firstly important for protein rational design and enzyme engineering.

Methodology

A computational approach, namely amino acid position conservation-mutation correlation analysis (CMCA), is developed to predict mutually mutative positions and find the evolutionary network in protein family. The amino acid position mutative function, which is the foundational equation of CMCA measuring the mutation of a residue at a position, is derived from the MSA (multiple structure alignment) database of protein evolutionary family. Then the position conservation correlation matrix and position mutation correlation matrix is constructed from the amino acid position mutative equation. Unlike traditional SCA (statistical coupling analysis) approach, which is based on the statistical analysis of position conservations, the CMCA focuses on the correlation analysis of position mutations.

Conclusions

As an example the CMCA approach is used to study the PDZ domain of protein family, and the results well illustrate the distantly allosteric mechanism in PDZ protein family, and find the functional mutative network among residues. We expect that the CMCA approach may find applications in protein engineering study, and suggest new strategy to improve bioactivities and physicochemical properties of enzymes.     

Estimates of statistical significance for comparison of individual positions in multiple sequence alignments

Background  

Profile-based analysis of multiple sequence alignments (MSA) allows for accurate comparison of protein families. Here, we address the problems of detecting statistically confident dissimilarities between (1) MSA position and a set of predicted residue frequencies, and (2) between two MSA positions. These problems are important for (i) evaluation and optimization of methods predicting residue occurrence at protein positions; (ii) detection of potentially misaligned regions in automatically produced alignments and their further refinement; and (iii) detection of sites that determine functional or structural specificity in two related families.     

Plasma Homocysteine,Vitamin B12 and Folate Levels in Multiple System Atrophy: A Case-Control Study

Background

Multiple system atrophy (MSA) is a neurodegenerative disease, and its pathological hallmark is the accumulation of α-synuclein proteins. Homocysteine (Hcy) is an intermediate amino acid generated during the metabolism of methionine. Hcy may contribute to the pathogenesis of neurodegenerative disorders. Vitamin B12 and folate are cofactors necessary for the methylation of homocysteine.

Methods

This study compared the levels of serum Hcy, vitamin B12 and folate in patients with MSA with those in healthy people to reveal the possible association between MSA and plasma levels of Hcy, vitamin B12 and folate. We enrolled 161 patients with MSA and 161 healthy people in this study. The association between MSA and the levels of Hcy, vitamin B12 and folate were analyzed using binary logistic regression.

Results

The mean level of Hcy in patients with MSA was significantly higher than that in healthy controls (16.23 ± 8.09 umol/l vs 14.04 ± 4.25 umol/l, p < 0.05). After adjusting for age, sex and medical history, the odds ratio for Hcy was 1.07 (95% CI = 1.01–1.13, p < 0.05) for patients with MSA. Vitamin B12 and folate levels were not significantly different between patients with MSA and controls.

Conclusion

Our data suggest that higher levels of Hcy may be associated with an increased risk for MSA.     

Abnormal Pulmonary Function and Respiratory Muscle Strength Findings in Chinese Patients with Parkinson’s Disease and Multiple System Atrophy–Comparison with Normal Elderly

Background

There have been limited comparative data regarding the investigations on pulmonary and respiratory muscle function in the patients with different parkinsonism disorders such as Parkinson’s disease (PD) and multiple system atrophy (MSA) versus normal elderly. The present study is aiming to characterize the performance of pulmonary function and respiratory muscle strength in PD and MSA, and to investigate the association with severity of motor symptoms and disease duration.

Methods

Pulmonary function and respiratory muscle strength tests were performed in 30 patients with PD, 27 with MSA as well as in 20 age-, sex-, height-, weight-matched normal elderly controls. All the patients underwent United Parkinson’s disease rating scale (UPDRS) or united multiple system atrophy rating scale (UMSARS) separately as diagnosed.

Results

Vital capacity, forced expiratory volume in 1 second and forced vital capacity decreased, residual volume and ratio of residual volume to total lung capacity increased in both PD and MSA groups compared to controls (p<0.05). Diffusing capacity was decreased in the MSA group, compared with PD and normal elderly control groups (p<0.05). Respiratory muscle strength was lower in both PD and MSA groups than in controls (p<0.05). The values representing spirometry function and respiratory muscle strength were found to have a negative linear correlation with mean score of UPDRS-III in PD and mean score of UMSARS-I in MSA. Respiratory muscle strength showed a negative linear correlation with the mean score of UMSARS-II and disease duration in MSA patients.

Conclusions

These findings suggest that respiratory dysfunction is involved in PD and MSA. Respiratory muscle strength is remarkably reduced, and some of the parameters correlate with disease duration and illness severity. The compromised respiratory function in neurodegenerative disorders should be the focus of further researches.     

An evolutionary relationship between Stearoyl-CoA Desaturase (SCD) protein sequences involved in fatty acid metabolism

Background:

Stearoyl-CoA desaturase (SCD) is a key enzyme that converts saturated fatty acids (SFAs) to monounsaturated fatty acids (MUFAs) in fat biosynthesis. Despite being crucial for interpreting SCDs’ roles across species, the evolutionary relationship of SCD proteins across species has yet to be elucidated. This study aims to present this evolutionary relationship based on amino acid sequences.

Methods:

Using Multiple Sequence Alignment (MSA) and phylogenetic construction methods, a hypothetical evolutionary relationship was generated between the stearoyl-CoA desaturase (SCD) protein sequences between 18 different species.

Results:

SCD protein sequences from Homo sapiens, Pan troglodytes (chimpanzee), and Pongo abelii (orangutan) have the lowest genetic distances of 0.006 of the 18 species studied. Capra hircus (goat) and Ovis aries (Sheep) had the next lowest genetic distance of 0.023. These farm animals are 99.987% identical at the amino acid level.

Conclusions:

The SCD proteins are conserved in these 18 species, and their evolutionary relationships are similar. Key Words: Phylogenetic analysis, Stearoyl-CoA desaturase (SCD) proteins, Multiple sequence alignment     

Wavelets filtering for classification of very noisy electron microscopic single particles images- application on structure determination of VP5-VP19C recombinant

Background

Images of frozen hydrated [vitrified] virus particles were taken close-to-focus in an electron microscope containing structural signals at high spatial frequencies. These images had very low contrast due to the high levels of noise present in the image. The low contrast made particle selection, classification and orientation determination very difficult. The final purpose of the classification is to improve the signal-to-noise ratio of the particle representing the class, which is usually the average. In this paper, the proposed method is based on wavelet filtering and multi-resolution processing for the classification and reconstruction of this very noisy data. A multivariate statistical analysis (MSA) is used for this classification.

Results

The MSA classification method is noise dependant. A set of 2600 projections from a 3D map of a herpes simplex virus -to which noise was added- was classified by MSA. The classification shows the power of wavelet filtering in enhancing the quality of class averages (used in 3D reconstruction) compared to Fourier band pass filtering. A 3D reconstruction of a recombinant virus (VP5-VP19C) is presented as an application of multi-resolution processing for classification and reconstruction.

Conclusion

The wavelet filtering and multi-resolution processing method proposed in this paper offers a new way for processing very noisy images obtained from electron cryo-microscopes. The multi-resolution and filtering improves the speed and accuracy of classification, which is vital for the 3D reconstruction of biological objects. The VP5-VP19C recombinant virus reconstruction presented here is an example, which demonstrates the power of this method. Without this processing, it is not possible to get the correct 3D map of this virus.

     

Reducing the false positive rate in the non-parametric analysis of molecular coevolution

Background  

The strength of selective constraints operating on amino acid sites of proteins has a multifactorial nature. In fact, amino acid sites within proteins coevolve due to their functional and/or structural relationships. Different methods have been developed that attempt to account for the evolutionary dependencies between amino acid sites. Researchers have invested a significant effort to increase the sensitivity of such methods. However, the difficulty in disentangling functional co-dependencies from historical covariation has fuelled the scepticism over their power to detect biologically meaningful results. In addition, the biological parameters connecting linear sequence evolution to structure evolution remain elusive. For these reasons, most of the evolutionary studies aimed at identifying functional dependencies among protein domains have focused on the structural properties of proteins rather than on the information extracted from linear multiple sequence alignments (MSA). Non-parametric methods to detect coevolution have been reported to be especially susceptible to produce false positive results based on the properties of MSAs. However, no formal statistical analysis has been performed to definitively test the differential effects of these properties on the sensitivity of such methods.     

Proinflammatory Cytokines Are Elevated in Serum of Patients with Multiple System Atrophy

Background

Despite several lines of evidence from preclinical and post-mortem studies suggesting that inflammation is involved in Multiple System Atrophy (MSA), no previous studies have measured peripheral indices of inflammation in MSA patients.

Methods

We measured C-reactive protein, interleukin (IL)-6, soluble IL-2 receptor and tumor necrosis factor (TNF)-α in blood samples from MSA patients (n = 14) and healthy controls (n = 40).

Results

IL-6 and TNF-α were significantly elevated in MSA patients compared to healthy controls. After controlling for the potentially confounding effects of age, gender, and somatic co-morbidities, a diagnosis of MSA was still significantly associated with high levels of TNF-α. Higher TNF-α levels were associated with less severe motor symptoms and earlier disease stage.

Conclusions

Our findings are in line with the hypothesis that inflammation might be involved at an early stage of MSA pathophysiology.     

Accurate Simulation and Detection of Coevolution Signals in Multiple Sequence Alignments

Background

While the conserved positions of a multiple sequence alignment (MSA) are clearly of interest, non-conserved positions can also be important because, for example, destabilizing effects at one position can be compensated by stabilizing effects at another position. Different methods have been developed to recognize the evolutionary relationship between amino acid sites, and to disentangle functional/structural dependencies from historical/phylogenetic ones.

Methodology/Principal Findings

We have used two complementary approaches to test the efficacy of these methods. In the first approach, we have used a new program, MSAvolve, for the in silico evolution of MSAs, which records a detailed history of all covarying positions, and builds a global coevolution matrix as the accumulated sum of individual matrices for the positions forced to co-vary, the recombinant coevolution, and the stochastic coevolution. We have simulated over 1600 MSAs for 8 protein families, which reflect sequences of different sizes and proteins with widely different functions. The calculated coevolution matrices were compared with the coevolution matrices obtained for the same evolved MSAs with different coevolution detection methods. In a second approach we have evaluated the capacity of the different methods to predict close contacts in the representative X-ray structures of an additional 150 protein families using only experimental MSAs.

Conclusions/Significance

Methods based on the identification of global correlations between pairs were found to be generally superior to methods based only on local correlations in their capacity to identify coevolving residues using either simulated or experimental MSAs. However, the significant variability in the performance of different methods with different proteins suggests that the simulation of MSAs that replicate the statistical properties of the experimental MSA can be a valuable tool to identify the coevolution detection method that is most effective in each case.     

Kalign – an accurate and fast multiple sequence alignment algorithm

Background  

The alignment of multiple protein sequences is a fundamental step in the analysis of biological data. It has traditionally been applied to analyzing protein families for conserved motifs, phylogeny, structural properties, and to improve sensitivity in homology searching. The availability of complete genome sequences has increased the demands on multiple sequence alignment (MSA) programs. Current MSA methods suffer from being either too inaccurate or too computationally expensive to be applied effectively in large-scale comparative genomics.     

The partitioned LASSO-patternsearch algorithm with application to gene expression data

Background

The generation of multiple sequence alignments (MSAs) is a crucial step for many bioinformatic analyses. Thus improving MSA accuracy and identifying potential errors in MSAs is important for a wide range of post-genomic research. We present a novel method called MergeAlign which constructs consensus MSAs from multiple independent MSAs and assigns an alignment precision score to each column.

Results

Using conventional benchmark tests we demonstrate that on average MergeAlign MSAs are more accurate than MSAs generated using any single matrix of sequence substitution. We show that MergeAlign column scores are related to alignment precision and hence provide an ab initio method of estimating alignment precision in the absence of curated reference MSAs. Using two novel and independent alignment performance tests that utilise a large set of orthologous gene families we demonstrate that increasing MSA performance leads to an increase in the performance of downstream phylogenetic analyses.

Conclusion

Using multiple tests of alignment performance we demonstrate that this novel method has broad general application in biological research.     

Analysis of <Emphasis Type="Italic">COQ2</Emphasis>gene in multiple system atrophy

Background

Loss of function COQ2 mutations results in primary CoQ10 deficiency. Recently, recessive mutations of the COQ2 gene have been identified in two unrelated Japanese families with multiple system atrophy (MSA). It has also been proposed that specific heterozygous variants in the COQ2 gene may confer susceptibility to sporadic MSA. To assess the frequency of COQ2 variants in patients with MSA, we sequenced the entire coding region and investigated all exonic copy number variants of the COQ2 gene in 97 pathologically-confirmed and 58 clinically-diagnosed MSA patients from the United States.

Results

We did not find any homozygous or compound heterozygous pathogenic COQ2 mutations including deletion or multiplication within our series of MSA patients. In two patients, we identified two heterozygous COQ2 variants (p.S54W and c.403?+?10G?>?T) of unknown significance, which were not observed in 360 control subjects. We also identified one heterozygous carrier of a known loss of function p.S146N substitution in a severe MSA-C pathologically-confirmed patient.

Conclusions

The COQ2 p.S146N substitution has been previously reported as a pathogenic mutation in primary CoQ10 deficiency (including infantile multisystem disorder) in a recessive manner. This variant is the third primary CoQ10 deficiency mutation observed in an MSA case (p.R387X and p.R197H). Therefore it is possible that in the heterozygous state it may increase susceptibility to MSA. Further studies, including reassessing family history in patients of primary CoQ10 deficiency for the possible occurrence of MSA, are now warranted to resolve the role of COQ2 variation in MSA.

     

Urinary Dysfunction in Progressive Supranuclear Palsy Compared with Other Parkinsonian Disorders

Background

Autonomic urinary dysfunction affects patients with progressive supranuclear palsy (PSP); however, the severity and prevalence of urinary dysfunctions in these patients compared with those observed in patients with Parkinson’s disease (PD) and multiple system atrophy (MSA) are unknown.

Objective

We compared urinary dysfunction characteristics in patients with PSP, PD, and MSA.

Patients and Methods

Forty-seven patients who satisfied the probable or possible criteria of the National Institute for Neurological Diseases and Stroke and Society for PSP were assessed using the urinary symptoms questionnaire and the urodynamic study at Chiba and Toho Universities (n = 26 and 21, respectively). The results were compared with those of patients with PD and MSA (n = 218 and 193, respectively).

Results

The mean disease duration of PSP and the mean age were 2.97 ± 0.26 and 71.4 ± 0.88 years, respectively. The mini-mental state examination and frontal assessment battery scores were 22.6 ± 0.70 and 10.7 ± 0.49, respectively. Urinary storage and voiding symptoms were observed in 57% and 56% of patients with PSP, respectively. Detrusor overactivity in the urodynamic study was detected in 81% of patients with PSP, which was slightly more than that found in patients with PD (69%) and MSA (67%); however, this was not statistically significant. Postvoid residual volume in patients with PSP was significantly more than that in patients with PD (P < 0.01), but was equivalent to that in patients with MSA.

Conclusions

The present study demonstrated that patients with PSP experienced various urinary dysfunctions. Urinary storage dysfunction in patients with PSP was not different from that in patients with PD or MSA, whereas urinary voiding dysfunction in patients with PSP was milder than that in patients with MSA and more severe than that in patients with PD. These features should be taken into account for the differentiation of PSP from PD and MSA.     

Involvement of Peripheral Nerves in the Transgenic PLP-α-Syn Model of Multiple System Atrophy: Extending the Phenotype

Multiple system atrophy (MSA) is a fatal, rapidly progressive neurodegenerative disease with (oligodendro-)glial cytoplasmic α-synuclein (α-syn) inclusions (GCIs). Peripheral neuropathies have been reported in up to 40% of MSA patients, the cause remaining unclear. In a transgenic MSA mouse model featuring GCI-like inclusion pathology based on PLP-promoter driven overexpression of human α-syn in oligodendroglia motor and non-motor deficits are associated with MSA-like neurodegeneration. Since α-syn is also expressed in Schwann cells we aimed to investigate whether peripheral nerves are anatomically and functionally affected in the PLP-α-syn MSA mouse model.

Results

To this end, heat/cold as well as mechanical sensitivity tests were performed. Furthermore, in vivo and ex vivo nerve conduction and the G-ratios of the sciatic nerve were analyzed, and thermosensitive ion channel mRNA expression in dorsal root ganglia (DRG) was assessed. The presence of human α-syn in Schwann cells was associated with subtle behavioral impairments. The G-ratio of the sciatic nerve, the conduction velocity of myelinated and unmyelinated primary afferents and the expression of thermosensitive ion channels in the sensory neurons, however, were similar to wildtype mice.

Conclusion

Our results suggest that the PNS appears to be affected by Schwann cell α-syn deposits in the PLP-α-syn MSA mouse model. However, there was no consistent evidence for functional PNS perturbations resulting from such α-syn aggregates suggesting a more central cause of the observed behavioral abnormalities. Nonetheless, our results do not exclude a causal role of α-syn in the pathogenesis of MSA associated peripheral neuropathy.     

Cortical Thickness,Surface Area and Volume Measures in Parkinson's Disease,Multiple System Atrophy and Progressive Supranuclear Palsy

Objective

Parkinson''s disease (PD), Multiple System Atrophy (MSA) and Progressive Supranuclear Palsy (PSP) are neurodegenerative diseases that can be difficult to distinguish clinically. The objective of the current study was to use surface-based analysis techniques to assess cortical thickness, surface area and grey matter volume to identify unique morphological patterns of cortical atrophy in PD, MSA and PSP and to relate these patterns of change to disease duration and clinical features.

Methods

High resolution 3D T1-weighted MRI volumes were acquired from 14 PD patients, 18 MSA, 14 PSP and 19 healthy control participants. Cortical thickness, surface area and volume analyses were carried out using the automated surface-based analysis package FreeSurfer (version 5.1.0). Measures of disease severity and duration were assessed for correlation with cortical morphometric changes in each clinical group.

Results

Results show that in PSP, widespread cortical thinning and volume loss occurs within the frontal lobe, particularly the superior frontal gyrus. In addition, PSP patients also displayed increased surface area in the pericalcarine. In comparison, PD and MSA did not display significant changes in cortical morphology.

Conclusion

These results demonstrate that patients with clinically established PSP exhibit distinct patterns of cortical atrophy, particularly affecting the frontal lobe. These results could be used in the future to develop a useful clinical application of MRI to distinguish PSP patients from PD and MSA patients.     

Improvement in accuracy of multiple sequence alignment using novel group-to-group sequence alignment algorithm with piecewise linear gap cost

Background  

Multiple sequence alignment (MSA) is a useful tool in bioinformatics. Although many MSA algorithms have been developed, there is still room for improvement in accuracy and speed. In the alignment of a family of protein sequences, global MSA algorithms perform better than local ones in many cases, while local ones perform better than global ones when some sequences have long insertions or deletions (indels) relative to others. Many recent leading MSA algorithms have incorporated pairwise alignment information obtained from a mixture of sources into their scoring system to improve accuracy of alignment containing long indels.     

Mesenchymal stem cells in a transgenic mouse model of multiple system atrophy: immunomodulation and neuroprotection

Background

Mesenchymal stem cells (MSC) are currently strong candidates for cell-based therapies. They are well known for their differentiation potential and immunoregulatory properties and have been proven to be potentially effective in the treatment of a large variety of diseases, including neurodegenerative disorders. Currently there is no treatment that provides consistent long-term benefits for patients with multiple system atrophy (MSA), a fatal late onset α-synucleinopathy. Principally neuroprotective or regenerative strategies, including cell-based therapies, represent a powerful approach for treating MSA. In this study we investigated the efficacy of intravenously applied MSCs in terms of behavioural improvement, neuroprotection and modulation of neuroinflammation in the (PLP)-αsynuclein (αSYN) MSA model.

Methodology/Principal Findings

MSCs were intravenously applied in aged (PLP)-αSYN transgenic mice. Behavioural analyses, defining fine motor coordination and balance capabilities as well as stride length analysis, were performed to measure behavioural outcome. Neuroprotection was assessed by quantifying TH neurons in the substantia nigra pars compacta (SNc). MSC treatment on neuroinflammation was analysed by cytokine measurements (IL-1α, IL-2, IL-4, IL-5, IL-6, IL-10, IL-17, GM-CSF, INFγ, MCP-1, TGF-β1, TNF-α) in brain lysates together with immunohistochemistry for T-cells and microglia.Four weeks post MSC treatment we observed neuroprotection in the SNc, as well as downregulation of cytokines involved in neuroinflammation. However, there was no behavioural improvement after MSC application.

Conclusions/Significance

To our knowledge this is the first experimental approach of MSC treatment in a transgenic MSA mouse model. Our data suggest that intravenously infused MSCs have a potent effect on immunomodulation and neuroprotection. Our data warrant further studies to elucidate the efficacy of systemically administered MSCs in transgenic MSA models.     

Symptoms and Quality of Life in Late Stage Parkinson Syndromes: A Longitudinal Community Study of Predictive Factors

Background

Palliative care is increasingly offered earlier in the cancer trajectory but rarely in Idiopathic Parkinson''s Disease(IPD), Progressive Supranuclear Palsy(PSP) or Multiple System Atrophy(MSA). There is little longitudinal data of people with late stage disease to understand levels of need. We aimed to determine how symptoms and quality of life of these patients change over time; and what demographic and clinical factors predicted changes.

Methods

We recruited 82 patients into a longitudinal study, consenting patients with a diagnosis of IPD, MSA or PSP, stages 3–5 Hoehn and Yahr(H&Y). At baseline and then on up to 3 occasions over one year, we collected self-reported demographic, clinical, symptom, palliative and quality of life data, using Parkinson''s specific and generic validated scales, including the Palliative care Outcome Scale (POS). We tested for predictors using multivariable analysis, adjusting for confounders.

Findings

Over two thirds of patients had severe disability, over one third being wheelchair-bound/bedridden. Symptoms were highly prevalent in all conditions - mean (SD) of 10.6(4.0) symptoms. More than 50% of the MSA and PSP patients died over the year. Over the year, half of the patients showed either an upward (worsening, 24/60) or fluctuant (8/60) trajectory for POS and symptoms. The strongest predictors of higher levels of symptoms at the end of follow-up were initial scores on POS (AOR 1.30; 95%CI:1.05–1.60) and being male (AOR 5.18; 95% CI 1.17 to 22.92), both were more predictive than initial H&Y scores.

Interpretation

The findings point to profound and complex mix of non-motor and motor symptoms in patients with late stage IPD, MSA and PSP. Symptoms are not resolved and half of the patients deteriorate. Palliative problems are predictive of future symptoms, suggesting that an early palliative assessment might help screen for those in need of earlier intervention.     

Global variation in the relationship between avian phylogenetic diversity and functional distance is driven by environmental context and constraints

Aim

If evolutionary distance is akin to evolutionary chance, then it follows that species assemblages that are distantly related will also be more disparate in terms of their traits, features and the niches they occupy. Yet, studies have found that the total phylogenetic distance of an assemblage, known as phylogenetic diversity, is an unreliable surrogate for functional diversity. We investigate global variation in the relationship between Faith's phylogenetic diversity (PD) and mean pairwise functional distance (MPFD) across latitude and the influence of migratory species on both these aspects of diversity.

Location

Global.

Time period

Present day.

Major taxa studied

Birds.

Methods

We measure PD and MPFD for over 9000 species of bird across more than 17,000 globally distributed assemblages. We obtain standardised effect sizes for both indices by simulating assemblage composition under an ecologically informed null model. We employ path analysis to characterize variation in the relationship between PD's and MPFD across latitude, elevation and with proportion of migratory species.

Results

Globally, assemblages that were phylogenetically diverse tended to be less functionally dispersed than expected; however, this relationship showed considerable variation across latitude decreasing with distance from the equator. The proportion of migratory species in an assemblage was found to be an important predictor of functional diversity, with migrant rich assemblages generally showing less functional diversity than expected. We identify the Andes and Hengduan Mountains as regions of exceptional bird functional diversity.

Main conclusions

The relationship between phylogenetic diversity and function diversity is context specific, varying across environmental gradients such as latitude, and influenced by ecological phenomena such as migration. Thus, care should be taken using phylogenetic diversity as a proxy for functional diversity, particularly in clades with sparse functional data. Instead, we recommend that studies consider how phylogenetic diversity's surrogacy for functional diversity may be impacted by environmental context and evaluate empirical observations against biogeographically constrained and ecologically informed null models.