Molecular-crowding effects on single-molecule RNA folding/unfolding thermodynamics and kinetics

Macromolecular crowding effects on protein dynamics

Macromolecular crowding experiments bridge the gap between in-vivo and in-vitro studies by mimicking some of the cellular complexities like high viscosity and limited space, while still manageable for experiments and analysis. Macromolecular crowding impacts all biological processes and is a focus of contemporary research. Recent reviews have highlighted the effect of crowding on various protein properties. One of the essential characteristics of protein is its dynamic nature; however, how protein dynamics get modulated in the crowded milieu has been largely ignored. This article discusses how protein translational, rotational, conformational, and solvation dynamics change under crowded conditions, summarizing key observations in the literature. We emphasize our research on microsecond conformational and water dynamics in crowded milieus and their impact on enzymatic activity and stability. Lastly, we provided our outlook on how this field might move forward in the future.

Role of Associated Water Dynamics on Protein Stability and Activity in Crowded Milieu

Macromolecular crowding bridges in vivo and in vitro studies by simulating cellular complexities such as high viscosity and limited space while maintaining the experimental feasibility. Over the last two decades, the impact of macromolecular crowding on protein stability and activity has been a significant topic of study and discussion, though still lacking a thorough mechanistic understanding. This article investigates the role of associated water dynamics on protein stability and activity within crowded environments, using bromelain and Ficoll-70 as the model systems. Traditional crowding theory primarily attributes protein stability to entropic effects (excluded volume) and enthalpic interactions. However, our recent findings suggest that water structure modulation plays a crucial role in a crowded environment. In this report, we strengthen the conclusion of our previous study, i.e., rigid-associated water stabilizes proteins via entropy and destabilizes them via enthalpy, while flexible water has the opposite effect. In the process, we addressed previous shortcomings with a systematic concentration-dependent study using a single-domain protein and component analysis of solvation dynamics. More importantly, we analyze bromelain's hydrolytic activity using the Michaelis-Menten model to understand kinetic parameters like maximum velocity (Vmax) achieved by the system and the Michaelis-Menten coefficient (KM). Results indicate that microviscosity (not the bulk viscosity) controls the enzyme-substrate (ES) complex formation, where an increase in the microviscosity makes the ES complex formation less favorable. On the other hand, flexible associated water dynamics were found to favor the rate of product formation significantly from the ES complex, while rigid associated water hinders it. This study improves our understanding of protein stability and activity in crowded environments, highlighting the critical role of associated water dynamics.

Nonequilibrium relaxation of soft responsive colloids

Stimuli-responsive macromolecules display large conformational changes during their dynamics, sometimes switching between states. Such a multi-stability is useful for the development of soft functional materials. Here, we introduce a mean-field dynamical density functional theory for a model of responsive colloids to study the nonequilibrium dynamics of a colloidal dispersion in time-dependent external fields, with a focus on the coupling of translational and conformational dynamics during their relaxation. Specifically, we consider soft Gaussian particles with a bimodal size distribution between two confining walls with time-dependent (switching-on and off) external gravitational and osmotic fields. We find a rich relaxation behavior of the systems in excellent agreement with particle-based Brownian dynamics computer simulations. In particular, we find time-asymmetric relaxations of integrated observables (wall pressures, mean size, and liquid center-of-mass) for activation/deactivation of external potentials, respectively, which are tunable by the ratio of translational and conformational diffusion time scales. Our work thus paves the way for studying the nonequilibrium relaxation dynamics of complex soft matter with multiple degrees of freedom and hierarchical relaxations.

Macromolecular crowding has opposite effects on two critical sub-steps of transcription initiation

Transcription initiation, the first step in gene expression, has been studied extensively in dilute buffer, a condition which fails to consider the crowded environment in live cells. Recent reports indicate the kinetics of promoter escape is altered in crowded conditions for a consensus bacterial promoter. Here, we use a real-time fluorescence enhancement assay to study the kinetics of unwound bubble formation and promoter escape for three separate promoters. We find that the effect of crowding on transcription initiation is complex, with lower rates of unwound bubble formation, higher rates of promoter escape, and large variations depending on promoter identity. Based on our results, we suggest that altered conditions of crowding inside a live cell can trigger global changes.

The disordered N-terminal tail of SARS-CoV-2 Nucleocapsid protein forms a dynamic complex with RNA

The SARS-CoV-2 Nucleocapsid (N) protein is responsible for condensation of the viral genome. Characterizing the mechanisms controlling nucleic acid binding is a key step in understanding how condensation is realized. Here, we focus on the role of the RNA binding domain (RBD) and its flanking disordered N-terminal domain (NTD) tail, using single-molecule Förster Resonance Energy Transfer and coarse-grained simulations. We quantified contact site size and binding affinity for nucleic acids and concomitant conformational changes occurring in the disordered region. We found that the disordered NTD increases the affinity of the RBD for RNA by about 50-fold. Binding of both nonspecific and specific RNA results in a modulation of the tail configurations, which respond in an RNA length-dependent manner. Not only does the disordered NTD increase affinity for RNA, but mutations that occur in the Omicron variant modulate the interactions, indicating a functional role of the disordered tail. Finally, we found that the NTD-RBD preferentially interacts with single-stranded RNA and that the resulting protein:RNA complexes are flexible and dynamic. We speculate that this mechanism of interaction enables the Nucleocapsid protein to search the viral genome for and bind to high-affinity motifs.

Modulation of Aβ 16-22 aggregation by glucose

The self-assembly of amyloid-beta (Aβ) peptides into fibrillar structures in the brain is a signature of Alzheimer's disease. Recent studies have reported correlations between Alzheimer's disease and type-2 diabetes. Structurally, hyperglycemia induces covalent protein crosslinkings by advanced glycation end products (AGE), which can affect the stability of Aβ oligomers. In this work, we leverage physics-based coarse-grained molecular simulations to probe alternate thermodynamic pathways that affect peptide aggregation propensities at varying concentrations of glucose molecules. Similar to previous experimental reports, our simulations show a glucose concentration-dependent increase in Aβ aggregation rates, without changes in the overall secondary structure content. We discovered that glucose molecules prefer partitioning onto the aggregate-water interface at a specific orientation, resulting in a loss of molecular rotational entropy. This effectively hastens the aggregation rates, as peptide self-assembly can reduce the available surface area for peptide-glucose interactions. This work introduces a new thermodynamic-driven pathway, beyond chemical cross-linking, that can modulate Aβ aggregation.

Molecular crowding effect in Hantzch pyridine synthesis in polyethylene glycol aqueous solution

In a molecular crowding environment, the kinetics and thermodynamics differ from those in a diluted solution. Although the molecular crowding effect has been extensively investigated, its fundamental kinetics and thermodynamics remain unclear. In this study, we investigated the change in the rate constant (k) of the Hantzch pyridine reaction in a molecular crowding environment using polyethylene glycol (PEG). While the k value increased to a PEG concentration (CPEG) of 10 vol%, a decreasing trend was observed for CPEG > 20 vol%. This intriguing behavior was analyzed based on the increase in reactant activity due to volume exclusion and the decrease in water activity due to osmotic pressure. Volume exclusion and osmotic pressure had opposing effects on the reaction, which were positive for volume exclusion and negative for osmotic pressure. We found that k decreased when the negative effect of the osmotic pressure surpassed the volume exclusion effect.

The 3' terminal region of Zika virus RNA contains a conserved G-quadruplex and is unfolded by human DDX17

Zika virus (ZIKV) infection remains a worldwide concern, and currently no effective treatments or vaccines are available. Novel therapeutics are an avenue of interest that could probe viral RNA-human protein communication to stop viral replication. One specific RNA structure, G-quadruplexes (G4s), possess various roles in viruses and all domains of life, including transcription and translation regulation and genome stability, and serves as nucleation points for RNA liquid-liquid phase separation. Previous G4 studies on ZIKV using a quadruplex forming G-rich sequences Mapper located a potential G-quadruplex sequence in the 3' terminal region (TR) and was validated structurally using a 25-mer oligo. It is currently unknown if this structure is conserved and maintained in a large ZIKV RNA transcript and its specific roles in viral replication. Using bioinformatic analysis and biochemical assays, we demonstrate that the ZIKV 3' TR G4 is conserved across all ZIKV isolates and maintains its structure in a 3' TR full-length transcript. We further established the G4 formation using pyridostatin and the BG4 G4-recognizing antibody binding assays. Our study also demonstrates that the human DEAD-box helicases, DDX3X132-607 and DDX17135-555, bind to the 3' TR and that DDX17135-555 unfolds the G4 present in the 3' TR. These findings provide a path forward in potential therapeutic targeting of DDX3X or DDX17's binding to the 3' TR G4 region for novel treatments against ZIKV.

External field-driven property localization in liquids of responsive macromolecules

We explore theoretically the effects of external potentials on the spatial distribution of particle properties in a liquid of explicitly responsive macromolecules. In particular, we focus on the bistable particle size as a coarse-grained internal degree of freedom (DoF, or "property"), σ, that moves in a bimodal energy landscape, in order to model the response of a state-switching (big-to-small) macromolecular liquid to external stimuli. We employ a mean-field density functional theory (DFT) that provides the full inhomogeneous equilibrium distributions of a one-component model system of responsive colloids (RCs) interacting with a Gaussian pair potential. For systems confined between two parallel hard walls, we observe and rationalize a significant localization of the big particle state close to the walls, with pressures described by an exact RC wall theorem. Application of more complex external potentials, such as linear (gravitational), osmotic, and Hamaker potentials, promotes even stronger particle size segregation, in which macromolecules of different size are localized in different spatial regions. Importantly, we demonstrate how the degree of responsiveness of the particle size and its coupling to the external potential tune the position-dependent size distribution. The DFT predictions are corroborated by Brownian dynamics simulations. Our study highlights the fact that particle responsiveness can be used to localize liquid properties and therefore helps to control the property- and position-dependent function of macromolecules, e.g., in biomedical applications.

Ionic liquid-caged nucleic acids enable active folding-based molecular recognition with hydrolysis resistance

Beyond storage and transmission of genetic information in cellular life, nucleic acids can perform diverse interesting functions, including specific target recognition and biochemical reaction acceleration; the versatile biopolymers, however, are acutely vulnerable to hydrolysis-driven degradation. Here, we demonstrate that the cage effect of choline dihydrogen phosphate permits active folding of nucleic acids like water, but prevents their phosphodiester hydrolysis unlike water. The choline-based ionic liquid not only serves as a universal inhibitor of nucleases, exceptionally extending half-lives of nucleic acids up to 6 500 000 times, but highly useful tasks of nucleic acids (e.g. mRNA detection of molecular beacons, ligand recognition of aptamers, and transesterification reaction of ribozymes) can be also conducted with well-conserved affinities and specificities. As liberated from the function loss and degradation risk, the presence of undesired and unknown nucleases does not undermine desired molecular functions of nucleic acids without hydrolysis artifacts even in nuclease cocktails and human saliva.

Dynamics of G-Quadruplex Formation under Molecular Crowding

G-quadruplex (G4) structures assemble from guanine-rich DNA sequences and are believed to regulate several key cellular processes. G4 formation and conformational interconversions are well-established to occur dynamically in vitro. However, a clear understanding of G4 formation dynamics in cells as well as under conditions mimicking the cellular environment is missing. To fill this gap, we have investigated the G4 dynamics in molecularly crowded solutions, thus replicating the effect of the excluded volume present in cells. The results show that the volume exclusion exerted by large crowding agents accelerates the rate of G4 formation by at least an order of magnitude, leading to significant G4 stabilization. Extrapolation from our experimental data predicts crowding-induced G4 stabilization by more than 3 kcal/mol, under crowding levels found in the cellular environment. Such effects are likely to be important for G4-driven regulatory functions.

Molecular Crowding Facilitates Ribozyme-Catalyzed RNA Assembly

Catalytic RNAs or ribozymes are considered to be central to primordial biology. Most ribozymes require moderate to high concentrations of divalent cations such as Mg2+ to fold into their catalytically competent structures and perform catalysis. However, undesirable effects of Mg2+ such as hydrolysis of reactive RNA building blocks and degradation of RNA structures are likely to undermine its beneficial roles in ribozyme catalysis. Further, prebiotic cell-like compartments bounded by fatty acid membranes are destabilized in the presence of Mg2+, making ribozyme function inside prebiotically relevant protocells a significant challenge. Therefore, we sought to identify conditions that would enable ribozymes to retain activity at low concentrations of Mg2+. Inspired by the ability of ribozymes to function inside crowded cellular environments with <1 mM free Mg2+, we tested molecular crowding as a potential mechanism to lower the Mg2+ concentration required for ribozyme-catalyzed RNA assembly. Here, we show that the ribozyme-catalyzed ligation of phosphorimidazolide RNA substrates is significantly enhanced in the presence of the artificial crowding agent polyethylene glycol. We also found that molecular crowding preserves ligase activity under denaturing conditions such as alkaline pH and the presence of urea. Additionally, we show that crowding-induced stimulation of RNA-catalyzed RNA assembly is not limited to phosphorimidazolide ligation but extends to the RNA-catalyzed polymerization of nucleoside triphosphates. RNA-catalyzed RNA ligation is also stimulated by the presence of prebiotically relevant small molecules such as ethylene glycol, ribose, and amino acids, consistent with a role for molecular crowding in primordial ribozyme function and more generally in the emergence of RNA-based cellular life.

Excluded Volume and Weak Interactions in Crowded Solutions Modulate Conformations and RNA Binding of an Intrinsically Disordered Tail

The cellular milieu is a solution crowded with a significant concentration of different components (proteins, nucleic acids, metabolites, etc.). Such a crowded environment affects protein conformations, dynamics, and interactions. Intrinsically disordered proteins and regions are particularly sensitive to these effects. Here, we investigate the impact on an intrinsically disordered tail that flanks a folded domain, the N-terminal domain, and the RNA-binding domain of the SARS-CoV-2 nucleocapsid protein. We mimic the crowded environment of the cell using polyethylene glycol (PEG) and study its impact on protein conformations using single-molecule Förster resonance energy transfer. We found that high-molecular-weight PEG induces a collapse of the disordered N-terminal tail, whereas low-molecular-weight PEG induces a chain expansion. Our data can be explained by accounting for two opposing contributions: favorable interactions between the protein and crowder molecules and screening of excluded volume interactions. We further characterized the interaction between protein and RNA in the presence of crowding agents. While for all PEG molecules tested, we observed an increase in the binding affinity, the trend is not monotonic as a function of the degree of PEG polymerization. This points to the role of nonspecific protein-PEG interactions on binding in addition to the entropic effects due to crowding. To separate the enthalpic and entropic components of the effects, we investigated the temperature dependence of the association constants in the absence and presence of crowders. Finally, we compared the effects of crowding across mutations in the disordered region and found that the threefold difference in association constants for two naturally occurring variants of the SARS-CoV-2 nucleocapsid protein is reduced to almost identical affinities in the presence of crowders. Overall, our data provide new insights into understanding and modeling the contribution of crowding effects on disordered regions, including the impact of interactions between proteins and crowders and their interplay when binding a ligand.

Quantification of macromolecule crowding at single-molecule level

Macromolecule crowding has a prominent impact on a series of biochemical processes in the cell. It is also expected to promote macromolecular complexation induced by excluded volume effects, which conflicts with recent advances in the thermodynamic interaction between inert, synthetic polymers, and nucleic acids. Along this line, a method combining high-resolution magnetic tweezers and extended crowder-oxDNA model was applied to resolve these discrepancies by systematically studying the kinetics and thermodynamics of the folding-unfolding transition for an individual DNA hairpin in a crowded environment. More specifically, from the magnetic tweezers-based experiments, the linear dependence of the critical force of the DNA hairpin on the polyethylene glycol (PEG) concentration was demonstrated, which is consistent with the results based on the crowder-oxDNA model in which the Lennard-Jones potential was adopted to express the interaction between the crowders and the DNA hairpin. These consistencies highlight that the excluded volume effects are mainly responsible for the interaction between PEG and the DNA hairpin, which is different from the interaction between dextran and the DNA hairpin. In the meantime, the dependence of the folding rate on the molecule weight of PEG, which was different from fluorescence resonance energy transfer-based results, was identified. The proposed method opens a path to detect the interaction between an inert, synthetic molecule, and the DNA hairpin, which is important to accurately mimic the cytosolic environments using mixtures of different inert molecules.

Crowded environments tune the fold-switching in metamorphic proteins

Metamorphic proteins such as circadian clock protein KaiB and human chemokine XCL1 play vital roles in regulating biological processes, including gene expression, circadian clock and innate immune responses, and perform distinct functions in living cell by switching different structures in response to cellular environment stimuli. However, it is unclear how complex and crowded intracellular environments affect conformational rearrangement of metamorphic proteins. Here, the kinetics and thermodynamics of two well-characterized metamorphic proteins, circadian clock protein KaiB and human chemokine XCL1, were quantified in physiologically relevant environments by using NMR spectroscopy, indicating that crowded agents shift equilibrium towards the inactive form (ground-state KaiB and Ltn10-like state XCL1) without disturbing the corresponding structures, and crowded agents have predominantly impact on the exchange rate of XCL1 that switches folds on timescales of seconds, but have slightly impact on the exchange rate of KaiB that switches folds on timescales of hours. Our data shed light on how metamorphic proteins can respond immediately to the changed crowded intracellular conditions that induced by environmental cues and then execute different functions in living cell, and it also enhances our understanding of how environments enrich the sequence-structure-function paradigm.

Macromolecular crowding impact on anti-CRISPR AcrIIC3/NmeCas9 complex: Insights from scaled particle theory, molecular dynamics, and elastic networks models

The coupling of Cas9 and its inhibitor AcrIIC3, both from the bacterium Neisseria meningitidis (Nme), form a homodimer of the (NmeCas9/AcrIIC3)₂ type. This coupling was studied to assess the impact of their interaction with the crowders in the following environments: (1) homogeneous crowded, (2) heterogeneous, and (3) microheterogeneous cytoplasmic. For this, statistical thermodynamic models based on the scaled particle theory (SPT) were used, considering the attractive and repulsive protein-crowders contributions and the stability of the formation of spherocylindrical homodimers and the effects of changes in the size of spherical dimers were estimated. Studies based on models of dynamics, elastic networks, and statistical potentials to the formation of complexes NmeCas9/AcrIIC3 using PEG as the crowding agent support the predictions from SPT. Macromolecular crowding stabilizes the formation of the dimers, being more significant when the attractive protein-crowder interactions are weaker and the crowders are smaller. The coupling is favored towards the formation of spherical and compact dimers due to crowding addition (excluded-volume effects) and the thermodynamic stability of the dimers is markedly dependent on the size of the crowders. These results support the experimental mechanistic proposal of inhibition of NmeCas9 mediated by AcrIIC3.

Liquid structure of bistable responsive macromolecules using mean-field density-functional theory

Macromolecular crowding typically applies to biomolecular and polymer-based systems in which the individual particles often feature a two-state folded/unfolded or coil-to-globule transition, such as found for proteins and peptides, DNA and RNA, or supramolecular polymers. Here, we employ a mean-field density functional theory (DFT) of a model of soft and bistable responsive colloids (RCs) in which the size of the macromolecule is explicitly resolved as a degree of freedom living in a bimodal 'Landau' energy landscape (exhibiting big and small states), thus directly responding to the crowding environment. Using this RC-DFT we study the effects of self-crowding on the liquid bulk structure and thermodynamics for different energy barriers and softnesses of the bimodal energy landscape, in conditions close to the coil-to-globule transition. We find substantial crowding effects on the internal distributions, a complex polydispersity behavior, and quasi-universal compression curves for increasing (generalized) packing fractions. Moreover, we uncover distinct signatures of bimodal versus unimodal behavior in the particle compression. Finally, the analysis of the pair structure - derived from the test particle route - reveals that the microstructure of the liquid is quite inhomogeneous due to local depletion effects, tuneable by particle softness.

Stability of DNA and RNA hairpins: a comparative study based on ox-DNA

Advances in single-molecule experiments on macromolecular crowding urgently need an efficient simulation method to resolve their discrepancies quantitatively. Ox-DNA model has been since reworked to treat the thermodynamics and mechanical properties of DNA/RNA hairpin at a stretching force. In hopping experiments, the critical forces of RNA hairpins at different temperatures are greater than those of DNA hairpins, in addition, the Gibbs free energy at a fixed temperature required to convert an RNA hairpin into a single-stranded molecule at zero force is obviously greater than that of DNA hairpin and gradually decreases by increasing the temperature. As far as force-ramping experiments are concerned, the first-rupture forces of RNA/DNA hairpins corresponding to the maximum probability density linearly pertain to the force-loading rate, with those of RNA hairpins being greater. The extended ox-DNA model could potentially identify the interaction between biologically inert polymer and RNA/DNA hairpins in crowded environments.

GB1 Dimerization in Crowders: A Multiple Resolution Approach

In-cell protein-protein association, which is crucial in enzyme catalysis and polymerization, occurs in an environment that is highly heterogeneous and crowded. The crowder molecules exclude the reactant molecules from occupying certain regions of the cell, resulting in changes in the reaction thermodynamics and kinetics. Recent studies, both experiment and simulations, revealed that the nature of the interaction between crowder and protein species, in particular the soft interactions, plays an important role in crowder induced effects on protein association. To this end, from a simulation perspective, it is important to decipher the level of structural resolution in a protein-crowder model that can faithfully capture the influence of crowding on protein association. Here, we investigate the dimerization of model system GB1 in the presence of lysozyme crowders at two structural resolutions. The lower resolution model assumes both protein and crowder species as spherical beads, similar to the analytical scaled particle theory model, whereas the higher resolution model retains residue specific structural details for protein and crowder species. From the higher resolution model, it is found that GB1 dimer formation is destabilized in the presence of lysozyme crowders, and the destabilization is more for the side-by-side dimer compared to the domain-swapped dimer, in qualitative agreement with experimental findings. However, the low resolution CG model predicts stabilization of the dimers in the presence of the lysozyme crowder, similar to the SPT model. Our results indicate a nontrivial role of the choice of model resolution in computer simulation studies investigating crowder induced effects.

Intrinsic Dynamics of the ClpXP Proteolytic Machine Using Elastic Network Models

ClpXP complex is an ATP-dependent mitochondrial matrix protease that binds, unfolds, translocates, and subsequently degrades specific protein substrates. Its mechanisms of operation are still being debated, and several have been proposed, including the sequential translocation of two residues (SC/2R), six residues (SC/6R), and even long-pass probabilistic models. Therefore, it has been suggested to employ biophysical-computational approaches that can determine the kinetics and thermodynamics of the translocation. In this sense, and based on the apparent inconsistency between structural and functional studies, we propose to apply biophysical approaches based on elastic network models (ENM) to study the intrinsic dynamics of the theoretically most probable hydrolysis mechanism. The proposed models ENM suggest that the ClpP region is decisive for the stabilization of the ClpXP complex, contributing to the flexibility of the residues adjacent to the pore, favoring the increase in pore size and, therefore, with the energy of interaction of its residues with a larger portion of the substrate. It is predicted that the complex may undergo a stable configurational change once assembled and that the deformability of the system once assembled is oriented, to increase the rigidity of the domains of each region (ClpP and ClpX) and to gain flexibility of the pore. Our predictions could suggest under the conditions of this study the mechanism of the interaction of the system, of which the substrate passes through the unfolding of the pore in parallel with a folding of the bottleneck. The variations in the distance calculated by molecular dynamics could allow the passage of a substrate with a size equivalent to ∼3 residues. The theoretical behavior of the pore and the stability and energy of binding to the substrate based on ENM models suggest that in this system, there are thermodynamic, structural, and configurational conditions that allow a possible translocation mechanism that is not strictly sequential.

Structurally Different Yet Functionally Similar: Aptamers Specific for the Ebola Virus Soluble Glycoprotein and GP1,2 and Their Application in Electrochemical Sensing

The Ebola virus glycoprotein (GP) gene templates several mRNAs that produce either the virion-associated transmembrane protein or one of two secreted glycoproteins. Soluble glycoprotein (sGP) is the predominant product. GP1 and sGP share an amino terminal sequence of 295 amino acids but differ in quaternary structure, with GP1 being a heterohexamer with GP2 and sGP a homodimer. Two structurally different DNA aptamers were selected against sGP that also bound GP1,2. These DNA aptamers were compared with a 2'FY-RNA aptamer for their interactions with the Ebola GP gene products. The three aptamers have almost identical binding isotherms for sGP and GP1,2 in solution and on the virion. They demonstrated high affinity and selectivity for sGP and GP1,2. Furthermore, one aptamer, used as a sensing element in an electrochemical format, detected GP1,2 on pseudotyped virions and sGP with high sensitivity in the presence of serum, including from an Ebola-virus-infected monkey. Our results suggest that the aptamers interact with sGP across the interface between the monomers, which is different from the sites on the protein bound by most antibodies. The remarkable similarity in functional features of three structurally distinct aptamers suggests that aptamers, like antibodies, have preferred binding sites on proteins.

Crowding-induced phase separation and gelling by co-condensation of PEG in NPM1-rRNA condensates

The crowdedness of the cell calls for adequate intracellular organization. Biomolecular condensates, formed by liquid-liquid phase separation of intrinsically disordered proteins and nucleic acids, are important organizers of cellular fluids. To underpin the molecular mechanisms of protein condensation, cell-free studies are often used where the role of crowding is not investigated in detail. Here, we investigate the effects of macromolecular crowding on the formation and material properties of a model heterotypic biomolecular condensate, consisting of nucleophosmin (NPM1) and ribosomal RNA (rRNA). We studied the effect of the macromolecular crowding agent poly(ethylene glycol) (PEG), which is often considered an inert crowding agent. We observed that PEG could induce both homotypic and heterotypic phase separation of NPM1 and NPM1-rRNA, respectively. Crowding increases the condensed concentration of NPM1 and decreases its equilibrium dilute phase concentration, although no significant change in the concentration of rRNA in the dilute phase was observed. Interestingly, the crowder itself is concentrated in the condensates, suggesting that co-condensation rather than excluded volume interactions underlie the enhanced phase separation by PEG. Fluorescence recovery after photobleaching measurements indicated that both NPM1 and rRNA become immobile at high PEG concentrations, indicative of a liquid-to-gel transition. Together, these results provide more insight into the role of synthetic crowding agents in phase separation and demonstrate that condensate properties determined in vitro depend strongly on the addition of crowding agents.

PEG mediated destabilization of holo α-lactalbumin probed by in silico and in vitro studies: deviation from excluded volume effect

Crowded and confined macromolecular milieus surround proteins, and both are stabilizing if the nature of the interaction between crowder and proteins are considered hard-core repulsive interactions. However, non-specific chemical interactions between a protein and its surroundings also play a significant role and the sum effect of both hard-core repulsion and soft interaction balances the overall effect of crowding/confinement. Previous studies showing the effect of polyethylene glycol (PEG) on protein and nucleic acid may be interpreted as either primarily excluded volume effect or, in some cases, chemical effect by changing solvent properties. In case of destabilizing interactions, charge-charge and hydrophobic contact have to gain more attention. For instance, in vitro and in vivo studies using protein as crowding agent revealed the destabilization of proteins induced by charge-charge interactions. To investigate the effect of PEG 10 kDa on holo α-lactalbumin (holo α-LA), structure and thermal stability of the protein were measured at different pH values using several techniques. Structural characterization by Trp-fluorescence, near-UV CD and far-UV measurements at different pH values clearly shows perturbation of tertiary and secondary structure of holo α-LA by PEG 10 kDa. Furthermore, the dynamic light scattering measurement shows that the protein is homogeneous under all experimental conditions. Analysis of the heat-induced denaturation profile in the presence of the crowder shows destabilization of the protein in terms of T_m (midpoint of denaturation) and ΔG_D⁰ (Gibbs free energy change at 25 °C). To evaluate the interaction of PEG 10 kDa with holo α-LA and stability of PEG-α-LA complex, docking and molecular dynamic simulation were carried out for 100 ns.Communicated by Ramaswamy H. Sarma.

The biophysics of disordered proteins from the point of view of single-molecule fluorescence spectroscopy

Intrinsically disordered proteins (IDPs) and regions (IDRs) have emerged as key players across many biological functions and diseases. Differently from structured proteins, disordered proteins lack stable structure and are particularly sensitive to changes in the surrounding environment. Investigation of disordered ensembles requires new approaches and concepts for quantifying conformations, dynamics, and interactions. Here, we provide a short description of the fundamental biophysical properties of disordered proteins as understood through the lens of single-molecule fluorescence observations. Single-molecule Förster resonance energy transfer (FRET) and fluorescence correlation spectroscopy (FCS) provides an extensive and versatile toolbox for quantifying the characteristics of conformational distributions and the dynamics of disordered proteins across many different solution conditions, both in vitro and in living cells.

Crowding-induced morphological changes in synthetic lipid vesicles determined using smFRET

Lipid vesicles are valuable mesoscale molecular confinement vessels for studying membrane mechanics and lipid-protein interactions, and they have found utility among bio-inspired technologies, including drug delivery vehicles. While vesicle morphology can be modified by changing the lipid composition and introducing fusion or pore-forming proteins and detergents, the influence of extramembrane crowding on vesicle morphology has remained under-explored owing to a lack of experimental tools capable of capturing morphological changes on the nanoscale. Here, we use biocompatible polymers to simulate molecular crowding in vitro, and through combinations of FRET spectroscopy, lifetime analysis, dynamic light scattering, and single-vesicle imaging, we characterize how crowding regulates vesicle morphology. We show that both freely diffusing and surface-tethered vesicles fluorescently tagged with the DiI and DiD FRET pair undergo compaction in response to modest concentrations of sorbitol, polyethylene glycol, and Ficoll. A striking observation is that sorbitol results in irreversible compaction, whereas the influence of high molecular weight PEG-based crowders was found to be reversible. Regulation of molecular crowding allows for precise control of the vesicle architecture in vitro, with vast implications for drug delivery and vesicle trafficking systems. Furthermore, our observations of vesicle compaction may also serve to act as a mechanosensitive readout of extramembrane crowding.

An in vitro cytomimetic of in-cell RNA folding

To discover the cytomimetic that accounts for cytoplasmic crowding and sticking on RNA stability, we conducted a two-dimensional scan of mixtures of artificial crowding and sticking agents, PEG10k and M-PERTM. As our model RNA, we investigate the fourU RNA thermometer motif of Salmonella, a hairpin-structured RNA that regulates translation by unfolding and exposing its RBS in response to temperature perturbations. We found that the addition of artificial crowding and sticking agents leads to a stabilization and destabilization of RNA folding, respectively, through the excluded volume effect and surface interactions. FRET-labels were added to the fourU RNA and Fast Relaxation Imaging (FReI), fluorescence microscopy coupled to temperature-jump spectroscopy, probed differences between folding stability of RNA inside single living cells and in vitro. Our results suggest that the cytoplasmic environment affecting RNA folding is comparable to a combination of 20% v/v M-PERTM and 150 g/L PEG10k.

Synergism in the Molecular Crowding of Ligand-Induced Riboswitch Folding: Kinetic/Thermodynamic Insights from Single-Molecule Spectroscopy

Conformational dynamics in riboswitches involves ligand binding and folding of RNA, each of which can be influenced by excluded volume effects under "crowded" in vivo cellular conditions and thus incompletely characterized by in vitro studies under dilute buffer conditions. In this work, temperature-dependent single-molecule fluorescence resonance energy transfer (FRET) spectroscopy is used to characterize the thermodynamics of (i) cognate ligand and (ii) molecular crowders (PEG, polyethylene glycol) on folding of the B. subtilis LysC lysine riboswitch. With the help of detailed kinetic analysis, we isolate and study the effects of PEG on lysine binding and riboswitch folding steps individually, from which we find that PEG crowding facilitates riboswitch folding primarily via a surprising increase in affinity for the cognate ligand. This is furthermore confirmed by temperature-dependent studies, which reveal that PEG crowding is not purely entropic and instead significantly impacts both enthalpic and entropic contributions to the free energy landscape for folding. The results indicate that PEG molecular crowding/stabilization of the lysine riboswitch is more mechanistically complex and requires extension beyond the conventional picture of purely repulsive solvent-solute steric interactions arising from excluded volume and entropy. Instead, the current experimental FRET data support an alternative multistep mechanism, whereby PEG first entropically crowds the unfolded riboswitch into a "pre-folded" conformation, which in turn greatly increases the ligand binding affinity and thereby enhances the overall equilibrium for riboswitch folding.

Sensitive and Low-Bias Transcriptome Sequencing Using Agarose PCR

Transcriptome sequencing has emerged as an important research tool for exploring the mysteries of life at the single-cell level. However, its wide application is limited by the bias associated with the amplification reactions which is essential for library building of trace RNA. In this study, low-melting-point agarose was added to the amplification reactions to take advantage of its molecular crowding effect and polymer cross-linked structure to improve the sensitivity of the reactions and reduce bias. To further evaluate the performance of the method, it was applied to transcriptome sequencing of microregion samples from brain tissue sections of mice with Parkinson's disease at the single cell level. The results showed that agarose PCR had better performance than in-tube PCR. Further application of agarose PCR to transcriptome library sequencing could obtain data closer to that of unamplified. With the addition of low melting point agarose, the sensitivity of the amplification reaction was significantly increased, while homogeneity was increased by approximately 2-fold. Not only that, but this work also provides 11% sensitivity improvement for spatial transcriptomic study on Parkinson's disease-associated gene detection. The agarose PCR provides a new tool for efficient and homogeneous amplification of trace samples and can be widely used for spatial transcriptome library sequencing and studies.

Comparative study of SARS-CoV-2 infection in different cell types: Biophysical-computational approach to the role of potential receptors

Cellular susceptibility to SARS-CoV-2 infection in the respiratory tract has been associated with the ability of the virus to interact with potential receptors on the host membrane. We have modeled viral dynamics by simulating various cellular systems and artificial conditions, including macromolecular crowding, based on experimental and transcriptomic data to infer parameters associated with viral growth and predict cell susceptibility. We have accomplished this based on the type, number and level of expression of the angiotensin-converting enzyme 2 (ACE2), transmembrane serine 2 (TMPRSS2), basigin2 (CD147), FURIN protease, neuropilin 1 (NRP1) or other less studied candidate receptors such as heat shock protein A5 (HSPA5) and angiotensin II receptor type 2 (AGTR2). In parallel, we studied the effect of simulated artificial environments on the accessibility to said proposed receptors. In addition, viral kinetic behavior dependent on the degree of cellular susceptibility was predicted. The latter was observed to be more influenced by the type of proteins and expression level, than by the number of potential proteins associated with the SARS CoV-2 infection. We predict a greater theoretical propensity to susceptibility in cell lines such as NTERA-2, SCLC-21H, HepG2 and Vero6, and a lower theoretical propensity in lines such as CaLu3, RT4, HEK293, A549 and U-251MG. An important relationship was observed between expression levels, protein diffusivity, and thermodynamically favorable interactions between host proteins and the viral spike, suggesting potential sites of early infection other than the lungs. This research is expected to stimulate future quantitative experiments and promote systematic investigation of the effect of crowding presented here.

Kineto-Mechanistic Investigation of Effect of Macromolecular Crowding on the Breathing of DNA Bubble

Macromolecular crowding along with hydrogen bonding or stacking interactions and hydration reportedly has enormous repercussions on elementary biochemical processes, such as the folding of proteins or nucleic acids involving the stability of DNA base pairing. By using the mismatch-induced DNA bubble as a mesoscopic model, the complex interplay of macromolecular crowding on the dynamical fluctuations at the bubble region within the thermodynamic limit has been monitored using single-molecule fluorescence resonance energy transfer (sm-FRET). These single-molecule experimental results have been further corroborated using physical models such as "scaled particle theory" (SPT) and "Gaussian cloud model" (GCM), to predict the biological activity of DNA. The two-state fluctuation of the DNA bubble has been visualized as a function of the nature, size, and concentration of the crowder. The influence of crowders on the DNA conformation has been investigated with the help of the m-factor, the eccentricity, and the kinetic and thermodynamic parameters without any prior assumption. The clear effect of crowding on the dynamics of such a simple biomolecular system emphasizes the power of single-molecule methods and the dependency of the radius of gyration of the co-solute as well as the preferential interaction with the crowder on the distinct conformational states adopted by the bubble. This study provides an idea and hypothesizes the preferential propensity of the DNA bubble to adopt a conformation with the single-stranded domains being far apart, independent of the crowder size, that may be beneficial for efficient recognition by proteins for an uninterrupted procession of the biological process of the central dogma.

A pressure-jump study on the interaction of osmolytes and crowders with cubic monoolein structures

Many vital processes that take place in biological cells involve remodeling of lipid membranes. These processes take place in a milieu that is packed with various solutes, ranging from ions and small organic osmolytes to proteins and other macromolecules, occupying about 30% of the available volume. In this work, we investigated how molecular crowding, simulated with the polymer polyethylene glycol (PEG), and the osmolytes urea and trimethylamine-N-oxide (TMAO) affect the equilibration of cubic monoolein structures after a phase transition from a lamellar state induced by an abrupt pressure reduction. In absence of additives, swollen cubic crystallites form after the transition, releasing excess water over several hours. This process is reflected in a decreasing lattice constant and was monitored with small angle X-ray scattering. We found that the osmotic pressure exerted by PEG and TMAO, which are displaced from narrow inter-bilayer spaces, accelerates the equilibration. When the radius of gyration of the added PEG was smaller than the radius of the water channels of the cubic phase, the effect became more pronounced with increasing molecular weight of the polymers. As the release of hydration water from the cubic structures is accompanied by an increasing membrane curvature and a reduction of the interface between lipids and aqueous phase, urea, which has a slight affinity to reside near membrane surfaces, stabilized the swollen crystallites and slowed down the equilibration dynamics. Our results support the view that cellular solutes are important contributors to dynamic membrane processes, as they can accelerate dehydration of inter-bilayer spaces and promote or counteract membrane curvature.

An anionic ligand snap-locks a long-range interaction in a magnesium-folded riboswitch

The archetypical transcriptional crcB fluoride riboswitch from Bacillus cereus is an intricately structured non-coding RNA element enhancing gene expression in response to toxic levels of fluoride. Here, we used single molecule FRET to uncover three dynamically interconverting conformations appearing along the transcription process: two distinct undocked states and one pseudoknotted docked state. We find that the fluoride anion specifically snap-locks the magnesium-induced, dynamically docked state. The long-range, nesting, single base pair A40-U48 acts as the main linchpin, rather than the multiple base pairs comprising the pseudoknot. We observe that the proximally paused RNA polymerase further fine-tunes the free energy to promote riboswitch docking. Finally, we show that fluoride binding at short transcript lengths is an early step toward partitioning folding into the docked conformation. These results reveal how the anionic fluoride ion cooperates with the magnesium-associated RNA to govern regulation of downstream genes needed for fluoride detoxification of the cell.

Vesicle encapsulation stabilizes intermolecular association and structure formation of functional RNA and DNA

During the origin of life, encapsulation of RNA inside vesicles is believed to have been a defining feature of the earliest cells (protocells). The confined biophysical environment provided by membrane encapsulation differs from that of bulk solution and has been shown to increase activity as well as evolutionary rate for functional RNA. However, the structural basis of the effect on RNA has not been clear. Here, we studied how encapsulation of the hairpin ribozyme inside model protocells affects ribozyme kinetics, ribozyme folding into the active conformation, and cleavage and ligation activities. We further examined the effect of encapsulation on the folding of a stem-loop RNA structure and on the formation of a triplex structure in a pH-sensitive DNA switch. The results indicate that encapsulation promotes RNA-RNA association, both intermolecular and intramolecular, and also stabilizes tertiary folding, including the docked conformation characteristic of the active hairpin ribozyme and the triplex structure. The effects of encapsulation were sufficient to rescue the activity of folding-deficient mutants of the hairpin ribozyme. Stabilization of multiple modes of nucleic acid folding and interaction thus enhanced the activity of encapsulated nucleic acids. Increased association between RNA molecules may facilitate the formation of more complex structures and cooperative interactions. These effects could promote the emergence of biological functions in an "RNA world" and may have utility in the construction of minimal synthetic cells.

Modulating internal transition kinetics of responsive macromolecules by collective crowding

Packing and crowding are used in biology as mechanisms to (self-)regulate internal molecular or cellular processes based on collective signaling. Here, we study how the transition kinetics of an internal "switch" of responsive macromolecules is modified collectively by their spatial packing. We employ Brownian dynamics simulations of a model of Responsive Colloids, in which an explicit internal degree of freedom-here, the particle size-moving in a bimodal energy landscape self-consistently responds to the density fluctuations of the crowded environment. We demonstrate that populations and transition times for the two-state switching kinetics can be tuned over one order of magnitude by "self-crowding." An exponential scaling law derived from a combination of Kramers' and liquid state perturbation theory is in very good agreement with the simulations.

Effects of Molecular Crowders on Single-Molecule Nucleic Acid Folding: Temperature-Dependent Studies Reveal True Crowding vs Enthalpic Interactions

Biomolecular folding in cells can be strongly influenced by spatial overlap/excluded volume interactions (i.e., "crowding") with intracellular solutes. As a result, traditional in vitro experiments with dilute buffers may not accurately recapitulate biomolecule folding behavior in vivo. In order to account for such ubiquitous excluded volume effects, biologically inert polyethylene glycol (PEG) and polysaccharides (dextran and Ficoll) are often used as in vitro crowding agents to mimic in vivo crowding conditions, with a common observation that high concentrations of these polymers stabilize the more compact biomolecule conformation. However, such an analysis can be distorted by differences in polymer interactions with the folded vs unfolded conformers, requiring temperature-dependent analysis of the thermodynamics to reliably assess competing enthalpic vs entropic contributions and thus the explicit role of excluded volume. In this work, temperature-controlled single-molecule fluorescence resonance energy transfer (smFRET) is used to characterize the thermodynamic interaction between nucleic acids and common polymer crowders PEG, dextran, and Ficoll. The results reveal that PEG promotes secondary and tertiary nucleic acid folding by simultaneously increasing the folding rate while decreasing the unfolding rate, with temperature-dependent studies confirming that the source of PEG stabilization is predominantly entropic and consistent with a true excluded volume crowding mechanism. By way of contrast, neither dextran nor Ficoll induces any significant concentration-dependent change in nucleic acid folding stability at room temperature, but instead, stabilization effects gradually appear with a temperature increase. Such a thermal response indicates that both folding enthalpies and entropies are impacted by dextran and Ficoll. A detailed thermodynamic analysis of the kinetics suggests that, instead of true entropic molecular crowding, dextran and Ficoll associate preferentially with the unfolded vs folded nucleic acid conformer as a result of larger solvent accessible surface area, thereby skewing the free energy landscapes through both significant entropic/enthalpic contributions that compete and fortuitously cancel near room temperature.

Discovery of highly reactive self-splicing group II introns within the mitochondrial genomes of human pathogenic fungi

Fungal pathogens represent an expanding global health threat for which treatment options are limited. Self-splicing group II introns have emerged as promising drug targets, but their development has been limited by a lack of information on their distribution and architecture in pathogenic fungi. To meet this challenge, we developed a bioinformatic workflow for scanning sequence data to identify unique RNA structural signatures within group II introns. Using this approach, we discovered a set of ubiquitous introns within thermally dimorphic fungi (genera of Blastomyces, Coccidioides and Histoplasma). These introns are the most biochemically reactive group II introns ever reported, and they self-splice rapidly under near-physiological conditions without protein cofactors. Moreover, we demonstrated the small molecule targetability of these introns by showing that they can be inhibited by the FDA-approved drug mitoxantrone in vitro. Taken together, our results highlight the utility of structure-based informatic searches for identifying riboregulatory elements in pathogens, revealing a striking diversity of reactive self-splicing introns with great promise as antifungal drug targets.

The RNA chaperone StpA enables fast RNA refolding by destabilization of mutually exclusive base pairs within competing secondary structure elements

In bacteria RNA gene regulatory elements refold dependent on environmental clues between two or more long-lived conformational states each associated with a distinct regulatory state. The refolding kinetics are strongly temperature-dependent and especially at lower temperatures they reach timescales that are biologically not accessible. To overcome this problem, RNA chaperones have evolved. However, the precise molecular mechanism of how these proteins accelerate RNA refolding reactions remains enigmatic. Here we show how the RNA chaperone StpA of Escherichia coli leads to an acceleration of a bistable RNA's refolding kinetics through the selective destabilization of key base pairing interactions. We find in laser assisted real-time NMR experiments on photocaged bistable RNAs that the RNA chaperone leads to a two-fold increase in refolding rates at low temperatures due to reduced stability of ground state conformations. Further, we can show that upon interaction with StpA, base pairing interactions in the bistable RNA are modulated to favor refolding through the dominant pseudoknotted transition pathway. Our results shed light on the molecular mechanism of the interaction between RNA chaperones and bistable RNAs and are the first step into a functional classification of chaperones dependent on their biophysical mode of operation.

Comparative study of the interaction of ivermectin with proteins of interest associated with SARS-CoV-2: A computational and biophysical approach

The SARS-CoV-2 pandemic has accelerated the study of existing drugs. The mixture of homologs called ivermectin (avermectin-B1a [HB1a] + avermectin-B1b [HB1b]) has shown antiviral activity against SARS-CoV-2 in vitro. However, there are few reports on the behavior of each homolog. We investigated the interaction of each homolog with promising targets of interest associated with SARS-CoV-2 infection from a biophysical and computational-chemistry perspective using docking and molecular dynamics. We observed a differential behavior for each homolog, with an affinity of HB1b for viral structures, and of HB1a for host structures considered. The induced disturbances were differential and influenced by the hydrophobicity of each homolog and of the binding pockets. We present the first comparative analysis of the potential theoretical inhibitory effect of both avermectins on biomolecules associated with COVID-19, and suggest that ivermectin through its homologs, has a multiobjective behavior.

Leveraging lifetime information to perform real-time 3D single-particle tracking in noisy environments

A microscopy platform that leverages the arrival time of individual photons to enable 3D single-particle tracking of fast-moving (translational diffusion coefficient of ≃3.3 µm²/s) particles in high-background environments is reported here. It combines a hardware-based time-gating module, which enables the rate of photon processing to be as high as 100 MHz, with a two-photon-excited 3D single-particle tracking confocal microscope to enable high sample penetration depth. Proof-of-principle experiments where single quantum dots are tracked in solutions containing dye-stained cellulose, are shown with tracking performance markedly improved using the hardware-based time-gating module. Such a microscope design is anticipated to be of use to a variety of communities who wish to track single particles in cellular environments, which commonly have high fluorescence and scattering background.

Fluorine NMR Spectroscopy Enables to Quantify the Affinity Between DNA and Proteins in Cell Lysate

The determination of the binding affinity quantifying the interaction between proteins and nucleic acids is of crucial interest in biological and chemical research. Here, we have made use of site-specific fluorine labeling of the cold shock protein from Bacillus subtilis, BsCspB, enabling to directly monitor the interaction with single stranded DNA molecules in cell lysate. High-resolution 19 F NMR spectroscopy has been applied to exclusively report on resonance signals arising from the protein under study. We have found that this experimental approach advances the reliable determination of the binding affinity between single stranded DNA molecules and its target protein in this complex biological environment by intertwining analyses based on NMR chemical shifts, signal heights, line shapes and simulations. We propose that the developed experimental platform offers a potent approach for the identification of binding affinities characterizing intermolecular interactions in native surroundings covering the nano-to-micromolar range that can be even expanded to in cell applications in future studies.

Effect of Molecular Crowding on Complexation of Metal Ions and 8-Quinolinol-5-Sulfonic Acid

The effect of molecular crowding on macromolecular reactions has been revealed by many researchers. In this study, we investigate the complexation of metal ions (Zn, Co, and Cd) with 8-quinolinol-5-sulfonic acid as a model of small-molecular reactions in molecular crowding. The complexation constants for 1:1, 1:2, and total complexation in the presence of polyethylene glycol (PEG, a molecular crowding reagent) are evaluated based on the increase in the reactant activity by volume exclusion and the decrease in the water activity due to the change in osmotic pressure. All complexation constants are enhanced by increasing the concentration of PEG. Its mechanisms differ for 1:1, 1:2, and total complexation. The 1:1 complexation is promoted only by the influence of the water activity, while the reactant and water activities influence the increase in the 1:2 complexation constant. Increasing the molecular weight of PEG further increases the complexation constants, as dehydration of the complex is promoted by a higher hydration number of PEG. Because this study gives the fundamental knowledge for the protein-metal interaction, in which solvation is an important factor, in molecular crowding, it provides new insights into molecular crowding studies and should attract the attention of a broad spectrum of biochemistry researchers.

Effects of chain length of polyethylene glycol molecular crowders on a mutant Tetrahymena group I ribozyme lacking large peripheral module

While current group I ribozymes use several distinct strategies to function under conditions of low Mg²⁺ concentration (≤ 3 mM), a deletion mutant of the Tetrahymena ribozyme (ΔP5 ribozyme) is virtually inactive with 3 mM Mg²⁺ due to removal of the large peripheral module, P5abc, supporting the active conformation of the core module. We investigated the molecular crowding effects of synthetic polyethylene glycols (PEGs) on the activity of the ΔP5 ribozyme. Among PEG molecules with different chain lengths, PEG600 improved the activity of the ΔP5 ribozyme most effectively in the presence of 3 mM Mg²⁺.

Liquid-liquid phase separation as a common organizing principle of intracellular space and biomembranes providing dynamic adaptive responses

This work is devoted to the phenomenon of liquid-liquid phase separation (LLPS), which has come to be recognized as fundamental organizing principle of living cells. We distinguish separation processes with different dimensions. Well-known 3D-condensation occurs in aqueous solution and leads to membraneless organelle (MLOs) formation. 2D-films may be formed near membrane surfaces and lateral phase separation (membrane rafts) occurs within the membranes themselves. LLPS may also occur on 1D structures like DNA and the cyto- and nucleoskeleton. Phase separation provides efficient transport and sorting of proteins and metabolites, accelerates the assembly of metabolic and signaling complexes, and mediates stress responses. In this work, we propose a model in which the processes of polymerization (1D structures), phase separation in membranes (2D structures), and LLPS in the volume (3D structures) influence each other. Disordered proteins and whole condensates may provide membrane raft separation or polymerization of specific proteins. On the other hand, 1D and 2D structures with special composition or embedded IDRs can nucleate condensates. We hypothesized that environmental change may trigger a LLPS which can propagate within the cell interior moving along the cytoskeleton or as an autowave. New phase propagation quickly and using a low amount of energy adjusts cell signaling and metabolic systems to new demands. Cumulatively, the interconnected phase separation phenomena in different dimensions represent a previously unexplored system of intracellular communication and regulation which cannot be ignored when considering both physiological and pathological cell processes.

Modeling protein association from homogeneous to mixed environments: A reaction-diffusion dynamics approach

Protein-protein association in vivo occur in a crowded and complex environment. Theoretical models based on hard-core repulsion predict stabilization of the product under crowded conditions. Soft interactions, on the contrary, can either stabilize or destabilize the product formation. Here we modeled protein association in presence of crowders of varying size, shape, interaction potential and used different mixing parameters for constituent crowders to study the influence on the association reaction. It was found that size is a more dominant factor in crowder-induced stabilization than the shape. Furthermore, in a mixture of crowders having different sizes but identical interaction potential, the change of free energy is additive of the free energy changes produced by individual crowders. However, the free energy change is not additive if two crowders of same size interact via different interaction potentials. These findings provide a systematic understanding of crowding influences in heterogeneous medium.

Comparative study of polymer looping kinetics in passive and active environments

Intra-chain looping in complex environments is significant in advancing our understanding of biological processes in life. We adopt Langevin dynamics simulations to perform a comparative study of polymer looping kinetics in passive and active environments. From the analysis of looping quantities, including looping-unlooping times and looping probabilities, we unraveled the intriguing effects of active crowder size, activity and crowding. Firstly, we figured out the phase diagram involving a novel facilitation-inhibition transition in the parameter space of active crowder size and active force, and the two-fold roles of activity are clarified. In particular, we find that active particles of a size comparable to the polymer monomer are most favorable for facilitated looping, while those with a similar size to the polymer gyration radius impede the looping most seriously. Secondly, the underlying looping mechanisms in different active crowder size regimes are rationalized by the interplay among diffusion, polymer conformational change and the free-energy barrier. For small active crowders, activity significantly promotes end-to-end distance diffusion, which dominantly facilitates both looping and unlooping processes. In the case of moderate active crowders, the polymer chain suffers from prominent swelling, and thus inevitable inhibited looping will occur. For large active crowders, activity induces a counterintuitive non-cage effect on the looping kinetics, through yielding a higher effective temperature and larger unlooping free-energy barrier. This is in sharp contrast to the caging phenomena observed in passive media. Lastly, the volume-fraction dependence of the looping quantities in an active bath demonstrates dramatic discrepancies from that in a passive bath, which highlights the contrasting effects of activity and crowding.

Pathway regulation mechanism revealed by cotranslational folding of villin headpiece subdomain HP35

Cotranslational folding is one of the most important features of protein folding in vivo. Although many studies have shown that the folding pathways of cotranslational folding are different from free folding in vitro, the detailed mechanism of folding dynamics is lacking. Here we combine all-atom molecular simulations with an ideal ribosome tunnel model to investigate the cotranslational folding of villin headpiece subdomain HP35. By comparing the folding dynamics between cotranslational folding and free folding, we found that cotranslational folding tends to fold along the pathway that is easier to fold into native state in the latter. In addition, the roles of the ribosome tunnel and sequential folding are analyzed separately. Our results show that the ribosome can prevent the untimely folding of the C segment of HP35 to reduce the non-native interactions, while the translation speed can regulate the amounts of native and non-native interactions and the balance between them. Overall, these results give insights into the general mechanisms of cotranslational protein folding.

Smaller molecules crowd better: Crowder size dependence revealed by single-molecule FRET studies and depletion force modeling analysis

The cell is an extremely crowded environment, which is known to have a profound impact on the thermodynamics, functionality, and conformational stability of biomolecules. Speculations from recent theoretical molecular dynamics studies suggest an intriguing size dependence to such purely entropic crowding effects, whereby small molecular weight crowders under constant enthalpy conditions are more effective than larger crowders on a per volume basis. If experimentally confirmed, this would be profoundly significant, as the cellular cytoplasm is also quite concentrated in smaller molecular weight solutes such as inorganic ions, amino acids, and various metabolites. The challenge is to perform such studies isolating entropic effects under isoenthalpic conditions. In this work, we first present results from single-molecule FRET spectroscopy (smFRET) on the molecular size-dependent crowding stabilization of a simple RNA tertiary motif (the GAAA tetraloop-tetraloop receptor), indeed providing evidence in support of the surprising notion in the crowding literature that "smaller is better." Specifically, systematic smFRET studies as a function of crowder solute size reveal that smaller molecules both significantly increase the RNA tertiary folding rate and, yet, simultaneously decrease the unfolding rate, predicting strongly size-dependent stabilization of RNA tertiary structures under crowded cellular conditions. The size dependence of these effects has been explored via systematic variation of crowder size over a broad range of molecular weights (90-3000 amu). Furthermore, corresponding temperature dependent studies indicate the systematic changes in the folding equilibrium to be predominantly entropic in origin, i.e., consistent with a fundamental picture of entropic molecular crowding without additional enthalpic interactions. Most importantly, all trends in the single-molecule crowding data can be quantitatively recapitulated by a simple analytic depletion force model, whereby excluded volume interactions represent the major thermodynamic driving force toward folding. Our study, thus, not only provides experimental evidence and theoretical support for small molecule crowding but also predicts further enhancement of crowding effects for even smaller molecules on a per volume basis.

Quantify the combined effects of temperature and force on the stability of DNA hairpin

OxDNA, as a successful coarse-grain model, has been applied to reproduce the thermodynamic and mechanical properties of both single- and double-stranded DNA. In current simulation, oxDNA is extended to explore the combined effects of temperature and force on the stability of DNA hairpin and its free energy landscape. Simulations were carried out at different forces and temperatures, at each temperature, a 18-base-pair DNA hairpin dynamically transited between folded state and unfolded state, and the separation between two states is consistent with the full contour length of single-stranded DNA in the unfolded state. Two methods were used to identify the critical force of DNA hairpin at each temperature and the critical forces obtained from two methods were consistent with each other and gradually decreased with the increasing temperature from 300 K to 326 K. The critical force at 300 K is reasonably consistent with the single molecule result of DNA hairpin with the same stem length. The two-state free energy landscape can be elucidated from the probability distribution of DNA hairpin extension and its dependence on the force and temperature is totally different. The increasing temperature not only reduces the free energy barrier, but also alters the position of transition point along the extension coordinate, resulting in the reduction of folding distance and the extension of unfolding distance, but their sum is not obviously dependent on the temperature. Generally, an assumption that the location of transition state in two-state energy landscape is independent of the stretching force is used to analyze the data of the single molecule experiment, but current simulation results indicate that effects of stretching forces on the location of transition state in two-state energy landscape are dependent on temperature. At relatively high temperature, stretching force can also change the location of transition state in the free energy landscape.

Salt-Dependent RNA Pseudoknot Stability: Effect of Spatial Confinement

Macromolecules, such as RNAs, reside in crowded cell environments, which could strongly affect the folded structures and stability of RNAs. The emergence of RNA-driven phase separation in biology further stresses the potential functional roles of molecular crowding. In this work, we employed the coarse-grained model that was previously developed by us to predict 3D structures and stability of the mouse mammary tumor virus (MMTV) pseudoknot under different spatial confinements over a wide range of salt concentrations. The results show that spatial confinements can not only enhance the compactness and stability of MMTV pseudoknot structures but also weaken the dependence of the RNA structure compactness and stability on salt concentration. Based on our microscopic analyses, we found that the effect of spatial confinement on the salt-dependent RNA pseudoknot stability mainly comes through the spatial suppression of extended conformations, which are prevalent in the partially/fully unfolded states, especially at low ion concentrations. Furthermore, our comprehensive analyses revealed that the thermally unfolding pathway of the pseudoknot can be significantly modulated by spatial confinements, since the intermediate states with more extended conformations would loss favor when spatial confinements are introduced.