Diffusion-based determination of protein homodimerization on reconstituted membrane surfaces.

The transient interactions between cellular components, particularly on membrane surfaces, are critical in the proper function of many biochemical reactions. For example, many signaling pathways involve dimerization, oligomerization, or other types of clustering of signaling proteins as a key step in the signaling cascade. However, it is often experimentally challenging to directly observe and characterize the molecular mechanisms such interactions-the greatest difficulty lies in the fact that living cells have an unknown number of background processes that may or may not participate in the molecular process of interest, and as a consequence, it is usually impossible to definitively correlate an observation to a well-defined cellular mechanism. One of the experimental methods that can quantitatively capture these interactions is through membrane reconstitution, whereby a lipid bilayer is fabricated to mimic the membrane environment, and the biological components of interest are systematically introduced, without unknown background processes. This configuration allows the extensive use of fluorescence techniques, particularly fluorescence fluctuation spectroscopy and single-molecule fluorescence microscopy. In this review, we describe how the equilibrium diffusion of two proteins, K-Ras4B and the PH domain of Bruton's tyrosine kinase (Btk), on fluid lipid membranes can be used to determine the kinetics of homodimerization reactions. [BMB Reports 2021; 54(3): 157-163].

INTRODUCTION

The cellular membrane is the primary theater of signal transduction, which begins at the receptor binding at the membrane surface (1, 2). Transient protein-protein and proteinlipid interactions on membrane surfaces constitute the key steps in these signaling cascade reactions (2-5). These interactions are often combinatorial in design, involving reactions such as homo- and hetero-oligomerization (6-10). However, the molecular details of many such systems still remain mysterious. This is largely because mechanistic investigations, which must constitute quantitative measurements, on self-assembled lipid bilayer structures can be a prodigious technical challenge. On experimental systems utilizing live cells, there are added complications due to the copious number of additional background processes which may or may not participate the process under observation. Fortunately, many advances in membrane reconstitution as well as imaging and spectroscopic techniques have made it possible to probe quantitative information such as kinetic rate constants of membrane association and dissociation, catalytic rates, and diffusion constants (11-14). In particular, supported lipid bilayers (SLBs) stand out in its straightforwardness in preparation, experimental implementations, and data analysis (15-21). SLBs are planar lipid bilayers that can be as large as square centimeters in area. They are typically prepared by the rupturing small unilamellar vesicles, which then fuse with each other into a single bilayer on a solid support such as glass coverslips (16, 22).

One way to detect protein-protein interactions on membrane surfaces is by measuring diffusion, which is sensitive to the size of the proteins as well as their conformations, orientations, and overall interaction with the membrane. In addition to the obvious advantages of SLBs in allowing precise control over their chemical environment, another important advantage is that the SLB lipids and other molecules embedded within are fluid and exhibit unencumbered Brownian diffusion (17, 23, 24). This is in stark contrast with the membrane diffusion in live cells where diffusion is usually anomalous (24). The well-defined diffusion behavior allows the diffusion measurements performed on SLBs to be interpreted without much ambiguity. Thus, SLBs have been demonstrated in its usefulness in many works examining molecular interactions on membrane surfaces. Two such examples are K-Ras4B GTPase (25) and the pleckstrin homology (PH) domain of Bruton’s tyrosine kinase (Btk) in their dimerization behavior (26). Homodimerization is perhaps the simplest possible protein-protein interaction and is ubiquitous in signal transduction. Nevertheless, the accurate determination of its kinetics is still a significant technical challenge, particularly if it is dependent on membranes. This is the case for both K-Ras4B and Btk, whose dimerization has never been observed in solution even though they were thought to be critical in their function. In this review, the diffusion-based investigation of their dimerization behavior of K-Ras4B and Btk using fluorescence spectroscopy and imaging will be described.

TWO-DIMENSIONAL DIFFUSION ON MEMBRANES

While it can be reasonably assumed that the diffusion coefficient of a molecule tethered to a membrane would be dependent on its overall size, their theoretical relationship is not necessarily straightforward. In a three-dimensional space, particle diffusion follows Stokes-Einstein relation, where the diffusion coefficient D is inversely proportional to the hydrodynamic radius r of the particle,

Where kB is Boltzmann’s constant, T is the absolute temperature, f is the friction, and η is the viscosity of the medium. In socalled Stokes’ Paradox, however, no simple relationship between the particle size and the diffusion coefficient in two dimensions. In specific case of a cylinder embedded in a viscous two-dimensional medium (such as a lipid bilayer), its diffusion is described by Saffman-Delbruck model, in which diffusion coefficient is proportional to the natural log of inverse radius of the cylinder (27). Experiments have mixed results with respect to this model’s ability to accurately predict the diffusion from the protein size (28, 29). Nevertheless, the model’s assumptions are incompatible with many membrane proteins if they deviate from the cylindrical shape with heights mismatching that of the membrane. Thus, the applicability of Saffman-Delbruck model is limited in cases such as peripheral membrane proteins that are only associated with one leaflet of the membrane to varying degrees.

For peripheral membrane proteins, many factors influence the observed diffusion: [1] the number of lipid anchors attached to the protein, if any, [2] the degree of immersion and interaction (including nonspecific ones such as Coulomb interactions) of the protein with the membrane, [3] the structure and orientation of the protein with respect to the membrane, among others (30). Without a detailed knowledge of these factors, it is in principle impossible to predict the diffusion behavior of the peripheral membrane proteins. A useful alternative is to experimentally calibrate the diffusion coefficients as a function of the molecular size of the proteins. In an in vitro study by Knight et al. (31), multimers of pleckstrin homology (PH) domains of general receptor of phosphoinositides 1 (Grp1), which selectively binds to phosphatidylinositol (3, 4, 5)-trisphosphate (PIP3) lipids were engineered, and their diffusion on supported lipid bilayers were measured as a function of the lipid binding stoichiometry. It was found that diffusion coefficient is inversely proportional to the number of lipid anchors–in other words, the PH multimers followed Stokes-Einstein relation (Equation 2):

Where fN is the additive friction introduced by the Nth lipid anchor.

DIFFUSION MEASUREMENTS BY FLUORESCENCE

Fluorescence spectroscopy and microscopy are by far the most common methods to measure diffusion on biological membranes. Among various techniques, fluorescence correlation spectroscopy (FCS) and single-molecule tracking (SMT) with surface-selective total internal reflection microscopy (TIRFM) offer high precision and reliability and have been used in a number of studies seeking to establish molecular interactions (Fig. 1), including the K-Ras 4B and Btk PH domain studies.

Fig. 1

(A) In fluorescence correlation spectroscopy, a confocal laser illumination spot is used to monitor the intensity fluctuation due to fluorescence molecules diffusing in and out of the focus. The auto-correlation function (ACF) of the resulting time trace is used to compute the particle number (N, from the y-intercept) and residence time (τd, from the half-point of the ACF decay) of the fluorescent molecule in the laser focus, which is related to the diffusion coefficient by w2 = 4Dτd. In this example, a GFP-labeled protein on a supported lipid bilayer is excited by a blue (488 nm) laser. (B) In total internal relfection fluorescence microscopy (TIRFM) single molecule tracking, the displacement of the particle between each frame are reconstructed into particle trajectories (example trajectory shown in inset). The step size distribution may be used to calculate the diffusion coefficient and the relative population of multiple diffusing species, following the step size equation P(r) (Equations 5 and 6 for single and two species, respectively). The fit to a single species step size distribution is shown. (C) Left: theoretical binding curves for a dimerization reaction for a wide range of dissociation constants, right: simulated FCS data using experimentally determined diffusion coefficients for monomers and dimers for the same dissociation constants.

Fluorescence correlation spectroscopy

FCS is a type of fluorescence fluctuation spectroscopy (FFS), which are techniques based on the time-dependent fluorescence intensity fluctuations. It was invented in more than thirty years ago (32), but it has not been very useful until the widespread use of very stable lasers in more recent years. FCS is now commonly used to characterize biological molecules in solutions, membranes, and in live cells (33, 34). It relies on observing fluorescent molecules in a laser focus, which is typically close to diffraction-limited and very small focus volume (∼10−15 L). Therefore, the number of molecules observed is very small— which is why sometimes they are considered “single-molecule” techniques (35). The time-dependent intensity fluctuations are due to fluorescent molecules diffusing in and out of the focus, and therefore contain information about the average number of the molecules within the focus volume and their diffusivity (Fig. 1A). To extract these information, the autocorrelation function (ACF) of the fluorescence intensity fluctuation trace is computed, and it is fit to a autocorrelation model with Gaussian focus G(τ). For a two-dimensional Brownian diffusion, the model is given by (Equation 3)

Where N is the number of the particles in the focus, D is the diffusion coefficient, τ is the average focus residence time of the molecule, and w is the radius of the Gaussian focus. Therefore, the absolute values of concentration and diffusion coefficient may be obtained with calibrated focus radius w.

An added benefit in using FCS to detect multimerization of fluorescent molecules is that it is particularly sensitive to brighter species. If there are mixed species present in FCS, the autocorrelation model is given by (Equation 4)

Where Bi is the molecular brightness of the ith species and (25). Note that the overall autocorrelation is dependent on the square of Bi, so brighter species will contribute non-linearly to the overall autocorrelation, and further outweighs the contribution by N. For example, if the sample is an equimolar mixture of a monomer with brightness with B and dimer of brightness 2B, the detected autocorrelation will be mostly due to the dimer even though the concentrations of monomers and dimers are the same. This “exaggerated” presence of brighter species is particularly advantageous in detecting and measuring the multimeric species where the multimerization affinities are very weak and there is only a small fraction of dimers in a mostly monomeric system.

Single-molecule tracking

In SMT, the motions of single fluorophores are directly visualized and captured by a high quantum efficiency camera (4, 23, 30, 36-38). By using extrinsic fluorophores of high quantum yield and photostability, long-lived, high signal-to-noise ratio trajectories of single molecules can be recorded (Fig. 1B). In the data analysis, the molecule position is identified, and by linking them the particles between teach time frame, a trajectory for each molecule is constructed (39). Thousands of such trajectories can be collected in a single movie, and the diffusion coefficient D of the particles can be obtained by analyzing the step size distribution P(r) given by (Equation 5)

Where r is the particle displacement and t is the time interval (11, 23, 40). If there are two distinct diffusing populations with diffusion coefficients D1 and D2, the distribution is (Equation 6)

Where α is the relative population of the diffusing species with diffusion coefficient D1. Thus, provided that the camera has a sufficient frame rate to resolve the step sizes associated with diffusion coefficients for monomeric and dimeric species, SMT is capable of resolving the diffusion coefficients and relative populations of monomers and dimers. A practical consideration in the implementation of SMT experiments the surface density of fluorophores must be kept low as to spatially resolve the diffusing species. Stoichiometric under-labeling of fluorophores (for example, if one in 106-109 molecules are labeled with a fluorophore on a bilayer with a surface density of 1,000 molecules/µm2, there are < 100 fluorescent molecules per field of view in a typical EMCCD camera) is frequently used to achieve this condition (36, 41).

DIMERIZATION OF RAS

K-Ras4B is one of four isoforms of Ras, which is a central molecular switch in a number of signaling pathways such as the mitogen-activated protein kinase (MAPK) and phosphoinositide 3-kinases (PI3K)-AKT pathways (42, 43). A small globular protein that is tethered to the membrane by a lipid anchor, Ras—and especially K-Ras4B—is of a paramount clinical importance in that its oncogenic mutations, where it is unable to return to the GDP-bound “off” state, are found in ∼30% of all human cancers (44). Despite concerted efforts to develop therapeutic strategies to inhibit Ras, all conventional, small molecule-based competitive inhibitors to date have failed (45). This has led to an interest in novel strategies to target Ras by modulating its spatial organization, as Ras is observed in dimers and clusters on cellular surfaces (46-49). In particular, it was speculated that its dimerization drove the assembly—this was supported by the consistent dimeric structures in crystals, as well as the fact that some interacting partners of Ras do so in pairs (for example, SOS inter-acts with two Ras molecules) (50-52). However, Ras dimers are never observed in solution, and despite many attempts to detect them on membranes no consensus exists regarding the nature of Ras dimers—whether they are intrinsic properties of Ras, or if so, which specific interactions are responsible—and even whether they exist at all (41, 53-56).

Whether natively expressed and fully processed K-Ras4B can homodimerize on its own on membrane surfaces without extrinsic factors was investigated by using SLB-based reconstitution, FCS, and SMT (25, 57). Two fluorescent labeling strategies were used: eGFP fused to K-Ras4B and an extrinsically labeled guanosine nucleotides that binds to K-Ras4B, and both displayed identical results. A definitive answer to this question would involve a measurement of the two-dimensional dissociation constant, , of the dimerization reaction, (Equation 7)

Where and are equilibrium surface densities of Ras monomers and dimers, respectively. The basic strategy to determine the dissociation constant was to measure the change in diffusion, induced by dimerization, as a function of surface density: in other words, obtain a titration curve where the readout for the surface density of dimers was the decrease in average diffusion. The theoretical surface density titration is shown in Fig. 1C, left panel. A number of titration curves for a range of Kd are plot-ted; a dimerization reaction with Kd of 200 molecules/µm2 indicates a fairly strong dimerization affinity, whereas Kd of 1.25 × 104 molecules/µm2 would represent a dimerization affinity so weak that it is only significant at a protein surface density approaching the close-packed state.

This strategy requires the knowledge of protein surface density and the diffusion coefficients of the monomers and dimers. While the surface density and diffusion coefficient of monomers may be readily measured by FCS or SMT, it is less straightforward to measure the diffusion coefficient of dimers of unknown dissociation constant, especially if it is weak. This problem was addressed by deliberately inducing dimerization using a strong cross-linker, in this case the Ras binding domain (RBD) of c-Raf fused to leucine zipper (LeuZ) motif (RBD-LeuZ). RBD selectively binds to GTP-bound Ras, and LeuZ is a constitutively dimeric motif. Therefore, the introduction of LeuZ to GTP-bound Ras results in crosslinking of Ras via LeuZ, and the diffusion coefficient of the dimeric Ras could be determined. In this case, the diffusion coefficient a dimer was approximately half that of a monomer, akin to the system constructed by PH domain multimers (31).

Because it was expected that Ras would have a very weak dimerization affinity, it was important to establish the detection limit of the measurement. Using the experimentally determined values of the diffusion coefficients of monomers and dimers, theoretical FCS titration measurements for a range of dissociation constants were simulated (Fig. 1C, right). The simulation shows that, due to the nonlinear contribution of brighter species to the overall measurement, FCS provides a wide dynamic range from strong to very weak dimerization reactions. The brightness effect is so prominent that even very weak dimerization reactions with Kd is in the surface density range of close-packed proteins would be detected by FCS within the experimentally accessible surface density (∼1,500 molecules/µm2).

The diffusion measurements indicated that, by itself, K-Ras4B remained unambiguously monomeric. Fig. 2A shows representative FCS and SMT step-size distribution data. In FCS (Fig. 2A, left), GTP-bound K-Ras4B does not change significantly with respect to increasing surface density, unlike K-Ras4B crosslinked by RBD-LeuZ. The estimated Kd for GTP-bound K-Ras was 1.25 × 105 µm2, which is close to the close-packed density (58) for the protein—indicating effectively no dimerization. The same was observed with SMT step size distributions: they are independent of the overall surface density, unless there is the crosslinker present (Fig. 2A, right). The same was observed in all conditions tested, such as different membrane lipid compositions including cholesterol, activation states, and chemical environments. These results suggested that higher order multimerization of Ras observed on cells are not driven by its intrinsic dimerization capability, but rather by other factors, such as scaffolding factors or downstream proteins, may be involved.

Fig. 2

(A) Left: FCS data for the dimerization of GTP-bound K-Ras4B (GTP) and GTP-bound K-Ras4B with the RBD-LeuZ crosslinker (GTP + crosslinker) on supported lipid bilayers, and the estimated dissociation constants for each. The lack of change in diffusion as a function of surface density without the crosslinker suggests that K-Ras4B does not dimerize on its own. Right: step size distributions of K-Ras4B shows similar results: without the crosslinker, there is no change in diffusion as the surface density is increased. (B) FCS measurements for the wild type Btk PH-TH domains (left) and the mutant lacking key dimerization residues (right). Substantial dimerization is only observed for the wild type at a high PIP3 surface density.

DIMERIZATION OF BRUTON’S TYROSINE KINASE (Btk) PLECKSTRIN HOMOLOGY (PH) DOMAIN

Bruton’s tyrosine kinase (Btk) is a protein kinase that is critical in activation of B-cells (59, 60). Conventionally, its activation depends on the recruitment to the plasma membrane by interaction with phosphatidylinositol (3, 4, 5)-trisphosphate (PIP3) lipids, which are produced upon receptor stimulation. The pleckstrin homology domain fused with a tec homology domain (PH-TH) of Btk selectively binds two PIP3 headgroups, and results in activation of Btk (61-64). Crystal structures show a consistent dimer structures of the PH-TH domain, and evidence suggested that trans-autophosphorylation (a Btk molecule phosphorylates the other in a dimer) mediated by dimerization may be involved in activation of Btk (65, 66). However, as it was the case for Ras, PH-TH domain dimers were never observed in solution. This led to a speculation that the interaction with the lipid ligand, PIP3, in the membrane may be required to conformational changes that allow dimerization.

The membrane-dependent dimerization was studied with FCS on supported lipid bilayers. The eGFP-fused PH-TH domain of Btk (an alternative, cysteine-labeled Alexa Fluor 647 also resulted in identical behavior) was introduced to SLBs containing a trace amount (0.005%) of fluorescently labeled lipids (Texas Red-DHPE) and 1-4% PIP3, and which recruited the protein to the membrane. The diffusion of both eGFP-Btk PH-TH domain and Texas Red-DHPE were simultaneously measured by two-color FCS. The titration curve was obtained by plotting the relative diffusion (PH-TH diffusion divided by lipid diffusion) as a function of the protein surface density. This had an effect of correcting for variations in the measured diffusion arising from microscopic heterogeneity in the membrane environment.

The left panel of Fig. 2B shows the FCS titration measurements for the wild type Btk PH-TH domain. There is a robust decrease in the overall diffusion as the surface density is increased on a 4% PIP3 membrane, consistent with dimerization. It was also found that when the PIP3 density lower at 1%, there is a substantially less change in diffusion, indicating smaller fraction of dimers at the same protein surface density. To confirm that this observed change in diffusion is due to dimerization, the same experiment was performed with a mutant construct whose putative dimerization interactions (65) had been removed. The resulting titration curve (Fig. 2B, right panel) shows that the diffusion is no longer modulated by the increased surface density, indicating that dimerization predicted by the crystal structure is responsible for the observed changes in diffusion. This observation, combined with adsorption kinetics measurements by total internal fluorescence microscopy, led to a Btk activation model in which each Btk PH-TH domain interacts with two PIP3 lipids then dimerizes and activates through trans-autophosphorylation. It was also suggested that such an elegant combinatorial design in protein-protein and proteinlipid interactions may grant switch-like PIP3 sensitivity to the Btk activation scheme, which is a common theme in cellular signal transduction (6, 67-69).

CONCLUSIONS

Dimerization and higher order multimerization of membrane proteins are a recurring motif in the cellular signal transduction, and a possible protein-protein interaction target for therapeutic interventions. However, it has been challenging to observe complexes and assemblies with traditional in vitro strategies, as they fall apart when as soon as membranes are solubilized. They are better preserved and observed on cells, but the platform suffers a different problem that it is difficult to systematically quantify molecular interactions on living cells. Reconstitution strategies such as those using supported lipid bilayers have had shown a great promise, particularly when combined with quantitative fluorescence microscopy such as FCS and SMT. Using diffusion as a readout for molecular interactions allow robust and precise measurement of the dissociation constant. In the case of K-Ras, it was unambiguously shown that it lacks intrinsic ability to dimerize on membranes, suggesting that the clusters observed on cells are assembled by other extrinsic factors. For Btk PH-TH domain, dimerization was never observed in solution even though it was suspected that it was critical in its activation. Its diffusion measurements revealed that the dimerization contact predicted by the crystal structures are respon-sible for membrane-mediated dimerization. Combined with membrane adsorption measurements, the entirety of the protein’s interaction kinetics were elucidated. While there are limitations to this method in that it can only be applied to reconstituted systems in homogeneous lipid bilayers, these results de-monstrate that diffusion-based investigation is a promising strategy for the mechanistic characterization of molecular interactions on membranes. With the continuing advances in microscopy, data analysis algorithms, and membrane fabrication technology, the utility of this method may be expanded to live cell contexts and transmembrane receptors, opening an opportunity to unprecedented understanding the cellular signaling phenomena.