Here, we show the use of traditional dark-field microscopy to monitor the dynamics of gold nanorods (AuNRs) on cell membrane. The location and orientation of single AuNRs are detected using ImageJ and MATLAB, and the diffusive states of AuNRs are characterized by single particle tracking analysis.
Analyzing the diffusional dynamics of nanoparticles on cell membrane plays a significant role in better understanding the cellular uptake process and provides a theoretical basis for the rational design of nano-medicine delivery. Single particle tracking (SPT) analysis could probe the position and orientation of individual nanoparticles on cell membrane, and reveal their translational and rotational states. Here, we show how to use traditional dark-field microscopy to monitor the dynamics of gold nanorods (AuNRs) on live cell membrane. We also show how to extract the location and orientation of AuNRs using ImageJ and MATLAB, and how to characterize the diffusive states of AuNRs. Statistical analysis of hundreds of particles show that single AuNRs perform Brownian motion on the surface of U87 MG cell membrane. However, individual long trajectory analysis shows that AuNRs have two distinctly different types of motion states on the membrane, namely long-range transport and limited-area confinement. Our SPT methods can be potentially used to study the surface or intracellular particle diffusion in different biological cells and can become a powerful tool for investigations of complex cellular mechanisms.
The dynamics of nanoparticles (NPs) on the membrane is closely associated with the cellular uptake process, which is essential for the understanding of cell functions, viral or bacterial infections and the development of artificial nanomedical delivery systems1,2. Single particle tracking (SPT) technique is a robust tool for characterizing the heterogeneous behaviors of NPs3,4. In general, cell membrane is fluidic, which means that the components such as proteins and lipids can move laterally in the plasma membrane plane5,6,7. The spatiotemporal complexity of membrane organization and structure may lead to spatiotemporal heterogeneity of the interaction between NPs and membrane. Hence, direct visualization of the movement of NPs on the membrane requires both high spatial and temporal resolution.
Single particle tracking microscopy that monitors the localization of individual particles in living cells with a spatial resolution of tens of nanometers and a time resolution of milliseconds has been well developed to study the NPs or membrane molecules dynamics8,9. Fluorescence-based microscopic imaging techniques have become valuable tools for observing NPs/molecules in living cell environment9,10,11,12. For example, total internal reflection fluorescence microscopy, which images thin layers (~100 nm) of the sample at the substrate/solution interface with a high spatiotemporal resolution has been widely used in studies of membrane molecules dynamics13,14. However, the inherent disadvantages of single fluorophores, such as low intensity and rapid irreversible photobleaching reduce the accuracy and duration of tracking13. Therefore, non-fluorescent plasmonic NPs, which replace the fluorescent probes, have attracted more and more attention in long-term imaging studies due to their unique optical characteristics15. Based on the scattering signals of plasmonic NP probes, several kinds of optical microscopic imaging technologies have been used to study the mechanism of biological processes, such as dark-field microscopy (DFM)16, interferometric scattering (iSCAT) microscopy17 and differential interference contrast microscopy (DICM)18. In addition, the motion and rotation dynamic of AuNRs can be obtained using DFM and DICM18,19,20,21,22. Typically, in an SPT experiment, the motion of the object is recorded by the optical microscope, and then analyzed by SPT analysis methods3. The time-resolved trajectories and orientational angles generated by individual NPs are normally stochastic and heterogeneous, so it is necessary to present abundant dynamic information with various analysis methods.
Here, we provide an integrated protocol that monitors the dynamics of AuNRs on cell membrane using DFM, extracts the location and orientation of AuNRs with ImageJ and MATLAB and characterizes the diffusion of AuNRs with SPT analysis methods. As a demonstration, we show here how to use the SPT protocol to visualize dynamics of unmodified AuNRs (CTAB-AuNRs, synthesized by cetyltrimethylammonium ammonium bromide molecule as protective agent) on U87 MG cell membrane. It has been demonstrated that CTAB-AuNRs can adsorb proteins in biological environment, move on cell membrane and then enter cells2,20,22. U87 MG cell is the most common and most malignant tumor of the central nervous system, and its membrane receptors are abnormally expressed. The membrane receptors can interact with proteins on AuNRs, which influence the dynamics of AuNRs. Our protocol is generally applicable to other SPT experiments in the field of biology.
1. Cell culture
2. Microscope slide preparation
NOTE: U87 MG cells of third to tenth generation with high activity are used in SPT experiments.
3. Performing single particle tracking experiments with darkfield microscopy (Figure 1).
4. Data acquisition
5. Data analysis
NOTE: A systematic and robust data analysis framework is essential for the performance and efficiency of SPT analysis methods. The custom software written in MATLAB is used (https://github.com/fenggeqd/JOVE-2020/tree/master/Analysis_parameters). A graphing and analysis software (see Table of Materials) is used for drawing the plots.
Measures | Definition | Physical meaning | |
Displacement | Changes in position of objects | ||
Step size | Distance between two adjacent points | ||
Velocity | Speed of objects motion | ||
Rg | Moving range of objects in a specific time interval | ||
Ta | Motion direction of objects between two adjacent points | ||
MSD | Average moving distance of objects in a specific time interval | ||
Dt | Diffusion ability of objects | ||
α | Normal diffusion (α~1) | ||
Polar angle | 3D orientation information of objects | ||
σ | Degree of dispersion of the polar angle data set |
Table 1: Three types of parameters used for analysis. These include trajectory related parameters (displacement, step size, velocity, Rg and Ta), MSD parameters (MSD, Dt and α) and rotation related parameters (polar angle and rotational lability).
In the protocol, the unmodified 40 x 85 nm CTAB-AuNRs were used. As shown in Figure 2B, its longitudinal plasmonic maximum at is ~650 nm (red region) and transverse resonance is at 520 nm (green region). Previous literatures have revealed that the optical properties (such as LSPR intensity) of plasmonic AuNRs will change significantly with their diameter20,22. In Figure 2C, the scattering intensity from CTAB-AuNRs on U87 cell membrane, showed typical Gaussian distribution with narrow width and was consistent with that of CTAB-AuNRs on glass, which indicates that CTAB-AuNRs tracked in this experiment were well monodispersed.
The dynamic of CTAB-AuNRs on membrane was monitored by DFM (12 fps). As shown in Figure 7, more than 500 trajectories (the trajectory length exceeds 300 frames) were obtained. Additionally, approximately 40 trajectories per cell can be generated. The Rg of all 500 trajectories showed a small-value distribution, the mean Rg was 0.5 µm (SD=0.6 µm), and max-displacement was also more distributed at small values (Figure 8). All 500 trajectories can be divided into two groups: long-range diffusion (Rg>0.5 µm, black trajectories) and confined diffusion (Rg<0.5 µm, red trajectories).
According to position and orientation of CTAB-AuNRs, several parameters are calculated for data analysis. There are many mays to analyze and present parameters during SPT analysis. The ensemble-time averaged MSD analysis shown in Figure 9A shows that CTAB-AuNRs normally diffuse with α~1. However, the density distribution of Dt-α, obtained from all trajectories, reveals that the dynamic of AuNRs showed a heterogeneous distribution with superdifusion motion, Brownian motion and subdiffusion motion (Figure 9B).
In addition to multi-particles statistical analysis, single-particle (multi-frames) statistical analysis was also performed. Figure 10 shows two representative long-term trajectories: confined (Rg=0.34 µm) and moving (Rg=1.48 µm). The polar-mapping trajectory showed that both the motion range and polar angle of confined AuNRs were smaller than those of moving AuNRs, which means that stronger the interaction between AuNRs and membrane, the more serve the confinement of its translational and rotational dynamics. Correspondingly, the histograms of Ta and polar angle of the two AuNRs showed different modes.
To further study the motion behavior of AuNRs, the time-series parameters were analyzed, including trajectory, displacement, polar angle and Rg. As shown in Figure 11, different patterns were presented over time. Compared with the long-distance diffusion AuNR, the rotation of the confined AuNR was relatively limited, and its polar angles were more concentrated, distributed in the range of 10° ~ 40°. However, although the overall Rg of the moving AuNR was larger than the overall Rg of confined AuNR, the time-series Rg of moving AuNR was smaller. Combining the time-series distribution of trajectories points and displacement, it can be inferred that confined AuNR moves in a limited area, but its motion range in a short time is larger than that of the moving AuNR confined on the membrane. Overall, the spatiotemporal heterogeneous distribution of translational and rotational dynamics is caused by the spatial and temporal complexity of membrane organization and structure.
Figure 1: Schematic diagram of the optical path of darkfield microscopy, microscope slide preparation, and the dual-channel difference method for calculating polar angle (θ) of AuNRs. Scale bar in image is 2 μm. Please click here to view a larger version of this figure.
Figure 2: Characterization of CTAB-AuNRs. (A) TEM images of CTAB@AuNRs. (B) LSPR spectrums of CTAB@AuNRs. (C) Distribution of single particle intensity: CTAB@AuNRs on glass (left panel), CTAB@AuNRs on plasmon membrane (right panel). Please click here to view a larger version of this figure.
Figure 3: Pretreatment of darkfield images. Time-series color images were converted to “8 bit” mode and their contrast were adjusted. Please click here to view a larger version of this figure.
Figure 4: Extraction of single long-term trajectory with ImageJ. Please click here to view a larger version of this figure.
Figure 5: Extraction of multi-trajectories with Fiji. Please click here to view a larger version of this figure.
Figure 6: Extraction of R/G values of AuNRs with MATLAB. Please click here to view a larger version of this figure.
Figure 7: All 527 trajectories of CTAB@AuNRs on U87 MG cell membrane tracked with Fiji. There were two groups: red trajectories (Rg<0.5 μm) and black trajectories (Rg>0.5 μm). Please click here to view a larger version of this figure.
Figure 8: Multi-particles statistical analysis of all 527 trajectories. Distribution of Rg (A) and maximum displacement (B). Please click here to view a larger version of this figure.
Figure 9: MSD analysis of all 527 trajectories. (A) Plot of ensemble-time averaged MSD versus τ. (B) Density plot in the α– log (Dt) plane. The color code represents the local density of trajectories between 0 (blue) and 1 (red). Please click here to view a larger version of this figure.
Figure 10: Single-particle statistical analysis. (A-B) Two representative long-term polar-mapping trajectories which were tracked with ImageJ. The color code represents the polar angle value between 0 (green) and 90 (red). Distribution of turning angle (C-D) and polar angle (E-F). Please click here to view a larger version of this figure.
Figure 11: Time-series analysis. (A,C) Time-mapping trajectories. The color from red to blue represents the direction of time. (B,D) Time-series of parameters: displacement (blue), polar angle (green) and radius of gyration (Rg, red). Please click here to view a larger version of this figure.
The presented protocol is used to study the dynamics of AuNRs on cell membrane. The protocol consists of four parts, including microscopic imaging, data extraction, dynamic parameters calculation and data analysis methods, and each part is flexible and universal. Therefore, there are many possible future applications, for instance, studying movement of NP-linked membrane molecules on membrane, endocytosis dynamics of NP-labeled receptors, dynamic analysis of intracellular NPs and vesicle-coated NPs transportation along microtubules, and so on.
The basic step is to use DFM to image the movement of AuNRs on the cell membrane. Currently, many imaging technologies for the analysis of cell membrane dynamics have been developed, of which fluorescence microscopy has been widely used24. However, the photobleaching properties of fluorescent probes lead to the inability to track particle dynamics for long periods of time. Here, light-stable plasmonic AuNRs were used. Several imaging techniques have been developed based on the scattering characteristics of plasmonic NPs15. Particularly, DFM with oblique illumination mode only collects scattered light from the samples, which makes it have a high signal-to-noise ratio.
Both the translational and the rotational dynamics of AuNRs are obtained and analyzed. Most studies focus on investigating the positional fluctuations of the NPs on the membrane, and ignoring the effects of orientational changes25. In addition, the cell membrane plane and the imaging plane (x-y plane) are largely the same in most cases, and the difference in the interaction between the AuNR and the cell membrane is more notable in the z direction. Hence, even without azimuth angles, the polar angle is valuable for the analysis of the interaction between the AuNRs and cell membrane. Based on DFM and unique optical properties of AuNRs, dual-channel polarization DFM is used to obtain the orientation information of AuNRs (I1 + I2 sin2 θ)19,26. Nevertheless, the calculation accuracy of polar angles is affected by the maximum and minimum detection intensities which are strongly influenced by surroundings. Hence, some modifications need to be made. Here, (R-G) was used to calculate the polar angle, and the G value serving as an internal reference could efficiently decrease the systematic errors and increase the measurement accuracy.
The critical steps are data extraction and dynamic parameter calculation, which is essential for analyzing the dynamics of NP on the membrane. The trajectories of AuNRs on membrane are extracted with high spatial resolution27,28 from sequence analysis using ImageJ/Fiji. The RG values and polar angles are calculated with MATLAB (https://github.com/fenggeqd/JOVE-2020). Several parameters are similarly calculated using MATLAB according to positions and polar angles (Table 1). In SPT analysis, there are many ways to perform parameters analysis and display. Hence, it is necessary to use as many analysis methods and visualization tools as possible to analyze the dynamics of NP, and then to summarize and extract a novel point of view from many analysis results. Normally, there are always some trade-off to make, which is not an easy work.
The dynamic state and motion behavior of NPs can be analyzed and presented systematically from top to bottom using comprehensive analytical methods and multivariate presentation ways. Using multi-particles statistical analysis, we can obtain the ensemble motion state of NPs across the whole cell membrane plane. The dynamic studies on individual NPs/molecules based on fluorescence imaging29 mostly use this overall statistical method because only a large number of short-term trajectories can be monitored. However, the ensemble averaging method will smooth and ignore individual details. As a supplement, the single-particle statistical analysis is used to quantitatively present the dynamic state of individual NPs. Meanwhile, it is important to perform time-series parameter analysis of individual or segment trajectories, which can provide motion behaviors of NPs as a function of time.
In summary, we propose an integrated protocol for studying the diffusion dynamics of AuNRs on live cell membrane with SPT method. The dynamic of AuNRs was monitored with single nanoparticle DFM, extracted using ImageJ and MATLAB, and characterized by comprehensive analytical methods. AuNRs did abnormal diffusion motion on U87 MG cell membrane, and its motion shows spatiotemporal heterogeneity. This protocol can be potentially used for the studies of other type of complex biological systems.
The authors have nothing to disclose.
This work was supported by the National Natural Science Foundation of China with grant numbers of 21425519, 91853105 and 21621003.
CTAB coated gold nanorods(CTAB-AuNRs) | Nanoseedz | NR-40-650 | 85 nm * 40 nm |
Color CMOS camera | Olympus | DP74 | Japan |
Coverslips | Citoglas | z10212222C | 22*22 mm |
Dark-field microscopy | Nikon | 80i | upright microscope |
Fetal bovine serum (FBS) | Gibco | 10099141 | |
Fiji | National Institutes of Health (NIH) | 2.0.0-rc-69/1.52 p | a distribution of ImageJ |
Grooved glass slide | Sail brand | 7103 | Single concave |
Image J | National Institutes of Health (NIH) | 1.52 j | |
MATLAB | MathWorks | R2019b | |
MATLAB Code | https://github.com/fenggeqd/JOVE-2020 | ||
Minimum essential medium (MEM) | Gibco | 10-010-CVR | with phenol red |
Minimum essential medium (MEM) | Gibco | 51200038 | no phenol red |
Origin | OriginLab | Origin Pro 2018C | |
Penicillin-streptomycin | Gibco | 15140122 | |
Plastic cell culture dishes | Falcon | 353002 | |
Plastic cell culture dishes | Falcon | 353001 | 35*10 mm |
U87 MG cell | American Type Culture Collection | ATCC HTB-14 | a human primary glioblastoma cell line |