WO2011128703A1 - Ion mobility method and apparatus for identifying and/or characterising a sample compound - Google Patents
Ion mobility method and apparatus for identifying and/or characterising a sample compound Download PDFInfo
- Publication number
- WO2011128703A1 WO2011128703A1 PCT/GB2011/050756 GB2011050756W WO2011128703A1 WO 2011128703 A1 WO2011128703 A1 WO 2011128703A1 GB 2011050756 W GB2011050756 W GB 2011050756W WO 2011128703 A1 WO2011128703 A1 WO 2011128703A1
- Authority
- WO
- WIPO (PCT)
- Prior art keywords
- collision cross
- compound
- section
- ion mobility
- theoretical
- Prior art date
Links
Classifications
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16C—COMPUTATIONAL CHEMISTRY; CHEMOINFORMATICS; COMPUTATIONAL MATERIALS SCIENCE
- G16C10/00—Computational theoretical chemistry, i.e. ICT specially adapted for theoretical aspects of quantum chemistry, molecular mechanics, molecular dynamics or the like
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16C—COMPUTATIONAL CHEMISTRY; CHEMOINFORMATICS; COMPUTATIONAL MATERIALS SCIENCE
- G16C20/00—Chemoinformatics, i.e. ICT specially adapted for the handling of physicochemical or structural data of chemical particles, elements, compounds or mixtures
- G16C20/20—Identification of molecular entities, parts thereof or of chemical compositions
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N27/00—Investigating or analysing materials by the use of electric, electrochemical, or magnetic means
- G01N27/62—Investigating or analysing materials by the use of electric, electrochemical, or magnetic means by investigating the ionisation of gases, e.g. aerosols; by investigating electric discharges, e.g. emission of cathode
- G01N27/622—Ion mobility spectrometry
- G01N27/623—Ion mobility spectrometry combined with mass spectrometry
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N27/00—Investigating or analysing materials by the use of electric, electrochemical, or magnetic means
- G01N27/62—Investigating or analysing materials by the use of electric, electrochemical, or magnetic means by investigating the ionisation of gases, e.g. aerosols; by investigating electric discharges, e.g. emission of cathode
- G01N27/64—Investigating or analysing materials by the use of electric, electrochemical, or magnetic means by investigating the ionisation of gases, e.g. aerosols; by investigating electric discharges, e.g. emission of cathode using wave or particle radiation to ionise a gas, e.g. in an ionisation chamber
-
- 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
Definitions
- This invention relates generally to a method and apparatus for identifying and/or characterising a sample compound. More specifically, although not exclusively, this invention relates to such a method for differentiating isomeric or isobaric compounds.
- MS mass spectrometry
- MS/MS tandem MS
- substituted molecules such as hydroxylated metabolites which produce identical MS/MS spectra. This is problematic if the exact position of the hydroxylation is required.
- Ion mobility (IM) has the ability to separate isomeric or isobaric species, such as hydroxylated metabolites, rapidly (msec) based on differences in their collision cross-section ( ⁇ ; physical size, and shape) in the gas-phase, thus providing specific information on ionic configuration and therefore, the position of the hydroxyl moiety.
- ⁇ collision cross-section
- the present invention provides methods and apparatus that are particularly suited for identification of such compounds where specific uses of these advances are made. More specifically, the methods and apparatus of the present invention enable more accurate characterisation and/or identification of such compounds.
- One aspect of the invention provides a method of identifying a sample compound, the method comprising the steps of: (i) proving models of an ensemble of possible or expected structures for each of two or more known compounds; (ii) calculating a theoretical collision cross-section for each modelled structure; (iii) averaging the calculated theoretical collision cross-sections for each known compound to provide a theoretical collision cross-section value for each known compound; (iv) measuring a collision cross-section value for the sample compound using an ion mobility technique or cell; and (v) comparing the measured collision cross-section value with the calculated theoretical collision cross-section values to identify which of the two or more known compounds the sample compound most closely resembles.
- a second aspect of the invention provides a method of characterising the structure of a sample compound, the method comprising the steps of: (i) providing models of an ensemble of possible or expected structures for each of two or more known compounds; (ii) calculating a theoretical collision cross-section for each modelled structure; (iii) averaging the calculated theoretical collision cross-sections for each known compound to provide a theoretical collision cross-section value for each known compound; (iv) measuring a collision cross-section value for the sample compound using an ion mobility technique or cell; and (v) comparing the measured collision cross-section value with the theoretical collision cross-section values to characterise the structure of a sample compound.
- the two or more known compounds may comprise two or more known isomeric or isobaric compounds, for example metabolites, e.g. hydroxylated metabolites.
- the ion mobility technique or cell preferably comprises a travelling wave ion mobility technique or cell.
- the method may further comprise the step of analysing a sample compound, for example a further sample compound, that may be from the same sample, e.g. using a mass spectrometer, for example to produce a mass spectrum.
- This analysing step may be carried out before the comparing step, e.g. in order to identify and/or further characterise the sample compound.
- the step of providing models may further comprise checking simplified molecular input line entry specification strings and/or protonates / deprotonates ionisation centres, for example according to physiological pH. Additionally or alternatively, the step of providing models may comprise generating three dimensional structures and may further comprise energy minimising the three dimensional structures.
- Another aspect of the invention provides a computer program element, which may comprise computer readable program code means, e.g. for causing a processor to execute a procedure, for example to implement one or more of the above steps.
- the computer program element may be embodied on a computer readable medium.
- a further aspect of the invention provides a computer readable medium, preferably having a program stored thereon, for example where the program is to make a computer execute a procedure, e.g. to implement one or more of the above steps.
- a yet further aspect of the invention provides an analytical apparatus suitable for carrying out a method as described above and/or comprising a program element as described above and/or a computer readable medium as described above, which apparatus may be specifically adapted for carrying out the method.
- the analytical instrument may comprise a mass spectrometer.
- a yet further aspect of the invention provides a retrofit kit for adapting an analytical apparatus or mass spectrometer to provide an apparatus or mass spectrometer as described above.
- the kit may comprise a program element as described above and/or a computer readable medium as described above.
- Figure 1 illustrates a parent drug ondansetron and three hydroxylated metabolites GR60661 , GR90315 and GR63418;
- Figure 2 illustrates the ensemble of hydroxylated metabolite GR60661 structures created by molecular dynamics and structure refinement, wherein the structures differ by degree rotation around the two rotatable bonds;
- Figure 4 illustrates the UPLC/MS separation of ondansetron and hydroxylated metabolites GR60661 , GR63418 and GR90315
- Figure 5 illustrates the collision induced dissociation MS/MS spectrum for the hydroxylated metabolite GR60661 m/z 310.1 ;
- Figure 6 is a table showing the comparison of the collision cross-section ( ⁇ , in A 2 ) values of ondansetron and hydroxylated metabolites GR60661 , GR63418 and GR90315 obtained using travelling wave ion mobility measurements, WatersTM Projection Approximation (PA) and the MOBCALTM Trajectory Method (TM) calculations;
- PA WatersTM Projection Approximation
- TM MOBCALTM Trajectory Method
- Figure 7A is an overlay of ion mobility arrival time distributions for the parent compounds ondansteron and metabolites GR90315, GR63418 and GR60661 , in order of ion mobility elution (travelling wave 1 100m/sec, 3.3 mbar nitrogen);
- Figure 7B is an overlay of ion mobility arrival time distributions for metabolites GR63418 and GR60661 (travelling wave 1 100m/sec, 3.3 mbar nitrogen);
- Figure 8 shows a DriftScope 2.1 plot of drift time (msec) verses retention time (mins) of the LC-IMS-MS analysis of ondansetron, GR60661 , GR63418 and GR90315 and assigned travelling wave collision cross-section values;
- Figure 9 is a table showing the comparison of the collision cross-section ( ⁇ , in A 2 ) values of methylated imidazole neutral loss fragment of ondansetron and hydroxylated metabolites GR60661 , GR63418 and GR90315 obtained using travelling wave ion mobility measurements and the MOBCALTM Trajectory Method calculations.
- the ⁇ value of the parent compound ondansetron and the three hydroxylated metabolites GR60661 , GR63418 and GR90315 were analysed by UPLC-IMS-MS. All three hydroxylated metabolites showed different chromatographic elution time. However, the MS and MS/MS spectra are all identical, thus hampering identification.
- the theoretically derived ⁇ values for hydroxylated metabolites GR60661 , GR63418 and GR90315 are 1 12.5 A2, 1 1 1 1 .8 A2 and 1 10.9 A2 respectively. All reported values based on the Projection Approximation and Trajectory Method.
- Travelling wave ion mobility derived collision cross-section values are in very close agreements with the theoretically derived values, thus allowing accurate identification of the hydroxylated metabolites, even those with a collision cross-section difference of less than 1 A2.
- No additional structural characterization techniques, such as NMR or X-ray crystallography were required for metabolite structural determination. This is a very time efficient and effective means of identifying drug metabolites directly from biological matrices. Additionally it can be utilised in a situation where chemically synthesised standards may not be available.
- the invention provides a new methodology which combines molecular modelling, theoretical ⁇ calculations and travelling wave ion mobility derived ⁇ measurements, to successfully indentify a number of drug metabolites, which cannot be differentiated by traditional MS/MS techniques.
- In-silico molecular modelling is used to create an ensemble (approximately 10,000) of possible metabolite structures (in this case study hydroxylated) from the parent compound ondansetron.
- the theoretical ⁇ values were calculated for the ten lowest energy structures.
- the theoretical ⁇ values are averaged for each structure (metabolite) therefore, providing an average of all the possible structures/conformations the molecule may "explore" in the gas-phase.
- the travelling wave ⁇ values are then derived for each metabolite by LC-IMS-MS.
- Synapt G2 HDMS was operated in ion mobility mode (positive ion mode).
- the travelling wave ion mobility cell was operated at a pressure of 3.3 mbar nitrogen with a travelling wave speed of 1000 m/sec and a pulse amplitude of 40 V.
- Acquity UPLC was operated at a flow rate of 600ul_/min.
- Solvent A was 5 mM ammonium acetate.
- Solvent B was 1 00% acetonitrile.
- Chromatographic separation was performed on a 2.1 x 50 mm BEH C1 8 1 .7 um column utilising a 5 minute gradient, starting at 5% solvent B, increasing to 30% solvent B, over 3 minutes.
- Simplified Molecular Input Line Entry Specification (SMILES) strings for all the compounds were transformed into 3D structures utilizing the software package Molecular Operating Environment (MOE). The following steps were performed: First, wash smiles. This checks the smiles string and protonates/deprotonates ionisation centers according to physiological pH . Second, generate 3-D structures (converts 2-D into 3-D). Third, energy minimise the 3-D structures. Once the 3D compounds were generated, the conformations were explored with the "Systematic Search" protocol in MOE. This is suitable since ondansetron only has two rotatable bonds (3 if you count the -OH in the metabolites) so, the combinatorial explosion is controlled.
- MOE Molecular Operating Environment
- Theoretical ⁇ values were calculated using the open source software program MOBCALTM. Briefly, the MOBCALTM output is based on three different models/algorithms which calculate the theoretical ⁇ of any molecule possessing a 3- dimensional coordinate file. The models/algorithms are the Projection Approximation (PA), Exact Hard Sphere Scattering (EHSS) and the Trajectory Method (TM), however, only the MOBCALTM Trajectory Method ⁇ value will be reported. The WatersTM ⁇ algorithm, based in DriftScope 2.1 was also utilized. The values obtained from the MOBCALTM and the WatersTM ⁇ algorithm were then compared to the travelling wave derived ⁇ values.
- PA Projection Approximation
- EHSS Exact Hard Sphere Scattering
- TM Trajectory Method
- Travelling wave ion mobility calibration of the Synapt G2 HDMS instrument was undertaken using five singly charged polyalanine ions (AAA, AAAA, AAAAA, AAAAAA and AAAAAAA) whose absolute ⁇ values are known from standard drift- tube studies.
- the polyalanine ions and their respective absolute ⁇ used for T-Wave ion mobility calibration were as follows: - AAA (m/z 232.1 ), 89.0 A 2 , AAAA (m/z 303.1 ), 102.9 A 2 , AAAAA (m/z 374.1 ) 1 15.0 A 2 , AAAAAA (m/z 445.2) 127.0 A 2 and AAAAAAA (m/z 516.2) 140.5 A 2 .
- the drift tube ⁇ values were corrected by multiplying by the square root of the reduced mass and dividing by the charge state of the calibrant species to provide a term related directly to 1 /K.
- the measured travelling wave drift times were finally corrected for their mass-independent and mass-dependent flight times between the travelling wave ion mobility cell and the ToF analyser.
- the derived coefficients of a and b obtained from the calibration curve were used to calculate the ⁇ value of ondansetron and related hydroxylated metabolites from their measured travelling wave drift times obtained under identical operating conditions to those of the calibrant species polyalanine.
- Figure 5 shows the MS/MS spectra of hydroxylated metabolite GR60661 .
- the fragmentation spectra for metabolites GR63418 and GR90315 are identical to that presented in Figure 5. Therefore, one cannot identify the hydroxylated metabolites based purely on their MS/MS fragmentation pattern.
- the present invention provides a new approach to metabolite structural elucidation and characterisation by generating an ensemble of energy minimised metabolite structures in-silico as described above in the Molecular Modelling section. This allows for the theoretical ⁇ calculation of the metabolite structures, which can then be compared to the travelling wave ion mobility measured ⁇ value, based on the LC- IMS-MS experiment.
- Figure 6 shows the comparison of the MOBCALTM Trajectory Method, the WatersTM PA algorithm derived ⁇ values for the ten lowest energy structures of ondansetron, GR60661 , GR63418 and GR90315 and a comparison to the travelling wave derived ⁇ values.
- Potential candidate structures can also be created and optimised for fragment ions.
- the metabolite fragment ion m/z 228.1 corresponding to the methylated imidazole neutral loss species was selected due to the absence of any rotatable bonds and almost planar structure.
- ⁇ value GR60661 and GR63418, it is also possible to differentiate them based on the fragment ion (m/z 228.1 ) conformation, both T-wave and theoretically derived.
- FIG. 8 shows the LC-IMS-MS data, with annotated travelling wave, displayed in DriftScope 2.1 , as a function of drift time verses chromatographic retention time.
- the invention allows travelling wave ion mobility combined with in-silico molecular modelling of hydroxylated metabolite structures can be successfully implemented to positively identify the hydroxylated metabolites based on their collisional cross- sections.
- This technique can be rapidly utilised to supplement all traditional methodologies of metabolite identification, such as MS/MS.
- This IM-MS based approach can enable a researcher to identify metabolites without the need for metabolite synthesis and/or NMR structural elucidation.
- this method can not identify all metabolites present, since some metabolites will possess an identical collision cross-sections (both theoretically and instrument ion mobility derived), this method can be used to "rule-out" potential structures on the basis of their theoretical Q:T-Wave derived ⁇ comparison.
Landscapes
- Physics & Mathematics (AREA)
- Chemical & Material Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Health & Medical Sciences (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Engineering & Computer Science (AREA)
- Computing Systems (AREA)
- Theoretical Computer Science (AREA)
- General Health & Medical Sciences (AREA)
- Analytical Chemistry (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Bioinformatics & Computational Biology (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Biochemistry (AREA)
- General Physics & Mathematics (AREA)
- Chemical Kinetics & Catalysis (AREA)
- Electrochemistry (AREA)
- Crystallography & Structural Chemistry (AREA)
- Molecular Biology (AREA)
- Toxicology (AREA)
- Other Investigation Or Analysis Of Materials By Electrical Means (AREA)
Abstract
A method of and apparatus for identifying and/or characterising a sample that may incorporate two or more isomeric or isobaric compounds such as hydroxylated metabolites. The method involves modelling an ensemble of possible structures for each of two or more known isomeric or isobaric compounds, calculating a theoretical collision cross-section for each modelled structure, averaging the calculated values for each known compound to provide a theoretical collision cross-section value for each known compound. A travelling wave ion mobility cell is used to measuring a collision cross-section value for the sample compound and the measured value is then compared with the theoretical values to identify which of the two or more known compounds the sample compound most closely resembles.
Description
ION MOBILITY METHOD AND APPARATUS FOR IDENTIFYING
AND/OR CHARACTERISING A SAMPLE COMPOUND
This invention relates generally to a method and apparatus for identifying and/or characterising a sample compound. More specifically, although not exclusively, this invention relates to such a method for differentiating isomeric or isobaric compounds.
As an analytical technique, mass spectrometry (MS) has limited capabilities for separating isomeric or isobaric species or for providing structural conformational information. Some structural information can be gained through tandem MS (MS/MS) techniques. There are classes of substituted molecules, such as hydroxylated metabolites which produce identical MS/MS spectra. This is problematic if the exact position of the hydroxylation is required. Ion mobility (IM) has the ability to separate isomeric or isobaric species, such as hydroxylated metabolites, rapidly (msec) based on differences in their collision cross-section (Ω; physical size, and shape) in the gas-phase, thus providing specific information on ionic configuration and therefore, the position of the hydroxyl moiety. The combination of MS with IM provides an extremely powerful analytical tool.
In recent years, the ability to measure a protein or protein/protein complex's collision cross-section (Ω) value has proven to be a very powerful feature and has significantly added to the advance of structural biology. This however, cannot be said about small molecule and metabolite identification. There has been very little addition to the literature in this area of research. However, recent research has begun to investigate, successfully, the separative powers of travelling wave ion mobility for small ruthenium anticancer structural isomers or isobars.
The present invention provides methods and apparatus that are particularly suited for identification of such compounds where specific uses of these advances are made. More specifically, the methods and apparatus of the present invention enable more accurate characterisation and/or identification of such compounds.
One aspect of the invention provides a method of identifying a sample compound, the method comprising the steps of: (i) proving models of an ensemble of possible or expected structures for each of two or more known compounds; (ii) calculating a theoretical collision cross-section for each modelled structure; (iii) averaging the
calculated theoretical collision cross-sections for each known compound to provide a theoretical collision cross-section value for each known compound; (iv) measuring a collision cross-section value for the sample compound using an ion mobility technique or cell; and (v) comparing the measured collision cross-section value with the calculated theoretical collision cross-section values to identify which of the two or more known compounds the sample compound most closely resembles.
A second aspect of the invention provides a method of characterising the structure of a sample compound, the method comprising the steps of: (i) providing models of an ensemble of possible or expected structures for each of two or more known compounds; (ii) calculating a theoretical collision cross-section for each modelled structure; (iii) averaging the calculated theoretical collision cross-sections for each known compound to provide a theoretical collision cross-section value for each known compound; (iv) measuring a collision cross-section value for the sample compound using an ion mobility technique or cell; and (v) comparing the measured collision cross-section value with the theoretical collision cross-section values to characterise the structure of a sample compound.
The two or more known compounds may comprise two or more known isomeric or isobaric compounds, for example metabolites, e.g. hydroxylated metabolites.
The ion mobility technique or cell preferably comprises a travelling wave ion mobility technique or cell.
The method may further comprise the step of analysing a sample compound, for example a further sample compound, that may be from the same sample, e.g. using a mass spectrometer, for example to produce a mass spectrum. This analysing step may be carried out before the comparing step, e.g. in order to identify and/or further characterise the sample compound.
The step of providing models may further comprise checking simplified molecular input line entry specification strings and/or protonates / deprotonates ionisation centres, for example according to physiological pH. Additionally or alternatively, the step of providing models may comprise generating three dimensional structures and may further comprise energy minimising the three dimensional structures.
Another aspect of the invention provides a computer program element, which may comprise computer readable program code means, e.g. for causing a processor to execute a procedure, for example to implement one or more of the above steps. The computer program element may be embodied on a computer readable medium.
A further aspect of the invention provides a computer readable medium, preferably having a program stored thereon, for example where the program is to make a computer execute a procedure, e.g. to implement one or more of the above steps.
A yet further aspect of the invention provides an analytical apparatus suitable for carrying out a method as described above and/or comprising a program element as described above and/or a computer readable medium as described above, which apparatus may be specifically adapted for carrying out the method. The analytical instrument may comprise a mass spectrometer.
A yet further aspect of the invention provides a retrofit kit for adapting an analytical apparatus or mass spectrometer to provide an apparatus or mass spectrometer as described above. The kit may comprise a program element as described above and/or a computer readable medium as described above.
Embodiments of the invention will now be described by way of example only with reference to the accompanying drawings in which:
Figure 1 illustrates a parent drug ondansetron and three hydroxylated metabolites GR60661 , GR90315 and GR63418;
Figure 2 illustrates the ensemble of hydroxylated metabolite GR60661 structures created by molecular dynamics and structure refinement, wherein the structures differ by degree rotation around the two rotatable bonds;
Figure 3 is a graph showing a travelling wave ion mobility calibration curve of drift-time (msec) verses adjusted collision cross-section for five polyalanine ions, wherein a power curve of the form y=axb is fitted to the data;
Figure 4 illustrates the UPLC/MS separation of ondansetron and hydroxylated metabolites GR60661 , GR63418 and GR90315;
Figure 5 illustrates the collision induced dissociation MS/MS spectrum for the hydroxylated metabolite GR60661 m/z 310.1 ;
Figure 6 is a table showing the comparison of the collision cross-section (Ω, in A2) values of ondansetron and hydroxylated metabolites GR60661 , GR63418 and GR90315 obtained using travelling wave ion mobility measurements, Waters™ Projection Approximation (PA) and the MOBCAL™ Trajectory Method (TM) calculations;
Figure 7A is an overlay of ion mobility arrival time distributions for the parent compounds ondansteron and metabolites GR90315, GR63418 and GR60661 , in order of ion mobility elution (travelling wave 1 100m/sec, 3.3 mbar nitrogen);
Figure 7B is an overlay of ion mobility arrival time distributions for metabolites GR63418 and GR60661 (travelling wave 1 100m/sec, 3.3 mbar nitrogen);
Figure 8 shows a DriftScope 2.1 plot of drift time (msec) verses retention time (mins) of the LC-IMS-MS analysis of ondansetron, GR60661 , GR63418 and GR90315 and assigned travelling wave collision cross-section values; and
Figure 9 is a table showing the comparison of the collision cross-section (Ω, in A2) values of methylated imidazole neutral loss fragment of ondansetron and hydroxylated metabolites GR60661 , GR63418 and GR90315 obtained using travelling wave ion mobility measurements and the MOBCAL™ Trajectory Method calculations.
Experiments were performed using a hybrid quadrupole/ion mobility/orthogonal acceleration time-of-flight Synapt G2 HDMS instrument. Ion mobility separation was performed at a pressure of 3.2 mbar nitrogen with a wave velocity of 1000 m/sec and a fixed travelling wave amplitude of 40 V was used. Theoretical Ω values were calculated using MOBCAL™ and the Waters™ Ω algorithm (as described in WO2010/1 19289 and WO2010/1 1923, the contents of which are incorporated herein by reference) and compared to the travelling wave derived Ω values. Ion mobility calibration was undertaken using five singly charged polyalanine ions. Three- dimensional metabolite conformations were explored with the "Systematic Search" protocol in MOE. Ten thousand combinations of possible values for variables were
explored. All conformations were minimised with the MMFF94 forcefield. The ten conformations with the lowest energy were selected for each metabolite.
The Ω value of the parent compound ondansetron and the three hydroxylated metabolites GR60661 , GR63418 and GR90315 were analysed by UPLC-IMS-MS. All three hydroxylated metabolites showed different chromatographic elution time. However, the MS and MS/MS spectra are all identical, thus hampering identification. The theoretically derived Ω values for hydroxylated metabolites GR60661 , GR63418 and GR90315 are 1 12.5 A2, 1 1 1 .8 A2 and 1 10.9 A2 respectively. All reported values based on the Projection Approximation and Trajectory Method. Travelling wave ion mobility derived collision cross-section values are in very close agreements with the theoretically derived values, thus allowing accurate identification of the hydroxylated metabolites, even those with a collision cross-section difference of less than 1 A2. No additional structural characterization techniques, such as NMR or X-ray crystallography were required for metabolite structural determination. This is a very time efficient and effective means of identifying drug metabolites directly from biological matrices. Additionally it can be utilised in a situation where chemically synthesised standards may not be available.
The invention provides a new methodology which combines molecular modelling, theoretical Ω calculations and travelling wave ion mobility derived Ω measurements, to successfully indentify a number of drug metabolites, which cannot be differentiated by traditional MS/MS techniques. In-silico molecular modelling is used to create an ensemble (approximately 10,000) of possible metabolite structures (in this case study hydroxylated) from the parent compound ondansetron. The theoretical Ω values were calculated for the ten lowest energy structures. The theoretical Ω values are averaged for each structure (metabolite) therefore, providing an average of all the possible structures/conformations the molecule may "explore" in the gas-phase. The travelling wave Ω values are then derived for each metabolite by LC-IMS-MS.
Instrumentation
Synapt G2 HDMS was operated in ion mobility mode (positive ion mode). The travelling wave ion mobility cell was operated at a pressure of 3.3 mbar nitrogen with a travelling wave speed of 1000 m/sec and a pulse amplitude of 40 V.
Acquity UPLC was operated at a flow rate of 600ul_/min. Solvent A was 5 mM ammonium acetate. Solvent B was 1 00% acetonitrile. Chromatographic separation was performed on a 2.1 x 50 mm BEH C1 8 1 .7 um column utilising a 5 minute gradient, starting at 5% solvent B, increasing to 30% solvent B, over 3 minutes.
Molecular Modelling
Simplified Molecular Input Line Entry Specification (SMILES) strings for all the compounds were transformed into 3D structures utilizing the software package Molecular Operating Environment (MOE). The following steps were performed: First, wash smiles. This checks the smiles string and protonates/deprotonates ionisation centers according to physiological pH . Second, generate 3-D structures (converts 2-D into 3-D). Third, energy minimise the 3-D structures. Once the 3D compounds were generated, the conformations were explored with the "Systematic Search" protocol in MOE. This is suitable since ondansetron only has two rotatable bonds (3 if you count the -OH in the metabolites) so, the combinatorial explosion is controlled.
Possible conformations for the cyclohexane ring were also considered. In total, around 1 0,000 combinations of possible values for the variables explored were generated. All the conformations were energy minimised with the MMFF94 forcefield and the ten conformations (Figure 2) with the lowest energy were selected for each metabolite.
Energy minimized candidate structures were also created for the fragment ions relating to the neutral loss of the methylated imidazole structure (m/z 228.1 , C1402NH13, for hydroxylated metabolites GR60661 , GR9031 5, GR63418 and m/z 212.1 , Ci4ONHi3, for ondansetron). The basis set HF/6-31 G++ was used for structure refinement.
Theoretical Collision Cross-section Calculations
Theoretical Ω values were calculated using the open source software program MOBCAL™. Briefly, the MOBCAL™ output is based on three different models/algorithms which calculate the theoretical Ω of any molecule possessing a 3- dimensional coordinate file. The models/algorithms are the Projection Approximation (PA), Exact Hard Sphere Scattering (EHSS) and the Trajectory Method (TM),
however, only the MOBCAL™ Trajectory Method Ω value will be reported. The Waters™ Ω algorithm, based in DriftScope 2.1 was also utilized. The values obtained from the MOBCAL™ and the Waters™ Ω algorithm were then compared to the travelling wave derived Ω values.
Travelling Wave Ion Mobility Calibration
Travelling wave ion mobility calibration of the Synapt G2 HDMS instrument was undertaken using five singly charged polyalanine ions (AAA, AAAA, AAAAA, AAAAAA and AAAAAAA) whose absolute Ω values are known from standard drift- tube studies. The polyalanine ions and their respective absolute Ωε used for T-Wave ion mobility calibration were as follows: - AAA (m/z 232.1 ), 89.0 A2, AAAA (m/z 303.1 ), 102.9 A2, AAAAA (m/z 374.1 ) 1 15.0 A2, AAAAAA (m/z 445.2) 127.0 A2 and AAAAAAA (m/z 516.2) 140.5 A2.
For the travelling wave calibration, the drift tube Ω values were corrected by multiplying by the square root of the reduced mass and dividing by the charge state of the calibrant species to provide a term related directly to 1 /K. The measured travelling wave drift times were finally corrected for their mass-independent and mass-dependent flight times between the travelling wave ion mobility cell and the ToF analyser. The adjusted Ω values were plotted against corrected drift time (msec) values determined using the T-wave device and the data fitted using an empirically determined power-law form, y = axb (Figure x). The derived coefficients of a and b obtained from the calibration curve were used to calculate the Ω value of ondansetron and related hydroxylated metabolites from their measured travelling wave drift times obtained under identical operating conditions to those of the calibrant species polyalanine.
It was found that the three hydroxylated metabolites can all be separated by UPLC. However, in the absence of any chemical standard as a reference, and one was analysing metabolites directly from a biological matrix, it would be very challenging to identify each metabolite based on the MS or the MS/MS spectra alone. Figure 5 shows the MS/MS spectra of hydroxylated metabolite GR60661 . The fragmentation spectra for metabolites GR63418 and GR90315 are identical to that presented in Figure 5. Therefore, one cannot identify the hydroxylated metabolites based purely on their MS/MS fragmentation pattern.
The present invention provides a new approach to metabolite structural elucidation and characterisation by generating an ensemble of energy minimised metabolite structures in-silico as described above in the Molecular Modelling section. This allows for the theoretical Ω calculation of the metabolite structures, which can then be compared to the travelling wave ion mobility measured Ω value, based on the LC- IMS-MS experiment. Figure 6 shows the comparison of the MOBCAL™ Trajectory Method, the Waters™ PA algorithm derived Ω values for the ten lowest energy structures of ondansetron, GR60661 , GR63418 and GR90315 and a comparison to the travelling wave derived Ω values.
Observed in Figure 6 is the substantial agreement of the theoretically derived Ω values and the travelling wave derived Ω values. The trend of the theoretical and travelling wave derived Ω values also match, in that metabolite GR60661 possesses the largest Ω value and GR90315, the smallest Ω value.
This demonstrates that it is possible with the methods of the present invention to generate accurate theoretical Ω values for metabolites, compare them to gas-phase Ω measurements and use this comparison to assign and identify the position of the hydroxylation site within the metabolite. In the example shown herein, structures GR60661 and GR63418 only differ in theoretical Ω value of 0.2 A2 (1 1 1 .4 and 1 1 1 .2 respectively, based on the PA method).
This small difference in theoretical collision cross-section can be measured and characterised using the Synapt G2 instrument. This is demonstrated in Figure 7A & 7B where the arrival time distributions of all metabolites have been overlaid. The small arrival time distribution difference of metabolites GR60661 and GR63418 is reproducible and measureable (Figure 7B). Additionally, the addition of the hydroxyl moiety results in a measurable increase in collision cross-section of 3-4 A2, depending on which metabolite is being considered.
IMS-MS of fragment ions
Potential candidate structures can also be created and optimised for fragment ions. The metabolite fragment ion m/z 228.1 , corresponding to the methylated imidazole neutral loss species was selected due to the absence of any rotatable bonds and almost planar structure. Considering the two closest structurally related metabolites, in terms of Ω value, GR60661 and GR63418, it is also possible to differentiate them based on the fragment ion (m/z 228.1 ) conformation, both T-wave and theoretically
derived. By first considering the T-wave derived Ω values, then the arrival time distributions are reproducibly offset by a value of "l OOusec (equating to a 0.1 A2 difference) with GR63418 always preceding GR60661 (Figure 9). The trend of the MOBCAL™ trajectory method derived Ω values follows that of the instrument derived values, albeit, a larger difference of 1 .0 A2 between GR63418 and GR60661 .
Based on this limited sample set, it would be challenging to predict the position of the hydroxyl moiety of the meta (GR63418) and para-hydroxylated (GR60661 ) metabolites by comparing theoretical and T-wave derived Ω values. The situation is more clear when considering the ortho-hydroxylated metabolite (GR90315). The T- wave and theoretically derived values correlate very well (Figure 9).
Using carbon dioxide as the ion mobility drift gas
As demonstrated previously by Hill (Anal Chem 2000) the ATDs of compounds under investigation can change dramatically depending on which drift gases are used. Herein, it is demonstrated that using the more polarizable carbon dioxide as the drift gas, the metabolite GR63418 has a longer drift time than GR60661 . This is opposite to what is observed when using the less polar nitrogen as the drift gas.
Benefits
The technique of the present invention is rapid, since it does not require lengthy chemical synthesis pathways for production of metabolites standards. Utilising LC- IMS-MS, it is possible to determine the Ω values of the metabolites directly from biological matrices. Additionally, and importantly, this method does not require skilled interpretation of NMR spectra to derive metabolite structures. Figure 8 shows the LC-IMS-MS data, with annotated travelling wave, displayed in DriftScope 2.1 , as a function of drift time verses chromatographic retention time.
The invention allows travelling wave ion mobility combined with in-silico molecular modelling of hydroxylated metabolite structures can be successfully implemented to positively identify the hydroxylated metabolites based on their collisional cross- sections. This technique can be rapidly utilised to supplement all traditional methodologies of metabolite identification, such as MS/MS.
This IM-MS based approach can enable a researcher to identify metabolites without the need for metabolite synthesis and/or NMR structural elucidation. In cases where this method can not identify all metabolites present, since some metabolites will possess an identical collision cross-sections (both theoretically and instrument ion mobility derived), this method can be used to "rule-out" potential structures on the basis of their theoretical Q:T-Wave derived Ω comparison.
It will be appreciated by those skilled in the art that any number of combinations of the aforementioned features and/or those shown in the appended drawings provide clear advantages over the prior art and are therefore within the scope of the invention described herein.
Claims
1 . A method of identifying a sample compound, the method comprising the steps of:
i. providing models of an ensemble of possible structures for each of two or more known compounds;
ii. calculating a theoretical collision cross-section for each modelled structure;
iii. averaging the calculated theoretical collision cross-sections for each known compound to provide a theoretical collision cross-section value for each known compound;
iv. measuring a collision cross-section value for the sample compound using an ion mobility technique; and
v. comparing the measured collision cross-section value with the theoretical collision cross-section values to identify which of the two or more known compounds the sample compound most closely resembles.
2. A method of characterising the structure of a sample compound, the method comprising the steps of:
i. providing models of an ensemble of possible structures for each of two or more known compounds;
ii. calculating a theoretical collision cross-section for each modelled structure;
iii. averaging the calculated theoretical collision cross-sections for each known compound to provide a theoretical collision cross-section value for each known compound;
iv. measuring a collision cross-section value for the sample compound using an ion mobility technique; and
v. comparing the measured collision cross-section value with the theoretical collision cross-section values to characterise the structure of a sample compound.
3. Method according to claim 1 or claim 2, wherein the two or more known compounds comprise two or more known isomeric or isobaric compounds.
4. Method according to any preceding claim, wherein the ion mobility technique comprises a travelling wave ion mobility technique.
5. Method according to any preceding claim further comprising the step of analysing a further sample compound from the same sample using a mass spectrometer to produce a mass spectrum before the comparing step in order to identify or further characterise the sample compound.
6. Method according to any preceding claim, wherein the step of providing models further comprises checking simplified molecular input line entry specification strings and protonates / deprotonates ionisation centres according to physiological pH.
7. Method according to any preceding claim, wherein the step of providing models further comprises generating three dimensional structures.
8. Method according to claim 7, wherein the step of providing models further comprises energy minimising the three dimensional structures.
9. Method according to any preceding claim, wherein the isomeric or isobaric compound structures comprise hydroxylated metabolites.
10. A computer program element comprising computer readable program code means for causing a processor to execute a procedure to implement the method of any one of claims 1 to 9.
1 1 . A computer program element according to claim 10 embodied on a computer readable medium.
12. A computer readable medium having a program stored thereon, where the program is to make a computer execute a procedure to implement the method of any one of claims 1 to 9.
13. An analytical apparatus suitable for carrying out a method according to any one of claims 1 to 9 or comprising a program element according to claim 10 or claim 1 1 or a computer readable medium according to claim 12.
14. An analytical apparatus according to claim 13, wherein the apparatus is specifically adapted for carrying out a method according to any one of claims 1 to 9.
15. A mass spectrometer suitable for carrying out a method according to any one of claims 1 to 9 or comprising a program element according to claim 10 or claim 1 1 or a computer readable medium according to claim 12.
16. A mass spectrometer according to claim 15, wherein the apparatus is specifically adapted for carrying out a method according to any one of claims 1 to 9.
17. A retrofit kit for adapting an analytical apparatus to provide an apparatus according to claim 13 or claim 14, the kit comprising a program element according to claim 10 or claim 1 1 or a computer readable medium according to claim 12.
18. A retrofit kit for adapting a mass spectrometer to provide a mass spectrometer according to claim 15 or claim 16, the kit comprising a program element according to claim 10 or claim 1 1 or a computer readable medium according to claim 12.
Priority Applications (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
EP11724438.4A EP2558850B1 (en) | 2010-04-15 | 2011-04-15 | Ion mobility method and apparatus for identifying a sample compound |
US13/640,848 US20130218478A1 (en) | 2010-04-15 | 2011-04-15 | Method And System Of Identifying A Sample By Analysing A Mass Spectrum By The Use Of A Bayesian Inference Technique |
Applications Claiming Priority (6)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
GBGB1006312.1A GB201006312D0 (en) | 2010-04-15 | 2010-04-15 | Method and apparatus for identifying and/or characterising a sample compound |
GB1006312.1 | 2010-04-15 | ||
US32784310P | 2010-04-26 | 2010-04-26 | |
US61/327,843 | 2010-04-26 | ||
GB1008414.3 | 2010-05-20 | ||
GBGB1008414.3A GB201008414D0 (en) | 2010-05-20 | 2010-05-20 | Sites of metabolic substitution; investigating metobolite structures utilising ion mobility and molecular modelling |
Publications (1)
Publication Number | Publication Date |
---|---|
WO2011128703A1 true WO2011128703A1 (en) | 2011-10-20 |
Family
ID=44246954
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
PCT/GB2011/050756 WO2011128703A1 (en) | 2010-04-15 | 2011-04-15 | Ion mobility method and apparatus for identifying and/or characterising a sample compound |
Country Status (3)
Country | Link |
---|---|
US (1) | US20130218478A1 (en) |
EP (1) | EP2558850B1 (en) |
WO (1) | WO2011128703A1 (en) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2015136273A1 (en) * | 2014-03-10 | 2015-09-17 | Micromass Uk Limited | Using theoretical collision cross section ("ccs") in sample identification |
GB2527632A (en) * | 2014-03-10 | 2015-12-30 | Micromass Ltd | Unknown identification using theroretical collision cross section |
US20170200594A1 (en) * | 2014-05-30 | 2017-07-13 | Micromass Uk Limited | Hybrid Mass Spectrometer |
EP3479110A4 (en) * | 2016-06-30 | 2020-02-26 | Waters Technologies Corporation | Chemical class compositions from collision cross-section fragment ions |
DE102022112760A1 (en) | 2021-06-02 | 2022-12-08 | Bruker Scientific Llc | EVALUATION USING PHYSICAL-CHEMICAL PROPERTIES FOR STRUCTURE IDENTIFICATION IN ION SPECTROMETRY |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US10115576B2 (en) | 2013-12-12 | 2018-10-30 | Waters Technologies Corporation | Method and an apparatus for analyzing a complex sample |
CN108918645B (en) * | 2018-06-29 | 2021-10-15 | 广州禾信仪器股份有限公司 | Isomeride body spectrum obtaining method and isomeride identification method |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2006040622A1 (en) * | 2004-10-14 | 2006-04-20 | Highchem, Ltd. | Determination of molecular structures using tandem mass spectrometry |
WO2010011923A1 (en) | 2008-07-24 | 2010-01-28 | Gesturetek, Inc. | Enhanced detection of circular engagement gesture |
WO2010119289A2 (en) | 2009-04-15 | 2010-10-21 | Micromass Uk Limited | A method and system of estimating the cross-sectional area of a molecule for use in the prediction of ion mobility |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8036867B2 (en) * | 2003-10-14 | 2011-10-11 | Verseon | Method and apparatus for analysis of molecular configurations and combinations |
-
2011
- 2011-04-15 WO PCT/GB2011/050756 patent/WO2011128703A1/en active Application Filing
- 2011-04-15 EP EP11724438.4A patent/EP2558850B1/en active Active
- 2011-04-15 US US13/640,848 patent/US20130218478A1/en not_active Abandoned
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2006040622A1 (en) * | 2004-10-14 | 2006-04-20 | Highchem, Ltd. | Determination of molecular structures using tandem mass spectrometry |
WO2010011923A1 (en) | 2008-07-24 | 2010-01-28 | Gesturetek, Inc. | Enhanced detection of circular engagement gesture |
WO2010119289A2 (en) | 2009-04-15 | 2010-10-21 | Micromass Uk Limited | A method and system of estimating the cross-sectional area of a molecule for use in the prediction of ion mobility |
Non-Patent Citations (3)
Title |
---|
HILL, ANAL CHEM, 2000 |
KNAPMAN T W ET AL: "Considerations in experimental and theoretical collision cross-section measurements of small molecules using travelling wave ion mobility spectrometry-mass spectrometry", INTERNATIONAL JOURNAL OF MASS SPECTROMETRY, ELSEVIER SCIENCE PUBLISHERS, AMSTERDAM, NL, vol. 298, no. 1-3, 4 October 2009 (2009-10-04), pages 17 - 23, XP027492678, ISSN: 1387-3806, [retrieved on 20091004], DOI: DOI:10.1016/J.IJMS.2009.09.011 * |
TAO ET AL: "A Collision Cross-Section Database of Singly-Charged Peptide Ions", JOURNAL OF THE AMERICAN SOCIETY FOR MASS SPECTROMETRY, ELSEVIER SCIENCE INC, US, vol. 18, no. 7, 26 June 2007 (2007-06-26), pages 1232 - 1238, XP022131654, ISSN: 1044-0305, DOI: DOI:10.1016/J.JASMS.2007.04.003 * |
Cited By (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
GB2527632B (en) * | 2014-03-10 | 2018-05-30 | Micromass Ltd | Unknown identification using theroretical collision cross section |
WO2015136272A1 (en) * | 2014-03-10 | 2015-09-17 | Micromass Uk Limited | Unknown identification using theoretical collision cross section |
WO2015136274A1 (en) * | 2014-03-10 | 2015-09-17 | Micromass Uk Limited | Theoretical collision cross section ("ccs") in experimental design |
GB2527632A (en) * | 2014-03-10 | 2015-12-30 | Micromass Ltd | Unknown identification using theroretical collision cross section |
GB2527879A (en) * | 2014-03-10 | 2016-01-06 | Micromass Ltd | Using theoretical collision cross section ("CCS") in sample identification |
WO2015136273A1 (en) * | 2014-03-10 | 2015-09-17 | Micromass Uk Limited | Using theoretical collision cross section ("ccs") in sample identification |
GB2527879B (en) * | 2014-03-10 | 2018-09-12 | Micromass Ltd | Using theoretical collision cross section ("CCS") in sample identification |
US20170200594A1 (en) * | 2014-05-30 | 2017-07-13 | Micromass Uk Limited | Hybrid Mass Spectrometer |
US10600627B2 (en) * | 2014-05-30 | 2020-03-24 | Micromass Uk Limited | Hybrid mass spectrometer |
US10886116B2 (en) * | 2014-05-30 | 2021-01-05 | Micromass Uk Limited | Hybrid mass spectrometer |
EP3479110A4 (en) * | 2016-06-30 | 2020-02-26 | Waters Technologies Corporation | Chemical class compositions from collision cross-section fragment ions |
US11333628B2 (en) | 2016-06-30 | 2022-05-17 | Waters Technologies Corporation | Chemical class compositions from collision cross-section fragment ions |
DE102022112760A1 (en) | 2021-06-02 | 2022-12-08 | Bruker Scientific Llc | EVALUATION USING PHYSICAL-CHEMICAL PROPERTIES FOR STRUCTURE IDENTIFICATION IN ION SPECTROMETRY |
Also Published As
Publication number | Publication date |
---|---|
US20130218478A1 (en) | 2013-08-22 |
EP2558850A1 (en) | 2013-02-20 |
EP2558850B1 (en) | 2017-03-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
EP2558850B1 (en) | Ion mobility method and apparatus for identifying a sample compound | |
Scheubert et al. | Computational mass spectrometry for small molecules | |
US11573201B2 (en) | Method and system for identifying structure of compound | |
Wang et al. | Mass spectrometry-based protein identification by integrating de novo sequencing with database searching | |
CN108508078B (en) | Method for identifying elemental composition of molecular species | |
Mullard et al. | A new strategy for MS/MS data acquisition applying multiple data dependent experiments on Orbitrap mass spectrometers in non-targeted metabolomic applications | |
CN107807198A (en) | For the method for the single isotopic mass for identifying various molecules | |
Ieritano et al. | How hot are your ions in differential mobility spectrometry? | |
Tugizimana et al. | Mass spectrometry in untargeted liquid chromatography/mass spectrometry metabolomics: Electrospray ionisation parameters and global coverage of the metabolome | |
Burns et al. | Platform dependencies in bottom-up hydrogen/deuterium exchange mass spectrometry | |
WO2013104004A1 (en) | Comprehensive interference treatment for icp-ms analysis | |
US20190018928A1 (en) | Methods for Mass Spectrometry-Based Structure Determination of Biomacromolecules | |
Donohoe et al. | Ion mobility spectrometry-hydrogen deuterium exchange mass spectrometry of anions: Part 1. Peptides to proteins | |
Kwan et al. | Molecular characterization of the surface excess charge layer in droplets | |
Ieritano et al. | The charge-state and structural stability of peptides conferred by microsolvating environments in differential mobility spectrometry | |
Ham | Even electron mass spectrometry with biomolecule applications | |
US20230251224A1 (en) | Method and system for identifying structure of compound | |
US11094399B2 (en) | Method, system and program for analyzing mass spectrometoric data | |
Olson et al. | Production of reliable MALDI spectra with quality threshold clustering of replicates | |
Yates III | Mass spectrometry as an emerging tool for systems biology | |
EP3523818B1 (en) | System and method for real-time isotope identification | |
Consta et al. | Advances in modeling the stability of noncovalent complexes in charged droplets with applications in electrospray ionization-ms experiments | |
WO2019175568A1 (en) | Methods and systems for analysis | |
Zhou | Computational analysis of LC-MS/MS data for metabolite identification | |
Knyazev et al. | Monte Carlo/RRKM/Classical Trajectories Modeling of Collisional Excitation and Dissociation of n-Butylbenzene Ion in Multipole Collision Cells of Tandem Mass Spectrometers |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
121 | Ep: the epo has been informed by wipo that ep was designated in this application |
Ref document number: 11724438 Country of ref document: EP Kind code of ref document: A1 |
|
REEP | Request for entry into the european phase |
Ref document number: 2011724438 Country of ref document: EP |
|
WWE | Wipo information: entry into national phase |
Ref document number: 2011724438 Country of ref document: EP |
|
NENP | Non-entry into the national phase |
Ref country code: DE |
|
WWE | Wipo information: entry into national phase |
Ref document number: 13640848 Country of ref document: US |