Effect of macromolecular crowding on reaction rates: a computational and theoretical study

Size Sensitivity of Metabolite Diffusion in Macromolecular Crowds

Metabolites play crucial roles in cellular processes, yet their diffusion in the densely packed interiors of cells remains poorly understood, compounded by conflicting reports in existing studies. Here, we employ pulsed-gradient stimulated-echo NMR and Brownian/Stokesian dynamics simulations to elucidate the behavior of nano- and subnanometer-sized tracers in crowded environments. Using Ficoll as a crowder, we observe a linear decrease in tracer diffusivity with increasing occupied volume fraction, persisting─somewhat surprisingly─up to volume fractions of 30-40%. While simulations suggest a linear correlation between diffusivity slowdown and particle size, experimental findings hint at a more intricate relationship, possibly influenced by Ficoll's porosity. Simulations and numerical calculations of tracer diffusivity in the E. coli cytoplasm show a nonlinear yet monotonic diffusion slowdown with particle size. We discuss our results in the context of nanoviscosity and discrepancies with existing studies.

Kinetic trapping organizes actin filaments within liquid-like protein droplets

Several actin-binding proteins (ABPs) phase separate to form condensates capable of curating the actin network shapes. Here, we use computational modeling to understand the principles of actin network organization within VASP condensate droplets. Our simulations reveal that the different actin shapes, namely shells, rings, and mixture states are highly dependent on the kinetics of VASP-actin interactions, suggesting that they arise from kinetic trapping. Specifically, we show that reducing the residence time of VASP on actin filaments reduces degree of bundling, thereby promoting assembly of shells rather than rings. We validate the model predictions experimentally using a VASP-mutant with decreased bundling capability. Finally, we investigate the ring opening within deformed droplets and found that the sphere-to-ellipsoid transition is favored under a wide range of filament lengths while the ellipsoid-to-rod transition is only permitted when filaments have a specific range of lengths. Our findings highlight key mechanisms of actin organization within phase-separated ABPs.

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.

Nonmonotonic Modulation of the Protein-Ligand Recognition Event by Inert Crowders

The ubiquitous event of a protein recognizing small molecules or ligands at its native binding site is crucial for initiating major biological processes. However, how a crowded environment, as is typically represented by a cellular interior, would modulate the protein-ligand search process is largely debated. Excluded volume-based theory suggests that the presence of an inert crowder would reinforce a steady stabilization and enhancement of the protein-ligand recognition process. Here, we counter this long-held perspective via the molecular dynamics simulation and Markov state model of the protein-ligand recognition event in the presence of inert crowders. Specifically, we demonstrate that, depending on concentration, even purely inert crowders can exert a nonmonotonic effect via either stabilizing or destabilizing the protein-ligand binding event. Analysis of the kinetic network of binding pathways reveals that the crowders would either modulate precedent non-native on-pathway intermediates or would devise additional ones in a multistate recognition event across a wide range of concentrations. As an important insight, crowders gradually shift the relative transitional preference of these intermediates toward a native-bound state, with ligand residence time at the binding pocket dictating the trend of nonmonotonic concentration dependence by simple inert crowders.

Molecular Crowders Can Induce Collapse in Hydrophilic Polymers via Soft Attractive Interactions

A comprehensive understanding of protein folding and biomolecular self-assembly in the intracellular environment requires obtaining a microscopic view of the crowding effects. The classical view of crowding explains biomolecular collapse in such an environment in terms of the entropic solvent excluded volume effects subjected to hard-core repulsions exerted by the inert crowders, neglecting their soft chemical interactions. In this study, the effects of nonspecific, soft interactions of molecular crowders in regulating the conformational equilibrium of hydrophilic (charged) polymers are examined. Using advanced molecular dynamics simulations, collapse free energies of an uncharged, a negatively charged, and a charge-neutral 32-mer generic polymer are computed. The strength of the polymer-crowder dispersion energy is modulated to examine its effect on polymer collapse. The results show that the crowders preferentially adsorb and drive the collapse of all three polymers. The uncharged polymer collapse is opposed by the change in solute-solvent interaction energy but is overcompensated by the favorable change in the solute-solvent entropy as observed in hydrophobic collapse. However, the negatively charged polymer collapses with a favorable change in solute-solvent interaction energy due to reduction in the dehydration energy penalty as the crowders partition to the polymer interface and shield the charged beads. The collapse of a charge-neutral polymer is opposed by the solute-solvent interaction energy but is overcompensated by the solute-solvent entropy change. However, for the strongly interacting crowders, the overall energetic penalty decreases since the crowders interact with polymer beads via cohesive bridging attractions to induce polymer collapse. These bridging attractions are found to be sensitive to the binding sites of the polymer, since they are absent in the negatively charged or uncharged polymers. These interesting differences in thermodynamic driving forces highlight the crucial role of the chemical nature of the macromolecule as well as of the crowder in determining the conformational equilibria in a crowded milieu. The results emphasize that the chemical interactions of the crowders should be explicitly considered to account for the crowding effects. The findings have implications in understanding the crowding effects on the protein free energy landscapes.

Molecular Basis for Actin Polymerization Kinetics Modulated by Solution Crowding

Actin polymerization drives cell movement and provides cells with structural integrity. Intracellular environments contain high concentrations of solutes, including organic compounds, macromolecules, and proteins. Macromolecular crowding has been shown to affect actin filament stability and bulk polymerization kinetics. However, the molecular mechanisms behind how crowding influences individual actin filament assembly are not well understood. In this study, we investigated how crowding modulates filament assembly kinetics using total internal reflection fluorescence (TIRF) microscopy imaging and pyrene fluorescence assays. The elongation rates of individual actin filaments analyzed from TIRF imaging depended on the type of crowding agent (polyethylene glycol, bovine serum albumin, and sucrose) as well as their concentrations. Further, we utilized all-atom molecular dynamics (MD) simulations to evaluate the effects of crowding molecules on the diffusion of actin monomers during filament assembly. Taken together, our data suggest that solution crowding can regulate actin assembly kinetics at the molecular level.

Brownian Motion in Optical Tweezers, a Comparison between MD Simulations and Experimental Data in the Ballistic Regime

The four most popular water models in molecular dynamics were studied in large-scale simulations of Brownian motion of colloidal particles in optical tweezers and then compared with experimental measurements in the same time scale. We present the most direct comparison of colloidal polystyrene particle diffusion in molecular dynamics simulations and experimental data on the same time scales in the ballistic regime. The four most popular water models, all of which take into account electrostatic interactions, are tested and compared based on yielded results and resources required. Three different conditions were simulated: a freely moving particle and one in a potential force field with two different strengths based on 1 pN/nm and 10 pN/nm. In all cases, the diameter of the colloidal particle was 50 nm. The acquired data were compared with experimental measurements performed using optical tweezers with position capture rates as high as 125 MHz. The experiments were performed in pure water on polystyrene particles with a 1 μm diameter in special microchannel cells.

Monte Carlo Simulation of Static and Dynamic Properties of Linear Polymer in a Crowded Environment

In this paper, we investigate the static and dynamic properties of linear polymer in the presence of obstacles. A Monte Carlo (MC) simulation method in two dimensions with a bond fluctuation model (BFM) was used to achieve this goal. To overcome the entropic barrier, we put the middle monomer of the polymer in the middle of the pore, which is placed between ordered and disordered obstacles. We probed the static properties of the polymer by calculating the mean square of the radius of gyration and the mean square end-to-end distance of the polymer, and we found that the scaling exponents of both the mean square end-to-end distance $〈R^2〉$ and the mean square radius of gyration $〈R_g^2〉$ as a function of the polymer length $N$ vary with the area fraction of crowding agents, $\varphi$ . The dynamic properties have also been studied by exploring the translocation of the polymer. Our current research shows that the escape time $\tau$ increases as $\varphi$ increases. Moreover, in the power-law relation of escape time $\tau$ as a function of polymer length $N$ , the scaling exponent ( $\alpha$ ) changes with $\varphi$ . Furthermore, the study has shown that the translocation of the polymer favors the disordered barriers.

Topological digestion drives time-varying rheology of entangled DNA fluids

Understanding and controlling the rheology of polymeric complex fluids that are pushed out-of-equilibrium is a fundamental problem in both industry and biology. For example, to package, repair, and replicate DNA, cells use enzymes to constantly manipulate DNA topology, length, and structure. Inspired by this feat, here we engineer and study DNA-based complex fluids that undergo enzymatically-driven topological and architectural alterations via restriction endonuclease (RE) reactions. We show that these systems display time-dependent rheological properties that depend on the concentrations and properties of the comprising DNA and REs. Through time-resolved microrheology experiments and Brownian Dynamics simulations, we show that conversion of supercoiled to linear DNA topology leads to a monotonic increase in viscosity. On the other hand, the viscosity of entangled linear DNA undergoing fragmentation displays a universal decrease that we rationalise using living polymer theory. Finally, to showcase the tunability of these behaviours, we design a DNA fluid that exhibits a time-dependent increase, followed by a temporally-gated decrease, of its viscosity. Our results present a class of polymeric fluids that leverage naturally occurring enzymes to drive diverse time-varying rheology by performing architectural alterations to the constituents.

Understanding and controlling the rheology of polymeric complex fluids is of fundamental importance in both industry and biology. Here, Michieletto et al. show how to achieve time-dependent rheology of DNA solutions via enzymatically-driven architectural alterations by restriction endonucleases.

Impact of Inert Crowders on Host-Guest Recognition Process

Biological environments typically contain high concentrations (300-400 mg/mL) of different macromolecules at volume fractions as large as 30%-40%. Biomolecular recognition processes, a ubiquitous biological phenomena, occurring in such crowded heterogeneous media would differ significantly compared to the dilute buffer solutions. Here we quantify the potential impact of inert crowders on prototypical host-guest recognition process by explicit-solvent molecular dynamics (MD) simulations in atomic resolution. We demonstrate that the crowders, when smaller in size, would facilitate the binding process of the guest molecule by decreasing the free energy barrier for binding via excluded volume effect and desolvation of the host receptor. However, the extent of crowder-induced stabilization of a host-guest complex is found to be significantly higher when the guest molecule is sterically constricted to approach the host along a centrosymmetric direction, compared to its unrestricted, freely diffusive movement. A kinetic analysis of the recognition process reveals that the origin of a relatively stronger crowder impact during constricted movement of guest molecule lies in the relatively enhanced residence time of the guest inside the host by crowders. Together, our results suggest that the extent of impact of crowding on recognition processes would be contingent upon the presence or absence of constriction on ligand movement.

Atomistic Simulation of Lysozyme in Solutions Crowded by Tetraethylene Glycol: Force Field Dependence

The behavior of biomolecules in crowded environments remains largely unknown due to the accuracy of simulation models and the limited experimental data for comparison. Here we chose a small crowder of tetraethylene glycol (PEG-4) to investigate the self-crowding of PEG-4 solutions and molecular crowding effects on the structure and diffusion of lysozyme at varied concentrations from dilute water to pure PEG-4 liquid. Two Amber-like force fields of Amber14SB and a99SB-disp were examined with TIP3P (fast diffusivity and low viscosity) and a99SB-disp (slow diffusivity and high viscosity) water models, respectively. Compared to the Amber14SB protein simulations, the a99SB-disp model yields more coordinated water and less PEG-4 molecules, less intramolecular hydrogen bonds (HBs), more protein-water HBs, and less protein-PEG HBs as well as stronger interactions and more hydrophilic and less hydrophobic contacts with solvent molecules. The a99SB-disp model offers comparable protein-solvent interactions in concentrated PEG-4 solutions to that in pure water. The PEG-4 crowding leads to a slow-down in the diffusivity of water, PEG-4, and protein, and the decline in the diffusion from atomistic simulations is close to or faster than the hard sphere model that neglects attractive interactions. Despite these differences, the overall structure of lysozyme appears to be maintained well at different PEG-4 concentrations for both force fields, except a slightly large deviation at 370 K at low concentrations with the a99SB-disp model. This is mainly attributed to the strong intramolecular interactions of the protein in the Amber14SB force field and to the large viscosity of the a99SB-disp water model. The results indicate that the protein force fields and the viscosity of crowder solutions affect the simulation of biomolecules under crowding conditions.

Hydrodynamics of immiscible binary fluids with viscosity contrast: a multiparticle collision dynamics approach

We present a multiparticle collision dynamics (MPC) implementation of layered immiscible fluids A and B of different shear viscosities separated by planar interfaces. The simulated flow profile for imposed steady shear motion and the time-dependent shear stress functions are in excellent agreement with our continuum hydrodynamics results for the composite fluid. The wave-vector dependent transverse velocity auto-correlation functions (TVAF) in the bulk-fluid regions of the layers decay exponentially, and agree with those of single-phase isotropic MPC fluids. In addition, we determine the hydrodynamic mobilities of an embedded colloidal sphere moving steadily parallel or transverse to a fluid-fluid interface, as functions of the distance from the interface. The obtained mobilities are in good agreement with hydrodynamic force multipoles calculations, for a no-slip sphere moving under creeping flow conditions near a clean, ideally flat interface. The proposed MPC fluid-layer model can be straightforwardly implemented, and it is computationally very efficient. Yet, owing to the spatial discretization inherent to the MPC method, the model can not reproduce all hydrodynamic features of an ideally flat interface between immiscible fluids.

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.

Oleate-modified hyaluronan: Controlling the number and distribution of side chains by varying the reaction conditions

In this work, low molecular weight hyaluronan was chemically modified by oleoyl moieties utilising mixed anhydrides methodology. The activation of oleic acid with benzoyl chloride in organic solvents miscible with water was followed by NMR spectroscopy. The product selectivity correlates with the solvent's Hildebrand solubility parameter. Furthermore, the effect of the solvent for the mixed anhydride formation was elucidated by density functional theory (DFT) and showed that the reactions are faster in acetonitrile or alcohols than in hexane. Furthermore, the solvent demonstrated to control the substituent distribution pattern along HA chain during esterification. An even distribution of substituents was observed in reactions performed in water mixed with ethers. The substituent distribution pattern clearly influenced the aggregation behaviour of amphiphilic HA, controlling the stability of the delivery system, while increasing the encapsulation capacity.

New Insights into the Functions of Nucleic Acids Controlled by Cellular Microenvironments

The right-handed double-helical B-form structure (B-form duplex) has been widely recognized as the canonical structure of nucleic acids since it was first proposed by James Watson and Francis Crick in 1953. This B-form duplex model has a monochronic and static structure and codes genetic information within a sequence. Interestingly, DNA and RNA can form various non-canonical structures, such as hairpin loops, left-handed helices, triplexes, tetraplexes of G-quadruplex and i-motif, and branched junctions, in addition to the canonical structure. The formation of non-canonical structures depends not only on sequence but also on the surrounding environment. Importantly, these non-canonical structures may exhibit a wide variety of biological roles by changing their structures and stabilities in response to the surrounding environments, which undergo vast changes at specific locations and at specific times in cells. Here, we review recent progress regarding the interesting behaviors and functions of nucleic acids controlled by molecularly crowded cellular conditions. New insights gained from recent studies suggest that nucleic acids not only code genetic information in sequences but also have unknown functions regarding their structures and stabilities through drastic structural changes in cellular environments.

Growth of Colloidal Nanocrystals by Liquid-Like Coalescence*

We here describe, model, and predict the growth kinetics of amine-capped PbS colloidal nanoparticles in the absence of supersaturation. The particles grow by coalescence rather than by Ostwald ripening. A comparison of different models indicates that the effective activation energy of coalescence (67.65 kJ mol^-1 ) is associated with two terms: a term proportional to the contact area between the ligand shells of two colliding particles, and a constant term. Our Brownian dynamics simulations show (i) how the remarkably low activation energy (or large rate constants) are most likely due to the large difference in size between the particles and their mean free path of diffusion, and (ii) how the low polydispersity is the likely result of the suppression of collision rates between rare populations due to crowding. The model successfully predicts the growth kinetics of nanoparticles, therefore enabling the precise control of the average particle size without the need of supersaturation.

Concentration sensing in crowded environments

Signal transduction within crowded cellular compartments is essential for the physiological function of cells. Although the accuracy with which receptors can probe the concentration of ligands has been thoroughly investigated in dilute systems, the effect of macromolecular crowding on the inference of concentration remains unclear. In this work, we develop an algorithm to simulate reversible reactions between reacting Brownian particles. Our algorithm facilitates the calculation of reaction rates and correlation times for ligand-receptor systems in the presence of macromolecular crowding. Using this method, we show that it is possible for crowding to increase the accuracy of estimated ligand concentration based on receptor occupancy. In particular, we find that crowding can enhance the effective association rates between small ligands and receptors to a degree sufficient to overcome the increased chance of rebinding due to caging by crowding molecules. For larger ligands, crowding decreases the accuracy of the receptor's estimate primarily by decreasing the microscopic association and dissociation rates.

Quantifying the roles of space and stochasticity in computer simulations for cell biology and cellular biochemistry

Most of the fascinating phenomena studied in cell biology emerge from interactions among highly organized multimolecular structures embedded into complex and frequently dynamic cellular morphologies. For the exploration of such systems, computer simulation has proved to be an invaluable tool, and many researchers in this field have developed sophisticated computational models for application to specific cell biological questions. However, it is often difficult to reconcile conflicting computational results that use different approaches to describe the same phenomenon. To address this issue systematically, we have defined a series of computational test cases ranging from very simple to moderately complex, varying key features of dimensionality, reaction type, reaction speed, crowding, and cell size. We then quantified how explicit spatial and/or stochastic implementations alter outcomes, even when all methods use the same reaction network, rates, and concentrations. For simple cases, we generally find minor differences in solutions of the same problem. However, we observe increasing discordance as the effects of localization, dimensionality reduction, and irreversible enzymatic reactions are combined. We discuss the strengths and limitations of commonly used computational approaches for exploring cell biological questions and provide a framework for decision making by researchers developing new models. As computational power and speed continue to increase at a remarkable rate, the dream of a fully comprehensive computational model of a living cell may be drawing closer to reality, but our analysis demonstrates that it will be crucial to evaluate the accuracy of such models critically and systematically.

Interplay of reactive interference and crowding effects in the diffusion-influenced reaction kinetics

We investigate the interplay of reactive interference and crowding effects in the irreversible diffusion-influenced bimolecular reactions of the type A+B→P+B by using the Brownian dynamics simulation method. It is known that the presence of nonreactive crowding agents retards the reaction rate when the volume fraction of the crowding agents is large enough. On the other hand, a high concentration of B is known to increase the reaction rate more than expected from the mass action law, although the B's may also act as crowders. Therefore, it would be interesting to see which effect dominates when the number density of B as well as the number density of the crowders increases. We will present an approximate theory that provides a reasonable account for the Brownian dynamics simulation results.

Screening and characterization of an α-L-fucosidase from Bacteroides fragilis NCTC9343 for synthesis of fucosyl-N-acetylglucosamine disaccharides

Fucosyl-N-acetylglucosamine disaccharides are present in many biologically important oligosaccharides, such as human milk oligosaccharides, Lewis carbohydrate antigens, and glycans on cell-surface glycoconjugate receptors, and thus have vast potential for infant formulas, prebiotics, and pharmaceutical applications. In this work, in order to screen biocatalysts for enzymatic synthesis of fucosyl-N-acetylglucosamine disaccharides, we performed sequence analysis of 12 putative and one known α-L-fucosidases of Bacteroides fragilis NCTC9343 and constructed a phylogenetic tree of the nine GH29 α-L-fucosidases. After that, five GH29A α-L-fucosidases were cloned, and four of them were successfully heterogeneous expressed and screened for transglycosylation activity, and a GH29A α-L-fucosidase (BF3242) that synthesized a mix of Fuc-α-1,3/1,6-GlcNAc disaccharides using pNPαFuc as donor and GlcNAc as acceptor was characterized. The effects of initial substrate concentration, pH, temperature, and reaction time on its transglycosylation activity were studied in detail. Under the optimum conditions of 0.05 U/mL enzyme, 20 mM pNPαFuc, and 500 mM GlcNAc in sodium buffer (pH 7.5) at 37 °C for 45 min, BF3242 efficiently synthesized Fuc-α-1,3/1,6-GlcNAc at a maximum yield of 79.0% with the ratio of 0.48 for 1,3/1,6. The molecular dynamics simulation analysis revealed that Loop-4 (His220-Ser245) in the putative 3D model of BF3242 displayed significant changes throughout the thermal simulations, might being responsible for the changes in the ratio of two regioisomeric products at different temperatures. This work provided not only a potential synthetic tool for enzymatic synthesis of fucosyl-N-acetylglucosamine disaccharides but also a possibility for the formation of regioisomeric products in glycosidase-catalyzed transglycosylation. KEY POINTS: • Sequence analysis of α-L-fucosidases of Bacteroides fragilis NCTC9343 • Obtainment of an α-L-fucosidase with high transglycosylation activity • Explanation why temperature affected the ratio of two regioisomeric products.

Cascade Kinetics in an Enzyme-Loaded Aqueous Two-Phase System

Macromolecular crowding plays a critical role in the kinetics of enzymatic reactions. Dynamic compartmentalization of biological components in living cells due to liquid-liquid phase separation represents an important cell regulatory mechanism that can increase enzyme concentration locally and influence the diffusion of substrates. In the present study, we probed partitioning of two enzymes (horseradish-peroxidase and urate-oxidase) in a poly(ethylene glycol)-dextran aqueous two-phase system (ATPS) as a function of salt concentration and ion position in the Hofmeister series. Moreover, we investigated enzymatic cascade reactions and their kinetics within the ATPS, which revealed a strong influence of the ion hydration stemming from the background electrolyte on the partitioning coefficients of proteins following the Hofmeister series. As a result, we were able to realize cross-partitioning of two enzymes because of different protein net charges at a chosen pH. Our study reveals a strong dependency of the enzyme activity on the substrate type and crowding agent interaction on the final kinetics of enzymatic reactions in the ATPS and therefore provides substantial implications en route toward dynamic regulation of reactivity in synthetic protocells.

Rotational and Translational Diffusion of Proteins as a Function of Concentration

Atomistic simulations of three different proteins at different concentrations are performed to obtain insight into protein mobility as a function of protein concentration. We report on simulations of proteins from diluted to the physiological water concentration (about 70% of the mass). First, the viscosity was computed and found to increase by a factor of 7-9 going from pure water to the highest protein concentration, in excellent agreement with in vivo nuclear magnetic resonance results. At a physiological concentration of proteins, the translational diffusion is found to be slowed down to about 30% of the in vitro values. The slow-down of diffusion found here using atomistic models is slightly more than that of a hard sphere model that neglects the electrostatic interactions. Interestingly, rotational diffusion of proteins is slowed down somewhat more (by about 80-95% compared to in vitro values) than translational diffusion, in line with experimental findings and consistent with the increased viscosity. The finding that rotation is retarded more than translation is attributed to solvent-separated clustering. No direct interactions between the proteins are found, and the clustering can likely be attributed to dispersion interactions that are stronger between proteins than between protein and water. Based on these simulations, we can also conclude that the internal dynamics of the proteins in our study are affected only marginally under crowding conditions, and the proteins become somewhat more stable at higher concentrations. Simulations were performed using a force field that was tuned for dealing with crowding conditions by strengthening the protein-water interactions. This force field seems to lead to a reproducible partial unfolding of an α-helix in one of the proteins, an effect that was not observed in the unmodified force field.

Rapid and room temperature detection of single nucleotide variation with enhanced discrimination by crowding assisted allele specific extension

In this paper we report the kinetics based detection of single nucleotide variation (SNV) at room temperature by allele specific extension with different concentrations and types of crowding agents. In general, the crowding conditions enhanced the specificity in the detection of SNV.

Macromolecular crowding effects on transcription and translation are regulated by free magnesium ion

Cell-free metabolic engineering is an emerging and promising alternative platform for the production of fuels and chemicals. In recent years, macromolecular crowding effect, which is an important function in living cells but ignored in cell-free systems, has been transferred to cell-free protein synthesis (CFPS). However, inhibitory effects of crowding agents on CFPS were frequently observed, and the mechanism is unclear. In this study, free Mg²⁺ was found to be a key factor that can regulate the macromolecular crowding effect on in vitro transcription, in vitro translation, and coupled transcript/translation. Addition of crowding agents (20% of Ficoll-70 or Ficoll-400) enhanced in vitro transcription at an index of free Mg²⁺ concentration (IFMC) below 2 mM but inhibited the transcription when the IFMC was higher than 2 mM. Similarly, Ficoll-400 enhanced in vitro translation and coupled transcription/translation at a lower IFMC (0.1-2 mM) and inhibited the reactions at higher IFMC (>2 mM). Based on the results, CFPS systems could be further optimized by adjusting the content of crowding agents and the IFMC. Besides, the results also indicate that macromolecular crowding effect is important for maintaining the efficiency of in vivo transcription and translation which occur at a low intracellular IFMC (<1 mM).

Measuring the effects of macromolecular crowding on antibody function with biolayer interferometry

Biotherapeutic proteins are commonly dosed at high concentrations into the blood, which is an inherently complex, crowded solution with substantial protein content. The effects of macromolecular crowding may lead to an appreciable level of non-specific hetero-association in this physiological environment. Therefore, developing a method to characterize the diverse consequences of non-specific interactions between proteins under such non-ideal, crowded conditions, which deviate substantially from those commonly employed for in vitro characterization, is vital to achieving a more complete picture of antibody function in a biological context. In this study, we investigated non-specific interactions between human serum albumin (HSA) and two monoclonal antibodies (mAbs) by static light scattering and determined these interactions are both ionic strength-dependent and mAb-dependent. Using biolayer interferometry (BLI), we assessed the effect of HSA on antigen binding by mAbs, demonstrating that these non-specific interactions have a functional impact on mAb:antigen interactions, particularly at low ionic strength. While this effect is mitigated at physiological ionic strength, our in vitro data support the notion that HSA in the blood may lead to non-specific interactions with mAbs in vivo, with a potential impact on their interactions with antigen. Furthermore, the BLI method offers a high-throughput advantage compared to orthogonal techniques such as analytical ultracentrifugation and is amenable to a greater variety of solution conditions compared to nuclear magnetic resonance spectroscopy. Our study demonstrates that BLI is a viable technology for examining the impact of non-specific interactions on specific biologically relevant interactions, providing a direct method to assess binding events in crowded conditions.

Single-molecule visualization of the effects of ionic strength and crowding on structure-mediated interactions in supercoiled DNA molecules

DNA unwinding is an important cellular process involved in DNA replication, transcription and repair. In cells, molecular crowding caused by the presence of organelles, proteins, and other molecules affects numerous internal cellular structures. Here, we visualize plasmid DNA unwinding and binding dynamics to an oligonucleotide probe as functions of ionic strength, crowding agent concentration, and crowding agent species using single-molecule CLiC microscopy. We demonstrate increased probe–plasmid interaction over time with increasing concentration of 8 kDa polyethylene glycol (PEG), a crowding agent. We show decreased probe–plasmid interactions as ionic strength is increased without crowding. However, when crowding is introduced via 10% 8 kDa PEG, interactions between plasmids and oligos are enhanced. This is beyond what is expected for normal in vitro conditions, and may be a critically important, but as of yet unknown, factor in DNA’s proper biological function in vivo. Our results show that crowding has a strong effect on the initial concentration of unwound plasmids. In the dilute conditions used in these experiments, crowding does not impact probe–plasmid interactions once the site is unwound.

Anomalous transport in the soft-sphere Lorentz model

The sensitivity of anomalous transport in crowded media to the form of the inter-particle interactions is investigated through computer simulations. We extend the highly simplified Lorentz model towards realistic natural systems by modeling the interactions between the tracer and the obstacles with a smooth potential. We find that the anomalous transport at the critical point happens to be governed by the same universal exponent as for hard exclusion interactions, although the mechanism of how narrow channels are probed is rather different. The scaling behavior of simulations close to the critical point confirm this exponent. Our result indicates that the simple Lorentz model may be applicable to describing the fundamental properties of long-range transport in real crowded environments.

Polydispersity and negative charge are key modulators of extracellular matrix deposition under macromolecular crowding conditions

Macromolecular crowding is a biophysical phenomenon that stems from the volume excluded by macromolecules, as they undergo steric repulsion and electrostatic interactions. The excluded volume depends on the shape, size, charge and polydispersity of the molecules. Although theoretical/computational models have been used to assess the influence of macromolecular crowding in biological media, real-time experiments are scarce. Herein, we evaluated the influence of hydrodynamic radius, charge and polydispersity of (a) various concentrations of different crowders (carrageenan, Ficoll™ and dextran sulphate); (b) various molecular weights of different crowders (70, 400 and 100 kDa of Ficoll™ and 10, 100 and 500 kDa of dextran sulphate) and (c) various cocktails of the same crowders (cocktails of various concentrations of different molecular weights Ficoll™ and dextran sulphate) on extracellular matrix deposition in human dermal fibroblast culture. The use of crowding cocktails with different molecular weight/concentrations of Ficoll™ or dextran sulphate molecules led to increased polydispersity and enhanced collagen type I deposition in comparison to their mono-domain counterparts. Carrageenan, however, induced the highest deposition of collagen type I due to its negative charge and inherent polydispersity. Our data contribute to a better understanding of the influence of the biophysical properties of the crowders on extracellular matrix deposition in vitro. STATEMENT OF SIGNIFICANCE: Macromolecular crowding is a biophysical phenomenon that accelerates and enhances extracellular matrix deposition in cell culture systems. Herein, we demonstrate that negatively charged and polydispersed macromolecules or cocktails of macromolecules, as opposed to neutral and monodomain macromolecules, induce highest extracellular matrix deposition in human dermal fibroblast cultures.

Liquid-liquid phase separation in artificial cells

Liquid-liquid phase separation (LLPS) in biology is a recently appreciated means of intracellular compartmentalization. Because the mechanisms driving phase separations are grounded in physical interactions, they can be recreated within less complex systems consisting of only a few simple components, to serve as artificial microcompartments. Within these simple systems, the effect of compartmentalization and microenvironments upon biological reactions and processes can be studied. This review will explore several approaches to incorporating LLPS as artificial cytoplasms and in artificial cells, including both segregative and associative phase separation.

Effect of the Organization of Rhodopsin on the Association between Transducin and a Photoactivated Receptor

After photoactivation, rhodopsin (R), a G-protein-coupled receptor, rapidly activates multiple transducin G-proteins (G) in an initial amplification step of phototransduction. G-protein activation requires diffusion-mediated association with an active rhodopsin (R*) at the rod disk membrane. Different organizations of R within the membrane have been revealded by several microscopy studies, including static and freely diffusing situations. However, it is unclear how such different scenarios influence the activation rate of G proteins. Through Monte Carlo simulations, we study the association reaction between a photoactivated rhodopsin and transducin under different reported receptor organizations including (a) R monomers diffusing freely, (b) R forming static dispersed crystalline domains made of rows of dimers, and (c) R arranged in static tracks formed by two adjacent rows of dimers. A key parameter in our simulations is the probability of binding following a collision ( p). For high p, the association rate between R* and G is higher in the freely diffusive system than in the static organizations, but for low collision efficiencies, the static organizations can result in faster association rates than the mobile system. We also observe that with low p, increasing the concentration of R increases the association rate significantly in the dispersed crystals configuration and just slightly in the free diffusive system. In summary, the lateral organization of rhodopsin influences the association rate between R* and G in a manner strongly dependent on the collision efficiency.

Discrete and Continuum Models for the Salt in Crowded Environments of Suspended Charged Particles

Electrostatic forces greatly affect the overall dynamics and diffusional activities of suspended charged particles in crowded environments. Accordingly, the concentration of counter- or co-ions in a fluid-''the salt"-determines the range, strength, and order of electrostatic interactions between particles. This environment fosters engineering routes for controlling directed assembly of particles at both the micro- and nanoscale. Here, we analyzed two computational modeling schemes that considered salt within suspensions of charged particles, or polyelectrolytes: discrete and continuum. Electrostatic interactions were included through a Green's function formalism, where the confined fundamental solution for Poisson's equation is resolved by the general geometry Ewald-like method. For the discrete model, the salt was considered as regularized point-charges with a specific valence and size, while concentration fields were defined for each ionic species for the continuum model. These considerations were evolved using Brownian dynamics of the suspended charged particles and the discrete salt ions, while a convection-diffusion transport equation, including the Nernst-Planck diffusion mechanism, accounted for the dynamics of the concentration fields. The salt/particle models were considered as suspensions under slit-confinement conditions for creating crowded "macro-ions", where density distributions and radial distribution functions were used to compare and differentiate computational models. Importantly, our analysis shows that disparate length scales or increased system size presented by the salt and suspended particles are best dealt with using concentration fields to model the ions. These findings were then validated by novel simulations of a semipermeable polyelectrolyte membrane, at the mesoscale, from which ionic channels emerged and enable ion conduction.

Pressure and cosolvent modulation of the catalytic activity of amyloid fibrils

We report on the effects of pressure and cosolvents on the catalytic activity of a designed amyloid fibril by applying a high-pressure stopped-flow methodology with rapid spectroscopic detection. FTIR spectroscopic data revealed a remarkable pressure and temperature stability of the fibrillar catalyst. The activity is further enhanced by osmolytes and macromolecular crowding.

Product interactions and feedback in diffusion-controlled reactions

Steric or attractive interactions among reactants or between reactants and inert crowders can substantially influence the total rate of a diffusion-influenced reaction in the liquid phase. However, the role of the product species, which has typically different physical properties than the reactant species, has been disregarded so far. Here we study the effects of reactant-product and product-product interactions as well as asymmetric diffusion properties on the rate of diffusion-controlled reactions in the classical Smoluchowski-setup for chemical transformations at a perfect catalytic sphere. For this, we solve the diffusion equation with appropriate boundary conditions coupled by a mean-field approach on the second virial level to account for the particle interactions. We find that all particle spatial distributions and the total rate can change significantly, depending on the diffusion and interaction properties of the accumulated products. Complex competing and self-regulating (homeostatic) or self-amplifying effects are observed for the system, leading to both decrease and increase in the rates, as the presence of interacting products feeds back to the reactant flux and thus the rate with which the products are generated.

Macromolecular Crowding Induces Spatial Correlations That Control Gene Expression Bursting Patterns

Recent superresolution microscopy studies in E. coli demonstrate that the cytoplasm has highly variable local concentrations where macromolecular crowding plays a central role in establishing membrane-less compartmentalization. This spatial inhomogeneity significantly influences molecular transport and association processes central to gene expression. Yet, little is known about how macromolecular crowding influences gene expression bursting-the episodic process where mRNA and proteins are produced in bursts. Here, we simultaneously measured mRNA and protein reporters in cell-free systems, showing that macromolecular crowding decoupled the well-known relationship between fluctuations in the protein population (noise) and mRNA population statistics. Crowded environments led to a 10-fold increase in protein noise even though there were only modest changes in the mRNA population and fluctuations. Instead, cell-like macromolecular crowding created an inhomogeneous spatial distribution of mRNA ("spatial noise") that led to large variability in the protein production burst size. As a result, the mRNA spatial noise created large temporal fluctuations in the protein population. These results highlight the interplay between macromolecular crowding, spatial inhomogeneities, and the resulting dynamics of gene expression, and provide insights into using these organizational principles in both cell-based and cell-free synthetic biology.

Making microenvironments: A look into incorporating macromolecular crowding into in vitro experiments, to generate biomimetic microenvironments which are capable of directing cell function for tissue engineering applications

Biomimetic microenvironments are key components to successful cell culture and tissue engineering in vitro. One of the most accurate biomimetic microenvironments is that made by the cells themselves. Cell-made microenvironments are most similar to the in vivo state as they are cell-specific and produced by the actual cells which reside in that specific microenvironment. However, cell-made microenvironments have been challenging to re-create in vitro due to the lack of extracellular matrix composition, volume and complexity which are required. By applying macromolecular crowding to current cell culture protocols, cell-made microenvironments, or cell-derived matrices, can be generated at significant rates in vitro. In this review, we will examine the causes and effects of macromolecular crowding and how it has been applied in several in vitro systems including tissue engineering.

The use of a genetically encoded molecular crowding sensor in various biological phenomena

We evaluated usability of a previously developed genetically encoded molecular crowding sensor in various biological phenomena. Molecular crowding refers to intracellular regions that are occupied more by proteins and nucleotides than by water molecules and is thought to have a strong effect on protein function. To evaluate intracellular molecular crowding, usually the diffusion coefficient of a probe is used because it is related to mobility of the surrounding molecular crowding agents. Recently, genetically encoded molecular crowding sensors based on Förster resonance energy transfer were reported. In the present study, to evaluate the usability of a genetically encoded molecular crowding sensor, molecular crowding was monitored during several biological events. Changes in molecular crowding during stem cell differentiation, cell division, and focal adhesion development and difference in molecular crowding in filopodia locations were examined. The results show usefulness of the genetically encoded molecular crowding sensor for understanding the biological phenomena relating to molecular crowding.

Facilitated diffusion in the presence of obstacles on the DNA

Biological functions of DNA depend on the sequence-specific binding of DNA-binding proteins to their corresponding binding sites. Binding of these proteins to their binding sites occurs through a facilitated diffusion process that combines three-dimensional diffusion in the cytoplasm with one-dimensional diffusion (sliding) along the DNA. In this work, we use a lattice model of facilitated diffusion to study how the dynamics of binding of a protein to a specific site (e.g., binding of an RNA polymerase to a promoter or of a transcription factor to its operator site) is affected by the presence of other proteins bound to the DNA, which act as 'obstacles' in the sliding process. Different types of these obstacles with different dynamics are implemented. While all types impair facilitated diffusion, the extent of the hindrance depends on the type of obstacle. As a consequence of hindrance by obstacles, more excursions into the cytoplasm are required for optimal target binding compared to the case without obstacles.

Does metabolite channeling accelerate enzyme-catalyzed cascade reactions?

Metabolite or substrate channeling is a direct transfer of metabolites from one enzyme to the next enzyme in a cascade. Among many potential advantages of substrate channeling, acceleration of the total reaction rate is considered as one of the most important and self-evident. However, using a simple model, supported by stochastic simulations, we show that it is not always the case; particularly at long times (i.e. in steady state) and high substrate concentrations, a channeled reaction cannot be faster, and can even be slower, than the original non-channeled cascade reaction. In addition we show that increasing the degree of channeling may lead to an increase of the metabolite pool size. We substantiate that the main advantage of channeling likely lies in protecting metabolites from degradation or competing side reactions.

Hyperbranched polyglycerols: recent advances in synthesis, biocompatibility and biomedical applications

Hyperbranched polyglycerol is one of the most widely studied biocompatible dendritic polymer and showed promising applications. Here, we summarized the recent advancements in the field.

Ion specificities of artificial macromolecules

Artificial macromolecules are well-defined synthetic polymers, with a relatively simple structure as compared to naturally occurring macromolecules. This review focuses on the ion specificities of artifical macromolecules. Ion specificities are influenced by solvent-mediated indirect ion-macromolecule interactions and also by direct ion-macromolecule interactions. In aqueous solutions, the role of water-mediated indirect ion-macromolecule interactions will be discussed. The addition of organic solvents to aqueous solutions significantly changes the ion specificities due to the formation of water-organic solvent complexes. For direct ion-macromolecule interactions, we will discuss specific ion-pairing interactions for charged macromolecules and specific ion-neutral site interactions for uncharged macromolecules. When the medium conditions change from dilute solutions to crowded environments, the ion specificities can be modified by either the volume exclusion effect, the variation of dielectric constant, or the interactions between ions, macromolecules, and crowding agents.

Population dynamics, information transfer, and spatial organization in a chemical reaction network under spatial confinement and crowding conditions

We investigate, via Brownian dynamics simulations, the reaction dynamics of a generic, nonlinear chemical network under spatial confinement and crowding conditions. In detail, the Willamowski-Rossler chemical reaction system has been "extended" and considered as a prototype reaction-diffusion system. Our results are potentially relevant to a number of open problems in biophysics and biochemistry, such as the synthesis of primitive cellular units (protocells) and the definition of their role in the chemical origin of life and the characterization of vesicle-mediated drug delivery processes. More generally, the computational approach presented in this work makes the case for the use of spatial stochastic simulation methods for the study of biochemical networks in vivo where the "well-mixed" approximation is invalid and both thermal and intrinsic fluctuations linked to the possible presence of molecular species in low number copies cannot be averaged out.

Dry/Wet Cycling and the Thermodynamics and Kinetics of Prebiotic Polymer Synthesis

The endoergic nature of protein and nucleic acid assembly in aqueous media presents two questions that are fundamental to the understanding of life’s origins: (i) how did the polymers arise in an aqueous prebiotic world; and (ii) once formed in some manner, how were they sufficiently persistent to engage in further chemistry. We propose here a quantitative resolution of these issues that evolved from recent accounts in which RNA-like polymers were produced in evaporation/rehydration cycles. The equilibrium N_m + N_n ↔ N_m+n + H₂O is endoergic by about 3.3 kcal/mol for polynucleotide formation, and the system thus lies far to the left in the starting solutions. Kinetic simulations of the evaporation showed that simple Le Châtelier’s principle shifts were insufficient, but the introduction of oligomer-stabilizing factors of 5–10 kcal/mol both moved the process to the right and respectively boosted and retarded the elongation and hydrolysis rates. Molecular crowding and excluded volume effects in present-day cells yield stabilizing factors of that order, and we argue here that the crowded conditions in the evaporites generate similar effects. Oligomer formation is thus energetically preferred in those settings, but the process is thwarted in each evaporation step as diffusion becomes rate limiting. Rehydration dissipates disordered oligomer clusters in the evaporites, however, and subsequent dry/wet cycling accordingly “ratchets up” the system to an ultimate population of kinetically trappedthermodynamically preferred biopolymers.

Quantifying the Influence of the Crowded Cytoplasm on Small Molecule Diffusion

Cytosolic crowding can influence the thermodynamics and kinetics of in vivo chemical reactions. Most significantly, proteins and nucleic acid crowders reduce the accessible volume fraction, ϕ, available to a diffusing substrate, thereby reducing its effective diffusion rate, Deff, relative to its rate in bulk solution. However, Deff can be further hindered or even enhanced, when long-range crowder/diffuser interactions are significant. To probe these effects, we numerically estimated Deff values for small, charged molecules in representative, cytosolic protein lattices up to 0.1 × 0.1 × 0.1 μm(3) in volume via the homogenized Smoluchowski electro-diffusion equation. We further validated our predictions against Deff estimates from ϕ-dependent analytical relationships, such as the Maxwell-Garnett (MG) bound, as well as explicit solutions of the time-dependent electro-diffusion equation. We find that in typical, moderately crowded cell cytoplasm (ϕ ≈ 0.8), Deff is primarily determined by ϕ; in other words, diverse protein shapes and heterogeneous distributions only modestly impact Deff. However, electrostatic interactions between diffusers and crowders, particularly at low electrolyte ionic strengths, can substantially modulate Deff. These findings help delineate the extent that cytoplasmic crowders influence small molecule diffusion, which ultimately may shape the efficiency and timing of intracellular signaling pathways. More generally, the quantitative agreement between computationally expensive solutions of the time-dependent electro-diffusion equation and its comparatively cheaper homogenized form suggest that the latter is a broadly effective model for diffusion in wide-ranging, crowded biological media.