Interaction of Polyelectrolytes with Proteins: Quantifying the Role of Water

Quantitative Prediction of Protein-Polyelectrolyte Binding Thermodynamics: Adsorption of Heparin-Analog Polysulfates to the SARS-CoV-2 Spike Protein RBD

Interactions of polyelectrolytes (PEs) with proteins play a crucial role in numerous biological processes, such as the internalization of virus particles into host cells. Although docking, machine learning methods, and molecular dynamics (MD) simulations are utilized to estimate binding poses and binding free energies of small-molecule drugs to proteins, quantitative prediction of the binding thermodynamics of PE-based drugs presents a significant obstacle in computer-aided drug design. This is due to the sluggish dynamics of PEs caused by their size and strong charge-charge correlations. In this paper, we introduce advanced sampling methods based on a force-spectroscopy setup and theoretical modeling to overcome this barrier. We exemplify our method with explicit solvent all-atom MD simulations of the interactions between anionic PEs that show antiviral properties, namely heparin and linear polyglycerol sulfate (LPGS), and the SARS-CoV-2 spike protein receptor binding domain (RBD). Our prediction for the binding free-energy of LPGS to the wild-type RBD matches experimentally measured dissociation constants within thermal energy, k B T, and correctly reproduces the experimental PE-length dependence. We find that LPGS binds to the Delta-variant RBD with an additional free-energy gain of 2.4 k B T, compared to the wild-type RBD, due to the additional presence of two mutated cationic residues contributing to the electrostatic energy gain. We show that the LPGS-RBD binding is solvent dominated and enthalpy driven, though with a large entropy-enthalpy compensation. Our method is applicable to general polymer adsorption phenomena and predicts precise binding free energies and reconfigurational friction as needed for drug and drug-delivery design.

Interaction between the Polyelectrolytes Unfractionated Heparin and Universal Heparin Reversal Agents

The interaction of unfractionated heparin (UFH) with universal heparin reversal agent 7 (UHRA-7) is investigated. UHRA-7 is composed of a hyperbranched polyglycerol core onto which an array of methylated tris(2-aminoethylamine) (Me-TREN) charged groups is grafted, which in turn are shielded with a layer of small chain poly(ethylene glycol) methyl ether (mPEG) chains. This system has previously been shown to be biocompatible and to be effective at neutralizing heparin. The binding constant Kb was determined from isothermal titration calorimetry experiments, at temperatures ranging from 278 to 323 K and salt concentrations ranging from 0.06 to 0.20 M NaCl. The data were analyzed in terms of a number of different theoretical models to determine the contribution of counterion release and water release to driving the interaction between UFH and UHRA-7. With the support of NMR and molecular dynamics simulation data, a model of the interaction between UFH and UHRA-7 is proposed. The binding of heparin and universal heparin reversal agent 7 is mainly due to charge-charge interactions between the negatively charged units on UFH with positively charged Me-TREN and mPEG chains.

Driving Forces in the Formation of Biocondensates of Highly Charged Proteins: A Thermodynamic Analysis of the Binary Complex Formation

A thermodynamic analysis of the binary complex formation of the highly positively charged linker histone H1 and the highly negatively charged chaperone prothymosin α (ProTα) is detailed. ProTα and H1 have large opposite net charges (-44 and +53, respectively) and form complexes at physiological salt concentrations with high affinities. The data obtained for the binary complex formation are analyzed by a thermodynamic model that is based on counterion condensation modulated by hydration effects. The analysis demonstrates that the release of the counterions mainly bound to ProTα is the main driving force, and effects related to water release play no role within the limits of error. A strongly negative Δcp (=-0.87 kJ/(K mol)) is found, which is due to the loss of conformational degrees of freedom.

Recent advances in wolfberry polysaccharides and whey protein-based biopolymers for regulating the diversity of gut microbiota and its mechanism: A review

Imbalances in gut microbiota diversity are associated with various health issues, including obesity and related disorders. There is a growing interest in developing synergistic biopolymers based on wolfberry polysaccharides and whey protein to address these problems due to their potential health benefits. This review explores recent advances in understanding how functional foods based on Lycium barbarum polysaccharides (LBP) and whey protein (WP) influence gut microbiota diversity and their underlying mechanisms. We examine the impact of these biopolymers on microbial composition and functionality, focusing on their roles in improving health by regulating gut microbiota. The combined effects of WP and LBP significantly enhance gut microbiome metabolic activities and taxonomic diversity, offering promising avenues for treating obesity and related disorders.

Polymers for Disrupting Protein-Protein Interactions: Where Are We and Where Should We Be?

Protein-protein interactions (PPIs) are central to the cellular signaling and regulatory networks that underlie many physiological and pathophysiological processes. It is challenging to target PPIs using traditional small molecule or peptide-based approaches due to the frequent lack of well-defined binding pockets at the large and flat PPI interfaces. Synthetic polymers offer an opportunity to circumvent these challenges by providing unparalleled flexibility in tuning their physiochemical properties to achieve the desired binding properties. In this review, we summarize the current state of the field pertaining to polymer-protein interactions in solution, highlighting various polyelectrolyte systems, their tunable parameters, and their characterization. We provide an outlook on how these architectures can be improved by incorporating sequence control, foldability, and machine learning to mimic proteins at every structural level. Advances in these directions will enable the design of more specific protein-binding polymers and provide an effective strategy for targeting dynamic proteins, such as intrinsically disordered proteins.

Interactions between Hyaluronic Acid and Chitosan by Isothermal Titration Calorimetry: The Effect of Ionic Strength, pH, and Polymer Molecular Weight

Hyaluronic acid (HA)/chitosan (CHI) complex coacervates have recently gained interest due to the pH-dependent ionization and semiflexibility of the polymers as well as their applicability in tissue engineering. Here, we apply isothermal titration calorimetry (ITC) to understand the apparent thermodynamics of coacervation for HA/CHI as a function of the pH, ionic strength, and chain length. We couple these ITC experiments with the knowledge of the charge states of HA and CHI from potentiometric titration to understand the mechanistic aspects of complex formation. Our data demonstrate that the driving force for the complex coacervation of HA and CHI is entropic in nature and this driving force decreased with increasing ionic strength. We also observed a decrease in the stoichiometry for ion-pairing with increasing ionic strength, which we suggest is a consequence of the changing degree of ionization for HA at higher ionic strengths. An increase in the strength of interactions with pH was hypothesized to also be a result of changes in the degree of ionization of HA, though stronger interactions were observed at the lowest pH tested, likely due to contributions from hydrogen bonding between HA and CHI.

Exploring the Origins of Association of Poly(acrylic acid) Polyelectrolyte with Lysozyme in Aqueous Environment through Molecular Simulations and Experiments

This study provides a detailed picture of how a protein (lysozyme) complexes with a poly(acrylic acid) polyelectrolyte (PAA) in water at the atomic level using a combination of all-atom molecular dynamics simulations and experiments. The effect of PAA and temperature on the protein's structure is explored. The simulations reveal that a lysozyme's structure is relatively stable except from local conformational changes induced by the presence of PAA and temperature increase. The effect of a specific thermal treatment on the complexation process is investigated, revealing both structural and energetic changes. Certain types of secondary structures (i.e., α-helix) are found to undergo a partially irreversible shift upon thermal treatment, which aligns qualitatively with experimental observations. This uncovers the origins of thermally induced aggregation of lysozyme with PAA and points to new PAA/lysozyme bonds that are formed and potentially enhance the stability in the complexes. As the temperature changes, distinct amino acids are found to exhibit the closest proximity to PAA, resulting into different PAA/lysozyme interactions; consequently, a different complexation pathway is followed. Energy calculations reveal the dominant role of electrostatic interactions. This detailed information can be useful for designing new biopolymer/protein materials and understanding protein function under immobilization of polyelectrolytes and upon mild denaturation processes.

Entropy Drives Interpolymer Association in Water: Insights into Molecular Mechanisms

Interpolymer association in aqueous solutions is essential for many industrial processes, new materials design, and the biochemistry of life. However, our understanding of the association mechanism is limited. Classical theories do not provide molecular details, creating a need for detailed mechanistic insights. This work consolidates previous literature with complementary isothermal titration calorimetry (ITC) measurements and molecular dynamics (MD) simulations to investigate molecular mechanisms to provide such insights. The large body of ITC data shows that intermolecular bonds, such as ionic or hydrogen bonds, cannot drive association. Instead, polymer association is entropy-driven due to the reorganization of water and ions. We propose a unifying entropy-driven association mechanism by generalizing previously suggested polyion association principles to include nonionic polymers, here termed polydipoles. In this mechanism, complementary charge densities of the polymers are the common denominators of association, for both polyions and polydipoles. The association of the polymers results mainly from two processes: charge exchange and amphiphilic association. MD simulations indicate that the amphiphilic assembly alone is enough for the initial association. Our proposed mechanism is a step toward a molecular understanding of the formation of complexes between synthetic and biological polymers under ambient or biological conditions.

Effects of Control Factors on Protein-Polyelectrolyte Complex Coacervation

Protein-polyelectrolyte complex coacervation is of particular interest for mimicking intracellular phase separation and organization. Yet, the challenge arises from regulating the coacervation due to the globular structure and anisotropic distributed charges of protein. Herein, we fully investigate the different control factors and reveal their effects on protein-polyelectrolyte coacervation. We prepared mixtures of BSA (bovine serum albumin) with different cationic polymers, which include linear and branched polyelectrolytes covering different spacer and charge groups, chain lengths, and polymer structures. With BSA-PDMAEMA [poly(N,N-dimethylaminomethyl methacrylate)] as the main investigated pair, we find that the moderate pH and ionic strength are essential for the adequate electrostatic interaction and formation of coacervate droplets. For most BSA-polymer mixtures, excess polyelectrolytes are required to achieve the full complexation, as evidenced by the deviated optimal charge mixing ratios from the charge stoichiometry. Polymers with longer chains or primary amine groups and a branched structure endow a strong electrostatic interaction with BSA and cause a bigger charge ratio deviation associated with the formation of solid-like coacervate complexes. Nevertheless, both the liquid- and solid-like coacervates hardly interrupt the BSA structure and activity, indicating the safe encapsulation of proteins by the coacervation with polyelectrolytes. Our study validates the crucial control of the diverse factors in regulating protein-polyelectrolyte coacervation, and the revealed principles shall be instructive for establishing other protein-based coacervations and boosting their potential applications.

Protein separation by sequential selective complex coacervation

In food manufacturing and particular biomedical products selected proteins are often required. Obtaining the desired proteins in a pure form from natural resources is therefore important, but often very challenging. Herein, we design a sequential coacervation process that allows to efficiently isolate and purify proteins with different isoelectric points (pIs) from a mixed solution, namely Bovine Serum Albumin (BSA, pI = 4.9) and Peroxidase from Horseradish (HRP, pI = 7.2). The key to separation is introducing a suitable polyelectrolyte that causes selective complex coacervation at appropriate pH and ionic strength. Specifically, polyethyleneimine (PEI), when added into the mixture at pH 6.0, produces a coacervation which exclusively contains BSA, leading to a supernatant solution containing 100 % HRP with a purity of 91 %. After separating the dilute and dense phases, BSA is recovered by adding poly(acrylic acid) (PAA) to the concentrated phase, which displaces BSA from the complex because it interacts more strongly with PEI. The supernatant phase after this step contains approximately 75 % of the initial amount of BSA with a purity of 99 %. Our results confirm that coacervation under well-defined conditions can be selective, enabling separation of proteins with adequate purity. Therefore, the established approach demonstrates a facile and sustainable strategy with potential for protein separation at industrial scale.

Respiratory viruses interacting with cells: the importance of electrostatics

The COVID-19 pandemic has rekindled interest in the molecular mechanisms involved in the early steps of infection of cells by viruses. Compared to SARS-CoV-1 which only caused a relatively small albeit deadly outbreak, SARS-CoV-2 has led to fulminant spread and a full-scale pandemic characterized by efficient virus transmission worldwide within a very short time. Moreover, the mutations the virus acquired over the many months of virus transmission, particularly those seen in the Omicron variant, have turned out to result in an even more transmissible virus. Here, we focus on the early events of virus infection of cells. We review evidence that the first decisive step in this process is the electrostatic interaction of the spike protein with heparan sulfate chains present on the surface of target cells: Patches of cationic amino acids located on the surface of the spike protein can interact intimately with the negatively charged heparan sulfate chains, which results in the binding of the virion to the cell surface. In a second step, the specific interaction of the receptor binding domain (RBD) within the spike with the angiotensin-converting enzyme 2 (ACE2) receptor leads to the uptake of bound virions into the cell. We show that these events can be expressed as a semi-quantitative model by calculating the surface potential of different spike proteins using the Adaptive Poison-Boltzmann-Solver (APBS). This software allows visualization of the positive surface potential caused by the cationic patches, which increased markedly from the original Wuhan strain of SARS-CoV-2 to the Omicron variant. The surface potential thus enhanced leads to a much stronger binding of the Omicron variant as compared to the original wild-type virus. At the same time, data taken from the literature demonstrate that the interaction of the RBD of the spike protein with the ACE2 receptor remains constant within the limits of error. Finally, we briefly digress to other viruses and show the usefulness of these electrostatic processes and calculations for cell-virus interactions more generally.

Sulfated Hyperbranched and Linear Polyglycerols Modulate HMGB1 and Morphological Plasticity in Neural Cells

The objective of this study was to establish if polyglycerols with sulfate or sialic acid functional groups interact with high mobility group box 1 (HMGB1), and if so, which polyglycerol could prevent loss of morphological plasticity in excitatory neurons in the hippocampus. Considering that HMGB1 binds to heparan sulfate and that heparan sulfate has structural similarities with dendritic polyglycerol sulfates (dPGS), we performed the experiments to show if polyglycerols can mimic heparin functions by addressing the following questions: (1) do dendritic and linear polyglycerols interact with the alarmin molecule HMGB1? (2) Does dPGS interaction with HMGB1 influence the redox status of HMGB1? (3) Can dPGS prevent the loss of dendritic spines in organotypic cultures challenged with lipopolysaccharide (LPS)? LPS plays a critical role in infections with Gram-negative bacteria and is commonly used to test candidate therapeutic agents for inflammation and endotoxemia. Pathologically high LPS concentrations and other stressful stimuli cause HMGB1 release and post-translational modifications. We hypothesized that (i) electrostatic interactions of hyperbranched and linear polysulfated polyglycerols with HMGB1 will likely involve sites similar to those of heparan sulfate. (ii) dPGS can normalize HMGB1 compartmentalization in microglia exposed to LPS and prevent dendritic spine loss in the excitatory hippocampal neurons. We performed immunocytochemistry and biochemical analyses combined with confocal microscopy to determine cellular and extracellular locations of HMGB1 and morphological plasticity. Our results suggest that dPGS interacts with HMGB1 similarly to heparan sulfate. Hyperbranched dPGS and linear sulfated polymers prevent dendritic spine loss in hippocampal excitatory neurons. MS/MS analyses reveal that dPGS-HMGB1 interactions result in fully oxidized HMGB1 at critical cysteine residues (Cys23, Cys45, and Cys106). Triply oxidized HMGB1 leads to the loss of its pro-inflammatory action and could participate in dPGS-mediated spine loss prevention. LPG-Sia exposure to HMGB1 results in the oxidation of Cys23 and Cys106 but does not normalize spine density.

Removal of nonspecific binding proteins is required in co-immunoprecipitation with nuclear proteins

Whether protein samples should be pretreated to remove nonspecific binding proteins in co-immunoprecipitation (CO-IP) is controversial. In this work, nonspecific binding of proteins to agarose beads was found to be greater than that to magnetic beads. The nonspecific binding was increased with the decrease of ion concentrations but reduced by Nonidet P40. Western blot indicated that p65 and β-actin were present as nonspecifically bound protein to the beads. p53 and β-actin were present in the CO-IP precipitates of nuclear proteins but pretreatment cleared the nonspecifically pulled down p53 and β-actin. These data suggest that magnetic beads are better for CO-IP, but preclearing is necessary to minimize false positive regardless of which bead is used, particularly for nuclear proteins.

Interaction of Heparin with Proteins: Hydration Effects

We present a thermodynamic investigation of the interaction of heparin with lysozyme in the presence of potassium glutamate (KGlu). The binding constant K_b is measured by isothermal titration calorimetry (ITC) in a temperature range from 288 to 310 K for concentrations of KGlu between 25 and 175 mM. The free energy of binding ΔG_b derived from K_b is strongly decreasing with increasing concentration of KGlu, whereas the dependence of ΔG_b on temperature T is found to be small. The decrease of ΔG_b can be explained in terms of counterion release: Binding of lysozyme to the strong polyelectrolyte heparin liberates approximately three of the condensed counterions of heparin, thus increasing the entropy of the system. The dependence of ΔG_b on T, on the other hand, is traced back to a change of hydration of the protein and the polyelectrolyte upon complex formation. This dependence is quantitatively described by the parameter Δw that depends on T and vanishes at a characteristic temperature T₀. A comparison of the complex formation in the presence of KGlu with the one in the presence of NaCl demonstrates that the parameters related to hydration are changed considerably. The characteristic temperature T₀ in the presence of KGlu solutions is considerably smaller than that in the presence of NaCl solutions. The change of specific heat Δc_p is found to become more negative with increasing salt concentration: This finding agrees with the model-free analysis by the generalized van't Hoff equation. The entire analysis reveals a small but important change of the free energy of binding by hydration. It shows that these ion-specific Hofmeister effects can be modeled quantitatively in terms of a characteristic temperature T₀ and a parameter describing the dependence of Δc_p on salt concentration.

Denaturation of proteins: electrostatic effects vs. hydration

The unfolding transition of proteins in aqueous solution containing various salts or uncharged solutes is a classical subject of biophysics. In many cases, this transition is a well-defined two-stage equilibrium process which can be described by a free energy of transition ΔG u and a transition temperature T m. For a long time, it has been known that solutes can change T m profoundly. Here we present a phenomenological model that describes the change of T m with the solute concentration c s in terms of two effects: (i) the change of the number of correlated counterions Δn ci and (ii) the change of hydration expressed through the parameter Δw and its dependence on temperature expressed through the parameter dΔc p/dc s. Proteins always carry charges and Δn ci describes the uptake or release of counterions during the transition. Likewise, the parameter Δw measures the uptake or release of water during the transition. The transition takes place in a reservoir with a given salt concentration c s that defines also the activity of water. The parameter Δn ci is a measure for the gain or loss of free energy because of the release or uptake of ions and is related to purely entropic effects that scale with ln c s. Δw describes the effect on ΔG u through the loss or uptake of water molecules and contains enthalpic as well as entropic effects that scale with c s. It is related to the enthalpy of transition ΔH u through a Maxwell relation: the dependence of ΔH u on c s is proportional to the dependence of Δw on temperature. While ionic effects embodied in Δn ci are independent of the kind of salt, the hydration effects described through Δw are directly related to Hofmeister effects of the various salt ions. A comparison with literature data underscores the general validity of the model.

Interaction of Linear Polyelectrolytes with Proteins: Role of Specific Charge-Charge Interaction and Ionic Strength

We present a thermodynamic study of the interaction of synthetic, linear polyelectrolytes with bovine serum albumin (BSA). All polyelectrolytes are based on poly(allyl glycidyl ether) which has been modified by polymer-analogous reaction with anionic (-SO₃Na), cationic (-NH₃Cl or -NHMe₂Cl) or zwitterionic groups (-NMe₂(CH₂)₃SO₃). While the anionic polymer shows a very weak interaction, the zwitterionic polymer exhibits no interaction with BSA (pI = 4.7) under the applied pH = 7.4, ionic strength (I = 23–80 mM) and temperature conditions (T = 20–37 °C). A strong binding, however, was observed for the polycations bearing primary amino or tertiary dimethyl amino groups, which could be analysed in detail by isothermal titration calorimetry (ITC). The analysis was done using an expression which describes the free energy of binding, ΔG_b, as the function of the two decisive variables, temperature, T, and salt concentration, c_s. The underlying model splits ΔG_b into a term related to counterion release and a term related to water release. While the number of released counter ions is similar for both systems, the release of bound water is more important for the primary amine compared to the tertiary N,N-dimethyl amine presenting polymer. This finding is further traced back to a closer contact of the polymers’ protonated primary amino groups in the complex with oppositely charged moieties of BSA as compared to the bulkier protonated tertiary amine groups. We thus present an investigation that quantifies both driving forces for electrostatic binding, namely counterion release and change of hydration, which contribute to a deeper understanding with direct impact on future advancements in the biomedical field.