$$\rightleftharpoonup{xx}$$
$$\longleftharp{xx}$$,
$$\longrightharp{xx}$$,
The development of therapeutic antibodies has shown remarkable success in the treatment of multiple tumor types1. The recent advancements of antibody-drug conjugates (ADC) as delivery mechanisms for cytotoxic compounds has expanded the requirements for understanding the dynamics of antibody:receptor interactions at the cell surface2. Following the successful targeting of an antibody to the cell surface receptor, these complexes can induce similar aggregation patterns to those observed in ligand:antibody interactions3. Alterations in receptor aggregation can induce changes to the membrane and result in the internalization of the receptor and its removal from the cell surface. In the context of an antibody-drug conjugate, this process subsequently releases the cytotoxic payload into internalized endosomes and subsequently the cytoplasm, resulting in effective cell killing.
Confocal microscopy has provided an effective means of visualizing these important interactions of antibodies and their target receptors4. To explore the changes of aggregation of target molecule at the cell surface this protocol utilizes post processing of confocal microscopy images via a spatial image correlation spectroscopy (ICS) technique 5,6,7.
The foundation of image correlation spectroscopy is the observation that spatial fluorescence intensity fluctuations share a relationship to the density and aggregation state of the labeled structures. This relationship is established following the calculation of a spatial autocorrelation function of a captured image5.
All variants of image correlation spectroscopy require the calculation of an image autocorrelation. This is followed by fitting this function to a two-dimensional Gaussian curve for the extraction of quantitative aggregation state parameters contained within the image. In simple terms the calculation of an image autocorrelation involves comparing all the possible pixel pairs contained within an image and calculating the likelihood that both equally as bright as each other. This is visualized as a function of the distance and directions of pixel separation8.
The theoretical framework for the image correlation spectroscopy was established and defined by Petersen and Wiseman et al.5,6. In this protocol, the autocorrelation calculations are performed in Fiji/ImageJ as well as a spreadsheet application, the basis for the intensity fluctuation spatial autocorrelation function can be described as (Eq 1):
where F represents the Fourier transform; F−1 the inverse Fourier transform; F* its complex conjugate; and the spatial lag variables ε and η. In spatial ICS, as described in this protocol, the autocorrelation function can be calculated using a 2D fast Fourier transform algorithm7,9. The autocorrelation at zero spatial separation otherwise known as zero-lag, g11(0,0), provides the inverse mean number of particles present per beam area of the microscope. It can be obtained by fitting the spatial autocorrelation function to a two-dimensional Gaussian function (Eq 2):
As the pixels captured within an image are contained within a set area and these measurements do not extend to infinity, the term g∞ is used as an offset to account for long-range spatial correlations contained within the image. For molecular-sized aggregates, ω is the point-spread function of the microscope and described by the full-width at half-maximum of the spatial autocorrelation function. The area contained within the point spread function of the instrument can be calibrated through the use of sub-resolution fluorescent beads.
For the image correlation spectroscopy protocol described herein, the autocorrelation and mathematical functions required to complete ICS are performed using the open-source imaging-processing platform, Fiji10, a distribution of the ImageJ program11,12. Fiji/ImageJ utilizes the preinstalled fast Fourier transformation in the FFT Math function. This function reduces the compute time required of this calculation by reducing the range of data by a factor of two in each dimension13. As the 2D autocorrelation function is approximately symmetrical in x,y axis, a single line profile plot through the autocorrelation image can be used to measure the raw autocorrelation as a function of the spatial lag. Any zero-lag noise is removed prior to further calculation, with the resulting autocorrelation amplitude (peak value, g(0)T) corrected for background with the expression (Eq 3):
where Ib is the mean intensity from a background region excluding the cell. The cluster density, or density of fluorescent objects, is defined by (Eq 4):
In the protocol described herein, we further simply the calculation of cluster density (CD), with the assumption based on the observation that a normalized autocorrelation function will decay to a value approaching 1.0 with increasing spatial lag. With maximum spatial lag, there is no longer any correlation of fluorescence intensity values and thus without a correlation the calculations at this region are computing the value of an intensity multiplied by this intensity which is subsequently divided by the square of that intensity, which by definition is equal to 1.0. Thus, cluster density from a normalized autocorrelation function can be computed by subtracting 1.0 from the normalized autocorrelation function prior to taking its reciprocal (Eq 5):
Further calibration of the beam area can be performed to quantitate the number of clusters contained within the area of the point spread function of the instrument. This calibration must be performed using the same optical conditions used during the image correlation spectroscopy analysis.