|
|
||||||||
1 Departments of Chemistry, 2 Biochemistry, and 3 Physics, Michigan State University, East Lansing, Michigan 48824, USA
Reprint requests to: J. Throck Watson, Departments of Chemistry and Biochemistry, Michigan State University, East Lansing, MI 48824, USA; e-mail: watsonj{at}msu.edu; fax: (517) 353-9334.
(RECEIVED March 15, 2005; FINAL REVISION April 18, 2005; ACCEPTED April 27, 2005)
| Abstract |
|---|
|
|
|---|
Keywords: secondary structure; hydrogen/deuterium exchange; mass spectrometry; collision-induced dissociation; autonomous folding unit; structured peptides; bovine pancreatic trypsin inhibitor (BPTI)
Abbreviations: H-bond, hydrogen bond HDX, hydrogen/deuterium exchange NMR, nuclear magnetic resonance MS, mass spectrometry MALDI-TOF, matrix-assisted laser-desorption/ionization time-offlight CID-MS/MS, collision-induced dissociation tandem mass spectrometry BPTI, bovine pancreatic trypsin inhibitor TCEP, tris-(2-carboxyethyl)phosphine hydrochloride TFA, trifluoroacetic acid GdnHCl, guanidinium hydrochloride HPLC, high-performance liquid chromatography CD, circular dichroism FA, formic acid
Article and publication are at http://www.proteinscience.org/cgi/doi/10.1110/ps.051458905.
| Introduction |
|---|
|
|
|---|
-helices,
-hairpins,
-meanders, etc.) (definitions of the three types of structure discussed in this worksecondary structure, tertiary structure, and side-chain structureare provided in the Supplemental Material) from naturally occurring proteins exhibit little secondary structure in aqueous solution (Montelione and Scheraga 1989), probably owing to competition between wateramide and amideamide hydrogen bonds (Hbonds 1) (Harrington and Schellman 1956; Klotz and Franzen 1962). Therefore, long-range interactions are usually invoked to account for the stabilization of short-range structure in protein folding. However, the mechanisms by which interactions between peptides stabilize the structure within peptides have not been studied systematically. It would be useful, for example, to assess the correlations between secondary, side-chain, and tertiary structure in pairs of covalently linked peptides and to characterize quantitatively the stabilization of local structure conferred by various factors, e.g., amino acid sequence, specific side-chain packing, steric, and/or hydrophobic exclusion of competing water molecules, etc. Hydrogen/deuterium exchange (HDX), as monitored by mass spectrometry (MS) and nuclear magnetic resonance (NMR) spectroscopy, has been used widely to study conformations and conformational dynamics of proteins in various conditions ranging from isolated states to noncovalent interactions (Englander 2000; Kaltashov and Eyles 2002). However, HDX has been rarely applied to study oligopeptide systems (Carulla et al. 2000). We chose to measure HDX using MS because it requires significantly less material than NMR and tolerates poorly soluble peptides. Collision-induced dissociation tandem mass spectrometry (CID-MS/MS) can assess the deuterium incorporation of individual amides from the b-series of CID fragment ions (Deng et al. 1999; Kim et al. 2001; Hoerner et al. 2004). Recently, HDX-CID-MS/MS results were reported in a study of two oligopeptides and their homodimers formed by noncovalent interactions using electrospray ionization with an ion-trap mass spectrometer (Cai and Dass 2005).
The current understanding of protein folding suggests three connected hypotheses that can be tested by HDX on pairs of covalently linked peptides:
These hypotheses will be tested by comparing the structure in monomeric and dimeric peptides with HDX and circular dichroism (CD). To enforce the association of designated peptides, we restrict our study here to (1) individual monocysteinyl peptides that are linked by a disulfide bond in the native structure of a folded protein and (2) all possible covalent dimers of these peptides.
We illustrate application of our method using two cysteinyl peptides and their three possible disulfide-linked pairs, denoted P
, P
, P
, P
, and P
. The model peptides mimic selected segments of the sequence of bovine pancreatic trypsin inhibitor (BPTI); specifically, the monomers, P
and P
, correspond to residues 4358 (NNFKSAEDCMRTAGGA) and 2033 (RYFYNAKAGLCQTF), respectively, which form the disulfide-linked C-terminal
-helix and central
-hairpin, respectively, in native BPTI (Fig. 1A
, Supplemental Material). Because P
and P
each have a cysteine, they can be disulfide-bonded to form the P
and P
homodimers and the P
heterodimer. The P
heterodimer corresponds to the (Cys30Cys51) disulfide species of intact BPTI, which accumulates in the onedisulfide ensemble and seems to be partially folded (Weissman and Kim 1991; Staley and Kim 1994). Oas and Kim have shown that the P
heterodimer has stable secondary structure by CD and 1D NMR, whereas the P
and P
monomers do not (Oas and Kim 1988). The P
and P
dimers have not been studied structurally. However, it was reported recently that the homodimer of another peptide model for the BPTI
-sheet folds into a four-stranded anti-parallel
-sheet when the two monomers are connected by an "optimized" synthetic linker, although homodimers formed using "unoptimized" linkers and the monomer peptide do not form stable structure (Carulla et al. 2002). Throughout this paper, the residues are numbered according to their position in native BPTI (e.g., Cys30, Ser47, etc.).
|
| Results and Discussion |
|---|
|
|
|---|
and P
peptides showed average mass shifts of 10.2 and 8.6 Da, respectively, from their 15 and 13 exchangeable hydrogens. Hence, P
and P
show some protection against HDX under our experimental conditions (see Materials and Methods). This protection may reflect "flickering" local structure, as has been observed in other experiments (Goodman and Kim 1989; Carulla et al. 2000) and in our CD and CID-MS/MS experiments reported below. Much smaller mass shifts of 6.8 and 3.5 Da were measured for the subunits of the P
and P
dimers, respectively; the native-like P
dimer showed slightly smaller mass shifts of 5.5 Da (
subunit) and 2.9 Da (
subunit). Hence, all three dimers showed significantly more protection than their component monomers, consistent with hypothesis H1.
|
and P
(Fig. 1A
and P
(Fig. 1B
(Fig. 1C
, P
, and P
) relative to that in their component monomers (P
and P
). The clearest evidence of secondary structure is visible in P
(Fig. 1C
exhibits two minima near those expected for helical proteins (208 nm and 222 nm) (Woody et al. 1996); at 60°C, these wells are replaced by one well at 204 nm characteristic of random-coil structure (Woody et al. 1996). The summed CD spectrum of the P
and P
peptides resembles the spectrum of P
at 60°C, but shifted to lower wavelengths; this contrasts with results from the earlier study (Oas and Kim 1988), which found little deviation. Comparison of P
and P
shows that P
at 0°C has more
-helical structure than P
at 60°C, which has more
-helical structure than P
at 0°C (judging by ellipticity at 222 nm). Similarly, P
at 0°C has significantly more secondary structure than its corresponding monomer P
at 0°C (judging by ellipticity at 217 nm) (Woody et al. 1996). These results are consistent with hypothesis H2 that local secondary structure content is stabilized upon the association of two peptides, presumably by the exclusion of water from attacking the amideamide H-bonds. However, the increase in secondary structure in P
seems small and local, not global; the P
CD spectrum is not fully
-helical even at 0°C, and the difference spectrum of P
and P
does not exhibit the double minima at 208 nm and 222 nm expected for
-helical structure (Fig. 1A
against HDX (which rivals that of the native-like P
) likely stems from the association of the two P
monomers connected by the disulfide bond, rather than from enhanced secondary structure within each, which supports hypothesis H1. By contrast, the CD spectrum of the P
peptide at 0°C (but not at 60°C) suggests the presence of side-chain structure; specifically, the positive peak at 230 nm likely corresponds to the CD signal of one or more ordered tyrosine side chains (Sreerama et al. 1999). This suggests that P
may also have stable tertiary structure, similar to that in P
; further evidence is cited below.
To identify the position(s) of the enhanced HDX protection induced by the dimerization in the
-unit, we used CID-MS/MS to probe the HDX of individual amide groups in P
, P
, and P
. These measurements were done using direct infusion electrospray ionization on an LCQDeca ion trap mass spectrometer. HDX-CID of b ions has been reported to be reliable for assessing deuterium content at the single amino acid level, whereas that of y ions is unreliable (Deng et al. 1999; Kim et al. 2001; Hoerner et al. 2004). A series of b ions (b3b10) was formed during the CID process for the P
subunit; however, only y ions were produced for the P
subunit. The masses of the b3b10 ions were used to assess deuterium incorporation in residues Tyr23Leu29 of the P
subunit in P
, P
, and P
(Table 2
). The results indicate that the amide hydrogens of Tyr23, Asn24, Lys26, and Ala27 exchange readily in P
(0.61.1 Da), but are protected in P
(0.10.3 Da) and even more so in P
(0.0 Da for Lys26, 0.1 Da for Ala27). By contrast, the Ala25 amide hydrogen exchanges equally well and moderately in P
, P
, and P
(0.4 Da). Finally, the Gly28 and Leu29 amide hydrogens exchange weakly in P
(0.2 Da) and less so in P
and P
(0.1 and 0.0 Da, respectively). Due to the missing fragment of the b1, b2, and b11b13 ions, deuterium incorporations for Tyr21, Phe22, and Cys30Phe33 are not available. However, the 2.0-Da mass shift of the b3 ion of P
indicates that Tyr21 and Phe22 amide hydrogens exchange readily in the monomer. In addition, the 1.0-Da decrease in the mass shift of the b3 ion in P
and P
suggests that one of these two amide hydrogens is protected in the dimers. Similarly, segment Cys30Phe33 of the
sequence shows less protection in P
(2.3 Da) than in P
and P
(0.9 Da and 1.0 Da, respectively).
|
subunit are generally consistent with the solvent exposure and H-bonding of the native structure (see Fig. 1B
-chain adopts its native structure, the amide hydrogens of Phe22, Asn24, Gly28, Gln31, and Phe33 should be well protected because they are all involved in strong backbone (amideamide) H-bonds; however, the Phe33Arg20 H-bond is unlikely to be maintained in the peptide model, owing to the well-known flexibility of chain termini. Consistent with these predictions, significant protection for Asn24 and Gly28 is observed in P
and P
(Table 2
and P
are consistent with significant protection of Phe22 and Gln31 against HDX, which are mutually H-bonded in the native structure. The amide hydrogen of Leu29 is well protected in P
and P
(Table 2
, suggesting that P
adopts local structure in its
hairpin (residues 510). Finally, the Lys26 and Ala27 amide hydrogens form a bifurcated H-bond with the side-chain oxygen of Asn24 in the native structure of BPTI. These hydrogens are well protected in P
and even more so in P
, consistent with a structured Asn24 side chain. These two hydrogens are chiefly responsible for the difference between the protection of the P
subunit against HDX in P
and P
. Combined with the CD results, these data suggest significant native secondary and side-chain structure in the P
homodimer, as well as the P
heterodimer (Oas and Kim 1988).
The absence of protection against HDX also can be a cross-check of native structure; the amide hydrogens of Ala25, Cys30, and Thr32 should exchange readily because they are solvent-exposed and not H-bonded in the native structure. The increased structure in P
and P
does not provide enhanced protection to Ala25 (Table 2
). Notably, Ala25 does not fully exchange even in P
; this effect remains to be explained, although slow exchange at this position has been noted in BPTI variants (Schulman and Kim 1994), and may reflect a stable
-turn at residues Asn24Leu29. The mass shifts of b14 and b10 from P
and P
suggests that at least one residue in the segment Cys30Phe33 is unprotected. Finally, the amide hydrogens of Tyr21 and Tyr23 should be protected in the P
dimer owing to intermolecular H-bonds with Phe22 and Asn20 from the P
subunit. Consistent with this prediction, no protection is observed for Tyr21 and Tyr23 in P
, but modest and roughly equal protection of Tyr23 is observed in both P
and P
(Table 2
). This protection of Tyr23 in the P
dimer suggests stable intermolecular H-bonding (and, thus, tertiary structure) in P
, in addition to the secondary and side-chain structure cited above.
| Conclusions |
|---|
|
|
|---|
Applied to peptide models from BPTI, our results provide evidence for the three hypotheses outlined in the introduction. First, P
shows some protection against HDX, but the P
dimer has almost as much protection as the native-like P
heterodimer, despite having significantly less secondary structure than P
(as assessed by CD). This supports hypothesis H1 that the association of two peptides can inhibit HDX by lowering the local activity of water independently of secondary structure. Second, the P
dimer exhibits more secondary structure than the P
monomer, even at 60°C; analogous results pertain to P
and P
. This supports hypothesis H2 that the decreased local activity of water stabilizes local secondary structure. Finally, the P
dimer does not appear to be fully structured, judging from its nonhelical CD spectrum; by contrast, the data suggest that the P
and P
dimers have stable secondary, side-chain, and even tertiary structure. Thus, our results are consistent with hypothesis H3 that stable secondary, side-chain, and even nonnative tertiary structure may develop when two sufficiently nonpolar peptides are linked covalently. The difference between the P
and P
dimers may result from the hydrophobicity of the two peptides; P
has only two large hydrophobic residues (Phe45 and Met52) of 16 residues, whereas P
has four (Tyr21, Phe22, Tyr23, and Leu29) of 14 residues. Therefore, the P
dimer is better able to form a hydrophobic core around its disulfide bond, further shielding its amideamide H-bonds from attack by water molecules. This is consistent with the global HDX results; P
exchanges 27% (3.5/13) of its amide hydrogens, whereas P
exchanges 45% (6.8/15).
| Materials and methods |
|---|
|
|
|---|
and P
peptides were synthesized by the Genomics Technology Support Facility (GTSF) at Michigan State University and purified by reverse-phase high-performance liquid chromatography (HPLC) using a Vydac C18 analytical column (catalog no. 218TP54). All other reagents were of the highest purity commercially available and used without further purification. Dibasic sodium phosphate (Na2HPO4) and sodium sulfate (Na2SO4) were purchased from Spectrum Chemical Mfg. Corp. and Columbus Chemical Industries, Inc., respectively. Guanidine hydrochloride (GdnHCl) was purchased from Invitrogen. Acetonitrile UV and acetontrile (HPLC grade) were purchased from Honeywell International, Inc. and EMD Chemicals, respectively. The tris-(2-carboxyethyl) phosphine (TCEP) hydrochloride solutions were freshly prepared before use. All other chemicals were purchased from Sigma.
All pH/pD values were measured using a Beckman
40 pH meter or ColorpHast pH indicator strips (pH 014; EM Science). pD values were computed by adding 0.4 to the corresponding pH readings. The standard buffer in this study contains 200 mM Na2SO4, 10 mM NaH2PO4, adjusted to pH 6.0. For the monomeric P
and P
peptides, the buffer also contained 1 mM reduced dithiothreitol to prevent dimerization of the monomers. No noticeable oxidation was detected after the experiment, as checked by reverse-phase HPLC.
Preparation of the peptide dimers
The P
and P
homodimers and the P
heterodimer were prepared as described elsewhere (Oas and Kim 1988). Briefly, roughly equal amounts of the P
and P
monomers were mixed in a solution containing 5 M GdnHCl and 200 mM Tris-HCl at pH 8.0 and were air-oxidized for 48 h at 25°C. The dimers were separated and purified on a Vydac C18 analytical column using a water/acetonitrile gradient in the presence of 0.1% trifluoroacetic acid (TFA). Each dimer was collected manually and dried under reduced pressure for further use.
Circular dichroism
CD spectra were obtained at 0°C on a Jasco J-810 CD spectropolarimeter (Jasco, Inc.) using a thermostatted 1-mm pathlength cell. The samples were dissolved in the standard buffer. The concentrations were determined by UV absorbance at 205 nm (Scopes 1974) and were typically 0.15 mg/mL.
Hydrogen/deuterium exchange (HDX)
The HDX experiments were carried out in a cold room, typically at 4°C. Each peptide sample was equilibrated with the standard buffer (1 nmol/µL) for about half an hour before each experiment. HDX was initiated by diluting 1 µL of aqueous protein solution with 19 µL of standard buffer made with D2O at pD 6.4. After 15 sec, 18 µL of aqueous TCEP solution (500 mM in 2M GdnHCl [pH 2.5], 4°C) were added to quench the exchange by decreasing the pH to 2.5 and to reduce the disulfide bond. The pH was maintained at 2.5 in all subsequent steps. Incubation of the sample with D2O for 15 sec should be long enough to exchange exposed amide hydrogens with deuteriums (Bai et al. 1993), but short enough to avoid significant exchange with H-bonded/protected amide hydrogens.
Analysis by MALDI-TOF mass spectrometry
Mass spectra were acquired on a PerSpective Biosystems Voyager DE STR instrument (PE Biosystems) in the positive ion reflectron mode of operation. The accelerating voltage was 20 kV; the grid voltage and guide wire voltage were 76% and 0.04% of the accelerating voltage, respectively. The mirror voltage ratio was 1.12, whereas the extraction delay time was 150 nsec. The matrix was
-cyano-4-hydroxycinnaminic acid in a 5 mg/mL solution containing acetonitrile, ethanol, and H2O (1:1:1, v:v:v); the final pH was adjusted to 2.5 with TFA. The measurement was done using next-spot external calibration with the P
and P
monomers as the calibrants, for which the calculated monoisotopic [M+H]+ masses are 1671.7219 Da and 1681.8160 Da, respectively. The average mass of each peptide was measured experimentally by determining the centroid of the distribution of its isotope peak intensities using Data Explorer TM 4.0. Ion current from 256 laser shots was accumulated for each spectrum.
Immediately after quenching the HDX process, the sample was desalted with cold water containing 0.1% TFA (pH 2.5) using a ZipTip (Millipore), eluted with matrix solution, and spotted on a chilled MALDI target. The target was immediately placed in a desiccator under moderate vacuum to dry the spot within 12 min. The target was immediately transferred to the MALDI instrument. A highly deuterated standard was prepared by dissolving the sample in D2O at pH 7.0 with 8 M D4-urea and incubating at room temperature for over 48 h. Treating the highly deuterated standard by the same protocol indicated that the extent of back-exchange during the experimental procedures was 35%. The uncertainty of the mass shift values is ±0.1 Da, expressed as the standard deviation of results from triplicate runs.
Analysis by CID tandem mass spectrometry
Mass spectra were acquired on a Finnigan LCQDeca ion-trap mass spectrometer (ThermoQuest) with an ion-spray voltage of 3 kV, a capillary voltage of 24 V, a tube lens offset of 55 V, and a capillary temperature of 150°C. Data-dependent MS/MS conditions were set with a default collision energy of 35%, a default charge state of 2, and an isolation width of 5 (m/z). Scans were taken over the range m/z 2351500. Typically, 100 scans were accumulated per spectrum.
Immediately after quenching the HDX process, the sample was desalted with cold water containing 0.1% formic acid (FA) (pH 2.5) using a ZipTip and eluted with 100 µL of a water/acetonitrile solution (1:1, v:v) preadjusted to pH 2.5 with FA.The processed sample was infused into the LCQDeca immediately with a prechilled (ice) 100-µL syringe at a flow rate of 20 µL/min; the syringe was kept buried in ice throughout the experiment. The extent of back-exchange during the infusion-CID-MS/MS procedures was 25%, as determined from analysis of the highly deuterated standard prepared as described above. Mass spectra were obtained with the operating software Xcalibur, and the centroid values were calculated with the Magtran software (Zhang and Marshall 1998).
Electronic supplemental material
Three items were provided as supplementary material: definitions of structural terms; Figure 1A
, ribbon diagram of bovine pancreatic trypsin inhibitor (BPTI), with peptides P
and P
highlighted in red and blue, respectively; and Figure 1B
, ribbon diagram of the natively structured P
peptide.
| Footnotes |
|---|
| Acknowledgments |
|---|
| References |
|---|
|
|
|---|
Cai, X. and Dass, C. 2005. Structural characterization of methionine and leucine enkephalins by hydrogen/deuterium exchange and electrospray ionization tandem mass spectrometry. Rapid Commun. Mass Spectrom. 19: 18.[CrossRef][Medline]
Carulla, N., Woodward, C., and Barany, G. 2000. Synthesis and characterization of a
-hairpin peptide that represents a "core module" of bovine pancreatic trypsin inhibitor (BPTI). Biochemistry 39: 79277937.[CrossRef][Medline]
. 2002. BetaCore, a designed water soluble four-stranded antiparallel
-sheet protein. Protein Sci. 11: 15391551.
Deng, Y., Pan, H., and Smith, D.L. 1999. Selective isotope labeling demonstrates that hydrogen exchange at individual peptide amide linkages can be determined by collision-induced dissociation mass spectrometry. J. Am. Chem. Soc. 121: 19661967.[CrossRef]
Englander, S.W. 2000. Protein folding intermediates and pathways studied by hydrogen exchange. Annu. Rev. Biophys. Biomol. Struct. 29: 213238.[CrossRef][Medline]
Goodman, E.M. and Kim, P.S. 1989. Folding of a peptide corresponding to the
-helix in bovine pancreatic trypsin inhibitor. Biochemistry 28: 43434347.[CrossRef][Medline]
Harrington, W.F. and Schellman, J.A. 1956. Evidence for the instability of hydrogen-bonded peptide structures in water, based on studies of ribonuclease and oxidized ribonuclease. C. R. Trav. Lab. Carlsberg [Chim.] 30: 2143.
Hoerner, J.K., Xiao, H., Dobo, A., and Kaltashov, I.A. 2004. Is there hydrogen scrambling in the gas phase? Energetic and structural determinants of proton mobility within protein ions. J. Am. Chem. Soc. 126: 77097717.[CrossRef][Medline]
Kaltashov, I.A. and Eyles, S.J. 2002. Studies of biomolecular conformations and conformational dynamics by mass spectrometry. Mass Spectrom. Rev. 21: 3771.[CrossRef][Medline]
Kim, M.Y., Maier, C.S., Reed, D.J., and Deinzer, M.L. 2001. Site-specific amide hydrogen/deuterium exchange in E. coli thioredoxins measured by electrospray ionization mass spectrometry. J. Am. Chem. Soc. 123: 98609866.[CrossRef][Medline]
Klotz, I.M. and Franzen, J.S. 1962. Hydrogen bonds between model peptide groups in solution. J. Am. Chem. Soc. 84: 34613466.[CrossRef]
Lewis, P.N., Go, N., Go, M., Kotelchuck, D., and Scheraga, H.A. 1970. Helix probability profiles of denatured proteins and their correlation with native structures. Proc. Natl. Acad. Sci. 65: 810815.
Montelione, G.T. and Scheraga, H.A. 1989. Formation of local structures in protein folding. Acc. Chem. Res. 22: 7076.[CrossRef]
Oas, T.G. and Kim, P.S. 1988. A peptide model of a protein folding intermediate. Nature 336: 4248.[CrossRef][Medline]
Scheraga, H.A. 1971. Theoretical and experimental studies of conformations of polypeptides. Chem. Rev. 71: 195217.[CrossRef][Medline]
Schulman, B.A. and Kim, P.S. 1994. Hydrogen exchange in BPTI variants that do not share a common disulfide bond. Protein Sci. 3: 22262232.[Abstract]
Scopes, R.K. 1974. Measurement of protein by spectrophotometry at 205 nm. Anal. Biochem. 59: 277282.[CrossRef][Medline]
Sreerama, N., Manning, M.C., Powers, M.E., Zhang, J.X., Goldenberg, D.P., and Woody, R.W. 1999. Tyrosine, phenylalanine, and disulfide contributions to the circular dichroism of proteins: Circular dichroism spectra of wild-type and mutant bovine pancreatic trypsin inhibitor. Biochemistry 38: 1081410822.[CrossRef][Medline]
Staley, J.P. and Kim, P.S. 1994. Formation of a native-like subdomain in a partially folded intermediate of bovine pancreatic trypsin inhibitor. Protein Sci. 3: 18221832.[Abstract]
Taniuchi, H. and Anfinsen, C.B. 1971. Simultaneous formation of two alternative enzymology active structures by complementation of two overlapping fragments of staphylococcal nuclease. J. Biol. Chem. 246: 22912301.
Weissman, J.S. and Kim, P.S. 1991. Reexamination of the folding of BPTI: Predominance of native intermediates. Science 253: 13861393.
Woody, R.W., Sugeta, H., and Kodama, T.S. 1996. Circular dichroism of proteins: Recent developments in analysis and prediction. Tanpakushitsu Kakusan Koso 41: 5669.[Medline]
Zhang, Z. and Marshall, A.G. 1998. A universal algorithm for fast and automated charge state deconvolution of electrospray mass-to-charge ratio spectra. J. Am. Soc. Mass Spectrom. 9: 225233.[CrossRef][Medline]
![]()
CiteULike
Connotea
Del.icio.us
Digg
Reddit
Technorati What's this?
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |