Abstract

Bordetella pertussis (Bp) is the causative agent of pertussis, a vaccine preventable disease occurring primarily in children. In recent years, there has been increased reporting of pertussis. Current pertussis vaccines are acellular and consist of Bp proteins including the major virulence factor pertussis toxin (Ptx), a 5-subunit exotoxin. Variation in Ptx subunit amino acid (AA) sequence could possibly affect the immune response. A blind comparative mass spectrometric (MS) analysis of commercially available Ptx as well as the chemically modified toxoid (Ptxd) from licensed vaccines was performed to assess peptide sequence and AA coverage variability as well as relative amounts of Ptx subunits. Qualitatively, there are similarities among the various sources based on AA percent coverages and MS/MS fragmentation profiles. Additionally, based on a label-free mass spectrometry-based quantification method there is differential relative abundance of the subunits among the sources.

1. Introduction

Bordetella pertussis (Bp), a gram-negative bacterium, is the causative agent of pertussis, commonly known as whooping cough [1]. Pertussis is a highly communicable respiratory infection that in the absence of vaccination primarily afflicts infants and children. Transmission occurs via direct contact with discharges from respiratory mucous membranes of infected persons. Persons who contract pertussis develop a persistent cough that includes the characteristic high-pitched “whoop” sound that can last up to eight weeks [1, 2]. Pertussis is a vaccine preventable disease but according to the World Health Organization (WHO), there are an estimated 50 million cases of pertussis per year worldwide with approximately 300,000 leading to death [3, 4]. The majority of the high mortality rates come from developing countries. In the United States acellular pertussis vaccines are administered to infants and young children [5] and in 2006 were recommended for adolescents and adults [68]. Acellular vaccines consist of up to five Bp components including the major virulence factor pertussis toxin (Ptx), filamentous hemagluttinin, pertactin, fimbrae 2, and/or fimbrae 3 [3]. Pertussis vaccines are generally formulated in combination with vaccine immunogens for other diseases, such as tetanus and diphtheria [9].

Ptx, a five-subunit protein complex (S1–S5) is a major bacterial exotoxin, surface adhesin, and porin, is involved in the establishment of infection in humans and has been shown to be an important component of the vaccine [10]. Ptx, a 105 kDa A-B protein complex, comprises the protomer (A) and oligomer (B). The A component consists of S1 ( 22 kDa), which contains ADP ribosylation enzymatic activity and is the most immunoreactive region of the toxin. The B portion consists of S2, S3, S4, and S5, decreasing in size from 20 kDa to 10 kDa, respectively, and is associated with adherence to human host cells [11, 12]. From a cellular perspective, the B portion binds to cells inducing the enzymatic activity of S1 and its transversion into the cell, followed by inactivation of G proteins by S1 and subsequent initiation of an intracellular cascade of signaling events which includes factors that regulate the immune response [11]. Importantly, the Ptxs in current vaccines have been inactivated (toxoids; Ptxd) by covalent crosslinkages formed using agents such as formaldehyde or glutaraldehyde [13, 14].

Recently, the U.S. and several other countries have reported an increase in pertussis cases. Factors such as waning immunity after natural infection or immunization, improved molecular-based diagnostics, increased surveillance, or changes in circulating strain antigenicity have all been hypothesized as contributing factors [3] and have led researchers to re-examine pertussis vaccine components and preparation strategies. Changes in Ptx subunit protein sequences as well as the amounts present in the vaccine formulation could affect the overall immune response. Moreover, identification of commonalities as well as potential disparities of various Ptx and Ptxd sources on a proteomic level could provide insight into changes in population immune responses to the vaccine.

Mass spectrometry (MS), a powerful and sensitive analytical tool has been used to identify, differentiate, and characterize proteins in complex mixtures [15]. In particular, MS techniques, such as liquid chromatography-electrospray ionization (LC-ESI) tandem MS (MS/MS) or matrix-assisted laser desorption ionization time-of-flight (MALDI-TOF) MS/MS, which vary in their rapidity and sensitivity of peptide detection as a result of enzymatic cleavage, are commonly used in proteomic analysis [15, 16]. Proteins can be proteolytically cleaved by enzymes such as trypsin and the peptides generated are ionized and separated based on their mass-to-charge ratios. The peptides and the proteins, from which they originated, are identified based on comparisons between the MS/MS fragmentation patterns and protein databases. [17]. For instance, Bottero et al. used MALDI-TOF MS/MS as a means to characterize enriched Bp membrane proteins, which included the Ptx S1, from different strains [18]. Alternatively, Tumala and colleagues have used liquid chromatography MS/MS (LC-MS/MS) to accurately measure the molecular weights of Ptx subunits as well as to identify novel chemical modification sites [19].

In this study, a blind comparative mass spectrometric analysis of Ptx and Ptxd from multiple sources was performed to assess potential sequence and amino acid (AA) coverage variability as well as relative amounts of pertussis toxin subunits. This first-phase study demonstrates the sensitivity of mass spectrometry in ascertaining differences at the peptide sequence level and may provide a basis for future correlative immunoproteomic studies.

2. Methods

2.1. Pertussis Toxin and Toxoid Sources

The commercially available Ptx and Ptxd were purchased from Protein Express, Inc. (Cincinnati, OH, USA), LIST Biologicals, Inc. (Campbell, CA, USA), Sanofi Pasteur Pharmaceutical (Paris, France), or GSK Pharmaceutical (London, UK). The sources were blindly assigned a letter, in which Ptxs are represented as A or B, while the chemically inactivated Ptxd samples are denoted as C or D.

2.2. Pertussis Toxin and Toxoid Digestions

The Ptx and Ptxd sources (10  g) were vortexed and precipitated with 100% acetonitrile (ACN), for 60 minutes, and then dried using a vacuum centrifuge. Trypsin (5  g) (Promega Corporation, Madison, WI, USA) was added and the volume was brought to 100  L with tryptic digestion buffer (50 mM ammonium bicarbonate adjusted to pH 8.5 with NaOH, 0.6 mM CaCl2, 0.1% (w/v) potassium azide (KN3), and 3 mg/L phenol red). The samples were incubated at for 3, 6, and 24 hours. Digested samples (250 ng total) were subsequently prepared for either MALDI-TOF MS/MS or LC-MS/MS. Additionally, in separate experiments, samples were treated with an acid cleavable detergent, RapiGest (RG; Waters Corporation, Milford, MA, USA) to assess the efficiency of tryptic digestion. 0.1% RG in tryptic digestion buffer was added to Ptx or Ptxd (10  g) samples. The samples were next incubated at for 5 minutes followed by cooling at room temperature with the subsequent addition of trypsin (5  g). Trypsin digestions were performed as stated above. However, prior to MS analysis, the RG was inactivated by cleavage with 0.175 M HCL at for 30 minutes, and the supernatant (obtained by centrifugation for 15 minutes at 17,300 g) was transferred to a new vial.

2.3. MS Sample Preparation

Samples for MALDI-TOF MS/MS were resuspended in an equal volume of 0.1% trifluoroacetic acid (TFA), and the peptides were separated using a C18 zip-tip minicolumn (Millipore Corporation, Billerica, MA, USA) and eluted in two steps with 3  L of 30% ACN in 0.1% TFA and then 3  L of 80% ACN in 0.1% TFA. The eluates were mixed with 10  L -cyano-4-hydroxycinnamic acid ( -CHCA) (Sigma Aldrich, St. Louis, MO, USA) and 0.5  L was next spotted on to a stainless steel plate. A protein mix diluted at 1 : 10 was spotted on the plate in the appropriate wells and used for instrument calibration purposes. For LC-MS/MS, samples were simply resuspended in equal volumes of 0.1% formic acid (FA).

2.4. Mass Spectrometry
2.4.1. MALDI-TOF MS/MS

Mass spectra of each sample well were obtained by scanning from 800 to 4000 m/z in MS positive-ion reflector mode on the Applied Biosystems 4800 Proteomics Analyzer (Framingham, MA, USA). The instrument uses a nd:YAG laser at 337 nm with a 200 MHz repetition rate, and each spectrum is an average of 1000 laser shots. In MS/MS mode, the CID spectra were generated using air as the collision gas at 1 kV and an average of 500 laser shots for MS/MS acquisition. The instrument was tuned and calibrated for MS analysis using a mixture of five peptides: des-Arg1-Bradykinin (m/z 904.47), angiotensin I (m/z 1,296.69), Glu1-fibrinopeptide B (m/z 1,570.68), ACTH (1–17)(m/z 2093.08), and ACTH (18–39)(m/z 2,465.20). For MS/MS mass calibration, known fragment ion masses observed in the Glu1-fibrinopeptide B were used. The instrument was programmed for automated MS/MS acquisition using data dependent acquisition (DDA) mode. Selected peptides were subjected to MS/MS analysis to confirm their identity through sequence determination.

All MS/MS spectra were automatically submitted to a Mascot Server through Mascot Distiller (Matrix Science Inc., London, UK; version 2.2.1.0) and searched against the entire NCBInr database or a modified NCBInr database created to search “Bordetella” or “pertussis” proteins with deamidation, oxidation and carbamidomethylation assigned as variable modifications. Tryptic constraints were applied, allowing up to two missed cleavages. Peptide searches were performed with an initial tolerance on mass measurement of 100 ppm in MS mode and 150 ppm in MS/MS mode. All matched peptides were verified manually and the searches combined and processed by Scaffold (Proteome Software Inc., Portland, OR, USA; version 2.03.01).

2.4.2. Nano-LC ESI-MS/MS (LTQ-Orbitrap)

Protein identification was achieved by using nanoflow liquid chromatography, data-dependant tandem mass spectrometry, and database searching. A 365  m by 75  m fused silica pulled needle capillary (New Objective, Inc., Woburn, MA) was filled in house with 10 cm of 5  m Symmetry 300 reverse phase packing material (Waters Inc., Bedford, MA). Protein digests were directly loaded onto the analytical column without using a trap column and introduced using an Eksigent autosampler. The gradient was delivered by an Eksigent 2D nanoLC system (Eksigent Technologies, Inc, Dublin, CA) where the mobile phase solvents consisted of: (solvent A) 0.2% formic acid (Thermo Scientific, Rockford, IL), 0.005% trifluoroacetic acid (Sigma-Aldrich) in water (Burdick and Jackson, Muskegon, MI), (solvent B) 0.2% formic acid, and 0.005% trifluoroacetic acid in acetonitrile (Burdick and Jackson). The gradient flow was set at 400 nL/min, and the profile consisted of a hold at 5% B for 5 minutes followed by a ramp to 30% B over 100 minutes then a ramp up to 90% B in 5 minutes and a hold at 90% for 2 minutes before returning to 5% B in 2 minutes and re-equilibration at 5% B for 20 minutes. Peptides were nano-electrosprayed into an LTQ Orbitrap tandem mass spectrometer (Thermo Scientific, San Jose, CA). The voltage applied to the nano-electrospray source was 2.0 kV. The mass spectrometer was programmed to perform data-dependant acquisition by scanning the mass range from m/z 400 to 1600 at a nominal resolution setting of 60,000 for parent ion acquisition in the Orbitrap. Then tandem mass spectra of doubly-charged and higher-charge states were acquired for the top 10 most intense ions. Singly-charged ions were not considered for MS/MS. All tandem mass spectra were recorded using the linear ion trap. This DDA process cycled continuously throughout the duration of the gradient.

All tandem mass spectra were extracted from the raw data file using Mascot Distiller (Matrix Science, London, UK; version 2.2.1.0) and searched using Mascot (version 2.2.0). Mascot was setup to search using the entire NCBInr database or a modified NCBInr database created to search “Bordetella” or “pertussis” recognized proteins in which trypsin is used as the digestion agent. Mascot was searched with a fragment ion tolerance mass of 0.80 Da and a parent ion tolerance of 200 ppm. Variable modifications for the mascot search were deamidation, oxidation, and carbamidomethylation. Scaffold was used to validate MS/MS-based peptide and protein identifications. Peptide identifications were accepted if they could be established at greater than 95.0% probability as specified by the Peptide Prophet algorithm [20]. Protein identifications were accepted if they could be established at greater than 99.0% probability and contained at least 2 identified peptides [21]. With these stringent parameters of Peptide Prophet and Protein Prophet, the false discovery rate was zero.

2.4.3. Nano-LC ESI-MS/MS (Q-Tof Premier)

LC-MS/MS was carried out using a nanoAcquity ultra-performance liquid chromatography (UPLC) coupled to a QTof Premier MS system (Waters Corporation). LC separation was performed using a Symmetry C18 trapping column and a BEH C18 column (100  m I.D. 100 mm long with 1.7  m packing), at a flow rate of 1.2  L/min. Standard water/ACN gradients containing 0.1% formic acid in both solvents were used for elution, with the ACN increasing linearly from 1% to 50% in 50 minutes, and an overall cycle time of 90 minutes. The QTof utilized an (or Protein Expression) method, which involved acquiring data-independent alternating low- and high-energy scans over the m/z range 50–1990 in 0.6 second, along with lockmass data on a separate channel to obtain accurate mass measurement. The data obtained using the method were further processed using ProteinLynx Global Server v2.3, which provided statistically validated peptide and protein identification using a SwissProt-Trembl protein database, along with relative protein quantification analysis. Root mean square mass accuracies were typically within 8 ppm for the MS data, and within 15 ppm for MS/MS data. Also the false positive rate was 4%. Lastly, the relative quantification assessed in this study used an internal standard from a known protein and is considered a surrogate.

3. Results

A mass spectrometric analysis was performed to assess protein sequence homogeneity of Ptx and Ptxd. MALDI-TOF MS/MS, and LC-MS/MS approaches revealed comparable to varying AA percent composition profiles (percent is based on the number of AA (within tryptic peptides) generated divided by the total number of AA in the protein) among Ptx and Ptxd (Table 1) as well as a compilation of detected peptides (Table 2). Generally, Ptx A and B digested with trypsin generated higher AA coverages among the nontreated oligomeric subunits compared to their chemically modified Ptxd counterparts regardless of the MS system used. Additionally, trypsin digests of Ptx B generated the highest average among S2–S5 compared to all samples. Also, for both Ptxd samples there was low to no detection of S2–S5 peptides.

However, AA coverage reports for S1 on average between Ptx and Ptxd were fairly comparable with 45 and 44%, respectively. Probing further, LC-MS/MS instruments proved to be more fruitful in terms of yielding a higher-sequence coverage and an increase in peptide detection. Specifically, S1 gave a higher-percent coverage by LC-MS/MS instruments resulting in a 1.7-fold higher-overall average as compared to the MALDI-TOF mass-spectrometric approach.

Tryptic digestions supplemented with RG resulted in an increase in sequence coverage in all samples and differentially within the subunits (Table 3). RG had a greater effect on the Ptxd giving rise to the greatest increase in sequence coverage. Ptxd D resulted in a slightly greater overall fold change as compared to Ptxd C in the presence of RG, while fold changes between all subunits were lower in both Ptx A and Ptx B Collectively for both Ptx and Ptxd, S1, S3, and S5 had marked gains in the presence of RG with a nearly 1.5, 2.5, and 10-fold increase, respectively.

As the sequence coverage increased in the presence of RG, the number of peptides detected by LC-MS/MS increased as well and to a lesser extent in MALDI-TOF (data not shown). Since the LC-MS/MS LTQ Orbitrap instrument gave the optimum sequence coverage (in particular for the toxins) in the presence of RG, a comparative assessment of peptide detection for Ptx subunits among all sources was compiled (Table 4). Table 4 reveals common and unique peptides that were detected among Ptx and Ptxd sources analyzed. More importantly, an in-depth visual scan analysis of MS/MS spectra of the common peptides revealed similar fragmentation profiles (based on y- and b-ion distribution and relative intensity of fragment ions), suggestive of Ptx subunit sequence homogeneity as illustrated in Figure 1.

In addition to the determination of sequence homogeneity among sources, the relative amounts of protein as measured by relative abundance (fmol) of subunits from Ptx and Ptxd was preliminarily explored (Table 5). An LC-MS/MS approach was used, in which the Q-Tof data were employed for this assessment in which differential abundance was revealed among sources and the subunits. The absolute abundance of Ptx B was 7.2-fold higher among all subunits compared to Ptx A. Ptx B S4 gave the highest fold change at nearly 19 in comparison to Ptx A. With respect to Ptxd, the abundance was nearly the same for all of the subunits, with Ptxd D S1 revealing the highest fold change versus Ptxd C.

4. Discussion

We initiated a study using mass spectrometry to determine if AA sequence differences exist between Ptx and Ptxd when examined from different sources. Using proteins (i.e., Ptx) from strains that may not share the exact same AA sequence or that have undergone different inactivation treatments could result in differences in immune effectiveness. Protein sequence variation between strains may be small (i.e., single AA addition, deletion, or substitution) or potentially more extensive with the addition or deletion of immunogenically relevant motifs or domain structures. As a result, protein sequence variance could lead to subtle structural, conformational, and/or functional differences that can be detected only through the use of sensitive technologies.

Previous MS investigations have been used for identification, differentiation, and characterization of protein complexes. For example, Nandakumar et al. [22] utilized MS techniques for identification of bacterial proteins. Protein mixtures were subject to tryptic in-solution digestion and analyzed by reverse phase nano LC-MS/MS followed by database analysis. Other investigators have used MALDI-TOF MS/MS for the analysis of peptides generated from the enzymatic cleavage of Clostridium botulinum C3 coenzyme [23]. Also, Yoder and colleagues analyzed tryptic digests of Vaccinia virus using MALDI-TOF, as well as a quadrupole ion trap mass spectrometer and a quadrupole TOF to identify proteins associated with infection [24]. Additionally, nano-LC-MS/MS was used to identify AA variability between two-bovine surfactant proteins [25]. In our study, an analytical and proteomic evaluation was conducted to assess any composition (i.e., AA sequence of the protein components) and formulation (i.e., how much of a specific component) disparities among Ptx and Ptxd.

Based on AA coverages and peptide detection alone, the two Ptxs and two Ptxds appear to have similar sequences, an indication that they were derived from the same origin, although different purification schemes and inactivation methods most likely were employed. The fact that the Ptx and Ptxd sources shared similar tryptic digest patterns, in terms of the peptides being generated and detected as well as spectrum profiles, suggests homogeneity. The presence of RG improved protein digestion resulting in higher-percent coverages generally across source as well as Ptx subunits. In addition, RG treatment appeared to have exposed the smallest-sized S5, resulting in an increase in coverage for this subunit. In doing so, additional peptides were detected that once again were shown to be similar among Ptx and Ptxd sources based on y- and b-ion distribution and relative intensities of the fragmented peptides. Moreover, studies by Sekura et al. [26] revealed, by amino acid analysis, that Ptx purified from Bp strain 165 was indistinguishable from Ptx prepared from the well-characterized Bp Tahoma strain. In addition, Tummala et al. confirmed Ptx sequence homogeneity among their Bp culture supernatant purified Ptx as compared with previously reported sequence information [19]. These data are confirmatory and strengthen the MS results in this study for Ptx sequence homogeneity. However, although the data suggest protein sequence homogeneity, if the sequence is indeed different but not present in the database then, upon searching, the generated peptides would not be detected resulting in a lower peptide map coverage.

Ptx is a major virulence factor and in its natural state has physiological inhibitory effects [11]. The subunits are highly immunoreactive, including the dominant S1, making it an ideal component for vaccines. To reduce the cellular lethality induced by Ptx, the Ptx is inactivated or detoxified in its final vaccine form. Inactivation occurs by chemical modification (i.e., formaldehyde or glutaraldehyde) via the formation of covalent crosslinkings on targeted AA (usually lysines). Interestingly in Ptxd samples, the oligomeric containing subunits percent AA coverage was lower. The S2–S5 subunits, rich in lysine are the subunits that are known to be heavily modified due to the chemical action exhibited by formaldehyde and/or glutaraldehyde [13, 27, 28]. As a result of modification, these peptides and subsequent MS-MS profiles would thus not be detected using traditional selected database searching tools resulting in lower AA coverages.

Although formaldehyde, for instance, is an ideal reagent for chemical crosslinking, it also immobilizes the antigen (i.e., Ptx) and alters its conformation preventing a natural interaction between the antibody and antigen. In other words, formaldehyde reduces the accessibility of antibodies to their epitopes by creating steric hindrance or by altering the epitope site itself. In Ptxd samples, inactivation by formaldehyde appears to have a severe impact on peptide detection, due to crosslinking of lysine on the Ptx subunits [13, 2729]. These subunits are most likely masked preventing trypsin, whose natural recognition sites are lysine and arginine, to effectively cleave the protein resulting in a lower AA coverage. On a cellular level, it is sensible for these subunits to be inactivated heavily, since in vivo, they make up the oligomeric complex that adheres to the human host cells which upon binding to the appropriate receptor initiates Ptx function. Blocking these cellular activities, in essence, inhibits or diminishes further downstream enzymatic activity by S1, which would allow for progression of normal immune responses by the host towards Bp. Besides, S1 has consistently had higher AA coverage reporting among all sources, which could be explained by the fact that S1 lacks lysines. Although S1 is deficient in lysine, free cellular lysine in the presence of formaldehyde would still create crosslinking bridges that lead to the detoxification [28, 29] but to a lesser extent compared to S2–S5. As a result of minimal chemical bombardment and thus potentially incomplete inactivation by formaldehyde on the toxic S1 component, residual toxicity may remain, with varying levels among Ptxd sources, could lead and be an explanative factor for differential immunity, albeit speculative.

5. Conclusions

In conclusion, this study incorporated multiple mass-spectrometric approaches to identify protein sequence similarities among Ptx and Ptxd sources, based on AA coverage profiles, peptide detection, and fragmentation profiles. Lower-percent coverage among Ptxd oligomeric subunits as compared to S1 could be explained most likely due to chemical modifications. As a result of these modifications, inactivation of Ptx subunits, as evident by failed detection of peptides, could lead to diminished immuno-cellular interactions between human-generated antibodies against chemically modified Bp antigens. In addition to the protein components, changes in the formulation, including the concentration, ratio and abundance of the Ptx subunits could translate into differences with respect to immune responses and vaccine efficacy. Moreover, information from this MS study could prove useful in improving immunological tests that use Ptx as an examining factor for vaccine evaluation and even pertussis diagnosis.

Future studies such as de novo sequencing or protein sequence tag analysis may be considered to validate and strengthen the current finding that the protein sequences among all the sources are indeed homogeneous. In addition, incorporation of other enzymes in the presence of RG to expand the AA coverage and peptide detection may also be performed. Also, a comparative quantitative MS analysis assessing the functionality of Ptx and Ptxd sources may be able to determine the effect, if any, that Ptx inactivation, Ptx function, and potential residual toxicity play on vaccine efficacy such as immune responses. Lastly, in vivo animal studies may be conducted to examine the correlation between Ptx protein sequence homogeneity and immunoreactivity.

Disclaimer

The references in this paper to any specific commercial products, process, service, manufacturer, or company do not constitute an endorsement or a recommendation by the U.S. Government or the Centers for Disease Control and Prevention. The findings and conclusions in this paper are those of the authors and do not necessarily represent the views of CDC.