Statistical discovery of site inter-dependencies in sub-molecular hierarchical protein structuring

Background

Much progress has been made in understanding the 3D structure of proteins using methods such as NMR and X-ray crystallography. The resulting 3D structures are extremely informative, but do not always reveal which sites and residues within the structure are of special importance. Recently, there are indications that multiple-residue, sub-domain structural relationships within the larger 3D consensus structure of a protein can be inferred from the analysis of the multiple sequence alignment data of a protein family. These intra-dependent clusters of associated sites are used to indicate hierarchical inter-residue relationships within the 3D structure. To reveal the patterns of associations among individual amino acids or sub-domain components within the structure, we apply a k-modes attribute (aligned site) clustering algorithm to the ubiquitin and transthyretin families in order to discover associations among groups of sites within the multiple sequence alignment. We then observe what these associations imply within the 3D structure of these two protein families.

Results

The k-modes site clustering algorithm we developed maximizes the intra-group interdependencies based on a normalized mutual information measure. The clusters formed correspond to sub-structural components or binding and interface locations. Applying this data-directed method to the ubiquitin and transthyretin protein family multiple sequence alignments as a test bed, we located numerous interesting associations of interdependent sites. These clusters were then arranged into cluster tree diagrams which revealed four structural sub-domains within the single domain structure of ubiquitin and a single large sub-domain within transthyretin associated with the interface among transthyretin monomers. In addition, several clusters of mutually interdependent sites were discovered for each protein family, each of which appear to play an important role in the molecular structure and/or function.

Conclusions

Our results demonstrate that the method we present here using a k- modes site clustering algorithm based on interdependency evaluation among sites obtained from a sequence alignment of homologous proteins can provide significant insights into the complex, hierarchical inter-residue structural relationships within the 3D structure of a protein family.

     

MTHFR and MTRR genotype and haplotype analysis and colorectal cancer susceptibility in a case-control study from the Czech Republic

Polymorphic variants in genes involved in one-carbon metabolism, in particular of dietary folate, may modulate the risk for colorectal cancer through aberrant DNA-methylation and altered nucleotide synthesis and repair. In the present study, we have assessed the association of six polymorphisms and relative haplotypes in the MTHFR gene (rs1801133 and rs1801131) and in the MTRR gene (rs1801394, rs1532268, rs162036, and rs10380) with the risk for colorectal cancer in 666 patients and 1377 controls from the Czech Republic. We found that the 677 C>T polymorphism in the MTHFR gene significantly decreased the risk for colorectal cancer in homozygous carriers of the variant allele (OR, 0.58; 95% CI, 0.39-0.87). Also, we noted a significantly different distribution of genotypes between cases and controls for the 66A>G polymorphism in the MTRR gene. In particular, homozygous carriers of the G-containing allele of this polymorphism were at an increased risk for colorectal cancer (OR, 1.39; 95% CI, 1.04-1.85). Haplotype analysis of the two MTHFR polymorphisms showed a moderate difference in the distribution of the TA haplotype between cases and controls. In comparison to the most common haplotype (CA), the TA haplotype was associated with a decreased risk for colorectal cancer (OR, 0.84; 95% CI, 0.71-0.99). No difference in the distribution between cases and controls was observed for the haplotypes based on the four polymorphisms in the MTRR gene. The present study suggests that the 677TT genotype and the TA haplotype in the MTHFR gene may also have a role in colorectal cancer risk in the Czech population, indicating the importance of genes involved in folate metabolism with respect to cancer risk. For MTRR, additional studies on larger populations are needed to clarify the possible role of variation in this gene in colorectal carcinogenesis.     

Cell cycle phosphorylation of mitotic exit network (MEN) proteins

Phosphorylation of proteins is an important mechanism used to regulate most cellular processes. Recently, we completed an extensive phosphoproteomic analysis of the core proteins that constitute the Saccharomyces cerevisiae centrosome. Here, we present a study of phosphorylation sites found on the mitotic exit network (MEN) proteins, most of which are associated with the cytoplasmic face of the centrosome. We identified 55 sites on Bfa1, Cdc5, Cdc14 and Cdc15. Eight sites lie in cyclin-dependent kinase motifs (Cdk, S/T-P), and 22 sites are completely conserved within fungi. More than half of the sites were found in centrosomes from mitotic cells, possibly in preparation for their roles in mitotic exit. Finally, we report phosphorylation site information for other important cell cycle and regulatory proteins.Key words: in vivo phosphorylation, yeast centrosome, mitotic exit network (MEN), cell cycle, protein kinase, Cdk (cyclin-dependent kinase)/Cdc28, Plk1 (polo-like kinase)/Cdc5Reversible protein phosphorylation leads to changes in targeting, structure and stability of proteins and is used widely to modulate biochemical reactions in the cell. We are interested in phosphoregulation of centrosome duplication and function in the yeast Saccharomyces cerevisiae. Centrosomes nucleate microtubules and, upon duplication during the cell cycle, form the two poles of the bipolar mitotic spindle used to segregate replicated chromosomes into the two daughter cells. Timing and spatial cues are highly regulated to ensure that elongation of the mitotic spindle and separation of sister chromatids occur prior to progression into late telophase and initiation of mitotic exit. The mitotic exit network (MEN) regulates this timing through a complex signaling cascade activated at the centrosome that triggers the end of mitosis, ultimately through mitotic cyclin-dependent kinase (Cdk) inactivation (reviewed in ref. ¹).The major components of the MEN pathway (Fig. 1) are a Ras-like GTPase (Tem1), an activator (Lte1) with homology to nucleotide exchange factors, a GTPase-activating protein (GAP) complex (Bfa1/Bub2), several protein kinases [Cdc5 (Plk1 in humans), Cdc15 and Dbf2/Mob1] and Cdc14 phosphatase (reviewed in ref. ^2–5). Tem1 initiates the signal for the MEN pathway when switched to a GTP-active state. Prior to activation at anaphase, it is held at the centrosome in an inactive GDP-bound state by an inhibiting GAP complex, Bfa1/Bub2.⁶ The Bfa1/Bub2 complex and the inactive Tem1 are localized at the mother centrosome destined to move into the budded cell upon chromosome segregation, whereas the activator Lte1 is localized at the tip of the budded cell. These separate localizations ensure that Lte1 and Tem1 only interact in late anaphase, when the mitotic spindle elongates.^7,8 Lte1 has been thought to activate Tem1 as a nucleotide exchange factor, although more recent evidence suggests that it may instead affect Bfa1 localization.⁹ In addition, full activation of Tem1 is achieved through Cdc5 phosphorylation of the negative regulator Bfa1 ¹⁰ and potentially through phosphorylation of Lte1. GTP-bound Tem1 is then able to recruit Cdc15 to the centrosome, allowing for Dbf2 activation.³ The final step in the MEN pathway is release of Cdc14 from the nucleolus, which is at least partially due to phosphorylation by Dbf2¹¹ an leads to mitotic cyclin degradation and inactivation of the mitotic kinase.²Open in a separate windowFigure 1Schematic representation of the MEN proteins and pathway. MEN protein localization is shown within a metaphase cell when mitotic exit is inhibited and in a late anaphase cell when mitotic exit is initiated. Primary inhibition and activation events are described below the cells.Recently, we performed a large-scale analysis of phosphorylation sites on the 18 core yeast centrosomal proteins present in enriched centrosomal preparations.¹² In total, we mapped 297 sites on 17 of the 18 proteins and described their cell cycle regulation, levels of conservation and demonstrated defects in centrosome assembly and function resulting from mutating selected sites. MEN proteins were also identified in the centrosome preparations. This was expected, because Nud1, one of the 18 core centrosome components, is known to recruit several MEN proteins to the centrosome¹³ as part of its function in mitotic exit.^14,15 As phosphorylation is essential to several steps in the MEN pathway, beginning with recruitment of Bfa1/Bub2 by phosphorylated Nud1,¹⁵ we were interested in mapping in vivo phosphorylation sites on the MEN proteins associated with centrosomes and identifying when they occur during the cell cycle.We combined centrosome enrichment with mass spectrometry analysis to examine phosphorylation from asynchronously growing cells.¹² Centrosomes were also prepared from cells arrested in G₁ and mitosis¹² to monitor potentially cell cycle-regulated sites. We obtained significant coverage of a number of the MEN proteins, several of which have human homologs (​and33, column 1), of which eight sites lie within Cdk/Cdc28 motifs [S/T(P)], (23 Mob1 and Dbf2 are known phosphoproteins²⁴ for which we observed peptide coverage but no phosphorylation. Surprisingly, we did not detect phosphorylation on Bub2 despite the high peptide coverage; it is possible that the mitotic centrosome preparations (using a Cdc20 depletion protocol) affect the phosphorylation state of Bub2, as Bub2 is required for mitotic exit arrest in cdc20 mutants.²⁵ Additionally, specific phosphorylation sites have not been mapped on Bub2, suggesting that modifications on this protein may be difficult to observe by mass spectrometry. Lte1 does not localize to the centrosome, and we did not recover Lte1 peptides in our preparations. Many phosphorylation events on MEN proteins were observed in mitotic centrosomal preparations, most likely in preparation for their subsequent role in exit from mitosis (

7 ? resolution in protein two-dimensional-crystal X-ray diffraction at Linac Coherent Light Source

Membrane proteins arranged as two-dimensional crystals in the lipid environment provide close-to-physiological structural information, which is essential for understanding the molecular mechanisms of protein function. Previously, X-ray diffraction from individual two-dimensional crystals did not represent a suitable investigational tool because of radiation damage. The recent availability of ultrashort pulses from X-ray free-electron lasers (XFELs) has now provided a means to outrun the damage. Here, we report on measurements performed at the Linac Coherent Light Source XFEL on bacteriorhodopsin two-dimensional crystals mounted on a solid support and kept at room temperature. By merging data from about a dozen single crystal diffraction images, we unambiguously identified the diffraction peaks to a resolution of 7 Å, thus improving the observable resolution with respect to that achievable from a single pattern alone. This indicates that a larger dataset will allow for reliable quantification of peak intensities, and in turn a corresponding increase in the resolution. The presented results pave the way for further XFEL studies on two-dimensional crystals, which may include pump–probe experiments at subpicosecond time resolution.     

Infection of human THP-1 cells with dormant Mycobacterium tuberculosis

Dormant, non-replicating Mycobacterium tuberculosis H37Rv strain cultured in hypoxic conditions was used to infect THP-1 cells. CFUs counting, Kinyoun staining and electron microscopy showed that dormant bacilli infected THP-1 cells at a rate similar to replicating M. tuberculosis, but failed to grow during the first 6 days of infection. The absence of growth was specific to the intracellular compartment, as demonstrated by efficient growth in liquid medium. Quantification of β-actin mRNA recovered from infected cells showed that, in contrast with log-phase bacteria, infection with dormant bacilli determined a reduced THP-1 cell death. Gene expression of intracellular non-replicating bacteria showed a pattern typical of a dormant state. Intracellular dormant bacteria induced the activation of genes associated to a proinflammatory response in THP-1 cells. Though, higher levels of TNFα, IL-1β and IL-8 mRNAs compared to aerobic H37Rv infected cells were not paralleled by increased cytokine accumulation in the supernatants. Moreover, dormant bacilli induced a higher expression of inducible cox-2 gene, accompanied by increased PGE2 secretion. Overall, our data describe a new model of in vitro infection using dormant M. tuberculosis that could provide the basis for understanding how non-replicating bacilli survive intracellularly and influence the maintenance of the hypoxic granuloma.     

Frequency of FMR1 premutation in individuals with ataxia and/or tremor and/or parkinsonism

A late onset neurological syndrome in carriers of premutation in FMR1 gene was recently described. The condition was named fragile-X-associated tremor/ataxia syndrome (FXTAS) and includes intentional tremor, cerebellar ataxia, parkinsonism, and cognitive deficit. We ascertained the contribution of FMR1 premutation to the phenotypes ataxia, tremor and/or parkinsonism. Sixty-six men over 45 years old presenting these symptoms, isolated or combined, were tested. Also, 74 normal men, randomly chosen in the population, formed the control group. In the patient group, no premutation carrier was found, which is in agreement with other observed frequencies reported elsewhere (0-5% variation). No significant differences were found when comparing gray zone allele frequencies among target and control groups. The FXTAS contribution in patients with phenotypic manifestations of FXTAS was 15/748 (2%). The presence of gray zone alleles is not correlated with FXTAS occurrence.     

Polymorphic dinucleotide microsatellite loci in the clonal ascidian Diplosoma listerianum: predominance of compound and highly interrupted imperfect repeats

We report the isolation of the first dinucleotide microsatellite loci from the clonal ascidian, Diplosoma listerianum. Most repeats were compound and highly interrupted, with flanking sequences containing many small interspersed repetitive regions. Consequently, most primers detected polymerase chain reaction products outside the expected size range, and only five out of 15 primers detected polymorphic single‐locus markers. Characterization of five variable loci from two UK populations revealed three to seven alleles per locus, with observed heterozygosity of 0.00–0.86 and expected heterozygosity of 0.10–0.66. Three loci showed significant heterozygote deficits either because of inbreeding, population substructure or the presence of null alleles.     

Importance of estimating matrix quality for modeling species distribution in complex tropical landscapes: a test with Atlantic forest small mammals

Information to guide decision making is especially urgent in human dominated landscapes in the tropics, where urban and agricultural frontiers are still expanding in an unplanned manner. Nevertheless, most studies that have investigated the influence of landscape structure on species distribution have not considered the heterogeneity of altered habitats of the matrix, which is usually high in human dominated landscapes. Using the distribution of small mammals in forest remnants and in the four main altered habitats in an Atlantic forest landscape, we investigated 1) how explanatory power of models describing species distribution in forest remnants varies between landscape structure variables that do or do not incorporate matrix quality and 2) the importance of spatial scale for analyzing the influence of landscape structure. We used standardized sampling in remnants and altered habitats to generate two indices of habitat quality, corresponding to the abundance and to the occurrence of small mammals. For each remnant, we calculated habitat quantity and connectivity in different spatial scales, considering or not the quality of surrounding habitats. The incorporation of matrix quality increased model explanatory power across all spatial scales for half the species that occurred in the matrix, but only when taking into account the distance between habitat patches (connectivity). These connectivity models were also less affected by spatial scale than habitat quantity models. The few consistent responses to the variation in spatial scales indicate that despite their small size, small mammals perceive landscape features at large spatial scales. Matrix quality index corresponding to species occurrence presented a better or similar performance compared to that of species abundance. Results indicate the importance of the matrix for the dynamics of fragmented landscapes and suggest that relatively simple indices can improve our understanding of species distribution, and could be applied in modeling, monitoring and managing complex tropical landscapes.     

An easy and reliable method for establishment and maintenance of leaf and root cell cultures ofArabidopsis thaliana

Cell suspension cultures are useful for a wide range of biochemical and physiological studies, yet their production can be technically demanding and often unreliable. Here we describe a protocol for producing Arabidopsis cell suspension cultures that is reliable and easy to use.