Techniques for estimating genetically variable peptides and semi-continuous likelihoods from massively parallel sequencing data.

Forensic genetic investigations typically rely on analysis of DNA for attribution purposes. There are times, however, when the amount and/or the quality of the DNA is limited, and thus little or no information can be obtained regarding the source of the sample. An alternative biochemical target that also contains genetic signatures is protein. One class of genetic signatures is protein polymorphisms that are a direct consequence of simple/single/short nucleotide polymorphisms (SNPs) in DNA. However, to interpret protein polymorphisms in a forensic context, certain complexities must be understood and addressed. These complexities include: 1) SNPs can generate 0, 1, or arbitrarily many polymorphisms in a polypeptide; and 2) as an object of expression that is modulated by alleles, genes and interactions with the environment, proteins may be present or absent in a given sample. To address these issues, a novel approach was taken to generate the expected protein alleles in a reference sample based on whole genome (or exome) sequence data and assess the significance of the evidence using a haplotype-based semi-continuous likelihood algorithm that leverages whole proteome data. Converting the genomic information into the proteomic information allows for the zero-to-many relationship between SNPs and GVPs to be abstracted away. When viewed as a haplotype, many GVPs that correspond to the same SNP is equivalent to many SNPs in perfect linkage disequilibrium (LD). As long as the likelihood formulation correctly accounts for LD, the correspondence between the SNP and the proteome can be safely neglected. Tests were performed on simulated samples, including single-source and two-person mixtures, and the power of using a classical semi-continuous likelihood versus one that has been adapted to neglect drop-out was compared. Additionally, summary statistics and a rudimentary set of decision guidelines were introduced to help identify mixtures from protein data.

Optimization of proteomics sample preparation for forensic analysis of skin samples.

We present an efficient protein extraction and in-solution enzymatic digestion protocol optimized for mass spectrometry-based proteomics studies of human skin samples. Human skin cells are a proteinaceous matrix that can enable forensic identification of individuals. We performed a systematic optimization of proteomic sample preparation for a protein-based human forensic identification application. Digestion parameters, including incubation duration, temperature, and the type and concentration of surfactant, were systematically varied to maximize digestion completeness. Through replicate digestions, parameter optimization was performed to maximize repeatability and increase the number of identified peptides and proteins. Final digestion conditions were selected based on the parameters that yielded the greatest percent of peptides with zero missed tryptic cleavages, which benefit the analysis of genetically variable peptides (GVPs). We evaluated the final digestion conditions for identification of GVPs by applying MS-based proteomics on a mixed-donor sample. The results were searched against a human proteome database appended with a database of GVPs constructed from known non-synonymous single nucleotide polymorphisms (SNPs) that occur at known population frequencies. The aim of this study was to demonstrate the potential of our proteomics sample preparation for future implementation of GVP analysis by forensic laboratories to facilitate human identification. SIGNIFICANCE: Genetically variable peptides (GVPs) can provide forensic evidence that is complementary to traditional DNA profiling and be potentially used for human identification. An efficient protein extraction and reproducible digestion method of skin proteins is a key contributor for downstream analysis of GVPs and further development of this technology in forensic application. In this study, we optimized the enzymatic digestion conditions, such as incubation time and temperature, for skin samples. Our study is among the first attempts towards optimization of proteomics sample preparation for protein-based skin identification in forensic applications such as touch samples. Our digestion method employs RapiGest (an acid-labile surfactant), trypsin enzymatic digestion, and an incubation time of 16 h at 37 °C.

Trace explosives sampling for security applications (TESSA) study: Evaluation of procedures and methodology for contact sampling efficiency.

The detection of trace amounts of explosive materials is critical to the security at mass transit centers (e.g., airports and railway stations). In a typical screening process, a trap is used to probe a surface of interest to collect and transfer particulate residue to a detector for analysis. The collection of residues from the surface being probed is widely viewed as the limiting step in this process. A multi-institutional study was performed to establish a methodology for the evaluation of sampling media collection efficiencies. Dry deposited residues of 1,3,5-trinitroperhydro-1,3,5-triazine (RDX), C-4 (an RDX-based explosive), and pentaerythritol tetranitrate (PETN) were harvested from acrylonitrile butadiene styrene (ABS) plastic, ballistic nylon (NYL), and uncoated aluminum surfaces using muslin, Texwipe cotton, and stainless-steel mesh traps. Transfer and collection efficiencies of the sample media were calculated based on liquid chromatography-mass spectrometry analysis. Dry transfer efficiencies (DTE%) to all tested surfaces were greater than 75%, with transfer to ABS plastic being the lowest. Collection efficiency (CE%) varied significantly across the traps and the surfaces, yet some conclusions can be drawn; nylon had the lowest CE% for all cases (∼10%), and the stainless steel mesh had the lowest CE% for the evaluated traps (∼20%). Though the testing parameters have been standardized among the participants to establish a framework for an independent comparison of contact sampling media and surfaces, substantial variations in the DTE% and the CE% were observed, suggesting that other variables can affect contact sampling.

Fractionation of DNA and protein from individual latent fingerprints for forensic analysis.

Human touch samples represent a significant portion of forensic DNA casework. Yet, the generally low abundance of genetic material combined with the predominantly extracellular nature of DNA in these samples makes DNA-based forensic analysis exceptionally challenging. Human proteins present in these same touch samples offer an abundant and environmentally-robust alternative. Proteogenomic methods, using protein sequence variants arising from nonsynonymous DNA mutations, have recently been applied to forensic analysis and may represent a viable option looking forward. However, DNA analysis remains the gold standard and any proteomics-based methods would need to consider how DNA could be co-extracted from samples without significant loss. Herein, we describe a simple workflow for the collection, enrichment and fractionation of DNA and protein in latent fingerprint samples. This approach ensures that DNA collected from a latent fingerprint can be analyzed by traditional DNA casework methods, while protein can be proteolytically digested and analyzed via standard liquid chromatography-tandem mass spectrometry-based proteomics methods from the same touch sample. Sample collection from non-porous surfaces (i.e., glass) is performed through the application of an anionic surfactant over the fingermark. The sample is then split into separate DNA and protein fractions following centrifugation to enrich the protein fraction by pelleting skin cells. The results indicate that this workflow permits analysis of DNA within the sample, yet highlights the challenge posed by the trace nature of DNA in touch samples and the potential for DNA to degrade over time. Protein deposited in touch samples does not appear to share this limitation, with robust protein quantities collected across multiple human donors. The quantity and quality of protein remains robust regardless of fingerprint age. The proteomic content of these samples is consistent across individual donors and fingerprint age, supporting the future application of genetically variable peptide (GVP) analysis of touch samples for forensic identification.

An algorithm for random match probability calculation from peptide sequences.

For the past three decades, forensic genetic investigations have focused on elucidating DNA signatures. While DNA has a number of desirable properties (e.g., presence in most biological materials, an amenable chemistry for analysis and well-developed statistics), DNA also has limitations. DNA may be in low quantity in some tissues, such as hair, and in some tissues it may degrade more readily than its protein counterparts. Recent research efforts have shown the feasibility of performing protein-based human identification in cases in which recovery of DNA is challenged; however, the methods involved in assessing the rarity of a given protein profile have not been addressed adequately. In this paper an algorithm is proposed that describes the computation of a random match probability (RMP) resulting from a genetically variable peptide signature. The approach described herein explicitly models proteomic error and genetic linkage, makes no assumptions as to allelic drop-out, and maps the observed proteomic alleles to their expected protein products from DNA which, in turn, permits standard corrections for population structure and finite database sizes. To assess the feasibility of this approach, RMPs were estimated from peptide profiles of skin samples from 25 individuals of European ancestry. 126 common peptide alleles were used in this approach, yielding a mean RMP of approximately 10-2.

Artificial fingerprints for cross-comparison of forensic DNA and protein recovery methods.

Quantitative genomic and proteomic evaluation of human latent fingerprint depositions represents a challenge within the forensic field, due to the high variability in the amount of DNA and protein initially deposited. To better assess recovery techniques for touch depositions, we present a method to produce simple and customizable artificial fingerprints. These artificial fingerprint samples include the primary components of a typical latent fingerprint, specifically sebaceous fluid, eccrine perspiration, extracellular DNA, and proteinaceous epidermal skin material (i.e., shed skin cells). A commercially available emulsion of sebaceous and eccrine perspiration material provides a chemically-relevant suspension solution for fingerprint deposition, simplifying artificial fingerprint production. Extracted human genomic DNA is added to accurately mimic the extracellular DNA content of a typical latent print and comparable DNA yields are recovered from the artificial prints relative to human prints across surface types. Capitalizing on recent advancements in the use of protein sequence identification for human forensic analysis, these samples also contain a representative quantity of protein, originating from epidermal skin cells collected from the fingers and palms of volunteers. Proteomic sequencing by liquid chromatography-tandem mass spectrometry (LC-MS/MS) analysis indicates a high level of protein overlap between artificial and latent prints. Data are available via ProteomeXchange with identifier PXD015445. By including known quantities of DNA and protein into each artificial print, this method enables total DNA and protein recovery to be quantitatively assessed across different sample collection and extraction methods to better evaluate extraction efficiency. Collectively, these artificial fingerprint samples are simple to make, highly versatile and customizable, and accurately represent the biochemical composition and biological signatures of human fingerprints.

Preferential cleavage of N-N hydrazone bonds for sequencing bis-arylhydrazone conjugated peptides by electron transfer dissociation.

Electron transfer dissociation (ETD) was used to sequence bis-arylhydrazone (BAH)-cross-linked peptides through preferential cleavage of the hydrazone bond. On average, 58% of the observed ETD product ion abundance was accounted for by fragment ions due to selective cleavage of the N12-N13 hydrazone bond. Dissociation of the N12-N13 hydrazone bond yielded the two constituent peptides, one an even-electron product ion termed Lalpha12, the other an odd-electron radical ion termed Lbeta11(*), which allowed each peptide to be individually sequenced by MS/MS methods and the site of cross-linking to be identified. The proposed pathway for the dissociation of the hydrazone bond involves transfer of the electron directly to the protonated hydrazone functionality and subsequent rearrangement to yield the Lalpha12 and Lbeta11(*) products. Collision induced dissociation (CID) of the even-electron Lalpha12 product yielded a series of b- and y-type ions; CID of the odd-electron Lbeta11(*) product resulted in a wide range of fragment ions including a-, b-, c-, y-, and z-type ions.

Infrared multiphoton dissociation of small-interfering RNA anions and cations.

Infrared multiphoton dissociation (IRMPD) on a linear ion trap mass spectrometer is applied for the sequencing of small interfering RNA (siRNA). Both single-strand siRNAs and duplex siRNA were characterized by IRMPD, and the results were compared with that obtained by traditional ion trap-based collision induced dissociation (CID). The single-strand siRNA anions were observed to dissociate via cleavage of the 5' P-O bonds yielding c- and y-type product ions as well as undergo neutral base loss. Full sequence coverage of the siRNA anions was obtained by both IRMPD and CID. While the CID mass spectra were dominated by base loss ions, accounting for approximately 25% to 40% of the product ion current, these ions were eliminated through secondary dissociation by increasing the irradiation time in the IRMPD mass spectra to produce higher abundances of informative sequence ions. With longer irradiation times, however, internal ions corresponding to cleavage of two 5' P-O bonds began to populate the product ion mass spectra as well as higher abundances of [a - Base] and w-type ions. IRMPD of siRNA cations predominantly produced c- and y-type ions with minimal contributions of [a - Base] and w-type ions to the product ion current; the presence of only two complementary series of product ions in the IRMPD mass spectra simplified spectral interpretation. In addition, IRMPD produced high abundances of protonated nucleobases, [G + H](+), [A + H](+), and [C + H](+), which were not detected in the CID mass spectra due to the low-mass cut-off associated with conventional CID in ion traps. CID and IRMPD using short irradiation times of duplex siRNA resulted in strand separation, similar to the dissociation trends observed for duplex DNA. With longer irradiation times, however, the individual single-strands underwent secondary dissociation to yield informative sequence ions not obtained by CID.

Top-down protein fragmentation by infrared multiphoton dissociation in a dual pressure linear ion trap.

Infrared multiphoton dissociation (IRMPD) was implemented in a novel dual pressure linear ion trap for rapid top-down proteomics. The high pressure cell provided improved trapping and isolation efficiencies while the isotopic profiles of 10+ charged ions could be resolved by mass analysis in the low pressure cell that enabled effective top down protein identification. Striking differences between IRMPD in the low pressure cell and CID in the high pressure cell were observed for proteins ranging from 8.6 to 29 kDa. Because of secondary dissociation, IRMPD yielded product ions in significantly lower charge states as compared to CID, thus facilitating more accurate mass identification and streamlining product ion assignment. This outcome was especially useful for database searching of larger proteins (approximately 29 kDa) as IRMPD substantially improved protein identification and scoring confidence. Also, IRMPD showed an increased selectivity toward backbone cleavages N-terminal to proline and C-terminal to acidic residues (especially for the lowest charge states), which could be useful for a priori spectral predictions and enhanced database searching for protein identification.

Reduction of chemical noise in electrospray ionization mass spectrometry by supplemental IR activation.

Supplemental infrared (IR) activation was applied to reduce background chemical noise and increase analyte ion signal in a linear ion trap mass spectrometer. Peptides, proteins, and small molecules were all introduced by electrospray ionization, and when regions of chemical noise were isolated and subjected to IR irradiation, protonated analyte molecules were observed in the product ion mass spectra. By isolating the entire mass range (e.g., m/z 400-2000) and then irradiating all ions in the trap, supplemental IR activation increased the signal of singly protonated peptides by almost 70% and by 40%-55% for the lower charge states of cytochrome c. This increase in analyte ion signal was less dramatic for the higher charge states of peptides and proteins. The chemical noise present in the mass spectra is attributed to incomplete desolvation of the electrospray, as the abundance of the protonated peptides observed upon supplemental IR activation of the chemical noise decreased with higher inlet capillary temperatures. Collision activation was not as effective for desolvating the ions present in the chemical noise.

Infrared multiphoton dissociation of peptide cations in a dual pressure linear ion trap mass spectrometer.

A dual pressure linear ion trap mass spectrometer was modified to permit infrared multiphoton dissociation (IRMPD) in each of the two cells-the first a high pressure cell operated at nominally 5 x 10(-3) Torr and the second a low pressure cell operated at nominally 3 x 10(-4) Torr. When IRMPD was performed in the high pressure cell, most peptide ions did not undergo significant photodissociation; however, in the low pressure cell peptide cations were efficiently dissociated with less than 25 ms of IR irradiation regardless of charge state. IRMPD of peptide cations allowed the detection of low m/z product ions including the y(1) fragments and immonium ions which are not typically observed by ion trap collision induced dissociation (CID). Photodissociation efficiencies of approximately 100% and MS/MS (tandem mass spectrometry) efficiencies of greater than 60% were observed for both multiply and singly protonated peptides. In general, higher sequence coverage of peptides was obtained using IRMPD over CID. Further, greater than 90% of the product ion current in the IRMPD mass spectra of doubly charged peptide ions was composed of singly charged product ions compared to the CID mass spectra in which the abundances of the multiply and singly charged product ions were equally divided. Highly charged primary product ions also underwent efficient photodissociation to yield singly charged secondary product ions, thus simplifying the IRMPD product ion mass spectra.

Ultraviolet photodissociation mass spectrometry of bis-aryl hydrazone conjugated peptides.

Ultraviolet photodissociation (UVPD) at 355 nm was used to rapidly identify peptides which had been chemically conjugated through bis-aryl hydrazone (BAH) moieties. The two biomolecules of interest were separately tagged to introduce either an aldehyde or a hydrazine and then conjugated together through these functional groups to from the UV-chromogenic BAH-group. In a mock mixture of peptides, UVPD was used to screen for the BAH-conjugated peptides in direct infusion ESI-UVPD-MS and online LC-UVPD-MS methods by comparing the abundances of the ions with the laser off and with the laser on. Only the BAH-conjugated peptides were observed to photodissociate upon exposure to UV irradiation, thus affording excellent selectivity for the pinpointing the relevant conjugated peptides in a complex mixture of nonconjugated peptides. UVPD analysis of conjugated model peptides indicated that the UVPD efficiencies of these species were charge state dependent. BAH-conjugated peptides that had a mobile proton which could protonate the basic BAH-moiety underwent more efficient photodissociation than the peptide ions with sequestered protons. Ultraviolet photodissociation of BAH-cross-linked peptides also yielded more diagnostic sequence ions than CID to unambiguously locate the site of conjugation.

Chromogenic cross-linker for the characterization of protein structure by infrared multiphoton dissociation mass spectrometry.

We have developed a new IR chromogenic cross-linker (IRCX) to aid in rapidly distinguishing cross-linked peptides from unmodified species in complex mixtures. By incorporating a phosphate functional group into the cross-linker, one can take advantage of its unique IR absorption properties, affording selective infrared multiphoton dissociation (IRMPD) of the cross-linked peptides. In a mock mixture of unmodified peptides and IRCX-cross-linked peptides (intramolecularly and intermolecularly cross-linked), only the peptides containing the IRCX modification were shown to dissociate upon exposure to 50 ms of 10.6-microm radiation. LC-IRMPD-MS proved to be an effective method to distinguish the cross-linked peptides in a tryptic digest of IRCX-cross-linked ubiquitin. A total of four intermolecular cross-links and two dead-end modifications were identified using IRCX and LC-IRMPD-MS. IRMPD of these cross-linked peptides resulted in secondary dissociation of all primary fragment ions containing the chromophore, producing a series of unmodified b- or y-type ions that allowed the cross-linked peptides to be sequenced without the need for collision-induced dissociation.

Impact of proline and aspartic acid residues on the dissociation of intermolecularly crosslinked peptides.

The dissociation of intermolecularly crosslinked peptides was evaluated for a series of peptides with proline or aspartic acid residues positioned adjacent to the crosslinking sites (lysine residues). The peptides were crosslinked with either disuccinimidyl suberate (DSS) or disuccinimidyl L-tartrate (DST), and the influence of proline and aspartic acid residues on the fragmentation patterns were investigated for precursor ions with and without a mobile proton. Collisionally activated dissociation (CAD) spectra of aspartic acid-containing crosslinked peptide ions, doubly-charged with both protons sequestered, were dominated by cleavage C-terminal to the Asp residue, similar to that of unmodified peptides. The proline-containing crosslinked peptides exhibited a high degree of internal ion formation, with the resulting product ions having an N-terminal proline residue. Upon dissociation of the doubly-charged crosslinked peptides, twenty to fifty percent of the fragment ion abundance was accounted for by multiple cleavage products. Crosslinked peptides possessing a mobile proton yielded almost a full series of b- and y-type fragment ions, with only proline-directed fragments still observed at high abundances. Interestingly, the crosslinked peptides exhibited a tendency to dissociate at the amide bond C-terminal to the crosslinked lysine residue, relative to the N-terminal side. One could envision updating computer algorithms to include these crosslinker specific product ions--particularly for precursor ions with localized protons--that provide complementary and confirmatory information, to offer more confident identification of both the crosslinked peptides and the location of the crosslink, as well as affording predictive guidelines for interpretation of the product-ion spectra of crosslinked peptides.

Protocol for the thermodynamic analysis of some proteins using an H/D exchange- and mass spectrometry-based technique.

SUPREX (stability of unpurified proteins from rates of H/D exchange) is a new H/D exchange- and mass spectrometry-based technique for the measurement of protein folding free energies (i.e., DeltaG values) and protein folding m values (i.e., deltaDeltaG/delta[denaturant]). Robust protocols for the acquisition and analysis of SUPREX data have been established and shown to be useful for the analysis of a number of different protein systems. Here we report on the SUPREX behavior of a special class of proteins that are not amenable to conventional SUPREX analyses using previously established protocols. This class of proteins includes protein systems that require an extended time to reach a folding/unfolding equilibrium in chemical denaturant-induced equilibrium unfolding experiments. As part of this work we use ubiquitin as a model system to highlight the complications that can arise in the conventional SUPREX analysis of such protein systems, and we describe a modified SUPREX protocol that can be used to eliminate these complications.