Crowding Promotes the Switch from Hairpin to Pseudoknot Conformation in Human Telomerase RNA

Deeper Insights into Mixed Crowding through Enzyme Activity, Dynamics, and Crowder Diffusion

Given the fact that the cellular interior is crowded by many different kinds of macromolecules, it is important that in vitro studies be carried out in the presence of mixed crowder systems. In this regard, we have used binary crowders formed by the combination of some of the commonly used crowding agents, namely, Ficoll 70, Dextran 70, Dextran 40, and PEG 8000 (PEG 8), to study how these affect enzyme activity, dynamics, and crowder diffusion. The enzyme chosen is AK3L1, an isoform of adenylate kinase. To investigate its dynamics, we have carried out three single point mutations (A74C, A132C, and A209C) with the cysteine residues being labeled with a coumarin-based solvatochromic probe [CPM: (7-diethylamino-3-(4-maleimido-phenyl)-4-methylcoumarin)]. Both enzyme activity and dynamics decreased in the binary mixtures as compared with the sum of the individual crowders, suggesting a reduction in excluded volume (in the mixture). To gain deeper insights into the binary mixtures, fluorescence correlation spectroscopy studies were carried out using fluorescein isothiocyanate-labeled Dextran 70 and tetramethylrhodamine-labeled AK3L1 as the diffusion probes. Diffusion in binary mixtures was observed to be much more constrained (relative to the sum of the individual crowders) for the labeled enzyme as compared to the labeled crowder showing different environments being faced by the two species. This was further confirmed during imaging of the phase-separated droplets formed in the binary mixtures having PEG as one of the crowding agents. The interior of these droplets was found to be rich in crowders and densely packed, as shown by confocal and digital holographic microscopy images, with the enzymes predominantly residing outside these droplets, that is, in the relatively less crowded regions. Taken together, our data provide important insights into various aspects of the simplest form of mixed crowding, that is, composed of just two components, and also hint at the enhanced complexity that the cellular interior presents toward having a detailed and comprehensive understanding of the same.

Crowder titrations enable the quantification of driving forces for macromolecular phase separation

Macromolecular solubility is an important contributor to the driving forces for phase separation. Formally, the driving forces in a binary mixture comprising a macromolecule dissolved in a solvent can be quantified in terms of the saturation concentration, which is the threshold macromolecular concentration above which the mixture separates into coexisting dense and dilute phases. In addition, the second virial coefficient, which measures the effective strength of solvent-mediated intermolecular interactions provides direct assessments of solvent quality. The sign and magnitude of second virial coefficients will be governed by a combination of solution conditions and the nature of the macromolecule of interest. Here, we show, using a combination of theory, simulation, and in vitro experiments, that titrations of crowders, providing they are true depletants, can be used to extract the intrinsic driving forces for macromolecular phase separation. This refers to saturation concentrations in the absence of crowders and the second virial coefficients that quantify the magnitude of the incompatibility between macromolecules and the solvent. Our results show how the depletion-mediated attractions afforded by crowders can be leveraged to obtain comparative assessments of macromolecule-specific, intrinsic driving forces for phase separation.

Computational Approaches to Predict Protein-Protein Interactions in Crowded Cellular Environments

Investigating protein-protein interactions is crucial for understanding cellular biological processes because proteins often function within molecular complexes rather than in isolation. While experimental and computational methods have provided valuable insights into these interactions, they often overlook a critical factor: the crowded cellular environment. This environment significantly impacts protein behavior, including structural stability, diffusion, and ultimately the nature of binding. In this review, we discuss theoretical and computational approaches that allow the modeling of biological systems to guide and complement experiments and can thus significantly advance the investigation, and possibly the predictions, of protein-protein interactions in the crowded environment of cell cytoplasm. We explore topics such as statistical mechanics for lattice simulations, hydrodynamic interactions, diffusion processes in high-viscosity environments, and several methods based on molecular dynamics simulations. By synergistically leveraging methods from biophysics and computational biology, we review the state of the art of computational methods to study the impact of molecular crowding on protein-protein interactions and discuss its potential revolutionizing effects on the characterization of the human interactome.

Attractive crowding effect on passive and active polymer looping kinetics

Loop formation in complex environments is crucially important to many biological processes in life. In the present work, we adopt three-dimensional Langevin dynamics simulations to investigate passive and active polymer looping kinetics in crowded media featuring polymer-crowder attraction. We find polymers undergo a remarkable coil-globule-coil transition, highlighted by a marked change in the Flory scaling exponent of the gyration radius. Meanwhile, looping time as a function of the crowder's volume fraction demonstrates an apparent non-monotonic alteration. A small number of crowders induce a compact structure, which largely facilitates the looping process. While a large number of crowders heavily impede end-to-end diffusion, looping kinetics is greatly inhibited. For a self-propelled chain, we find that the attractive crowding triggers an unusual activity effect on looping kinetics. Once a globular state is formed, activity takes an effort to open the chain from the compact structure, leading to an unexpected activity-induced inhibition of looping. If the chain maintains a coil state, the dominant role of activity is to enhance diffusivity and, thus, speed up looping kinetics. The novel conformational change and looping kinetics of both passive and active polymers in the presence of attractive crowding highlight a rather distinct scenario that has no analogy in a repulsive crowding counterpart. The underlying mechanism enriches our understanding of the crucial role of attractive interactions in modulating polymer structure and dynamics.

RNA folding pathways from all-atom simulations with a variationally improved history-dependent bias

Atomically detailed simulations of RNA folding have proven very challenging in view of the difficulties of developing realistic force fields and the intrinsic computational complexity of sampling rare conformational transitions. As a step forward in tackling these issues, we extend to RNA an enhanced path-sampling method previously successfully applied to proteins. In this scheme, the information about the RNA's native structure is harnessed by a soft history-dependent biasing force promoting the generation of productive folding trajectories in an all-atom force field with explicit solvent. A rigorous variational principle is then applied to minimize the effect of the bias. Here, we report on an application of this method to RNA molecules from 20 to 47 nucleotides long and increasing topological complexity. By comparison with analog simulations performed on small proteins with similar size and architecture, we show that the RNA folding landscape is significantly more frustrated, even for relatively small chains with a simple topology. The predicted RNA folding mechanisms are found to be consistent with the available experiments and some of the existing coarse-grained models. Due to its computational performance, this scheme provides a promising platform to efficiently gather atomistic RNA folding trajectories, thus retain the information about the chemical composition of the sequence.

Crowder titrations enable the quantification of driving forces for macromolecular phase separation

Macromolecular solubility is an important contributor to the driving forces for phase separation. Formally, the driving forces in a binary mixture comprising a macromolecule dissolved in a solvent can be quantified in terms of the saturation concentration, which is the threshold macromolecular concentration above which the mixture separates into coexisting dense and dilute phases. Additionally, the second virial coefficient, which measures the effective strength of solvent-mediated intermolecular interactions provides direct assessments of solvent quality. The sign and magnitude of second virial coefficients will be governed by a combination of solution conditions and the nature of the macromolecule of interest. Here, we show, using a combination of theory, simulation, and in vitro experiments, that titrations of crowders, providing they are true depletants, can be used to extract the intrinsic driving forces for macromolecular phase separation. This refers to saturation concentrations in the absence of crowders and the second virial coefficients that quantify the magnitude of the incompatibility between macromolecules and the solvent. Our results show how the depletion-mediated attractions afforded by crowders can be leveraged to obtain comparative assessments of macromolecule-specific, intrinsic driving forces for phase separation.

SIGNIFICANCE

Phase separation has emerged as a process of significant relevance to sorting macromolecules into distinct compartments, thereby enabling spatial and temporal control over cellular matter. Considerable effort is being invested into uncovering the driving forces that enable the separation of macromolecular solutions into coexisting phases. At its heart, this process is governed by the balance of macromolecule-solvent, inter-macromolecule, and solvent-solvent interactions. We show that the driving forces for phase separation, including the coefficients that measure interaction strengths between macromolecules, can be extracted by titrating the concentrations of crowders that enable macromolecules to phase separate at lower concentrations. Our work paves the way to leverage specific categories of measurements for quantitative characterizations of driving forces for phase separation.

Odd-even disparity in the population of slipped hairpins in RNA repeat sequences with implications for phase separation

Low-complexity nucleotide repeat sequences, which are implicated in several neurological disorders, undergo liquid-liquid phase separation (LLPS) provided the number of repeat units, n, exceeds a critical value. Here, we establish a link between the folding landscapes of the monomers of trinucleotide repeats and their propensity to self-associate. Simulations using a coarse-grained Self-Organized Polymer (SOP) model for (CAG)n repeats in monovalent salt solutions reproduce experimentally measured melting temperatures, which are available only for small n. By extending the simulations to large n, we show that the free-energy gap, ΔGS, between the ground state (GS) and slipped hairpin (SH) states is a predictor of aggregation propensity. The GS for even n is a perfect hairpin (PH), whereas it is a SH when n is odd. The value of ΔGS (zero for odd n) is larger for even n than for odd n. As a result, the rate of dimer formation is slower in (CAG)30 relative to (CAG)31, thus linking ΔGS to RNA-RNA association. The yield of the dimer decreases dramatically, compared to the wild type, in mutant sequences in which the population of the SH decreases substantially. Association between RNA chains is preceded by a transition to the SH even if the GS is a PH. The finding that the excitation spectrum-which depends on the exact sequence, n, and ionic conditions-is a predictor of self-association should also hold for other RNAs (mRNA for example) that undergo LLPS.

Many Functions of Telomerase Components: Certainties, Doubts, and Inconsistencies

A growing number of studies have evidenced non-telomeric functions of "telomerase". Almost all of them, however, investigated the non-canonical effects of the catalytic subunit TERT, and not the telomerase ribonucleoprotein holoenzyme. These functions mainly comprise signal transduction, gene regulation and the increase of anti-oxidative systems. Although less studied, TERC (the RNA component of telomerase) has also been shown to be involved in gene regulation, as well as other functions. All this has led to the publication of many reviews on the subject, which, however, are often disseminating personal interpretations of experimental studies of other researchers as original proofs. Indeed, while some functions such as gene regulation seem ascertained, especially because mechanistic findings have been provided, other ones remain dubious and/or are contradicted by other direct or indirect evidence (e.g., telomerase activity at double-strand break site, RNA polymerase activity of TERT, translation of TERC, mitochondrion-processed TERC). In a critical study of the primary evidence so far obtained, we show those functions for which there is consensus, those showing contradictory results and those needing confirmation. The resulting picture, together with some usually neglected aspects, seems to indicate a link between TERT and TERC functions and cellular stemness and gives possible directions for future research.

Thermodynamics of unfolding mechanisms of mouse mammary tumor virus pseudoknot from a coarse-grained loop-entropy model

Pseudoknotted RNA molecules play important biological roles that depend on their folded structure. To understand the underlying principles that determine their thermodynamics and folding/unfolding mechanisms, we carried out a study on a variant of the mouse mammary tumor virus pseudoknotted RNA (VPK), a widely studied model system for RNA pseudoknots. Our method is based on a coarse-grained discrete-state model and the algorithm of PK3D (pseudoknot structure predictor in three-dimensional space), with RNA loops explicitly constructed and their conformational entropic effects incorporated. Our loop entropy calculations are validated by accurately capturing previously measured melting temperatures of RNA hairpins with varying loop lengths. For each of the hairpins that constitutes the VPK, we identified alternative conformations that are more stable than the hairpin structures at low temperatures and predicted their populations at different temperatures. Our predictions were validated by thermodynamic experiments on these hairpins. We further computed the heat capacity profiles of VPK, which are in excellent agreement with available experimental data. Notably, our model provides detailed information on the unfolding mechanisms of pseudoknotted RNA. Analysis of the distribution of base-pairing probability of VPK reveals a cooperative unfolding mechanism instead of a simple sequential unfolding of first one stem and then the other. Specifically, we find a simultaneous "loosening" of both stems as the temperature is raised, whereby both stems become partially melted and co-exist during the unfolding process.

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.

Computer-aided comprehensive explorations of RNA structural polymorphism through complementary simulation methods

While RNA folding was originally seen as a simple problem to solve, it has been shown that the promiscuous interactions of the nucleobases result in structural polymorphism, with several competing structures generally observed for non-coding RNA. This inherent complexity limits our understanding of these molecules from experiments alone, and computational methods are commonly used to study RNA. Here, we discuss three advanced sampling schemes, namely Hamiltonian-replica exchange molecular dynamics (MD), ratchet-and-pawl MD and discrete path sampling, as well as the HiRE-RNA coarse-graining scheme, and highlight how these approaches are complementary with reference to recent case studies. While all computational methods have their shortcomings, the plurality of simulation methods leads to a better understanding of experimental findings and can inform and guide experimental work on RNA polymorphism.

Crowding-induced polymer trapping in a channel

In this work, we report an intriguing phenomenon: crowding-induced polymer trapping in a channel. Using Langevin dynamics simulations and analytical calculations, we find that for a polymer confined in a channel, crowding particles can push a polymer into the channel corner through inducing an effective polymer-corner attraction due to the depletion effect. This phenomenon is referred to as polymer trapping. The occurrence of polymer trapping requires a minimum volume fraction of crowders, ϕ^{*}, which scales as ϕ^{*}∼(a_{c}/L_{p})^{1/3} for a_{c}≫a_{m} and ϕ^{*}∼(a_{c}/L_{p})^{1/3}(a_{c}/a_{m})^{1/2} for a_{c}≪a_{m}, where a_{c} is the crowder diameter, a_{m} is the monomer diameter, and L_{p} is the polymer persistence length. For DNA, ϕ^{*} is estimated to be around 0.25 for crowders with a_{c}=2nm. We find that ϕ^{*} also strongly depends on the shape of the channel cross section, and ϕ^{*} is much smaller for a triangle channel than a square channel. The polymer trapping leads to a nearly fully stretched polymer conformation along a channel corner, which may have practical applications, such as full stretching of DNA for the nanochannel-based genome mapping technology.

Polymer looping kinetics in active heterogeneous environments

A typical biological environment is usually featured by crowding and heterogeneity, leading to complex reaction kinetics of the immersed macromolecules. In the present work, we adopt Langevin dynamics simulations to systematically investigate polymer looping kinetics in active heterogeneous media crowded with a mixture of mobile active particles and immobile obstacles. For comparison, a parallel study is also performed in the passive heterogeneous media. We explicitly analyze the change of looping time and looping probability with the variation of obstacle ratio, volume fraction and crowder size. We reveal the novel phenomena of inhibition-facilitation transition of the looping rate induced by heterogeneity, crowdedness and activity. In addition, our results demonstrate a very non-trivial crowder size effect on the looping kinetics. The underlying mechanism is rationalized by the interplay of polymer diffusion, conformational change and looping free-energy barrier. The competing effect arising from active particles and obstacles on structural and dynamical properties of the polymer yields a consistent scenario for our observations. Lastly, the non-exponential kinetics of the looping process is also analyzed. We find that both activity and crowding can strengthen the heterogeneity degree of the looping kinetics.

Influence of solvent quality on depletion potentials in colloid-polymer mixtures

As first explained by the classic Asakura-Oosawa (AO) model, effective attractive forces between colloidal particles induced by depletion of nonadsorbing polymers can drive demixing of colloid-polymer mixtures into colloid-rich and colloid-poor phases, with practical relevance for purification of water, stability of foods and pharmaceuticals, and macromolecular crowding in biological cells. By idealizing polymer coils as effective penetrable spheres, the AO model qualitatively captures the influence of polymer depletion on thermodynamic phase behavior of colloidal suspensions. In previous work, we extended the AO model to incorporate aspherical polymer conformations and showed that fluctuating shapes of random-walk coils can significantly modify depletion potentials [W. K. Lim and A. R. Denton, Soft Matter 12, 2247 (2016); J. Chem. Phys. 144, 024904 (2016)]. We further demonstrated that the shapes of polymers in crowded environments sensitively depend on solvent quality [W. J. Davis and A. R. Denton, J. Chem. Phys. 149, 124901 (2018)]. Here, we apply Monte Carlo simulation to analyze the influence of solvent quality on depletion potentials in mixtures of hard-sphere colloids and nonadsorbing polymer coils, modeled as ellipsoids whose principal radii fluctuate according to random-walk statistics. We consider both self-avoiding and non-self-avoiding random walks, corresponding to polymers in good and theta solvents, respectively. Our simulation results demonstrate that depletion of polymers of equal molecular weight induces much stronger attraction between colloids in good solvents than in theta solvents and confirm that depletion interactions are significantly influenced by aspherical polymer conformations.

Polymer brush-induced depletion interactions and clustering of membrane proteins

We investigate the effect of mobile polymer brushes on proteins embedded in biological membranes by employing both Asakura-Oosawa type of theoretical model and coarse-grained molecular dynamics simulations. The brush polymer-induced depletion attraction between proteins changes non-monotonically with the size of brush. The depletion interaction, which is determined by the ratio of the protein size to the grafting distance between brush polymers, increases linearly with the brush size as long as the polymer brush height is shorter than the protein size. When the brush height exceeds the protein size, however, the depletion attraction among proteins is slightly reduced. We also explore the possibility of the brush polymer-induced assembly of a large protein cluster, which can be related to one of many molecular mechanisms underlying recent experimental observations of integrin nanocluster formation and signaling.

Shape changes and cooperativity in the folding of the central domain of the 16S ribosomal RNA

Both the small and large subunits of the ribosome, the molecular machine that synthesizes proteins, are complexes of ribosomal RNAs (rRNAs) and a number of proteins. In bacteria, the small subunit has a single 16S rRNA whose folding is the first step in its assembly. The central domain of the 16S rRNA folds independently, driven either by Mg²⁺ ions or by interaction with ribosomal proteins. To provide a quantitative description of ion-induced folding of the ∼350-nucleotide rRNA, we carried out extensive coarse-grained molecular simulations spanning Mg²⁺ concentration between 0 and 30 mM. The Mg²⁺ dependence of the radius of gyration shows that globally the rRNA folds cooperatively. Surprisingly, various structural elements order at different Mg²⁺ concentrations, indicative of the heterogeneous assembly even within a single domain of the rRNA. Binding of Mg²⁺ ions is highly specific, with successive ion condensation resulting in nucleation of tertiary structures. We also predict the Mg²⁺-dependent protection factors, measurable in hydroxyl radical footprinting experiments, which corroborate the specificity of Mg²⁺-induced folding. The simulations, which agree quantitatively with several experiments on the folding of a three-way junction, show that its folding is preceded by formation of other tertiary contacts in the central junction. Our work provides a starting point in simulating the early events in the assembly of the small subunit of the ribosome.

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.

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.

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.

RNA contributions to the form and function of biomolecular condensates

Biomolecular condensation partitions cellular contents and has important roles in stress responses, maintaining homeostasis, development and disease. Many nuclear and cytoplasmic condensates are rich in RNA and RNA-binding proteins (RBPs), which undergo liquid-liquid phase separation (LLPS). Whereas the role of RBPs in condensates has been well studied, less attention has been paid to the contribution of RNA to LLPS. In this Review, we discuss the role of RNA in biomolecular condensation and highlight considerations for designing condensate reconstitution experiments. We focus on RNA properties such as composition, length, structure, modifications and expression level. These properties can modulate the biophysical features of native condensates, including their size, shape, viscosity, liquidity, surface tension and composition. We also discuss the role of RNA-protein condensates in development, disease and homeostasis, emphasizing how their properties and function can be determined by RNA. Finally, we discuss the multifaceted cellular functions of biomolecular condensates, including cell compartmentalization through RNA transport and localization, supporting catalytic processes, storage and inheritance of specific molecules, and buffering noise and responding to stress.

Direct observation of effect of crowding induced macromolecular hydration on molecular breathing in the stem of Fork-DNA by single-molecule FRET microspectroscopy

The perpetually changing cellular conditions, nucleotide sequence, and environmental effects including osmotic stress have multiple effects on DNA, leading to several conformational alternations and subsequently influencing their activity, too. In this work, single-molecule FRET microspectroscopy has been employed to monitor the breathing dynamics as an effect of molecular crowding in the stem region of Fork-DNA. The structural integrity greatly alters with the presence or absence of nucleotide overhangs and on the nature and concentration of the crowding agent, thus affecting the stability of the stem region and hence the forked DNA. The multiple hydrogen bonds and hydrophobic interactions between the polynucleotide strands appear to be altered with osmotic crowding. This induces increased flexibility in the double helix and allows DNA to breath. The conformational alternation of the DNA happens in nanometer resolution, that is been monitored by the change in the FRET efficiency between the dyes attached to two different strands of the DNA. The nature and molecular weight of crowding agents control the degree of spatial breathing in the stem of Fork-DNA. These constant fluctuations between the entropically favorable partially folded structures to an enthalpically favorable folded structure are not only valuable for elucidating nucleic acid structure but might play an important role in enzyme kinetics.

Understanding DNA interactions in crowded environments with a coarse-grained model

Nucleic acid interactions under crowded environments are of great importance for biological processes and nanotechnology. However, the kinetics and thermodynamics of nucleic acid interactions in a crowded environment remain poorly understood. We use a coarse-grained model of DNA to study the kinetics and thermodynamics of DNA duplex and hairpin formation in crowded environments. We find that crowders can increase the melting temperature of both an 8-mer DNA duplex and a hairpin with a stem of 6-nt depending on the excluded volume fraction of crowders in solution and the crowder size. The crowding induced stability originates from the entropic effect caused by the crowding particles in the system. Additionally, we study the hybridization kinetics of DNA duplex formation and the formation of hairpin stems, finding that the reaction rate k_on is increased by the crowding effect, while k_off is changed only moderately. The increase in k_on mostly comes from increasing the probability of reaching a transition state with one base pair formed. A DNA strand displacement reaction in a crowded environment is also studied with the model and we find that rate of toehold association is increased, with possible applications to speeding up strand displacement cascades in nucleic acid nanotechnology.

Molecular crowding and RNA catalysis

RNA enzymes or ribozymes catalyze some of the most important reactions in biology and are thought to have played a central role in the origin and evolution of life on earth. Catalytic function in RNA has evolved in crowded cellular environments that are different from dilute solutions in which most in vitro assays are performed. The presence of molecules such as amino acids, polypeptides, alcohols, and sugars in the cell introduces forces that modify the kinetics and thermodynamics of ribozyme-catalyzed reactions. Synthetic molecules are routinely used in in vitro studies to better approximate the properties of biomolecules under in vivo conditions. This review discusses the various forces that operate within simulated crowded solutions in the context of RNA structure, folding, and catalysis. It also explores ideas about how crowding could have been beneficial to the evolution of functional RNAs and the development of primitive cellular systems in a prebiotic milieu.

Crowding-Activity Coupling Effect on Conformational Change of a Semi-Flexible Polymer

The behavior of a polymer in a passive crowded medium or in a very dilute active bath has been well studied, while a polymer immersed in an environment featured by both crowding and activity remains an open problem. In this paper, a systematic Langevin simulation is performed to investigate the conformational change of a semi-flexible chain in a concentrated solution packed with spherical active crowders. A very novel shrinkage-to-swelling transition is observed for a polymer with small rigidity. The underlying phase diagram is constructed in the parameter space of active force and crowder size. Moreover, the variation of the polymer gyration radius demonstrates a non-monotonic dependence on the dynamical persistence length of the active particle. Lastly, the activity-crowding coupling effect in different crowder size baths is clarified. In the case of small crowders, activity strengthens the crowding-induced shrinkage to the chain. As crowder size increases, activity turns out to be a contrasting factor to crowding, resulting in a competitive shrinkage and swelling. In the large size situation, the swelling effect arising from activity eventually becomes dominant. The present study provides a deeper understanding of the unusual behavior of a semi-flexible polymer in an active and crowded medium, associated with the nontrivial activity-crowding coupling and the cooperative crowder size effect.

Comparative study of the crowding-induced collapse effect in hard-sphere, flexible polymer and rod-like polymer systems

A systematic Langevin simulation is performed to study the crowding-induced collapse effect on a probed chain in three typical systems: hard sphere (HS), flexible polymer and rod-like polymer.

Optimal molecular crowding accelerates group II intron folding and maximizes catalysis

Unlike in vivo conditions, group II intron ribozymes are known to require high magnesium(II) concentrations ([Mg2+]) and high temperatures (42 °C) for folding and catalysis in vitro. A possible explanation for this difference is the highly crowded cellular environment, which can be mimicked in vitro by macromolecular crowding agents. Here, we combined bulk activity assays and single-molecule Förster Resonance Energy Transfer (smFRET) to study the influence of polyethylene glycol (PEG) on catalysis and folding of the ribozyme. Our activity studies reveal that PEG reduces the [Mg2+] required, and we found an "optimum" [PEG] that yields maximum activity. smFRET experiments show that the most compact state population, the putative active state, increases with increasing [PEG]. Dynamic transitions between folded states also increase. Therefore, this study shows that optimal molecular crowding concentrations help the ribozyme not only to reach the native fold but also to increase its in vitro activity to approach that in physiological conditions.

Understanding chain looping kinetics in polymer solutions: crowding effects of microviscosity and collapse

A theoretical framework based on a generalized Langevin equation with fractional Gaussian noise is presented to describe the looping kinetics of chains in polymer solutions. Particular attention is paid to quantitatively revealing crowding effects on the loop formation rate in terms of microviscosity and collapse. By the aid of empirical relations for these two crowding associated physical quantities, we explicitly investigate the relationship between the looping rate and polymer concentration, the degree of polymerization, and system parameters. According to our analysis, the dependence of the looping rate on the crowder volume fraction exhibits three typical regimes: monotonic decreasing, a non-monotonic trend and monotonic increasing. We reveal that these non-trivial behaviors can be attributed to the competition between the two opposing factors of viscosity-associated inhibition and collapse-induced facilitation of loop formation. We apply our theory to analyze the kinetics of single-stranded DNA hairpin base pairing in polyethylene glycol solutions. The theoretical results can reproduce the experimental data on the closing rate of hairpins quantitatively to a certain degree with reasonable fitting parameters. The unexpected increase of the closing rate upon the addition of increasing amounts of polymer is well rationalised. Such good agreements clearly demonstrate the validity of our theory, appropriately addressing the very role of crowding effects in the relevant kinetics.

Influence of solvent quality on conformations of crowded polymers

The structure and function of polymers in confined environments, e.g., biopolymers in the cytoplasm of a cell, are strongly affected by macromolecular crowding. To explore the influence of solvent quality on conformations of crowded polymers, we model polymers as penetrable ellipsoids, whose shape fluctuations are governed by the statistics of self-avoiding walks, appropriate for a polymer in a good solvent. Within this coarse-grained model, we perform Monte Carlo simulations of mixtures of polymers and hard-nanosphere crowders, including trial changes in polymer size and shape. Penetration of polymers by crowders is incorporated via a free energy cost predicted by polymer field theory. To analyze the impact of crowding on polymer conformations in different solvents, we compute the average polymer shape distributions, radius of gyration, volume, and asphericity over ranges of the polymer-to-crowder size ratio and crowder volume fraction. The simulation results are accurately predicted by a free-volume theory of polymer crowding. Comparison of results for polymers in good and theta solvents indicates that excluded-volume interactions between polymer segments significantly affect crowding, especially in the limit of crowders much smaller than polymers. Our approach may help to motivate future experimental studies of polymers in crowded environments, with possible relevance for drug delivery and gene therapy.

Frictional Effects on RNA Folding: Speed Limit and Kramers Turnover

We investigated frictional effects on the folding rates of a human telomerase hairpin (hTR HP) and H-type pseudoknot from the Beet Western Yellow Virus (BWYV PK) using simulations of the Three Interaction Site (TIS) model for RNA. The heat capacity from TIS model simulations, calculated using temperature replica exchange simulations, reproduces nearly quantitatively the available experimental data for the hTR HP. The corresponding results for BWYV PK serve as predictions. We calculated the folding rates ( k_F) from more than 100 folding trajectories for each value of the solvent viscosity (η) at a fixed salt concentration of 200 mM. By using the theoretical estimate (∝ √N where N is the number of nucleotides) for folding free energy barrier, k_F data for both the RNAs are quantitatively fit using one-dimensional Kramers's theory with two parameters specifying the curvatures in the unfolded basin and the barrier top. In the high-friction regime (η ≳ 10^-5 Pa·s), for both HP and PK, k_F values decrease as 1/η, whereas in the low friction regime, k_F values increase as η increases, leading to a maximum folding rate at a moderate viscosity (∼10^-6 Pa·s), which is the Kramers turnover. From the fits, we find that the speed limit to RNA folding at water viscosity is between 1 and 4 μs, which is in accord with our previous theoretical prediction as well as results from several single molecule experiments. Both the RNA constructs fold by parallel pathways. Surprisingly, we find that the flux through the pathways could be altered by changing solvent viscosity, a prediction that is more easily testable in RNA than in proteins.

Effects of flanking regions on HDV cotranscriptional folding kinetics

Hepatitis delta virus (HDV) ribozyme performs the self-cleavage activity through folding to a double pseudoknot structure. The folding of functional RNA structures is often coupled with the transcription process. In this work, we developed a new approach for predicting the cotranscriptional folding kinetics of RNA secondary structures with pseudoknots. We theoretically studied the cotranscriptional folding behavior of the 99-nucleotide (nt) HDV sequence, two upstream flanking sequences, and one downstream flanking sequence. During transcription, the 99-nt HDV can effectively avoid the trap intermediates and quickly fold to the cleavage-active state. It is different from its refolding kinetics, which folds into an intermediate trap state. For all the sequences, the ribozyme regions (from 1 to 73) all fold to the same structure during transcription. However, the existence of the 30-nt upstream flanking sequence can inhibit the ribozyme region folding into the active native state through forming an alternative helix Alt1 with the segments 70-90. The longer upstream flanking sequence of 54 nt itself forms a stable hairpin structure, which sequesters the formation of the Alt1 helix and leads to rapid formation of the cleavage-active structure. Although the 55-nt downstream flanking sequence could invade the already folded active structure during transcription by forming a more stable helix with the ribozyme region, the slow transition rate could keep the structure in the cleavage-active structure to perform the activity.

Sequence-Dependent Three Interaction Site Model for Single- and Double-Stranded DNA

We develop a robust coarse-grained model for single- and double-stranded DNA by representing each nucleotide by three interaction sites (TIS) located at the centers of mass of sugar, phosphate, and base. The resulting TIS model includes base-stacking, hydrogen bond, and electrostatic interactions as well as bond-stretching and bond angle potentials that account for the polymeric nature of DNA. The choices of force constants for stretching and the bending potentials were guided by a Boltzmann inversion procedure using a large representative set of DNA structures extracted from the Protein Data Bank. Some of the parameters in the stacking interactions were calculated using a learning procedure, which ensured that the experimentally measured melting temperatures of dimers are faithfully reproduced. Without any further adjustments, the calculations based on the TIS model reproduce the experimentally measured salt and sequence-dependence of the size of single-stranded DNA (ssDNA), as well as the persistence lengths of poly(dA) and poly(dT) chains. Interestingly, upon application of mechanical force, the extension of poly(dA) exhibits a plateau, which we trace to the formation of stacked helical domains. In contrast, the force-extension curve (FEC) of poly(dT) is entropic in origin and could be described by a standard polymer model. We also show that the persistence length of double-stranded DNA, formed from two complementary ssDNAs, is consistent with the prediction based on the worm-like chain. The persistence length, which decreases with increasing salt concentration, is in accord with the Odijk-Skolnick-Fixman theory intended for stiff polyelectrolyte chains near the rod limit. Our model predicts the melting temperatures of DNA hairpins with excellent accuracy, and we are able to recover the experimentally known sequence-specific trends. The range of applications, which did not require adjusting any parameter after the initial construction based solely on PDB structures and melting profiles of dimers, attests to the transferability and robustness of the TIS model for ssDNA and dsDNA.

Lipid vesicles chaperone an encapsulated RNA aptamer

The organization of molecules into cells is believed to have been critical for the emergence of living systems. Early protocells likely consisted of RNA functioning inside vesicles made of simple lipids. However, little is known about how encapsulation would affect the activity and folding of RNA. Here we find that confinement of the malachite green RNA aptamer inside fatty acid vesicles increases binding affinity and locally stabilizes the bound conformation of the RNA. The vesicle effectively ‘chaperones’ the aptamer, consistent with an excluded volume mechanism due to confinement. Protocellular organization thereby leads to a direct benefit for the RNA. Coupled with previously described mechanisms by which encapsulated RNA aids membrane growth, this effect illustrates how the membrane and RNA might cooperate for mutual benefit. Encapsulation could thus increase RNA fitness and the likelihood that functional sequences would emerge during the origin of life.

So far little is known about how encapsulation affects the activity and folding of RNA, which is of interest for understanding the origin of cellular life. Here the authors show that encapsulation of functional RNA in vesicles increases RNA activity and improves RNA folding through a biophysical confinement effect.

Markov-state model for CO2 binding with carbonic anhydrase under confinement

Enzyme immobilization with a nanostructure material can enhance its stability and facilitate reusability. However, the apparent activity is often compromised due to additional diffusion barriers and complex interactions with the substrates and solvent molecules. The present study elucidates the effects of the surface hydrophobicity of nano-confinement on CO₂ diffusion to the active site of human carbonic anhydrase II (CA), an enzyme that is able to catalyze CO₂ hydration at extremely high turnover rates. Using the Markov-state model in combination with coarse-grained molecular dynamics simulations, we demonstrate that a hydrophobic cage increases CO₂ local density but hinders its diffusion towards the active site of CA under confinement. By contrast, a hydrophilic cage hinders CO₂ adsorption but promotes its binding with CA. An optimal surface hydrophobicity can be identified to maximize both the CO₂ occupation probability and the diffusion rate. The simulation results offer insight into understanding enzyme performance under nano-confinement and help us to advance broader applications of CA for CO₂ absorption and recovery.

Predicting 3D structure and stability of RNA pseudoknots in monovalent and divalent ion solutions

RNA pseudoknots are a kind of minimal RNA tertiary structural motifs, and their three-dimensional (3D) structures and stability play essential roles in a variety of biological functions. Therefore, to predict 3D structures and stability of RNA pseudoknots is essential for understanding their functions. In the work, we employed our previously developed coarse-grained model with implicit salt to make extensive predictions and comprehensive analyses on the 3D structures and stability for RNA pseudoknots in monovalent/divalent ion solutions. The comparisons with available experimental data show that our model can successfully predict the 3D structures of RNA pseudoknots from their sequences, and can also make reliable predictions for the stability of RNA pseudoknots with different lengths and sequences over a wide range of monovalent/divalent ion concentrations. Furthermore, we made comprehensive analyses on the unfolding pathway for various RNA pseudoknots in ion solutions. Our analyses for extensive pseudokonts and the wide range of monovalent/divalent ion concentrations verify that the unfolding pathway of RNA pseudoknots is mainly dependent on the relative stability of unfolded intermediate states, and show that the unfolding pathway of RNA pseudoknots can be significantly modulated by their sequences and solution ion conditions.

RNA Structural Dynamics As Captured by Molecular Simulations: A Comprehensive Overview

With both catalytic and genetic functions, ribonucleic acid (RNA) is perhaps the most pluripotent chemical species in molecular biology, and its functions are intimately linked to its structure and dynamics. Computer simulations, and in particular atomistic molecular dynamics (MD), allow structural dynamics of biomolecular systems to be investigated with unprecedented temporal and spatial resolution. We here provide a comprehensive overview of the fast-developing field of MD simulations of RNA molecules. We begin with an in-depth, evaluatory coverage of the most fundamental methodological challenges that set the basis for the future development of the field, in particular, the current developments and inherent physical limitations of the atomistic force fields and the recent advances in a broad spectrum of enhanced sampling methods. We also survey the closely related field of coarse-grained modeling of RNA systems. After dealing with the methodological aspects, we provide an exhaustive overview of the available RNA simulation literature, ranging from studies of the smallest RNA oligonucleotides to investigations of the entire ribosome. Our review encompasses tetranucleotides, tetraloops, a number of small RNA motifs, A-helix RNA, kissing-loop complexes, the TAR RNA element, the decoding center and other important regions of the ribosome, as well as assorted others systems. Extended sections are devoted to RNA-ion interactions, ribozymes, riboswitches, and protein/RNA complexes. Our overview is written for as broad of an audience as possible, aiming to provide a much-needed interdisciplinary bridge between computation and experiment, together with a perspective on the future of the field.

Theory and Modeling of RNA Structure and Interactions with Metal Ions and Small Molecules

In addition to continuous rapid progress in RNA structure determination, probing, and biophysical studies, the past decade has seen remarkable advances in the development of a new generation of RNA folding theories and models. In this article, we review RNA structure prediction models and models for ion-RNA and ligand-RNA interactions. These new models are becoming increasingly important for a mechanistic understanding of RNA function and quantitative design of RNA nanotechnology. We focus on new methods for physics-based, knowledge-based, and experimental data-directed modeling for RNA structures and explore the new theories for the predictions of metal ion and ligand binding sites and metal ion-dependent RNA stabilities. The integration of these new methods with theories about the cellular environment effects in RNA folding, such as molecular crowding and cotranscriptional kinetic effects, may ultimately lead to an all-encompassing RNA folding model.

Combining Diffusion NMR and Small-Angle Neutron Scattering Enables Precise Measurements of Polymer Chain Compression in a Crowded Environment

The effect of particles on the behavior of polymers in solution is important in a number of important phenomena such as the effect of "crowding" proteins in cells, colloid-polymer mixtures, and nanoparticle "fillers" in polymer solutions and melts. In this Letter, we study the effect of spherical inert nanoparticles (which we refer to as "crowders") on the diffusion coefficient and radius of gyration of polymers in solution using pulsed-field-gradient NMR and small-angle neutron scattering (SANS), respectively. The diffusion coefficients exhibit a plateau below a characteristic polymer concentration, which we identify as the overlap threshold concentration c^{⋆}. Above c^{⋆}, in a crossover region between the dilute and semidilute regimes, the (long-time) self-diffusion coefficients are found, universally, to decrease exponentially with polymer concentration at all crowder packing fractions, consistent with a structural basis for the long-time dynamics. The radius of gyration obtained from SANS in the crossover regime changes linearly with an increase in polymer concentration, and must be extrapolated to c^{⋆} in order to obtain the radius of gyration of an individual polymer chain. When the polymer radius of gyration and crowder size are comparable, the polymer size is very weakly affected by the presence of crowders, consistent with recent computer simulations. There is significant chain compression, however, when the crowder size is much smaller than the polymer radius gyration.

Predicting Molecular Crowding Effects in Ion-RNA Interactions

We develop a new statistical mechanical model to predict the molecular crowding effects in ion-RNA interactions. By considering discrete distributions of the crowders, the model can treat the main crowder-induced effects, such as the competition with ions for RNA binding, changes of electrostatic interaction due to crowder-induced changes in the dielectric environment, and changes in the nonpolar hydration state of the crowder-RNA system. To enhance the computational efficiency, we sample the crowder distribution using a hybrid approach: For crowders in the close vicinity of RNA surface, we sample their discrete distributions; for crowders in the bulk solvent away from the RNA surface, we use a continuous mean-field distribution for the crowders. Moreover, using the tightly bound ion (TBI) model, we account for ion fluctuation and correlation effects in the calculation for ion-RNA interactions. Applications of the model to a variety of simple RNA structures such as RNA helices show a crowder-induced increase in free energy and decrease in ion binding. Such crowding effects tend to contribute to the destabilization of RNA structure. Further analysis indicates that these effects are associated with the crowder-ion competition in RNA binding and the effective decrease in the dielectric constant. This simple ion effect model may serve as a useful framework for modeling more realistic crowders with larger, more complex RNA structures.

Non-Brownian diffusion in lipid membranes: Experiments and simulations

The dynamics of constituents and the surface response of cellular membranes-also in connection to the binding of various particles and macromolecules to the membrane-are still a matter of controversy in the membrane biophysics community, particularly with respect to crowded membranes of living biological cells. We here put into perspective recent single particle tracking experiments in the plasma membranes of living cells and supercomputing studies of lipid bilayer model membranes with and without protein crowding. Special emphasis is put on the observation of anomalous, non-Brownian diffusion of both lipid molecules and proteins embedded in the lipid bilayer. While single component, pure lipid bilayers in simulations exhibit only transient anomalous diffusion of lipid molecules on nanosecond time scales, the persistence of anomalous diffusion becomes significantly longer ranged on the addition of disorder-through the addition of cholesterol or proteins-and on passing of the membrane lipids to the gel phase. Concurrently, experiments demonstrate the anomalous diffusion of membrane embedded proteins up to macroscopic time scales in the minute time range. Particular emphasis will be put on the physical character of the anomalous diffusion, in particular, the occurrence of ageing observed in the experiments-the effective diffusivity of the measured particles is a decreasing function of time. Moreover, we present results for the time dependent local scaling exponent of the mean squared displacement of the monitored particles. Recent results finding deviations from the commonly assumed Gaussian diffusion patterns in protein crowded membranes are reported. The properties of the displacement autocorrelation function of the lipid molecules are discussed in the light of their appropriate physical anomalous diffusion models, both for non-crowded and crowded membranes. In the last part of this review we address the upcoming field of membrane distortion by elongated membrane-binding particles. We discuss how membrane compartmentalisation and the particle-membrane binding energy may impact the dynamics and response of lipid membranes. This article is part of a Special Issue entitled: Biosimulations edited by Ilpo Vattulainen and Tomasz Róg.

Reconstituting Intracellular Vesicle Fusion Reactions: The Essential Role of Macromolecular Crowding

Intracellular vesicle fusion is mediated by SNAREs and Sec1/Munc18 (SM) proteins. Despite intensive efforts, the SNARE-SM mediated vesicle fusion reaction has not been faithfully reconstituted in biochemical assays. Here, we present an unexpected discovery that macromolecular crowding is required for reconstituting the vesicle fusion reaction in vitro. Macromolecular crowding is known to profoundly influence the kinetic and thermodynamic behaviors of macromolecules, but its role in membrane transport processes such as vesicle fusion remains unexplored. We introduced macromolecular crowding agents into reconstituted fusion reactions to mimic the crowded cellular environment. In this crowded assay, SNAREs and SM proteins acted in concert to drive efficient membrane fusion. In uncrowded assays, by contrast, SM proteins failed to associate with the SNAREs and the fusion rate decreased more than 30-fold, close to undetectable levels. The activities of SM proteins were strictly specific to their cognate SNARE isoforms and sensitive to biologically relevant mutations, further supporting that the crowded fusion assay accurately recapitulates the vesicle fusion reaction. Using this crowded fusion assay, we also showed that the SNARE-SM mediated fusion reaction can be modulated by two additional factors: NSF and α-SNAP. These findings suggest that the vesicle fusion machinery likely has been evolutionarily selected to function optimally in the crowded milieu of the cell. Accordingly, macromolecular crowding should constitute an integral element of any reconstituted fusion assay.

How do metal ions direct ribozyme folding?

Ribozymes, which carry out phosphoryl-transfer reactions, often require Mg(2+) ions for catalytic activity. The correct folding of the active site and ribozyme tertiary structure is also regulated by metal ions in a manner that is not fully understood. Here we employ coarse-grained molecular simulations to show that individual structural elements of the group I ribozyme from the bacterium Azoarcus form spontaneously in the unfolded ribozyme even at very low Mg(2+) concentrations, and are transiently stabilized by the coordination of Mg(2+) ions to specific nucleotides. However, competition for scarce Mg(2+) and topological constraints that arise from chain connectivity prevent the complete folding of the ribozyme. A much higher Mg(2+) concentration is required for complete folding of the ribozyme and stabilization of the active site. When Mg(2+) is replaced by Ca(2+) the ribozyme folds, but the active site remains unstable. Our results suggest that group I ribozymes utilize the same interactions with specific metal ligands for both structural stability and chemical activity.

Unexpected Swelling of Stiff DNA in a Polydisperse Crowded Environment

We investigate the conformations of DNA-like stiff chains, characterized by contour length (L) and persistence length (lp), in a variety of crowded environments containing monodisperse soft spherical (SS) and spherocylindrical (SC) particles, a mixture of SS and SC, and a milieu mimicking the composition of proteins in the Escherichia coli cytoplasm. The stiff chain, whose size modestly increases in SS crowders up to ϕ ≈ 0.1, is considerably more compact at low volume fractions (ϕ ≤ 0.2) in monodisperse SC particles than in a medium containing SS particles. A 1:1 mixture of SS and SC crowders induces greater chain compaction than the pure SS or SC crowders at the same ϕ, with the effect being highly nonadditive. We also discover a counterintuitive result that the polydisperse crowding environment, mimicking the composition of a cell lysate, swells the DNA-like polymer, which is in stark contrast to the size reduction of flexible polymers in the same milieu. Trapping of the stiff chain in a fluctuating tube-like environment created by large-sized crowders explains the dramatic increase in size and persistence length of the stiff chain. In the polydisperse medium, mimicking the cellular environment, the size of the DNA (or related RNA) is determined by L/lp. At low L/lp, the size of the polymer is unaffected, whereas there is a dramatic swelling at an intermediate value of L/lp. We use these results to provide insights into recent experiments on crowding effects on RNA and also make testable predictions.

Free energy landscape and multiple folding pathways of an H-type RNA pseudoknot

How RNA sequences fold to specific tertiary structures is one of the key problems for understanding their dynamics and functions. Here, we study the folding process of an H-type RNA pseudoknot by performing a large-scale all-atom MD simulation and bias-exchange metadynamics. The folding free energy landscapes are obtained and several folding intermediates are identified. It is suggested that the folding occurs via multiple mechanisms, including a step-wise mechanism starting either from the first helix or the second, and a cooperative mechanism with both helices forming simultaneously. Despite of the multiple mechanism nature, the ensemble folding kinetics estimated from a Markov state model is single-exponential. It is also found that the correlation between folding and binding of metal ions is significant, and the bound ions mediate long-range interactions in the intermediate structures. Non-native interactions are found to be dominant in the unfolded state and also present in some intermediates, possibly hinder the folding process of the RNA.

Polymer Looping Is Controlled by Macromolecular Crowding, Spatial Confinement, and Chain Stiffness

We study by extensive computer simulations the looping characteristics of linear polymers with varying persistence length inside a spherical cavity in the presence of macromolecular crowding. For stiff chains, the looping probability and looping time reveal wildly oscillating patterns as functions of the chain length. The effects of crowding differ dramatically for flexible versus stiff polymers. While for flexible chains the looping kinetics is slowed down by the crowders, for stiffer chains the kinetics turns out to be either decreased or facilitated, depending on the polymer length. For severe confinement, the looping kinetics may become strongly facilitated by crowding. Our findings are of broad impact for DNA looping in the crowded and compartmentalized interior of living biological cells.

Effects of macromolecular crowding on the collapse of biopolymers

Experiments show that macromolecular crowding modestly reduces the size of intrinsically disordered proteins even at a volume fraction (ϕ) similar to that in the cytosol, whereas DNA undergoes a coil-to-globule transition at very small ϕ. We show using a combination of scaling arguments and simulations that the polymer size R̅(g)(ϕ) depends on x=R̅(g)(0)/D, where D is the ϕ-dependent distance between the crowders. If x≲O(1), there is only a small decrease in R̅(g)(ϕ) as ϕ increases. When x≫O(1), a cooperative coil-to-globule transition is induced. Our theory quantitatively explains a number of experiments.

Kinetics of polymer looping with macromolecular crowding: effects of volume fraction and crowder size

The looping of polymers such as DNA is a fundamental process in the molecular biology of living cells, whose interior is characterised by a high degree of molecular crowding. We here investigate in detail the looping dynamics of flexible polymer chains in the presence of different degrees of crowding. From the analysis of the looping-unlooping rates and the looping probabilities of the chain ends we show that the presence of small crowders typically slows down the chain dynamics but larger crowders may in fact facilitate the looping. We rationalise these non-trivial and often counterintuitive effects of the crowder size on the looping kinetics in terms of an effective solution viscosity and standard excluded volume. It is shown that for small crowders the effect of an increased viscosity dominates, while for big crowders we argue that confinement effects (caging) prevail. The tradeoff between both trends can thus result in the impediment or facilitation of polymer looping, depending on the crowder size. We also examine how the crowding volume fraction, chain length, and the attraction strength of the contact groups of the polymer chain affect the looping kinetics and hairpin formation dynamics. Our results are relevant for DNA looping in the absence and presence of protein mediation, DNA hairpin formation, RNA folding, and the folding of polypeptide chains under biologically relevant high-crowding conditions.