When using this technique, it is essential that the bead solution is properly diluted and the beads are chosen based on desired diameter, PA gel substrate stiffness, and the size scale of the phenomena that is explored in the desired experiment.
Caution should be taken when diluting the bead solution prior to functionalizing top glass cover slips. The spacing between beads on the gel surface can be altered by changing the dilution factor of the colloid solution. Diluting the solution too much will reduce the number of beads that couple to the cover slip and ultimately reduce the number of beads in the gel. Diluting the solution too little can yield too high a density of beads in the gel such that they are in contact and cannot be distinguished from one another. While it is often impossible to completely avoid the presence of nanoparticle aggregation, the application of ultrasonic waves to the diluted solution effectively minimizes nanoparticle aggregation. We have found that a 10,000 fold dilution of the colloid solution utilized in this protocol gives a bead concentration of 3.64 x 108 beads/ml, and yields a desirable particle density of 0.05 beads/μm2 for 0.1 μm diameter beads. A 1:20,000 dilution (1.82 x 108 beads/ml) yields an average of 52 μm2 per particle. Figure 2 shows images of the surface of gels with different dilutions. It should be noted that there is variability in the number of beads that transfer to the gel. This is due to the variability associated with removing the bead solution from the cover slip with compressed air during functionalization. In general, we have found that increasing the dilution factor by at least a factor of two consistently results in at least a twofold increase in the area per bead.
It is also critical to choose a bead size (diameter) that will result in full immersion of the bead inside the gel upon PA polymerization. Gel stiffness affects the location of beads as shown in Figure 5A. For 0.1 μm diameter beads, the bead depth increases with decreasing stiffness. This could be a result of the inverse relationship between pore size and PA gel stiffness13. For all of the three stiffness gels tested, the beads are consistently localized inside the gel, and within 1 μm of the gel surface as illustrated in Figure 8A. The average location of 1 μm beads within the gel also varies with stiffness. They are embedded deeper (approximately 1.25 μm) into the 1 kPa gels. For the 10 and 40 kPa gels, the location of the top of 1 μm diameter beads ranges from approximately 0.5 μm above to 0.5 μm below the gel (Figure 5) surface, which means the beads occasionally protrude through the top of the gel (Figure 8B). This is detrimental because it introduces topography for the cells, and topography is known to affect cell-substrate interactions14. Thus, it is necessary to be mindful of choosing beads for this technique when preparing gels of a specific rigidity. Of the two bead sizes utilized to characterize this method, 0.1 µm diameter proves the optimum size. Other bead sizes may be utilized in with this technique, but should be calibrated using confocal microscopy as we describe here, to understand the precise depth to which the beads are localized during gel polymerization.
This technique improves spatial resolution of beads in PA gel by increasing the number of distinguishable particles within the field of view. Localizing the beads to a layer near the top of the gel does not significantly affect the composite stiffness of the gel. Isostress composite theory is utilized to compute the composite elastic modulus, Ec, of only the 1 μm-thick layer of gel containing beads as

where Ef is the elastic modulus of the polystyrene beads (3 GPa), Em is the elastic modulus of the gel matrix (10 kPa), and vf and vm are the volume fractions of the beads and gel matrix, respectively. Because the beads comprise only 0.03% of the volume, Ec = 10.0025 kPa is not significantly different than the elastic modulus of an isotropic gel of stiffness 10 kPa. An established method for fabricating PA gels reduces the polymerization time to 30 sec, thus minimizing the effect of gravity on the distribution of beads during polymerization15. Our method relies on the longer (~30 min) polymerization time for the beads to polymerize into the gel such that they do not protrude the gel surface. Because localizing the beads to a single layer does not significantly affect the composite gel stiffness, and the longer polymerization time aids in the transfer of beads into the gel matrix, we do not incorporate any measure to reduce polymerization time in our technique.
Our technique demonstrates how to localize the beads to a known depth inside PA gel, which improves the measurement of substrate deformation, and thus the calculation of cell traction forces. A theoretical approach to quantify the cellular traction forces has been previously reported16,17. This approach treats the gel substrate as a semi-infinite solid in elastic half space. The Boussinesq solution is used to solve the inverse problem of calculating traction forces from a displacement field. We employ a simplified version of the Boussinesq integral to quantify the error in displacement measurement due to beads distal from the gel surface18. The equation for the displacement in a given direction based on an applied shear stress, τ, over a given radius, R, is given by

where ux is the displacement in the x direction, μ is the shear modulus, and z is the depth of the beads from the gel surface. For various locations of beads throughout the depth of the gel, corresponding displacements are shown in Table 2.
For the case shown here, the radius is estimated to be an average focal adhesion size of 1 μm. If the beads are assumed to be coplanar with the cells (z = 0), the displacement is given by 0.75 τ/μ. If the beads are only 1 μm below the gel surface, the displacement is given by 0.24 τ/μ. This represents greater than a threefold reduction in displacement from a bead coplanar with the surface to a location 1 μm below. This introduces error in the calculation of surface traction forces based on the motion of embedded beads whose precise locations within the gel are unknown. The variation in traction defined by τ/μ, as a function of the depth of the beads in the gel for both the proposed and traditional fabrication methods is significant. Figure 9 shows the distribution of forces (given by the ratio τ/μ) based on Boussinesq theory that result from computation utilizing the distribution of bead depths, z, within the optical depth of field when beads are uniformly distributed in the gel and localized near the surface. The most accurate measurement of cell force would result if beads were on the same plane as the cells, where z = 0 and τ/μ = 0.0133. When compared to this case, τ/μ increases by 118%, 80%, and 50% for 1 kPa, 10 kPa, and 40 kPa gels, respectively, when beads are localized near the surface using the technique we propose. When beads are uniformly distributed throughout the gel depth, τ/μ increases an additional 60%, 80%, and 150% beyond the increase induced by localizing beads near the surface, for 1 kPa, 10 kPa, and 40 kPa gels, respectively. This represents a total of nearly 200% error in traction calculations when the beads are distributed throughout the gel depth. Thus, localizing the beads to a known depth inside the gel improves accuracy in displacement measurements and cell traction forces.
The technique described here presents a new way to prepare PA gel substrates for traction force microscopy applications such that the beads are localized near the gel surface. Confining the beads to a known depth within PA gel provides new information about the actual displacement of beads upon cell force application, creating a more accurate picture of how the cells are deforming their underlying substrate.