Energy landscapes of homopolymeric RNAs revealed by deep unsupervised learning.

Conformational dynamics of RNA plays important roles in a variety of cellular functions such as transcriptional regulation, catalysis, scaffolding, and sensing. Recently, RNAs with low-complexity sequences have been shown to phase separate and form condensate phases similar to lowcomplexity protein domains. The affinity for phase separation and the material characteristics of RNA condensates are strongly dependent on sequence composition and patterning. We hypothesize that differences in the affinities for RNA phase separation can be uncovered by studying sequence-dependent conformational dynamics of single RNA chains. To this end, we have employed atomistic simulations and deep dimensionality reduction techniques to map temperature-dependent conformational free energy landscapes for 20 base-long homopolymeric RNA sequences: poly(U), poly(G), poly(C), and poly(A). The energy landscapes of homopolymeric RNAs reveal a plethora of metastable states with qualitatively different populations stemming from differences in base chemistry. Through detailed analysis of base, phosphate, and sugar interactions, we show that experimentally observed temperature-driven shifts in metastable state populations align with experiments on RNA phase transitions. Specifically, we find that the thermodynamics of unfolding of homopolymeric RNA follows the poly(G) > poly(A) > poly(C) > poly(U) order of stability, mirroring the propensity of RNA to form condensates. To conclude, this work shows that at least for homopolymeric RNA sequences the single-chain conformational dynamics contains sufficient information for predicting and quantifying condensate forming affinities of RNAs. Thus, we anticipate that atomically detailed studies of temeprature -dependent energy landscapes of RNAs will be a useful guide for understanding the propensity of various RNA molecules to form condensates.

Four-Dimensional Mesoscale Liquid Model of Nucleus Resolves Chromatin's Radial Organization.

Recent advances chromatin capture, imaging techniques, and polymer modeling have dramatically enhanced quantitative understanding of chromosomal folding. However, the dynamism inherent in genome architectures due to physical and biochemical forces and their impact on nuclear architecture and cellular functions remains elusive. While imaging of chromatin in four dimensions is becoming more common, there is a conspicuous lack of physics-based computational tools appropriate for revealing the forces that shape nuclear architecture and dynamics. To this end, we have developed a multiphase liquid model of the nucleus, which can resolve chromosomal territories, compartments, and nuclear lamina using a physics-based and data-informed free-energy function. The model enables rapid hypothesis-driven prototyping of nuclear dynamics in four dimensions, thereby facilitating comparison with whole nucleus imaging experiments. As an application, we model the Drosophila nucleus and map phase diagram of various possible nuclear morphologies. We shed light on the interplay of adhesive and cohesive interactions which give rise to distinct radial organization seen in conventional, inverted, and senescent nuclear architectures. The results also show the highly dynamic nature of the radial organization, the disruption of which leads to significant variability in domain coarsening dynamics and consequently variability of chromatin architecture. The model also highlights the impact of oblate nuclear geometry and heterochromatin-subtype interactions on the global chromatin architecture and local asymmetry of chromatin compartments.

Protein thermal sensing regulates physiological amyloid aggregation.

To survive, cells must respond to changing environmental conditions. One way that eukaryotic cells react to harsh stimuli is by forming physiological, RNA-seeded subnuclear condensates, termed amyloid bodies (A-bodies). The molecular constituents of A-bodies induced by different stressors vary significantly, suggesting this pathway can tailor the cellular response by selectively aggregating a subset of proteins under a given condition. Here, we identify critical structural elements that regulate heat shock-specific amyloid aggregation. Our data demonstrates that manipulating structural pockets in constituent proteins can either induce or restrict their A-body targeting at elevated temperatures. We propose a model where selective aggregation within A-bodies is mediated by the thermal stability of a protein, with temperature-sensitive structural regions acting as an intrinsic form of post-translational regulation. This system would provide cells with a rapid and stress-specific response mechanism, to tightly control physiological amyloid aggregation or other cellular stress response pathways.

Determinants of viscoelasticity and flow activation energy in biomolecular condensates.

The form and function of biomolecular condensates are intimately linked to their material properties. Here, we integrate microrheology with molecular simulations to dissect the physical determinants of condensate fluid phase dynamics. By quantifying the timescales and energetics of network relaxation in a series of heterotypic viscoelastic condensates, we uncover distinctive roles of sticker motifs, binding energy, and chain length in dictating condensate dynamical properties. We find that the mechanical relaxation times of condensate-spanning networks are determined by both intermolecular interactions and chain length. We demonstrate, however, that the energy barrier for network reconfiguration, termed flow activation energy, is independent of chain length and only varies with the strengths of intermolecular interactions. Biomolecular diffusion in the dense phase depends on a complex interplay between viscoelasticity and flow activation energy. Our results illuminate distinctive roles of chain length and sequence-specific multivalent interactions underlying the complex material and transport properties of biomolecular condensates.

Phase separation of multicomponent peptide mixtures into dehydrated clusters with hydrophilic cores.

Phase separation of biomolecules underlies the formation and regulation of various membraneless condensates in cells. How condensates function reliably while surrounded by heterogeneous and dynamic mixtures of biomolecular components with specific and nonspecific interactions is yet to be understood. Studying multicomponent biomolecular mixtures with designer peptides has recently become an attractive avenue for learning about physicochemical principles governing cellular condensates. In this work, we employed long-timescale atomistic simulations of multicomponent tripeptide mixtures with all residue substitutions to illuminate the nature of direct and water-mediated interactions in a prototypical cellular condensate environment. We find that peptide mixtures form clusters with inverse hydrophobic order. Most multivalent and charged residues are localized in the cluster's core, with a large fraction of nonaromatic hydrophobic residues remaining on the surface. This inverse hydrophobic order in peptide clusters is partly driven by the expulsion of nonspecifically bound water molecules following peptide cluster growth. The growth of clusters is also accompanied by the formation of increasing numbers of specific water-mediated interactions between polar and charged residues. While the present study focused on the condensation of short peptide motifs, the general findings and analysis techniques should be helpful for future studies on larger peptides and protein condensation.

Temperature sensitive contact modes allosterically gate TRPV3.

TRPV Ion channels are sophisticated molecular sensors designed to respond to distinct temperature thresholds. The recent surge in cryo-EM structures has provided numerous insights into the structural rearrangements accompanying their opening and closing; however, the molecular mechanisms by which TRPV channels establish precise and robust temperature sensing remain elusive. In this work we employ molecular simulations, multi-ensemble contact analysis, graph theory, and machine learning techniques to reveal the temperature-sensitive residue-residue interactions driving allostery in TRPV3. We find that groups of residues exhibiting similar temperature-dependent contact frequency profiles cluster at specific regions of the channel. The dominant mode clusters on the ankyrin repeat domain and displays a linear melting trend while others display non-linear trends. These modes describe the residue-level temperature response patterns that underlie the channel's functional dynamics. With network analysis, we find that the community structure of the channel changes with temperature. And that a network of high centrality contacts connects distant regions of the protomer to the gate, serving as a means for the temperature-sensitive contact modes to allosterically regulate channel gating. Using a random forest model, we show that the contact states of specific temperature-sensitive modes are indeed predictive of the channel gate's state. Supporting the physical validity of these modes and networks are several residues identified with our analyses that are reported in literature to be functionally critical. Our results offer high resolution insight into thermo-TRP channel function and demonstrate the utility of temperature-sensitive contact analysis.

Re-entrant transitions of locally stiff RNA chains in the presence of polycations leads to gelated architectures.

The liquid-liquid phase separation of protein and nucleic acid mixtures drives the formation of numerous membraneless compartments in cells. Temperature variation is commonly used for mapping condensate phase diagrams, which often display unique upper critical temperatures. Recent report on peptide-RNA mixtures has shown the existence of lower and upper critical solution temperatures, highlighting the importance of temperature-dependent solvent and ion-mediated forces. In the present work, we employ residue-level coarse-grained models of RNA and polycation peptide chains for simulating temperature-induced re-entrant transitions and shedding light on the role played by mobile ions, temperature-dependent dielectric permittivity, and local chain stiffness. We show that differences in bending rigidity can significantly modulate condensate topology leading to the formation of gelated or fibril like architectures. The study also finds that temperature dependence of water permittivity is generally sufficient for recapitulating experimentally observed closed loop and LCST phase diagrams of highly charged protein-RNA mixtures. However, we find that similar-looking closed-loop phase diagrams can correspond to vastly different condensate topologies.

Solution Structure Ensembles of the Open and Closed Forms of the ∼130 kDa Enzyme I via AlphaFold Modeling, Coarse Grained Simulations, and NMR.

Large-scale interdomain rearrangements are essential to protein function, governing the activity of large enzymes and molecular machineries. Yet, obtaining an atomic-resolution understanding of how the relative domain positioning is affected by external stimuli is a hard task in modern structural biology. Here, we show that combining structural modeling by AlphaFold2 with coarse-grained molecular dynamics simulations and NMR residual dipolar coupling data is sufficient to characterize the spatial domain organization of bacterial enzyme I (EI), a ∼130 kDa multidomain oligomeric protein that undergoes large-scale conformational changes during its catalytic cycle. In particular, we solve conformational ensembles for EI at two different experimental temperatures and demonstrate that a lower temperature favors sampling of the catalytically competent closed state of the enzyme. These results suggest a role for conformational entropy in the activation of EI and demonstrate the ability of our protocol to detect and characterize the effect of external stimuli (such as mutations, ligand binding, and post-translational modifications) on the interdomain organization of multidomain proteins. We expect the ensemble refinement protocol described here to be easily transferrable to the investigation of the structure and dynamics of other uncharted multidomain systems and have assembled a Google Colab page (https://potoyangroup.github.io/Seq2Ensemble/) to facilitate implementation of the presented methodology elsewhere.

Protein Conformational Dynamics Underlie Selective Recognition of Thermophilic over Mesophilic Enzyme I by a Substrate Analogue.

Substrate selectivity is an important preventive measure to decrease the possibility of cross interactions between enzymes and metabolites that share structural similarities. In addition, understanding the mechanisms that determine selectivity towards a particular substrate increases the knowledge base for designing specific inhibitors for target enzymes. Here, we combine NMR, molecular dynamics (MD) simulations, and protein engineering to investigate how two substrate analogues, allylicphosphonate (cPEP) and sulfoenolpyruvate (SEP), recognize the mesophilic (eEIC) and thermophilic (tEIC) homologues of the receptor domain of bacterial Enzyme I, which has been proposed as a target for antimicrobial research. Chemical Shift Perturbation (CSP) experiments show that cPEP and SEP recognize tEIC over the mesophilic homologue. Combined Principal Component Analysis of half-microsecond-long MD simulations reveals that incomplete quenching of a breathing motion in the eEIC-ligand complex destabilizes the interaction and makes the investigated substrate analogues selective toward the thermophilic enzyme. Our results indicate that residual protein motions need to be considered carefully when optimizing small molecule inhibitors of EI. In general, our work demonstrates that protein conformational dynamics can be exploited in the rational design and optimization of inhibitors with subfamily selectivity.

Temperature-Sensitive Contact Modes Allosterically Gate TRPV3.

TRPV Ion channels are sophisticated molecular sensors designed to respond to distinct temperature thresholds. The recent surge in cryo-EM structures has provided numerous insights into the structural rearrangements accompanying their opening and closing; however, the molecular mechanisms by which TRPV channels establish precise and robust temperature sensing remain elusive. In this work we employ molecular simulations, multi-ensemble contact analysis, graph theory, and machine learning techniques to reveal the temperature-sensitive residue-residue interactions driving allostery in TRPV3. We find that groups of residues exhibiting similar temperature-dependent contact frequency profiles cluster at specific regions of the channel. The dominant mode clusters on the ankyrin repeat domain and displays a linear melting trend while others display non-linear trends. These modes describe the residue-level temperature response patterns that underlie the channel's functional dynamics. With network analysis, we find that the community structure of the channel changes with temperature. And that a network of high centrality contacts connects distant regions of the protomer to the gate, serving as a means for the temperature-sensitive contact modes to allosterically regulate channel gating. Using a random forest model, we show that the contact states of specific temperature-sensitive modes are indeed predictive of the channel gate's state. Supporting the physical validity of these modes and networks are several residues identified with our analyses that are reported in literature to be functionally critical. Our results offer high resolution insight into thermo-TRP channel function and demonstrate the utility of temperature-sensitive contact analysis.

Multiscale Modeling of Protein-RNA Condensation in and Out of Equilibrium.

A vast number of intracellular membraneless bodies also known as biomolecular condensates form through a liquid-liquid phase separation (LLPS) of biomolecules. To date, phase separation has been identified as the main driving force for a membraneless organelles such as nucleoli, Cajal bodies, stress granules, and chromatin compartments. Recently, the protein-RNA condensation is receiving increased attention, because it is closely related to the biological function of cells such as transcription, translation, and RNA metabolism. Despite the multidisciplinary efforts put forth to study the biophysical properties of protein-RNA condensates, there are many fundamental unanswered questions regarding the mechanism of formation and regulation of protein-RNA condensates in eukaryotic cells. Major challenges in studying protein-RNA condensation stem from (i) the molecular heterogeneity and conformational flexibility of RNA and protein chains and (ii) the nonequilibrium nature of transcription and cellular environment. Computer simulations, bioinformatics, and mathematical models are uniquely positioned for shedding light on the microscopic nature of protein-RNA phase separation. To this end, there is an urgent need for innovative models with the right spatiotemporal resolution for confronting the experimental observables in a comprehensive and physics-based manner. In this chapter, we will summarize the currently emerging research efforts, which employ atomistic and coarse-grained molecular models and field theoretical models to understand equilibrium and nonequilibrium aspects of protein-RNA condensation.

Temperature-sensitive contacts in disordered loops tune enzyme I activity.

Homologous enzymes with identical folds often exhibit different thermal and kinetic behaviors. Understanding how an enzyme sequence encodes catalytic activity at functionally optimal temperatures is a fundamental problem in biophysics. Recently it was shown that the residues that tune catalytic activities of thermophilic/mesophilic variants of the C-terminal domain of bacterial enzyme I (EIC) are largely localized within disordered loops, offering a model system with which to investigate this phenomenon. In this work, we use molecular dynamics simulations and mutagenesis experiments to reveal a mechanism of sequence-dependent activity tuning of EIC homologs. We find that a network of contacts in the catalytic loops is particularly sensitive to changes in temperature, with some contacts exhibiting distinct linear or nonlinear temperature-dependent trends. Moreover, these trends define structurally clustered dynamical modes and can distinguish regions that tend toward order or disorder at higher temperatures. Assaying several thermophilic EIC mutants, we show that complementary mesophilic mutations to the most temperature-sensitive positions exhibit the most enhanced activity, while mutations to relatively temperature insensitive positions exhibit the least enhanced activities. These results provide a mechanistic explanation of sequence-dependent temperature tuning and offer a computational method for rational enzyme modification.

BIAPSS: A Comprehensive Physicochemical Analyzer of Proteins Undergoing Liquid-Liquid Phase Separation.

The liquid-liquid phase separation (LLPS) of biomolecules is a phenomenon which is nowadays recognized as the driving force for the biogenesis of numerous functional membraneless organelles and cellular bodies. The interplay between the protein primary sequence and phase separation remains poorly understood, despite intensive research. To uncover the sequence-encoded signals of protein capable of undergoing LLPS, we developed a novel web platform named BIAPSS (Bioinformatics Analysis of LLPS Sequences). This web server provides on-the-fly analysis, visualization, and interpretation of the physicochemical and structural features for the superset of curated LLPS proteins.

RNA chain length and stoichiometry govern surface tension and stability of protein-RNA condensates.

Proteomic studies have shown that cellular condensates are frequently enriched in diverse RNA molecules, which is suggestive of mechanistic links between phase separation and transcriptional activities. Here, we report a systematic experimental and computational study of thermodynamic landscapes and interfacial properties of protein-RNA condensates. We have studied the affinity of protein-RNA condensation as a function of variable RNA sequence length and RNA-protein stoichiometry under different ionic environments and external crowding. We have chosen the PolyU sequences for RNA and arginine/glycine-rich intrinsically disordered peptide (RGG) for proteins as a model system of RNA-protein condensates, which we then investigate through in vitro microscopy measurements and coarse-grained molecular dynamics simulations. We find that crowding and RNA chain length can have a major stabilizing effect on the condensation. We also find that the RNA-protein charge ratio is a crucial variable controlling stability, interfacial properties, and the reentrant phase behavior of RGG-RNA mixtures.

Programmable viscoelasticity in protein-RNA condensates with disordered sticker-spacer polypeptides.

Liquid-liquid phase separation of multivalent proteins and RNAs drives the formation of biomolecular condensates that facilitate membrane-free compartmentalization of subcellular processes. With recent advances, it is becoming increasingly clear that biomolecular condensates are network fluids with time-dependent material properties. Here, employing microrheology with optical tweezers, we reveal molecular determinants that govern the viscoelastic behavior of condensates formed by multivalent Arg/Gly-rich sticker-spacer polypeptides and RNA. These condensates behave as Maxwell fluids with an elastically-dominant rheological response at shorter timescales and a liquid-like behavior at longer timescales. The viscous and elastic regimes of these condensates can be tuned by the polypeptide and RNA sequences as well as their mixture compositions. Our results establish a quantitative link between the sequence- and structure-encoded biomolecular interactions at the microscopic scale and the rheological properties of the resulting condensates at the mesoscale, enabling a route to systematically probe and rationally engineer biomolecular condensates with programmable mechanics.

The interplay of chromatin phase separation and lamina interactions in nuclear organization.

The genetic material of eukaryotes is segregated into transcriptionally active euchromatin and silent heterochromatin compartments. The spatial arrangement of chromatin compartments evolves over the course of cellular life in a process that remains poorly understood. The latest nuclear imaging experiments reveal a number of dynamical signatures of chromatin that are reminiscent of active multiphase liquids. This includes the observations of viscoelastic response, coherent motions, Ostwald ripening, and coalescence of chromatin compartments. There is also growing evidence that liquid-liquid phase separation of protein and nucleic acid components is the underlying mechanism for the dynamical behavior of chromatin. To dissect the organizational and dynamical implications of chromatin's liquid behavior, we have devised a phenomenological field-theoretic model of the nucleus as a multiphase condensate of liquid chromatin types. Employing the liquid chromatin model of the Drosophila nucleus, we have carried out an extensive set of simulations with an objective to shed light on the dynamics and chromatin patterning observed in the latest nuclear imaging experiments. Our simulations reveal the emergence of experimentally detected mesoscale chromatin channels and spheroidal droplets which arise from the dynamic interplay of chromatin type to type interactions and intermingling of chromosomal territories. We also quantitatively reproduce coherent motions of chromatin domains observed in displacement correlation spectroscopy measurements which are explained within the framework of our model by phase separation of chromatin types operating within constrained intrachromosomal and interchromosomal boundaries. Finally, we illuminate the role of heterochromatin-lamina interactions in the nuclear organization by showing that these interactions enhance the mobility of euchromatin and indirectly introduce correlated motions of heterochromatin droplets.

Solvent Exposure and Ionic Condensation Drive Fuzzy Dimerization of Disordered Heterochromatin Protein Sequence.

Proteins with low complexity, disordered sequences are receiving increasing attention due to their central roles in the biogenesis and regulation of membraneless organelles. In eukaryotic organisms, a substantial fraction of disordered proteins reside in the nucleus, thereby facilitating the formation of nuclear bodies, nucleolus, and chromatin compartmentalization. The heterochromatin family of proteins (HP1) is an important player in driving the formation of gene silenced mesoscopic heterochromatin B compartments and pericentric regions. Recent experiments have shown that the HP1a sequence of Drosophila melanogaster can undergo liquid-liquid phase separation under both in vitro and in vivo conditions, induced by changes of the monovalent salt concentration. While the phase separation of HP1a is thought to be the mechanism underlying chromatin compartmentalization, the molecular level mechanistic picture of salt-driven phase separation of HP1a has remained poorly understood. The disordered hinge region of HP1a is seen as the driver of salt-induced condensation because of its charge enriched sequence and post-translational modifications. Here, we set out to decipher the mechanisms of salt-induced condensation of HP1a through a systematic study of salt-dependent conformations of single chains and fuzzy dimers of disordered HP1a hinge sequences. Using multiple independent all-atom simulations with and without enhanced sampling, we carry out detailed characterization of conformational ensembles of disordered HP1a chains under different ionic conditions using various polymeric and structural measures. We show that the mobile ion release, enhancement of local transient secondary structural elements, and side-chain exposure to solvent are robust trends that accompany fuzzy dimer formation. Furthermore, we find that salt-induced changes in the ensemble of conformations of HP1a disordered hinge sequence fine-tune the inter-chain vs. self-chain interactions in ways that favor fuzzy dimer formation under low salt conditions in the agreement with condensation trends seen in experiments.

Sequence-encoded and composition-dependent protein-RNA interactions control multiphasic condensate morphologies.

Multivalent protein-protein and protein-RNA interactions are the drivers of biological phase separation. Biomolecular condensates typically contain a dense network of multiple proteins and RNAs, and their competing molecular interactions play key roles in regulating the condensate composition and structure. Employing a ternary system comprising of a prion-like polypeptide (PLP), arginine-rich polypeptide (RRP), and RNA, we show that competition between the PLP and RNA for a single shared partner, the RRP, leads to RNA-induced demixing of PLP-RRP condensates into stable coexisting phases-homotypic PLP condensates and heterotypic RRP-RNA condensates. The morphology of these biphasic condensates (non-engulfing/ partial engulfing/ complete engulfing) is determined by the RNA-to-RRP stoichiometry and the hierarchy of intermolecular interactions, providing a glimpse of the broad range of multiphasic patterns that are accessible to these condensates. Our findings provide a minimal set of physical rules that govern the composition and spatial organization of multicomponent and multiphasic biomolecular condensates.

Liquid-liquid phase separation driven compartmentalization of reactive nucleoplasm.

The nucleus of eukaryotic cells harbors active and out of equilibrium environments conducive to diverse gene regulatory processes. On a molecular scale, gene regulatory processes take place within hierarchically compartmentalized sub-nuclear bodies. While the impact of nuclear structure on gene regulation is widely appreciated, it has remained much less clear whether and how gene regulation is impacting nuclear order itself. Recently, the liquid-liquid phase separation emerged as a fundamental mechanism driving the formation of biomolecular condensates, including membrane-less organelles, chromatin territories, and transcriptional domains. The transience and environmental sensitivity of biomolecular condensation are strongly suggestive of kinetic gene-regulatory control of phase separation. To better understand kinetic aspects controlling biomolecular phase-separation, we have constructed a minimalist model of the reactive nucleoplasm. The model is based on the Cahn-Hilliard formulation of ternary protein-RNA-nucleoplasm components coupled to non-equilibrium and spatially dependent gene expression. We find a broad range of kinetic regimes through an extensive set of simulations where the interplay of phase separation and reactive timescales can generate heterogeneous multi-modal gene expression patterns. Furthermore, the significance of this finding is that heterogeneity of gene expression is linked directly with the heterogeneity of length-scales in phase-separated condensates.

Hybrid Thermophilic/Mesophilic Enzymes Reveal a Role for Conformational Disorder in Regulation of Bacterial Enzyme I.

Conformational disorder is emerging as an important feature of biopolymers, regulating a vast array of cellular functions, including signaling, phase separation, and enzyme catalysis. Here we combine NMR, crystallography, computer simulations, protein engineering, and functional assays to investigate the role played by conformational heterogeneity in determining the activity of the C-terminal domain of bacterial Enzyme I (EIC). In particular, we design chimeric proteins by hybridizing EIC from thermophilic and mesophilic organisms, and we characterize the resulting constructs for structure, dynamics, and biological function. We show that EIC exists as a mixture of active and inactive conformations and that functional regulation is achieved by tuning the thermodynamic balance between active and inactive states. Interestingly, we also present a hybrid thermophilic/mesophilic enzyme that is thermostable and more active than the wild-type thermophilic enzyme, suggesting that hybridizing thermophilic and mesophilic proteins is a valid strategy to engineer thermostable enzymes with significant low-temperature activity.