METHOD FOR PREDICTING ATHEROSCLEROTIC RISK
BACKGROUND OF THE INVENTION
Cross-Reference to Related Application
This application is a continuation-in-part of of U.S. Patent Application Serial Number 07/454,069 filed December 21 1989.
Statement Regarding Federally Sponsored Research
Funding for work described herein was provided by the Federal Government under a grant from the Department of Healt and Human Services. The Government may have certain rights in this invention.
Field of the Invention
The present invention relates to a method for predicting atherosclerotic risk in a living patient.
Prior Art
The plasma lipoproteins are complexes in which the lipids and proteins occur in a relatively fixed ratio. They carry water-insoluble lipids, such as cholesterol and cholesterol esters for eventual cellular utilization between various organs via the blood, in a form with a relatively small and constant
particle diameter and weight. While all cells require cholesterol for growth, excess accumulation of cholesterol by cells is known to result in a disease state referred to as atherosclerosis. It is also known that total serum cholesterol can be correlated with the incidence of atherosclerosis.
Human plasma lipoproteins occur in four major classes that differ in density as well as particle size as shown in the table below.
They are physically distinguished by their relative rates of flotation in high gravitational fields in the ultracentrifuge. All four lipoprotein classes have densities less than 1.21 g ml-1, whereas the other plasma proteins, such as albumin and -globulin, have densities in the range of 1.33 to 1.35 g ml -1 The characteristic flotation rates in Svedberg flotation units (Sf) of the lipoproteins are determined in an NaCl medium of density 1.063 g ml- 1at 26°C, in which lipoproteins float upward and simple proteins sediment.
The plasma lipoproteins contain varying proportions of protein and different types of lipid. The very low-density lipoproteins contain four different types of polypeptide chains having distinctive amino acid sequences. The high-density lipoproteins have two different types of polypeptide chains, of molecular weight 17,500 and 28,000. The polypeptide chains of the plasma lipoproteins are believed to be arranged on the surface of the molecules, thus conferring hydrophilic properties. However, in the very low-density lipoproteins and chylomicrons, there is insufficient protein to cover the surface; presumably the polar heads of the phospholipid components also contribute hydrophilic groups on the surface, with the nonpolar triacylglycerols in the interior. Biochemistry, Lehninger, Worth Publishers, Inc., New York, 1975, p.301.
The different lipoprotein classes contain varying amounts of cholesterol. A total serum cholesterol measurement is an average of the amount that each lipoprotein class contributes to the total serum lipoprotein.
It has long been hypothesized that the concentration of specific lipoprotein classes are predominantly responsible for the development of atherosclerosis. Studies indicate that low density lipoproteins are responsible for the accumulation of cholesterol in cells and high density lipoproteins are responsible for removing excess cholesterol from cells.
Notwithstanding the desirability of isolating the various lipoprotein classes in blood plasma; a fast, accurate, and reproducible technique for use in clinical laboratories has not heretofore been developed. The method most often used relies on ultracentrifugation techniques, which are time consuming and expensive.
Bell et al, FEBS LETTERS, Vol. 219, no. 1, July 1987; reported single-pulse and Hahn spin-echo 500 MHz H-1 NMR spectra of human blood plasma and isolated chylomicrons, very low density lipoproteins (VLDL), low density lipoproteins (LDL) and high density lipoproteins (HDL). They made specific assignments for the resonances of individual lipoproteins in the CH2 and CH3 (fatty acid), and NMe3 + (phospholipid choline head group) regions of the spectra of plasma (0.8-1.3 and 3.25ppm, respectively). Measurements were not made on whole plasma and no attempt was made to determine the concentration of HDL or LDL or to associate them with atherosclerotic risk.
Thus, a need exists for a fast, accurate, and reproducible method to predict atherosclerotic risk in living patients.
SUMMARY OF THE INVENTION
The present invention is a method for predicting atherosclerotic risk in a living patient. In accordance with the present invention, a sample of a patient's bodily fluid is subj ected to proton nuclear magnetic spectroscopy to generate a
water-suppressed proton nuclear magnetic resonance spectrum. A curve resolution procedure is performed on the methyl or methylene resonance envelope to allow classifying of the resonances as resulting from very low density lipoproteins, low density lipoproteins, high density lipoproteins and chylomicrons. The amplitude of the high density lipoprotein resonance is measured and a mean normalized value calculated. The measured amplitude is then classified into a normal or abnormal category compared to a predetermined standard for which an abnormal spectrum indicates a high atherosclerotic risk.
The VLDL, LDL, HDL and chylomicron-cholesterol values are obtained using standard purified lipoprotein preparations. Standard curves are prepared in which peak height is plotted versus concentration for the methyl or methylene resonance as well as peak area versus concentration for the methyl or methylene resonance. Test sample results are then compared to the standard curves to obtain the lipoprotein concentrations in the test sample. The triglyceride level is determined by adding the peak areas of the chylomicron and VLDL components. A risk index is determined using the peak heights and peak areas of the resolved components.
Accordingly, an object of the present invention is to provide a method for predicting atherosclerotic risk in a living patient.
Another object of the present invention is to provide the concentrations, expressed in mg/dl, of the VLDL, LDL and HDL components as well as the chylomicron and triglyceride levels.
Other objects and advantages of the invention will become apparent from the description of the drawings and the invention, which follow.
BRIEF DESCRIPTION OF THE DRAWINGS
FIG. 1 is a typical 360 MHz NMR spectrum of the non-water components (water suppressed) of a plasma sample from a healthy control obtained in accordance with the present invention;
FIG. 2 is an NMR spectrum for the same plasma sample from which the spectrum of Fig. 1 was obtained, using the same equipment and pulse frequency, except without water suppression;
FIG. 3 is an expanded view of the methyl and methylene region of the Fig. 1 sample;
FIG. 4 is an NMR spectrum of the methyl and methylene regions of a plasma sample of a patient with a high atherosclerotic risk;
FIG. 5 schematically illustrates the apparatus used to perform the method of the present invention;
FIG. 6 shows the results of a study performed using the method of the present invention; and
FIG. 7 shows an NMR spectrum of the methylene region and the four lipoprotein components as obtained using the curve resolution procedure.
DESCRIPTION OF THE PREFERRED EMBODIMENTS
At the outset, the invention is described in its broadest overall aspects with a more detailed description following.
The present invention is a method for predicting atherosclerotic risk in a living patient.
The test for high atherosclerotic risk, which is an object of the present invention, will typically be performed in vitro. The process operates on any lipid-containing body fluid, blood, or bone marrow plasma. Whole blood, serum, or plasma may be used. While the test may be performed on any such lipid-containing body fluid, work to date has focused on blood plasma and thus in a preferred embodiment of the present invention plasma is used. The test sample need not be fasting.
In accordance with the present invention, a sample of patient's whole blood or plasma is subjected to proton nuclear magnetic resonance spectroscopy to generate a nuclear magnetic
resonance spectrum. Since components of the NMR spectrum which have significant predictive value may be masked by other materials, such as the water signal, steps are taken to eliminate the masking, thereby producing an informative NMR spectrum.
A curve resolution procedure is performed on the methyl or methylene resonance envelope to allow classifying of the resonances as resulting from very low density lipoproteins, low density lipoproteins, high density lipoproteins and chylomicrons. FIG. 7 shows a methylene resonance NMR spectrum 1 before the curve resolution procedure and the spectrum as resolved into the four components; chylomicrons 2, VLDL 3, LDL 4 and HDL 5. The amplitude of the high density lipoprotein resonance, the amplitude of the low density lipoprotein resonance and the amplitude of the very low density lipoprotein are then measured and a mean normalized value calculated for each. The measured amplitudes are then classified into a normal or an abnormal category compared to a predetermined standard for which an abnormal spectrum indicates a high atherosclerotic risk, i.e., high HDL is a sign of low risk, high LDL is a sign of high risk.
In a preferred embodiment of the present invention, a water-suppressed proton NMR spectrum of a preprandial or postprandial plasma sample (or serum sample) is obtained and the chemical shift values of VLDL, LDL, HDL and chylomicrons in
the methyl or methylene resonance are identified using a curve resolution program. The area, line-width, and height of the VLDL, LDL, HDL and chylomicron components are then determined. The heights of the peaks are measured with a ruler or computer from the center of the baseline noise to the top of the peak.
Next, an LDL-cholesterol value is obtained using a standard purified LDL preparation. Standard curves are prepared in which peak height is plotted versus concentration for the methyl or methylene resonance as well as peak area versus concentration for the methyl or methylene resonance. Test sample results are then compared to the standard curves to obtain the LDL-cholesterol concentration in the test sample.
An HDL-cholesterol value is then obtained using a standard purified HDL preparation. Standard curves are prepared in which peak height is plotted versus concentration for the methyl or methylene resonance as well as peak area versus concentration for the methyl or methylene resonance. Test sample results are then compared to the standard curves to obtain the HDL-cholesterol concentration in the test sample.
A VLDL-cholesterol value is then obtained using a standard purified HDL preparation. Standard curves are prepared in which peak height is plotted versus concentration for the methyl or methylene resonance as well as peak area versus
concentration for the methyl or methylene resonance. Test sample results are then compared to the standard curves to obtain the VLDL-cholesterol concentration in the test sample.
A risk index is then determined using the resolved components of either the methyl or methylene resonance using one of the following equations:
e) determine Height HDL relative normalized intensity using CH3 resonances.
In a preferred embodiment, values for the risk index of 3.9 + 4.5 at a proton frequency of 360 MHz (8.45T) or 400 MHz (9.40T) indicates a normal atherosclerotic risk or a healthy/normal individual, and a value of 11.0 ± 4.1 indicates a high atherosclerotic risk.
In addition, the concentration of VLDL, LDL, HDL chylomicrons in mg/dl is determined. Values for each lipoprotein are obtained using a standard purified lipoprotein preparation. Standard curves are prepared in which the known
lipoprotein concentration is plotted against peak area for each lipoprotein fraction using either the methyl or methylene resonance. Test sample results are then compared to the standard curves to determine the concentration in mg/dl of each lipoprotein fraction; VLDL, LDL, HDL and chylomicrons. The triglyceride concentration is then determined by adding the peak area of the chylomicron component to the peak area of the VLDL component. The peak area is measured by a standard mathematical method using a computer program to carry out the calculations.
FIG. 1 shows water suppressed proton spectrum of a healthy control, and FIG. 2 shows a proton spectrum of the same sample without water suppression. The truncated resonance line of water is denoted A in FIG. 2. The resonance lines between 2 and 3 ppm (part per million of resonance frequency) arise from the methyl and methylene groups of the lipoprotein lipids. An expanded view of this region of the proton spectrum is shown in FIG. 3 for a normal control. FIG. 4 shows an NMR spectrum of the methyl and methylene regions of a plasma sample for a patient with a high atherosclerotic risk.
In a preferred embodiment of this invention, proton NMR spectroscopy is performed on a human blood plasma sample with the water signal is suppressed. The water suppressed proton NMR spectrum obtained is dominated by resonances of plasma lipoprotein lipids. Without water suppression, these non-water resonances are virtually overwhelmed by the water. Signal
averaging allows observation of resonances associated with non-water body fluid components, at high magnetic fields, even in the presence of water resonance. However, modern NMR spectrometers can almost completely suppress the water proton resonance. The water suppressed proton NMR spectrum of plasma is essentially plasma lipoproteins and a few low molecular weight molecules . The plasma protein protons are obscured because they comprise a broad smear of unresolved resonances. The sharper resonances of the more mobile lipoprotein protons are superimposed on this broad background.
Accordingly, in its preferred embodiment, the present invention uses one of a number of conventional water suppression techniques, i.e., techniques for suppression of the water proton NMR signal. Numerous techniques have been devised to suppress the water proton NMR signal in other contexts. These may be broadly divided into two categories: (1) those that attempt avoiding excitement of the water proton signal, e.g., rapid scan correlation spectroscopy and the selective excitation technique, and (2) those that arrange for the water proton magnetization to be extremely small at the time the observed radio frequency (rf) pulse is applied, e.g., the inversion recovery technique and saturation. These and other solvent suppression techniques are described by P.J. Hore in "Solvent Suppression in Fourier Transform Nuclear Magnetic Resonance", Journal of Magnetic Resonance, 55:283-300 (1983) and the reference footnoted therein.
Although the water suppression technique is preferred when using conventional NMR apparatus due to their inability to distinguish between the signal of the solvent protons and those of the moiety or species of interest, a sufficiently sensitive apparatus would eliminate the need for water suppression altogether.
Any conventional modern NMR spectrometer may be used in the practice of the present invention. In the preferred embodiments, however, an NMR spectrometer with a magnet at constant field strength is used. The NMR signal is Fourier transformed and the mean normalized HDL resonance amplitude for proton resonances of methyl or methylene groups are the NMR parameters of interest.
FIG. 5 illustrates a nuclear magnetic resonance (NMR) spectrometer 2 which is capable of performing proton NMR spectroscopy and which is preferably, but not necessarily, of the type that suppresses the NMR signal of water. The spectrometer 2 is adapted for examination of a sample 4 which is human blood plasma contained within a test tube 6. The spectrometer 2 contains means 8 for selecting the HDL NMR resonance line in the NMR spectrum of the sample 4 and measuring the mean normalized amplitude of the resonance so selected.
The spectrometer 2 also is of conventional construction and includes, in addition to all its other structure, a means 10 for storing a value or range of values. In the preferred embodiment, the mean normalized HDL resonance amplitude is compared with a value or range of values which represents the value or range of values to be expected from normal patients, i.e., patients who are not at a high risk of developing atherosclerosis. In accordance with the invention, the spectrometer 2 also includes means 12 for classifying the measured HDL resonance amplitudes as normal or abnormal, i.e., high atherosclerotic risk, based upon the stored information. This may be done by comparison, subtraction, or any other appropriate mathematical operation.
In the preferred embodiment, the selecting and measuring means 8 is pre-adjusted to measure the mean normalized HDL resonance amplitude. This may include suppressing the water signal from the NMR spectrum of the sample 4, or may alternatively be done directly where the spectrometer 2 is sensitive enough to do so.
Typical spectrometers that can perform the method of the present invention are the Bruker AM-360 and the Bruker AM-500. Of course, others skilled in the art will know of similar equipment to perform the method of the present invention.
Correct sample preparation and execution is essential to carry out a successful measurement on plasma. Blood is collected in tubes containing 70 ul of a solution of 15% Na2- EDTA and is maintained at 4°C until centrifugation. Plasma is separated and stored at 4°C until NMR analysis. Plasma samples are never frozen because freezing destroys lipoprotein lipid structural integrity. Samples which show any visible sign of hemolysis are excluded.
The spectra are obtained at preferably 20-22°C and most preferably at 21° C. A relatively broad range of proton frequencies may be employed, e.g., 60 MHz and higher, however, 360 MHz is the most preferred frequency. If cost is not a factor, 500 MHz may be the preferred frequency. The samples are shimmed individually on the area of the proton free induction decay until the full width at half height of the water resonance is 4 Hz of less. Careful shimming is of course an assumed component of good NMR laboratory technique.
The present invention is further illustrated by the following nonlimiting examples.
Example 1
The method of the present invention was applied to a group of 17 patients undergoing traditional lipid profile analysis. Blood was collected in non-siliconized vacutainer tubes
containing 70 ul of a solution of 15% Na2EDTA and maintained a 4°C until centrifugation. Plasma was separated and stored at 4º C until NMR analysis. Plasma samples were not frozen because freezing destroys lipoprotein lipid structural integrity. Samples which showed any visible sign of hemolysis were excluded.
All spectra were obtained at 21ºC using an 8.45 T Bruker AM spectrometer operating at 360 MHz for proton (H-1) NMR. Al studies were carried out in 5 mm OD sample tubes (Wilmad, Vineland, New Jersey; #507PP or #528PP). Each sample, containing 0.6 ml plasma, was shimmed individually on the area of the proton free induction decay (FID) until the full-width at half-height (FWHH) of the water resonance was 4 Hz or less. An internal quality control was found in the linewidth of the EDTA resonances. If all was well with the sample preparation and shimming, the linewidth (FWHH) of the EDTA resonances (without exponential broadening) had to be 2 Hz or less and was often between 1.0 - 1.5 Hz. In order to accomplish this, most H-1 probes require detuning to avoid radiation damping. The probe was detuned until the 90° radio-frequency pulse became 20 msec. In the 8.45 T spectrometer, this resulted in probe detuning of about 2 MHz. The sample was spun during shimming of the Z shim coils and during data acquisition. Our H-1 spectra were acquired using presaturation to suppress water and an inversion-recovery sequence to null any lactate methyl protons present. The presaturation pulse was 4.0 sec, with a delay of about 0.8 sec between the 180° and 90° pulse. Eight
FIDs were signal averaged and then Fourier transformed following multiplication by an exponential resulting in 2 Hz line-broadening. The portion of the spectrum from 0.5 to 1.6 ppm was phased so that the baseline level at the edges of the plot was the same. This resulted in defective phasing of other (non-plotted) portions of the spectra.
Results of the H-1 NMR spectroscopic evaluation are shown in FIG. 6. These results were obtained by analyzing the data using the equation:
risk value = Heights (LDL + HDL)
Height HDL
Example 2
The method of the present invention was applied to a plasma sample. Blood was collected in non-siliconized vacutainer tubes containing 70 ul of a solution of 15% Na2EDTA and maintained at 4°C until centrifugation. Plasma was separated and stored at 4º C until NMR analysis. Plasma samples were not frozen because freezing destroys lipoprotein lipid structural integrity. Samples which showed any visible sign of hemolysis were excluded.
All water-suppressed proton NMR spectra were obtained at 16-18°C using an 8.45 T Bruker AM spectrometer operating at 360 MHz with zero-filling to 16K or 32K data points. All studies were carried out in 5 mm OD sample tubes (Wilmad, Vineland, New
Jersey; #507PP or #528PP). Each sample, containing 0.6 ml plasma, was shimmed individually on the area of the proton free induction decay (FID) until the full-width at half-height (FWHH) of the water resonance was 4 Hz or less.
An internal quality control was found in the linewidth of the EDTA resonances. If all was well with the sample preparation and shimming, the linewidth (FWHH) of the EDTA resonances (without exponential broadening) had to be 2 Hz or less and was often between 1.0 - 1.5 Hz. In order to accomplish this, most H-1 probes require detuning to avoid radiation damping. The probe was detuned until the 90° radio-frequency pulse became 20 msec. In the 8.45 T spectrometer, this resulted in probe detuning of about 2 MHz. The sample was spun during shimming of the Z shim coils and during data acquisition. Our H-1 spectra were acquired using presaturation to suppress water and an inversion-recovery sequence to null any lactate methyl protons present. The presaturation pulse was 4.0 sec, with a delay of about 0.8 sec between the 180° and 90° pulse.
Eight FIDs were signal averaged and then Fourier transformed following multiplication by an exponential resulting in 2 Hz line-broadening. The portion of the spectrum from 0.5 to 1.6 ppm was phased so that the baseline level at the edges of the plot was the same. This resulted in defective phasing of other (non-plotted) portions of the spectra.
The chemical shift values of VLDL, LDL, HDL and chylomicrons in the methyl and methylene resonances were identified. A curve resolution program, such as NMRtwo from NMRi of Syracuse, NY, was used to resolve the component residues using initial spectrum and chemical shift values of the methyl and methylene resonances of the VLDL, LDL, HDL and chylomicrons. The area, line-width, and height of each resolved component was obtained.
Next, an LDL-cholesterol value was obtained using a standard purified LDL preparation. Standard curves were prepared in which peak height was plotted versus concentration for the methyl and methylene resonances as well as peak area versus concentration for the methyl and methylene resonances. Test sample results were then compared to the standard curves to obtain the LDL-cholesterol concentration in the test sample.
An HDL-cholesterol value was then obtained using a standard purified HDL preparation. Standard curves were prepared in which peak height was plotted versus concentration for the methyl and methylene resonances as well as peak area versus concentration for the methyl and methylene resonances. Test sample results were then compared to the standard curves to obtain the HDL-cholesterol concentration in the test sample .
A VLDL-cholesterol value was obtained using a standard purified HDL preparation. Standard curves were prepared in which peak height is plotted versus concentration for the
methyl or methylene resonance as well as peak area versus concentration for the methyl or methylene resonance. Test sample results were then compared to the standard curves to obtain the VLDL-cholesterol concentration in the test sample.
A risk index was then determined using the resolved components of either the methyl or methylene resonance using one of the following equations:
a) Heights (LDL + HDL)
Height HDL b) Areas (VLDL + HDL + LDL)
Area HDL c) Heights (VLDL + HDL + LDL)
Height HDL d) Area (HDL + LDL)
Area LDL e) determine Height HDL relative normalized intensity using CH3 resonances.
The concentrations of VLDL, LDL, HDL and chylomicrons, expressed in mg/dl, was determined by measuring the peak area of each resolved component using the methylene or methyl resonance spectrum. Standard curves were prepared where known concentrations of the lipoprotein fractions were plotted against the peak area and these were compared to the test results to determine the concentration of VLDL, LDL, HDL and chylomicrons in the test sample. The triglyceride concentration was then determined by adding the peak area of the chylomicron
component to the peak area of the VLDL component. The peak area was measured by a standard mathematical method using a computer program to carry out the calculations.
The invention may be embodied in other specified forms without departing from the spirit of essential characteristics thereof. The present embodiments are therefore to be considered in all respects as illustrative and not restrictive, the scope of the invention being indicated by the appended claims rather than by the foregoing description, and all changes which come within the meaning and range or equivalency of the claims are therefore intended to be embraced therein.