Multiplexed, high-throughput measurements of cell contraction and endothelial barrier function.

Vascular leakage, protein exudation, and edema formation are events commonly triggered by inflammation and facilitated by gaps that form between adjacent endothelial cells (ECs) of the vasculature. In such paracellular gap formation, the role of EC contraction is widely implicated, and even therapeutically targeted. However, related measurement approaches remain slow, tedious, and complex to perform. Here, we have developed a multiplexed, high-throughput screen to simultaneously quantify paracellular gaps, EC contractile forces, and to visualize F-actin stress fibers, and VE-cadherin. As proof-of-principle, we examined barrier-protective mechanisms of the Rho-associated kinase inhibitor, Y-27632, and the canonical agonist of the Tie2 receptor, Angiopoietin-1 (Angpt-1). Y-27632 reduced EC contraction and actin stress fiber formation, whereas Angpt-1 did not. Yet both agents reduced thrombin-, LPS-, and TNFα-induced paracellular gap formation. This unexpected result suggests that Angpt-1 can achieve barrier defense without reducing EC contraction, a mechanism that has not been previously described. This insight was enabled by the multiplex nature of the force-based platform. The high-throughput format we describe should accelerate both mechanistic studies and the screening of pharmacological modulators of endothelial barrier function.

Introduction:

Across several major diseases including acute lung injury, atherosclerosis, tumor neovascularization, hypertension, diabetes mellitus, and coronary artery disease, the endothelial barrier becomes disrupted through the formation of paracellular gaps [1-8]. Since the pioneering work of Guido, Majno, and George Palade in the 1960s [9, 10], paracellular gaps in response to inflammatory mediators have been considered to arise by contractile forces that act to pull intercellular junctions apart [11-22]. These forces are generated via acto-myosin interactions and manifest upon the EC substrate as tractions [12, 18–21, 23, 24].

Endothelial barrier function is routinely measured through changes in electrical resistance or macromolecular permeability across a confluent monolayer [25]. While these measurements are fast, they cannot simultaneously quantify underlying biophysical forces. As a result, further validation requires follow-up measurements of EC contraction using single cells [21], cell doublets [20, 23], micropatterned cell clusters [12, 18], and continuous cell monolayers [19, 24]. These follow-up measurements are slow, complex, and low-throughput, limiting the pace of mechanistic insights and the breadth of pharmacological investigation.

Here, we have developed “mini-XPerT”, a force-based high-throughput assay that combines monolayer traction microscopy [19, 26] to measure EC contraction, Xpress Permeability Testing (XPerT) [27] to measure paracellular gaps, and immunofluorescence imaging to measure structural and morphological changes. These measurements are miniaturized in 96-well plates using the approach of contractile force screening [28, 29]. Mini-XPerT utilizes elastic silicone substrates that are tunable in Young’s modulus (0.3–150kPa) [29], thereby enabling studies spanning the entire known (patho)physiological range of substrate stiffness. Additional practical advantages of mini-XPerT include the following: it is non-invasive, non-destructive, and suitable for both short- and long-term assessment of barrier function.

Applying mini-XPerT, we have studied monolayers of primary human dermal microvascular ECs subjected to acute and progressive mediators of pathological permeability. We have compared barrier protective mechanisms of two well-known and mechanistically distinct barrier-protective agents, the Rho-kinase inhibitor Y-26732 [30] and the canonical Tie2 receptor agonist Angiopoietin-1 (Angpt-1) [31]. The multiplexed nature of mini-XPerT unveiled evidence to suggest that Angpt-1 achieves barrier protection without diminishing contractile forces. Taken together, these results reveal that the endothelium can protect itself by: 1) reorganizing the actin cytoskeleton – a contraction-dependent mechanism – as in the case of Y-27632 or 2) by accumulating VE-cadherin to fortify cell-cell junctions – a contraction-independent mechanism – as in the case of Angpt-1.

Material and Methods:

While our methodologies are generalizable across adherent cell types for a wide spectrum of experimental conditions, we focused here on a representative endothelial cell type, a commonly studied range of substrate stiffness values, and a commonly used extracellular matrix ligand, collagen I.

Cell culture:

All experiments were performed using primary dermal human microvascular endothelial cells (Lonza, Basel, Switzerland) at passages 4–7 in EGM2-MV medium (Lonza, Basel, Switzerland) containing 4.7% Fetal Bovine Serum. The cells were plated for 2–3 days prior to experimentation.

Fabrication of multi-well plates:

Soft, elastic (Young’s Modulus = 0.3 or 3kPa), and deformable substrates were prepared using NuSil® 8100 (NuSil Silicone Technologies, Carpinteria, CA) in 96-well plates [29]. Embedded in the substrate surface are fluorescent beads (diameter = ~600nm) whose displacements enable traction force calculations. Ligated to the substrate surface is biotinylated collagen I (collagen I, Advanced Biomatrix, Carlsbad, CA; 10x PBS, Corning Life Sciences, Tewksbury, MA; EZ-Link NHS-LC-LC-Biotin, Thermo Fisher Scientific, Waltham, MA) whose binding with FITC-avidin enables the identification of paracellular gaps. Detailed information of plate preparation, ligation, and cell seeding is provided in the supplemental material.

Measurement of cellular contractile forces – monolayer traction microscopy[32]:

The local contraction that an endothelial cell exerts upon their substrate is called traction [33]. To measure tractions, we utilized an inverted epi-fluorescence microscope (DMI 6000B, Leica Inc., Germany) equipped with a heated chamber (37°C), a monochrome camera (Leica DFC365 FX), and a motorized stage. We recorded spatial images of substrate-embedded fluorescent beads at 10x magnification. Based on the bead displacements (resolution = ~15μm) relative to a cell-free image, together with knowledge of substrate stiffness and thickness, we computed tractions using the approach of Fourier Transform Traction Cytometry [34], modified to the case of cell monolayers [32]. From each traction map, we calculated the root-mean squared traction (RMST) value and report this as a measure of average contraction in the monolayer. On a well-by-well basis, we first computed the ratio of the RMST after vs. before stimulation. We then normalized all values of a given treatment group to the mean of the corresponding vehicle control group and further calculated the average of each treatment group.

This permits us to compare differences across treatment conditions in each 96-well plate as well as across plates.

Measurement of paracellular gaps –Express permeability testing (XPerT)[27]:

At the end of drug treatment, FITC-conjugated avidin dissolved in PBS (Thermo Fisher Scientific, Waltham, MA) was added to each well and incubated for 210s. The cellular monolayer acts as a mask; Avidin binds to biotinylated collagen I only in areas of paracellular gaps [27]. When further visualized using an inverted epi-fluorescence microscope (DMI 6000B, Leica Inc., Germany, 10x magnification), the FITC label localizes the gap and thereby permits quantitation of both its location and overall incidence in a well. The raw FITC images were converted into 8-bit binary images; the same threshold value was applied to all images in the manuscript. Gaps were computed from the areas in white and expressed as the percentage of the total area in a given image. We then normalized all values of a given treatment group to the mean of the corresponding vehicle control group and further calculated the average of each treatment group. To perform these analyses in high-throughput, a custom-macro was written using Fiji software. The macro is provided in the supplemental material.

Quantification of F-actin stress fiber expression and orientation:

At the end of drug treatment, the cells were fixed with 3% Formalin (Sigma-Aldrich, St. Louis, MO) for 5 minutes and then washed three times with PBS. Alexa Fluor 480 phalloidin (Thermo Fisher Scientific, Waltham, MA) at a concentration of 1:300 and Hoechst dye (Thermo Fisher Scientific, Waltham, MA) at a concentration of 1:2000 was added together to each well for 30 minutes. The wells were washed twice with PBS and visualized using an inverted epi-fluorescence microscope (DMI 6000B, Leica Inc., Germany, 20x magnification). The Alexa 480 labels the actin stress fibers while the Hoechst dye labels the nuclei. From each image, we chose at random, six different regions. In these regions, we quantified F-actin orientation (in terms of anisotropy) using an open source plugin called Fibriltool[35] for Fiji. We also quantified corresponding F-actin expression (in terms of mean-fluorescent intensity) in those same regions. Across all images for a given treatment group, we normalized all values to the mean of the corresponding vehicle control group and further calculated the average of each treatment group.

Evaluation of VE-Cadherin expression:

At the end of drug treatment, the cells were fixed with 3% Formalin for 5 minutes, washed three times with PBS, and then permeabilized with 0.1% saponin (Sigma-Aldrich, St. Louis, MO) for 5 minutes. A mouse-derived monoclonal anti-VE-cadherin antibody (BD Biosciences, Franklin Lakes, NJ) at a concentration of 1:250 was added to each well for 2 hours. The wells were washed twice with PBS and a secondary antibody conjugated with FITC (Thermo Fisher Scientific, Waltham, MA) at a concentration of 1:500 and Hoechst dye (Thermo Fisher Scientific, Waltham, MA) at a concentration of 1:2000 was added together for an additional 30 minutes. When visualized using an inverted epi-fluorescence microscope (DMI 6000B, Leica Inc., Germany, 20x magnification), the FITC-conjugated antibody labels VE-Cadherin expression at the cell-cell junction and the Hoechst dye labels nuclei.

Measurement of trans-endothelial electrical resistance (TEER)[36]:

An electric cell-substrate impedance sensing system (ECIS, Applied BioPhysics Inc., Troy, NY) was used to measure trans-endothelial electrical resistance of EC monolayers. Briefly, cells were seeded on collagen I coated electrode arrays (8W10E+, Applied BioPhysics Inc., Troy, NY) and grown to confluency (24–48h), determined by achieving electric resistance of >1500 ohms. Cells were then incubated with indicated drugs or EBM-2 medium alone as a vehicle control with raw resistance continually measured over time at an amplitude of 4,000 Hz. For every well, raw resistance values were divided by the correspondent time-zero raw resistance value to obtain the normalized TEER.

Drug treatments:

EBM-2 basal medium served as vehicle for drug treatments, and as the vehicle control. Angiopoietin-1, CD14 and LPS binding protein were purchased from R&D Systems (Minneapolis, MN). Thrombin was purchased from EMD Millipore (Billerica, MA). LPS serotype O111:B4 was purchased from Sigma-Aldrich (St. Louis, MO). Y-27632 was purchased from EMD Biosciences (LaJolla, CA).

Statistics:

Data are presented as mean ± SEM. For gaps, ECIS measurements, F-actin intensity and anisotropy measurements, statistical significance was tested by two-tailed unpaired t test with Welch’s correction. For traction measurements, the more appropriate non-parametric two-tailed Mann-Whitney U Test for data deviating from Gaussian distribution was used. A p value less than 0.05 was considered significant.

Results and Discussion:

A key downstream effector of EC barrier disruption is the Rho-associated protein kinase, ROCK [37-39]. Activation of ROCK signaling promotes actin stress fiber formation, myosin light chain phosphorylation, endothelial cell contraction and, ultimately, paracellular gap formation [15, 30, 37–43]. Conversely, inhibition of ROCK signaling via the pharmacological agent Y-27632 promotes endothelial barrier defense [12, 30].

We chose Y-27632 as an exemplar to evaluate mini-XPerT (Fig. 1). Based on typically used concentrations [12, 43–47], and supportive dose-response measurements (Fig. S1), we picked a sub-maximal dose of 5μM for further evaluation. Individual wells of a single 96-well plate (Young’s Modulus = 3kPa) were pre-treated with either vehicle or Y-27632 for 20 min, and then stimulated further with the barrier disruptive agent, thrombin (1U/ml, 30min). Paracellular gaps and monolayer contraction were significantly enhanced by thrombin and diminished by pre-treatment with Y-27632 (Fig. 1B-E). Moreover, these force changes correlated with F-actin cytoskeletal changes, as was revealed by simultaneous measurements of the F-actin cytoskeleton (Fig. S2). Specifically, while thrombin enhanced F-actin expression and orientation, these same quantities were significantly reduced by pre-treatment with Y-27632. Finally, the barrier-protective effects of Y-27632 was consistent with ECIS measurements (Fig. S3).

Figure 1:

Mini-XPerT – simultaneous measurements of barrier function and cellular contractile forces in a high-throughput format.

(A) Silicone-based elastic substrates (Young’s Modulus = 0.3kPa or 3kPa) were prepared in a 96-well plate format. Embedded in the substrate are fluorescent beads (diameter = ~400nm; shown in red) whose displacement enables the computation of monolayer contractile forces (shown as white arrows). Ligated to the substrate is biotinylated collagen I (shown as black triangles) whose binding with FITC-avidin (shown in green) enables the identification of paracellular gaps. (B) Representative gap images for every experimental condition, (C) corresponding traction force maps, and (D-E) averaged values. While gap formation and contraction were increased with thrombin (1U/ml, 30min) it was significantly reduced by pre-treatment for 20 min with the ROCK inhibitor, Y-27632 (5μM). In (D), plotted is the gap area normalized to the vehicle group. In (E), on a well-by-well basis, we quantified the ratio of the average monolayer contraction after vs. before stimulation. We further normalized this ratio to the vehicle group. Each treatment group comprises of measurements performed over n=7–16 individual wells of a 96-well plate. Plotted is the mean and standard error. * indicates p<0.05 while ns indicates p>0.05.

Next, we examined biophysical mechanisms of barrier protection for the Tie2 receptor agonist, Angpt-1. In cells, tissue, organs, and animals, Angpt-1 has been demonstrated to counteract permeability and weakening of the endothelial barrier [31]. While roles for actin stress fiber rearrangements, junctional remodeling, and mechanotransduction have been implicated [48-53], in the absence of direct mechanical measurements, the underlying biophysical mechanisms remain unclear. Based on typically used concentrations [54–5936, 49, 60], and additional dose-response measurements (Fig. S4), we picked a sub-maximal dose of 300ng/ml for further evaluation. Individual wells of a single 96-well plate (Young’s Modulus = 3kPa) were co-stimulated with either vehicle or Angpt-1 together with the barrier disruptive agent, thrombin (1U/ml; 30min) (Fig. 2). In accordance with previous findings [36], Angpt-1 co-stimulation reduced thrombin-induced gap formation (Fig. 2A). While this barrier protective effect, including time-course and efficacy, was recapitulated by ECIS (Fig. 2C-D), it was not correlated with a reduction in cell contractile force (Fig. 2B) or in the orientation and expression of the F-actin cytoskeleton (Fig. S2). However, it was associated with an increase in VE-Cadherin expression at the cell-cell junctions (Fig. S5).

Figure 2:

Angpt-1 reduces thrombin-induced gap formation but does not reduce cellular contractile force.

ECs in a 96-well plate were incubated with either vehicle, vehicle + thrombin (1U/ml), Angpt-1 (300ng/ml), or Angpt-1 + thrombin for 30 min. (A) While paracellular gaps were reduced by Angpt-1 co-treatment (D) cellular contractile forces were not. For gaps, we report changes normalized to vehicle. For contraction, quantities have been normalized to t=0min values and further to the vehicle group. Each treatment group comprises of measurements performed over n=7–16 individual wells of a 96-well plate. In all bar plots, shown is the mean and standard error. * indicates p<0.05. (E) ECIS measurements confirm mini-XPerT findings. Each curve is averaged over n=3–5 individual wells. The grey line marks the time point of 30 minutes post treatment; this time point corresponds to gap measurements with mini-XPerT. (F) Shown are corresponding resistance values. Plotted is the mean and standard error. * indicates p<0.05 while ns indicates p>0.05.

To demonstrate the versatility of mini-XPerT, we examined other physiologically-relevant endothelial leak provocateurs that are commonly utilized in barrier function studies: 1) Gram-negative endotoxin lipopolysaccharides (LPS; 100ng/ml, 4h) [49, 61] and 2) the canonical host cytokine, tumor necrosis factor alpha (TNFα; 20ng/ml, 4h) [49, 62]. Notably, their barrier disruptive activity manifests after several hours (e.g. [49, 63]), as opposed to thrombin, whose impact is acute, in the order of minutes (e.g. Fig. 2C). We evaluated these triggers of permeability on soft (0.3kPa) and stiff (3kPa) substrates with and without Y-27632 pre-treatment or Angpt-1 co-stimulation (Fig. 3). Paracellular gaps were significantly enhanced by TNFα/LPS stimulation and were largely reduced by the additional presence of Y-27632 or Angpt-1. Correspondingly, Y-27632 reduced cell contractile forces (Fig. S6) while Angpt-1 did not (Fig. S7). Taken together, these findings reveal the generality of the barrier protective mechanism of Y-27632 and Angpt-1.

Figure 3:

Barrier protection by Y-27632 is correlated with a reduction of cellular contractile forces whereas barrier protection by Angpt-1 is not.

Plotted is gaps versus contraction for the indicated conditions. For gaps, we report changes normalized to vehicle. For contraction, quantities have been normalized to t=0min values and then to the vehicle group. Each treatment group comprises of measurements performed over n=7–16 individual wells of a 96-well plate. Plotted is the mean and standard error. Symbols in red, blue and green correspond to the LPS/TNFα/thrombin, Y-27632 and Angpt-1 treatments, respectively. The ellipses highlight contractile force dependent (blue) or independent (green) reduction of thrombin-, LPS, and TNFα-induced gap formation. A lack of concordance for the EC response to TNFα on the soft substrate remains unclear (Fig. S8).

If not through contractile force reduction, how then does Angpt-1 promote barrier defense? This question is of significant interest. Given that Angpt-1 has been shown to counteract endothelial barrier disruption triggered by numerous, unrelated permeability mediators ranging from VEGF to anthrax lethal toxin to thrombin [36, 49], the most widely invoked explanation is that Angpt-1 signaling modulates one or more common downstream effectors of barrier function [31]. We evaluated the two major candidates for this conserved action: 1) F-actin cytoskeletal rearrangement, and 2) VE-Cadherin accumulation at the cell-cell junction. We discovered a dominant role for the latter (Fig. S5), consistent with the postulated role for Tie2 signaling for VE-cadherin stabilization in vivo [49, 64]. Unexpectedly, Angpt-1 did not prevent agonist-induced enhancement of F-actin (Fig. S2) or contractile force (Fig. 2, Fig. S2).

In vascular biology, the dogma for EC barrier function has long held that junctional mechanisms are typically coupled to contractile-force dependent F-actin rearrangements. The present results, in contrast, are the first to our knowledge, that these two cellular processes as they relate to Angpt-1 induced barrier protection. If actin-cytoskeletal rearrangements can be uncoupled from junctional pro-barrier mechanisms in response to physiological stimuli, this could have implications for other cell types that exert barrier function such as squamous epithelium in the skin, columnar epithelium in the GI tract, and cuboidal epithelium in ducts and acinar structures of secretory organs.

Our data does not rule out the possibility that Angpt-Tie2 signaling is reorganizing forces on a much smaller scale, for example, in the peri-junctional region where the Rho family GTPase Rac1 has been shown to act in ECs [65]. Although beyond the resolution of the current technique, how such localized force reorganization could impact VE-cadherin accumulation at the cell-cell junction remains unclear.

We propose three important future directions for mini-XPerT. First, by combining mini-XPerT with monolayer stress microscopy [19], we can also determine intercellular stresses and its consequent effects on gap formation. This will require improved spatial resolution for traction measurements as well as clearly defined monolayer boundaries that can be specified through micropatterning procedures [12, 18, 19]. Second, by imaging gap formation in real-time, we can even more precisely correlate force and gap dynamics. To this end, we are currently evaluating cellular based markers to be able to delineate gaps in real time. Finally, it would be ideal to combine mini-XPerT with concomitant measurements of molecular signaling.