Thermodynamics of Macromolecular Association in Heterogeneous Crowding Environments: Theoretical and Simulation Studies with a Simplified Model

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.

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

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

GB1 Dimerization in Crowders: A Multiple Resolution Approach

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

Macromolecular crowding effects on electrostatic binding affinity: Fundamental insights from theoretical, idealized models

The crowded cellular environment can affect biomolecular binding energetics, with specific effects depending on the properties of the binding partners and the local environment. Often, crowding effects on binding are studied on particular complexes, which provide system-specific insights but may not provide comprehensive trends or a generalized framework to better understand how crowding affects energetics involved in molecular recognition. Here, we use theoretical, idealized molecules whose physical properties can be systematically varied along with samplings of crowder placements to understand how electrostatic binding energetics are altered through crowding and how these effects depend on the charge distribution, shape, and size of the binding partners or crowders. We focus on electrostatic binding energetics using a continuum electrostatic framework to understand effects due to depletion of a polar, aqueous solvent in a crowded environment. We find that crowding effects can depend predictably on a system's charge distribution, with coupling between the crowder size and the geometry of the partners' binding interface in determining crowder effects. We also explore the effect of crowder charge on binding interactions as a function of the monopoles of the system components. Finally, we find that modeling crowding via a lowered solvent dielectric constant cannot account for certain electrostatic crowding effects due to the finite size, shape, or placement of system components. This study, which comprehensively examines solvent depletion effects due to crowding, complements work focusing on other crowding aspects to help build a holistic understanding of environmental impacts on molecular recognition.

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

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

Spermidine strongly increases the fidelity of Escherichia coli CRISPR Cas1-Cas2 integrase

Site-selective CRISPR array expansion at the origin of bacterial adaptive immunity relies on recognition of sequence-dependent DNA structures by the conserved Cas1-Cas2 integrase. Off-target integration of a new spacer sequence outside canonical CRISPR arrays has been described in vitro However, this nonspecific integration activity is rare in vivo Here, we designed gel assays to monitor fluorescently labeled protospacer insertion in a supercoiled 3-kb plasmid harboring a minimal CRISPR locus derived from the Escherichia coli type I-E system. This assay enabled us to distinguish and quantify target and off-target insertion events catalyzed by E. coli Cas1-Cas2 integrase. We show that addition of the ubiquitous polyamine spermidine or of another polyamine, spermine, significantly alters the ratio between target and off-target insertions. Notably, addition of 2 mm spermidine quenched the off-target spacer insertion rate by a factor of 20-fold, and, in the presence of integration host factor, spermidine also increased insertion at the CRISPR locus 1.5-fold. The observation made in our in vitro system that spermidine strongly decreases nonspecific activity of Cas1-Cas2 integrase outside the leader-proximal region of a CRISPR array suggests that this polyamine plays a potential role in the fidelity of the spacer integration also in vivo.

Computational approaches to macromolecular interactions in the cell

Structural modeling of a cell is an evolving strategic direction in computational structural biology. It takes advantage of new powerful modeling techniques, deeper understanding of fundamental principles of molecular structure and assembly, and rapid growth of the amount of structural data generated by experimental techniques. Key modeling approaches to principal types of macromolecular assemblies in a cell already exist. The main challenge, along with the further development of these modeling approaches, is putting them together in a consistent, unified whole cell model. This opinion piece addresses the fundamental aspects of modeling macromolecular assemblies in a cell, and the state-of-the-art in modeling of the principal types of such assemblies.

A global map of the protein shape universe

Proteins are involved in almost all functions in a living cell, and functions of proteins are realized by their tertiary structures. Obtaining a global perspective of the variety and distribution of protein structures lays a foundation for our understanding of the building principle of protein structures. In light of the rapid accumulation of low-resolution structure data from electron tomography and cryo-electron microscopy, here we map and classify three-dimensional (3D) surface shapes of proteins into a similarity space. Surface shapes of proteins were represented with 3D Zernike descriptors, mathematical moment-based invariants, which have previously been demonstrated effective for biomolecular structure similarity search. In addition to single chains of proteins, we have also analyzed the shape space occupied by protein complexes. From the mapping, we have obtained various new insights into the relationship between shapes, main-chain folds, and complex formation. The unique view obtained from shape mapping opens up new ways to understand design principles, functions, and evolution of proteins.

Proteins are the major molecules involved in almost all cellular processes. In this work, we present a novel mapping of protein shapes that represents the variety and the similarities of 3D shapes of proteins and their assemblies. This mapping provides various novel insights into protein shapes including determinant factors of protein 3D shapes, which enhance our understanding of the design principles of protein shapes. The mapping will also be a valuable resource for artificial protein design as well as references for classifying medium- to low-resolution protein structure images of determined by cryo-electron microscopy and tomography.

Recent advances in carrier mediated nose-to-brain delivery of pharmaceutics

Central nervous system (CNS) disorders (e.g., multiple sclerosis, Alzheimer's disease, etc.) represent a growing public health issue, primarily due to the increased life expectancy and the aging population. The treatment of such disorders is notably elaborate and requires the delivery of therapeutics to the brain in appropriate amounts to elicit a pharmacological response. However, despite the major advances both in neuroscience and drug delivery research, the administration of drugs to the CNS still remains elusive. It is commonly accepted that effectiveness-related issues arise due to the inability of parenterally administered macromolecules to cross the Blood-Brain Barrier (BBB) in order to access the CNS, thus impeding their successful delivery to brain tissues. As a result, the direct Nose-to-Brain delivery has emerged as a powerful strategy to circumvent the BBB and deliver drugs to the brain. The present review article attempts to highlight the different experimental and computational approaches pursued so far to attain and enhance the direct delivery of therapeutic agents to the brain and shed some light on the underlying mechanisms involved in the pathogenesis and treatment of neurological disorders.

In Aqua Veritas: The Indispensable yet Mostly Ignored Role of Water in Phase Separation and Membrane-less Organelles

Despite the common practice of presenting structures of biological molecules on an empty background and the assumption that interactions between biological macromolecules take place within the inert solvent, water represents an active component of various biological processes. This Perspective addresses indispensable, yet mostly ignored, roles of water in biological liquid-liquid phase transitions and in the biogenesis of various proteinaceous membrane-less organelles. We point out that changes in the structure of water reflected in the changes in its abilities to donate and/or accept hydrogen bonds and participate in dipole-dipole and dipole-induced dipole interactions in the presence of various solutes (ranging from small molecules to synthetic polymers and biological macromolecules) might represent a driving force for the liquid-liquid phase separation, define partitioning of various solutes in formed phases, and define the exceptional ability of intrinsically disordered proteins to be engaged in the formation of proteinaceous membrane-less organelles.

Molecular simulations of cellular processes

It is, nowadays, possible to simulate biological processes in conditions that mimic the different cellular compartments. Several groups have performed these calculations using molecular models that vary in performance and accuracy. In many cases, the atomistic degrees of freedom have been eliminated, sacrificing both structural complexity and chemical specificity to be able to explore slow processes. In this review, we will discuss the insights gained from computer simulations on macromolecule diffusion, nuclear body formation, and processes involving the genetic material inside cell-mimicking spaces. We will also discuss the challenges to generate new models suitable for the simulations of biological processes on a cell scale and for cell-cycle-long times, including non-equilibrium events such as the co-translational folding, misfolding, and aggregation of proteins. A prominent role will be played by the wise choice of the structural simplifications and, simultaneously, of a relatively complex energetic description. These challenging tasks will rely on the integration of experimental and computational methods, achieved through the application of efficient algorithms. Graphical abstract.

Role of Macromolecular Crowding on Stability and Iron Release Kinetics of Serum Transferrin

The macromolecular crowding influences the structural stability and functional properties of transferrin (Tf). The equilibrium as well as kinetic studies of Tf at different concentrations of crowding agents (dextran 40, dextran 70, and ficoll 70) and at a fixed concentration of dextran 40 under different concentrations of NaCl at pH 7.4 and 5.6 (±1) revealed that (i) the crowder environment increases the diferric-Tf (Fe₂Tf) stability against iron loss and overall denaturation of the protein, (ii) both in the absence and presence of crowder, the presence of salt promotes the loss of iron and overall denaturation of Fe₂Tf which is due to ionic screening of electrostatic interactions, (iii) the crowder environment retards iron release from monoferric N-lobe of Tf (Fe_NTf) by increasing enthalpic barrier, (iv) the retardation of iron release by crowding is enthalpically dominated than the entropic one, (v) both in the absence and presence of crowder, the presence of salt accelerates the iron release from Fe_NTf due to ionic screening of electrostatic interactions and anion binding to KISAB sites, and (vi) the crowders environment is unable to diminish (a) the salt-induced destabilization of Fe₂Tf against the loss of iron and overall denaturation and (b) the anion effect and ionic screening of diffusive counterions responsible to promote iron release from Fe_NTf.

Crowding in Cellular Environments at an Atomistic Level from Computer Simulations

The effects of crowding in biological environments on biomolecular structure, dynamics, and function remain not well understood. Computer simulations of atomistic models of concentrated peptide and protein systems at different levels of complexity are beginning to provide new insights. Crowding, weak interactions with other macromolecules and metabolites, and altered solvent properties within cellular environments appear to remodel the energy landscape of peptides and proteins in significant ways including the possibility of native state destabilization. Crowding is also seen to affect dynamic properties, both conformational dynamics and diffusional properties of macromolecules. Recent simulations that address these questions are reviewed here and discussed in the context of relevant experiments.

Explicit Incorporation of Hard and Soft Protein-Protein Interactions into Models for Crowding Effects in Protein Mixtures. 2. Effects of Varying Hard and Soft Interactions upon Prototypical Chemical Equilibria

Previously derived approximate analytical relations for the activity coefficient of each solute in a mixture of up to three spherical solutes in a highly nonideal solution interacting via square well potentials of mean force (Hoppe, T.; Minton, A. P. J Phys Chem B. 2016, 120, 11866-11872) were used to explore the effect of heterogeneity in volume occupancy and intermolecular interactions upon prototypical schemes representing solubility, partitioning, conformational isomerization, and self-association in crowded solutions. Results generally indicate that all of the equilibria explored are exquisitely sensitive to variations in both volume occupancy and intermolecular interaction and have important implications for the design and execution of more detailed simulations of complex media.

Role of solvent properties of water in crowding effects induced by macromolecular agents and osmolytes

Dipolarity/polarizability of water in polymer mixtures may be additive, reduced or enhanced depending on the composition of the mixture.