EP3559658B1 - Détermination automatique de temps de rétention prévu et de fenêtre de temps de rétention optimale prévue - Google Patents

Détermination automatique de temps de rétention prévu et de fenêtre de temps de rétention optimale prévue Download PDF

Info

Publication number
EP3559658B1
EP3559658B1 EP17885088.9A EP17885088A EP3559658B1 EP 3559658 B1 EP3559658 B1 EP 3559658B1 EP 17885088 A EP17885088 A EP 17885088A EP 3559658 B1 EP3559658 B1 EP 3559658B1
Authority
EP
European Patent Office
Prior art keywords
interest
compound
samples
xic
sample
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
EP17885088.9A
Other languages
German (de)
English (en)
Other versions
EP3559658A1 (fr
EP3559658A4 (fr
Inventor
Gordana Ivosev
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
DH Technologies Development Pte Ltd
Original Assignee
DH Technologies Development Pte Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by DH Technologies Development Pte Ltd filed Critical DH Technologies Development Pte Ltd
Publication of EP3559658A1 publication Critical patent/EP3559658A1/fr
Publication of EP3559658A4 publication Critical patent/EP3559658A4/fr
Application granted granted Critical
Publication of EP3559658B1 publication Critical patent/EP3559658B1/fr
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N30/00—Investigating or analysing materials by separation into components using adsorption, absorption or similar phenomena or using ion-exchange, e.g. chromatography or field flow fractionation
    • G01N30/02—Column chromatography
    • G01N30/86—Signal analysis
    • G01N30/8624—Detection of slopes or peaks; baseline correction
    • G01N30/8644—Data segmentation, e.g. time windows
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N30/00—Investigating or analysing materials by separation into components using adsorption, absorption or similar phenomena or using ion-exchange, e.g. chromatography or field flow fractionation
    • G01N30/02—Column chromatography
    • G01N30/86—Signal analysis
    • G01N30/8675—Evaluation, i.e. decoding of the signal into analytical information
    • G01N30/8679—Target compound analysis, i.e. whereby a limited number of peaks is analysed
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N30/00—Investigating or analysing materials by separation into components using adsorption, absorption or similar phenomena or using ion-exchange, e.g. chromatography or field flow fractionation
    • G01N30/02—Column chromatography
    • G01N30/62—Detectors specially adapted therefor
    • G01N30/72—Mass spectrometers
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N30/00—Investigating or analysing materials by separation into components using adsorption, absorption or similar phenomena or using ion-exchange, e.g. chromatography or field flow fractionation
    • G01N30/02—Column chromatography
    • G01N30/86—Signal analysis
    • G01N30/8624—Detection of slopes or peaks; baseline correction
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N30/00—Investigating or analysing materials by separation into components using adsorption, absorption or similar phenomena or using ion-exchange, e.g. chromatography or field flow fractionation
    • G01N30/02—Column chromatography
    • G01N30/86—Signal analysis
    • G01N30/8693—Models, e.g. prediction of retention times, method development and validation
    • H—ELECTRICITY
    • H01—ELECTRIC ELEMENTS
    • H01J—ELECTRIC DISCHARGE TUBES OR DISCHARGE LAMPS
    • H01J49/00—Particle spectrometers or separator tubes
    • H01J49/0027—Methods for using particle spectrometers
    • H01J49/0036—Step by step routines describing the handling of the data generated during a measurement
    • H—ELECTRICITY
    • H01—ELECTRIC ELEMENTS
    • H01J—ELECTRIC DISCHARGE TUBES OR DISCHARGE LAMPS
    • H01J49/00—Particle spectrometers or separator tubes
    • H01J49/004—Combinations of spectrometers, tandem spectrometers, e.g. MS/MS, MSn
    • H01J49/0045—Combinations of spectrometers, tandem spectrometers, e.g. MS/MS, MSn characterised by the fragmentation or other specific reaction
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N30/00—Investigating or analysing materials by separation into components using adsorption, absorption or similar phenomena or using ion-exchange, e.g. chromatography or field flow fractionation
    • G01N30/02—Column chromatography
    • G01N30/04—Preparation or injection of sample to be analysed
    • G01N2030/042—Standards
    • G01N2030/047—Standards external
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N30/00—Investigating or analysing materials by separation into components using adsorption, absorption or similar phenomena or using ion-exchange, e.g. chromatography or field flow fractionation
    • G01N30/02—Column chromatography
    • G01N30/88—Integrated analysis systems specially adapted therefor, not covered by a single one of the groups G01N30/04 - G01N30/86

Definitions

  • the teachings herein relate to determining the expected retention time and the optimal retention time window of the chromatographic peak for a compound of interest from standard samples or from monitoring a large number of other compounds of interest in addition to the compound of interest in a large number of samples. More particularly the teachings herein relate to systems and methods for automatically identifying the actual chromatographic peak of the compound of interest from standard samples or from monitoring a large number of other compounds of interest in addition to the compound of interest in a large number of samples. Once the actual chromatographic peak of the compound of interest is found, the expected retention time and the optimal retention time window can be computed.
  • Mass spectrometry is an analytical technique for detection and quantification of chemical compounds based on the analysis of mass-to-charge ratios (m/z) of ions formed from those compounds.
  • MS mass-to-charge ratios
  • LC liquid chromatography
  • a fluid sample under analysis is passed through a column filled with a solid adsorbent material (typically in the form of small solid particles, e.g., silica). Due to slightly different interactions of components of the mixture with the solid adsorbent material (typically referred to as the stationary phase), the different components can have different transit (elution) times through the packed column, resulting in separation of the various components.
  • a solid adsorbent material typically in the form of small solid particles, e.g., silica
  • the effluent exiting the LC column can be continuously subjected to mass spectrometric analysis to generate an extracted ion chromatogram (XIC), which can depict detected ion intensity (a measure of the number of detected ions, total ion intensity or of one or more particular analytes) as a function of retention time.
  • XIC extracted ion chromatogram
  • the LC effluents can be subjected to tandem mass spectrometry (or mass spectrometry/mass spectrometry MS/MS) for the identification of product ions corresponding to the peaks in the XIC.
  • the precursor ions can be selected based on their mass/charge ratio to be subjected to subsequent stages of mass analysis.
  • the selected precursor ions can be fragmented (e.g., via collision induced dissociation), and the fragmented ions (product ions) can be analyzed via a subsequent stage of mass spectrometry.
  • peaks present in an XIC There are often multiple peaks present in an XIC.
  • a compound of interest can also be called an analyte of interest. This selection process, however, may be cumbersome, as it requires that the user specify the expected retention time for each compound of interest. Moreover, it can in some cases result in misidentification of the peak corresponding to the compound of interest.
  • sMRM scheduled multiple reaction monitoring
  • Tandem mass spectrometry or MS/MS involves ionization of one or more compounds of interest from a sample, selection of one or more precursor ions of the one or more compounds, fragmentation of the one or more precursor ions into product ions, and mass analysis of the product ions.
  • Tandem mass spectrometry can provide both qualitative and quantitative information.
  • the product ion spectrum can be used to identify a molecule of interest.
  • the intensity of one or more product ions can be used to quantitate the amount of the compound present in a sample.
  • a large number of different types of experimental methods or workflows can be performed using a tandem mass spectrometer.
  • Three broad categories of these workflows are, targeted acquisition, information dependent acquisition (IDA) or data dependent acquisition (DDA), and data independent acquisition (DIA).
  • one or more transitions of a precursor ion to a product ion are predefined for a compound of interest.
  • the one or more transitions are interrogated during each time period or cycle of a plurality of time periods or cycles.
  • the mass spectrometer selects and fragments the precursor ion of each transition and performs a targeted mass analysis for the product ion of the transition.
  • a mass spectrum is produced for each transition.
  • Targeted acquisition methods include, but are not limited to, multiple reaction monitoring (MRM) and selected reaction monitoring (SRM).
  • a user can specify criteria for performing targeted or untargeted mass analysis of product ions while a sample is being introduced into the tandem mass spectrometer. For example, in an IDA method a precursor ion or mass spectrometry (MS) survey scan is performed to generate a precursor ion peak list. The user can select criteria to filter the peak list for a subset of the precursor ions on the peak list. MS/MS is then performed on each precursor ion of the subset of precursor ions. A product ion spectrum is produced for each precursor ion. MS/MS is repeatedly performed on the precursor ions of the subset of precursor ions as the sample is being introduced into the tandem mass spectrometer.
  • MS mass spectrometry
  • sMRM scheduled MRM
  • DIA methods the third broad category of tandem mass spectrometry. These DIA methods have been used to increase the reproducibility and comprehensiveness of data collection from complex samples. DIA methods can also be called non-specific fragmentation methods.
  • the actions of the tandem mass spectrometer are not varied among MS/MS scans based on data acquired in a previous precursor or product ion scan. Instead a precursor ion mass range is selected. A precursor ion mass selection window is then stepped across the precursor ion mass range. All precursor ions in the precursor ion mass selection window are fragmented and all of the product ions of all of the precursor ions in the precursor ion mass selection window are mass analyzed.
  • the precursor ion mass selection window used to scan the mass range can be very narrow so that the likelihood of multiple precursors within the window is small.
  • This type of DIA method is called, for example, MS/MS ALL .
  • a precursor ion mass selection window of about 1 amu is scanned or stepped across an entire mass range.
  • a product ion spectrum is produced for each 1 amu precursor mass window.
  • a product ion spectrum for the entire precursor ion mass range is produced by combining the product ion spectra for each mass selection window.
  • the time it takes to analyze or scan the entire mass range once is referred to as one scan cycle. Scanning a narrow precursor ion mass selection window across a wide precursor ion mass range during each cycle, however, is not practical for some instruments and experiments.
  • a larger precursor ion mass selection window, or selection window with a greater width is stepped across the entire precursor mass range.
  • This type of DIA method is called, for example, SWATH TM acquisition.
  • the precursor ion mass selection window stepped across the precursor mass range in each cycle may have a width of 5-25 amu, or even larger.
  • all the precursor ions in each precursor ion mass selection window are fragmented, and all of the product ions of all of the precursor ions in each mass selection window are mass analyzed.
  • the cycle time can be significantly reduced in comparison to the cycle time of the MS/MS ALL method.
  • Systems and methods are disclosed for identifying actual XIC peaks of compounds of interest from samples so that more accurate expected retention times and more accurate expected retention time windows can be calculated.
  • an actual XIC peak of a compound of interest is identified using a plurality of standard samples.
  • the quantity of the compound of interest is known in each of the plurality of standard samples and this information is used to identify an actual XIC peak of a compound of interest.
  • the system includes a separation device, an ion source device, a mass spectrometer, and a processor.
  • the separation device separates a compound of interest over time from each sample of a plurality of standard samples.
  • the ion source device ionizes and transforms the separated compound of interest into an ion beam.
  • the mass spectrometer mass analyzes the ion beam for at least one ion of the compound of interest, producing a plurality of intensity measurements at a plurality of different acquisition times for the at least one ion of the compound of interest for each sample of the plurality of standard samples.
  • the mass spectrometer can perform MS producing intensity measurements for precursor ions or can perform MS/MS producing intensity measurements for product ions.
  • the processor calculates an XIC for at least one ion of the compound of interest for each sample of the plurality of standard samples from the plurality of intensity measurements, producing an XIC for each sample.
  • the processor finds one or more XIC peaks in each XIC.
  • the processor identifies an XIC peak in each XIC of each sample as an actual peak of the compound of interest, if a ratio of the intensity of the XIC peak in that sample and the intensity of another XIC peak in another sample matches a ratio of the known quantity of the compound of interest in the two samples.
  • an actual XIC peak of a compound of interest is identified using information about other compounds of interest in a plurality of samples. It is known that the XIC peaks of compounds of interest in the same samples have a similar distribution of retention times across those samples, and it is known that the mean of these distributions are correlated in that they differ from a known retention of their compound of interest according to a monotonic function of the acquisition time.
  • This system also includes a separation device, an ion source device, a mass spectrometer, and a processor.
  • the separation device separates a plurality of compounds of interest over time from each sample of a plurality of samples.
  • the ion source device ionizes and transforms the separated compounds of interest into an ion beam.
  • the mass spectrometer mass analyzes the ion beam for at least one ion of each compound of interest of the plurality of compounds of interest, producing a plurality of intensity measurements at a plurality of different acquisition times for the at least one ion of each compound of interest of the plurality of compounds of interest.
  • the processor calculates an XIC for at least one ion of each compound of interest of the plurality of compounds of interest for each sample, producing an XIC for each ion of each compound of interest for each sample.
  • the processor finds one or more XIC peaks in each XIC.
  • the processor calculates for a first XIC peak of a first compound of interest a retention time for each sample and a first distribution function that describes the distribution of the retention times of the first XIC peak across samples.
  • the processor calculates for a second XIC peak of a second compound of interest a retention time for each sample and a second distribution function that describes the distribution of the retention times of the second XIC peak across samples.
  • the processor identifies the first XIC peak as an actual peak, if the first distribution function and the second distribution function have the same width and a first difference between the mean of the first distribution function and the known retention time of the first compound of interest and a second difference between the mean of the second distribution function and the known retention time of the second compound of interest are described by a monotonic function of the plurality of acquisition times.
  • FIG. 1 is a block diagram that illustrates a computer system 100, upon which embodiments of the present teachings may be implemented.
  • Computer system 100 includes a bus 102 or other communication mechanism for communicating information, and a processor 104 coupled with bus 102 for processing information.
  • Computer system 100 also includes a memory 106, which can be a random access memory (RAM) or other dynamic storage device, coupled to bus 102 for storing instructions to be executed by processor 104.
  • Memory 106 also may be used for storing temporary variables or other intermediate information during execution of instructions to be executed by processor 104.
  • Computer system 100 further includes a read only memory (ROM) 108 or other static storage device coupled to bus 102 for storing static information and instructions for processor 104.
  • a storage device 110 such as a magnetic disk or optical disk, is provided and coupled to bus 102 for storing information and instructions.
  • Computer system 100 may be coupled via bus 102 to a display 112, such as a cathode ray tube (CRT) or liquid crystal display (LCD), for displaying information to a computer user.
  • a display 112 such as a cathode ray tube (CRT) or liquid crystal display (LCD)
  • An input device 114 is coupled to bus 102 for communicating information and command selections to processor 104.
  • cursor control 116 is Another type of user input device, such as a mouse, a trackball or cursor direction keys for communicating direction information and command selections to processor 104 and for controlling cursor movement on display 112.
  • This input device typically has two degrees of freedom in two axes, a first axis (i.e., x) and a second axis (i.e., y), that allows the device to specify positions in a plane.
  • a computer system 100 can perform the present teachings. Consistent with certain implementations of the present teachings, results are provided by computer system 100 in response to processor 104 executing one or more sequences of one or more instructions contained in memory 106. Such instructions may be read into memory 106 from another computer-readable medium, such as storage device 110. Execution of the sequences of instructions contained in memory 106 causes processor 104 to perform the process described herein. Alternatively, hard-wired circuitry may be used in place of or in combination with software instructions to implement the present teachings. Thus, implementations of the present teachings are not limited to any specific combination of hardware circuitry and software.
  • Non-volatile media includes, for example, optical or magnetic disks, such as storage device 110.
  • Volatile media includes dynamic memory, such as memory 106.
  • Precursor ion mass selection media includes coaxial cables, copper wire, and fiber optics, including the wires that comprise bus 102.
  • Computer-readable media include, for example, a floppy disk, a flexible disk, hard disk, magnetic tape, or any other magnetic medium, a CD-ROM, digital video disc (DVD), a Blu-ray Disc, any other optical medium, a thumb drive, a memory card, a RAM, PROM, and EPROM, a FLASH-EPROM, any other memory chip or cartridge, or any other tangible medium from which a computer can read.
  • Various forms of computer readable media may be involved in carrying one or more sequences of one or more instructions to processor 104 for execution.
  • the instructions may initially be carried on the magnetic disk of a remote computer.
  • the remote computer can load the instructions into its dynamic memory and send the instructions over a telephone line using a modem.
  • a modem local to computer system 100 can receive the data on the telephone line and use an infra-red transmitter to convert the data to an infra-red signal.
  • An infra-red detector coupled to bus 102 can receive the data carried in the infra-red signal and place the data on bus 102.
  • Bus 102 carries the data to memory 106, from which processor 104 retrieves and executes the instructions.
  • the instructions received by memory 106 may optionally be stored on storage device 110 either before or after execution by processor 104.
  • instructions configured to be executed by a processor to perform a method are stored on a computer-readable medium.
  • the computer-readable medium can be a device that stores digital information.
  • a computer-readable medium includes a compact disc read-only memory (CD-ROM) as is known in the art for storing software.
  • CD-ROM compact disc read-only memory
  • the computer-readable medium is accessed by a processor suitable for executing instructions configured to be executed.
  • MS mass spectrometry
  • LC liquid chromatography
  • peaks present in an XIC there are often multiple peaks present in an XIC.
  • a peak in the XIC that is closest to an expected retention time for a compound of interest, or simply a peak exhibiting the highest intensity within the expected retention time is selected as the peak corresponding to the compound of interest.
  • This selection process may be cumbersome, as it requires that the user specify the expected retention time for each compound of interest. Moreover, it can in some cases result in misidentification of the peak corresponding to the compound of interest.
  • XIC peaks of compounds of interest are distinguished from other XIC peaks using some additional information or pattern found from performing mass spectrometry experiments on a plurality of samples.
  • an XIC peak of a compound of interest is distinguished from other XIC peaks using the additional information provided from performing mass spectrometry experiments on a plurality of quality control (QC) or standard samples for the compound of interest.
  • Standard samples are, for example, a set of samples with known properties, such as a known concentration or known intensity ratio.
  • the varying known quantities of the compound of interest across the standard samples are used to distinguish XIC peaks of the compound of interest from other XIC peaks in the standard samples.
  • an XIC peak of a compound of interest is distinguished from other XIC peaks using the additional information provided from performing mass spectrometry experiments on a plurality of samples that include a plurality of other known compounds of interest.
  • distributions of the retention times of the other known compounds of interest across the plurality of samples are used distinguish XIC peaks of the compound of interest from other XIC peaks in the plurality of samples.
  • results from a sample with an unknown amount of a compound of interest are compared to results from samples with known amounts of a compound of interest.
  • the samples with known amounts of a compound of interest are quality control (QC) or standard samples and include different known quantities of the compound of interest.
  • this information about different known quantities of the compound of interest in different standard samples is used to improve the expected retention time and expected retention time window.
  • These different known quantities of the compound of interest appear as different XIC peak intensities in the results from the different standard samples.
  • Figure 2 is an exemplary plot 200 of the product ion XICs calculated from results of performing LC-MS/MS analysis on 15 different standard samples for the compound Muron, in accordance with various embodiments.
  • Plot 200 illustrates the problem that can be encountered when determining expected retention time and an expected retention time window from actual samples.
  • an initial expected retention time and an initial expected retention time window are parameters entered by a user for a particular experiment. If, for example, an initial expected retention time 210 of 7.5 minutes and an initial expected retention time window 220 between 5 minutes and 10 minutes is entered by a user of the LC-MS/MS system, then a problem is encountered.
  • Initial expected retention time window 220 includes two peaks for each sample: the actual peak and a background peak.
  • the parameter or option can be to select the peak closest to the initial expected retention time or the peak with the largest intensity. In this case, both selections would, in most samples, result in the selection of the wrong peak.
  • the peaks at initial expected retention time 210, 7.5 minutes are the background peaks.
  • the peaks at about 9 minutes are the actual peaks representing the compound Muron.
  • samples 10 through 15 the actual peaks representing the compound Muron do have the larger intensities.
  • the background peaks have a larger intensity.
  • Plot 200 shows how user defined parameters for the initial expected retention time, the initial expected retention time window, and the option for selecting peaks can lead to an incorrect result.
  • Figure 3 is an exemplary plot 300 of the product ion XICs calculated from the results of performing LC-MS/MS analysis on 15 different standard samples for the compound Naled, in accordance with various embodiments.
  • Plot 300 also illustrates the problem that can be encountered when determining expected retention time and an expected retention time window from actual samples.
  • an initial expected retention time window 320 includes actual peaks of Naled at about 11 minutes and background peaks at about 16 minutes. The distance between the actual peaks and background peaks of Naled is larger than it is for Muron in Figure 2 , but the peaks are still close enough to introduce an ambiguity depending on the size of the initial expected retention time window selected.
  • Figure 3 shows another problem that can occur.
  • Some samples may contain no peaks for the compound of interest. For example, samples 1 and 3 of plot 300 have no peaks for Naled. These samples, however, still include peaks of the background.
  • Plot 300 also, however, illustrates the solution to this problem. Note that the actual peaks of Naled at a retention time of about 11 minutes vary in intensity across samples. This intensity variation or intensity pattern is also known, since the samples are standard samples and the quantity of Naled in a standard sample is known. In contrast, the background peaks at a retention time of about 16 minutes do not significantly vary in intensity across samples or do not vary according the known quantity of Naled in the sample.
  • an actual peak of a compound of interest is found in a standard sample by comparing the intensity of that peak to the intensity of at least one other peak in another standard sample. If the intensity ratio of the two peaks is found to match the ratio of the quantities of the compound of interest in the two different standard samples, then both peaks are found to be actual peaks of the compound of interest.
  • sample 10 includes background peak 331 and actual peak 332, and sample 15 includes background peak 341 and actual peak 342.
  • the ratio of the quantities of Naled in sample 15 as compared to sample 10 is known to be 2, or sample 15 has twice as much Naled as sample 10. This is the additional information that is known from the standard samples.
  • the intensity of actual peak 342 of sample 15 is compared to the intensity of actual peak 332 of sample 10, the ratio of the two intensities is also found to be 2, matching the ratio of Naled in the two samples.
  • actual peak 342 of sample 15 and actual peak 332 of sample 10 are found to be peaks of the compound of interest Naled.
  • the expected retention time and the expected retention time are modified based on these peaks.
  • the smallest retention time shift between any two standard samples that provides the intensity pattern of at least one matching peak across samples that matches expected intensity pattern is found. Those peaks are then identified as target peaks and their retention time and retention shift across samples (acquisition times) are used to predict the retention time of the compound of interest in the unknown samples.
  • the expected retention time window is also determined as well as the uncertainty of the solution.
  • the uncertainty measure is used to bring problematic cases to the user's attention and to ask for assistance in order to reduce errors in the final results and also to reduce the time needed for review or correction, for example.
  • Such embodiments utilize known information or features of a subset of samples. In various embodiments, this information can also be combined with measured ion ratios of all peaks, in case of multiple fragments per compound of interest.
  • a matrix of all possible combinations of peaks between two samples is created to find a new or modified expected retention time from the standard samples.
  • the new or modified expected retention time is, therefore, found from all samples or just a subset of samples.
  • the new or modified expected retention time window is calculated from the distribution of retention times across all samples or just a subset of samples. For example, a histogram of retention times or shifts from the expected retention time can be created across all samples or just a subset of samples. The width of the distribution of the histogram or some multiple of the width of the distribution can be used to provide the new or modified expected retention time window.
  • Figure 4 is a plot 400 of a histogram of the actual retention times found for the actual peaks of the compound Naled across the 15 standard samples, in accordance with various embodiments.
  • the actual peaks of Naled ions are found by comparing the peak intensities of peaks in different samples and then comparing the intensity ratio of two peaks to the ratio of the quantities of Naled known to be in the two different samples, as described above.
  • the histogram of plot 400 is plotted by counting the number of samples that have an actual peak with a retention within one of five different retention time ranges.
  • the histogram of plot 400 has a Gaussian distribution or shape 410 and a mean retention time 420 of about 10.73 minutes. If the original expected retention time of Naled is 10.5 minutes, the histogram shows that the retention is shifted about 0.23 minutes in the standard samples. The new or modified expected retention time is then 10.73 minutes.
  • the width of the Gaussian distribution can be selected as the new or modified expected retention time window.
  • a multiple such as three times the 0.5 minute range, can be used.
  • the new or modified expected retention time window can be from 10 to 11.5 minutes.
  • Figure 5 is a different view 500 of the XICs of Figure 3 that are calculated from the results of performing LC-MS/MS analysis on 15 different standard samples for the compound Naled, in accordance with various embodiments.
  • View 500 shows the original expected retention time window 320 that include the actual peaks and the background peaks.
  • View 500 also shows a new or modified expected retention time window 530 from 10 to 11.5 minutes that is determined from the histogram of Figure 4 , for example.
  • Figure 5 shows how the expected retention time window is significantly improved in going from window 320 to window 530.
  • New or modified expected retention time window 530 is now large enough to find a Naled ion peak in a sample where the quantity is unknown, but is small enough to avoid detecting the background peaks.
  • the fact that all compounds of interest in a plurality of samples have a similar Gaussian distribution of actual retention times across the plurality of samples is the additional information that is used to distinguish actual XIC peaks of a compound of interest from background peaks.
  • all compounds of interest have an expected retention time.
  • the distributions of the actual retention times of the compounds of interest are plotted across samples, they are all found to have a similar Gaussian distribution.
  • the shift of the mean of the Gaussian distribution from the expected retention time for each compound of interest is substantially the same.
  • a peak of a compound of interest is found in a plurality of samples by calculating the distribution of the retention time of the peak across the plurality of samples and then comparing the distribution to one or more distributions of the retention times of peaks of other compounds of interest across the same plurality of samples.
  • a peak of a compound of interest is found if the distribution of the retention time of the peak across the plurality of samples matches a distribution of the retention time of a peak of at least one different compound of interest across the plurality of samples in both its shape and in the amount its mean is shifted or differs from the expected retention time.
  • Figure 6 is a series 600 of plots of the distribution of retention times of particular XIC peaks of compounds of interest across a plurality of samples, in accordance with various embodiments.
  • Plot 610 shows a distribution 611 of the retention times of ion peak 1 of compound 1 across a plurality of samples. Peak 1 is one of the peaks found in the XICs calculated for compound 1 for each of the plurality of samples. An XIC is calculated based on the known m/z of an ion of the known compound.
  • the distribution 611 is found, for example, by plotting a histogram of the number of samples that include ion peak 1 at each retention time or retention time range. If the original expected retention time of compound 1 is 3 minutes, then distribution 611 shows a 2 minute shift 612 of the retention time. In other words, the mean of distribution 611 has a shift 612 of 2 minutes from the original expected retention time of 3 minutes.
  • peak 1 is an actual peak, by comparing distribution 611 and shift 612 to the distribution and shift of an ion peak of a different known compound of interest in the plurality of samples.
  • plot 630 shows a distribution 631 of the retention times of ion peak 1 of compound 2 across the plurality of samples. The mean of distribution 631 is also shifted from the original expected retention time of 10 minutes by a shift 632 of 2 minutes.
  • distribution 611 and distribution 631 both have the shape of a Gaussian distribution. Actual peaks are known to have retention times that follow a Gaussian distribution across samples. Distribution 611 and distribution 631 have similar widths, so they are likely similar distributions. Finally, shift 612 of distribution 11 and shift 632 of distribution 631 have the same magnitude and direction. As a result, peak 1 of compound 1 is found to be an actual peak of compound 1 and not a background peak.
  • the retention time distribution of an ion peak of compound of interest is compared to the retention time distributions of ion peaks of more than one different compound of interest.
  • the vast majority of ion peaks found for these compounds can be assumed to be actual peaks of the compound of interest.
  • the largest peak found in an XIC or the peak closest to the expected retention time usually is the actual peak.
  • the peak found can be a background peak or an incorrect interfering peak.
  • plot 640 shows a distribution 641 of the retention times of ion peak 1 of compound 3 across the plurality of samples.
  • the mean of distribution 641 is also shifted from the original expected retention time of 10 minutes by a shift 642 of 2 minutes.
  • the direction of shift 642 is incorrect. This means that peak 1 of compound 3 is actually an incorrect peak.
  • distribution 611 of peak 1 of compound 1 in plot 610 is compared to distribution 641 of peak 1 of compound 3 in plot 630, peak 1 of compound 1 will not be found to be an actual peak, which is an incorrect result.
  • comparing the retention time distribution 611 of peak 1 of compound 1 to other distributions can improve the certainty in determining whether or not peak 1 of compound 1 is an actual peak.
  • Compound 1 can also include a background peak.
  • Plot 620 shows a distribution 621 of the retention times of ion peak 2 of compound 1 across the plurality of samples. Peak 2 of compound 1 is a background peak. A true negative result can be found by comparing plots 620 and 630.
  • Distribution 621 has the shape of a uniform distribution, while distribution 631 has the shape of a Gaussian distribution.
  • Shift 622 of peak 2 of compound 1 has the same magnitude of shift 632 of peak 1 of compound 2. However, shift 622 and 632 differ in direction. This comparison shows that peak 2 of compound 1 is not an actual peak of compound 1.
  • this method of identifying the actual ion peak of compound of interest using ion peaks of other compounds of interest is used in scheduled multiple reaction monitoring (sMRM).
  • sMRM typically involves monitoring of large number of compounds of interest over relatively short time windows so that the number of peak features for XIC alignment are not sufficient for standard alignment approaches that utilize a larger number of matching reference peaks. For example, it is common to have one or two peaks observed in an XIC across multiple samples. Relating single-peak samples with two-peak samples introduces uncertainty. It is reasonable to expect that uncertain cases produce random patterns of retention time shift due to the random error of retention time detection. In contrast, compounds of interest with mostly correctly identified retention times show a similar pattern across the acquisition time (retention time shifts tend to be correlated, not necessarily identical, also the amount of shifting is a monotonic function of the acquisition time).
  • Figure 7 is a schematic diagram 700 showing a system for identifying an XIC peak of a compound of interest, in accordance with various embodiments.
  • the XIC peak is found, in various embodiments, using a plurality of standard samples.
  • the system of Figure 7 includes separation device 710, ion source device 720, mass spectrometer 730, and processor 740.
  • Separation device 710 separates a compound of interest over time from each sample of a plurality of standard samples. Each sample of the plurality of standard samples includes a different known quantity of the compound of interest. Each sample is obtained from sample plate 705, for example. Separation device 710 can perform separation techniques that include, but are not limited to, liquid chromatography, gas chromatography, capillary electrophoresis, or ion mobility.
  • Ion source device 720 ionizes and transforms the separated compound of interest for each sample of the plurality of standard samples into an ion beam.
  • Ion source device 720 can perform ionization techniques that include, but are not limited to, matrix assisted laser desorption/ionization (MALDI) or electrospray ionization (ESI).
  • MALDI matrix assisted laser desorption/ionization
  • ESI electrospray ionization
  • Mass spectrometer 730 receives the ion beam for each sample of the plurality of standard samples. Mass spectrometer 730 mass analyzes the ion beam for at least one ion of the compound of interest. A plurality of intensity measurements at a plurality of different acquisition times are produced for the at least one ion of the compound of interest for each sample of the plurality of standard samples.
  • mass spectrometer 730 is shown as a triple quadrupole device.
  • mass spectrometer 730 can include other types of mass spectrometry devices including, but not limited to, ion traps, orbitraps, time-of-flight (TOF) devices, ion mobility devices, or Fourier transform ion cyclotron resonance (FT-ICR) devices.
  • TOF time-of-flight
  • FT-ICR Fourier transform ion cyclotron resonance
  • mass spectrometer 730 can mass filter or select a precursor ion of the compound of interest by selecting its mass-to-charge ratio at each acquisition time of the plurality of different acquisition times. Mass spectrometer 730 can then measure the intensity of the precursor ion at each acquisition time of the plurality of different acquisition times.
  • the system of Figure 7 can perform LC-MS, for example.
  • mass spectrometer 730 can mass filter or select a precursor ion of the compound of interest by selecting its mass-to-charge ratio at each acquisition time of the plurality of different acquisition times. Mass spectrometer 730 can then fragment the precursor ion into one or more product ions. Finally, mass spectrometer 730 can then measure the intensity of at least one product ion of the compound of interest at each acquisition time of the plurality of different acquisition times.
  • the system of Figure 7 can perform LC-MS/MS, for example.
  • Processor 740 can be, but is not limited to, a computer, a microprocessor, the computer system of Figure 1 , or any device capable of sending and receiving control signals and data from mass spectrometer 730 and processing data.
  • Processor 740 is in communication with ion separation device 720, ion source device 720, and mass spectrometer 730.
  • Processor 740 receives the plurality of intensity measurements for the at least one ion of the compound of interest for each sample of the plurality of standard samples from mass spectrometer 730. Processor 740 calculates an XIC for the at least one ion of the compound of interest for each sample of the plurality of standard samples from the plurality of intensity measurements. An XIC is produced for each sample of the plurality of standard samples.
  • Processor 740 finds one or more XIC peaks in each XIC of each sample of the plurality of standard samples.
  • Processor 740 can, for example, use any conventional peak finding algorithm.
  • Processor 740 identifies an XIC peak of the one or more XIC peaks in each XIC of each sample of the plurality of standard samples as an actual peak of the compound of interest, if a ratio of the intensity of the XIC peak in each sample and the intensity of another XIC peak in another sample matches a ratio of the known quantity of the compound of interest in each sample and the known quantity of the compound of interest in another sample. In other words, processor 740 compares the ratio of XIC peak intensities in different samples to the ratio of known quantities of the compound of interest in the different samples to identify the XIC peak.
  • processor 740 identifies an XIC peak in two or more samples of the plurality of standard samples as an actual peak of the compound of interest and calculates a retention time for each XIC peak identified as an actual peak of the compound of interest. A retention time is produced for each of the two or more samples of the plurality of standard samples.
  • processor 740 calculates a distribution function for the distribution of retention times among the two or more samples of the plurality of standard samples using the retention times calculated for each of the two or more samples of the plurality of standard samples.
  • Figure 4 shows a distribution function 410 calculated for the distribution of retention times among two or more samples of the plurality of standard samples using the retention times calculated for each of the two or more samples of the plurality of standard samples.
  • processor 740 of Figure 7 calculates an expected retention time and an expected retention time window from the distribution function. For example, processor 740 calculates an expected retention time as the expected value of the distribution function.
  • the expected value of the distribution function can be estimated using various statistical methods for distribution parameter estimation.
  • the expected value can include the mean and/or the median. Both can be used if the distribution is binomial (double hump) or just one can be used if the distribution has a single hump, like in a Gaussian distribution.
  • the expected value can be found from a hypothesis test that distribution has an expected value equal to average value calculated from observation. If the test fails, then the distribution is not Gaussian-like and other modeling techniques can be applied. In practice, for single batch, it is most likely Gaussian and if it is not Gaussian, then it is sign of chromatographic system property change. A sign of chromatographic system property change can be used for signaling issue, or to automatically adjust scheduling windows and other dependent parameters of the acquisition in order to make sure all samples will have useful measurements (will not cut-off the peak of interest due to inappropriate scheduling window).
  • Figure 4 shows an exemplary mean 420 calculated from the distribution function 410 that can be used as a new or modified expected retention time.
  • processor 740 of Figure 7 calculates the width of the expected retention time window as the width of the distribution function or as a multiple of the width of the distribution function.
  • Figure 4 shows an exemplary distribution function 410 with a width of about 0.5 min. This width or a multiple of this width can be used as a new or modified expected retention time window.
  • the same system can be used to identify an extracted ion chromatogram (XIC) peak of a compound of interest using a plurality of samples that include a plurality of other compounds of interest.
  • XIC extracted ion chromatogram
  • Separation device 710 separates a plurality of compounds of interest over time from each sample of a plurality of samples.
  • Ion source device 720 ionizes and transforms the separating compounds of interest of the plurality of compounds of interest for each sample of the plurality of samples into an ion beam.
  • Mass spectrometer 730 receives the ion beam for each sample of the plurality of samples. Mass spectrometer 730 mass analyzes the ion beam for at least one ion of each compound of interest of the plurality of compounds of interest. A plurality of intensity measurements at a plurality of different acquisition times are produced for the at least one ion of each compound of interest of the plurality of compounds of interest for each sample of the plurality of samples.
  • Processor 740 receives the plurality of intensity measurements for the at least one ion of each compound of interest of the plurality of compounds of interest for each sample of the plurality of sample from mass spectrometer 730. Processor 740 also receives a known retention time for each compound of interest of the plurality of compounds of interest from a user, for example.
  • Processor 740 calculates an XIC for the at least one ion of each compound of interest of the plurality of compounds of interest for each sample of the plurality of samples from the plurality of intensity measurements. An XIC is produced for the at least one ion of each compound of interest of the plurality of compounds of interest for each sample of the plurality of samples. Processor 740 finds one or more XIC peaks in each XIC for the at least one ion of each compound of interest of the plurality of compounds of interest for each sample of the plurality of samples.
  • Processor 740 calculates for a first XIC peak of an ion of a first compound of interest of the plurality of compounds of interest a retention time for each sample of the plurality of samples and a first distribution function that describes the distribution of the retention times of the first XIC peak across the plurality of samples.
  • Plot 610 of Figure 6 shows an exemplary first distribution function 611 for the retention times of a first XIC peak across a plurality of samples.
  • Processor 740 of Figure 7 calculates for a second XIC peak of an ion of a first compound of interest of the plurality of compounds of interest a retention time for each sample of the plurality of samples and a second distribution function that describes the distribution of the retention times of the second XIC peak across the plurality of samples.
  • Plot 630 of Figure 6 shows an exemplary second distribution function 631 for the retention times of a second XIC peak across a plurality of samples.
  • processor 740 of Figure 7 identifies the first XIC peak as an actual peak of the first compound of interest, if the first distribution function and the second distribution function have the same width and a first difference between the mean of the first distribution function and the known retention time of the first compound of interest and a second difference between the mean of the second distribution function and the known retention time of the second compound of interest are described by a monotonic function of the plurality of acquisition times.
  • Exemplary first distribution function 611 of plot 610 of Figure 6 and exemplary second distribution function 631 of plot 630 of Figure 6 have the same width.
  • Exemplary first distribution function 611 has a first difference 612 between the mean of first distribution function 611 and the known retention time of the first compound of interest.
  • Exemplary second distribution function 631 has a second difference 632 between the mean of second distribution function 631 and the known retention time of the second compound of interest.
  • First difference 612 and second difference 632 are described by a monotonic function of the plurality of acquisition times. In fact, first difference 612 and second difference 632 are both 2 minutes.
  • a monotonic function is a function that is non-increasing or non-decreasing.
  • the first XIC peak of a first compound of interest that has first distribution function 611 can be identified as an actual peak using the second distribution function 631 of a second XIC peak of a second compound of interest.
  • an XIC peak can be identified by comparing its distribution to more than one distribution of other XIC peaks of other compounds of interest.
  • processor 740 of Figure 7 can further calculate for one or more other XIC peaks of ions of one or more other compounds of interest of the plurality of compounds of interest a retention time for each sample of the plurality of samples and a distribution function that describes the distribution of the retention times of the one or more other XIC peaks across the plurality of samples.
  • Processor 740 can then identify the first XIC peak as an actual peak of the first compound of interest, if the first distribution function and the distribution functions of the one or more other XIC peaks have the same width and a first difference between the mean of the first distribution function and the known retention time of the first compound of interest and the differences between the means of the distribution functions of the one or more other XIC peaks and the known retention times of the one or more other compounds of interest are described by a monotonic function of the plurality of acquisition times.
  • processor 740 calculates an expected retention time and an expected retention time window from the first distribution function. For example, processor 740 calculates an expected retention time as the mean of the first distribution function. Processor 740 can calculate the width of the expected retention time window as the width of the first distribution function or as a multiple of the width of the first distribution function.
  • the specification may have presented a method and/or process as a particular sequence of steps.
  • the method or process should not be limited to the particular sequence of steps described.
  • other sequences of steps may be possible. Therefore, the particular order of the steps set forth in the specification should not be construed as limitations on the claims.
  • the claims directed to the method and/or process should not be limited to the performance of their steps in the order written, and one skilled in the art can readily appreciate that the sequences may be varied and still remain within the scope of the various embodiments.

Landscapes

  • Chemical & Material Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Physics & Mathematics (AREA)
  • Biochemistry (AREA)
  • General Health & Medical Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Immunology (AREA)
  • Pathology (AREA)
  • Engineering & Computer Science (AREA)
  • Library & Information Science (AREA)
  • Chemical Kinetics & Catalysis (AREA)
  • Other Investigation Or Analysis Of Materials By Electrical Means (AREA)

Claims (13)

  1. Système pour identifier un pic de chromatogramme ionique extrait (XIC) d'un composé d'intérêt en utilisant une pluralité d'échantillons standards, comprenant :
    un dispositif de séparation (710) pour séparer un composé d'intérêt dans le temps de chaque échantillon d'une pluralité d'échantillons standards, dans lequel chaque échantillon de la pluralité d'échantillons standards comprend une quantité connue différente du composé d'intérêt ;
    un dispositif de source ionique (720) pour ioniser et transformer le composé d'intérêt séparé pour chaque échantillon de la pluralité d'échantillons standards en un faisceau ionique ;
    un spectromètre de masse (730) pour recevoir le faisceau ionique pour chaque échantillon de la pluralité d'échantillons standards et pour analyser en masse le faisceau ionique pour au moins un ion du composé d'intérêt, produisant une pluralité de mesures d'intensité à une pluralité de temps d'acquisition différents pour l'au moins un ion du composé d'intérêt pour chaque échantillon de la pluralité d'échantillons standards ;
    un processeur (740) en communication avec le spectromètre de masse pour
    recevoir la pluralité de mesures d'intensité de l'au moins un ion du composé d'intérêt pour chaque échantillon de la pluralité d'échantillons standards depuis le spectromètre de masse,
    calculer un XIC pour l'au moins un ion du composé d'intérêt pour chaque échantillon de la pluralité d'échantillons standards à partir de la pluralité de mesures d'intensité, produisant un XIC pour chaque échantillon de la pluralité d'échantillons standards,
    trouver un ou plusieurs pics XIC dans chaque XIC de chaque échantillon de la pluralité d'échantillons standards, et
    identifier un pic XIC de l'un ou plusieurs pics XIC dans chaque XIC de chaque échantillon de la pluralité d'échantillons standards en tant que pic réel du composé d'intérêt, si un rapport de l'intensité du pic XIC dans chaque échantillon et l'intensité d'un autre pic XIC dans un autre échantillon correspondent à un rapport entre la quantité connue du composé d'intérêt dans chaque échantillon et la quantité connue du composé d'intérêt dans l'autre échantillon.
  2. Système selon la revendication 1, dans lequel le processeur (740) est configuré pour identifier un pic XIC dans deux ou plusieurs échantillons de la pluralité d'échantillons standards comme un pic réel du composé d'intérêt et calculer un temps de rétention pour chaque pic XIC identifié comme un pic réel du composé d'intérêt, produisant un temps de rétention pour chacun des deux ou plusieurs échantillons de la pluralité d'échantillons standards.
  3. Système selon la revendication 2, dans lequel le processeur (740) est configuré pour calculer une fonction de distribution pour la distribution des temps de rétention parmi les deux ou plusieurs échantillons de la pluralité d'échantillons standards en utilisant les temps de rétention calculés pour chacun des deux ou plusieurs échantillons de la pluralité d'échantillons standards.
  4. Système selon la revendication 3, dans lequel le processeur (740) est configuré pour calculer un temps de rétention attendu et une fenêtre de temps de rétention attendu à partir de la fonction de distribution.
  5. Système selon la revendication 4, dans lequel le processeur (740) est configuré pour calculer un temps de rétention attendu comme valeur attendue de la fonction de distribution.
  6. Système selon la revendication 4, dans lequel le processeur (740) est configuré pour calculer une largeur de la fenêtre de temps de rétention attendu comme un multiple de la largeur de la fonction de distribution.
  7. Système selon la revendication 6, dans lequel le multiple de la largeur de la fonction de distribution est la largeur de la fonction de distribution.
  8. Système pour identifier un pic de chromatogramme ionique extrait (XIC) d'un composé d'intérêt en utilisant une pluralité d'échantillons qui incluent une pluralité d'autres composés d'intérêt, comprenant :
    un dispositif de séparation (710) pour séparer une pluralité de composés connus d'intérêt dans le temps de chaque échantillon d'une pluralité d'échantillons ;
    un dispositif de source ionique (720) pour ioniser et transformer les composés d'intérêt séparés de la pluralité de composés d'intérêt pour chaque échantillon de la pluralité d'échantillons en un faisceau ionique ;
    un spectromètre de masse (730) pour recevoir le faisceau ionique pour chaque échantillon de la pluralité d'échantillons et pour analyser en masse le faisceau ionique pour au moins un ion de chaque composé d'intérêt de la pluralité de composés d'intérêt, produisant une pluralité de mesures d'intensité à une pluralité de temps d'acquisition différents pour l'au moins un ion de chaque composé d'intérêt de la pluralité de composés d'intérêt pour chaque échantillon de la pluralité d'échantillons ;
    un processeur (740) en communication avec le spectromètre de masse pour
    recevoir la pluralité de mesures d'intensité pour au moins un ion de chaque composé d'intérêt de la pluralité de composés d'intérêt pour chaque échantillon de la pluralité d'échantillons depuis le spectromètre de masse et un temps de rétention connu pour chaque composé d'intérêt de la pluralité de composés d'intérêt depuis un utilisateur,
    calculer un XIC pour l'au moins un ion de chaque composé d'intérêt de la pluralité de composés d'intérêt pour chaque échantillon de la pluralité d'échantillons à partir de la pluralité de mesures d'intensité, produisant un XIC pour l'au moins un ion de chaque composé d'intérêt de la pluralité de composés d'intérêt pour chaque échantillon de la pluralité d'échantillons,
    trouver un ou plusieurs pics XIC dans chaque XIC pour l'au moins un ion de chaque composé d'intérêt de la pluralité de composés d'intérêt pour chaque échantillon de la pluralité d'échantillons,
    calculer pour un premier pic XIC d'un ion d'un premier composé d'intérêt de la pluralité de composés d'intérêt un temps de rétention pour chaque échantillon de la pluralité d'échantillons et une première fonction de distribution qui décrit la distribution des temps de rétention du premier pic XIC à travers la pluralité d'échantillons,
    calculer pour un deuxième pic XIC d'un ion d'un deuxième composé d'intérêt de la pluralité de composés d'intérêt un temps de rétention pour chaque échantillon de la pluralité d'échantillons et une deuxième fonction de distribution qui décrit la distribution des temps de rétention du deuxième pic XIC à travers la pluralité de échantillons, et
    identifier le premier pic XIC comme un pic réel du premier composé d'intérêt, si la première fonction de distribution et la deuxième fonction de distribution ont la même largeur et une première différence entre la moyenne de la première fonction de distribution et le temps de rétention connu du premier composé d'intérêt et une deuxième différence entre la moyenne de la deuxième fonction de distribution et le temps de rétention connu du deuxième composé d'intérêt sont décrites par une fonction monotone de la pluralité de temps d'acquisition.
  9. Système selon la revendication 8, dans lequel le processeur (740) est configuré pour
    calculer pour un ou plusieurs autres pics XIC d'ions d'un ou plusieurs autres composés d'intérêt de la pluralité de composés d'intérêt un temps de rétention pour chaque échantillon de la pluralité d'échantillons et une fonction de distribution qui décrit la distribution des temps de rétention d'un ou plusieurs autres pics XIC sur la pluralité d'échantillons et
    identifier le premier pic XIC comme un pic réel du premier composé d'intérêt, si la première fonction de distribution et les fonctions de distribution de l'un ou plusieurs autres pics XIC ont la même largeur et une première différence entre la moyenne de la première fonction de distribution et le temps de rétention connu du premier composé d'intérêt et les différences entre les moyennes des fonctions de distribution de l'un ou plusieurs autres pics XIC et les temps de rétention connus de l'un ou plusieurs autres composés d'intérêt sont décrites par une fonction monotone de la pluralité de temps d'acquisition.
  10. Système selon la revendication 8, dans lequel le processeur (740) est configuré pour calculer un temps de rétention attendu et une fenêtre de temps de rétention attendu à partir de la première fonction de distribution.
  11. Système selon la revendication 10, dans lequel le processeur (740) est configuré pour calculer un temps de rétention attendu comme la moyenne de la première fonction de distribution.
  12. Système selon la revendication 10, dans lequel le processeur (740) est configuré pour calculer la largeur de la fenêtre de temps de rétention attendu comme la largeur de la première fonction de distribution.
  13. Système selon la revendication 10, dans lequel le processeur (740) est configuré pour calculer la largeur de la fenêtre de temps de rétention attendu comme un multiple de la largeur de la première fonction de distribution.
EP17885088.9A 2016-12-20 2017-12-01 Détermination automatique de temps de rétention prévu et de fenêtre de temps de rétention optimale prévue Active EP3559658B1 (fr)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
US201662436633P 2016-12-20 2016-12-20
PCT/IB2017/057560 WO2018116039A1 (fr) 2016-12-20 2017-12-01 Détermination automatique de temps de rétention prévu et de fenêtre de temps de rétention optimale prévue

Publications (3)

Publication Number Publication Date
EP3559658A1 EP3559658A1 (fr) 2019-10-30
EP3559658A4 EP3559658A4 (fr) 2020-08-19
EP3559658B1 true EP3559658B1 (fr) 2022-11-16

Family

ID=62625988

Family Applications (1)

Application Number Title Priority Date Filing Date
EP17885088.9A Active EP3559658B1 (fr) 2016-12-20 2017-12-01 Détermination automatique de temps de rétention prévu et de fenêtre de temps de rétention optimale prévue

Country Status (3)

Country Link
US (1) US11474087B2 (fr)
EP (1) EP3559658B1 (fr)
WO (1) WO2018116039A1 (fr)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2021122784A1 (fr) 2019-12-17 2021-06-24 Roche Diagnostics Gmbh Procédé et dispositif de surveillance de transitions multiples
US20240183829A1 (en) * 2022-12-05 2024-06-06 Halliburton Energy Services, Inc. Predictive chromatograph peak detection

Family Cites Families (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP3801355B2 (ja) * 1998-07-16 2006-07-26 フロンティア・ラボ株式会社 未知試料検索方法、未知試料検索装置及び未知試料検索用記録媒体
EP1697024B1 (fr) * 2003-12-18 2017-08-30 DH Technologies Development Pte. Ltd. Procedes et appareil permettant d'ameliorer la detection d'echantillons fondee sur des ions au moyen d'une pre-separation et d'une amplification selectives
WO2005079263A2 (fr) 2004-02-13 2005-09-01 Waters Investments Limited Appareil et procede d'identification de pics dans des donnes de spectrometrie de masse/chromatographie liquide et de formation de spectres et de chromatogrammes
US20130080073A1 (en) * 2010-06-11 2013-03-28 Waters Technologies Corporation Techniques for mass spectrometry peak list computation using parallel processing
US9640373B2 (en) * 2011-06-29 2017-05-02 Shimadzu Corporation Method and system for processing analysis data
EP2798664B1 (fr) * 2011-12-29 2020-10-28 DH Technologies Development Pte. Ltd. Utilisation de données de spectrométrie de masse en fenêtre pour confirmation ou détermination de temps de maintien
EP2834835B1 (fr) * 2012-04-02 2018-11-28 Thermo Fisher Scientific (Bremen) GmbH Procédé et appareil pour la quantification améliorée par spectrométrie de masse
GB201306868D0 (en) * 2013-04-15 2013-05-29 Micromass Ltd A method of screening samples
US10121643B2 (en) * 2014-07-24 2018-11-06 Shimadzu Corporation Chromatography/mass spectrometry data processing device
US9804136B2 (en) * 2014-09-18 2017-10-31 Dionex Corporation Automated method of calibrating a chromatography system and analysis of a sample
WO2016125060A1 (fr) * 2015-02-05 2016-08-11 Dh Technologies Development Pte. Ltd. Détection de similarité basée sur la spectrométrie de masse par le biais d'une soustraction de courbes
CN108780071B (zh) * 2016-03-07 2021-06-08 株式会社岛津制作所 气相色谱质谱分析装置

Also Published As

Publication number Publication date
EP3559658A1 (fr) 2019-10-30
EP3559658A4 (fr) 2020-08-19
US11474087B2 (en) 2022-10-18
US20200096489A1 (en) 2020-03-26
WO2018116039A1 (fr) 2018-06-28

Similar Documents

Publication Publication Date Title
US12033839B2 (en) Data independent acquisition of product ion spectra and reference spectra library matching
US9343276B2 (en) Use of windowed mass spectrometry data for retention time determination or confirmation
US9583323B2 (en) Use of variable XIC widths of TOF-MSMS data for the determination of background interference in SRM assays
CN109526237B (zh) 使用swath方法扩展动态范围
EP4393003A1 (fr) Procédé d'amélioration d'informations dans la spectrométrie de masse de dda
EP3559658B1 (fr) Détermination automatique de temps de rétention prévu et de fenêtre de temps de rétention optimale prévue
WO2016125061A1 (fr) Balayage rapide de grandes fenêtres rf quadripolaires effectué pendant le basculement simultané de l'énergie de fragmentation
WO2020044161A1 (fr) Identification et notation de composés apparentés dans des échantillons complexes
US11953478B2 (en) Agnostic compound elution determination
EP3688788B1 (fr) Évaluation de la pureté de pic mrm avec des ms/ms sélectifs d'isotopes
US20240355604A1 (en) Data independent acquisition of product ion spectra and reference spectra library matching
WO2023199139A1 (fr) Optimisation de paramètres de traitement pour analyse ms/ms top/middle down
WO2024075058A1 (fr) Réduction de la complexité de données pour un alignement de temps de retention (rt) ultérieur
CN117242543A (zh) 线性定量动态范围延展方法
WO2023199138A1 (fr) Notation de spectres msms de protéine entière sur la base d'une note de pertinence de liaison
WO2022091047A1 (fr) Identification de composés par spectrométrie de masse

Legal Events

Date Code Title Description
STAA Information on the status of an ep patent application or granted ep patent

Free format text: STATUS: THE INTERNATIONAL PUBLICATION HAS BEEN MADE

PUAI Public reference made under article 153(3) epc to a published international application that has entered the european phase

Free format text: ORIGINAL CODE: 0009012

STAA Information on the status of an ep patent application or granted ep patent

Free format text: STATUS: REQUEST FOR EXAMINATION WAS MADE

17P Request for examination filed

Effective date: 20190717

AK Designated contracting states

Kind code of ref document: A1

Designated state(s): AL AT BE BG CH CY CZ DE DK EE ES FI FR GB GR HR HU IE IS IT LI LT LU LV MC MK MT NL NO PL PT RO RS SE SI SK SM TR

AX Request for extension of the european patent

Extension state: BA ME

DAV Request for validation of the european patent (deleted)
DAX Request for extension of the european patent (deleted)
A4 Supplementary search report drawn up and despatched

Effective date: 20200720

RIC1 Information provided on ipc code assigned before grant

Ipc: G01N 30/72 20060101ALI20200714BHEP

Ipc: H01J 49/00 20060101ALI20200714BHEP

Ipc: G01N 33/68 20060101ALI20200714BHEP

Ipc: G01N 30/86 20060101AFI20200714BHEP

GRAP Despatch of communication of intention to grant a patent

Free format text: ORIGINAL CODE: EPIDOSNIGR1

STAA Information on the status of an ep patent application or granted ep patent

Free format text: STATUS: GRANT OF PATENT IS INTENDED

RIC1 Information provided on ipc code assigned before grant

Ipc: G01N 30/04 20060101ALN20220509BHEP

Ipc: G01N 30/88 20060101ALN20220509BHEP

Ipc: G01N 30/72 20060101ALN20220509BHEP

Ipc: H01J 49/00 20060101ALI20220509BHEP

Ipc: G01N 33/68 20060101ALI20220509BHEP

Ipc: G01N 30/86 20060101AFI20220509BHEP

INTG Intention to grant announced

Effective date: 20220530

GRAS Grant fee paid

Free format text: ORIGINAL CODE: EPIDOSNIGR3

GRAA (expected) grant

Free format text: ORIGINAL CODE: 0009210

STAA Information on the status of an ep patent application or granted ep patent

Free format text: STATUS: THE PATENT HAS BEEN GRANTED

AK Designated contracting states

Kind code of ref document: B1

Designated state(s): AL AT BE BG CH CY CZ DE DK EE ES FI FR GB GR HR HU IE IS IT LI LT LU LV MC MK MT NL NO PL PT RO RS SE SI SK SM TR

REG Reference to a national code

Ref country code: GB

Ref legal event code: FG4D

REG Reference to a national code

Ref country code: CH

Ref legal event code: EP

REG Reference to a national code

Ref country code: DE

Ref legal event code: R096

Ref document number: 602017063839

Country of ref document: DE

REG Reference to a national code

Ref country code: IE

Ref legal event code: FG4D

REG Reference to a national code

Ref country code: AT

Ref legal event code: REF

Ref document number: 1532058

Country of ref document: AT

Kind code of ref document: T

Effective date: 20221215

REG Reference to a national code

Ref country code: LT

Ref legal event code: MG9D

REG Reference to a national code

Ref country code: NL

Ref legal event code: MP

Effective date: 20221116

REG Reference to a national code

Ref country code: AT

Ref legal event code: MK05

Ref document number: 1532058

Country of ref document: AT

Kind code of ref document: T

Effective date: 20221116

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: SE

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20221116

Ref country code: PT

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20230316

Ref country code: NO

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20230216

Ref country code: LT

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20221116

Ref country code: FI

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20221116

Ref country code: ES

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20221116

Ref country code: AT

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20221116

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: RS

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20221116

Ref country code: PL

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20221116

Ref country code: LV

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20221116

Ref country code: IS

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20230316

Ref country code: HR

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20221116

Ref country code: GR

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20230217

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: NL

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20221116

P01 Opt-out of the competence of the unified patent court (upc) registered

Effective date: 20230601

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: SM

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20221116

Ref country code: RO

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20221116

Ref country code: EE

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20221116

Ref country code: DK

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20221116

Ref country code: CZ

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20221116

REG Reference to a national code

Ref country code: CH

Ref legal event code: PL

REG Reference to a national code

Ref country code: DE

Ref legal event code: R097

Ref document number: 602017063839

Country of ref document: DE

REG Reference to a national code

Ref country code: BE

Ref legal event code: MM

Effective date: 20221231

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: SK

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20221116

Ref country code: LU

Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES

Effective date: 20221201

Ref country code: AL

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20221116

PLBE No opposition filed within time limit

Free format text: ORIGINAL CODE: 0009261

STAA Information on the status of an ep patent application or granted ep patent

Free format text: STATUS: NO OPPOSITION FILED WITHIN TIME LIMIT

26N No opposition filed

Effective date: 20230817

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: LI

Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES

Effective date: 20221231

Ref country code: IE

Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES

Effective date: 20221201

Ref country code: CH

Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES

Effective date: 20221231

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: SI

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20221116

Ref country code: BE

Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES

Effective date: 20221231

PGFP Annual fee paid to national office [announced via postgrant information from national office to epo]

Ref country code: FR

Payment date: 20230929

Year of fee payment: 7

PGFP Annual fee paid to national office [announced via postgrant information from national office to epo]

Ref country code: GB

Payment date: 20231012

Year of fee payment: 7

PGFP Annual fee paid to national office [announced via postgrant information from national office to epo]

Ref country code: DE

Payment date: 20231003

Year of fee payment: 7

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: HU

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT; INVALID AB INITIO

Effective date: 20171201

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: CY

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20221116

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: MK

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20221116

Ref country code: IT

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20221116

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: MC

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20221116

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: MC

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20221116

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: BG

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20221116

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: MT

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20221116