METHOD AND APPARATUS FOR DETERMINING LOCAL NISCOELASTTCTTV
Background of the Invention
In the study of colloidal suspensions, polymer solutions, and in cell and microbiology it is often desirable to have a direct observation in real time of the changes in the local viscoelasticity which affect the dynamics of objects such as colloidal particles, cells and organelles. For example, the cells transport and diffusion of macro-molecules both depend on the local environment. It is especially important to have high spatial resolution in complex environments in which local inhomogeneities are present.
Existing methods such as localized dynamic light scattering (Bar-Ziv et al., Phys. Rev. Lett. 78: 154-157 (1997)) require a complex setup. In most methods, determining the stiffness of a trapping beam requires a well-calibrated position detector (Nisscher et al., Selected topics in quantum electronics 2: 1066-1076 (1996)). Determining the characteristic time of the motion is done through a cumbersome fit calculation to the power spectrum of the position, which requires substantial time to collect. Digital video microscopy (Crocker et al., Phys. Rev. Lett. 82: 4352-4355 (1998)) is a good tool to study the dynamics but requires the processing of many digitized video images, which can not be done in real time and is often complicated by the image quality. A method which can be used to measure local changes in viscoelasticity at high spatial resolution and high temporal resolution would represent a significant advance in the art.
Summary of the Invention
The present invention relates, in one aspect, to an apparatus and method for determining local viscoelasticity in a fluid or a gel. The method employs a scanning optical microscope in combination with a force manipulating function. A particle is provided within the fluid or gel in which the local viscoelasticity determination is to be made. The particle is periodically modulated using the force manipulating function to set the particle in periodic motion about the focal region of the scanning optical microscope. The second harmonic of the signal from the focal region of the scanning optical microscope is determined with respect to the particle in motion. The phase delay of the second harmonic determined in this manner is used to calculate local viscosity.
Brief Description of the Drawings
Fig. 1 is a diagrammatic representation of a confocal embodiment of the present invention.
Fig. 2 is a diagrammatic representation of a confocal embodiment of the type employed in the experiments described in the Exemplification section set forth below.
Fig. 3 is a diagrammatic representation of measurements made at a constant laser power of 1.3 mw.
Detailed Description of the Invention
The present invention relates to a method for determining local viscoelasticity in a fluid or a gel. This invention finds application in a variety of fields including, for example, biological and chemical applications.
With respect to biological application, in preferred embodiments, the fluid is a biological fluid (i.e., a fluid produced by an organism), or a non-biological fluid which contains biological cells or biomolecules (i.e., molecules produced by an organism). A biological fluid may be an intracellular fluid, such as cytosol. Alternatively, a biological fluid may be an extracellular fluid such as serum. Examples of non-biological fluids which contain biological cells or biomolecules include cell culture medium, buffered saline solution, tissue culture medium, etc.
The apparatus used to determine local viscoelasticity in the method of the present invention is comprised of multiple elements. A scanning optical microscope provides the ability to probe the local environment in a fast and sensitive manner. Examples of scanning optical microscope formats useful in connection with the present invention include confocal microscopes as well as two photon microscopes.
A force manipulating function is also provided in connection with the methods of the present invention. The force manipulating function is used to sequester and manipulate a particle within the fluid or gel of interest. Examples of suitable force manipulating functions include, for example, photon force manipulating functions (e.g., laser trap or optical tweezers) and magnetic force manipulating functions. With respect to the former, the laser is preferably tuned to the infrared spectrum.
To determine local viscoelasticity, it is necessary to provide a particle which is suitable for manipulation in the fluid or gel using the force manipulating function. The required properties of the particle for use in connection with the methods of the invention vary depending upon the nature of the force manipulating function. For example, when using a laser trap as a force manipulating function, the particle must be transparent and have a refractive index higher than that of the surrounding medium. When using a magnetic force manipulating function, the particle must be a magnetic particle. It should be noted that suitable particles for use in connection with the present invention may be naturally ocurring in the fluid or gel of interest. For example, it is a matter or routine experimentation to determine whether an intracellular organelle, such as mitochondria or chloroplasts, can be visualized, trapped and set in motion using the force generating function of the present invention.
Once sequestered using the force generating function of the present invention, the particle is set into periodic motion by modulating the force generating function with an oscillation amplitude that is slightly larger than the diameter of the particle of interest. The resulting particle motion is periodic with a frequency equal to the driving frequency and has a phase lag because of the drag.
The second harmonic of the signal from the focal region of the scanning optical microscope is determined with respect to the particle in motion. The phase delay of this second harmonic is then determined with respect to the second harmonic of the period modulation of the force function. Such determinations are made, for example, using a phase sensitive detector. These values are then used to determine local viscoelasticity. Given the present disclosure, one of skill in the art could derive and solve local viscoelasticity equations without the use of undue experimentation. Many sets of equations could be derived to determine local viscoelasticity depending, for example, upon the degree of precision required. It is not possible within the context of a patent application to set forth all of the many possible variations which ultimately will yield a value for local viscoelasticity.
Following is an illustration of a solution route for a variant of a photonic force microscope in which confocal detection is employed. This simplifies alignment, improves stability and allows scanning imaging. In confocal detection we receive a signal from the
bead whenever it is in the focal region of the trapping beam. Therefore, in typical applications in which the bead is moving about in the trap, it is necessary to have a model to indicate when the position x of the particle is the same as the position/? of the trap. In the following, we develop a mathematical model for u =x —p and expect a signal whenever u = 0.
In optical tweezers, particles are trapped by the force exerted by the electric field of the trapping beam on the dipole moment induced in the trapped particle. The dipole moment μ induced by the optical field is given by μ = aE where a is the polarizability of the trapped particle. For particles much smaller than the wavelength of light is given by
a = 3ε g2 ~ g' V ε2 + 2εl
Fis the particle volume, ε the particle dielectric constant and ε \ the dielectric constant of the host medium. In one dimension the force on the particle is
so that the particle finds itself in a potential
The motion of a micron-sized particle in common fluids takes place at small
Reynolds number and viscous drag dominates inertial forces. The corresponding equation of motion of the particle in the trap, ignoring Brownian motion for the moment is dx _ + E = 0 dt where ; is the hydrodynamic drag coefficient. For a sphere far from any surface γ equals
6πηa where a is the radius of the sphere and η is the dynamic viscosity. Suppose the optical electric field of the trap is Gaussian, as occurs in a diffraction-limited beam. Then the field amplitude is proportional to exp(-^2 / 2σ2 ) so that the equation of motion is dx γ — + a(x-p(t))&φ { (*-p(t))2 = 0. dt 2σ2
If the tweezer is set into sinusoidal motion in a fluid flowing with velocity v then
p(t) = vt + a sin ωt . The corresponding equation for u is
These differential equations can be readily solved using the Runge-Kutta method, or can be approximated by linear differential equations with constant coefficients if the experimental conditions are such that the particle stays near the center of the trap and a parabolic approximation can be made for the particle potential. It is convenient to recast the equation of motion into dimensionless form = —V - AcosT
where U= u/σ is the relative bead position in units of σ, T=ω t is the normalized time, and A-alσ is the oscillation amplitude of the trap in units of beam width σ and N is a normalized velocity. In the linear case this equation reduces to
— + CU = -V - A cosT dT whose solution is
V A(CcosT + sinT) C l + C2
Zero crossings occur at
For the simple case V=Q, the second harmonic phase is then
φ2 = 2tan_1 C
Local viscoelasticity determinations of the type described above find application in a variety of biological contexts. When cells produce biological materials such as proteins the viscosity of the fluid around and within the cells changes. The greater the change, the greater the concentration of product. Disclosed is a novel scanning probe microscope capable of forming spatial images of viscosity distributions to allow, for example, selection of overproducers.
In a preferred embodiment, the instrument combines optical tweezers with confocal optical detection. In an optical tweezer the trapped particle appears as an object in the focal point of the laser beam. If the particle moves about in the trap, the signal detected confocally from it will vary. If the tweezers are set into slight sinusoidal spatial oscillation then the position of the particle in the trap can be accurately monitored by lock-in detection of the confocal signal. In regions of higher viscosity, the trapped particle will lag further in phase behind the oscillating beam. Such an instrument can also be used to map fluid flow distributions. In regions of faster flow, the particle will be further pushed to one side of the oscillating beam. This offset too can be monitored by confocal detection. This method of microscopy will address the need in biology and medicine for real-time visualization of fluid parameters in cellular systems such as cultures. Typical applications are:
• Selection of overproduces of bioproducts such as polysaccharides (eg Pullulan). Cells secrete these polysaccharides in a halo around their bodies. A real time map of the spatial extent of the viscous polysaccharide distribution would enable realtime screening.
• Selection of hybridomas from a background population of precursor cancer cells. Hybridomas are efficient antibody factories.
• Selection of secretory mutants for protein expression. One problem with protein expression in bacteria is that the product is often left inaccessible inside the bacteria, or is hindered by a toxic buildup of the expression product inside the cell.
• Analysis of mass transfer of gasses in cell populations.
• Cell sorting by refractive index signature. The trapping strength of the optical tweezers in the microscope increases with the refractive index of the trapped particles. Sporulating bacteria can be used to make antibiotics. The microscope could select these by the high refractive index of their spore bodies.
• Differentiation of live cells from dead cells, by examining differences in their extracellular environment.
Cell interiors can be explored by using the microscope to manipulate intracellular particles such as mitochondria and chloroplasts (especially important for selection of
plant cell lines). Refinements of the present invention include variants involving fluorescent microscopy and dielectric and magnetic particles. The present invention relates also to the three dimensional case with traps of Gaussian rather than quadratic profile, and include corrections due to the fact that the beads in use actually have radii too large to use in the Rayleigh (small bead) limit.
Even a continuous wave laser can elicit two-photon fluorescence in appropriately labeled beads. The fluorescence intensity is highest when the particle is at the center of the trap, so the method of second harmonic lock-in detection is also useful in this context. One significant advantage of two-photon fluorescence is that visible fluorescence can be excited with minimally damaging near infrared light. If regular single photon excitation were employed, the tweezer wavelength would have to be in the visible, where biological materials show considerable absorption and photodamage. It is true that a weak visible exciter beam could be made collinear with an infrared trapping beam but this will come at substantial alignment cost and increased system complexity. Another significant advantage is that two-photon fluorescence is proportional to the square of the intensity. This localizes the fluorescence to the focal volume of the trap, and eliminates the need for a physical pinhole in the confocal detection unit.
Magnetic beads may also be used in place of dielectric beads. In this context, a time varying magnetic field is used to force the beads into oscillation. A substantially lower laser power can be used since the laser no longer need to supply the trapping force. A feedback system would be designed to lock the probe laser to the bead under test by nulling out the fundamental signal in x, y and z, or by maximizing the average confocal signal. The second harmonic signal would give information on viscosity as described above.
Sample Stage
Detector
Figure 1
In one embodiment, the present invention relates to a device of the type shown in Figure 1. In this embodiment, a scanning mirror oscillates the trap back and forth and a lock-in amplifier measures the phase shift between the confocal signal and the second harmonic of the spatial oscillation frequency of the trap. At higher viscosities, this phase shift increases and as shown below a mathematical model can be developed to enable determination of the medium viscosity.
An important distinction between the present invention and prior art work is the associated use of 2nd harmonic phase detection (using for example a lock-in amplifier) to monitor the position of the probe particle relative to the center of the oscillating trap. In the usual use of lock-in detection the signal is detected at the same frequency (fundamental) as the sinusoidal oscillation of the stimulus. But in connection with the present invention, the bead comes into focus twice every cycle, so there will be no signal at the fundamental. Detection at twice the frequency of oscillation of the stimulus is required and, therefore, second harmonic detection is required.
Exemplifi cation
1. Theoretical Considerations
The equation of motion of a particle in a harmonic potential is: βx(t) + a - [x(t) - p(t)] = 0 (1) where x(t) is the particle position, p(t) is the position of the center of the potential well (the focused laser beam in our case), β and a are the drag coefficient and the trap stiffness respectively. We have neglected the inertial term and the Langevin (stochastic) force. For a potential moving in a sinusoidal fashion, this equation becomes:
x(t) + ^& = ?° smωt (2)
T T where po and ω are the amplitude and angular frequency of the laser motion,, and τ is the characteristic time scale of the motion: τ = - (3
Since the trap stiffness is proportional to the optical power P. we get
where C Q is a measure of the trapping efficiency, b is an invariant property of the combined particle-medium-optics system (independent of the laser power or the driving frequency) and has units of energy. The solution to equation 2 is:
XV) = (1 + ^2)1 • SinM + θ) (δ)
where the phase difference between the particle's position relative to the laser beam
θ — - tan 1(ωr) (6)
Given that relation 6 holds, a measurement of θ at a given ω constitutes a measurement of r.
2. Experimental Design
The spectrum of the back-scattered light signal from a laser beam moving with a frequency ω across a microparticle will exhibit peaks at the second harmonic since the beam traverses the particle twice in one cycle and since the confocality condition means that we detect a signal only when the focused laser beam coincides with the position of the particle.
The particle also exhibits Brownian motion, which appears as noise that hides the signal at 2ω. In order to account for that, we use a lock-in amplifier, which locks on to the second harmonic of the driving frequency by using a reference signal taken directly from the oscillating mirror. Figure 2 shows a schematic of the experimental setup. A 0.9'μτn diameter silica microsphere suspended in water is held in the optical tweezers and observed by a camera. lOOx. NA1.25 oil immersion objective lens is used to form the trap and capture the scattered light as well as for the imaging. The laser source is an Argon ion laser operating at 488nm. The lock-in amplifier measures the phase difference between the back-scattered light signal and the reference signal at 2ω. To find the relation between the measured phase angle ψ and the particle's phase lag θ, we need to look at the measurement procedure itself. A simple position detection setup is used to generate a reliable reference signal. It consists of a He-Ne laser, which is reflected off the back side of the moving mirror to a split photodiode (see fig.2 ). A difference amplifier connected to the outputs of the split photodiode produces the desired position signal of the mirror. The lock-in amplifier receives this sinusoidally varying voltage at a frequency ω as a reference. It then generates internally a pure sine wave signal at the second harmonic:
re/2ω = sin(2^f) (7)
The signal from the photomultiplier tube (PMT) can be written as:
. sig = - sin[2w(j; - tc) + ] + noise (8)
tc represents the times when the particle and the center of the trap coincide. We add π/2 to account for the fact that the signal is high at the coincidence times. The signal is negative due to the negative bias voltage of the PMT. Note that the lock-in detection treats contributions to the signal at other frequencies as noise. We can rewrite the signal at 2ω:
sig2ω = sin(2ωt + φ) (9)
with the phase angle φ given by:
φ = -2ωtc - (10)
φ is the measured quantity in our experiment. To find tc we need to solve the equation x(t) = p(t) - po sm' ωt for t. The solution is given from equation 5 by:
tc = - t&n' 1 (cot(θ)) (11) ω which gives us the relation between θ and φ:
Finally, using equation 6, we get the relation between φ, r and b which forms the basis of our measurement technique:
,φ + . ωb , . cot(^) = ωr = - (13)
3. Results
Figure 2 shows the results of measurements of φ at different ω's at a constant laser power of 1.3 mw (power at the sample). Similar measurements were made at different laser powers of w=218.6 rad/sec ( 34.8Hz ) . All measurements were made on the same particle at the same distance from the cover glass ~ 50μm. The data clearly show that relation 13 holds.
From the slope of the graph in figure 3 we get b = 1/ (slope -ω) = 2.20±0.0lμJ. From the slope in figure 2 we get r = 1.68 ± 0.01ms, which gives b = τ - P = 2.18±0.02 zJ. Each measurement of φ is a statistical average over a period of ~ 30sec and the error bars are derived from the standard deviation. The slopes were calculated using a linear regression algorithm that forces the curve to go through the origin. We further calculate a = 5.0 • 10~3pn/nm from our measurement of τ using a Stokes drag coefficient of β = 6πηa = 8.48- 10-9 kg/s, using the viscosity of water η = 1 • 10-3 kg /mj sec.
4. Summary
The results reported above demonstrate a new analytical method method which combines optical tweezers with confocal detection to probe the microscopic scale environment and dynamics of microparticles. The results confirm the validity of the point particle model in an oscillating harmonic potential. While the multiple measurements taken over a range of frequencies and optical powers are necessary to validate the measurement technique, actual measurements of τ require in principle only one measurement at a given frequency and power. This is important because it allows one to monitor time variations of these properties in a time varying medium (e.g., a polymer solution undergoing a phase transition). If this method is used in conjunction with another method to get α then β can be determined absolutely.
A significant advantage of this technique relative to prior art methods is that it can be used with a relatively low laser power, since the optical tweezers does not need to trap the particle with a strong force, but rather has to exert enough force to move the particle in its direction. We note that there is a direct relationship between the range of frequencies and optical powers one should use for this measurement to succeed. For a given optical power, as the frequency of the motion is increase there is a point where the particle's response becomes too slow to follow the trap. At this point the lock frequency breaks down because Brownian motion starts to dominate. At a very low frequency of motion the particle with move together with the trap, and the lock frequency will break down because the second harmonic component of the signal will be very weak and noisy. This does not mean that one could not use this method at very high or very low frequencies, it only means that the strength of the trap has to be adjusted appropriately at these extremes. Finally, this method can be used as a new confocal imaging technique by obtaining successive measurements of τ at different locations in the sample.