Abstract
The phase separation of intrinsically disordered proteins is emerging as an important mechanism for cellular organization. However, efforts to connect protein sequences to the physical properties of condensates, i.e., the molecular grammar, are hampered by a lack of effective approaches for probing high-resolution structural details. Using a combination of multiscale simulations and fluorescence lifetime imaging microscopy experiments, we systematically explored a series of systems consisting of diblock elastin-like polypeptides (ELP). The simulations succeeded in reproducing the variation of condensate stability upon amino acid substitution and revealed different microenvironments within a single condensate, which we verified with environmentally sensitive fluorophores. The interspersion of hydrophilic and hydrophobic residues and a lack of secondary structure formation result in an interfacial environment, which explains both the strong correlation between ELP condensate stability and interfacial hydrophobicity scales, as well as the prevalence of protein-water hydrogen bonds. Our study uncovers new mechanisms for condensate stability and organization that may be broadly applicable.
Introduction
Biological condensates are found in both the cytosol1–4 and nucleus,5–10 playing essential roles in a variety of cellular processes11,12 from stress response13 to genome organization.14,15 Similar to membrane-bound organelles, they assemble a collection of molecules to raise the efficiency of sophisticated tasks. The lack of a membrane barrier allows fast material exchange between condensates and the cellular environment, rendering the molecular composition and stability of condensates more prone to regulations by external signals. 16–18
Intrinsically disordered proteins (IDPs) that promote multivalent, promiscuous interactions are key drivers of condensate formation. 19–21 Multiple mechanisms, including electrostatic, cation-π, π-π, hydrogen bonding, and hydrophobic interactions, contribute to the affinity among various chemical groups.22–24 Above a threshold concentration, as predicted by the Flory-Huggins theory, 25 interactions among IDPs can drive liquid-liquid phase separation (LLPS) to produce a highly concentrated phase that nevertheless remains dynamic. The simplicity of theory, while insightful, may prove insufficient for a comprehensive understanding of biological condensates.
Much remains to be learned regarding the connection between amino acid sequences and protein phase behaviors, or the so-called “molecular grammar” of protein condensates.26–31 These systems often exhibit complex viscoelastic behaviors and substructures with a layered organization,32–35 defying the mean-field assumption in the Flory-Huggins theory. Several more advanced theories have been introduced to better account for long-lived structural features that might persist in the polymer network. The sticker and spacer model36,37 has been adopted by Pappu and coworkers to account for strong and specific interactions that might form physical cross-links among protein molecules.38 In the meantime, the block copolymer theory may explain the microphase separation that can lead to layered structures.39–45 Because of their inherent assumptions, different theories are more appropriate for some systems but not others. High-resolution structural characterizations of condensates could offer further insight into their organizational principles and the applicable theories.
Further decoding the condensate grammar could also benefit from studies of simpler systems. Natural proteins often utilize highly complex amino acid sequences, rendering the attribution of an individual residue’s contribution to collective phenomena challenging. On the other hand, elastin-like polypeptides (ELPs), which are composed of pentapeptide repeats of valine-proline-glycine-X-glycine, serve as excellent models for studying protein phase behaviors. ELPs undergo phase separation upon heating with a lower critical solution temperature (LCST).46–50 Importantly, the guest residue X can be modified to any amino acid except proline with modern engineering approaches, 51 enabling systematic characterization of the impact on condensate stability and organization upon introducing specific residues.
We combine multiscale simulations with fluorescence lifetime imaging microscopy (FLIM) to probe condensates formed by diblock ELPs and investigate the contribution of amino acid composition to condensate stability. The simulation approach allows the sampling of largescale conformational rearrangements while providing atomistic resolution to quantify the solvation environment of individual residues. It succeeds in reproducing the stability of condensates, and reveals the formation of tubular network phases similar to gyroids. Such structures deviate from weak micelles that have been proposed for diblock ELPs in solution, and result in heterogeneous microenvironments within the condensate, which we verify in-vitro with FLIM. In addition, we find that the condensate stability exhibits a striking correlation with a hydrophobicity scale derived from interfacial transfer free energy, supporting an interfacial microenvironment of the condensate interior. The chemical specificity of the microenvironment is dictated by the peptide sequence that prevents a complete microphase separation between hydrophobic and hydrophilic groups. As a result, all condensates remain highly solvated after phase separation, producing water molecules that maintain hydrogen bonding with the exposed peptide backbone which persists even in residues with hydrophobic side chains. Our study takes a significant stride at connecting amino acid sequences with the structural and physical properties of biological condensates.
Results
Multiscale simulations of diblock ELP condensates
Decoding the condensate grammar necessitates connecting protein sequences with collective physical properties. ELPs stand out because of their sequence simplicity and amenability for biological engineering, allowing systematic exploration and precise attribution of amino acid contributions. We focus on diblock ELPs of sequence (V-P-G-V-G)n-(V-P-G-X-G)n, which we abbreviate as VnXn (Fig. 1A). The guest residue in the first block is set as valine to promote phase separation, and we explore 20 systems in which every natural amino acid is substituted into the X position in the second block.
We adopt a multiscale approach to balance accuracy and efficiency for simulating ELP condensates (Fig. 1B). The first stage of this strategy is to simulate condensate formation using a coarse-grained force field, MOFF,52–55 with one bead per amino acid and implicit solvation. These simulations are highly efficient for system relaxation and promoting largescale conformational rearrangements that occur over slow timescales. We further converted these α-carbon based models to approximately four heavy atoms per bead resolution, and introduced explicit water and ions for simulations with the MARTINI force field.56 The latest version of this force field (MARTINI3) provides a more balanced set of parameters for protein-protein interactions and has been applied to study biological condensates.57–60 These simulations further relax the protein configurations and the partition of water and counterions in condensed and dilute phases. Finally, we carried out explicit solvent all-atom simulations for 250 ns to produce an accurate characterization of the chemical environment of the condensates.
Utilizing the simulated condensate conformations, we computed various quantities to benchmark against experimental measurements. While the critical temperature has been widely used as a measure for condensate stability, determining it computationally is expensive. As an alternative, we computed the surface tension, τ, using 100-µs-long MARTINI simulations performed with the NPNAT ensemble.61 As detailed in the Supplemental Theory in the Supporting information, an inverse relationship is expected between τ and the critical temperature, Tc, for systems exhibiting LCSTs. We further approximate Tc with the transition temperatures (Tt) of ELP sequences,46 which are the temperatures at which ELPs undergo an LCST transition at a specified solution condition. Tt was shown to be linearly proportional to TC.62,63 As expected, a negative correlation can be readily seen between computed surface tension and experimental Tt (Fig. 1C). This observed negative correlation between Tt and τ supports the simulation approach’s accuracy in reproducing the sequence-dependent changes in ELP phase behavior.
We note that the experimental values were determined using different ELP sequences from those simulated here, and the relation between the surface tension and Tc is likely nonlinear,64 both contributing to the imperfect agreement. In another study,65 we adopted the same multiscale approach to determine the dielectric constant for several ELP condensates. The simulated values agree well with those determined using FLIM experiments.
Microphase separation of ELP condensates
Upon validating the accuracy of the simulation protocol for reproducing collective properties of ELP condensates, we examined their structural organization. Hassouneh et al. argued that, in dilute solutions, ELP diblocks self-assemble into the so-called weak spherical micelles with dense cores and almost unstretched coronas.66 At phase separation conditions, whether similar structural features are preserved in condensates remains unclear.
We found that most ELP condensates undergo a microphase separation. Representative configurations for four typical systems from MARTINI simulations are shown in the top panels of Fig. 2, with the V blocks in blue and X in red. Overall, the microphase separation manages to orient the more hydrophobic blocks towards the condensate interior as opposed to the interface. The relative density of guest blocks near the center of the condensate positively correlates with their hydrophobicity, as can be seen in both MARTINI (Fig. S1 in the Supporting Information) and all-atom (Fig. S2 in the Supporting Information) simulations.
Surprisingly, microphase separation did not produce lamellar morphology as expected for block copolymers with equal volume fraction of the two blocks (Fig. S3 in the Supporting Information).39–45,67–69 In particular, the condensates appear to form gyroid-like structures, in which the V and X blocks form two interpenetrating networks (Fig. S4 in the Supporting Information). This morphology also differs from micelle-like structures seen in simplified hydrophobic-polar (HP) polymers.70,71 It promotes interfacial contacts while maintaining substantial self-interactions as well. Weak interfacial tension between different ELP blocks has also been noted by Hassouneh et al.66 These qualitative observations are insensitive to the polymer length and system size in our simulations (Fig. S5 in the Supporting Information). While longer diblock copolymers drive more prominent microphase separation, similar gyroid structures can also be observed. Notably, for more hydrophilic X blocks, the condensates begin to dissolve into a collection of micelle-like structures, consistent with the predictions by Hassouneh et al. in dilute solutions.66
Quantitative characterization of the condensate interior supports the presence of interpenetrating networks with substantial interfacial contacts as well. For example, significant contacts between V and X blocks can be readily seen in the inter-chain contact maps (Fig. 2, Fig. S6 in the Supporting Information). We further computed the radial distribution functions for amino acids in various blocks (Fig. 2, Fig. S7 in the Supporting Information). For all systems, we found that the probability of finding cross-block contacts V-X is often comparable to the intra-block contacts between X-X or V-V, depending on the relative hydrophobicity of V and X. For example, the more hydrophobic F blocks exhibit strong self-clustering, much more prominent than V-V blocks that are comparable to V-F blocks. On the other hand, clustering among V-V blocks and V-A blocks is more substantial than among A-A blocks. These results are again robust with respect to system setups (Fig. S5 in the Supporting Information).
To experimentally test the microphase separation behavior uncovered in simulations, we studied the micro-physicochemical properties of the V-end and X-end of the peptides. We constructed diblock peptides with the combination of 30 pentameric repeats of V block and X (A or G) block, namely V30A30 and V30G30 (Experimental Sequences Section in the Supporting Information). The amino-termini of V30A30 and V30G30 sequences were subsequently labeled with environmentally sensitive BODIPY or SBD fluorophores,72,73 whose lifetime could be measured to quantify the viscosity or polarity of the V-end (Fig. 3A, left panel).74 These probes have been reported to be only sensitive to single physicochemical properties.65,73 To avoid artifacts induced by fluorophore labeling, we usually used ELPs labeled with a low fraction of dyes. We also constructed A30V30 and G30V30 diblock peptides, wherein the viscosity or polarity of the A-end or the G-end could be measured by fluorophores that are attached at the amino-terminus (Fig. 3A, right panel).
Using FLIM, we found that the lifetime of BODIPY for the V-end (5.43 ns) was longer than that for the A-end (4.35 ns), suggesting that the V-end indeed has a higher microviscosity than the A-end (ηV = 2233.54 cp vs ηA= 969.57 cp). Accordingly, the lifetime of SBD was longer for the V-end (8.75 ns) than the A-end (7.00 ns), indicating that the micropolarity of the V-end was lower than the A-end (ɛV = 13.25 vs ɛA = 18.97). These observations could be largely attributed to the greater extent of dehydration at the V-end due to its higher local peptide density. We further showed that the observed differences are not results of possible artifacts arising from any subtle distinctions between the two sequences V30A30 and A30V30 (Experimental Characterization of ELP Condensates Section in the Supporting Information, Fig. S8-S9 in the Supporting Information). Similar results were observed using the V-G sequences. FLIM experiments revealed that the V-end was more viscous than the G-end (ηV = 2972.72 cp vs ηG= 1958.60 cp) and the V-end was less polar than the G-end (ɛV = 9.14 vs ɛG = 27.50). These experimental observations provided the first line of evidence to support the microphase separation, as suggested by the simulation results. It is worth noting that the reported values, although related, may not quantitatively represent the steady-state viscosity. This discrepancy arises from the slow relaxation timescale inherent in ELP condensates with viscoelastic properties.
To further investigate these phenomena in real-time in motion, optical tweezers were used to stretch the droplet to observe the viscosity changes of different blocks. A single diblock condensate was captured by two optical traps (trap 1 and 2) at the same position, and it was stretched as trap 1 moved away from the original position (Fig. 3C-D). We used the peptides V30A30 and A30V30, with BODIPY labeled at their amino-terminus, to study the microviscosity changes where the fluorescence intensity decreases as the viscosity decrease. We found that the normalized fluorescence intensity of the V-end decreases greater than the A-end, indicating that the viscosity change of the V-end is larger (Fig. 3E). We anticipate that during stretching, condensates deform and molecules may deviate from the preferred chemical environment established with microphase separation. Such deviations impact the hydrophobic segments more, leading to more significant reduction in their interactions with surrounding molecules and correspondingly their viscosity.
Frustration produces condensate interior with interfacial properties
The strong dependence of molecular organization on amino acid hydrophobicity suggests that the solvation environment of individual residues might be a determining factor for condensate stability. Indeed, as shown in the Supplemental Theory of the Supporting Information, the critical temperature is closely related to the free energy cost of transferring polymer beads from a solution state to a polymer-only environment. This transfer free energy is often used to quantify the hydrophobicity of amino acids.75–81 To explore their relationship more quantitatively, we compared the transition temperature for ELP condensates measured by Urry46 to several hydrophobicity scales. Fig. 4A shows that water-octanol and water-POPC-interface transfer free energies best correlate with the transition temperature. Meanwhile, other hydrophobicity measures, even including those derived specifically for disordered proteins, do not match as well. We note that our conclusion here is not dependent on the limited number of hydrophobicity measures examined in this study, and can be seen in a previous clustering of 98 different hydrophobicity scales.82
The correlation analysis of the Urry transition temperature with various hydrophobicity scales supports a unique chemical environment emerging from the collective behaviors of condensates. To better understand the microenvironment surrounding guest residues in X blocks, we computed the radial distribution functions using MARTINI simulations. As shown in Fig. 4B-C and Fig. S10 in the Supporting Information, the guest residues are immediately surrounded with hydrophilic glycine residues as the nearest neighbors. In the second shell, more hydrophobic valine and proline residues appear. These profiles reveal the interfacial nature of the chemical environment surrounding the guest residues consisting of both hydrophobic and hydrophilic amino acids. This interfacial solvation environment is consistent with the strong correlation of the Urry scale with water-POPC-interface transfer free energy. Therefore, ELP condensates do not simply feature a hydrophobic interior that separates from hydrophilic residues as seen for folded proteins, and hence the poor correlation shown in Fig. 4A with the burial propensity of amino acids (RSA1, RSA2).
The ELP peptide sequences partly dictate the interfacial environment. For example, the guest residues are juxtaposed between two glycine residues in the pentamer (Fig. 1A), creating the immediate shell of hydrophilic residues. However, the presence of valine in the pentamer prevents the movement of GXG motifs to the interface due to the clustering of hydrophobic residues that prefer the condensate interior. Therefore, while microphase separation occurs in ELP condensates, frustration remains in the system. Hydrophilic residues cannot completely separate from hydrophobic ones due to constraints imposed by the amino acid sequence, creating unique microenvironments. Such a microenvironment arises from the collective behavior of many proteins and can deviate from individual chains,89 explaining the less-than-ideal correlation between the Urry scale and the hydrophobicity scales optimized to reproduce the conformational ensemble of IDP monomers.
Solvation environment from atomistic simulations support interfacial properties
The interfacial environment uncovered in the previous section suggests the presence of significant solvation in the condensate interior. For condensates with lower surface tension, such as V5G5, the protein density near the center of the condensate is only around ∼ 30%. Even for the stable condensates with high surface tension, such as V5F5, the protein density near the center only reaches ∼ 70%. The relative density of protein, water, and ions in different systems generally follows the trend observed for the transition temperature (Fig. S11-S12 in the Supporting Information).
The large water mass fractions in the condensates prompt the question of the extent to which water is engaged in meeting the protein’s hydrogen bonding requirements. We determined the average number of hydrogen bonds per residue with ELP residues or water from extensive atomistic simulations. As shown in Fig. 5A, over 75% of the hydrogen bonds are between protein and water, in all of the condensates. Prior work on a different ELP condensate47 similarly indicated that water-protein hydrogen bonds constitute most of the hydrogen bonds found in the system. This finding highlights that ELPs, lacking noticeable secondary structures (Fig. S13 in the Supporting Information), are less capable of meeting its own hydrogen bonding needs, as compared to folded proteins.90 Therefore, high water concentration within the condensate is crucial to fulfilling the protein’s hydrogen bonding requirements.
Strikingly, given the importance of water molecules, the average number of hydrogen bonds per residue remained relatively constant across condensates, even when the amount of water inside the condensate varied by almost two folds. Therefore, we hypothesize that each water molecule in more hydrophobic condensates must engage more in hydrogen bonding than each water in hydrophilic condensates. To test this hypothesis, we determined the water hydrogen bond density, defined as the ratio of the number of protein-water hydrogen bonds divided by the total number of water molecules near the condensate (Fig. 5D and 5E). As shown in Fig. 5B and 5C, the water hydrogen bond density indeed negatively correlates with the condensate transition temperature, Tt, supporting our hypothesis that more hydrophobic condensates would have a higher hydrogen bond density.
The water hydrogen bond density also highlights an interfacial environment of blended hydrophobic and hydrophilic regions. For example, when we limit the analysis to water molecules near the guest residues, X, the hydrogen bond density now becomes positively correlated with Tt. Therefore, locally, protein-water interactions are strongly influenced by the chemistry of side chains, and water molecules have fewer hydrogen bonding opportunities with more hydrophobic residues. However, the opposite trend is observed for the hydrogen bond density near glycine residues that immediately surround the guest residues. Therefore, the exposed backbones retain water molecules even when hydrophobic side chains do not provide hydrogen bonding opportunities, and the inseparability between the two further contributes to the observed interfacial property of condensates.
As a complementary measure of the solvation environment of individual residues, we computed the relative solvent accessibility (RSA) for the guest amino acids using atomistic simulations (Fig. 6A-B). The RSA measures the solvent-accessible surface area normalized by the maximum possible solvent exposure. RSA values in ELP condensates are much lower than those computed from all-atom simulations of ELP monomers, consistent with the increase of polymer density upon phase separation. However, they are almost always higher than the values estimated for folded proteins due to many water molecules roaming inside the condensates.87 Therefore, while ELP condensates provide a mechanism to shield hydrophobic amino acids from solvent, they are less effective in doing so than folded proteins.
While the amino acid RSA values increase upon phase separation, unlike protein folding, the increase does not result from polymer collapse. We observe a significant expansion for almost all the studied peptides upon condensation, supporting the notion that water is a poor solvent for ELP (Fig. 6C). The expansion is most evident for more hydrophobic ELPs, the radii of gyration of which are now comparable to that of an ideal Gaussian chain. Notably, the sizes of the two separate peptide blocks are similar in most systems, supporting an interfacial environment that can accommodate both hydrophobic and hydrophilic residues (Fig. S14 in the Supporting Information).
Conclusions and discussion
We carried out multiscale simulations to elucidate the connection between protein sequences and emergent properties of ELP condensates. The multiscale approach overcomes challenges in sampling slow conformational rearrangements to produce equilibrium atomistic condensate configurations. The simulations successfully reproduced condensate stability variation upon amino acid substitution. While our study is performed at set salt concentration and temperature to isolate the contributions of amino acid hydrophobicity to condensate organization, future studies may consider implementing temperature50 or salt93 dependent models to explore how solution conditions affect the organization of ELP condensates.
We found that diblock ELPs undergo microphase separation, with the more hydrophobic blocks localizing towards the condensate interior. The more hydrophobic blocks increase the local protein density, resulting in higher viscosity and lower dielectric constants, as determined in FLIM experiments. In contrast to the interior of folded proteins, these blocks are significantly solvated with water occupying over 30% of mass density. Unlike the weak micelles formed in solution,66 ELP condensates exhibit gyroid-like morphologies that promotes contacts between blocks. Further studies on the thermodynamic stability of these morphologies and comparing with predictions from the self-consistent field theory shall provide more insights into the driving forces for their emergence.
Our study uncovered a remarkably simple grammar for determining the stability of ELP condensates. A single quantity specific to the guest residue suffices to predict the emergent properties of the condensate, as evidenced by the strong correlation of the transition temperatures with the transfer free energies. The lack of higher-order effects is striking, supporting a mean field approximation to account for the contribution of individual residues in the system. However, the hydrophobicity scale is itself sensitive to the local microenvironment that is sequence specific. Such a microenvironment arises from the collective behavior of many proteins, can deviate from that of individual chains, and is likely sensitive to the solution conditions,94 which are held constant in our study. Future work on systems with double amino acid substitutions or changes to salt concentration or temperature could elucidate the generality of the mean field interpretation and the additivity of individual contributions.
While the transferability of the extracted rule for condensate stability to other systems remains to be shown, there are reasons to be hopeful. For instance, we anticipate that the microphase separation that produces the interfacial properties is a general principle for condensate organization. For most IDPs, the hydrophobic and hydrophilic residues are often interspersed throughout the chain to prevent the collapse of individual proteins95,96 and, correspondingly, the complete phase separation as in the lamellar phase. In addition, the hydrophilic backbones are expected to be exposed due to a lack of secondary structure formation, contributing to the interfacial property in the presence of hydrophobic side chains. As such, we anticipate that future research will discover the broad presence of interfacial environments within microphase separated motifs, as observed here.
Methods
Molecular dynamics simulation details
To better understand the properties of biological condensates, we performed simulations on 40 diblock ELPs with the sequence (V-P-G-V-G)n-(V-P-G-X-G)n. We set n = 5 to balance computational efficiency and polymer topological effects. As shown in Fig. S5 in the Supporting Information, the results presented in the main text are qualitatively robust with respect to the system size and peptide length.
The multiscale simulations of ELP condensates started with an implicit solvent model, MOFF,52 and a similar homopolymer model, to generate initial configurations for each system. Next, we converted the α−carbon only configurations into those consistent with the MARTINI force field57 for coarse-grained, explicit solvent simulations. These simulations measured the coexistence between the condensate and solvent, lasted for 100 µs, and were conducted in the NPNAT ensemble. Finally, we converted the MARTINI configurations into atomistic structures for simulations with CHARMM36m97 force field, with protein-water interactions shifted to better characterize IDPs. Atomistic simulations were conducted in the NPNAT ensemble and lasted 250 ns. All simulations were conducted using the software GROMACS98 and are described in detail in the Condensate simulation details section of the Supporting Information.
In addition to the condensates, we performed atomistic simulations of ELP monomers. These simulations began from RoseTTAFold predictions of the peptide structure 99,100 and lasted for 1 µs in the NPT ensemble. Full details are available in the Monomer simulation details section of the Supporting Information.
Analysis of condensate simulations
Surface tension
In explicit solvent simulations, two coexisting phases are present, in which proteins form a dense slab in the x-y plane, and the solvent creates a dilute phase in the z-dimension.101,102 We controlled the pressure normal to the protein-solvent interface by performing simulations in the NPNAT ensemble.61 The surface tension (τ) was calculated from these simulations as
where Lz is the length of the Z-dimension and N is the number of protein droplets. PXX, PYY, and PZZ are the components of the pressure tensor along the three axes.
Clustering
Since the condensate can diffuse along the z-axis, i.e., the direction perpendicular to the interface, we aligned the simulated configurations using the position of the largest cluster of peptide molecules. To find the largest cluster, we first computed a center of mass contact matrix between peptide molecules with a distance cutoff of 4 nm. Then, using a depth first search algorithm,103 we identified the size and location of the largest protein cluster. The system was then shifted such that the z coordinate of the center of mass of the largest cluster is at zero. We implemented the alignment using the software MDAnalysis.104,105
Contact map
Protein contact maps were calculated using a cutoff distance of 0.6 nm. Specifically, two residues were defined as in contact if any pair of their coarse-grained beads were within 0.6 nm.
Radial distribution function
Radial distribution functions were calculated using the InterRDF function of MDAnalysis between 0.00001 nm and 2 nm. For analyzing the microphase separation of each condensate, we divided each protein into the V-substituted blocks and X-substituted blocks. This radial distribution function excluded intra-chain contacts to better capture how inter-chain contacts determine condensate organization. To analyze the origin of chain frustration within diblock phase separation, we computed the radial distribution function of each guest amino acid to the amino acids native to the sequence. Intra-chain contacts were included to fully capture the local environment experienced by the guest amino acid.
Relative solvent accessibility
The relative solvent accessibility (RSA) was measured based on the solvent accessible surface area (SASA) of guest amino acids determined from all atom simulations. We first calculated the SASA for each residue using GROMACS tools sasa, with a probe radius of 0.14 nm.106 We then averaged across all repeats of guest amino acids, and normalized the SASA by the maximum possible value for a given amino acid to compute RSA87 (see Table S1).
Radius of gyration
The estimated Rg from atomistic simulations were compared to the following analytical expression91
where N is the length of the polymer, b = 0.38nm is the bond length, and lp = 0.40nm is the persistence length.
Secondary structure
Secondary structure for protein condensates was calculated from atomistic simulations using the GROMACS do dssp command.107,108 Fig. S13 in the Supporting Information displays the ordered fraction of secondary structures, which includes α-helices, β-sheets, β-bridges, and turns.
Hydrogen bonding
Protein-protein and protein-water hydrogen bonds were determined using the hydrogenbonds package within MDAnalysis.104,105,109 A hydrogen bond was defined as having a donor heavy atom within 3.30 Å of an acceptor heavy atom and a donor heavy atom-hydrogen atom-acceptor heavy atom angle of 135.0◦ or larger. The nearby water count for each analysis selection was determined by counting the water molecules whose oxygen atoms were within 4.0 Å of the heavy atoms of this selection. Water hydrogen bond density for each residue selection for each frame was computed as the ratio of protein-water hydrogen bonds to the number of nearby water molecules (Fig. S15 in the Supporting Information).
Experimental methods
Fluorescence lifetime imaging microscopy (FLIM)
The proteins of interest were labeled with NHS ester fluorophore. We used ELPs with 1% BODIPY labels or 2.5% SBD labels to form condensates, which avoid the artifacts induced by fluorophores. Droplets were formed with the final concentration of 70 µM ELP in 2 M NaCl for V-A and 1.5 M NH4SO4 for V-G diblock, respectively. A drop of droplets containing solution was placed on a 0.17 mm coverslip with a 500 µm spacer. Images were acquired by Leica Falcon Fluorescence Microscope equipped with Wil pulse laser and 63X/0.12 oilimmersion objective. The BODIPY was excited at 488 nm and the SBD was excited at 448 nm. The fluorescence lifetime fitting and image analysis were performed in LAS X and Image J.
Optical tweezer experiments
The protocol describes the preparation and manipulation of a sample containing a BODIPY-labeled ELP protein using an optical tweezer, specifically the C-Trap from LULICKS. The goal is to form a single large condensate from small droplets and then stretch it using the optical tweezers while monitoring its fluorescence intensity. The first step is to mix the BODIPY-labeled ELP protein with NaCl solution to obtain a final concentration of 20 µM protein in 2 M NaCl. The sample is then loaded into the flow channel of the optical tweezer. The dual trap mode of the C-Trap is then used to move the two traps to the same position. The sample is flushed until a single large condensate of about 4 µm is caught by the dual trap via the fusion of small droplets. The traps are then moved along with the droplet to another flow channel containing 2 M NaCl buffer. The channel and pressure are turned off to avoid any additional force on the droplet. The next step is to stretch the droplet by moving trap 1 0.2 µm each time. Fluorescent images of each stretch are taken using a 488 nm excitation laser. The mean fluorescence intensity of the whole droplet is then quantified using Image J.
Acknowledgements
This work was supported by the National Institutes of Health (Grant R35GM133580, B.Z.), the Research Center for Industries of the Future (RCIF) at Westlake University (X.Z), PEW Biomedical Scholars Program 00033066 (X.Z.), and the National Science Foundation under CHE-1654415 (A.P.W.). A.L. acknowledges support from the National Science Foundation Graduate Research Fellowship Program (grant number 1745302). D.S. acknowledges support from the Hertz Foundation and the National Science Foundation Graduate Research Fellowship Program (grant number 2141064).
References
- (1)Germline P Granules Are Liquid Droplets That Localize by Controlled Dissolution/CondensationScience 324:1729–1732
- (2)Regulation of biomolecular condensates by interfacial protein clustersScience 373:1218–1224
- (3)A Liquid-to-Solid Phase Transition of the ALS Protein FUS Accelerated by Disease MutationCell 162:1066–1077
- (4)A Membraneless Organelle Associated with the Endoplasmic Reticulum Enables 3’UTR-Mediated Protein-Protein InteractionsCell 175:1492–1506
- (5)Coexisting Liquid Phases Underlie Nucleolar SubcompartmentsCell 165:1686–1697
- (6)Liquid droplet formation by HP1α suggests a role for phase separation in heterochromatinNature 547:236–240
- (7)Phase separation drives heterochromatin domain formationNature 547:241–245
- (8)Coactivator condensation at super-enhancers links phase separation and gene controlScience 361
- (9)Single-stranded nucleic acid binding and coacervation by linker histone H1Nat. Struct. Mol. Biol 29:463–471
- (10)On the stability and layered organization of protein-DNA condensatesBiophys. J 121:1727–1737
- (11)Biomolecular condensates: Organizers of cellular biochemistryNat. Rev. Mol. Cell Biol 18:285–298
- (12)Intrinsically disordered proteins in overcrowded milieu: Membrane-less organelles, phase separation, and intrinsic disorderCurr. Opin. Struct. Biol 44:18–30
- (13)Stress-Triggered Phase Separation Is an Adaptive, Evolutionarily Tuned ResponseCell 168:1028–1040
- (14)Multiscale Modeling of Genome Organization with Maximum Entropy OptimizationJ. Chem. Phys 155
- (15)Biomolecular Condensates in the NucleusTrends Biochem. Sci :961–977
- (16)A phase separation model predicts key features of transcriptional controlCell 169:13–23
- (17)Composition-dependent thermodynamics of intracellular phase separationNature 581:209–214
- (18)Partitioning of cancer therapeutics in nuclear condensatesScience 368:1386–1392
- (19)Liquid-Liquid Phase Separation in BiologyAnnu. Rev. Cell Dev. Biol 30:39–58
- (20)Protein Phase Separation: A New Phase in Cell BiologyTrends Cell Biol 28:420–435
- (21)Conformational Dynamics of Intrinsically Disordered Proteins Regulate Biomolecular Condensate ChemistryChem. Rev 122:6719–6748
- (22)Biomolecular Phase Separation: From Molecular Driving Forces to Macroscopic PropertiesAnnu. Rev. Phys. Chem 71:1–23
- (23)Comparative roles of charge, π, and hydrophobic interactions in sequence-dependent phase separation of intrinsically disordered proteinsProc. Natl. Acad. Sci. U.S.A 117:28795–28805
- (24)Molecular interactions contributing to FUS SYGQ LC-RGG phase separation and co-partitioning with RNA polymerase II heptadsNat. Struct. Mol. Biol 28:923–935
- (25)Thermodynamics of high polymer solutionsJ. Chem. Phys 10
- (26)Biomolecular Condensates: Sequence Determinants of Phase Separation, Microstructural Organization, Enzymatic Activity, and Material PropertiesJ. Phys. Chem. B 125:3441–3451
- (27)Physical Principles Underlying the Complex Biology of Intracellular Phase TransitionsAnnu. Rev. Biophys 49:107–133
- (28)A Molecular Grammar Governing the Driving Forces for Phase Separation of Prion-like RNA Binding ProteinsCell 174:688–699
- (29)Elucidating binding mechanisms and dynamics of intrinsically disordered protein complexes using NMR spectroscopyCurr. Opin. Struct. Biol 54:10–18
- (30)The disordered P granule protein LAF-1 drives phase separation into droplets with tunable viscosity and dynamicsProc. Natl. Acad. Sci. U.S.A 112:7189–7194
- (31)Learning the chemical grammar of biomolecular condensatesNat. Chem. Biol
- (32)Molecular Determinants for the Layering and Coarsening of Biological CondensatesAggregate
- (33)A conceptual framework for understanding phase separation and addressing open questions and challengesMol. Cell 82:2201–2214
- (34)Phase-separating RNA-binding proteins form heterogeneous distributions of clusters in subsaturated solutionsProc. Natl. Acad. Sci. U.S.A 119
- (35)Single fluorogen imaging reveals spatial inhomogeneities within biomolecular condensatesbioRxiv
- (36)Theory of Thermoreversible GelationMacromolecules 22:1988–1994
- (37)Thermoreversible Gelation in Solutions of Associative Polymers. 1. StaticsMacromolecules 31:1373–1385
- (38)Intrinsically disordered linkers determine the interplay between phase separation and gelation in multivalent proteinsElife 6:1–31
- (39)Theory of Microphase Separation in Block CopolymersMacromolecules 13:1602–1617
- (40)Block Copolymer Thermodynamics: Theory and ExperimentAnnu. Rev. Phys. Chem 41:525–557
- (41)Stable and unstable phases of a diblock copolymer meltPhys. Rev. Lett 72:2660–2663
- (42)Unifying weak- and strong-segregation block copolymer theoriesMacromolecules 29:1091–1098
- (43)The packing of soft materials: Molecular asymmetry, geometric frustration and optimal lattices in block copolymer meltsPhys. Rep 433:1–64
- (44)Design and application of nanoscale actuators using block-copolymersPolymers 2:454–469
- (45)Frustration in block copolymer assembliesJ. Phys. Condens. Matter 33
- (46)Physical chemistry of biological free energy transduction as demonstrated by elastic protein-based polymersJ. Phys. Chem. B 101:11007–11028
- (47)The liquid structure of elastineLife 6
- (48)Thermal Compaction of Disordered and Elastin-like Polypeptides: A Temperature-Dependent, Sequence-Specific Coarse-Grained Simulation ModelBiomacromolecules 21:3523–3538
- (49)Temperature, Hydrostatic Pressure, and Osmolyte Effects on Liquid–Liquid Phase Separation in Protein Condensates: Physical Chemistry and Biological ImplicationsChem. Eur. J 25:13049–13069
- (50)Temperature-Controlled Liquid–Liquid Phase Separation of Disordered ProteinsACS Cent. Sci 5
- (51)Recursive directional ligation by plasmid reconstruction allows rapid and seamless cloning of oligomeric genesBiomacromolecules 11:944–952
- (52)Consistent Force Field Captures Homologue-Resolved HP1 Phase SeparationJ. Chem. Theory Comput 17:3134–3144
- (53)Improving Coarse-Grained Protein Force Fields with Small-Angle X-ray Scattering DataJ. Phys. Chem. B 123:1026–1034
- (54)Maximum Entropy Optimized Force Field for Intrinsically Disordered ProteinsJ. Chem. Theory Comput 16:773–781
- (55)Unifying coarse-grained force fields for folded and disordered proteinsCurr. Opin. Struct. Biol 72:63–70
- (56)The MARTINI force field: Coarse grained model for biomolecular simulationsJ. Phys. Chem. B 111:7812–7824
- (57)Martini 3: a general purpose force field for coarse-grained molecular dynamicsNat. Methods 18:382–388
- (58)Improving Martini 3 for Disordered and Multidomain ProteinsJ. Chem. Theory Comput 18:2033–2041
- (59)Simulation of FUS Protein Condensates with an Adapted Coarse-Grained ModelJ. Chem. Theory Comput 17:525–537
- (60)Coacervate formation studied by explicit solvent coarse-grain molecular dynamics with the MartiniChem. Sci
- (61)Computer simulation of liquid/liquid interfaces. I. Theory and application to octane/waterJ. Chem. Phys 103:10252–10266
- (62)A Unified Model for de Novo Design of Elastin-like Polypeptides with Tunable Inverse Transition TemperaturesBiomacromolecules 14:2866–2872
- (63)Quantification of the Effects of Chain Length and Concentration on the Thermal Behavior of Elastin-like PolypeptidesBiomacromolecules 5:846–851
- (64)Theory of the interface between polymers or polymer solutions. I. Two components systemJ. Chem. Phys 62:490–499
- (65)Micropolarity governs the structural organization of biomolecular condensatesNat. Chem. Biol :1–9
- (66)Elastin-like Polypeptide Diblock Copolymers Self-Assemble into Weak MicellesMacromolecules 48:4183–4195
- (67)Block copolymers-designer soft materialsPhys. Today 52:32–38
- (68)Self-assembly of block copolymersChem. Soc. Rev 41:5969–5985
- (69)Diblock copolymer micelles in a dilute solutionMacromolecules 38:5330–5351
- (70)Model for disordered proteins with strongly sequence-dependent liquid phase behaviorJ. Chem. Phys 152
- (71)Analytical Formulation and Field-Theoretic Simulation of Sequence-Specific Phase Separation of Protein-Like Heteropolymers with Short- and Long-Spatial-Range InteractionsJ. Phys. Chem. B 126:9222–9245
- (72)AgHalo: A Facile Fluorogenic Sensor to Detect Drug-Induced Proteome StressAngew. Chemie - Int. Ed 56:8672–8676
- (73)A dual-functional BODIPY-based molecular rotor probe reveals different viscosity of protein aggregates in live cellsAggregate :2–8
- (74)An Optical Technique for Mapping Microviscosity Dynamics in Cellular OrganellesACS Nano 12:4398–4407
- (75)Improved coarse-grained model for studying sequence dependent phase separation of disordered proteinsProtein Sci
- (76)Experimentally determined hydrophobicity scale for proteins at membrane interfacesNat. Struct. Biol 3:842–848
- (77)Solvation energies of amino acid side chains and backbone in a family of host - Guest pentapeptidesBiochemistry 35:5109–5124
- (78)Physicochemical Basis of Amino Acid Hydrophobicity Scales: Evaluation of Four New Scales of Amino Acid Hydrophobicity Coefficients Derived from RP-HPLC of PeptidesAnal. Chem 67:1210–1219
- (79)The characterization of amino acid sequences in proteins by statistical methodsJ. Theor. Biol 21:170–201
- (80)Comparing the Polarities of the Amino Acids: Side-Chain Distribution Coefficients between the Vapor Phase, Cyclohexane, 1-Octanol, and Neutral Aqueous SolutionBiochemistry :1664–1670
- (81)A simple experimental model for hydrophobic interactions in proteinsJ. Biol. Chem 259:2910–2912
- (82)50 years of amino acid hydrophobicity scales: Revisiting the capacity for peptide classificationBiol. Res 49:1–19
- (83)A simple atomic-level hydrophobicity scale reveals protein interfacial structureJ. Mol. Biol 426:484–498
- (84)Accurate model of liquid – liquid phase behavior of intrinsically disordered proteins from optimization of single-chain propertiesProc. Natl. Acad. Sci. U.S.A
- (85)A Data-driven Hydrophobicity Scale for Predicting Liquid-Liquid Phase Separation of ProteinsJ. Phys. Chem. B 125:4046–4056
- (86)Hydrophobicity of Amino Acid Residues in Globular ProteinsScience 229
- (87)Maximum allowed solvent accessibilites of residues in proteinsPLoS ONE 8
- (88)A simple method for displaying the hydropathic character of a proteinJ. Mol. Biol 157:105–132
- (89)Phase behaviour of disordered proteins underlying low density and high permeability of liquid organellesNat. Chem 9
- (90)Hydrogen bonding in globular proteinsJ. Mol. Biol 226:1143–1159
- (91)Polymer scaling laws of unfolded and intrinsically disordered proteins quantified with single-molecule spectroscopyProc. Natl. Acad. Sci. U.S.A 109:16155–16160
- (92)Relation between single-molecule properties and phase behavior of intrinsically disordered proteinsProc. Natl. Acad. Sci. U.S.A 201804177
- (93)Salt-Dependent Conformational Changes of Intrinsically Disordered ProteinsJ. Phys. Chem. Lett 12:6684–6691
- (94)Surface Electrostatics Govern the Emulsion Stability of Biomolecular CondensatesNano Lett 22:612–621
- (95)Valence and patterning of aromatic residues determine the phase behavior of prion-like domainsScience 367:694–699
- (96)Hydropathy Patterning Complements Charge Patterning to Describe Conformational Preferences of Disordered ProteinsJ. Phys. Chem. Lett 11:3408–3415
- (97)CHARMM36m: an improved force field for folded and intrinsically disordered proteinsNat. Methods 14:71–73
- (98)GROMACS: A message-passing parallel molecular dynamics implementationComput. Phys. Commun 91:43–56
- (99)Accurate prediction of protein structures and interactions using a three-track neural networkScience 373:871–876
- (100)Rosetta in CASP4: Progress in ab initio protein structure predictionProteins 45:119–126
- (101)Sequence determinants of protein phase behavior from a coarse-grained modelPLoS Comput. Biol 14:1–23
- (102)Triple-point coexistence properties of the lennard-jones systemChem. Phys. Lett 51:155–159
- (103)Analyzing and Driving Cluster Formation in Atomistic SimulationsJ. Chem. Theory and Comput 13:1317–1327
- (104)MDAnalysis: A Python Package for the Rapid Analysis of Molecular Dynamics SimulationsProceedings of the 15th Python in Science Conference
- (105)MDAnalysis: a toolkit for the analysis of molecular dynamics simulationsJ. Comput. Chem 32:2319–2327
- (106)The Double Cubic Lattice Method: Efficient Approaches to Numerical Integration of Surface Area and Volume and to Dot Surface Contouring of Molecular AssembliesJ. Comput. Chem 16:273–284
- (107)A series of PDB related databases for everyday needsNucleic Acids Res 43:D364–D368
- (108)Dictionary of Protein Secondary Structure: Pattern Recognition of Hydrogen-Bonded and GeometricalBiopolymers 22:2577–2637
- (109)On the interaction of hyaluronic acid with synovial fluid lipid membranesPhys. Chem. Chem. Phys 21:9845–9857
Article and author information
Author information
Version history
- Preprint posted:
- Sent for peer review:
- Reviewed Preprint version 1:
- Reviewed Preprint version 2:
- Reviewed Preprint version 3:
Copyright
© 2023, Latham et al.
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.
Metrics
- views
- 670
- downloads
- 46
- citations
- 0
Views, downloads and citations are aggregated across all versions of this paper published by eLife.