EP1573308A1 - Systeme et procede de definition de profils dielectrophoretiques de particules - Google Patents

Systeme et procede de definition de profils dielectrophoretiques de particules

Info

Publication number
EP1573308A1
EP1573308A1 EP02794251A EP02794251A EP1573308A1 EP 1573308 A1 EP1573308 A1 EP 1573308A1 EP 02794251 A EP02794251 A EP 02794251A EP 02794251 A EP02794251 A EP 02794251A EP 1573308 A1 EP1573308 A1 EP 1573308A1
Authority
EP
European Patent Office
Prior art keywords
particles
velocity
frequency
particle
electric field
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.)
Withdrawn
Application number
EP02794251A
Other languages
German (de)
English (en)
Inventor
Richard Stanley Lee
Ron Pethig
Mark Stewart Talary
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.)
Aura Biosystems Inc
Original Assignee
Aura Biosystems Inc
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 Aura Biosystems Inc filed Critical Aura Biosystems Inc
Priority claimed from PCT/US2002/039931 external-priority patent/WO2004055505A1/fr
Publication of EP1573308A1 publication Critical patent/EP1573308A1/fr
Withdrawn legal-status Critical Current

Links

Definitions

  • This invention relates to dielectrophoresis.
  • Dielectrophoresis refers to the force experienced by particles suspended in a fluid medium in applied electric field gradients. Due to the electric field gradient, differences in the dielectric polarization between the particles and the fluid medium cause the particles to experience a dielectrophoretic force. This effect can be quantified in terms of the electromagnetic momentum balance via the Maxwell stress tensor, or in terms of the magnitude and distribution of the charges induced on and within the particle by the applied field. Particles, such as blood cells, experiencing strong DEP motion will typically experience a DEP force of about 10 "11 N, which is about 40 times greater than the gravitational settling force and about 2 x 10 5 times larger than the maximum Brownian diffusion force.
  • a particle's structural and physico-chemical properties can contribute towards its DEP response. This response can also depend on the frequency of the applied electric field. Due to these dependencies, variations in applied field frequencies and external environment can simultaneously probe different particle substructures and processes. For example, some fundamental electrical properties of cells, such as membrane capacitance, membrane resistance and cytoplasmic conductance affect their DEP response. These properties also reflect a cell's ability to maintain ion balances and are a measure of metabolic work and biological organization. Thus, DEP can provide a non-invasive method for determining the electrical properties of cell populations, down to the single cell level.
  • DEP has potential uses in a number of fields.
  • DEP can be used as a drug discovery tool, e.g., monitoring the dielectrophoretic response of a cell population to candidate chemical compounds.
  • Other potential applications include separating particle populations using their differing dielectrophoretic response.
  • the invention features an apparatus and method for monitoring the dielectrophoretic response of one or more particles by capturing video sequences of the particles under the influence of varying DEP forces, and using image analysis to determine the response of each particle to the field.
  • the image analysis includes determining particle velocities as a function of electric field frequency for a known field profile, and determining the dielectrophoretic force on each particle from the particle size and velocity.
  • the invention features a method that includes (i) sequentially applying an electric field at a plurality of frequencies to a medium in a chamber, the medium comprising one or more particles (e.g., cells) suspended in a fluid; (ii) tracking the location of the one or more particles in the chamber while applying the AC electric field; (iii) calculating a velocity of each of the one or more tracked particles at each electric field frequency from the locations; and (iv) determining a characteristic of the one or more particles from the calculated velocities.
  • Embodiments of the method can include one or more of the following features and/or features of other aspects.
  • Tracking the location of the one or more particles can include acquiring a series of images of the one or more particles at each frequency. Tracking the location can include identifying the one or more particles in the series of images at a frequency. Calculating the velocity of a particle can include normalizing the velocity based on the size of the particle. In some embodiments, calculating the velocity of a particle includes normalizing the velocity based on the location of the particle. Alternatively, or additionally, calculating the velocity of a particle can include normalizing the velocity based on the electric field strength at the particle's location. The electric field strength is determined using a computer model (e.g., a computer model that accounts for a finite thickness of the electrode elements).
  • a computer model e.g., a computer model that accounts for a finite thickness of the electrode elements.
  • Calculating the velocity of a particle can include normalizing the velocity based on the electric field gradient at the particle's location.
  • the electric field gradient can also be determined using a computer model.
  • the parametric fitting function can include at least three fitting parameters.
  • Determining the characteristic of the one or more particles can include determining a frequency or range of frequencies where the one or more particles experience a maximum positive dielectrophoretic force. Calculating a velocity of each of the one or more particles can include calculating multiple velocity values of the velocity at the frequency or in the range of frequencies where the one or more particles experience a maximum positive dielectrophoretic force.
  • Determining the characteristic of the one or more particles can include determining a frequency or range of frequencies where the one or more particles experience a maximum negative dielectrophoretic force. Calculating a velocity of each of the one or more particles can include calculating multiple velocity values of the velocity at the frequency or in the range of frequencies where the one or more particles experience a maximum negative dielectrophoretic force.
  • the applied electric field can be generated by a waveform comprising sinusoidal waveform and/or by a digitally synthesized waveform.
  • the one or more particles can include particles that are labeled with a fluorescent moiety.
  • the method can further include treating the medium with a chemical or physical agent and determining an effect of the treatment on the characteristic of the one or more particles as a function time.
  • the method can include treating a surface adjacent the medium with an agent to modify the interaction of one or more of the particles with the surface.
  • the agent can cause the surface to selectively adhere or repel one or more of the particles.
  • the method can include determining an effect of the treatment on the characteristic of the one or more particles as a function of time.
  • the method can include identifying any of the one or more particles that form a pearl chain.
  • the characteristic can be determined from particles that do not form a pearl chain.
  • the invention features a method, including: (i) sequentially applying an AC electric field at a plurality of frequencies to a medium in a chamber, the medium including a plurality of particles suspended in a fluid; (ii) measuring a velocity of each of the plurality of particles at each electric field frequency; and (iii) determining a dielectrophoretic cross-over frequency for the particles from the measured velocities.
  • Embodiments of the invention can include one or more of the features of other aspects.
  • the invention features a method that includes: (i) sequentially applying an AC electric field at a plurality of frequencies to a medium in a chamber, the medium including a plurality of particles suspended in a fluid; (ii) acquiring a series of images of the plurality of particles at each frequency; (iii) measuring a velocity of each of the plurality of particles at each electric field frequency; and (iv) measuring a size of each of the plurality of particles.
  • Embodiments of the invention can include one or more of the following features and/or features of other aspects.
  • the method can include determining a cross-over frequency for the particles from the measured velocities.
  • Measuring the velocity of each of the plurality of particles can include monitoring positions of the plurality of particles in each of the series of images. Measuring the velocity of each particle can include determining the change in position of each particle in between images.
  • Determining a cross-over frequency for the plurality of particles can include determining an individual cross-over frequency for each of the plurality of particles and statistically analyzing the individual cross-over frequencies.
  • the invention features a method, including: (i) sequentially applying an AC electric field at a plurality of frequencies to a medium in a chamber, the medium including a plurality of particles suspended in a fluid; (ii) electronically acquiring a series of images of the plurality of particles at each frequency; and (iii) determining a cross-over frequency for the particles from the images.
  • Embodiments of the invention can include one or more of the features of other aspects.
  • the invention features a method, including: (i)sequentially applying an AC electric field at a plurality of frequencies to a medium in a chamber, the medium including a plurality of particles suspended in a fluid; (ii) acquiring a series of images of the plurality of particles at each frequency; and (iii)determining a velocity for each particle from each series of images.
  • Embodiments of the invention can include one or more of the features of other aspects.
  • the invention features a system for performing dielectrophoretic analysis on a population of particles.
  • the system includes a chamber having a wall, a plurality of electrodes disposed on the wall of the chamber, a function generator configured to supply an AC voltage to the plurality of electrodes, a detector configured to acquire images of a sample disposed in the chamber; and an electronic controller in electrical communication with the detector.
  • the function generator sequentially supplies an AC voltage to the plurality of electrodes at a plurality of frequencies while the detector acquires a series of images of a sample in a chamber, the sample including one or more particles (e.g., cells) suspended in a fluid, and the electronic controller tracks the location of the one or more particles in the chamber, calculates a velocity of each of the one or more particles at each electric field frequency from the tracked locations, and determines a characteristic of the one or more particles from the calculated velocities.
  • the function generator sequentially supplies an AC voltage to the plurality of electrodes at a plurality of frequencies while the detector acquires a series of images of a sample in a chamber, the sample including one or more particles (e.g., cells) suspended in a fluid
  • the electronic controller tracks the location of the one or more particles in the chamber, calculates a velocity of each of the one or more particles at each electric field frequency from the tracked locations, and determines a characteristic of the one or more particles from the calculated velocities.
  • Embodiments of the invention can include one or more of the following features and/or features of other aspects.
  • the system can include an optical microscope positioned relative to the chamber to image the sample to an image plane.
  • the detector can be positioned at the image plane.
  • the electrodes can include polynomial electrodes and/or interdigitated electrodes.
  • the plurality of electrodes includes an array of electrode regions.
  • the array of electrode regions can include an 8 x 12 array of electrode regions (e.g., configured for use with a 96-well plate).
  • Embodiments of the invention may include one or more of the following advantages.
  • Embodiments can enable rapid and/or accurate measurement of physical properties (e.g., DEP cross-over frequency and particle size) of a particle or population of particles. Due to speed and accuracy of measurements, embodiments can be used to investigate the DEP properties of large numbers of particles in a single cycle of the DEP apparatus (e.g., a sufficiently large number of particles to reduce statistical variance to useful levels).
  • the DEP cross-over frequency can be determined for a population of about 50 particles in one experimental run lasting less than two minutes.
  • the DEP characteristics of 1000 or more particles can be determined in one experimental run of less than two minutes.
  • particle parameters can be determined with greater statistical accuracy, and for biological particles errors associated with biological variability can be minimized.
  • FIG. 1 is a schematic diagram showing an embodiment of a dielectrophoresis system.
  • FIG. 2 is a plan view of an electrode apparatus used in the dielectrophoresis system shown in FIG. 1.
  • FIG. 3(a) and FIG. 3(b) are schematic diagrams showing motion of particles under (a) a negative dielectrophoretic force, and (b) a positive dielectrophoretic force.
  • FIG. 4 is a flowchart of an algorithm used to analyze images acquired using a DEP system.
  • FIG. 5 is a schematic diagram showing pixels used to determine a figure of merit for the radius of a particle.
  • FIG. 6 is a plot showing a particle's position as a function of time while under the influence of a negative dielectrophoretic force.
  • FIG. 7 is a plot showing a curve fitted to velocity/frequency data according to a three parameter fitting function.
  • FIG 8 is a schematic diagram showing interdigitated electrodes.
  • DEP system 100 includes an optical microscope 110, which images the particles onto a CCD camera 120.
  • the output of CCD camera 120 is connected to a PC computer 130.
  • the particles are confined to a DEP electrode apparatus 150, which is positioned on microscope stage 114 relative to the optical microscope's objective lens 112.
  • the output of a function generator 140 is connected to an electrode array in DEP electrode apparatus 150.
  • DEP electrode apparatus is also connected to PC computer 130 via a function generator 140 so that the PC can control the application of signals to electrode apparatus 150. This connection allows PC computer 130 to monitor particle activity at electrodes in DEP electrode apparatus.
  • function generator 140 applies an AC voltage to microelectrodes in DEP electrode apparatus 150, generating a spatially non-uniform electric field in the fluid.
  • the electric field when applied at certain frequencies, causes the particles in the fluid to move.
  • An image capture card installed in PC computer 130 acquires images from CCD camera 120. PC computer 130 stores and analyzes the acquired images.
  • the image capture card and CCD camera 120 can be selected to provide a suitable image acquisition rate.
  • a suitable acquisition rate is sufficiently high that a particle's displacement between successive frames under the influence of a dielectrophoretic force is sufficiently small that a particle can be tracked from one frame to the next.
  • a frame rate suitable for many applications is about 10 Hz, although the frame rate can be slower (e.g., about 1 Hz or slower) or higher (e.g., about 100 Hz or higher).
  • the higher frame rates could be used, for example, to monitor rapid cell physiology events, such as changes in membrane ion channel activity after exposure to a chemical agent (e.g., events that take place over milliseconds or tens of milliseconds).
  • DEP electrode apparatus includes an electrode driving board 210, which includes connectors for interfacing with the function generator.
  • Electrode driving board 210 supports an electrode loading tray 220 and an electrode array 230.
  • Electrode loading tray 220 facilitates easily loading of electrode array 230 into and out of DEP electrode apparatus 150.
  • Electrode array 230 includes several electrode regions 240 supported by a transparent substrate. Each of the electrode regions includes one or more electrodes. The electrodes are designed to provide a desired electric field profile to the fluid between the electrodes when a voltage is applied to the electrodes.
  • Electrode array 230 can be fabricated using photolithographic techniques. The electrodes are etched from a layer of a conductive material disposed on a transparent substrate (e.g., a glass substrate). The conductive material can be transparent (e.g., indium tin oxide) or opaque (e.g., gold or chromium). h the described embodiment, electrode array 230 is covered by a removable coverslip. Electrode array 230 and the coverslip form a chamber, which substantially confines the fluid and particles to a region of the electrode array 230 proximate to electrode regions 240. Typically, a user pipettes the fluid and particles onto the top surface of electrode array 230 and then covers the exposed suspension with the coverslip to form the chamber.
  • a transparent substrate e.g., a glass substrate
  • the conductive material can be transparent (e.g., indium tin oxide) or opaque (e.g., gold or chromium).
  • electrode array 230 is covered by a removable coverslip. Electrode array 230
  • an electrode region includes four electrode elements 310, 320, 330, and 340.
  • the electrode elements are separated by an electrically insulated region 350.
  • Center 360 of electrically insulated area 350 is defined as the point equidistant from electrode element tips 311, 321, 331, and 341.
  • Each electrode element pair located opposite insulated region 350 is connected to a voltage source.
  • electrode elements 310 and 330 are connected to one voltage source
  • electrode elements 320 and 340 are connected to another voltage source.
  • Particles 370 are suspended in the fluid adjacent insulated region 350.
  • Electrode elements 310, 320, 330, and 340 form polynomial electrodes.
  • polynomial electrodes refer to electrodes designed on the assumption that the electrical potential at any point created by an electrode system of interest is defined by a polynomial that obeys Laplace's equation. By substituting this polynomial into Laplace's equation the corresponding equipotentials can therefore be determined, and these in turn can be used to define the electrode boundaries.
  • a first AC potential difference is applied across electrode electrodes 310 and 340 (indicated as Visin(ffit) in FIG. 3(a) and 3(b)), and another AC potential difference is applied across electrode elements 320 and 330 (indicated as V ⁇ cos( ⁇ rt) in FIG.
  • the two potentials have same frequency and amplitude but are 90 degrees out of phase.
  • the potential difference between the electrodes generates an electric field in the fluid adjacent insulated region 350.
  • particles 370 can move in response to the electric field. Referring specifically to FIG. 3(a), for some frequencies, the particles can move away from electrode elements 310, 320, 330, and 340. In such cases, particles 370 motion is caused by negative dielectrophoresis. Referring now to FIG. 3(b), at other frequencies, particles 370 can move towards the electrode elements. This motion is caused by positive dielectrophoresis.
  • the amplitude of the applied voltage can be varied as desired. Typically, a voltage of sufficient amplitude to induce an observable dielectrophoretic response under investigation is applied to the electrodes at a given sample rate for the image capture card. For example, for electrode geometries where opposite electrode tips are separated by 0.5 millimeters, a voltage of about five volts peak-to-peak is sufficient for many particles.
  • the applied AC voltage is a sinusoidally-varying voltage, although other pulse shapes can be employed. Examples of other pulse shapes include square- wave pulse shapes and saw tooth pulse shapes.
  • the applied voltage can be the superposition of multiple sinusoidally-varying voltage waveforms. The frequency and/or amplitude of the sinusoidal waveforms can be the same or different.
  • One example of an applied voltage formed from a supe ⁇ osition of multiple sinusoidally-varying waveforms is a voltage formed from a waveform having frequency ⁇ and waveforms having harmonic frequencies (e.g., 2 ⁇ », 3 ⁇ >, 4 ⁇ , ).
  • system 100 can include a waveform synthesizer, and the applied voltage can be formed by digitally synthesizing a desired waveform.
  • the applied voltage waveform can vary as desired during data acquisition. Where the waveform is non-sinusoidal, the frequency refers to the number of times the waveform repeats itself per unit time.
  • the AC frequency may be varied as desired while images of the particles are acquired.
  • personal computer 130 will cause function generator 140 to sequentially scan through a range of frequencies while acquiring images from CCD camera 120 at a number of frequencies within the range. Several images (e.g., 100 or more) are acquired for each of the frequencies.
  • the range of frequencies scanned and the increment between successive frequencies depends on the sample being investigated.
  • the frequency range scanned can be from about 0.1 MHz to 1 MHz in 0.1 MHz increments.
  • the corresponding frequency range can be from about 10 kHz to 0.1 MHz, and for a population of bacteria a frequency range extending from around 1 MHz to 10 MHz can be used. More generally, frequency ranges can extend from the Hz and kHz range (e.g., about 1,000 Hz, 100 Hz, 10 Hz or less) to the MHz range or greater (e.g., about 0.1 MHz, 1 MHz, 100 MHz or more).
  • the increment between successive acquisition frequencies usually depends on the frequency range.
  • the increment will typically be small enough to provide a sufficient number of data points within the frequency range, but sufficiently large to scan the entire frequency range of interest. For example, where the range is small (e.g., 1,000 Hz) the increment will be small enough to provide sufficient data points in the range (e.g., 100 Hz increment providing 10 data points). Conversely, where the range is large (e.g., 1 MHz), the increment can be similarly large (e.g., 0.1 MHz).
  • the range of frequencies can be scanned linearly or non-linearly. During a linear scan, the increment between successive acquisition frequencies remains constant. During a non-linear scan, however, this increment can vary.
  • An example of a non-linear frequency scan is one in which large increments are used to determine a sub-range of frequencies in which the particles exhibit some identifying behavior. Within this sub-range, the frequency increment is reduced to provide additional data points, which can more accurately characterize the particles.
  • DEP system 100 can be used to study any dielectric particle that can be suspended in an appropriate fluid or medium.
  • the particles may be biological particles.
  • the particles can be cells, or components of cells and/or microorganisms. Components of cells include proteins and DNA. Microorganism's include bacteria. Biological particles also include pathogens, such as viruses.
  • Particles can be polymeric.
  • the particles may include polystyrene microspheres.
  • Solid particles include polymer spheres or protein macromolecules.
  • Semi-solid particles include poly- acrylamide or agar gel particles.
  • Liquid particles include the dispersed phase in an emulsion, such as oil droplets in water or liquid particles in an aerosol.
  • Gaseous particles include the dispersed phase in a foam, such as gas bubbles in a liquid.
  • particles can be tagged with a fluorescent moiety. Particles can be tagged for use with fluorescent microscopy techniques. Particle size may vary. Typically, particles are large enough to be observed using optical microscopy (e.g., larger than about 0.5 microns in diameter, such as 1 micron or larger). In some embodiments, particles can be larger than about 1 millimeter in diameter. In some embodiments, however, particles may be used that are too small to be observed directly using optical microscopy. For example, a sample may consist of protein or virus particles that have been fluorescently labeled to aid detection by optical microscopy.
  • a sample can include one type of particle, or a mixture of particles.
  • a sample may be a mixture such as a blood sample, including red and white blood cells, platelets, in addition to other types of particle.
  • An example of a single particle sample is a colloidal dispersion of polystyrene spheres in deionized water.
  • the number of particles in a sample may vary.
  • a sample may include a sufficiently large number of particles to take account of biological or natural variability of their size or properties.
  • a sample may include more than 10 particles (e.g., more than about 20 particles, 50 particles, 100 particles, 1000 particles).
  • PC controller 130 After PC controller 130 has acquired a series of images (also referred to as frames), it implements a software algorithm which identifies the particles in each frame, tracks the particles from frame to frame, determines the particles velocities at each applied frequency, and computes the particles' DEP cross-over frequency.
  • the DEP cross-over frequency is the frequency which, when applied to the electrodes, produces no net force on the particle under consideration.
  • the algorithm eliminates background pixels from each frame. Depending on the number of particles in the frame, many pixels in each image can belong to one of two sets, a low light intensity set corresponding to the electrodes, or a high light intensity set corresponding to the electrically insulated region between the electrodes. Because the particles usually correspond to intermediate intensity pixels, these low intensity pixels and high intensity images can be excluded from further analysis. Excluding these pixels can speed up the algorithm.
  • the algorithm first takes an intensity histogram of the entire image. Each pixel has 24 bit color depth (i.e., red, green, and blue sub-pixels each having eight bit gray scale). Accordingly, each sub-pixel has an intensity value from 0 to 255.
  • the intensity of each pixel is computed by averaging the red/green/blue values.
  • a histogram array with 256 bins is then used to create an intensity histogram of the entire image.
  • this histogram has two distinct peaks, which are identified using a windowing technique on the ID histogram array.
  • the low-intensity peak corresponds to pixels imaging the electrode
  • a high- intensity peak corresponds to pixels imaging the electrically insulated region. All pixels in the image that fall below a low intensity threshold (e.g., corresponding to the highest intensity pixel corresponding to an electrode) are assigned a first intensity value (e.g., 0).
  • Those pixels having intensity above a high-intensity threshold e.g., corresponding to the lowest intensity pixel corresponding to the electrically insulated region
  • a second intensity value e.g., 255
  • the next step in the analysis is to setup a coordinate system within which particles can be tracked from frame-to-frame.
  • a particle's position is defined by the distance from the center of the particle to a user specified origin (e.g., center 360 of electrically insulated area 350).
  • a user specified origin e.g., center 360 of electrically insulated area 350.
  • the velocity of a given particle is expected to be proportional to the time-dependent distance between that particle and center 360 of electrically insulated area 350.
  • the algorithm determines the pixel coordinate of center 360 of electrically insulated area 350.
  • the algorithm determines pixel coordinates of tips 311, 321, 331, and 341 of electrode elements 310, 320, 330, and 340, respectively.
  • the algorithm starts out at the center vertical pixel on the left side of the image and looks for transition from the second intensity value to the first intensity value (e.g., from 0 to 255).
  • the last pixel having the second intensity value is defined as the left side of the circle.
  • a similar procedure is used to find tips 311, 321, and 331.
  • the center of the coordinate system is defined as the average of the location of the left and right tips (i.e., tip 341 and tip 321, respectively) for the horizontal coordinate and the average of the top and bottom tips (i.e., tip 311 and tip 321, respectively) for the vertical coordinate.
  • the algorithm After setting up a coordinate system, the algorithm proceeds to scan each frame for particles.
  • the algorithm can restrict the search to an area inside of the electrodes, and can ignore the background in this region. Accordingly, the algorithm identifies individual particles within an area of interest and provides a rough estimate of the size and position of these particles.
  • the algorithm finds particles by scanning the frame column-by-column and identifying pixels that have an intensity value different from the background intensity (e.g., background pixels that have previously been assigned the first intensity value).
  • the algorithm finds a pixel that has an intensity value different from the background intensity value
  • the algorithm looks for adjacent pixels that also have intensities different from the background. Accordingly, the algorithm identifies sets of adjacent pixels. If a set matches certain criterion (e.g., number of pixels, vertical and horizontal range of the set), then this set of pixels is determined to be a particle.
  • the size and position of the cell may also be roughly determined by the size of and position of the set of adjacent pixels.
  • the algorithm determines the size and position of each identified particle. In this step, the algorithm uses the previously determined rough estimate of the center and size of each particle to more accurately determine the center and size of each particle.
  • the algorithm fits a predetermined shape (e.g., a circle, an ellipse, or a polygon, such as a square or hexagon) to each particle.
  • a predetermined shape e.g., a circle, an ellipse, or a polygon, such as a square or hexagon
  • the algorithm can fit a circle to the dark part (i.e., low intensity) of an image of each cell's membrane.
  • the image of the cell includes three concentric circular regions, a dark circle sandwiched between two lighter circles.
  • the algorithm identifies the intensity of pixels corresponding to a set of pixels at a series of different azimuthal angles, ⁇ , measured from a horizontal axis originating at the estimated center of the cell.
  • the algorithm determines an overall figure of merit for a cell image, according to the intensity of the nearest neighbor pixels to the pixel located the estimated cell radius away from the estimated cell center.
  • the overall figure of merit is derived from the addition of two parameters.
  • the first parameter quantifies the darkness of the cell membrane at the edge of the fitted shape. The larger this first parameter is, the more distinct the cell image. This first parameter rapidly reduces in magnitude for cells that become out of focus or have a physically disrupted membrane (e.g., after bursting).
  • the second parameter quantifies the uniformity of the membrane darkness around the cell perimeter. The smaller this second parameter, the more non-uniform is the cell membrane. This second parameter is sensitive to very uneven, discontinuous or oddly shaped cell membranes.
  • the center coordinate and cell radius are modified slightly (e.g., by moving a pixel in the horizontal and/or vertical direction) and the two figure of merit parameters are computed again.
  • the algorithm determines a figure of merit sum for the cell, hi an iterative procedure, the algorithm modifies one or more of the estimated vertical position of the cell; the estimated horizontal position of the cell; and/or the estimated radius of the cell. Changes in one of these parameters that increase the figure of merit sum for all positions are rewarded by continuing to change that parameter in the direction that increased fitness. When further changes to any of the three parameters fail to yield increases in fitness, then the iteration cycle is complete, and the final estimate of the cell center position and radius are taken as the cell center position and cell radius for that frame.
  • the figure of merit and predetermined shape used to determine the figure of merit are selected according to the type of particle being studied. For example, while the above description refers specifically to a figure of merit and shape (i.e., circle) for when the particles being studied are cells, other shapes and figures of merit may be used.
  • the resulting data includes the position (e.g., an x, y coordinate) and size (e.g., radius) of each particle with respect to the coordinate system for each frame. It may be desirable to run the aforementioned steps more than once for the first frame (or first few frames) in an experiment in order to obtain initial estimates of the position of each cell.
  • the algorithm calculates each particle's velocity for a series of frames.
  • a series of frames corresponds to frames acquired for a particular frequency. Accordingly, by performing this analysis for each acquisition frequency, the analysis provides a particle velocity as a function of applied field frequency.
  • a particle's velocity is determined from the displacement of the particle's center position between different frames, and the time increment between successive frame acquisitions.
  • the velocity is computed by creating a particle tracking file that appends the tracked particle position and size for each particle in each frame using the previously described algorithm and then comparing the change in position corresponding to the frame rate.
  • the velocity can be expressed in terms of pixels per second, or, where the apparatus is calibrated and the imaged area known, velocity can be expressed as, e.g., microns per second.
  • a particle's dynamic response can be determined by studying its position as a function of time.
  • the particle's position is determined as its radial distance from center 360.
  • the particle experiences negative DEP force and moves toward center 360.
  • the particle accelerates in the first second under the influence of the dielectrophoretic force and reaches a steady-state velocity after about one second.
  • the particle's velocity is determined by performing a linear regression on this data, ignoring those data points acquired in the first second while the particle is accelerating.
  • a velocity is determined for one or more particles at each applied frequency.
  • the frames used to compute each particle's velocity at a frequency can be a subset of the frames acquired at that frequency. For example, frames acquired while the particle is not at a steady-state velocity can be ignored while determining particle velocity.
  • the algorithm also adjusts the velocity to compensate for the electrode geometry. For example, for the polynomial electrode shape described above, the velocity is adjusted by multiplying by a factor of (1/7?), where R is the radial position of the particle center in each frame. This compensation normalizes the particle velocities to account for variations in the electric field strength and/or gradient, which are both dependant on the particles position with respect to the electrode.
  • the velocity of each particle can be normalized with respect to variations of the electric field strength and electric field gradient between the electrode edges.
  • the electric field profile for a particular electrode geometry can be determined by computer modeling.
  • An example of computer software that can be used to model the electric field profile for different electrode geometries is Maxwell 3D Electromagnetic Field Simulator from Ansoft (Pittsburgh, PA).
  • the electric field profile is modeled accounting for the three-dimensional nature of the electrodes. In other words, the model accounts for the finite thickness of the electrodes.
  • a particle Under the influence of a dielectrophoretic force, a particle reaches its steady- state velocity when the dielectrophoretic force accelerating the particle is balanced by the viscous force of the fluid resisting the acceleration. Because both viscous force and the dielectrophoretic force depend on the particle size, particle velocity also depends on particle size. Accordingly, in some embodiments, the particle velocities can also be normalized with respect to each particle's size.
  • the algorithm determines the DEP cross-over frequency for each particle.
  • the DEP cross-over frequency is computed using velocity values for a given particle over a range of frequencies.
  • there should be at least one measurement (preferably several) at a frequency that provides a negative dielectrophoretic response for the particles well below the DEP cross-over frequency e.g., for many cells suspended in an aqueous medium of conductivity of around 40 mS/m, a 20 kHz frequency is sufficiently low
  • at least one measurement preferably several at a frequency that provides a positive dielectrophoretic response well above the DEP cross-over frequency (e.g., 500 kHz for most cells in a suspending medium conductivity around 40 mS/m).
  • pearl chains In some cases, some particles attract each other to form particle chains. Such particle chains are often referred to as “pearl chains” and the phenomenon is referred to as “pearl chaining", which is described in detail in Electromechanics of Particles , by T. B. Jones, Ch. 6 & 7, Cambridge University Press (1995).
  • the response of particles in pearl chains to an applied electric field can be different from free particles (i.e., particles that have not formed pearl chains). Accordingly, it can be advantageous to eliminate data corresponding to pearl chains from subsequent analysis.
  • System 100 can analyze how the DEP cross-over frequency of a particle is altered as a result of it interacting with other particles to form pearl chains. This is achieved by monitoring the separation distance between tracked particles. The user can define a minimum separation distance between adjacent tracked cells. When two particles are closer than this minimum separation distance, the particles are identified as a pearl chain composed of two particles, and the reported DEP cross-over frequency is reported separately from the data for single particles to allow for separate analysis. Multiples of particles attached together in pearl chains are similarly reported, with an indication of the number of particles tracked as part of that chain.
  • the time averaged DEP force ⁇ F(t)> acting on a cell can be described using the Maxwell stress tensor formulation [see, e.g., Sauer, FA., in: Chiabrera, A., Nicolini, C, Schwan, H.P. (Eds.) Interactions between Electromagnetic Fields and Cells, Plenum Publ. Co ⁇ . 1985, 181-202.], but for the common case, where dielectric losses are not too high, a more computationally convenient equivalent dipole moment formulation can be used.
  • the force acting on a particle in a field of magnitude E and gradient VE is given by:
  • Re denotes the real (in phase) component of the dipole moment, m( ⁇ ), induced in the cell.
  • the factors ⁇ and ⁇ m are the particle and suspending medium complex
  • ⁇ is the real component of the Clausius-Mossotti factor defined by equation 3, and which theoretically can have a value ranging from -0.5 to +1.0.
  • the algorthim then proceeds to fit a curve to the velocity vs. frequency data using a three parameter model derived as follows.
  • the real component of the Clausius Mossotti can be expressed in terms of the conductivities, ⁇ p and ⁇ , occasion, and permittivities, ⁇ p and ⁇ m , of the particle and medium, respectively, as:
  • ⁇ represents the characteristic Maxwell- Wagner relaxation time describing the polarization (charge accumulation) that occurs at the interface between the particle and its surrounding suspending medium for low volume fractions of suspended particles.
  • equation 11 reduces to approximately ⁇ !.
  • the DEP cross-over frequency is given by the fitting parameter equation 12 and the fitting parameter B.
  • DEP system 100 can be used to determine the DEP response of particles under varying chemical and/or physical conditions. Knowledge of a particle's DEP response can be used for variety of pu ⁇ oses, including, for example, characterizing the properties of specific particle types, analyzing sub-populations of particles in particle mixtures, and/or analyzing the behavior in response to exposure to different chemical and/or physical agents.
  • DEP system 100 can be used to study the response of a population of cells to a chemical agent.
  • DEP system 100 can be used as a tool to follow physiological changes in a cell membrane that accompany transmembrane signaling events (e.g., to study human T lymphocytes stimulated using phorbol myristate acetate and ionomycin).
  • Changes in the DEP response e.g., changes in the DEP cross-over frequency
  • changes in the DEP cross-over frequency of such cells after activation can be correlated to changes in the cells' membrane topography (e.g., reduction of membrane associated microvilli, blebs, and/or folding) via a theoretical model.
  • DEP system 100 can be used to identify agents that produce a specific response in a cell population for the pu ⁇ oses of drug discovery.
  • DEP system 100 can be used to determine parameters that are appropriate for the separation and/or identification of a particular type of particle from other particles, or to differentiate between different types of particles present in the dispersion.
  • DEP system 100 can be used to analyze components of a mixture of several particles.
  • DEP system 100 can also be used to determine particle adhesion to a substrate. Because the dielectrophoretic force on a particle can be determined by correlating the particle's position to the electric field profile, the amount of force required to dislodge the particle from a substrate can be determined when the particle starts to move. By modifying a substrate (e.g., by chemical and/or physical surface treatment), DEP system 100 can be used to monitor changes in a particle's interaction with the substrate as a result of the modification.
  • injecting and/or releasing the fluid into and out of the chamber can be performed manually or can be automated.
  • fluid is manually supplied to the chamber by the user.
  • fluid supply can be automated.
  • the chamber may include an inlet port and an exhaust port, for supplying sample fluid into and removing sample fluid from the chamber.
  • Tubing is connected between a fluid supply to the chamber.
  • a pump draws a predetermined volume of fluid from the fluid supply and supplies it to the chamber on demand.
  • the fluid sample is pumped out of the chamber through the exhaust port.
  • the chamber can be flushed with a cleansing sample (e.g., deionized water) before the next fluid sample is introduced.
  • personal computer 130 can control the introduction and removal of fluid samples from the chamber.
  • PC 130 can also be used to control the function generator. This can increase the level of automation in the system, because, once initiated, the system can scan the entire frequency range without additional user input.
  • the above-described embodiment includes four polynomial electrode elements, other numbers of electrode elements and/or other electrode geometries can also be used.
  • the sum of the phase difference between the voltage applied to adjacent electrodes should be 2 ⁇ .
  • the phase difference between each adjacent element is 7z 2 radians.
  • the phase difference between the voltage applied to adjacent electrodes should be 2 ⁇ /N radians.
  • the voltage can be applied to the electrodes to produce a traveling wave electric field.
  • interdigitated electrodes 800 includes 10 parallel electrode elements. Electrode elements 811, 812, 813, 814, and 815 are in electrical contact to a first bus line 810, while electrode elements 821, 822, 823, 824, and 825 are in electrical contact to bus line 820. During operation, a power source applies an AC voltage between bus line 810 and bus line 820. Accordingly, the potential difference between electrode elements 811-815 and 821-825 gives rise to an electric field between the electrodes.
  • the dimensions of the interdigitated electrode elements can be varied as desired.
  • the electrodes are about five or more microns wide (e.g., 10 microns, 20 microns, 50 microns, 100 microns) and about 20 or more microns long (e.g., 30 microns, 40 microns, 50 microns, 75 microns, 100 microns, 200 microns, 500 microns).
  • the separation between adjacent electrode elements can vary. In some embodiments, the separation between adjacent electrodes is greater than the electrode's width. The separation can be, for example, more than about 10 microns (e.g., 20 microns, 30 microns, 50 microns, 100 microns, 200 microns, or more).
  • interdigitated electrodes 800 includes 10 electrode elements, other embodiments can have fewer or more electrode elements (e.g., more than 20 elements, 50 elements, 100 elements).
  • electrode array 230 includes multiple electrode regions, each electrode region being configured to apply an electrode field to a different sample of particles, h particular, electrode array 230 can be configured for use with an array of samples, such as for use with a 96-well plate. In the case of a 96-well plate, electrode array 230 includes an 8 x 12 array of electrode regions, each region corresponding to a sample in the 96-well plate. Such embodiments can provide high sample throughput as the system can be configured to automatically study each of the samples in sequence, with minimal additional user input after initial configuration. In general, electrode array 230 can include any number of electrode regions (e.g., less than 96, or more than 96, such as 384).
  • DEP system 100 can also include devices with which to apply additional forces to the particles.
  • DEP system 100 can include devices for applying hydrodynamic force, centrifugal force, gravitational force, ultrasonic force, electrophoretic force, magnetic force, and/or an optical force to one or more of the particles in the dispersion.
  • an output beam from a laser can be focused through the microscope objective lens to provide an optical trap in the chamber.
  • a hydrodynamic force can be applied by flowing a liquid through the chamber at a predetermined rate.
  • An example of a device for applying an ultrasonic force is an ultrasonic transducer, which can be coupled to the chamber. To apply a centrifugal force, one could rotate the chamber at a rotational velocity sufficient to apply the desired force.
  • DEP system 100 can include additional probes to monitor various parameters in the chamber that can affect the dielectrophoretic response.
  • DEP system can include a conductivity sensor to monitor the conductivity of the sample in the chamber during dielectrophoretic study, hi some embodiments, DEP system 100 can include a thermometer to monitor the temperature of the sample.
  • Other parameters that might affect the dielectrophoretic response that can be monitored include dielectric permittivity, pH value, fluid viscosity, and/or the material or chemical treatment of the electrode and/or chamber walls.
  • Personal computer 130 can monitor additional parameters of the system that can affect the dielectrophoretic response of the particles. For example, by including a temperature probe on electrode array and connecting the probe to it, the PC can be used to monitor the temperature of the dispersion during DEP manipulation. PC 130 can monitor other parameters, such as those mentioned above.
  • images are detected using CCD camera 120, however, other types of detectors can also be used (e.g., CMOS imaging chips or optical waveguides).
  • the PC computer and image capture card can be replaced by a custom electronic processor designed to analyze signals from the detector using signal processing techniques.
  • optical microscope 110 is an upright transmission mode configured for bright field imaging.
  • optical microscopy can also be used, including, for example, dark field microscopy, fluorescence microscopy, and/or confocal microscopy. More generally, the above- described techniques can be applied in conjunction with and/or using non-microscopic modes of data capture, for example, CCD arrays, CMOS arrays, or using a plate reader.

Landscapes

  • Investigating Or Analysing Biological Materials (AREA)
  • Investigating Or Analyzing Materials By The Use Of Electric Means (AREA)

Abstract

La présente invention concerne un système utilisé pour effectuer l'analyse diélectrophorétique d'une population de particules, le système comprenant: une chambre (150) comportant un puits; une pluralité d'électrodes situées sur la paroi de la chambre; un générateur de fonctions (140) configuré pour appliquer à la pluralité d'électrodes, un courant alternatif (CA) à une pluralité de fréquences; un détecteur (120) configuré pour acquérir des images d'un échantillon placé dans la chambre; et un dispositif de commande électronique (130) électriquement relié au détecteur (120).
EP02794251A 2002-12-12 2002-12-12 Systeme et procede de definition de profils dielectrophoretiques de particules Withdrawn EP1573308A1 (fr)

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
PCT/US2002/039931 WO2004055505A1 (fr) 2002-12-12 2002-12-12 Systeme et procede de definition de profils dielectrophoretiques de particules

Publications (1)

Publication Number Publication Date
EP1573308A1 true EP1573308A1 (fr) 2005-09-14

Family

ID=34374705

Family Applications (1)

Application Number Title Priority Date Filing Date
EP02794251A Withdrawn EP1573308A1 (fr) 2002-12-12 2002-12-12 Systeme et procede de definition de profils dielectrophoretiques de particules

Country Status (3)

Country Link
EP (1) EP1573308A1 (fr)
JP (1) JP2006510020A (fr)
AU (1) AU2002359698A1 (fr)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP5849503B2 (ja) * 2011-08-01 2016-01-27 株式会社Ihi 粒子観察装置及びマニピュレート装置
JP5991236B2 (ja) * 2013-03-13 2016-09-14 ソニー株式会社 分取装置
KR102244778B1 (ko) * 2018-08-31 2021-04-27 연세대학교 원주산학협력단 유전영동힘에 반응하는 세포의 위치변화를 이용한 세포의 교차주파수 측정방법
KR102221953B1 (ko) * 2018-08-31 2021-03-04 연세대학교 원주산학협력단 유전영동 힘 인가변화 조절에 따른 세포의 전기분극 특성 분석법

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
See references of WO2004055505A1 *

Also Published As

Publication number Publication date
AU2002359698A1 (en) 2004-07-09
JP2006510020A (ja) 2006-03-23

Similar Documents

Publication Publication Date Title
US7063777B2 (en) Dielectrophoretic particle profiling system and method
US6264815B1 (en) Apparatus and method for testing using dielectrophoresis
FI94646C (fi) Mikro-organismien ja muiden hiukkasten dielektroforeettinen määritys
EP3410107B1 (fr) Dispositif d'analyse et dispositif de séparation
Vahey et al. High-throughput cell and particle characterization using isodielectric separation
US8262886B2 (en) Apparatus for analyzing characteristics of particulate with dielectrophoresis of particulate by applying angle-modulated wave and method for the same
US20110290649A1 (en) Device and method for measuring microspheres
Antunez-Vela et al. Simultaneous determination of linear and nonlinear electrophoretic mobilities of cells and microparticles
WO1994029694A1 (fr) Mesure de la taille et de la charge de particules dans des colloïdes a multicomposants
Zhang et al. Characterization of single-cell biophysical properties and cell type classification using dielectrophoresis model reduction method
Reichle et al. Electro-rotation in octopole micro cages
JP2007006858A (ja) 微生物検査チップおよび微生物検査方法
KR20060085299A (ko) 농도구배발생부가 구비된 유전영동 장치, 그를 이용한물질의 분리방법 및 물질 분리의 최적 조건을 탐색하는 방법
US20180333724A1 (en) Assay systems and methods for processing sample entities
Rempfer et al. Selective trapping of DNA using glass microcapillaries
Ziebert et al. A planar lipid bilayer in an electric field: membrane instability, flow field, and electrical impedance
EP1573308A1 (fr) Systeme et procede de definition de profils dielectrophoretiques de particules
Abd Samad et al. Dielectrophoresis velocities response on tapered electrode profile: Simulation and experimental
US20180292352A1 (en) Methods and systems for identifying a particle using dielectrophoresis
JP3734131B2 (ja) 微生物数測定装置および微生物数測定方法
JP2005224171A (ja) 誘電泳動活性解析方法とそのための装置
EP4389282A1 (fr) Capture de particules d'analyte dans un liquide
KR20200026171A (ko) 유전영동힘에 반응하는 세포의 궤적분석을 통한 세포의 cross-over frequency 측정방법
Sherif et al. An Influence of the Microfluidic Channel Height and Distribution of Dielectrophoretic Force on the Impedance Extraction in Microfluidic Systems
JP2009192479A (ja) 微粒子測定装置およびそれに用いる電極

Legal Events

Date Code Title Description
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

17P Request for examination filed

Effective date: 20050627

AK Designated contracting states

Kind code of ref document: A1

Designated state(s): AT BE BG CH CY CZ DE DK EE ES FI FR GB GR IE IT LI LU MC NL PT SE SI SK TR

AX Request for extension of the european patent

Extension state: AL LT LV MK RO

DAX Request for extension of the european patent (deleted)
STAA Information on the status of an ep patent application or granted ep patent

Free format text: STATUS: THE APPLICATION IS DEEMED TO BE WITHDRAWN

18D Application deemed to be withdrawn

Effective date: 20060704