Confinement and Crowding Effects on Folding of a Multidomain Y-Family DNA Polymerase

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.

Impact of Confinement within a Hydrogel Mesh on Protein Thermodynamic Stability and Aggregation Kinetics

Though protein stability and aggregation have been well characterized in dilute solutions, the influence of a confining environment that exists (e.g., in intercellular and tissue spaces and therapeutic formulations) on the protein structure is largely unknown. Herein, the effects of confinement on stability and aggregation were explored for proteins of different sizes, stability, and hydrophobicity when encapsulated in hydrophilic poly(ethylene glycol) hydrogels. Denaturation curves show linear correlations between confinement size (mesh size) and thermodynamic stability, i.e., unfolding free energy and surface area accessible for solvation (represented by m-value). Two clusters of protein types are identifiable from these correlations; the clusters are defined by differences in protein stability, surface area, and aggregation propensity. Proteins with higher stability, larger surface area, and lower aggregation propensity (e.g., lysozyme and myoglobin) are less affected by the confinement imposed by mesh size than proteins with lower stability, smaller surface area, and higher aggregation propensity (e.g., growth hormone and aldehyde dehydrogenase). According to aggregation kinetics measured by thioflavin T fluorescence, confinement in smaller mesh sizes resulted in slower aggregation rates than that in larger mesh sizes. Compared to that in buffer solution, the confinement of a hydrophobic protein (e.g., human insulin) in the hydrogels accelerates protein aggregation. Conversely, the confinement of a hydrophilic protein (e.g., human amylin) in the hydrogels decelerates or prevents aggregation, with the rates of aggregation inversely proportional to mesh size. These findings provide new insights into protein conformational stability in confined microenvironments relevant to various cellular, tissue, and therapeutics scenarios.

Unveiling Metastable Ensembles of GRB2 and the Relevance of Interdomain Communication during Folding

The folding process of multidomain proteins is a highly intricate phenomenon involving the assembly of distinct domains into a functional three-dimensional structure. During this process, each domain may fold independently while interacting with others. The folding of multidomain proteins can be influenced by various factors, including their composition, the structure of each domain, or the presence of disordered regions, as well as the surrounding environment. Misfolding of multidomain proteins can lead to the formation of nonfunctional structures associated with a range of diseases, including cancers or neurodegenerative disorders. Understanding this process is an important step for many biophysical analyses such as stability, interaction, malfunctioning, and rational drug design. One such multidomain protein is growth factor receptor-bound protein 2 (GRB2), an adaptor protein that is essential in regulating cell survival. GRB2 consists of one central Src homology 2 (SH2) domain flanked by two Src homology 3 (SH3) domains. The SH2 domain interacts with phosphotyrosine regions in other proteins, while the SH3 domains recognize proline-rich regions on protein partners during cell signaling. Here, we combined computational and experimental techniques to investigate the folding process of GRB2. Through computational simulations, we sampled the conformational space and mapped the mechanisms involved by the free energy profiles, which may indicate possible intermediate states. From the molecular dynamics trajectories, we used the energy landscape visualization method (ELViM), which allowed us to visualize a three-dimensional (3D) representation of the overall energy surface. We identified two possible parallel folding routes that cannot be seen in a one-dimensional analysis, with one occurring more frequently during folding. Supporting these results, we used differential scanning calorimetry (DSC) and fluorescence spectroscopy techniques to confirm these intermediate states in vitro. Finally, we analyzed the deletion of domains to compare our model outputs to previously published results, supporting the presence of interdomain modulation. Overall, our study highlights the significance of interdomain communication within the GRB2 protein and its impact on the formation, stability, and structural plasticity of the protein, which are crucial for its interaction with other proteins in key signaling pathways.

Perspectives on the landscape and flux theory for describing emergent behaviors of the biological systems

We give a review on the landscape theory of the equilibrium biological systems and landscape-flux theory of the nonequilibrium biological systems as the global driving force. The emergences of the behaviors, the associated thermodynamics in terms of the entropy and free energy and dynamics in terms of the rate and paths have been quantitatively demonstrated. The hierarchical organization structures have been discussed. The biological applications ranging from protein folding, biomolecular recognition, specificity, biomolecular evolution and design for equilibrium systems as well as cell cycle, differentiation and development, cancer, neural networks and brain function, and evolution for nonequilibrium systems, cross-scale studies of genome structural dynamics and experimental quantifications/verifications of the landscape and flux are illustrated. Together, this gives an overall global physical and quantitative picture in terms of the landscape and flux for the behaviors, dynamics and functions of biological systems.

Investigating the Conformational Dynamics of a Y-Family DNA Polymerase during Its Folding and Binding to DNA and a Nucleotide

During DNA polymerization, the Y-family DNA polymerases are capable of bypassing various DNA damage, which can stall the replication fork progression. It has been well acknowledged that the structures of the Y-family DNA polymerases have been naturally evolved to undertake this vital task. However, the mechanisms of how these proteins utilize their unique structural and conformational dynamical features to perform the translesion DNA synthesis are less understood. Here, we developed structure-based models to study the precatalytic DNA polymerization process, including DNA and nucleotide binding to DPO4, a paradigmatic Y-family polymerase from Sulfolobus solfataricus. We studied the interplay between the folding and the conformational dynamics of DPO4 and found that DPO4 undergoes first unraveling (unfolding) and then folding for accomplishing the functional "open-to-closed" conformational transition. DNA binding dynamically modulates the conformational equilibrium in DPO4 during the stepwise binding through different types of interactions, leading to different conformational distributions of DPO4 at different DNA binding stages. We observed that nucleotide binding induces modulation of a few contacts surrounding the active site of the DPO4-DNA complex associated with a high free energy barrier. Our simulation results resonate with the experimental evidence that the conformational change at the active site led by nucleotide is the rate-limiting step of nucleotide incorporation. In combination with localized frustration analyses, we underlined the importance of DPO4 conformational dynamics and fluctuations in facilitating DNA and nucleotide binding. Our findings offer mechanistic insights into the processes of DPO4 conformational dynamics associated with the substrate binding and contribute to the understanding of the "structure-dynamics-function" relationship in the Y-family DNA polymerases.

Protein assembly and crowding simulations

Proteins encounter frequent molecular interactions in biological environments. Computer simulations have become an increasingly important tool in providing mechanistic insights into how such interactions in vivo relate to their biological function. The review here focuses on simulations describing protein assembly and molecular crowding effects as two important aspects that are distinguished mainly by how specific and long-lived protein contacts are. On the topic of crowding, recent simulations have provided new insights into how crowding affects protein folding and stability, modulates enzyme activity, and affects diffusive properties. Recent studies of assembly processes focus on assembly pathways, especially for virus capsids, amyloid aggregation pathways, and the role of multivalent interactions leading to phase separation. Also, discussed are technical challenges in achieving increasingly realistic simulations of interactions in cellular environments.

A small molecule high throughput screening platform to profile conformational properties of nascent, ribosome-bound proteins

Genetic mutations cause a wide spectrum of human disease by disrupting protein folding, both during and after synthesis. Transient de-novo folding intermediates therefore represent potential drug targets for pharmacological correction of protein folding disorders. Here we develop a FRET-based high-throughput screening (HTS) assay in 1,536-well format capable of identifying small molecules that interact with nascent polypeptides and correct genetic, cotranslational folding defects. Ribosome nascent chain complexes (RNCs) containing donor and acceptor fluorophores were isolated from cell free translation reactions, immobilized on Nickel-NTA/IDA beads, and imaged by high-content microscopy. Quantitative FRET measurements obtained from as little as 0.4 attomole of protein/bead enabled rapid assessment of conformational changes with a high degree of reproducibility. Using this assay, we performed a pilot screen of ~ 50,000 small molecules to identify compounds that interact with RNCs containing the first nucleotide-binding domain (NBD1) of the cystic fibrosis transmembrane conductance regulator (CFTR) harboring a disease-causing mutation (A455E). Screen results yielded 133 primary hits and 1 validated hit that normalized FRET values of the mutant nascent peptide. This system provides a scalable, tractable, structure-based discovery platform for screening small molecules that bind to or impact the folding of protein substrates that are not amenable to traditional biochemical analyses.

Topological frustration leading to backtracking in a coupled folding-binding process

Intrinsically disordered proteins (IDPs) are abundant in all species. Their discovery challenges the traditional "sequence-structure-function" paradigm of protein science because IDPs play important roles in various biological processes without preformed folded structures. Bioinformatic analysis reveals that the intrinsically conformational disorder of IDPs as well as their conformational transition upon binding to their targets is encoded by their amino acid sequences. The rRNase domain of colicin E3 and the immunity protein Im3 are a pair of proteins involved in bacterial survival. While the N-terminal segment and the central segment of E3 make comparable intermolecular contacts with Im3 in the bound state, binding of E3 with Im3 is dominantly triggered by the central segment of E3. In this work, to further investigate the binding mechanism of disordered E3 with Im3, we performed systematic free energy and transition path analysis through coarse-grained molecular dynamics simulations. We observed backtracking of the N-terminal segment of E3 in the binding process, whose occurrence depends on salt concentration. Conformational analysis revealed that initial binding of the N-terminal segment of E3 to Im3 usually leads to misorientation of a central hairpin of E3 on Im3, which generates topological frustration and results in backtracking of the N-terminal segment. Our results not only provide deeper mechanistic insights into the coupled folding-binding process of the E3/Im3 complex, but also suggest that topological frustration could be present in the coupled folding-binding process of IDPs and play an important role in regulating the binding transition pathways.

Physics of biomolecular recognition and conformational dynamics

Biomolecular recognition usually leads to the formation of binding complexes, often accompanied by large-scale conformational changes. This process is fundamental to biological functions at the molecular and cellular levels. Uncovering the physical mechanisms of biomolecular recognition and quantifying the key biomolecular interactions are vital to understand these functions. The recently developed energy landscape theory has been successful in quantifying recognition processes and revealing the underlying mechanisms. Recent studies have shown that in addition to affinity, specificity is also crucial for biomolecular recognition. The proposed physical concept of intrinsic specificity based on the underlying energy landscape theory provides a practical way to quantify the specificity. Optimization of affinity and specificity can be adopted as a principle to guide the evolution and design of molecular recognition. This approach can also be used in practice for drug discovery using multidimensional screening to identify lead compounds. The energy landscape topography of molecular recognition is important for revealing the underlying flexible binding or binding-folding mechanisms. In this review, we first introduce the energy landscape theory for molecular recognition and then address four critical issues related to biomolecular recognition and conformational dynamics: (1) specificity quantification of molecular recognition; (2) evolution and design in molecular recognition; (3) flexible molecular recognition; (4) chromosome structural dynamics. The results described here and the discussions of the insights gained from the energy landscape topography can provide valuable guidance for further computational and experimental investigations of biomolecular recognition and conformational dynamics.

Stability of proteins encapsulated in Michael-type addition polyethylene glycol hydrogels

Degradable polyethylene glycol (PEG) hydrogels are excellent vehicles for sustained drug release due to their biocompatibility, tunable physical properties, and customizable degradation. However, protein therapeutics are unstable under physiological conditions and releasing degraded or inactive therapeutics can induce immunogenic effects. While controlling protein release from PEG hydrogels has been extensively investigated, few studies have detailed protein stability long-term or under stress conditions. Here, lysozyme and alcohol dehydrogenase (ADH) stability were explored upon encapsulation in PEG hydrogels formed through Michael-type addition. The stability and structure of the two model proteins were monitored by measuring the free energy of unfolding and fluoresce quenching when confined in a hydrogel and compared to PEG solution and buffer. Hydrogels destabilized lysozyme structure at low denaturant concentrations but prevented complete unfolding at high concentrations. ADH was stabilized as the confining mesh size approached the protein radius of gyration. Both proteins retained enzymatic activity within the hydrogels under stress conditions, including denaturant, high temperature, and agitation. Conjugation between lysozyme and PEG-acrylate was identified at long reaction times but no conjugation was observed in the time required for complete gelation. Studies of protein stability in PEG hydrogels, as the one detailed here, can lead to designer technologies for the improved formulation, storage, and delivery of protein therapeutics.

Effects of electrostatic interactions on global folding and local conformational dynamics of a multidomain Y-family DNA polymerase

The Y-family DNA polymerases specialize in translesion DNA synthesis, which is essential for replicating damaged DNA. The Y-family polymerases, which are made up of four stable domains, exhibit extensive distributions of charged residues, and are responsible for the tight formation of the protein-DNA complex. However, it is still unclear how the electrostatic interactions influence the conformational dynamics of the polymerases. Here, we focus on the case of a prototype Y-family DNA polymerase, Dpo4. Using coarse-grained models including a salt-dependent electrostatic potential, we investigate the effects of the electrostatic interactions on the folding process of Dpo4. Our simulations show that strong electrostatic interactions result in a three-state folding of Dpo4, consistent with the experimental observations. This folding process exhibits low cooperativity led by low salt concentration, where the individual domains fold one by one through one single pathway. Since the refined folding order of domains in multidomain proteins can shrink the configurational space, we suggest that the electrostatic interactions facilitate the Dpo4 folding. In addition, we study the local conformational dynamics of Dpo4 in terms of fluctuation and frustration analyses. We show that the electrostatic interactions can exaggerate the local conformational properties, which are in favor of the large-scale conformational transition of Dpo4 during the functional DNA binding. Our results underline the importance of electrostatic interactions in the conformational dynamics of Dpo4 at both the global and local scale, providing useful guidance in protein engineering at the multidomain level.

It is time to crowd your cell culture media - Physicochemical considerations with biological consequences

In vivo, the interior and exterior of cells is populated by various macromolecules that create an extremely crowded milieu. Yet again, in vitro eukaryotic cell culture is conducted in dilute culture media that hardly imitate the native tissue density. Herein, the concept of macromolecular crowding is discussed in both intracellular and extracellular context. Particular emphasis is given on how the physicochemical properties of the crowding molecules govern and determine kinetics, equilibria and mechanism of action of biochemical and biological reactions, processes and functions. It is evidenced that we are still at the beginning of appreciating, let alone effectively implementing, the potential of macromolecular crowding in permanently differentiated and stem cell culture systems.

Combined Effects of Confinement and Macromolecular Crowding on Protein Stability

Confinement and crowding have been shown to affect protein fates, including folding, functional stability, and their interactions with self and other proteins. Using both theoretical and experimental studies, researchers have established the independent effects of confinement or crowding, but only a few studies have explored their effects in combination; therefore, their combined impact on protein fates is still relatively unknown. Here, we investigated the combined effects of confinement and crowding on protein stability using the pores of agarose hydrogels as a confining agent and the biopolymer, dextran, as a crowding agent. The addition of dextran further stabilized the enzymes encapsulated in agarose; moreover, the observed increases in enhancements (due to the addition of dextran) exceeded the sum of the individual enhancements due to confinement and crowding. These results suggest that even though confinement and crowding may behave differently in how they influence protein fates, these conditions may be combined to provide synergistic benefits for protein stabilization. In summary, our study demonstrated the successful use of polymer-based platforms to advance our understanding of how in vivo like environments impact protein function and structure.

Investigating the trade-off between folding and function in a multidomain Y-family DNA polymerase

The way in which multidomain proteins fold has been a puzzling question for decades. Until now, the mechanisms and functions of domain interactions involved in multidomain protein folding have been obscure. Here, we develop structure-based models to investigate the folding and DNA-binding processes of the multidomain Y-family DNA polymerase IV (DPO4). We uncover shifts in the folding mechanism among ordered domain-wise folding, backtracking folding, and cooperative folding, modulated by interdomain interactions. These lead to ‘U-shaped’ DPO4 folding kinetics. We characterize the effects of interdomain flexibility on the promotion of DPO4–DNA (un)binding, which probably contributes to the ability of DPO4 to bypass DNA lesions, which is a known biological role of Y-family polymerases. We suggest that the native topology of DPO4 leads to a trade-off between fast, stable folding and tight functional DNA binding. Our approach provides an effective way to quantitatively correlate the roles of protein interactions in conformational dynamics at the multidomain level.