Stable Isotope Labeling by Amino Acids in Cell Culture (SILAC) and Quantitative Comparison of the Membrane Proteomes of Self-renewing and Differentiating Human Embryonic Stem Cells

Stable isotope labeling by amino acids in cell culture (SILAC) is a powerful quantitative proteomics platform for comprehensive characterization of complex biological systems. However, the potential of SILAC-based approaches has not been fully utilized in human embryonic stem cell (hESC) research mainly because of the complex nature of hESC culture conditions. Here we describe complete SILAC labeling of hESCs with fully preserved pluripotency, self-renewal capabilities, and overall proteome status that was quantitatively analyzed to a depth of 1556 proteins and 527 phosphorylation events. SILAC-labeled hESCs appear to be perfectly suitable for functional studies, and we exploited a SILAC-based proteomics strategy for discovery of hESC-specific surface markers. We determined and quantitatively compared the membrane proteomes of the self-renewing versus differentiating cells of two distinct human embryonic stem cell lines. Of the 811 identified membrane proteins, six displayed significantly higher expression levels in the undifferentiated state compared with differentiating cells. This group includes the established marker CD133/Prominin-1 as well as novel candidates for hESC surface markers: Glypican-4, Neuroligin-4, ErbB2, receptor-type tyrosine-protein phosphatase ζ (PTPRZ), and Glycoprotein M6B. Our study also revealed 17 potential markers of hESC differentiation as their corresponding protein expression levels displayed a dramatic increase in differentiated embryonic stem cell populations.Human embryonic stem cells (hESCs)¹ are stem cells derived from the blastocyst inner cell mass. They are pluripotent; thus they are able to differentiate into any human cell type. The self-renewal capacity and pluripotency make hESCs an ideal system to study the processes of cell development and differentiation. Moreover hESC research is highly relevant for regenerative medicine, which aims at replacing or restoring tissue damaged by disease or injury through transplantation of functional hESCs (1,2). However, factors responsible for maintaining the undifferentiated and pluripotent nature of hESCs are still largely unknown. Before hESCs can be used for transplantation into the human body, reliable and reproducible protocols for differentiating them into specific cell types are needed. To create such protocols we need to develop a thorough understanding of the mechanisms maintaining the undifferentiated pluripotent nature of hESCs and those guiding their differentiation into specific lineages.A number of factors involved in the maintenance of pluripotency have been described over the last few years (3). It has also been demonstrated that overexpression of some of these factors in somatic cells is sufficient to turn them into pluripotent stem cells very similar to hESCs (4–8). However, it is apparent that the processes occurring during such transformation are extremely complex. A large number of factors and pathways are involved in maintaining the pluripotent state and regulating self-renewal and differentiation. The process of specific hESC differentiation into distinct cell types is even less understood. Most current attempts to directionally differentiate hESCs are based on sequential application of empirically selected growth factors and consequent selection for markers expressed in the target cell types (9). A more systematic approach is needed to improve our understanding of the pathways that control the conversion of precursors into specific cell types, progressing toward the goal of reproducing these processes in vitro for the generation of functional cells and tissues for transplantation.Comprehensive quantitative analysis of the hESC proteome would mean an important advance in understanding the nature of “stemness,” pluripotency, and differentiation. Several studies targeting various aspects of the hESC proteome have already been reported (for reviews, see Refs. 10 and 11). The task, however, is so enormous that further detailed analysis and novel strategies are necessary and will be of high interest and importance. In this regard, MS-based quantitative proteomics and in particular stable isotope labeling by amino acids in cell culture (SILAC) may greatly facilitate the process of defining the mechanisms of hESC self-renewal and differentiation. With SILAC, the entire proteome of a given cell population is metabolically labeled by heavy, non-radioactive isotopic variants of amino acids, thus making it distinguishable by MS analysis (12). Thereafter two or more distinctly SILAC-labeled cell populations can be mixed and analyzed in one MS experiment that allows accurate quantitation of proteins from the different cellular states (13). This versatile strategy has been demonstrated to be very useful for comprehensive characterization of complex biological phenomena (14–21) including in-depth comparison of signaling pathways to identify control points determining cell fate of adult mesenchymal stem cells (22).Here we report a procedure for complete SILAC labeling of human ES cells. We show that these SILAC-encoded hESCs have preserved self-renewing undifferentiated status as well as pluripotent capabilities based on analysis of known markers. In addition, we further compared the overall proteomes and phosphoproteomes of SILAC-labeled hESCs and equivalent cells grown under conventional culture conditions. We next compared the membrane proteomes of undifferentiated and differentiated hESCs in a quantitative manner. Our analysis identified 811 membrane proteins, which to our knowledge is the largest data set of ES cell membrane proteome. This study also revealed 23 membrane proteins with large changes in their expression levels during the differentiation. Six of those cell surface molecules displayed more than 3-fold higher levels in the self-renewing cells, whereas the remaining 17 were identified as more abundant in the differentiated population. These may be useful as specific hESC markers for the corresponding ES cell state and help to shed light on the mechanisms for self-renewal and differentiation.     

A simple probabilistic model of multibody interactions in proteins

Protein structure prediction methods typically use statistical potentials, which rely on statistics derived from a database of know protein structures. In the vast majority of cases, these potentials involve pairwise distances or contacts between amino acids or atoms. Although some potentials beyond pairwise interactions have been described, the formulation of a general multibody potential is seen as intractable due to the perceived limited amount of data. In this article, we show that it is possible to formulate a probabilistic model of higher order interactions in proteins, without arbitrarily limiting the number of contacts. The success of this approach is based on replacing a naive table‐based approach with a simple hierarchical model involving suitable probability distributions and conditional independence assumptions. The model captures the joint probability distribution of an amino acid and its neighbors, local structure and solvent exposure. We show that this model can be used to approximate the conditional probability distribution of an amino acid sequence given a structure using a pseudo‐likelihood approach. We verify the model by decoy recognition and site‐specific amino acid predictions. Our coarse‐grained model is compared to state‐of‐art methods that use full atomic detail. This article illustrates how the use of simple probabilistic models can lead to new opportunities in the treatment of nonlocal interactions in knowledge‐based protein structure prediction and design. Proteins 2013; 81:1340–1350. © 2013 Wiley Periodicals, Inc.     

Coordinate 5' and 3' endonucleolytic trimming of terminally blocked blunt DNA double-strand break ends by Artemis nuclease and DNA-dependent protein kinase

Previous work showed that, in the presence of DNA-dependent protein kinase (DNA-PK), Artemis slowly trims 3′-phosphoglycolate-terminated blunt ends. To examine the trimming reaction in more detail, long internally labeled DNA substrates were treated with Artemis. In the absence of DNA-PK, Artemis catalyzed extensive 5′→3′ exonucleolytic resection of double-stranded DNA. This resection required a 5′-phosphate, but did not require ATP, and was accompanied by endonucleolytic cleavage of the resulting 3′ overhang. In the presence of DNA-PK, Artemis-mediated trimming was more limited, was ATP-dependent and did not require a 5′-phosphate. For a blunt end with either a 3′-phosphoglycolate or 3′-hydroxyl terminus, endonucleolytic trimming of 2–4 nucleotides from the 3′-terminal strand was accompanied by trimming of 6 nt from the 5′-terminal strand. The results suggest that autophosphorylated DNA-PK suppresses the exonuclease activity of Artemis toward blunt-ended DNA, and promotes slow and limited endonucleolytic trimming of the 5′-terminal strand, resulting in short 3′ overhangs that are trimmed endonucleolytically. Thus, Artemis and DNA-PK can convert terminally blocked DNA ends of diverse geometry and chemical structure to a form suitable for polymerase-mediated patching and ligation, with minimal loss of terminal sequence. Such processing could account for the very small deletions often found at DNA double-strand break repair sites.     

Amyloid‐beta 1‐40 is associated with alterations in NG2+ pericyte population ex vivo and in vitro

The population of brain pericytes, a cell type important for vessel stability and blood brain barrier function, has recently been shown altered in patients with Alzheimer's disease (AD). The underlying reason for this alteration is not fully understood, but progressive accumulation of the AD characteristic peptide amyloid‐beta (Aβ) has been suggested as a potential culprit. In the current study, we show reduced number of hippocampal NG2+ pericytes and an association between NG2+ pericyte numbers and Aβ1‐40 levels in AD patients. We further demonstrate, using in vitro studies, an aggregation‐dependent impact of Aβ1‐40 on human NG2+ pericytes. Fibril‐EP Aβ1‐40 exposure reduced pericyte viability and proliferation and increased caspase 3/7 activity. Monomer Aβ1‐40 had quite the opposite effect: increased pericyte viability and proliferation and reduced caspase 3/7 activity. Oligomer‐EP Aβ1‐40 had no impact on either of the cellular events. Our findings add to the growing number of studies suggesting a significant impact on pericytes in the brains of AD patients and suggest different aggregation forms of Aβ1‐40 as potential key regulators of the brain pericyte population size.     

ATM-dependent ERK signaling via AKT in response to DNA double-strand breaks

Ionizing radiation (IR) triggers many signaling pathways primarily originating from either damaged DNA or non-nuclear sources such as growth factor receptors. Thus, to study the DNA damage-induced signaling component alone by irradiation would be a challenge. To generate DNA double-strand breaks (DSBs) and minimize non-nuclear signaling, human cancer cells having bromodeoxyuridine (BrdU)—substituted DNA were treated with the photosensitizer Hoechst 33258 followed by long wavelength UV (UV-A) treatment (BrdU photolysis). BrdU photolysis resulted in well-controlled, dose-dependent generation of DSBs equivalent to radiation doses between 0.2–20 Gy, as determined by pulsed-field gel electrophoresis and accompanied by dose-dependent ATM (ser-1981), H2AX (ser-139), Chk2 (thr-68) and p53 (ser15) phosphorylation. Interestingly, low levels (≤2 Gy equivalents) of BrdU photolysis stimulated ERK phosphorylation whereas higher (>2 Gy eq.) resulted in ERK dephosphorylation. ERK phosphorylation was ATM-dependent whereas dephosphorylation was ATM-independent. The ATM-dependent increase in ERK phosphorylation was also seen when DSBs were generated by transfection of cells with an EcoRI expression plasmid or by electroporation of EcoRI enzyme. Furthermore, AKT was critical for transmitting the DSB signal to ERK. Altogether, our results show that low levels of DSBs trigger ATM- and AKT-dependent ERK pro-survival signaling and increased cell proliferation whereas higher levels result in ERK dephosphorylation consistent with a dose-dependent switch from pro-survival to anti-survival signaling.Key words: bromodeoxyuridine, DNA repair, MAP kinase, p53, KU-55933, U87 glioma cells     

Production of diaspores at the landscape level regulates local colonization: an experiment with a spore‐dispersed moss

Effective dispersal is crucial to species inhabiting transient substrates in order for them to be able to persist in a landscape. Bryophytes, pteridophytes, lichens and fungi all have wind‐dispersed small diaspores and can be efficiently dispersed if their diaspores reach air masses above canopy height. However, empirical data on dispersal over landscape scales are scarce. We investigated how the colonization of an acrocarpous clay‐inhabiting pioneer moss, Discelium nudum, varied between sites that differed in connectivity to potential dispersal sources at spatial scales from 1 to 20 km in a region in northern Sweden. We recorded the colonization on ?25 introduced clay heaps at each of 14 experimental sites some months after the dispersal period. The colonization rate ranged from 0–82% and had a statistically significant relationship with a proxy for potential habitats (amount of clay‐dominated soil) in a buffer of 20 km radius surrounding the experimental sites (and also weakly with the amount of substrate in a 10 km buffer). There were no significant relationships between colonization rate and connectivity at smaller scales (1 and 5 km). We made a rough estimate of the number of spores available for dispersal in a landscape, given the amount of clay‐dominated soil, by recording the number of Discelium nudum colonies in two 25 × 25 km landscapes. The estimated available spore numbers in the different 20 km buffers were of the same order of magnitude as the deposition densities at the experimental sites calculated from the colonization rates. The results suggest that the spores of species with scattered occurrences and small diaspores (25 μm) in open landscapes can be deposited over extensive areas, at rates high enough to drive colonization patterns. This also implies that regional connectivity may be more important than local connectivity for these kinds of species.     

Novel SLC19A3 Promoter Deletion and Allelic Silencing in Biotin-Thiamine-Responsive Basal Ganglia Encephalopathy

Background

Biotin-thiamine responsive basal ganglia disease is a severe, but potentially treatable disorder caused by mutations in the SLC19A3 gene. Although the disease is inherited in an autosomal recessive manner, patients with typical phenotypes carrying single heterozygous mutations have been reported. This makes the diagnosis uncertain and may delay treatment.

Methods and Results

In two siblings with early-onset encephalopathy dystonia and epilepsy, whole-exome sequencing revealed a novel single heterozygous SLC19A3 mutation (c.337T>C). Although Sanger-sequencing and copy-number analysis revealed no other aberrations, RNA-sequencing in brain tissue suggested the second allele was silenced. Whole-genome sequencing resolved the genetic defect by revealing a novel 45,049 bp deletion in the 5’-UTR region of the gene abolishing the promoter. High dose thiamine and biotin therapy was started in the surviving sibling who remains stable. In another patient two novel compound heterozygous SLC19A3 mutations were found. He improved substantially on thiamine and biotin therapy.

Conclusions

We show that large genomic deletions occur in the regulatory region of SLC19A3 and should be considered in genetic testing. Moreover, our study highlights the power of whole-genome sequencing as a diagnostic tool for rare genetic disorders across a wide spectrum of mutations including non-coding large genomic rearrangements.     

Modifications of the 7-Hydroxyl Group of the Transthyretin Ligand Luteolin Provide Mechanistic Insights into Its Binding Properties and High Plasma Specificity

Amyloid formation of the plasma protein transthyretin (TTR) has been linked to familial amyloid polyneuropathy and senile systemic amyloidosis. Binding of ligands within its natural hormone binding site can stabilize the tetrameric structure and impair amyloid formation. We have recently shown that the flavonoid luteolin stabilizes TTR in human plasma with a very high selectivity. Luteolin, however, is inactivated in vivo via glucuronidation for which the preferred site is the hydroxy group at position 7 on its aromatic A-ring. We have evaluated the properties of two luteolin variants in which the 7-hydroxy group has been exchanged for a chlorine (7-Cl-Lut) or a methoxy group (7-MeO-Lut). Using an in vitro model, based on human liver microsomes, we verified that these modifications increase the persistence of the drug. Crystal structure determinations show that 7-Cl-Lut binds similarly to luteolin. The larger MeO substituent cannot be accommodated within the same space as the chlorine or hydroxy group and as a result 7-MeO-Lut binds in the opposite direction with the methoxy group in position 7 facing the solvent. Both 7-Cl-Lut and 7-MeO-Lut qualify as high-affinity binders, but in contrast to luteolin, they display a highly non-specific binding to other plasma components. The binding of the two conformations and the key-interactions to TTR are discussed in detail. Taken together, these results show a proof-of-concept that the persistence of luteolin towards enzymatic modification can be increased. We reveal two alternative high-affinity binding modes of luteolin to TTR and that modification in position 7 is restricted only to small substituents if the original orientation of luteolin should be preserved. In addition, the present work provides a general and convenient method to evaluate the efficacy of TTR-stabilizing drugs under conditions similar to an in vivo environment.     

Tuning of Thioredoxin Redox Properties by Intramolecular Hydrogen Bonds

Thioredoxin-like proteins contain a characteristic C-x-x-C active site motif and are involved in a large number of biological processes ranging from electron transfer, cellular redox level maintenance, and regulation of cellular processes. The mechanism for deprotonation of the buried C-terminal active site cysteine in thioredoxin, necessary for dissociation of the mixed-disulfide intermediate that occurs under thiol/disulfide mediated electron transfer, is not well understood for all thioredoxin superfamily members. Here we have characterized a 8.7 kD thioredoxin (BC3987) from Bacillus cereus that unlike the typical thioredoxin appears to use the conserved Thr8 side chain near the unusual C-P-P-C active site to increase enzymatic activity by forming a hydrogen bond to the buried cysteine. Our hypothesis is based on biochemical assays and thiolate pK_a titrations where the wild type and T8A mutant are compared, phylogenetic analysis of related thioredoxins, and QM/MM calculations with the BC3987 crystal structure as a precursor for modeling of reduced active sites. We suggest that our model applies to other thioredoxin subclasses with similar active site arrangements.