Paul
Pullara‡
,
Ibraheem
Alshareedah‡
and
Priya R.
Banerjee
*
Department of Physics, University at Buffalo, Buffalo, NY 14260, USA. E-mail: prbanerj@buffalo.edu
First published on 23rd December 2021
Liquid–liquid phase separation (LLPS) of multivalent biopolymers is a ubiquitous process in biological systems and is of importance in bio-mimetic soft matter design. The phase behavior of biomolecules, such as proteins and nucleic acids, is typically encoded by the primary chain sequence and regulated by solvent properties. One of the most important physical modulators of LLPS is temperature. Solutions of proteins and/or nucleic acids have been shown to undergo liquid–liquid phase separation either upon cooling (with an upper critical solution temperature, UCST) or upon heating (with a lower critical solution temperature, LCST). However, many theoretical frameworks suggest the possibility of more complex temperature-dependent phase behaviors, such as an hourglass or a closed-loop phase diagram with concurrent UCST and LCST transitions. Here, we report that RNA–polyamine mixtures undergo a reentrant phase separation with temperature. Specifically, at low temperatures, RNA–polyamine mixtures form a homogenous phase. Increasing the temperature leads to the formation of RNA–polyamine condensates. A further increase in temperature leads to the dissolution of condensates, rendering a reentrant homogenous phase. This dual-response phase separation of RNA is not unique to polyamines but also observed with short cationic peptides. The immiscibility gap is controlled by the charge of the polycation, salt concentration, and mixture composition. Based on the existing theories of complex coacervation, our results point to a complex interplay between desolvation entropy, ion-pairing, and electrostatic interactions in dictating the closed-loop phase behavior of RNA–polycation mixtures.
Earlier theories such as those by Flory and Huggins focused on enthalpy-driven LLPS.11 Subsequent efforts pointed to the existence of other factors such as solvation energies that lead to an entropy-driven LLPS.12–16 These two broad categories of LLPS have an orthogonal dependence on temperature. Enthalpy-driven LLPS is often manifested in an upper critical solution temperature (UCST) phase behavior, where increasing temperature leads to the suppression of LLPS.17 In contrast, entropy-driven LLPS often has a lower critical solution temperature (LCST) phase behavior, where LLPS is promoted at high temperatures and suppressed at low temperatures.17 It is worth noting that Flory and Huggins developed their theory for neutral polymer systems. For charged polymers, however, the existence of long-range electrostatic interactions, specific ion-mediated interactions, and solvation effects can give rise to more complex behaviors that are observed both in biology and in synthetic polymer systems. Therefore, mapping the temperature-dependent phase behavior of biopolymeric systems can give significant insights into the nature of thermodynamic driving forces underlying LLPS.
Many natural disordered proteins undergo liquid–liquid phase separation with a UCST phase behavior. For example, Molliex and coworkers reported that the ribonucleoprotein hnRNPA1 has a UCST phase behavior.18 Other proteins such as DDX4 and TIA1 were also shown to undergo UCST phase transitions.17 The roles of charged, polar, and aromatic residues in dictating the UCST phase behavior of proteins have been extensively discussed.17,19–21 On the other hand, there are few protein systems found in the literature that can undergo an LCST phase transition. Amongst naturally occurring proteins, microtubule-associated protein tau has been shown to undergo either LCST or UCST transition depending on the salt conditions.22–24 Another protein that has been shown to undergo an LCST transition is the proteasomal shuttle factor UBQLN2.25 Among synthetic systems, a seminal work from Quiroz and Chilkoti showed that elastin-like polypeptides (ELPs) and resilin-like polypeptides (RPLs) can be programmed to undergo either an LCST or a UCST phase transition, respectively, depending on their primary sequence.26 The authors identified distinct classes of amino acids that increase the tendency for LCST (charge depleted and hydrophobic amino acids) and UCST (polar and aromatic amino acids) transitions. Creating hybrid sequences with blocks of different amino acid compositions can lead to the coexistence of UCST and LCST in the same system, giving rise to an hourglass phase diagram with a miscibility gap at intermediate temperatures.27,28 Further, the condensation and aggregation of synthetic systems such as charged polyelectrolytes have also been extensively studied against temperature and salt with UCST behavior being frequently observed.29–33 However, a closed-loop phase diagram in the temperature-composition plane with TLCST < TUCST has been a rare trait for biologically relevant systems, although it is commonly discussed in theoretical classifications of phase diagrams and computational studies of polypeptides.17,34
In this letter, we report that a homopolymeric RNA, poly(rU), displays a temperature-controlled dual-response phase behavior with concurrent UCST and LCST transitions in presence of small polycationic molecules such as polyamines and peptides. The LCST transition occurs at temperatures below the UCST transition, leading to a reentrant phase transition with an immiscibility gap that is close to the physiological temperature. Combined with the fact that RNA–polycation condensates undergo a reentrant phase separation as the mixture composition is varied, we reveal a closed-loop phase diagram of RNA–polycation mixtures in the temperature-composition plane. The immiscibility gap is controlled by the valence of the cation, the mixture composition, and the salt concentration. Our experiments indicate that the LCST transition is not strongly affected by the ionic strength of the mixture or the composition asymmetry while the UCST transition is strongly dependent on these factors. These results suggest that the UCST transition is a result of electrostatically driven inter-complex interactions35 while the LCST transition is predominantly driven by solvation-mediated effects. Our results point to a unique system where the LCST and the UCST transitions are not encoded by the primary sequence pattern/composition, but rather by the collective properties of the RNA–polycation complexes.
In addition to divalent cations, poly(rU) RNA has previously been shown to undergo phase separation with multivalent cationic small molecules such as spermine and cationic polypeptides.37–41 However, unlike divalent cations, spermine and cationic polypeptides mediate a reentrant LLPS of RNA as the polycation-to-RNA mixing ratio is varied isothermally.39,41,42 This suggests that the complexation and phase behavior of RNA with polyamines and short cationic peptides is distinct from RNA–divalent cation mixtures. In fact, Keating and coworkers reported that poly(rU)–spermine mixtures undergo liquid–liquid phase separation with an LCST as opposed to the UCST behavior observed here for divalent cations.37 However, the temperature range at which poly(rU)–spermine mixtures were tested was up to ∼37 °C. We hypothesized that since an LCST is present in poly(rU)–spermine mixtures, a UCST must be present as well since poly(rU)–spermine phase separation is likely driven by energetically-favored inter-complex interactions.35,43 Accordingly, we prepared poly(rU)–spermine mixtures at variable compositions and inspected their phase behavior at different temperatures ranging from 4 °C to 90 °C. Remarkably, we observed that poly(rU)–spermine mixtures undergo a dual-response phase separation with both LCST and UCST transitions (Fig. 2a; Movie S1, ESI†). The LCST transition is observed at temperatures lower than the UCST transition, creating a closed-loop two-phase coexistence regime in the temperature-composition plane. Both the upper and lower cloud-point temperatures (UCPT and LCPT, respectively) of the mixture are dependent on the spermine-to-RNA ratio (Fig. 2b). At low spermine concentrations (<1 mM), the UCPT rapidly increases with the spermine concentration. Subsequently, we observe a sharp decrease in the UCPT at higher spermine concentrations (>50 mM, Fig. 2b). At intermediate spermine concentrations, the UCPT was >90 °C and therefore undetectable in our experimental setup. In contrast, the LCPT showed less sensitivity to mixture composition variation, indicating that the LCST transition might be driven by solvent-mediated interactions and not purely dependent on inter-complex interactions. We note that UCST phase transitions have been typically observed in complex coacervates where phase separation is driven by oppositely charged macromolecules.44,45 More recently, experimental observations of an LCST transition have been documented for complex coacervates37,46,47 and theoretically discussed by Muthukumar and colleagues.48 However, to our knowledge, our results report the first experimental observation of a thermo-responsive closed-loop phase diagram in an RNA-based complex coacervate system.
We next asked whether this temperature-controlled dual-response phase behavior is generic to other cationic macromolecules. To test that, we inspected the phase behavior of poly(rU) RNA with a tetra-amino acid peptide Lysine4 (K4) which carries the same formal charge (+4) as spermine. Indeed, a similar temperature-dependent dual-response phase behavior was observed for poly(rU)-K4 mixtures as a function of mixture composition (Fig. 3a and b, Movie S2, ESI†). Previous work from our lab and others indicated that phase separation in oppositely charged polypeptide–RNA mixtures is tunable by the ionic strength of the medium.36,39,41,49 Therefore, we next asked how increasing NaCl concentration in the medium impacts the two-phase coexistence region. We observed that increasing [NaCl] from 0 to 50 mM resulted in a substantial shrinkage of the two-phase regime. The UCPT in particular, but not the LCPT, was observed to rapidly decrease as the salt concentration increased (Fig. 3b–d). Finally, at ∼75 mM NaCl, LLPS was abrogated and the mixture remained soluble at all temperatures tested (90 °C ≥ T ≥ 5 °C). Overall, these data suggest that the electrostatic interactions, such as the dipolar interactions proposed by Muthukumar and colleagues, between the RNA–peptide complexes are likely to drive the UCST transition.35
In addition to salt and mixture composition variation, we measured the UCPTs and LCPTs as a function of total polymer concentration (at a fixed RNA/polycation ratio), at different pH values (pH 6–8.5), and in two independent buffers (HEPES and MOPS) which have substantially lower degree of pH variation as a function of temperature as compared to the Tris–HCl buffer.50 These control measurements showed no substantial dependence of the cloud point temperatures (both UCPT and LCPT) on the aforementioned variables (Fig. S2, ESI†).
Several studies probed the phase behavior of RNA–peptide mixtures and analogous complex coacervates with temperature variation and did not report a temperature-controlled dual-response phase behavior.37,44–47,51 We hypothesized that the size of the closed-loop is strongly dependent on the charge/size of the cation molecule. To test this idea, we examined the phase behavior of poly(rU) mixtures with a smaller cationic molecule, spermidine (Sp+3), and with a larger cationic peptide, Lysine10 (K10). We observed that spermidine–poly(rU) mixtures show a narrower LLPS region as a function of temperature (Fig. 4a). Poly(rU) samples with K10, on the other hand, remained phase-separated at all temperatures tested (4–90 °C; Fig. 4b). This observation is consistent with previous reports52 where the authors noted high temperature stability of complex coacervates. Although poly(rU)–K10 mixtures displayed phase separation at all temperature tested, we found that poly(rU)–K10 droplets became substantially smaller at 90 °C and 4 °C as compared to intermediate temperatures when K10 concentration was increased (Fig. S3, ESI†). These data indicate the proximity of both LCST and UCST transitions in this case, albeit their occurrence is outside our experimentally accessible temperature (4–90 °C) window. In summary, our results suggest that small variations in the charge and/or size of the polycation strongly affect the dimensions of the closed-loop of RNA–polycation two-phase coexistence.
In this work, we show that both types of phase behaviors (LCST and UCST) can simultaneously occur in RNA mixtures with small cationic polypeptides and polyamines in the form of a closed-loop phase diagram. Both RNA and the peptides/polyamines used here are homopolymeric without any modular architecture or complex primary sequence features. This suggests that RNA–cation complexes have distinct solvation properties and inter-complex interactions depending on the identity and size of the polycation and the mixture composition. While our experiments with small divalent cations showed a UCST phase separation, larger polycations such as polyamines showed a reentrant phase separation with temperature, comprising both UCST and LCST transitions. Manipulating intermolecular interactions with salt led to the shrinkage of the immiscibility gap by lowering the UCST boundary and keeping the LCST transition almost unaffected. The effects of salt concentration and the mixture asymmetry on the UCST of the mixture that we observed in this study are consistent with a mean-field model that considers enthalpic interactions between polyelectrolyte complexes as a basis of LLPS.35 Increasing the size of the polycation led to the extension of the immiscibility gap beyond our experimental temperature range, although, signs of the proximity of both UCST and LCST transitions were present at the extreme hot and cold temperatures tested, respectively (Fig. S3, ESI†). The increase in UCST as a function of the polycation valence is also consistent with the model proposed by Adhikari and colleagues.35
What are the molecular parameters that drive LCST and UCST transitions of RNA solutions? Ruff and co-workers presented a conceptual framework for the thermo-responsive phase behavior of biomolecules.17 Although this review mainly focuses on the phase behavior of intrinsically disordered proteins (IDPs), the framework could be extended to nucleic acid phase separation. A UCST transition could be conceptualized utilizing simple Flory–Huggins formalism where the inter-polymer interactions along with the solvent–solvent and solvent–polymer interactions are considered in Flory's mean-field parameter, called χ.53 To account for the temperature dependence of phase separation, χ can be empirically formalized as where A (represents the entropic part), B, and C are constants.46 Depending on the values and signs of the coefficients associated with the enthalpic part (B & C), various shapes of the phase diagram are possible including UCST behavior (B > 0; C = 0), LCST behavior (B < 0; C = 0), or a closed-loop (B < 0; C < 0). Although it is physically intuitive to conceptualize the UCST behavior based on the Flory–Huggins model by considering weakened polymer–polymer interactions as the temperature is increased, the LCST behavior is much more subtle and requires further attention. Extensive work by Vause and Walker, by Tanaka, and others in the 1980s and 1990s suggested that favorable hydrogen-bonding interactions between the polymer chain and the water molecules at lower temperatures underlie the occurrence of an LCST transition.54,55 Subsequently, theoretical work integrating phase separation in the framework of the Flory–Huggins model and temperature-dependent solvation via H-bonding interactions successfully captured the closed-loop shape of the experimentally determined phase diagram of aqueous polymer mixtures.55,56
In addition to the temperature dependence of the χ parameter, further considerations are needed for describing the thermo-responsive phase behavior of complex coacervates. Several theories have suggested that enthalpy-driven phase separation of oppositely charged polyelectrolytes is a two-step condensation process where the second step is the phase separation and it occurs through inter-complex interactions.35,48,57–60 In a recent work by Muthukumar and coworkers, the authors proposed the formation of polyelectrolyte complexes with dipolar trains upon mixing of oppositely charged polyelectrolytes.35 These complexes undergo phase separation via dipole–dipole inter-complex attraction. This theory successfully recapitulates the known effects of salt and mixture composition on the phase separation behavior of oppositely charged polyelectrolytes, including the occurrence of reentrant phase transition upon varying mixture stoichiometry.41,61,62 The extension of this theory suggested the possibility of both UCST, LCST, and the coexistence of UCST and LCST based on the interplay between polymer–solvent interactions and inter-complex dipolar attraction.48 The authors argued that increasing temperature can increase the strength of dipolar inter-complex interactions due to the lowering of the dielectric constant and hence favor liquid–liquid phase separation at high temperatures (LCST behavior). This model is supported by the experimental observation of LCST behavior in the mixtures of oppositely charged strong polyelectrolytes (MW ∼150–200 kDa).46 In another scenario where polymer–solvent interactions also increase with temperature, the model predicts the coexistence of UCST and LCST with LCST being higher than UCST due to the dominance of dipolar interactions at extremely high temperatures.48 Interestingly, this theory does not discuss a closed-loop phase diagram since for that to occur, phase separation needs to be suppressed at high temperatures.
Based on the aforementioned discussion, we now describe a possible molecular mechanism that leads to the existence of a closed-loop phase diagram in the temperature-composition plane of a mixture of RNA and small polycations. We argue that phase separation in these systems is a two-step process that starts with the formation of polyelectrolyte complexes followed by the condensation of these complexes into mesoscopic condensates. The RNA–polycation complexes presumably form due to the entropic gain of counterion release.63–65 The complexation occurs in the form of exchanging the counterions (which are small) with the oppositely charged polycations (which are larger) and hence the entropic gain of releasing counterions is larger than the entropic cost of constraining the polycations. Since this process is entropically driven, it is expected to be weakened at lower temperatures since the entropic contribution to the free energy of the system becomes less relevant at low temperatures53 (Fig. 5a). It is also possible that such complexation is not favored at lower temperatures due to the enthalpic cost of reducing the polymer–water hydrogen bond interactions that oppose the entropy gain of counterion release.54,55,63 As the temperature is increased, the entropic gain of exchanging counterions with the polycations becomes larger, dominating the enthalpic cost of breaking the polymer–water interactions, and hence more RNA–polycation complexes form (Fig. 5b). At the lower critical solution temperature (LCST), the population of complexes becomes sufficient to drive phase separation through inter-complex interactions (Fig. 5c) as proposed in the dipolar interaction model35 or the electrostatic correlation energy model described by Zhang, Nguyen, and Shklovskii.43,58–60 We note that, for stable dipolar trains to occur, the oppositely charged polyelectrolytes need to have symmetric or nearly symmetric lengths,35 which is not the case for RNA–polycation mixtures since the polycations used in our study are significantly smaller than the RNA chain (average length of poly(rU) is ∼2200 nucleotides). Therefore, we argue that the correlation energy model is more appropriate to describe our system mainly due to the large asymmetry of size between the components. This condensation has an entropic cost due to the partitioning of the complexes into phase-separated liquid droplets as well as the differential partitioning of salt molecules within these two phases. At higher temperatures, the entropic cost of bringing the complexes together into condensates becomes higher than the enthalpic gain of inter-complex attractions and hence the droplets dissolve, leading to a homogenous solution of RNA–polycation complexes (Fig. 5d–f). Further, the valence of the polycations and the RNA may change at high temperatures due to different protonation states, which can also affect the enthalpic gain of complexation.66 When the size of the polycation is increased (such as K4-to-K10 substitution), the entropic gain of forming complexes via counterion release becomes much larger than the entropic cost of constraining the polycations (due to the differences in size) leading to a lower temperature threshold of complex formation and subsequently a lower temperature threshold for phase separation (lowering of LCST). This may explain why mixtures of RNA with larger polycations (such as peptides) only show UCST transitions within the testable experimental temperatures.67 At the same time, the increased size of the polycation leads to stronger electrostatic correlation energies (or the formation of more dipoles in the complexes) which subsequently leads to a higher UCST. These arguments are qualitatively consistent with the experimental results presented in this study (Fig. 2–4). Increasing salt concentration strongly affects the UCST since it affects the correlation attractions and may impact the complexation equilibrium.63
In summary, we present RNA–polycation mixtures as suitable model systems to study the interplay between solvent-mediated interactions and intermolecular electrostatic interactions underlying the temperature-dependent dual-response phase behavior of the system. Our experimental results indicate that both types of interactions are significant and give rise to the coexistence of an LCST transition and a UCST transition within the experimentally accessible temperature range (4–90 °C). Together with previous studies,35,37,41,48 our data provide experimental evidence of rich phase behavior of RNA–polycation mixtures that is profoundly modulated by the charge/size of the cation and the composition of the mixture. Further, our findings add to the rich literature of the phase behavior of oppositely charged polyelectrolytes and complement the above-mentioned seminal experimental and theoretical findings. Overall, these results provide new insights into the molecular driving forces that lead to temperature-dependent phase separation in multicomponent mixtures where phase separation is driven by heterotypic interactions. Further development of the theories of complex coacervates will provide a better and more general understanding of the molecular mechanisms underlying the UCST and LCST transitions in RNA–polycation mixtures.
Footnotes |
† Electronic supplementary information (ESI) available. See DOI: 10.1039/d1sm01557e |
‡ These authors contributed equally. |
This journal is © The Royal Society of Chemistry 2022 |