Mahdiye
Ijavi
a,
Robert W.
Style
a,
Leonidas
Emmanouilidis
b,
Anil
Kumar
c,
Sandro M.
Meier
bc,
Alexandre L.
Torzynski
a,
Frédéric H. T.
Allain
b,
Yves
Barral
b,
Michel O.
Steinmetz
c and
Eric R.
Dufresne
*a
aDepartment of Materials, ETH Zürich, Switzerland. E-mail: eric.dufresne@mat.ethz.ch
bInstitute of Biochemistry, ETH Zürich, Switzerland
cLaboratory of Biomolecular Research, Paul Scherrer Institute, Switzerland
First published on 15th December 2020
Phase separated macromolecules play essential roles in many biological and synthetic systems. Physical characterization of these systems can be challenging because of limited sample volumes, particularly for phase-separated proteins. Here, we demonstrate that a classic method for measuring the surface tension of liquid droplets, based on the analysis of the shape of a sessile droplet, can be effectively scaled down to measure the interfacial tension between a macromolecule-rich droplet phase and its co-existing macromolecule-poor continuous phase. The connection between droplet shape and surface tension relies on the density difference between the droplet and its surroundings. This can be determined with small sample volumes in the same setup by measuring the droplet sedimentation velocity. An interactive MATLAB script for extracting the capillary length from a droplet image is included in the ESI.
Liquid–liquid phase separation of macromolecular solutions has enjoyed intense interest in recent years thanks to the discovery that many membrane-less compartments within eukaryotic cells appear to be liquid droplets of phase-separated macromolecules.5,6 Examples include P granules,7 P bodies,8 Cajal bodies9 and nucleoli.10 Membraneless organelles have complex compositions, including a large number of distinct RNAs and proteins. The bulk of the proteins found in these droplets are intrinsically-disordered and lack globular structure.11–14 These proteins readily phase separate in vitro to form a dilute ‘supernatant’ phase and a protein-rich phase, which can take the form of liquid droplets, hydrogels, fibrils or aggregates.15 The tendency of intrinsically disordered proteins to phase separate is thought to underlie the formation of membraneless compartments in the cell.16–18 Precise physical characterisation of droplets in vitro can provide insights into the molecular scale interactions of their components, and potentially shed light on their physiology in vivo.
The essential challenge in characterizing the properties of these droplets is their size. Typically no more than a few microns across in vivo, larger droplets can be formed in vitro. However, the production and purification of protein is labor intensive and expensive, and total volumes of the protein-rich phase are typically limited to the μL scale. Such samples are too small for conventional mechanical measurements. For example, quantification of bulk mechancial properties through conventional rheometry requires at least 100 μL of material. In the last twenty years, a number of new approaches have been developed to scale-down rheological measurements.19 Among these, particle tracking microrheology is well-suited to measure the linear rheology of phase separated droplets.20,21
Similarly efficient and effective approaches have not yet been established to characterize the surface mechanical properties of droplets. Classic spinning drop measurements can work for droplets on the μL scale, but only if they are less dense than their surroundings. In principle, the spinning droplet method could be applied to a small droplet of the protein-depleted phase within a larger volume of the protein-rich phase, but that would require volumes of the protein-rich phase around 1000 μL. A number of methods have been introduced in recent years to estimate the surface tension of micron-scale phase separated droplets. The most popular method is based on measuring the time it takes for two droplets to fuse.22 In principle, this method can be very accurate if the rheology of the droplets is known. While fusion times are easy to calculate from first principles in bulk, they can be strongly modified when droplets are in contact with a surface, as contact lines can dramatically slow down the movement of any droplet.23 Imaging of fusion events in the absence of contact lines is very challenging. At droplet volume fractions low enough to clearly resolve the dynamics, fusion events are rare. Optical tweezers can be used to initiate fusion, but care has to be taken to avoid contributions of optical forces to the fusion time as well as photo- and thermal damage to the droplets.24 Measurements of the force required to statically deform a droplet can also be used to accurately determine the surface tension using optical tweezers25 or atomic force microsopy.26,27
Alternatively, the shape of droplets of an appropriate size contains sufficient information to accurately determine their surface tension. For conventional simple liquids, made of small molecules, this is done through the analysis of photographs of millimeter scale pendant or sessile droplets.28,29 Surface tension is most easily inferred from the shape of a droplet when its radius is larger than the capillary length, . Here, Δρ is the density difference between a droplet and its surroundings and g is the acceleration due to gravity. The tiny values of surface tension for phase-separated droplets of macromolecules are somewhat mitigated by a simultaneous reduction of the density difference (e.g.ref. 30). All together, this leads to reported values of the capillary length varying from
(10 mm) to
(100 mm).4 This corresponds to droplets with volumes of the order
(nL). By analyzing the shape of fused fluorescently tagged nucleoli from giant Xenopus oocytes, Feric et al. determined their capillary length to be of order
(10 mm).31 One challenge of this approach is that the density difference between the droplet and surrounding fluid needs to be known accurately. This can be inferred from sedimentation speeds of microscopic droplets when the viscosities of the droplet and surroundings are known.32
Here, we apply established methods of sessile droplet shape analysis to determine the surface tension of phase-separated macromolecular droplets with unknown density. Using a standard camera equipped with a telecentric lens, we quantify phase separated liquid droplets, made of polymers or proteins, with surface tensions ranging from 7 to 90 μN m−1. We determine the density difference between the droplet and continuous phases in the same setup by measuring the droplets' sedimentation velocity. Typical values of the density difference range from 30 to 150 kg m−3. This approach is attractive for the in vitro measurement of phase-separated proteins because it uses small sample volumes and the proteins do not need to be fluorescently tagged.
To image droplets, we use a typical optical setup for sessile droplet analysis,33,34 described in the Materials and methods. Droplets are illuminated by an extended light source and imaged with a telecentric lens. To limit artifacts due to refraction, the droplet and its co-existing continuous phase are held in a rectangular cuvette. A typical image of dextran-rich droplet in its co-existing dextran-poor continuous phase is shown in Fig. 1. A series of droplets of different sizes of the same material are shown as insets to Fig. 2. The droplet shape, rexp(z), is manually extracted from the acquired images in MATLAB, using the script included in the ESI.†
![]() | ||
Fig. 2 Sessile droplet shape analysis. Images of sessile dextran-rich droplets in a dextran-poor continuous phase. The scale bars are 200 μm. Dashed-green lines on the droplets shown best fit curves to the droplet shape. These images are superimposed over a scatter plot of the Neumann number and aspect ratio, H/d. The dashed-grey curve is Neumann number in the perfect case where the sessile droplet has a contact angle of 180°, calculated numerically using eqn (1)–(3). |
Large droplets, appropriate for surface tension measurements, described in Section 3, can be formed by direct pipetting of the droplet phase into its coexisting continuous phase. Alternatively, a droplet dispersion can be pipetted into the same continuous phase. The latter approach allows for a wider range of droplet sizes, including smaller droplets appropriate for sedimentation velocity measurements of Section 4.
To accurately determine the 3D shape of the droplet from such 2D images, droplets must be axisymmetric. As we will see in Section 3, the accuracy of the surface tension measurement is best for droplets with a large contact angle. Phase separated droplets of proteins often interact strongly with standard glass and plastic surfaces, spreading to very low contact angles and adopting non-axisymmetric shapes due to pinning. Therefore, the bottom of the cuvette needs to be very clean and functionalized to minimize spreading. A PEG-silane treated coverslip placed at the bottom of the cuvette is sufficient for the droplets studied here. The protocol for creating this is described in the Materials and methods.
![]() | (1) |
The coordinate system, including s, r, z, θ, is defined in Fig. 1. The left-hand side of the equation gives the contribution from gravity, which increases linearly from the top of the droplet. The right-hand side of the equation gives the contribution of surface tension to the droplet pressure, which is determined by the local curvature. The constant P* = 2γ/R0 is the pressure at the top of the droplet, where R0 is the local radius of curvature.
To simplify the equation, we scale r, s and z by the capillary length (indicated by over-bars) and rearrange to obtain
![]() | (2) |
![]() | (3) |
The last equation holds at the top of the drop, and represents the fact that the two principle surface curvatures are equal at this point.35
To extract the material properties of the droplet, we find the values of β and Lc that minimise the difference between the numerical solution, and the non-dimensionalised shape of the droplet, rexp/Lc. Examples of the fits to droplet shape superimposed on raw droplet images are found in the insets to Fig. 2. An interactive MATLAB script for fitting droplet shape is found in the ESI.†
Previous work has shown that sessile droplet measurements of the surface tension are only accurate when the Neumann number, Ne = R0H/Lc2 = H/(Lcβ) > 0.3, where H is the height of the droplet from the substrate.36 For smaller Ne, droplets are effectively spherical, and it is not possible to extract a meaningful value of γ. Practically speaking, Ne can be increased by increasing droplet size, or by increasing the contact angle of the droplet. A scatter plot of the Neumann number and droplet aspect ratio for a series of dextran-rich droplets in a dextran-poor continous phase are shown in Fig. 2.
The density difference between the two phases can easily be determined by measuring the sedimentation speed of small droplets in the same setup used for the sessile droplet shape experiments. Example data for the PEG-dextran system is shown in Fig. 3. Images of a sedimenting dextran-rich droplet, acquired at equally spaced time intervals are shown in the inset. As shown by the straight line through their centers, these droplets move at a steady sedimentation velocity. The sedimentation speed depends on the droplet size, and is plotted over a factor of eight in droplet radius in Fig. 3.
The sedimentation velocity, v, is determined by a balance of the buoyant force and viscous drag. For a spherical droplet of radius R, the buoyant force is 4πΔρgR3/3. When the droplet viscosity, ηd, is much larger than the viscosity of the continuous phase, ηc, the viscous drag is given by the familiar Stokes form, 6πηcRv. Balancing the forces and solving for the density difference, we find
![]() | (4) |
If the continuous phase viscosity is comparable or larger than the droplet viscosity, the previous equation is corrected by a factor that depends on the viscosity ratio:37
![]() | (5) |
Note, however, that this correction is never more than a 33% change in density.
Since macromolecule-rich fluids tend to be much more viscous than simple liquids, eqn (4) is usually quite accuruate for protein and coacervate droplets in their typical buffers. However, for systems with signficant macromolecule concentrations in both phases eqn (5) may be necessary. In such cases, independent measurements of the viscosities of both phases are suggested. For the PEG-dextran example of Fig. 3, ηc/ηd = 0.05, so that eqn (4) is accurate to <2%.
For fixed compositions of the droplet and continuous phases, the sedimentation velocity should scale like the square of the droplet radius. This can be verified by rearranging eqn (4) or (5). The dashed-line in Fig. 3 shows a good fit of the sedimentation speed to the R2 form over this range of droplet radii. Systematic deviations from this size-dependence can indicate that the droplets are not an appropriate size.
On the one hand, droplets need to be significantly bigger than the resolution of the imaging system. On the other hand, the droplets have to be small enough that they maintain a spherical shape while falling. To remain spherical, droplets require the capillary number, Ca = vηc/γ ≪ 1.38,39 Inserting this into eqn (4), this requirement becomes that R ≪ Lc. This criterion contradicts the requirement for accurate sessile droplet measurements, R ≳ Lc. Therefore, one should not use the same droplets for sedimentation velocity and surface tension measurements. For that reason, we prefer to pipette a premixed emulsion into the continuous phase. However, when multiple droplets fall near each other, hydrodynamic interactions can dramatically effect their sedimentation velocity.40,41 Therefore, sedimentation measurements should be limited to extremely dilute, well-separated droplets.
The sedimentation of large droplets can also be impacted by inertial effects, quantified by the Reynolds number, Re = ρcvR/ηc. Eqn (4) and (5) are only valid when Re ≪ 1. For the largest expected density differences (around 500 kg m−3) in the least viscous continuous phase (an aqueous continuous phase with no macromolecules, ηc ≈ 10−3 Pa s), inertial effects will be negligible for droplets less than 100 mm in radius.
Measured values of the surface tension and relative density of each of these droplets is shown in Fig. 4. The surface tensions of these systems, shown in Fig. 4a, range from about 7 to 90 μN m−1. For each measurement of the capillary length, the Neumann number was greater than 0.3, as shown in Fig. 7. The density differences between the two phases, shown in Fig. 4b, range from 30 to 150 kg m−3. The surface tension and density values for PEG/dextran agree well with literature data.4 To connect sedimentation speeds of the FUS and Bik1 droplets to density differences, we assumed that the continuous phase viscosity is the same as water, and that the inner phase is much more viscous, the latter is justified by the microrheology measurements in Fig. 5. For dextran droplets, we measured the viscosities of both phases with a rheometer, which were 150 ± 50 and 3.9 ± 0.1 mPa s, for the dextran-rich and dextran-poor phases respectively.
While the two protein systems have similar capillary lengths, around 160 μm, they have very different surface tensions and densities. The surface tension of Bik1 is 7 μN m−1, much smaller than the surface tension of FUS267, about 90 μN m−1. There are no measurements of the surface tension of these specific protein droplets. The surface tension of other phase-separated protein droplets measured through other means have been reported to range from 0.4–100 μN m−1.7,22,25,31,45 Bik1 droplets had a significantly smaller density mismatch than FUS267 droplets, about 29 kg m−3 and 135 kg m−3, respectively. Using a nominal protein density of 1300 kg m−3, these density differences imply respective protein volume fractions within the droplets of 10% and 45%. There are few measurements of the densities of phase-separated protein droplets. Feric et al.32 measured the density of the nucleolus (composed of a complex mixture of proteins and nucleic acids) using a similar sedimentation method, and found a value of 1140 kg m−3. This is similar to our measured value of FUS267 droplets, about 1135 kg m−3.
These values are placed in a broader context using an Ashby plot of surface tension and density in Fig. 4c. In addition to these three droplets, we include a selection of droplets of far-from-critical simple liquid combinations. This plot higlights an essential feature that contrasts macromolecular droplets (shaded light green), from simple liquids (shaded gray): macromolecular droplets tend to have very low surface tensions, even when they are far from critical. Contours of the capillary length, Lc, are indicated by the dashed lines in Fig. 4c. While simple liquids have capillary lengths from the millimeter to centimeter scales, our macromolecular droplets have capillary lengths between 100 μm and 1 mm.
To complete our mechanical characterization, we measured the rheological properties of the protein droplets using standard particle tracking microrheology,19 as described in the Materials and methods. Mean-squared displacements of ≈0.1 μm radius tracer particles embedded in each of the three droplet phases all show a linear dependence on time, as shown in the inset of Fig. 5. This confirms that each droplet behaves as a simple Newtonian liquid, with no signs of viscoelasticity. Since the mean-squared displacements are linear, the rheology of each of the droplets is simply characterized by their viscosity. The viscosities of the Bik1 and FUS267 droplets are 18.2 and 1.2 Pa s, respectively.
An Ashby diagram showing the surface tension and viscosity of droplets is shown in Fig. 5b. It helps to illustrate potential advantages of the sessile droplet method over the popular ‘fusion method’ for determing the surface tension. In the latter case, the time it takes for the droplets to coalesce at low Reynolds number is set by the ratio of the droplet size and the capillary velocity Vc = γ/ηd. Thus, with a knowledge of ηd, γ can be extracted. If the capillary velocity is too high, fusion events cannot be resolved. The dashed lines on Fig. 5b indicate contours of constant capillary velocity. While the fusion of Bik1 droplets in the lower right corner of this diagram can be easily resolved, microscopic fusion experiments become much more challenging as droplets move to the top left of the diagram. The capillary velocities of dextran and FUS267 droplets are a hundred-fold higher, around 100 mm s−1. For typical 10 mm droplets, resolution of fusion events would require frame rates around 100 frames per second.
Compared to the popular droplet fusion technique, the sessile droplet method offers several advantages. First, it can measure a broader range of droplets, because only the most viscous droplets will have slow enough fusion events to be resolved with standard imaging systems. Second, it is not subject to systematic errors associated with contact line sliding.23,46 Third, viscoelasticity is expected to change the dynamics of droplet fusion. Sessile droplet shape analysis can work with viscoelastic droplets, as long as they have no elastic memory at long times.
One challenge of the current approach is that the droplets should have as large a contact angle as possible. To achieve Ne > 0.3, one typically requires contact angles greater than 90°. Larger contact angles will also typically reduce pinning, which is important so that the droplets remain axisymmetric. While a simple PEG-ylated surface was sufficient for the cases studied here, different surface treatments may be needed for other droplets. In cases where appropriate surface treatments are not readily available, a similar shape analysis of pendant droplets could be more convenient.
To coat the coverslips, we immerse them using a glass staining rack in a constantly-stirred solution containing 2300 mg of PEG silane, 500 mL of toluene, and 800 μL of aqueous hydrochloric acid (hydrochloric acid fuming 37% from VWR, MilliporeSigma). After 18 hours at room temperature, the coverslips are rinsed once in toluene (toluene 99.7+% from Alfa Aesar), and then twice more with ethanol (ethanol, denatured with IPA from Alcosuisse), before being dried with clean, dry air. Finally coverslips are stored in a dry chamber with a dessicant until required for use.
The cuvette is fill with about 200 μL of the continuous phase, which is separated from the co-existing droplet phase by centrifugation (10 min at 14 rcf). We consistently use the denser phase for sedimentation experiments (enriched in dextran or protein). Droplet phase or a coexisting dilute mixture of droplets in the continuous phase are added from the top with a micropipette. For sedimentation experiments, a series of images is recorded as the droplets fall, typically around 3 fps. For sessile droplet experiments, one must wait for the droplet to reach its equilibrium shape after coming in contact with the surface, which can take several minutes and depends on the size of the droplet and its capillary velocity.
![]() | ||
Fig. 7 The Neumann number calculated for all the sessile droplets of the 3 systems. Note that they are all greater than 0.3. |
Samples are imaged with fluorescence on a Nikon eclipse Ti-U inverted microscope with an Oil Nikon Apo TIRF 100× objective and recorded with a ORCA Flash 4.0 v2 (digital CMOS C11440-22CU) Hamamatsu camera, controlled by μManager. For each protein droplet, 200 frames we acquired with 10 ms exposures at 300 ms intervals.
Images were analyzed in MATLAB. Particle centers were found by fitting to a Gaussian distribution. Particle locations were linked into trajectories using.47 Diffusion coefficents were determined by fitting the mean-squared displacement versus time. The viscosity was determined from the diffusion coefficient using the Stokes–Einstein equation.
D = kBT/6πRη, | (6) |
Footnote |
† Electronic supplementary information (ESI) available. See DOI: 10.1039/d0sm01319f |
This journal is © The Royal Society of Chemistry 2021 |