|
|
||||||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Methods |
Center for Experimental Therapeutics and Reperfusion Injury, Department of Anesthesiology, Perioperative and Pain Medicine, Brigham and Women's Hospital and Harvard Medical School, Boston, MA 02115
Published, JLR Papers in Press, February 1, 2005. DOI 10.1194/jlr.D400020-JLR200
The online version of this article (available at http://www.jlr.org) contains supplemental data. ![]()
1 Y. Lu and S. Hong contributed equally to this work. ![]()
2 To whom correspondence should be addressed. e-mail: cnserhan{at}zeus.bwh.harvard.edu
| ABSTRACT |
|---|
|
|
|---|
These results indicate that COCAD has many advantages for profiling and identification of LMs compared with the conventional dot-product algorithm.
Abbreviations: COCAD, cognoscitive-contrast-angle algorithm and databases; CRT, chromatographic retention time; ESI, electrospray ionization; HETE, hydroxy-eicosatetraenoic acid; LC-UV-MS/MS, liquid chromatography-ultraviolet-tandem mass spectrometry; LM, lipid mediator; LTB4, 5S,12R-dihydroxy-6Z,8E,10E,14Z-eicosatetraenoic acid; LXA4, 5S,6R,15S-trihydroxy-7E,9E,11Z,13E-trans-11-cis-eicosatetraenoic acid; LXB4, 5S,14R,15S-trihydroxy-6E,8Z,10E,12E-eicosatetraenoic acid; PGE2, 11
,15S- dihydroxy-9-oxo-prosta-5Z,13E-dien-1-oic acid; PGF2
, 9
,11
,15S-trihydroxy-prosta-5Z,13E-dien-1-oic acid
Supplementary key words cognoscitive-contrast-angle algorithm and database theoretical database and algorithm inflammation resolution eicosanoids resolvins neuroprotectins docosatrienes polyunsaturated fatty acid
| INTRODUCTION |
|---|
|
|
|---|
Although the physicochemical properties of individual LMs can give different chromatographic behaviors and mass spectra that are instrument- and chromatographic system-dependent, to facilitate profiling analyses between laboratories and assembling of system metabolomic maps, it is essential to first construct databases and search algorithms that permit accurate, prompt identification of LMs in complex biological matrices such as human blood, urine, and inflammatory exudates. To this end, one specific and sensitive tool for chemical identification in mediator-lipidomics is liquid chromatography-ultraviolet-tandem mass spectrometry (LC-UV-MS/MS) with atmospheric pressure ionization (710), especially electrospray ionization (ESI). Low-energy ionization of LMs with ESI generates primarily molecular (or pseudomolecular) ions for low/intermediate energy collision-induced dissociation (CID) MS/MS analysis and can avoid unwanted degradation products (7). MS/MS spectra are used extensively to identify and elucidate the structures of LMs derived from PUFAs (24, 79, 11). LC-UV-MS/MS can also provide more spectral information and criteria for those compounds with specific UV chromophores and assist in the elucidation of both known and novel LM structures. In this context, many LMs derived from PUFAs possess conjugated double bond systems that are critical components for their bioactions; each gives a characteristic UV spectrum (24, 79). Unesterified PUFAs frequently coexist in samples, and these compounds contain 1,4-cis-pentadiene subunits that do not give a specific UV chromophore signature. For example, leukotrienes contain a conjugated triene chromophore and lipoxins possess a conjugated tetraene chromophore.
Mediator-lipidomic databases containing MS/MS as well as UV spectra and chromatographic profiles with appropriate search algorithms are imperative for accurate and timely analysis of LMs (12). Mediator-lipidomic analysis is summarized in
Scheme 1.
|
Although GC-MS can be used in many cases for structure identification, the required high column temperatures limit its use because many of the PUFA-derived LMs are thermolabile (7). Search algorithms for GC-MS electron impact ionization spectra have been well studied (1315), as have those for both GC-MS chromatograms and mass spectra (16). The contrast-angle (or dot-product, = cosine of the contrast-angle) algorithm is widely used. The contrast-angle is the angle between two spectra represented as vectors composed of ordered peak intensities (13). Because databases and search algorithms for LMs on comprehensive LC-UV-MS/MS are not yet available, we initially used MassFrontierTM (ThermoFinnigan) GC-MS mass spectral commercial software to construct a LM database and search algorithm. The search algorithm for MassFrontier is dot-product and was developed by Stein and Scott (13, 14).
The dot-product algorithm in MassFrontier uses the intensities of all peaks in MS/MS as a vector for computation and is not concerned with the identities of ions, whether they are molecular ions, ions derived from molecular ions, or ions from interfering substances in the samples (1316). Using the ions generated from interfering background can decrease the percentage correct for best match of the search. Even the ions generated from unknown LMs contribute differently to the identification. For LMs, ions generated by cleavage of the carbon-carbon bond inside the carbon chain have greater diagnostic yield for determining the structures than the ions formed by the loss of water or CO2, which are common features of most PUFA-derived LMs (3, 4, 7, 8, 11).
The relationships between ESI-MS/MS spectra and LM structures have been studied extensively (24, 79). The identification and elucidation of the unknown structures of novel LMs by LC-UV-MS/MS are based on the relationships between structural features (such as functional groups and double bonds) and characteristics of the spectra and chromatograms (MS/MS ions, UV spectra, and LC retention times). On C18 reversed-phase LC, chromatographic retention times (CRTs) generally increase in the following order: trihydroxy LMs, dihydroxy LMs, monohydroxy LMs, and then the nonoxidized PUFA precursors (such as arachidonic acid) are eluted, depending on the mobile phase used. For positional isomers of hydroxy-containing LMs, the CRTs decrease as the hydroxy groups locate closer to the methyl group at the
-end of the long-chain PUFA. For instance, the CRT of 5S,14R,15S-trihydroxy-6E,8Z,10E,12E-eicosatetraenoic acid (LXB4) is usually shorter than that of 5S,6R,15S-trihydroxy-7E,9E,11Z,13E-trans-11-cis-eicosatetraenoic acid (LXA4). The CRT of native LXA4 (containing a 15S-hydroxy at carbon 15) is shorter than that of the aspirin-triggered 15-epi-LXA4 but greater than that of 11-trans-LXA4, which is the natural all-trans-containing isomer of LXA4 (17).
The band multiplicity and maximum absorbance wavelength (
max) of UV spectra are additional class- or series-specific signatures of LMs, such as the presence of an asymmetric singlet band with
max
235 nm for the conjugated diene present within the lipoxygenase pathways. These include monohydroxy-eicosatetraenoic acids (mono-HETEs); a triplet with
max
270 nm for the leukotrienes, such as the conjugated triene within LTB4 (5S,12R-dihydroxy-6Z,8E,10E,14Z-eicosatetraenoic acid); a triplet with
max
300 nm for conjugated tetraene, as present within LXA4 and LXB4; and, for example, an asymmetric singlet with
max
242 nm from the two conjugated dienes interrupted by a methylene group, as present in the double lipoxygenation product 5S,15S-di-HETE (17). For compounds without a class and/or series specific chromophore or with
max in vacuum UV range (actually undetectable on our LC-UV-MS/MS instrument) such as 1,4-cis-pentadiene-containing PUFAs as well as some prostaglandins, namely PGE2 (11
,15S-dihydroxy-9-oxo-prosta-5Z,13E-dien-1-oic acid) and PGF2
(9
,11
,15S-trihydroxy-prosta-5Z,13E-dien-1-oic acid), which do not possess conjugated double bond systems (17, 18), the UV chromophore component is not used for identification. When the conjugated tetraene of LXA4 isomers is in the all-trans geometry, the
max is shifted to 302 nm instead of 300 nm (17). The relationship between MS/MS spectra and LM structures is of interest and presented vide infra.
With current chemical analytical technologies, most LMs are identified manually by comparing the spectra and chromatographic behaviors acquired from sample tissues with those of authentic standards of known LMs; when authentic standards are not available, as in the case of novel LMs and their metabolites, basic chemical structures can be obtained on the basis of the relationship between structures and features of their spectra and chromatographic behaviors compared with those of synthetic and biogenic products prepared to assist in the assignment. We routinely identify LMs by comparing the unknown spectra (MS/MS, GC-MS, and UV spectra) and CRTs with those of authentic and synthetic standards if available or with a theoretical database that consists of virtual UV and MS/MS spectra and CRTs for discovering potentially novel LMs (3, 4, 12) if standards are not available. We initially developed a theoretical database and algorithm according to the relationships between LM structures and their spectral and chromatographic characteristics (12). The proposed structures of novel potential LMs in the theoretical databases were based on PUFA precursors and established biosynthetic pathways.
In this report, we focus on constructing mediator-lipidomic databases and search algorithms useful for the identification of LM structures using LC-UV-MS/MS with the following objectives: 1) to assemble a database using currently available mass spectral software; 2) to construct a cognoscitive-contrast-angle algorithm and databases (COCAD) to improve the identification of LMs using MS/MS ion identities that currently cannot be determined with available software; and 3) to develop a theoretical database and algorithm for assessing potentially novel and/or unknown structures of LMs and their metabolites in biologic matrices. It is quite meaningful to develop mediator-lipidomic databases and algorithms using ion trap mass spectrometers that are relatively cheaper and popular. Moreover, the fragmentation rules and patterns for CID spectra from triple-quadruple mass spectrometers, another popular MS instrument, are similar to those we encounter using the ion trap (7, 19, 20). The databases and algorithms described in this report should be a useful initiation that can be extended to other types of MS instruments.
| SYSTEMS AND METHODOLOGY |
|---|
|
|
|---|
We used a LUNA C18-2 (Phenomenex, Torrance, CA) column packed with 5 µm spherical C-18 particles. The dimension was 100 mm (or 150 mm) long and 2 mm inner diameter. The column was kept in a column heater (30°C). The LC system used a P-4000 quaternary LC pump (ThermoFinnigan). When the column was 100 mm long, the mobile phase (methanol-water-acetic acid, 65:35:0.01) flowed at 0.2 ml/min in C channel from 0 to 8 min, linearly ramped to methanol in D channel from 8.01 to 30 min, eluted as methanol for 5 min, and finally eluted as C again. When the column was 150 mm, the mobile phase (methanol-water-acetic acid, 65:35:0.01) flow rate was 0.2 ml/min through the C channel from 0 to 15 min, linearly ramped to methanol in D channel from 15.01 to 48 min, eluted as methanol for 8 min, and return eluted as C again.
UV spectra were acquired using a tandem photodiode-array UV detector (UV-6000; ThermoFinnigan) that scanned from 200 to 360 nm. Samples were injected manually with a 7125E injector equipped with a 100 µl loop (Rheodyne, Rohnert Park, CA). The injection volume was 20 µl, and the amount of each LM and related compounds in each set of analyses ranged from 0.5 to 100.0 ng.
The ion source conditions for the LCQ were set as follows: electrospray voltage, 4.3 kV; heating capillary, 230°C and 39 V; tube lens offset, 60 V; sheath N2 gas, 80 units (near 1.2 l/min); auxiliary N2 gas, 3 units (0.045 l/min). The ion trap analyzer typically scanned from m/z 200 to 800 in MS mode to determine the m/z of ESI-ion source-generated ions. The parameters for MS/MS were as follows: damping and collision gas, He (at
0.1 Pa); isolation width, m/z 1.5; collision time, 30 ms; normalized collision energy on base of the maximum available tickling voltage (5 V),
38% for prostaglandins and dihydroxy- and trihydroxy-containing acyclic eicosanoids, 42% for monohydroxy-containing eicosanoids, and 50% for PUFAs. The settings for data-dependent scan MS/MS were as follows: isolation width, m/z 1.5; normalized collision energy, 40%; minimum signal required, 100,000; minimum MSn signal required, 10,000; repeat count, 3; repeat duration, 0.5 min; exclusion duration, 1.5 min. Automated gain control was set at
2 x 107 for MS and MS/MS full scan.
Deuterium-labeled PGE2 (d4-PGE2), d4-LTB4, and d(n)-PUFA were routinely used in our laboratory as the internal standard to calculate recoveries during sample preparation and quantification. Deuterium-labeled eicosanoids were also routinely used to calibrate CRTs and account for potential daily variations in relative chromatographic behavior. CRTs for LMs were routinely within ± 0.6 min.
Generation of eicosanoids by polymorphonuclear leukocytes and stimulated spleen
Results from LC-UV-MS/MS for eicosanoids acquired from biological samples were used to demonstrate the performance of the databases and algorithms constructed. One sample set was obtained by incubations of rabbit polymorphonuclear leukocytes (
50 x 106 cells per incubation) isolated from rabbit whole blood and stimulated with ionophore A23187 (15 µM) and arachidonic acid (10 µM, 20 min, 37°C) as described (19). Samples were also analyzed from murine spleens (kindly provided by J. Aliberti and coworkers) prepared as described (20). The spleens were collected from mice injected with homogenate of Toxoplasma gondii (RH strain). Incubation of samples was stopped with two volumes of methanol, and samples were maintained at 20°C before isolation and further analyses. All samples were extracted using a solid-phase extraction technique (21). In short, samples rapidly preequilibrated at pH 3.5 were loaded onto solid-phase extraction cartridges (100 mg of C18 adsorbent per cartridge; Varian, Palo Alto, CA) and eluted with 8 ml of methyl formate (EM Science, Gibbstown, NJ). Samples were taken to dryness using a flow of nitrogen gas. Each sample was suspended in 20 µl of mobile phase C and injected for LC-UV-MS/MS analysis.
Software
For LMs and related compounds for which synthetic and authentic standards were available, the databases were constructed with MassFrontier (version 2.0; HighChem, Ltd., San Jose, CA) (22), a commercial software designed for use with positive-ion MS spectra generated by electron impact ionization. We implemented cognoscitive-contrast-angle and theoretical algorithms and developed databases using Microsoft Visual Studio 6.0 for mass spectra obtained in the negative-ion mode. All of the spectra generated were processed with Microsoft Excel 2000.
Tests to evaluate and optimize systems of LC-UV-MS/MS mediator-lipidomic databases and algorithms
Murine tissue was extracted along with different amounts of LMs, their precursors, and metabolites that were added for these experiments. Each of six different murine tissues (kidney, spleen, liver, brain, colon, and fat) was homogenized in cold methanol. Each homogenate was centrifuged (14,000 rpm, 20 min), and supernatants were spiked with 1, 2.5, or 5 ng of the following 34 LMs and the related compounds denoted in Table 1 (n = 23). The trihydroxy-containing eicosanoids LXA4, LXB4, 5S,12R,20-trihydroxy-6E,8Z,10E,14Z-eicosatetraenoic acid; the dihydroxy-containing eicosanoids 5,15-diHETE (5,15-dihydroxy-6E,8E,11Z, 13Z-eicosatetraenoic acid), LTB4, 5S,6S-diHETE (5S,6S-dihydroxy-7E,9E,11Z,14Z-eicosatetraenoic acid), 5S,6R-diHETE (5S,6R-dihydroxy-7E,9E,11Z,14Z-eicosatetraenoic acid), 12S,20-dihydroxyeicosa-5Z,8Z,10E,-14Z-tetraenoic acid; HETEs (monohydroxy eicosatetraenoic acids), 20-HETE (20-hydroxy-5Z,8Z,11Z,14Z-eicosatetraenoic acid), 15S-HETE (15S-hydroxy-5E,8Z,11Z,13Z-eicosatetraenoic acid), 12S-HETE (12S-hydroxy-5Z,8Z,10E,14Z-eicosatetraenoic acid), 11-HETE (11-hydroxy-5E,8Z,12Z,14Z-eicosatetraenoic acid), 8-HETE (8-hydroxy-5E,9Z,11Z,14Z-eicosatetraenoic acid), and 5S-HETE (5S-hydroxy-6E,8Z,11Z,14Z-eicosatetraenoic acid); the prostaglandin PGE2, 9-oxo-11
,15S-dihydroxy-prost-13E-en-1-oic acid, 9,15-dioxo-11
-hydroxy-prosta-5Z,13E-dien-1-oic acid, 9,15-dioxo-11a-hydroxy-prost-5Z-en-1-oic acid, 9
,15S,-dihydroxy-11-oxo-prosta-5Z,13E-dien-1-oic acid, PGF2
, 9,15-dioxo-11
-hydroxy-prosta-13E-en-1-oic acid, 6-oxo-9
,11
,15S-trihydroxy-prost-13-en-1-oic acid, 11-oxo-prosta-5Z,9,12E,14Z-tetraen-1-oic acid, 15S-hydroxy-11-oxo-prosta-5Z,9E,13E-trien-1-oic acid, and 9
-hydroxy-11-oxo-prosta-5Z,12E,14E-trien-1-oic acid; the precursors of LMs arachidonic acid (eicosa-5Z,8Z,11Z,14Z-tetraenoic acid), eicosapentaenoic acid [eicosa-5Z,8Z,11Z,14Z,17Z-pentaenoic acid (20:5,n-3)], docosahexaenoic acid [docosa-4Z,7Z,10Z,13Z,16Z,19Z-hexaenoic acid (22:6, n-3)], and linolenic acid [octadeca-9Z,12Z,15Z-trienoic acid (18:3, n-3)]; and others: 15-oxo-eicosa-5Z,8Z,11Z,13E-tetraenoic acid, 5S-hydroxy-6R-(S-glutathionyl)-7E,9E,11Z,14Z-eicosatetraenoic acid, 5S-hydroxy-6R-(S-cysteinylglycinyl)-7E,9E,11Z,14Z-eicosatetraenoic acid, 9
,11,15S-trihydroxythromba-5Z,13E-dien-1-oic acid, and hepoxilin A3 [8(R,S)-hydroxy-11S,12S-epoxy-eicosa-5Z,9E,14Z-trienoic acid]. Each sample was then concentrated by a stream of N2 gas, dissolved into 20 µl of mobile phase C, and injected into the instrument to acquire LC-UV-MS/MS chromatograms and spectra, which were then searched against the databases using the algorithms that we developed (vide infra). For each amount of one compound in six different tissues, one value of percentage correct for best match was obtained. The mean and SEM of percentage correct for best match were calculated from one compound group at one amount for each database/algorithm. The primary measure of performance in these experiments was the correct best match and compound identification.
|
| RESULTS AND DISCUSSION |
|---|
|
|
|---|
|
0.3 min. Thus, only the unknown CRT within ±2 x 0.3 min of the CRT (for 95% confidence intervals) in the databases was taken as a correct match. If the UV spectral pattern was unclear, the MS/MS and CRT were still searched to avoid potential errors in assignment. A standard LM or theoretical fragmentation/ion fragmentation pattern that fulfilled the above match criteria was then assigned to the unknown set. If the match was a "hit" only with UV and MS/MS spectra but not with CRT, the LM in the sample was likely to be a geometric isomer of a known LM.
UV patterns are simply presented in the database as
max and divided into six clusters:
301 nm,
278 nm,
270 nm,
242 nm,
235 nm, and vacuum UV. In contrast, each MS/MS spectrum has dozens of ions and therefore is a much larger data set than a UV
max. If the UV
max is matched first, then only the MS/MS subdatabase having LM standards with this UV
max, which is a fraction of the whole MS/MS database, needs to be searched. This type of narrowed search is faster than that through the entire MS/MS database. This will become more significant when more and more standard MS/MS data are entered into the database. We tested the route with MS/MS spectra before UV; the search results are the same as having MS/MS after UV data. If the UV pattern is unclear, we can just search the MS/MS and CRT to avoid the assignment errors.
Databases constructed with commercial mass spectral software containing LC-UV-MS/MS chromatograms and spectra measured from standard LMs
A mediator-lipidomic database using LC-UV-MS/MS spectra and chromatograms acquired from authentic LMs was constructed with GC-MS spectral software MassFrontier. The UV
max of standard LMs were written into the subdatabase names and the CRTs were written into the LM names, so that MassFrontier could handle the acquired UV spectral results and CRTs for the identification of the unknown LMs according to the logic diagram in Scheme 2.
The demonstration of the efficacy of this database was assessed with the identification of eicosanoids biosynthesized by rabbit leukocytes using LC-UV-MS/MS (Fig. 1)
. For peak I in the inset of Fig. 1A, a subdatabase "mTOz 351Vacuum UV" was selected with the molecular ion of interest and matched UV
max. The search showed that the MS/MS spectrum of peak I matched best with PGE2, with the highest matching score of 907 (Fig. 1B). Furthermore, the CRT of peak I and standard PGE2 matched. Therefore, peak I was assigned as PGE2, consistent with the manual identification and assessment of MS and UV spectral and LC chromatographic features of PGE2.
|
|
Typical chain-cut ions for LMs in MS/MS are formed by
-cleavage of the carbon-carbon bonds connecting to the carbon with a functional group directly attached (24, 7, 8). LMs readily undergo
-cleavage (7). We proposed the nomenclatures illustrated in the LXA4 structure presented in Scheme 3 to systematically name the segments formed via chain cut and chain plus peripheral cut without concern for hydrogen shifts occurring during mass spectrometric analysis of PUFA-derived products. Each of these LMs derived from PUFA has a carboxyl terminus and a methyl terminus. An Arabic number was used to designate the position of the functional group on the LM carbon chain where the cleavage occurs. The uppercase letter right after the number indicates the side of the functional group on which the cleavage occurs: C is for cleavage on the carboxyl side of the functional group, and M is for cleavage on the methyl side of the functional group. Each cleavage can directly generate two segments. The lowercase letter that follows indicates the side of the cleavage on which the segment forms: c is for a segment formed on the carboxyl side of the cleavage, and m is for a segment formed on the methyl side of the cleavage. Segments formed through ß-cleavage or
-cleavage toward the functional group are named by adding ß or
between the uppercase and lowercase letters. Because three hydroxy groups of LXA4 are located at C5, C6, and C15 of the 20 carbon linear chain,
-cleavages can occur on the bonds at C4-C5, C5-C6, C6-C7, C14-C15, and C15-C16 to generate 10 chain-cut segments (Scheme 3). All of the possible chain-cut, peripheral-cut, and chain-plus-peripheral-cut segments for LXA4 are listed in the table in Scheme 3.
A MS/MS ion detected from LM samples in the negative-ion mode generally is formed from a specific segment with the addition or subtraction of hydrogen(s) caused by hydrogen shift during the cleavage. The charge (z) of the LM negative ion is usually equal to 1; therefore, the mass-to-charge ratio (m/z) of a LM ion is usually equal to its mass (m). Previous reports (7, 8, 11) and our unpublished results (data not shown) demonstrated that the MS/MS fragmentation observes the empirical rules regarding the addition or subtraction of hydrogen(s) for the chain-cut segments to form the chain-cut MS/MS ions: for segment Cc, the detected MS/MS ions are Cc and Cc + H; for segment Cm, the detected MS/MS ions are Cm, Cm ± H, and Cm ± 2H; for segment Mc, the detected MS/MS ions are Mc and Mc H; and for segment Mm, the detected MS/MS ions are Mm, Mm ± H, and Mm ± 2H (referred to below as the fragmentation rules).
These rules were used to identify the segments for instrument-detected MS/MS ions. The interpreted ions, such as M H H2O, Cc, Cc + H, and Mc H, identified as the MS/MS ions detected in the mass spectrometer, are called virtual ions. One detected MS/MS ion can be interpreted as one or several virtual ions. Through the loss of H2O, CO2, NH3, and/or amino acids, the chain-cut ions can form chain-plus-peripheral-cut ions. For the chain-cut and chain-plus-peripheral-cut ions in the present report, we focused on those formed by
-cleavages. Those detected MS/MS ions uninterpretable via rule 1 and neutral loss are taken as unidentified ions.
Modification of MS/MS ion intensities according to their identities In nature, isotope 13C constitutes 1.1% of elemental carbon. If the intensity of an ion with mass M (ion charge z is 1) containing only carbons, hydrogens, and oxygen(s) is 100, the contribution of 13C to the intensity of the ion with mass (M + 1) (the ion with mass 1 unit higher) is statistically equal to 100 x 1.1% x the number of carbons in ion M, which is listed in Table 2.2 of the report by McLafferty and Turecek (23). The contribution of 13C to the ion (M + 2) is much less than 100 x 1.1% x carbon number. The 13C contribution of ion M to the intensities of ions (M + 1) and (M + 2) was subtracted in the computation of the matching score for the COCAD and theoretical algorithms. If ion M could not be identified, the carbon number in the ion (M + 1) or (M + 2) is used because M, (M + 1), and (M + 2) generally have the same carbon number.
Chain-cut ions are most informative and can be diagnostic for determining LM structures such as the position of functional groups and double bonds. Generally, there are more chain-plus-peripheral-cut ions than chain-cut ions, because the former are derived from the latter by one or several neutral losses. It is likely that some chain-plus-peripheral-cut ions have the same m/z even though they originate by different mechanisms and represent different structural features. Therefore, chain-cut ions are more specific than chain-plus-peripheral-cut ions for defining the LM structure. Peripheral-cut ions in MS/MS spectra are similar among LM isomers; therefore, they were not specific enough for differentiation of individual LM isomers.
According to the fragmentation rules, the nth MS/MS peak can be identified as one or several chain-cut (C) ions, peripheral-cut (P) ions, and/or chain-plus-peripheral-cut (CP) ions. The weighted intensity (yIn) of each identified ion is as follows:
![]() |
![]() |
![]() |
![]() |
represents the contribution of peripheral-cut ions to In' (
= 3 for peripheral-cut ions formed via the loss of one CO2 from each molecular ion,
= 10 for peripheral-cut ions formed via the loss of one H2O from each molecular ion, and
= 1 for other peripheral-cut ions formed via multiple loss of CO2 and/or H2O from each molecular ion). The assignment of
values is arbitrary and based on the observation of relative intensities of peripheral-cut ions in MS/MS spectra of LMs.
COCAD contrast angle
COCAD used a contrast-angle algorithm to match an MS/MS spectrum between sample and standards. The contrast angle is calculated as follows:
![]() |
![]() |
![]() |
![]() |
Uv is equal to Cv, CPv, or Pv for unknown spectrum to be identified; Sv is equal to Cv, CPv, or Pv for standard spectrum.
![]() |
![]() |
![]() |
![]() |
The COCAD contrast angle in equation 6 represents how well the spectrum of the sample matches the standard: if it is 0°, the two spectra match exactly; if it is 90°, the two spectra do not match at all. The smaller the contrast angle between 0° and 90°, the better the match (13, 24). The value is integrated and normalized from dot-products DC, DCP, and DP (equation 6). The numeric coefficient 10 in equation 6 was found to be the best value (2, 20, and 100 were also tested) that emphasizes the fingerprinting feature of chain-cut ions, because chain-cut ions are more important for determining the LM structure than are other types of ions. To normalize [(10 x DC + DCP + DP) ÷ (11 +
CP)] in equation 6 to be no more than 1, 11 was used in the denominator of equation 6, and
CP is equal to 1 if at least one MS/MS ion is identified as a chain-plus-peripheral-cut virtual ion or is equal to 0 if no such ion is identified. No chain-plus-peripheral-cut ion is identified in a few LM standard spectra. Therefore,
CP is introduced in equation 6 to normalize the COCAD contrast angle to 0 when matching these types of spectra against themselves. Unidentified ions were excluded for matching in equations 2 to 6.
According to the empirical fragmentation rules, for monohydroxy-containing LMs, the vth chain-cut ion can be Cc, Cc + 1, Cm 2, Cm 1, Cm, Cm + 1, Cm + 2, Mc 1, Mc, Mm 2, Mm 1, Mm, Mm + 1, or Mm + 2. The vth chain-plus-peripheral-cut ion can be Cc CO2, Cc CO2 + 1, Cm H2O 2, Cm H2O 1, Cm H2O, Cm H2O + 1, Cm H2O + 2, Mc H2O 1, Mc H2O, Mc CO2 1, Mc CO2, Mc H2O CO2 1, or Mc H2O CO2. The vth peripheral-cut ion can be M H CO2, M H H2O, or M H H2O CO2. For LMs having multiple functional groups, V is accordingly greater.
Integrating MS/MS spectra with UV spectra and chromatograms in COCAD to identify LMs
COCAD databases incorporate the LC-UV-MS/MS chromatograms and spectra measured from standard LMs and the segments formed by cleavage illustrated in Scheme 3. The features of UV spectra and chromatograms are represented with the UV
max and CRTs, respectively. The search of COCAD databases was conducted stepwise, as depicted in the logical diagram shown in Scheme 2, from UV
max, to MS/MS spectra, and then to CRTs. Using COCAD, we identified LXA4 in murine spleen stimulated during exposure to T. gondii (Fig. 2A, B)
. The search was narrowed to the subdatabase with UV
max 301 nm (Fig. 2A, middle inset) and molecular ion m/z 351, which were detected for the unknown peak II in the chromatogram (m/z 251 of MS/MS 351) for the stimulated spleen sample. The CRT (5.2 min) of peak II was then matched against the candidates with the highest MS/MS match scores (Fig. 2A, left inset). For comparison, searching was also conducted with the system developed with MassFrontier (Fig. 2C). Both demonstrated that LXA4 was the best match for the unknown peak II shown in Fig. 2A. However, COCAD and MassFrontier could not be used alone to identify the structures of unknowns without LC-UV-MS/MS chromatograms and spectra obtained for authentic and/or synthetic standards.
|
max, and CRTs predicted for potentially novel LMs. Searching against a theoretical database is also performed stepwise as described in Scheme 2, from UV
max, to MS/MS spectra, and then to CRTs. Equation 7 is the matching score for an unknown MS/MS spectrum compared with the virtual spectrum based on the segments and the fragmentation rules.
![]() |
The matching score in equation 7 summates the weighted intensities of all of the identified MS/MS peaks in the spectrum acquired from the sample. The numerator of the formula is composed of three parts:
![]() |
![]() |
![]() |
![]() |
-cleavage formed from the f th functional group and matched to the nth MS/MS peak. F is the total number of functional groups in one LM. f is counted from the carboxyl terminus of LM. For example, f is 1 for the 5-hydroxy, 2 for the 6-hydroxy, and 3 for the 15S-hydroxy group present in LXA4. F for LXA4 is 3 (Scheme 3).
![]() |
-cleavage. The smallest MS/MS ion detectable in ion trap tandem mass spectrometry is generally
m/z 95 for LMs with a molecular ion of
400 Da (25). To compensate for the bias caused by the inability to detect an MS/MS ion of m/z < 95, factors (f TC ÷ fAC)0.5, (f TCP ÷ fACP)0.5, and (TP ÷ AP)0.5 are used in equation 7. f TC (or f TCP) is the total number of chain-cut (or chain-plus-peripheral-cut) ions formed from the f th functional group via
-cleavage. Tp is the total number of peripheral-cut ions formed from one LM. fAC (fACP or AP) is the fraction of f TC (f TCP or TP) representing the ions within the MS/MS detecting range (m/z from 95 to the m/z of molecular ion). 20-HETE is a typical example. The ions formed from the segment 20Cm and 20Cm H2O of 20-HETE are m/z 31 and 13 according to the fragmentation rules, which are too small to be detected in ion trap MS/MS. Consequently, without the use of factors (f TC ÷ fAC)0.5 and (f TCP ÷ fACP)0.5, 20-HETE could have a lower matching score than other HETEs, even though the MS/MS spectrum is that of 20-HETE.
![]() |
The unknown MS/MS spectrum was assigned to the LM and/or potentially novel LM that has the highest matching score in addition to having the predicted UV spectral and chromatographic features (e.g.,
max and CRTs). This theoretical system was used successfully to identify 15-HETE in murine spleen (Fig. 3)
. Peak III displayed at CRT of 20.4 min on the chromatogram (at m/z 219 of MS/MS 319; left inset) and had a UV
max of 235 nm (Fig. 3A, middle inset). Therefore, the search on the theoretical system was narrowed to the subdatabase with molecular ion m/z 319, UV
max 235 nm, and CRT 21 min. In this case example, 15-HETE gave the highest matching score among all compounds in the subdatabase (Fig. 3B). The MS/MS peaks identified were annotated with the ion interpretation that also shows a fragmentation mechanism (Fig. 3A). Segments of 15-HETE that matched the MS/MS peaks according to the fragmentation rules are italicized (Fig. 3B). For comparison, MassFrontier was also used to identify peak III, which also identified it as 15-HETE (Fig. 3C).
|
Only the identified ions were used for the COCAD and theoretical databases. The test results in Table 1 demonstrate that it is not necessary to use all of the MS/MS ions to identify a LM or other compounds. Moreover, in some cases it is not obvious and cost-efficient to identify all MS/MS ions. Rigorous identification can be done using isotope labeling and derivatization, which may not be feasible or realistic for routine analyses of the biologically relevant samples. We can identify a LM with precision using one or several chain-cut ions per functional group along with peripheral-cut ions and chain-plus-peripheral-cut ions according to the MS/MS fragmentation rules. For example, we use at least four ions for monohydroxy-containing eicosanoids, six ions for dihydroxy-containing eicosanoids, and eight ions for trihydroxy-containing eicosanoids. Along these lines, John Beynon and coworkers originated the Eight Peak Index of Mass Spectra database for manual identification of compounds (as reviewed in 26). For different amounts of LMs added to the samples (1, 2.5, and 5 ng), the percentage identified correctly was not significantly different. Based on the test experiments using the theoretical database and algorithm, the average matching score for correct best match is 2.02 with a standard error of 0.71 (n = 212). Therefore, the threshold score for a 95% confidence interval of matching score for correct best match is 2.02 1.96 x 0.71 = 0.62.
The databases and search algorithms were developed on the basis of LC-UV-MS/MS data of LMs. The ion intensity patterns of MS/MS spectra generated from an ESI-triple quadrupole mass spectrometer are quite similar to those from an ESI-ion trap mass spectrometer, because the collision energy for both types of instruments is in the low-energy region [a few to 100 eV (laboratory kinetic energy of ions)] (7, 19, 20, 27). Therefore, the constants and algorithms described in this report may fit the CID spectra from triple quadrupole MS/MS without much modification. For high-collision-energy (102 to
103 eV) CID spectra generated via sector or spectrum time-of-flight/time-of-flight analyzer, the relative intensity patterns are quite different compared with the low-energy spectra, although many ions occurred for both energy situations (7, 8, 19, 20, 27); for example, the peripheral-cut ions are less abundant than the chain-cut ions. For ion trap and triple quadrupole MS/MS, peripheral-cut ions are more abundant than chain-cut ions. Our constants and algorithms give chain-cut ions more weight than peripheral-cut ions because chain-cut ions are more important to define LM structures. Therefore, they may still fit high-collision-energy CID spectra. Nevertheless, the constants and algorithms should be thoroughly tested and modified accordingly to fit other instruments.
The CRTs used in this set of experiments were obtained using specified chromatographic conditions (for a column of 100 mm length and some for 150 mm length; see Systems and Methodology for details) and because several fundamental issues were our initial focus. Hence, the databases and algorithms were programmed so that the new LC-UV-MS/MS data, including other chromatographic conditions, can be easily entered and used in the databases.
| CONCLUSIONS |
|---|
|
|
|---|
The three mediator-lipidomic systems described in this report were constructed with results from LC-UV-MS/MS obtained using specific LC conditions and parameters. Although these low-collision-energy CID spectra were quite similar to those acquired from triple quadruple mass spectrometers (7, 19, 20), the lipidomic systems need to be thoroughly tested and modified wherever necessary to be useful with results from triple quadruple mass spectrometers, not to mention those from intermediate- and high-collision-energy mass spectrometers. Theoretical databases and algorithms can be powerful when used appropriately to reduce significantly the tedious work that experienced biochemical researchers have to do when identifying the structures of novel LMs. However, they should consider other information or constraints, such as biosynthetic and metabolic pathways, stereochemistry, and bioactions, to ensure that the results are correct. Moreover, new LC results can be entered into the systems to extend the mediator-lipidomic systems for use with other LC conditions. Copies of the program and libraries can be made available to investigators on request.
| ACKNOWLEDGMENTS |
|---|
Manuscript received August 12, 2004 and in revised form January 12, 2005.
| REFERENCES |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
Y.-P. Sun, S. F. Oh, J. Uddin, R. Yang, K. Gotlinger, E. Campbell, S. P. Colgan, N. A. Petasis, and C. N. Serhan Resolvin D1 and Its Aspirin-triggered 17R Epimer: STEREOCHEMICAL ASSIGNMENTS, ANTI-INFLAMMATORY PROPERTIES, AND ENZYMATIC INACTIVATION J. Biol. Chem., March 30, 2007; 282(13): 9323 - 9334. [Abstract] |