Quantitative functions of Argonaute proteins in mammalian development.

Argonaute proteins (Ago1-4) are essential components of the microRNA-induced silencing complex and play important roles in both microRNA biogenesis and function. Although Ago2 is the only one with the slicer activity, it is not clear whether the slicer activity is a universally critical determinant for Ago2's function in mammals. Furthermore, functional specificities associated with different Argonautes remain elusive. Here we report that microRNAs are randomly sorted to individual Argonautes in mammals, independent of the slicer activity. When both Ago1 and Ago2, but not either Ago1 or Ago2 alone, are ablated in the skin, the global expression of microRNAs is significantly compromised and it causes severe defects in skin morphogenesis. Surprisingly, Ago3 is able to load microRNAs efficiently in the absence of Ago1 and Ago2, despite a significant loss of global microRNA expression. Quantitative analyses reveal that Ago2 interacts with a majority of microRNAs (60%) in the skin, compared with Ago1 (30%) and Ago3 (<10%). This distribution is highly correlated with the abundance of each Argonaute, as quantified by shotgun proteomics. The quantitative correlation between Argonautes and their associated microRNAs is conserved in human cells. Finally, we measure the absolute expression of Argonaute proteins and determine that their copy number is ~1.4 × 10(5) to 1.7 × 10(5) molecules per cell. Together, our results reveal a quantitative picture for microRNA activity in mammals.

Functional proteomics identifies targets of phosphorylation by B-Raf signaling in melanoma.

Melanoma and other cancers harbor oncogenic mutations in the protein kinase B-Raf, which leads to constitutive activation and dysregulation of MAP kinase signaling. In order to elucidate molecular determinants responsible for B-Raf control of cancer phenotypes, we present a method for phosphoprotein profiling, using negative ionization mass spectrometry to detect phosphopeptides based on their fragment ion signature caused by release of PO(3)(-). The method provides an alternative strategy for phosphoproteomics, circumventing affinity enrichment of phosphopeptides and isotopic labeling of samples. Ninety phosphorylation events were regulated by oncogenic B-Raf signaling, based on their responses to treating melanoma cells with MKK1/2 inhibitor. Regulated phosphoproteins included known signaling effectors and cytoskeletal regulators. We investigated MINERVA/FAM129B, a target belonging to a protein family with unknown category and function, and established the importance of this protein and its MAP kinase-dependent phosphorylation in controlling melanoma cell invasion into three-dimensional collagen matrix.

Dosage and temporal thresholds in microRNA proteomics.

MicroRNAs (miRNAs) modulate protein and mRNA expression through translational repression and/or mRNA decay. In this study, we combined SILAC-based proteomics and RNAseq to identify primary targets based on measurements of protein and mRNA repression and analysis of transcript 3'UTR sequences. The primary target set was used to compare different prediction algorithms, revealing higher stringency of selection by Targetscan and PITA compared with miRanda, at the expense of higher false negatives. A key finding was that significant and unexpected variations occurred in the kinetics of repression as well as the sensitivity to exogeneous miRNA concentration. Bimodal thresholds were observed, which distinguished responses to low (10 nm) versus high (50-100 nm) miRNA, as well as the onset of repression at early (12-18 h) versus late (36-48 h) times. Similar behavior was seen at the transcript level with respect to kinetics of repression. The differential thresholds were most strongly correlated with ΔΔG, the net free energy of miRNA-target interactions, which mainly reflected inverse correlations with ΔGopen, the free energy of forming 3'UTR secondary structures, at or nearby the miRNA seed matching sites. Thus, our working model is that protein binding or other competitive mechanisms variably interfere with the accessibility of miRISC to the transcript binding site. In addition, biphasic responses were observed in a subset of proteins that were partially down-regulated at early times, and further down-regulated at later times. Taken together, our findings provide evidence for varying modes of miRNA target repression, which lead to different thresholds of target responses with respect to kinetics and concentration, and predict that certain transcripts will show graded responses in sensitivity and fold-change under cellular conditions that lead to varying steady state miRNA levels.

A Phosphoproteomic Comparison of B-RAFV600E and MKK1/2 Inhibitors in Melanoma Cells.

Inhibitors of oncogenic B-RAF(V600E) and MKK1/2 have yielded remarkable responses in B-RAF(V600E)-positive melanoma patients. However, the efficacy of these inhibitors is limited by the inevitable onset of resistance. Despite the fact that these inhibitors target the same pathway, combination treatment with B-RAF(V600E) and MKK1/2 inhibitors has been shown to improve both response rates and progression-free survival in B-RAF(V600E) melanoma patients. To provide insight into the molecular nature of the combinatorial response, we used quantitative mass spectrometry to characterize the inhibitor-dependent phosphoproteome of human melanoma cells treated with the B-RAF(V600E) inhibitor PLX4032 (vemurafenib) or the MKK1/2 inhibitor AZD6244 (selumetinib). In three replicate experiments, we quantified changes at a total of 23,986 phosphosites on 4784 proteins. This included 1317 phosphosites that reproducibly decreased in response to at least one inhibitor. Phosphosites that responded to both inhibitors grouped into networks that included the nuclear pore complex, growth factor signaling, and transcriptional regulators. Although the majority of phosphosites were responsive to both inhibitors, we identified 16 sites that decreased only in response to PLX4032, suggesting rare instances where oncogenic B-RAF signaling occurs in an MKK1/2-independent manner. Only two phosphosites were identified that appeared to be uniquely responsive to AZD6244. When cells were treated with the combination of AZD6244 and PLX4032 at subsaturating concentrations (30 nm), responses at nearly all phosphosites were additive. We conclude that AZD6244 does not substantially widen the range of phosphosites inhibited by PLX4032 and that the benefit of the drug combination is best explained by their additive effects on suppressing ERK1/2 signaling. Comparison of our results to another recent ERK1/2 phosphoproteomics study revealed a surprising degree of variability in the sensitivity of phosphosites to MKK1/2 inhibitors in human cell lines, revealing unexpected cell specificity in the molecular responses to pathway activation.

N- and O-glycosylation analysis of etanercept using liquid chromatography and quadrupole time-of-flight mass spectrometry equipped with electron-transfer dissociation functionality.

Etanercept is a highly glycosylated therapeutic Fc-fusion protein that contains multiple N- and O-glycosylation sites. An in-depth characterization of the glycosylation of etanercept was carried out using liquid chromatography/mass spectrometry (LC/MS) methods in a systematic approach in which we analyzed the N- and O-linked glycans and located the occupied O-glycosylation sites. Etanercept was first treated with peptide N-glycosidase F to release the N-glycans. The N-glycan pool was labeled with a 2-aminobenzamide (2-AB) fluorescence tag and separated using ultraperformance liquid chromatography-hydrophilic interaction liquid chromatography (UPLC-HILIC). Preliminary structures were assigned using Glycobase. These assignments, which included monosaccharide sequence and linkage information, were confirmed by exoglycosidase array digestions of aliquots of the N-glycan pool. The removal of the N-glycans from etanercept facilitated the selective characterization of O-glycopeptides and enabled the O-glycans to be identified. These were predominantly of the core 1 subtype (HexHexNAc O-structure) attached to Ser/Thr residues. α2→3,6,8,9 sialidase was used to remove the sialic acid residues on the O-glycans allowing the use of an automated LC/MS(E) protocol to identify the O-glycopeptides. Electron-transfer dissociation (ETD) was then used to pinpoint the 12 occupied O-glycosylation sites. The determination of N- and O-glycans and O-glycosylation sites in etanercept provides a basis for future studies addressing the biological importance of specific protein glycosylations in the production of safe and efficacious biotherapeutics.

Characterisation of glycoproteins using a quadrupole time-of-flight mass spectrometer configured for electron transfer dissociation.

RATIONALE: Electron transfer dissociation (ETD) within ion trapping mass spectrometers has proven to be a useful tool for the characterisation of post-translational modifications. In this study, we describe the implementation of ETD upon a modified quadrupole time-of-flight (Q-ToF) system and methods for the analysis of glycoproteins. METHODS: Liquid chromatography electrospray ionisation mass spectrometry (LC/ESI-MS) was performed using a hybrid quadrupole/ion mobility/oa-ToF mass spectrometer equipped with ETD functionality. 1,4-Dicyanobenzene reagent anions necessary for the ETD reaction were generated from a glow discharge region located within the ESI source block. ETD reactions occurred in the stacked ring travelling wave ion guide (located after the quadrupole mass filter and prior to the oa-ToF mass analyser). LC/ETD was performed upon 'super-charged' tryptic glycopeptide ions produced from the recombinant monoclonal antibody trastuzumab. LC/ETD was also performed on ions from the smaller glycopeptides obtained from erythropoietin. RESULTS: ETD performed upon the quadruply 'super-charged' N-linked glycopeptide ions of trastuzumab and the triply charged O-linked glycopeptide ions of erythropoietin provided both glycosylation site assignments and full sequence information, respectively. Tandem mass (MS/MS) spectra employing collision-induced dissociation (CID) were dominated by oxonium product ions hampering full peptide sequence characterisation. CONCLUSIONS: LC/ETD on the Q-ToF system proved effective at characterising a number of different N-linked glyco-forms of the tryptic peptide, EEQYNSTYR, from trastuzumab as well as glyco-forms from the O-linked tryptic peptide, EASIPPDAASAAPLR, from erythropoietin. The data demonstrates that the glycopeptide site heterogeneity of trastuzumab and erythropoietin can be accurately characterised. In addition, the post-column mixing of the super-charging reagent, m-NBA, is an effective method to increase the precursor ion charge state and to improve ETD reaction efficiency.

Spectrum-to-spectrum searching using a proteome-wide spectral library.

The unambiguous assignment of tandem mass spectra (MS/MS) to peptide sequences remains a key unsolved problem in proteomics. Spectral library search strategies have emerged as a promising alternative for peptide identification, in which MS/MS spectra are directly compared against a reference library of confidently assigned spectra. Two problems relate to library size. First, reference spectral libraries are limited to rediscovery of previously identified peptides and are not applicable to new peptides, because of their incomplete coverage of the human proteome. Second, problems arise when searching a spectral library the size of the entire human proteome. We observed that traditional dot product scoring methods do not scale well with spectral library size, showing reduction in sensitivity when library size is increased. We show that this problem can be addressed by optimizing scoring metrics for spectrum-to-spectrum searches with large spectral libraries. MS/MS spectra for the 1.3 million predicted tryptic peptides in the human proteome are simulated using a kinetic fragmentation model (MassAnalyzer version2.1) to create a proteome-wide simulated spectral library. Searches of the simulated library increase MS/MS assignments by 24% compared with Mascot, when using probabilistic and rank based scoring methods. The proteome-wide coverage of the simulated library leads to 11% increase in unique peptide assignments, compared with parallel searches of a reference spectral library. Further improvement is attained when reference spectra and simulated spectra are combined into a hybrid spectral library, yielding 52% increased MS/MS assignments compared with Mascot searches. Our study demonstrates the advantages of using probabilistic and rank based scores to improve performance of spectrum-to-spectrum search strategies.

IsoformResolver: A peptide-centric algorithm for protein inference.

When analyzing proteins in complex samples using tandem mass spectrometry of peptides generated by proteolysis, the inference of proteins can be ambiguous, even with well-validated peptides. Unresolved questions include whether to show all possible proteins vs a minimal list, what to do when proteins are inferred ambiguously, and how to quantify peptides that bridge multiple proteins, each with distinguishing evidence. Here we describe IsoformResolver, a peptide-centric protein inference algorithm that clusters proteins in two ways, one based on peptides experimentally identified from MS/MS spectra, and the other based on peptides derived from an in silico digest of the protein database. MS/MS-derived protein groups report minimal list proteins in the context of all possible proteins, without redundantly listing peptides. In silico-derived protein groups pull together functionally related proteins, providing stable identifiers. The peptide-centric grouping strategy used by IsoformResolver allows proteins to be displayed together when they share peptides in common, providing a comprehensive yet concise way to organize protein profiles. It also summarizes information on spectral counts and is especially useful for comparing results from multiple LC-MS/MS experiments. Finally, we examine the relatedness of proteins within IsoformResolver groups and compare its performance to other protein inference software.

A simulated MS/MS library for spectrum-to-spectrum searching in large scale identification of proteins.

Identifying peptides from mass spectrometric fragmentation data (MS/MS spectra) using search strategies that map protein sequences to spectra is computationally expensive. An alternative strategy uses direct spectrum-to-spectrum matching against a reference library of previously observed MS/MS that has the advantage of evaluating matches using fragment ion intensities and other ion types than the simple set normally used. However, this approach is limited by the small sizes of the available peptide MS/MS libraries and the inability to evaluate the rate of false assignments. In this study, we observed good performance of simulated spectra generated by the kinetic model implemented in MassAnalyzer (Zhang, Z. (2004) Prediction of low-energy collision-induced dissociation spectra of peptides. Anal. Chem. 76, 3908-3922; Zhang, Z. (2005) Prediction of low-energy collision-induced dissociation spectra of peptides with three or more charges. Anal. Chem. 77, 6364-6373) as a substitute for the reference libraries used by the spectrum-to-spectrum search programs X!Hunter and BiblioSpec and similar results in comparison with the spectrum-to-sequence program Mascot. We also demonstrate the use of simulated spectra for searching against decoy sequences to estimate false discovery rates. Although we found lower score discrimination with spectrum-to-spectrum searches than with Mascot, particularly for higher charge forms, comparable peptide assignments with low false discovery rate were achieved by examining consensus between X!Hunter and Mascot, filtering results by mass accuracy, and ignoring score thresholds. Protein identification results are comparable to those achieved when evaluating consensus between Sequest and Mascot. Run times with large scale data sets using X!Hunter with the simulated spectral library are 7 times faster than Mascot and 80 times faster than Sequest with the human International Protein Index (IPI) database. We conclude that simulated spectral libraries greatly expand the search space available for spectrum-to-spectrum searching while enabling principled analyses and that the approach can be used in consensus strategies for large scale studies while reducing search times.

Quantifying the impact of chimera MS/MS spectra on peptide identification in large-scale proteomics studies.

A complicating factor for protein identification within complex mixtures by LC/MS/MS is the problem of "chimera" spectra, where two or more precursor ions with similar mass and retention time are co-sequenced by MS/MS. Chimera spectra show reduced scores due to unidentifiable fragment ions derived from contaminating parents. However, the extent of chimeras in LC/MS/MS data sets and their impact on protein identification workflows are incompletely understood. We report ChimeraCounter, a software program which detects chimeras in data sets collected on an Orbitrap/LTQ instrument. Evaluation of synthetic chimeras created from pairs of well-defined peptide MS/MS spectra reveal that chimeras reduce database search scores most significantly when contaminating fragment ion intensities exceed 20% of the targeted fragment ion intensities. In large-scale data sets, the identification rate for chimera MS/MS is 2-fold lower compared to nonchimera spectra. Importantly, this occurs in a manner which depends not on absolute precursor ion intensity, but on intensity relative to the median precursor intensity distribution. We further show that chimeras reduce the number of accepted peptide identifications by increasing false negatives while showing little increase in false positives. The results provide a framework for identifying chimeras and characterizing their contribution to the poorly understood false negative class of MS/MS.

A Case Study to Identify the Drug Conjugation Site of a Site-Specific Antibody-Drug-Conjugate Using Middle-Down Mass Spectrometry.

Middle-down mass spectrometry (MD MS) has emerged as a promising alternative to classical bottom-up approaches for protein characterization. Middle-level experiments after enzymatic digestion are routinely used for subunit analysis of monoclonal antibody (mAb)-related compounds, providing information on drug load distribution and average drug-to-antibody ratio (DAR). However, peptide mapping is still the gold standard for primary amino acid sequence assessment, post-translational modifications (PTM), and drug conjugation identification and localization. However, peptide mapping strategies can be challenging when dealing with more complex and heterogeneous mAb formats, like antibody-drug conjugates (ADCs). We report here, for the first time, MD MS analysis of a third-generation site-specific DAR4 ADC using different fragmentation techniques, including higher-energy collisional- (HCD), electron-transfer (ETD) dissociation and 213 nm ultraviolet photodissociation (UVPD). UVPD used as a standalone technique for ADC subunit analysis afforded, within the same liquid chromatography-MS/MS run, enhanced performance in terms of primary sequence coverage compared to HCD- or ETD-based MD approaches, and generated substantially more MS/MS fragments containing either drug conjugation or glycosylation site information, leading to confident drug/glycosylation site identification. In addition, our results highlight the complementarity of ETD and UVPD for both primary sequence validation and drug conjugation/glycosylation site assessment. Altogether, our results highlight the potential of UVPD for ADC MD MS analysis for drug conjugation/glycosylation site assessment, and indicate that MD MS strategies can improve structural characterization of empowered next-generation mAb-based formats, especially for PTMs and drug conjugation sites validation.

Large-Scale Examination of Factors Influencing Phosphopeptide Neutral Loss during Collision Induced Dissociation.

Collision-induced dissociation (CID) remains the predominant mass spectrometry-based method for identifying phosphorylation sites in complex mixtures. Unfortunately, the gas-phase reactivity of phosphoester bonds results in MS/MS spectra dominated by phosphoric acid (H3PO4) neutral loss events, suppressing informative peptide backbone cleavages. To understand the major drivers of H3PO4 neutral loss, we performed robust nonparametric statistical analysis of local and distal sequence effects on the magnitude and variability of neutral loss, using a collection of over 35,000 unique phosphopeptide MS/MS spectra. In contrast to peptide amide dissociation pathways, which are strongly influenced by adjacent amino acid side chains, we find that neutral loss of H3PO4 is affected by both proximal and distal sites, most notably basic residues and the peptide N-terminal primary amine. Previous studies have suggested that protonated basic residues catalyze neutral loss through direct interactions with the phosphate. In contrast, we find that nearby basic groups decrease neutral loss regardless of mobility class, an effect only seen by stratifying spectra by charge-mobility. The most inhibitory bases are those immediately N-terminal to the phosphate, presumably because of steric hindrances in catalyzing neutral loss. Further evidence of steric effects is shown by the presence of proline, which can dramatically reduce the presence of neutral loss when between the phosphate and a possible charge donor. In mobile proton spectra, the N-terminus is the strongest predictor of high neutral loss, with proximity to the N-terminus essential for peptides to exhibit the highest levels of neutral loss.

Identification of a family of fatty-acid-speciated sonic hedgehog proteins, whose members display differential biological properties.

Hedgehog (HH) proteins are proteolytically processed into a biologically active form that is covalently modified by cholesterol and palmitate. However, most studies of HH biogenesis have characterized protein from cells in which HH is overexpressed. We purified Sonic Hedgehog (SHH) from cells expressing physiologically relevant levels and showed that it was more potent than SHH isolated from overexpressing cells. Furthermore, the SHH in our preparations was modified with a diverse spectrum of fatty acids on its amino termini, and this spectrum of fatty acids varied dramatically depending on the growth conditions of the cells. The fatty acid composition of SHH affected its trafficking to lipid rafts as well as its potency. Our results suggest that HH proteins exist as a family of diverse lipid-speciated proteins that might be altered in different physiological and pathological contexts in order to regulate distinct properties of HH proteins.

The influence of adnectin binding on the extracellular domain of epidermal growth factor receptor.

The precise and unambiguous elucidation and characterization of interactions between a high affinity recognition entity and its cognate protein provides important insights for the design and development of drugs with optimized properties and efficacy. In oncology, one important target protein has been shown to be the epidermal growth factor receptor (EGFR) through the development of therapeutic anticancer antibodies that are selective inhibitors of EGFR activity. More recently, smaller protein derived from the 10th type III domain of human fibronectin termed an adnectin has also been shown to inhibit EGFR in clinical studies. The mechanism of EGFR inhibition by either an adnectin or an antibody results from specific binding of the high affinity protein to the extracellular portion of EGFR (exEGFR) in a manner that prevents phosphorylation of the intracellular kinase domain of the receptor and thereby blocks intracellular signaling. Here, the structural changes induced upon binding were studied by probing the solution conformations of full length exEGFR alone and bound to a cognate adnectin through hydrogen/deuterium exchange mass spectrometry (HDX MS). The effects of binding in solution were identified and compared with the structure of a bound complex determined by X-ray crystallography.ᅟ

The non-coding B2 RNA binds to the DNA cleft and active-site region of RNA polymerase II.

The B2 family of short interspersed elements is transcribed into non-coding RNA by RNA polymerase III. The ~180-nt B2 RNA has been shown to potently repress mRNA transcription by binding tightly to RNA polymerase II (Pol II) and assembling with it into complexes on promoter DNA, where it keeps the polymerase from properly engaging the promoter DNA. Mammalian Pol II is an ~500-kDa complex that contains 12 different protein subunits, providing many possible surfaces for interaction with B2 RNA. We found that the carboxy-terminal domain of the largest Pol II subunit was not required for B2 RNA to bind Pol II and repress transcription in vitro. To identify the surface on Pol II to which the minimal functional region of B2 RNA binds, we coupled multi-step affinity purification, reversible formaldehyde cross-linking, peptide sequencing by mass spectrometry, and analysis of peptide enrichment. The Pol II peptides most highly recovered after cross-linking to B2 RNA mapped to the DNA binding cleft and active-site region of Pol II. These studies determine the location of a defined nucleic acid binding site on a large, native, multi-subunit complex and provide insight into the mechanism of transcriptional repression by B2 RNA.

Wnt5a directs polarized calcium gradients by recruiting cortical endoplasmic reticulum to the cell trailing edge.

Wnt5a directs the assembly of the Wnt-receptor-actin-myosin-polarity (WRAMP) structure, which integrates cell-adhesion receptors with F-actin and myosin to form a microfilament array associated with multivesicular bodies (MVBs). The WRAMP structure is polarized to the cell posterior, where it directs tail-end membrane retraction, driving forward translocation of the cell body. Here we define constituents of the WRAMP proteome, including regulators of microfilament and microtubule dynamics, protein interactions, and enzymatic activity. IQGAP1, a scaffold for F-actin nucleation and crosslinking, is necessary for WRAMP structure formation, potentially bridging microfilaments and MVBs. Vesicle coat proteins, including coatomer-I subunits, localize to and are required for the WRAMP structure. Electron microscopy and live imaging demonstrate movement of the ER to the WRAMP structure and plasma membrane, followed by elevation of intracellular Ca2+. Thus, Wnt5a controls directional movement by recruiting cortical ER to mobilize a rear-directed, localized Ca2+ signal, activating actomyosin contraction and adhesion disassembly for membrane retraction.

Elemental and molecular evidence of soot- and char-derived black carbon inputs to New York City's atmosphere during the 20th century.

Soot black carbon (here expressed as GBC) is present in sediments of Central Park and Prospect Park Lakes, New York City (NYC), and peaks in the middle of the 20th Century at the highest values (1-3% dry weight) ever reported in urban lakes. During that period (approximately 1940-1970), the GBC represents up to 28% of the total organic carbon (OC). Radionuclide-normalized whole core inventories of accumulated GBC are similar in the two lakes which are separated by approximately 15 km, suggesting that emissions of fine soot particles may have accumulated homogeneously over at least the urban center of NYC. The distribution of polycyclic aromatic hydrocarbons (PAHs) in the sediments is decoupled from that of GBC. The highest levels of total PAHs correspond to peak coal use for space heating in NYC in the early 1900s. In contrast, GBC concentrations were highest in the mid 1900s, a period when oil combustion dominated local fossil fuel use and incineration of municipal solid waste (MSW) was common practice in NYC. Decreases in GBC levels observed in more recently deposited sediments are consistent with improvements in particle emissions control systems. Non-soot BC (char) was identified by a high carbon to nitrogen (C/N) ratio that persisted after correction for GBC. This likely tracer of MSW incineration was estimated to contribute an additional '35% of total organic carbon found in the sediments deposited during the peak period of combustion. The temporal trends of soot-BC observed in our lake cores do not agree with published historical reconstructions based on fuel consumption and estimated emission factors.