|
|
||||||||
Molecular Diagnostics and Genetics |
1 W.M. Keck FT-ICR Mass Spectrometry Laboratory, Mayo Proteomics Research Center, Rochester, MN.
Departments of2
Biochemistry and Molecular Biology and3
Hematology, Mayo Clinic College of Medicine, Rochester, MN.
4 Department of Chemistry, Virginia Commonwealth University, Richmond, VA.
aAddress correspondence to this author at: Medical Sciences Building 3-115, Mayo Clinic College of Medicine, 200 First Street SW, Rochester, MN 55905. Fax 507-284-9261; e-mail muddiman.david{at}mayo.edu.
| Abstract |
|---|
|
|
|---|
Methods: TTR was immunoprecipitated from serum by use of a polyclonal antibody and subsequently reduced with tris(2-carboxyethyl)phosphine. The purified TTR was then analyzed by fast-gradient liquid chromatographydual-electrospray ionization Fourier-transform ion-cyclotron-resonance (FT-ICR) mass spectrometry. DNA sequencing was performed on all samples used in this study.
Results: Because of the inherent limitations in achieving high mass measurement accuracy based on the most abundant isotopic mass, we applied a fitting procedure that allowed determination of monoisotopic mass. Wild-type TTR (mean molecular mass, 13 761 Da) and its associated variant forms could be distinguished because of the high molecular mass accuracy afforded by FT-ICR (
3 ppm) except for instances involving isobaric species or when isotopic distributions overlapped significantly. The [M + 11 H+]11+ charge state for all samples was used to determine the mass accuracies for both wild-type and variant forms of the protein. We correctly assigned seven of seven TTR variants. Moreover, using a combination of proteomic and genomic technologies, we discovered and characterized a previously unreported cis double mutation with a mass only 2 Da different from wild-type TTR. Furthermore, DNA sequencing of the TTR gene for all individuals in this study completely agreed with the intact protein measurements.
Conclusions: FT-ICR mass spectrometry has sufficient mass accuracy to identify genetic variants of immunoaffinity-purified TTR. We believe that 91% of known TTR variants could be detected by this technique.
| Introduction |
|---|
|
|
|---|
Mass spectrometry (MS) has played a crucial role in proteomics because of the soft ionization techniques of electrospray ionization (ESI)(8)(9) and matrix-assisted laser desorption/ionization(10)(11). Wada et al.(12) were the first to perform the analysis of tryptic fragments of TTR by MS. This spawned several other mass spectrometric approaches to analyze and catalog TTR for any variations that differ from the wild-type amino acid sequence(6)(13)(14)(15)(16)(17)(18)(19)(20).
Fourier-transform ion-cyclotron-resonance (FT-ICR) MS can achieve mass measurement accuracy [MMA; (observed mass theoretical mass)/theoretical mass x 106] of
3 ppm if it can accurately measure the cyclotron frequency of the ions of interest. Changes in ion cyclotron frequency occur as a result of variations in magnetic field strength, ion population (space charge), trapping potentials, and excitation conditions(21)(22)(23)(24). Coulombic interactions between the ions themselves are known to cause systematic shifts in the observed cyclotron frequency, otherwise known as space-charge effects. There are two categories of space-charge effects that influence the observed cyclotron frequency. The more pronounced of the two is referred to as "global" space-charge effect and is commonly taken into account in FT-ICR calibration laws(23)(24)(25). Systematic error also originates from "local" space-charge effects attributable to more subtle "indirect" coulombic mechanisms, as observed by Masselson et al.(26) and Flora et al.(27). To achieve high MMA, these space-charge effects must be quantitatively, routinely, and confidently accounted for. Internal calibration effectively compensates for space-charge effects. Because the internal mass calibrant (IMC) and analyte experience identical conditions within the ICR, the global space-charge component can be accounted for in the calibration equation(23)(24).
Hannis and Muddiman(22) were the first to report success in using a dual ESI source for FT-ICR MS, routinely detecting intact biological macromolecules at mass accuracies of 3 ppm. Flora et al.(27) demonstrated MMAs of
3 ppm (99% confidence interval of the mean) with the use of a dual ESI source for internally calibrating product ions formed from a 15mer oligonucleotide by the use of sustained off-resonance irradiation collision dissociation at 4.7 Tesla. Our group has shown the applicability of a dual ESI source(22)(27)(28)(29)(30)(31)(32)(33)(34), has surmounted the limitations of the first-generation source(35), and has successfully applied this second-generation source(34)(35)(36). Other groups have recently adopted this general strategy to achieve high MMA in FT-ICR(37)(38)(39).
Here we describe a "targeted" proteomics approach for the characterization and identification of TTR variants based on immunoprecipitation by rabbit anti-human TTR antibodies and accurate mass measurements. Both wild-type and variant TTR are isolated from the vast number of proteins in serum by immunoprecipitation with a polyclonal antibody. TTR variants are detected by accurate mass measurements using a dual ESI source for internal calibration.
| Materials and Methods |
|---|
|
|
|---|
with a Barnstead Nanopure Infinity ultrapure water system.
preparation of ttr protein extraction from serum
The Institutional Review Board of the Mayo Clinic granted approval for the study, and peripheral blood samples were drawn after informed consent was obtained from each participant.
TTR was purified according to the method of Bergen et al.(40) with slight modifications as described here. Rabbit anti-human TTR antibodies (Dako Corp.) were coupled to POROS AL (PerSeptive Biosystems) medium by the "recycling method" according to the manufacturers instructions. We mixed 50 µL of serum and 50 µL of phosphate-buffered saline (10 mmol/L phosphate, pH 7.4, 150 mmol/L NaCl) before adding 20 µL of resin suspension (anti-human TTR antibodies) in a HandeeTM MiniSpin column (Pierce Biotechnology). The solution was then incubated for
1.5 h to allow for binding of the TTR to the antibody resin. The column was then centrifuged and the resin washed twice with phosphate-buffered saline. The purified TTR protein was eluted from the resin by 20 µL of 0.1 mol/L glycine (pH 2.8) and collected in a microcentrifuge tube.
We reduced the TTR samples with TCEP, adding 2 µL of a 100 mmol/L solution to the purified sample (20 µL), which we then kept at room temperature for 15 min.
one-dimensional polyacrylamide gel electrophoresis
Immunoaffinity-purified TTR and diluted plasma were diluted 1:1 with Laemmli sample buffer [63 mmol/L Tris, 100 mL/L glycerol, 20 g/L sodium dodecyl sulfate (SDS), and 0.0025 g/L bromphenol blue; 2-mercaptoethanol was added to this solution to a final concentration of 50 mL/L]. The immunoaffinity-purified TTR was acidic; we therefore added 1 µL of 3 mol/L Tris (pH 8.3) to each sample before boiling all samples for 3 min and then applying them to a 15% Tris-HCl Criterion precast gel (Bio-Rad). Tris/glycine running buffer (25 mmol/L Tris, 192 mmol/L glycine, 1 g/L SDS, pH 8.3) was used to run the gel for 0.9 h at 200 V. The gels were washed three times with water and then stained with Bio-Safe Coomassie G250 stain (Bio-Rad). After destaining with water, the gels were scanned with a GS800 densitometer (Bio-Rad). Protein molecular mass markers were included for reference.
ms
Mass spectra were obtained with a modified ESI-FT-ICR mass spectrometer (IonSpec) equipped with an actively shielded 7-Tesla superconducting magnet (Cryomagnetics). The ESI source (Analytica, Branford) was modified to accept a heated metal capillary(41) and a new dual ESI source(35) and has been described elsewhere in detail(34)(35). Software for calculation of monoisotopic mass using an algorithm described by Senko et al.(42) was written in C++ in Linux. Theoretical isotopic distributions were generated with the Mercury algorithm(43).
NANO-LC-MS of ttr
The chromatography system consisted of two Shimadzu LC-10ADvp pumps with an SCL-10Avp controller (flow rate was adjusted by splitting to give
500 nL/min through the column). Solvent A was wateracetonitrile (95:5 by volume) containing 2 g/L formic acid; solvent B was wateracetonitrile (5:95 by volume) containing 2 g/L formic acid. A third LC-10ADvp pump [methanolwater (10:90 by volume) containing 0.1 g/L trifluoroacetic acid; 30 µL/min] was used to load the TTR sample on a C4 trap. The trap consisted of a LC Packings micro precolumn (0.3 x 5 mm) packed with Magic C4 (300 Å pore size, 5-µm bead diameter; Michrom Bioresources, Inc.).
The C4 trap was plumbed into a 10-port valve (Cheminert; VICI Valco Instruments). In position A the sample was injected and loaded on the C4 trap with the sample-loading pump. In position B, the trap was in-line with the nano-LC column and the analyte was eluted off the trap with an organic gradient. The gradient profile for eluting TTR from the C4 trap was 03 min held at 5% B, 36 min linearly ramped to 60% B, 69 min held at 60% B, and 910 min linearly ramped to 95% B.
Data acquisition began simultaneously with the start of the LC gradient. The total hexapole accumulation time for both the internal standard and apomyoglobin digest was 2000 ms; 250 ms were spent on the IMC emitter, and the remaining 1750 ms were spent on the eluate from the nano-LC column. The IMC consisted of a 2 µmol/L solution of angiotensin I, bradykinin, and/or glucagon, which was infused directly through a second emitter at a flow rate of 300 nL/min while the LC effluent was sprayed through the primary emitter during the LC-MS run. Each mass spectrum was internally calibrated by use of the monoisotopic peaks that directly flanked the [M + 11H+]11+ charge state of the wild-type and variant proteins.
gene sequencing
High-quality genomic DNA was isolated from peripheral blood buffy coat by use of DNAzol (Invitrogen). Further details can be found elsewhere(40).
| Results and Discussion |
|---|
|
|
|---|
The purity of the immunoaffinity-purified TTR relative to plasma is shown in the one-dimensional polyacrylamide gel (Fig. 1
, top left). The immunoaffinity-purified TTR was composed primarily of TTR with some higher-molecular-mass contaminants that did not affect the analysis. The arrow indicates the band corresponding to TTR at
13.7 kDa. A typical multiply charged spectrum of TTR immunoaffinity purified by our method is shown in Fig. 1
. The TTR spectrum overwhelms the underlying noise, which may be a result of the high-molecular-mass species shown in the gel. Additionally, a total-ion chromatogram indicated a single chromatographic peak that is specific to TTR [see Fig. 6 in Bergen et al.(40), which appears in this issue]. Although the bands on the gel indicate the presence of other proteins, they did not interfere with obtaining a "clean" spectrum of TTR that was easy to interpret (Fig. 1
).
|
The spectrum in Fig. 1
is a LC-dual-ESI-FT-ICR mass spectrum of the immunopurified TTR from a patient diagnosed with ATTR. DNA sequencing verified that a transversion mutation of the DNA sequence produced a G47E variant in the amino acid sequence (data not shown). The broadband spectrum illustrates that there are two forms of TTR present: wild type and a variant form (G47E) producing a +72 Da intact molecular mass shift. The scan shown in the inset (top right) in Fig. 1
further demonstrates the effectiveness of using a reducing agent that makes the G47E variant readily distinguishable from cysteine adducts. The simplification of the mass spectrum of TTR is afforded by reduction, which can also be seen in Bergen et al.(40).
We considered using TCEP for reduction because of the low pH matrix (pH 2.7) needed to elute the TTR from the antibody resin. Getz et al.(45), who compared the reduction achieved with TCEP and dithiothreitol (DTT) for use in protein biochemistry, showed compelling evidence for the use of TCEP over DTT. Han and Han(46) reported that TCEP is a more effective reducing agent than DTT at lower pH values. Fig. 1
is an example of a TTR sample that was reduced by use of TCEP before analysis; no residual adducts can be detected. An added advantage of using TCEP was the time needed for reduction of TTR. A total of 15 min was required, compared with the 1.5 h needed for reduction by DTT (data not shown). The short reduction time facilitates high-throughput screening of TTR variants.
The total analysis time for a single measurement, as depicted in Fig. 1
, is
10 min and is largely a result of reestablishing the gradient. However, in a clinical environment, column-switching techniques could be used, which would reduce this time. Furthermore, decoupling the LC concentration/desalting step would allow use of flow-injection analysis requiring <30 s per sample.
determination of variants by use of accurate mass measurements
Comparison of accurate mass measurements of proteins with theoretical masses can reveal the presence of genetic variants. Although the most abundant isotopic mass of a protein or peptide is readily detectable experimentally, it may differ by up to 1 Da from the theoretical mass of the most abundant isotopic species because of the experimental variability of ion abundances(47). Monoisotopic mass is preferred for comparison because it is not subject to this uncertainty and because it does not depend on the characteristics of the isotopic distribution, which varies with mass. However, the monoisotopic masses of molecules larger than 5 kDa have extremely low abundance and are therefore not readily detected. In the present study, experimental monoisotopic masses were determined by an isotope-fitting procedure using the
2 statistic as our metric for the best fit. This approach, first described by Senko et al.(42), aligns a theoretical isotopic distribution to the experimental isotopic distribution. Several possible distribution alignments are obtained by shifting the theoretical isotopic distribution in
1/z increments relative to the experimental distribution.
For each possible alignment, a
2 value between experimental and theoretical abundances was calculated. Only pairs of peaks in which the experimental abundance was >20% of the most abundant isotopic peak in the distribution were included. The alignment with the lowest
2 value was taken to be optimal. The monoisotopic mass for each isotopic peak was then calculated by an abundance-weighted average of differences between experimental mass and theoretical monoisotopic mass(48). This procedure will be referred to as targeted isotopic fitting.
Whereas Senko et al.(42) used a model amino acid, averagine, whose chemical composition was obtained by averaging the relative occurrence of each amino acid in a database, we used the known amino acid composition of the wild-type protein to produce the theoretical distribution. The specificity of the antibody used for immunopurification of TTR before analysis ensured that the chemical composition of the experimental species closely resembled the wild type. Genetic variants obviously change the chemical composition and will therefore change the isotopic distribution, but, at least in the case of substitutions or small insertions and deletions, this change does not affect the targeted isotopic fitting procedure.
An expanded region of the [M + 11H+]11+ charge state for a mass spectrum that was acquired by LC-dual-ESI-FT-ICR is shown in Fig. 2
. The 11+ charge state was selected as it was typically the most abundant charge state and was bracketed by the IMCs. In Fig. 2
, two separate experimental isotopic distributions are observed; one corresponding to wild type (MMA = 0.6 ppm), the other clearly indicating the presence of a variant. The filled circles represent the theoretical isotopic distribution for this multiply charged ion, whereas the open inverted triangles show the value of the
2 statistic for a variety of alignments. The arrow for each distribution represents the calculated monoisotopic m/z values. The calculated experimental monoisotopic mass is then compared with a table of theoretical monoisotopic masses for known TTR variants. In this case, the calculated neutral monoisotopic mass (13 828.914 Da) corresponded well to the theoretical neutral monoisotopic mass of the S77Y variant (13 828.920 Da), giving a MMA of 0.4 ppm. DNA sequencing data were also used to verify that this was a S77Y variant (data not shown).
|
The results for seven individuals who were screened and tested positive for TTR variants are shown in Table 1
. Indicated in Table 1
are the neutral monoisotopic masses and the resulting mass errors. Once we had determined that a variant was present, we performed DNA analysis to confirm the results (data not shown). The mean MMA achieved for all seven samples was 0.21 ± 1.2 ppm (mean ± confidence interval of the mean at the 95% confidence level) for variant forms. Six of the seven individuals were heterozygous; the mean MMA achieved for the wild-type protein was 0.48 ± 2.26 ppm (0.01 ± 0.03 Da mass error). In comparison, Bergen et al.(40), using a low-resolution single-quadrupole instrument, obtained an intraassay mass difference of 28.1 (0.3) Da [mean (SD)] for a Val94Ala variant measured 10 times the same day and an interassay comparison of 28.7 (0.2) Da for the same sample measured 5 times on day 5. When 49 wild-type TTR samples were analyzed on the same day, Bergen et al.(40) obtained a mean mass of 13 761.5 ± 0.2 Da (mean ± confidence interval of the mean) at the 95% confidence level using the same low-resolution MS. Only the FT-ICR instrumentation provided the requisite high resolving power and mass accuracy to ascertain the mass differences of TTR containing a
2 Da mass shift, whereas the instrument with lower resolving power detected variants with mass shifts
10 Da(40).
|
Of the >90 known TTR genetic variants, it is possible that 8 variants would give a false-negative outcome. Two of the eight are isobaric (i.e., I
L or Q
K) with the wild type and would not be distinguished by mass spectrometric screening methods. The remaining six differ from the wild-type molecular mass by 1 Da. Identifying individuals homozygous for any ±1 Da variant becomes a feasible task because of the high MMA achieved. However, a majority of the individuals who have been diagnosed with ATTR are heterozygous, suggesting that there would be two forms of the protein, a wild type as well as the protein variant (± 1 Da). A mass-resolving power (m/
m) of >300 000 (full-width half-maximum; FWHM) must be achieved to distinguish the overlapping isotopic distributions from the wild-type and variant TTR proteins (e.g., to resolve the monoisotopic mass of a +1 Da variant from the wild-type protein with one 13C incorporated into it
m would need to be <0.05 Da). The typical resolving power obtained using the experimental conditions stated earlier was
35 000FWHM.
The remaining 82 genetic variants (91%) can be unambiguously detected based on mass. However, some variants are isobaric with each other, and the location of the variant cannot be determined from intact accurate mass measurements alone. Of the 82 variants, 27 (33%) can be conclusively located. Determining the site of the remaining 55 isobaric variants (77%) could be important to correlate with disease progression.
detection of novel cis double mutation
The [M + 11H+]11+ charge state of a theoretical wild-type isotopic distribution and an experimental isotopic distribution are shown in Fig. 3A
. The single experimental isotopic distribution in this spectrum was first assumed to be wild type. However, the neutral monoisotopic mass was calculated to be 13 753.882 Da (72 ppm MMA when compared with the theoretical wild-type neutral monoisotopic mass of 13 752.887 Da), which is significantly outside of our achievable limits (± 3 ppm) using the dual ESI source for internal calibration and fitting the isotopic distribution to obtain monoisotopic mass. These data clearly suggested that this individual was not homozygous wild type. DNA sequencing (Fig. 3B
) was performed to determine whether something other than an aberration in our methodology could explain this poor MMA. Indeed, the DNA sequencing data showed that there are two variants present in this particular individual, a G6S (+30 Da shift in protein mass) and a V30A (28 Da shift in protein mass).
|
Because ATTR is a familial disease, it is important to know from whom the variant(s) originated. If a double mutation were to occur in trans (the maternal and paternal chromosome would each contain a single mutation), there would be two distinct isotopic distributions for each variant; however, this is not observed (Fig. 3A
). If the mutations were to occur in cis (both mutations reside on the same chromosome), the resulting protein net mass shift would be +2 Da. If an individual is heterozygous, two isotopic distributions would be observed: one representing the wild type and the other the cis variant. The resolving power necessary to distinguish the wild type and the cis (G6S and V30A) variant is not typically achieved, providing a composite isotopic distribution, as observed (Fig. 3A
). Additional work is continuing on computational methods to routinely determine the presence of significantly overlapping isotopic distributions. In particular, we are exploring the use of the figure-of-merit calculation from THRASH(48) in an examination of how the signal-to-noise ratio affects the limits of detection for admixtures of small mass shifts in proteins of appreciable molecular mass.
Although MMA did not uniquely identify the genetic variant present, the combination of DNA sequencing and mass spectrometric results clearly indicated that the individual is heterozygous wild type/cis (G6S and V30A). Fig. 3C
shows the pedigree for the patient analyzed, the proband individual (index case) presented with typical ATTR symptoms in his late twenties. The immediate family (mother and brother) tested negative for the genetic variants of TTR. The father could not be tested because he is deceased, but his death had been attributed to amyloidosis. Because it is known that the mother was wild type, it can be concluded that the cis (G6S and V30A) variant originated from the paternal chromosome.
Double mutations are very rare: there has been only one other report, by Ando et al.(49), of an individual with a double TTR mutation. This particular finding showed that the individual carried a heterozygous trans variant. Ours is the first report of a cis double mutation of TTR. This particular cis double mutation was detected by a combination of genomic and proteomic approaches.
In conclusion, we used a targeted proteomics approach to characterize TTR variants by use of LC-dual-ESI FT-ICR MS. The high mass-resolving power translates into high MMA when an isotopic fitting procedure is used. Genetic variants of TTR can be detected effectively by use of quadrupole MS platforms(40), but unambiguous results can be obtained for >91% of known TTR genetic variants when FT-ICR MS is used for detection. Identification of variants is achieved routinely and accurately because of the high MMA (<3 ppm) largely afforded by the high mass-resolving power, the dual ESI source, and the isotopic fitting procedure. Importantly, we could determine, with contributions from both protein analysis and DNA sequencing, that an individual was heterozygous wild type and cis (G6S and V30A) variant. This is the first report of an individual with a cis double variant.
| Acknowledgments |
|---|
| Footnotes |
|---|
| References |
|---|
|
|
|---|
The following articles in journals at HighWire Press have cited this article:
![]() |
G. L. Hortin, N. Seam, and G. T. Hoehn Bound Homocysteine, Cysteine, and Cysteinylglycine Distribution between Albumin and Globulins Clin. Chem., December 1, 2006; 52(12): 2258 - 2264. [Abstract] [Full Text] [PDF] |
||||
![]() |
G. L. Hortin The MALDI-TOF Mass Spectrometric View of the Plasma Proteome and Peptidome Clin. Chem., July 1, 2006; 52(7): 1223 - 1237. [Abstract] [Full Text] [PDF] |
||||
![]() |
L. Anderson Candidate-based proteomics in the search for biomarkers of cardiovascular disease J. Physiol., February 15, 2005; 563(1): 23 - 60. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. Kinter Toward Broader Inclusion of Liquid Chromatography-Mass Spectrometry in the Clinical Laboratory Clin. Chem., September 1, 2004; 50(9): 1500 - 1502. [Full Text] [PDF] |
||||
![]() |
H. R. Bergen III, S. R. Zeldenrust, M. L. Butz, D. S. Snow, P. J. Dyck, P. J. B. Dyck, C. J. Klein, J. F. O'Brien, S. N. Thibodeau, and D. C. Muddiman Identification of Transthyretin Variants by Sequential Proteomic and Genomic Analysis Clin. Chem., September 1, 2004; 50(9): 1544 - 1552. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |