Cells carry out highly sophisticated tasks such as self-repair, migration, differentiation and division. These behaviors emerge from collective interactions between biomolecules as well as complex biochemical reaction networks. The robustness of cellular processes is especially remarkable when considering high molecular diversity and abundance found in the cell. The cell interior is highly crowded with numerous protein species present at an estimated concentration around 3 mM (Milo, 2013). A daunting task cells face is to coordinate various intracellular processes which take place simultaneously in the confined three-dimensional space. Cells address this seemingly chaotic situation, at least partially, through compartmentalization.

Eukaryotic cells contain numerous compartments or organelles that consist of a specific set of biomolecules and perform specialized biological functions. In addition to membrane-bound organelles such as nucleus and mitochondria, cells also harbor those without encapsulating membranes. Examples fall into the latter category include, among others, nucleoli, nuclear speckles, paraspeckles, promyelocytic leukemia (PML) bodies in nucleus and cytoplasmic stress granules, germ granules and P bodies (Buchan and Parker, 2009; Chujo and Hirose, 2017; Mao et al., 2011). These membrane-less organelles, recently termed condensates, have several interesting biophysical properties (Banani et al., 2017; Shin and Brangwynne, 2017). First, even in the absence of physical barriers such as lipid membranes, they can still enrich a set of biomolecules within them. Specific marker proteins are recruited into condensates, enabling the visualization of condensates using fluorescence microscopy (Thul et al., 2017). Second, biomolecules in these organelles are not organized in a well-defined three-dimensional arrangement, which is an important feature distinguishing them from macromolecular complexes such as ribosome and RNA polymerases of which individual subunits make tight contacts with one another in a designated stoichiometry and orientation (Freeman Rosenzweig et al., 2017; Wilfling et al., 2020). Moreover, condensates such as nucleoli and stress granules tend to be much larger than typical macromolecular complexes, reaching up to a few microns in size. Third, internal components of condensates are not completely confined inside the organelles, as examined with experimental techniques such as fluorescence recovery after photobleaching. Although the degree of molecular mobility varies depending on the type of biomolecules and the physical state of condensates, components often exhibit diffusive mobility within the organelle and also undergo dynamic exchange with the surrounding cellular space (Chen and Huang, 2001; Mollet et al., 2008; Weidtkamp-Peters et al., 2008). These data collectively suggest that membrane-less compartments are not a mere aggregate of biomolecules. Rather, they represent a condensed state of living matter where constituent biomolecules are cohesively held together by intermolecular interactions that constantly break and reform.

Recently, a biophysical principle driving the formation of membrane-less organelles has begun to be uncovered. In 2009, P granules in Caenorhabditis elegans embryo were shown to exhibit liquid-like material properties (Brangwynne et al., 2009). When the one-cell embryo divides, P granules spatially segregate into one of daughter cells, a process implicated in germline specification. Quantitative image analysis of the segregation process revealed that the segregation was not a result of transport of intact P granules to one side of the cell, but rather a combined outcome of localized condensation and dissolution. Moreover, when cytoplasmic flow was induced, P granules that initially wet the surface of nuclei became detached and exhibited dripping as well as fusion, a hallmark of liquid states. A later work showed that the nucleoli of Xenopus laevis oocytes had round morphology and upon fusing one another, the dumbbell-like shape gradually relaxed back to round one (Brangwynne et al., 2011). These behaviors are again consistent with liquid-like nature displaying surface tension which drives the system toward the minimum of surface area. The presence of liquid-like organelles within cyto- or nucleoplasm immediately suggests that liquid-liquid phase separation can be a physiochemical principle for the biogenesis of membrane-less organelles. Numerous works afterwards have shown that phase separation indeed plays a key role in organizing intracellular space (Li et al., 2012; Molliex et al., 2015; Nott et al., 2015; Patel et al., 2015) and the list of biological processes involving phase separation is continuously expanding (Feric et al., 2016; Franzmann et al., 2018; Sabari et al., 2018; Strom et al., 2017; Su et al., 2016; Yasuda et al., 2020; Zeng et al., 2016; Zhang et al., 2018).

In this paper, we review fundamentals of phase separation processes with a primary focus on diverse phase diagrams exhibited by biomolecular phase separation. We begin our discussion with the regular solution model of a binary mixture composed of one type of solute dissolved in solvent. This simple model explains a key aspect of phase separation such as the concentration-dependent assembly. We then build up complexity by considering a ternary regular solution model with additional type of solute. Special emphasis is given to rich behaviors the ternary phase-separating systems can exhibit. We discuss how intermolecular interactions can define the specific shape of phase diagrams of the ternary mixture. In doing so, we motivate a model-based approach in analyzing experimental data and designing new experiments that can be useful in gaining new insights into the phase behaviors of cells, an inherently multicomponent system.

Phase separation is commonly encountered in everyday experience in mixtures such as oil-water, oil-vinegar, and metal alloy systems. The simplest model capturing the essence of phase separation is a lattice-based regular solution model for a binary mixture, also sometimes called Flory-Huggins model (Rubinstein and Colby, 2003). The model deals with a lattice of which individual site is filled with either solute or solvent species (Fig. 1A). The volume fraction of the solute, the number of solutes divided by the total number of lattice sites, is a major parameter that can be controlled as an input in the model. The phase behavior of the binary regular solution model is dictated by the free-energy of mixing, comprising both entropic as well as energetic contributions (Brangwynne et al., 2015; Posey et al., 2018). In the absence of any interactions between species, the solutes stay in a well-mixed state without any phase separation (Figs. 1A and 1B). This is because the system can minimize free energy by maximizing possible configurations, equivalently entropy, in the well-mixed state. Energetic interactions between species can either promote or impede mixing. In the regular solution model, interaction energies between only nearest neighbors are taken into consideration under mean-field approximation. When homotypic interactions are favored over heterotypic ones, the binary system can decrease its free energy through phase separation, dividing system volume into two regions with distinct compositions (Figs. 1A and 1B).

These conditions for phase separation can be conveniently expressed in the form of a phase diagram where the X and Y axis typically represent the volume fraction or concentration of phase-separating solutes and temperature, respectively. In the phase diagram of binary systems, a curved line called binodal or coexistence curve delimits a single-phase versus a two-phase region (Fig. 1C). Thus, for a given temperature and solute type (a fixed Y value in the phase diagram), there exists a concentration threshold, called saturation concentration, above which phase separation is induced (Fig. 1D). Phase separation divides the system volume into two distinct regions, a dense phase with a high concentration of solutes and a dilute phase with low concentration. The concentration of solutes within each phase is again specified in the phase diagram: the dilute phase concentration is fixed to saturation concentration and the concentration in the dense phase is given by the right-hand side of binodals (Fig. 1D). Importantly, as the system enters deeper into the two-phase region, the model states that solute concentrations in each phase do not change, but the volume of the dense phase increases. The degree of phase separation can be also modulated by changing either temperature or interaction strengths between species. For example, at higher temperature (or when homotypic interactions become less favorable), the saturation concentration increases, and thus higher concentrations of solutes are required for phase separation. Taken together, the regular solution theory provides several detailed features that can be used to test whether binary phase separation takes place in the system being probed.

Liquid-liquid phase separation of a binary mixture of protein and aqueous solution has been observed at least as early as the 1970s (Ishimoto and Tanaka, 1977). For lysozyme and lens crystallin proteins, temperature-dependent phase separation behaviors were observed, and their coexistence curves were quantitatively measured (Broide et al., 1991; Thomson et al., 1987). The phase behaviors of these proteins are studied in relation to Cataract diseases (Benedek, 1997). More recently, numerous studies have reported binary phase separation of proteins harboring intrinsically disordered regions (IDRs) (Brangwynne et al., 2015; Uversky, 2021). IDRs represent polypeptide segments lacking well-defined three-dimensional folding (Hong et al., 2020). Their sequences tend to be simple in composition, and thus IDRs are closely related to low-complexity domains (Tompa, 2012). Recently, phase separation of IDR-containing proteins (IDPs) has attracted massive attention since many components of condensates, particularly RNA- and DNA-binding proteins, are known to harbor IDRs (Kato et al., 2012; Uversky, 2017). Purified IDPs undergo phase separation to form dense liquid droplets (Lin et al., 2015; Molliex et al., 2015; Patel et al., 2015). Notably, several features predicted from the regular solution model have been indeed observed, including the presence of saturation concentrations and the characteristic temperature-dependence of phase separation (Burke et al., 2015; Molliex et al., 2015; Nott et al., 2015). In particular, an in-vitro work on LAF-1, a P-granule component, demonstrated that the LAF-1 droplets are in equilibrium with a saturated protein solution outside of droplets by showing that the protein concentrations in the dilute phase are fixed at the saturation concentration even when total protein concentrations vary (Elbaum-Garfinkle et al., 2015).

Weak multivalent interactions mediated by residues along polypeptide chains are thought to be important for phase separation of IDPs, as often modeled with the sticker-spacer framework (Choi et al., 2020) (Fig. 1E). Stickers are a segment of polypeptides exhibiting attractive interactions with other stickers, acting as a key element driving phase separation. Spacers link neighboring stickers and play a modulatory role in phase separation (Harmon et al., 2017). Specific residues corresponding to stickers vary depending on the type of interactions biomolecules exhibit. Diverse interaction modes including electrostatic (Nott et al., 2015; Pak et al., 2016), cation-pi (Wang et al., 2018), and hydrophobic interactions (Murthy et al., 2019; Simon et al., 2017) have been shown to play key roles in promoting phase separation (Dignon et al., 2020). The regular solution model suggests that altering intermolecular interactions can impact phase separation of solutes (Fig. 1F). Indeed, screening electrostatic interactions using salt ions leads to a decrease in droplet sizes and eventual dissolution for electrostatically-driven condensates (Elbaum-Garfinkle et al., 2015; Nott et al., 2015; Zamudio et al., 2019). Interactions between stickers are also highly affected by mutations and post-translational modifications (PTMs), which can ultimately cause a significant change in phase behaviors (Hofweber et al., 2018; Monahan et al., 2017; Qamar et al., 2018; Ryan et al., 2018; Snead and Gladfelter, 2019). Moreover, cells utilize PTMs of condensate components as a regulatory mechanism for condensate assembly (Rai et al., 2018). Notably, although the regular solution model is successful in describing the overall phase separation properties of binary mixtures, a more advanced theory in combination with computer simulation is often required for dissecting quantitative aspects of phase behaviors such as the sequence-specific locations of the coexistence region (Lin et al., 2016; Wei et al., 2017).

For living cells with numerous types of biomolecules, the binary regular solution model appears to be over-simplified. However, several works have shown that the model still provides a biophysical explanation for key features of condensate assembly. Nucleoli in C. elegans are shown to assemble only when the concentration of nucleolar protein FIB-1 reaches a concentration threshold, consistent with the presence of the saturation concentration (Weber and Brangwynne, 2015). A recent live cell reconstitution work also shows that stress granules assemble with a clear concentration threshold for G3BP (Sanders et al., 2020). DDX4 is one of model IDPs of which phase behaviors are extensively studied (Nott et al., 2015; 2016). When expressed in cells, the DDX4 IDR concentration in dilute phase is maintained at a similar level even though the total protein concentration increases over an order of magnitude (Klosin et al., 2020), an effect fully consistent with the binary regular solution model. This behavior suggests that condensates can function to buffer protein expression noise. Similar behaviors of concentration thresholding are also observed for light-inducible systems such as optoDroplets (Shin et al., 2017; Wei et al., 2020) and Corelets (Bracha et al., 2018). Since the degree of super-saturation can be precisely defined using external light, these systems are particularly useful to quantitatively study phase separation in living cells (Basu et al., 2020; Zhu et al., 2019).

Adding another type of solute into the system leads to drastically richer phase behaviors, as hinted by a ternary regular solution model (Deviri and Safran, 2021; Dignon et al., 2020; Meijering, 1950). With additional components in the system, the phase diagram now incorporates another axis for the concentration of the second solute, and becomes three-dimensional (Fig. 2A). Typically, concentrations of coexisting phases are described by a two-dimensional slice at a given temperature, and thus phase diagrams of ternary systems are usually plotted in the concentration-concentration plane. Here, we consider five different scenarios of inter-species interactions, each of which leads to phase diagrams of distinct shape and molecular partitioning.

The first scenario involves two types of solutes which exhibit strong heterotypic interactions but not homotypic ones (Fig. 2B). Since each solute does not undergo self-association, a solution with only one type of solute stays in a well-mixed state. However, when both types of solutes are present, attractive heterotypic attractions can lead to phase separation. Thus, a phase diagram for this type of system typically forms a closed loop detached from both axes. When the solute concentrations fall in the coexistence region enclosed by the loop, phase separation occurs to generate two coexisting phases: a dense and a dilute phase. The concentrations of solutes in each phase are generally given by tie-lines. The tie-line of the system in the first scenario has a positive slope, indicating that both solutes are enriched in the dense phase and depleted in the dilute phase. A mixture of polycation and polyanion, such as positively-charged polypeptides and RNA, is a classic example in this category (Aumiller and Keating, 2016; Banerjee et al., 2017; Ukmar-Godec et al., 2019). More recently, pairs of tandem protein repeats with interacting modules, such as PRM_X-SH3_Y, SUMO_X-SIM_Y, and pY_X-SH2_Y (X and Y are the number of repeats), have been shown to play important roles in intracellular phase separation (Banani et al., 2016; Li et al., 2012). These multivalent proteins are found in various cell signaling pathways including T cell activation (Su et al., 2016) and neuronal synapses (Chen et al., 2020). Concentrating reaction components within condensed liquid droplets is shown to enhance biochemical reactions such as actin assembly (Su et al., 2016) as well as enzymatic activity (Peeples and Rosen, 2021).

The second case is a mixture of a self-associating solute and the other one exhibiting neither homo- nor heterotypic interactions. The self-associating solute can phase-separate even in the absence of the other solute, but the latter can promote phase separation through excluded volume interactions (i.e., two types of solutes cannot occupy the same volume simultaneously). Thus, the inert solute acts as a crowder. Typically, the coexistence region in this case intersects one of the axes, corresponding to the self-associating solute, and tie-lines are oriented with a negative slope (Fig. 2C). Thus, the dense phase consists mostly of the phase-separating solute, devoid of crowders. Examples in this scenario include purified IDPs in the presence of crowding reagents such as PEG or dextran (Cai et al., 2019; Lin et al., 2015; Sabari et al., 2018). Increasing the concentration of crowders further promotes phase separation, consistent with a decrease in the saturation concentration (Ray et al., 2020; Zhang et al., 2020). For IDPs with weaker homotypic interactions, crowders are required for phase separation (Kanaan et al., 2020; Patel et al., 2015; Sabari et al., 2018). The phase diagram for these proteins can also be described with a similar coexistence region but detached from the axis for the self-associating solute. The cell interior is highly crowded with biomolecules that do not undergo any specific interactions with many others. The crowded intracellular environment is likely to play a similar role in enhancing biomolecular phase separation.

Adding attractive heterotypic interactions to the second scenario, particularly relevant in the context of phase separation control in the cell is a mixture of a self-associating scaffold solute and a client solute exhibiting heterotypic interactions but minimal homotypic ones. Although the detailed modulatory effects of clients on scaffold phase behavior vary depending on the exact nature of scaffold-client interactions (Ruff et al., 2021), the ternary model predicts that, unlike crowders, client solutes can function in two opposite manners. When a small amount of client solutes is around, phase separation of scaffolds is enhanced due to additional crosslinks clients provide (Fig. 2D). However, an excess of clients can inhibit phase separation since interactions between clients and scaffolds begin to interfere with scaffold-scaffold interactions. This effect has been observed in a mixture of RNA and RNA binding proteins. RNA binding proteins often harbor positively charged IDRs such as RGG motifs together with various RNA recognition motifs (Wiedner and Giudice, 2021). Interactions between these motifs and RNA can strengthen a network of inter-molecular interactions driving phase separation. Indeed, multiple works have shown that RNA promotes phase separation of RNA binding proteins such as FUS, hnRNPA1, and G3BP1 (Guillén-Boixet et al., 2020; Maharana et al., 2018; Molliex et al., 2015; Yang et al., 2020). However, further increasing RNA concentrations begins to adversely affect phase separation, leading to shrinkage of protein droplets and eventual dissolution (Maharana et al., 2018). This “re-entrant phase behavior” plays an important regulatory role in modulating biomolecular phase separation in living cells (Banerjee et al., 2017). A recent work shows that during transcription, initial RNA production can enhance phase separation of transcriptional machineries such as MED1 coactivator, which further promotes transcription (Henninger et al., 2021). However, local accumulation of RNA then causes dissolution of transcriptional condensates. This RNA-mediated feedback control mechanism helps explain the origin of long-observed transcriptional bursts, a stochastic RNA production interspersed between periods of inactivity. The solubilizing effect of RNA also seems to play a role in maintaining protein homeostasis. Prion-like RNA binding proteins can undergo aberrant phase transitions from liquid droplets to solid aggregates (Lin et al., 2015; Molliex et al., 2015; Murakami et al., 2015; Patel et al., 2015), a process implicated in the onset of various neurodegenerative diseases (Mathieu et al., 2020). Disease-causing mutations are often found to impair nuclear-cytoplasmic trafficking, leading to the accumulation of mutated proteins in cytoplasmic aggregates (Dormann et al., 2010). It has been shown that a high concentration of nuclear RNA can maintain prion-like RNA binding proteins in the soluble state, preventing unwanted aberrant transitions into solid aggregates (Maharana et al., 2018).

When both solutes exhibit self-association but no heterotypic interactions, the system can have a three-phase region in addition to single or two-phase ones. Since each solute can phase-separate on its own, there are two binodals attached to each axis (Fig. 2E). If these binodals are large enough to meet at the center of the phase diagram, a three-phase coexisting region can emerge: single dilute phase and two distinct dense phases, each enriched with one type of solute but not the other, can form. This case highlights a possibility of multiple immiscible phases, and thus is particularly relevant to situations within the cell where multiple condensates with distinct compositions coexist. Stress granules and P-bodies are cytoplasmic condensates often contacting one another. A partial overlap in their molecular compositions and interaction networks provides an explanation for their immiscibility and adhesion (Sanders et al., 2020). Condensates themselves often exhibit multiple phases internally as can be seen in nucleoli, paraspeckles, PML bodies and stress granules (Jain et al., 2016; Lafontaine et al., 2021; Lallemand-Breitenbach and de Thé, 2018; Yamazaki et al., 2018). In general, relative magnitudes of interfacial energy impact spatial arrangements of multiphase droplets (Lu and Spruijt, 2020). Indeed, recent work provided evidence for the role of surface tension on the internal organization of condensates (Feric et al., 2016). Nucleoli exhibit a core-shell architecture consisting of FIB1-rich DFC core and NPM1-rich GC shell. Purified FIB1 and NPM1 proteins undergo phase separation to form liquid droplets, but when these two different types of droplets are combined, they stay as two immiscible phases. Consistent with the role of surface tension, hydrophobic FIB1 droplets are engulfed by hydrophilic NPM1 droplets, recapitulating nucleolar substructures.

The last case involves two solutes exhibiting both attractive homotypic and heterotypic interactions (Fig. 2F). The binodals intersect both axes, meaning that each solute can phase-separate on its own to form condensed phases. However, unlike the fourth case above where the three-phase region appears, attractive heterotypic interactions lead to mixing of condensed phases, resulting in a single condensed phase enriched with both types of solutes. Similar to this behavior, intracellular condensates are thought to form by multiple different scaffolds undergoing cooperative co-condensation. Examples of this case include components of transcriptional condensates (Cho et al., 2018). Typically, transcription factors consist of a folded DNA binding domain and a transactivation domain (Brent and Ptashne, 1985). The latter is often structurally disordered, and interacts with coactivators such as BRD4 and mediator complexes to facilitate transcription of target genes. Transcription factors such as Myc, p53, Nanog and Sox2 phase-separate to form in-vitro droplets (Boija et al., 2018). IDRs from coactivators such as MED1 and BRD4 also undergo phase separation in-vitro (Sabari et al., 2018). When transcription factors and coactivator IDRs are mixed together, two proteins cooperatively form a single dense phase enriched with both components (Boija et al., 2018). Mutations reducing intermolecular interactions between transcription factors and coactivators abrogate co-condensation. The same mutations also decrease transcriptional outputs of a reporter gene, implying a synergetic effect of co-condensation of transcriptional machineries on gene expression.

For multiple cases discussed above, ternary regular solution models provide qualitative description of diverse phase behaviors observed in biomolecular phase separation. Comparing phase separation in the ternary system to those in the binary one, it is clear that adding just another type of solute significantly increases the diversity and variability in phase behaviors. Notably, different classes of phase diagrams are a direct result of changes in intermolecular interactions. This implies that modulating interactions can be a regulatory knob that cells use to control phase behaviors. For example, tuning heterotypic interactions can turn crowders into clients or two immiscible droplets into a single droplet. The ternary model also provides important insights regarding the noise buffering effect of phase separation. Unlike the simple buffering effect in the binary system discussed above (Klosin et al., 2020), the degree of noise buffering in the ternary system varies significantly depending on the type of phase diagrams and the orientation of concentration variability imposed on the system (Deviri and Safran, 2021). Indeed, recent work experimentally demonstrated that the dilute-phase concentrations of intracellular condensates are often not fixed but rather increase as the expression of protein components increases (Riback et al., 2020). This behavior highlights the importance of heterotypic interactions in the formation of intracellular condensates.

PHASE BEHAVIORS IN MULTICOMPONENT SYSTEMS

Living cells contain numerous components undergoing both specific and nonspecific interactions with one another. Considering diverse types of phase diagrams appeared in the ternary system, it is not unreasonable to expect even higher levels of diversity in intracellular phase behavior. At first sight, the complexity of multicomponent phase behaviors appears intractable. However, several works have begun to bridge the gap between simple systems consisting of a few components on one hand, and multicomponent mixtures such as living cells (Jacobs and Frenkel, 2013; 2017). To model an N-component system, these works took a statistical approach, originally pioneered by Sear and Cuesta (2003), where pair-wise interactions between molecules are randomly drawn from a single Gaussian distribution. A major conclusion of the model is that the width of Gaussian distribution, i.e., variability in intermolecular interactions, is an important parameter distinguishing between two different types of phase separation: condensation and demixing. The former involves two coexisting phases with nearly identical compositions but differing in total concentration of molecules. In the latter scenario, a few species demix from others, forming phases with distinct compositions. These works suggest that tuning a few intermolecular interactions can induce the formation of demixed phases, and thus theoretically support the presence of diverse condensates in living cells.

A key question concerns how stability, miscibility and composition of individual condensates with multiple components are defined. Recent patchy-colloid based modeling suggested that the molecular connectivity in the multicomponent dense phase is a key parameter determining the stability and composition of condensates (Espinosa et al., 2020). Thus, the valency of individual components and the type of interaction nodes each component makes with others can shape together the phase diagram of specific condensates. Using stress granules and P-bodies as a model system, recent experimental work further elaborated this concept and emphasized the role of competing interaction nodes in tuning the phase miscibility (Sanders et al., 2020).

Phase separation is a physicochemical phenomenon driving a well-mixed solution into coexisting phases with distinct compositions. For the past few years, extensive studies have been carried out on diverse phase behaviors of biomolecules, and uncovered the role of phase separation in intracellular organization. The idea of phase separation has promoted our understanding of mechanisms by which cells regulate the assembly, composition and substructures of condensates. However, it must also be noted that experimental evidence for phase separation within living cells is often descriptive and lacks quantitative nature (McSwiggen et al., 2019), leaving a gap between phase behaviors observed with purified protein components and those of intracellular condensates. The model-based quantitative analysis of experimental data from living cells can help bridge this gap by providing predictive features that can be tested. This approach can then be fed back into the modeling process, which can ultimately provide important insights in the way intracellular phases behave (Riback et al., 2020).

In this regard, we discussed recent progress on biomolecular phase separation in light of regular solution models of binary and ternary systems. We showed that even a simple binary model can explain key biophysical features of intracellular condensates such as the presence of saturation concentrations and noise buffering. The ternary regular solution model highlights that heterotypic interactions provided by extra components can give rise to much richer phase behaviors including the change in saturation concentrations, re-entrant transitions and composition controls. Notably, despite being highly useful in describing the overall behaviors of many phase-separating systems, regular solution models have a limited capacity to predict quantitative details such as precise locations of phase boundary. In general, future work will be necessary to shed more light on the relation between the sequences of proteins and phase behaviors they exhibit. For example, little is known about how sequence features determine the phase immiscibility in multiphase systems. More work is also needed to elucidate how crowded intracellular environment would modulate phase separation properties. The model-based studies of biomolecular phase separation will be essential to systematically probe full complexity of intracellular organization.

This work is supported by the National Research Foundation (NRF) of Korea (NRF-2019R1C1C1006477), the National Convergence Research of Scientific Challenges (NRF-2020M3F7A1094300), and Creative-Pioneering Researchers Program through Seoul National University.

The author has no potential conflicts of interest to disclose.