Charge Regulation Effects in Polyelectrolyte Adsorption.

Polyelectrolyte (PE) adsorption plays a pivotal role in tailoring surface properties, finding diverse applications across scientific and industrial domains. In aqueous environments, these polymers gain charge through the dissociation equilibrium of ionizable groups. Consequently, accurately representing their electrostatic interactions necessitates considering many-body couplings between charge and configuration-a phenomenon known as charge regulation (CR)-instead of relying on the constant charge approximation. Here, we show that CR effects can significantly enhance PE adsorption. In the vicinity of a like-charged interface, where short-range PE-surface attraction drives adsorption, both PEs and the surface undergo a reduction in charges during the adsorption process, especially in a semidilute PE solution. This reduction in charges leads to a noticeable decrease in the electrostatic repulsion barrier, facilitating the adsorption process. Conversely, near an oppositely charged surface, CR effects enhance the electrostatic charges of both PEs and surfaces coherently, suppressing the charge reversal of surface charges and ultimately promoting adsorption. These findings underscore the crucial CR effects in PE adsorption and challenge the direct use of electrostatic charges, typically estimated in bulk systems, for an accurate understanding of adsorption phenomena.

Collapse and expansion kinetics of a single polyelectrolyte chain with hydrodynamic interactions.

We investigate the collapse and expansion dynamics of a linear polyelectrolyte (PE) with hydrodynamic interactions. Using dissipative particle dynamics with a bead-spring PE model, long-range electrostatics, and explicit ions, we examine how the timescales of collapse tcol and expansion texp depend on the chain length N and obtain scaling relationships tcol ∼ Nα and texp ∼ Nß. For neutral polymers, we derive values of α = 0.94 ± 0.01 and ß = 1.97 ± 0.10. Interestingly, the introduction of electrostatic interaction markedly shifts α to α ≈ 1.4 ± 0.1 for salt concentrations within c = 10-4 to 10-2 M. A reduction in the ion-to-monomer size ratio noticeably reduces α. On the other hand, the expansion scaling remains approximately constant, ß ≈ 2, regardless of the salt concentration or ion size considered. We find ß > α for all conditions considered, implying that expansion is always slower than collapse in the limit of long polymers. This asymmetry is explained by distinct kinetic pathways of collapse and expansion processes.

Impact of Hydrodynamic Interactions on the Kinetic Pathway of Protein Folding.

Protein folding is a fundamental process critical to cellular function and human health, but it remains a grand challenge in biophysics. Hydrodynamic interaction (HI) plays a vital role in the self-organization of soft and biological materials, yet its role in protein folding is not fully understood despite folding occurring in a fluid environment. Here, we use the fluid particle dynamics method to investigate many-body hydrodynamic couplings between amino acid residues and fluid motion in the folding kinetics of a coarse-grained four-α-helices bundle protein. Our results reveal that HI helps select fast folding pathways to the native state without being kinetically trapped, significantly speeding up the folding kinetics compared to its absence. First, the directional flow along the protein backbone expedites protein collapse. Then, the incompressibility-induced squeezing flow effects retard the accumulation of non-native hydrophobic contacts, thus preventing the protein from being trapped in local energy minima during the conformational search of the native structure. We also find that the significance of HI in folding kinetics depends on temperature, with a pronounced effect under biologically relevant conditions. Our findings suggest that HI, particularly the short-range squeezing effect, may be crucial in avoiding protein misfolding.

Hydrodynamic Effects on the Collapse Kinetics of Flexible Polyelectrolytes.

Understanding the collapse kinetics of polyelectrolytes (PEs) is crucial for comprehending various biological and industrial phenomena. Despite occurring in an aqueous environment, previous computational studies have overlooked the influence of hydrodynamic interactions (HIs) facilitated by fluid motion. Here, we directly compute the Navier-Stokes equation to investigate the collapse kinetics of a highly charged flexible PE. Our findings reveal that HI accelerates PE collapse induced by hydrophobicity and multivalent salt. In the case of hydrophobicity, HI induces long-range collective motion of monomers, accelerating the coarsening of local clusters through either Brownian-coagulation-like or evaporation-condensation-like processes, depending on the strength of hydrophobicity with respect to electrostatic interaction. Regarding multivalent salt, HI does not affect the condensation dynamics of multivalent ions but facilitates quicker movement of local dipolar clusters along the PE, thereby expediting the collapse process. These results provide valuable insights into the underlying mechanisms of HI in PE collapse kinetics.

Mechanical Slowing Down of Network-Forming Phase Separation of Polymer Solutions.

Phase separation is a fundamental phenomenon leading to spatially heterogeneous material distribution, which is critical in nature, biology, material science, and industry. In ordinary phase separation, the minority phase always forms droplets. Contrary to this common belief, even the minority phase can form a network structure in viscoelastic phase separation (VPS). VPS can occur in any mixture with significant mobility differences between their components and is highly relevant to soft matter and biomatter. In contrast to classical phase separation, experiments have shown that VPS in polymer solutions lacks self-similar coarsening, resulting in the absence of a domain-coarsening scaling law. However, the underlying microscopic mechanism of this behavior remains unknown. To this end, we perform fluid particle dynamics simulations of bead-spring polymers, incorporating many-body hydrodynamic interactions between polymers through a solvent. We discover that polymers in the dense-network-forming phase are stretched and store elastic energy when the deformation speed exceeds the polymer dynamics. This self-generated viscoelastic stress mechanically interferes with phase separation and slows its dynamics, disrupting self-similar growth. We also highlight the essential role of many-body hydrodynamic interactions in VPS. The implications of our findings may hold importance in areas such as biological phase separation, porous material formation, and other fields where network structures play a pivotal role.

Impact of Inverse Squeezing Flow on the Self-Assembly of Oppositely Charged Colloidal Particles under Electric Field.

Hydrodynamic interactions (HIs) play a critical role in the self-organization of colloidal suspensions and biological solutions. However, their roles have remained elusive particularly for charged soft matter systems. Here we consider the role of HIs in the self-assembly of oppositely charged colloidal particles, which is a promising candidate for electrical tunable soft materials. We employ the fluid particle dynamics method to consider many-body HIs and the coupling between the colloid, ion, and fluid motions. We find that, under a constant electric field, oppositely charged colloidal particles form clusters and percolate into a gel network, unlike bundlelike aggregates aligned in the field direction observed by Brownian dynamics simulations neglecting HIs. We reveal that the cluster-forming tendency originates from the incompressibility-induced "inverse squeezing flow" effect that dramatically slows down the disaggregation of attached colloids. Our findings indicate that the HI selects a unique kinetic pathway to the nonequilibrium colloidal self-assembly.

Impact of Charge Regulation on Self-Assembly of Zwitterionic Nanoparticles.

Zwitterionic modification of colloids with weak acids and bases represents a promising strategy in creating functional materials with tunable properties and modeling the self-organization of charged proteins. However, accurate incorporation of the dynamic dissociation or association of ionization groups known as charge regulation (CR) is often intractable in theoretical and computational investigations since charge redistribution and configuration need to be evolved self-consistently. Using hybrid Monte Carlo and molecular dynamics simulations, we demonstrate that a dilute suspension of overall charge-neutral zwitterionic Janus nanoparticles shows a conformational transition from an open assembly of string or bundle to compact cluster along with the variation in pH. The behavior under CR is qualitatively different from the commonly employed constant charge condition where the transition is absent. The CR-induced clustering is due to the inhomogeneous and fluctuating charges localized near the equatorial boundary of the Janus particle. These features are enhanced particularly at low salt concentration and high electrostatic coupling strength. Our results indicate the critical role of charge regulation in the spatial self-organization of zwitterionic nanoparticles.

Accelerated simulation method for charge regulation effects.

The net charge of solvated entities, ranging from polyelectrolytes and biomolecules to charged nanoparticles and membranes, depends on the local dissociation equilibrium of individual ionizable groups. Incorporation of this phenomenon, charge regulation (CR), in theoretical and computational models requires dynamic, configuration-dependent recalculation of surface charges and is therefore typically approximated by assuming constant net charge on particles. Various computational methods exist that address this. We present an alternative, particularly efficient CR Monte Carlo method (CR-MC), which explicitly models the redistribution of individual charges and accurately samples the correct grand-canonical charge distribution. In addition, we provide an open-source implementation in the large-scale Atomic/Molecular Massively Parallel Simulator molecular dynamics (MD) simulation package, resulting in a hybrid MD/CR-MC simulation method. This implementation is designed to handle a wide range of implicit-solvent systems that model discreet ionizable groups or surface sites. The computational cost of the method scales linearly with the number of ionizable groups, thereby allowing accurate simulations of systems containing thousands of individual ionizable sites. By matter of illustration, we use the CR-MC method to quantify the effects of CR on the nature of the polyelectrolyte coil-globule transition and on the effective interaction between oppositely charged nanoparticles.

Conformation and Ionization Behavior of Charge-Regulating Polyelectrolyte Brushes in a Poor Solvent.

Understanding the structural response of weak polyelectrolyte brushes upon external stimuli is crucial for their applications ranging from modifying surface properties to the development of smart and intelligent materials. In this work, coarse-grained molecular dynamics simulations were carried out to investigate the conformation and ionization behavior of charge-regulating polyelectrolyte brushes under poor solvent conditions, using an implicit solvent model. The results show that, while the thickness of a sparse polyelectrolyte brush shows a similar behavior to that of a single chain, namely, a monotonic change as a function of solvent quality (modeled by an effective segment-segment attraction strength parameter) and solution pH, a dense polyelectrolyte brush exhibits more complex behavior. An unexpected reexpansion is observed when the effective segment-segment attraction strength is further increased, especially in the case of a high pH. In the latter case, strong attraction in polymer segments promotes the formation of large, interchain, cylindrical aggregates, leading to an increase in brush thickness.

Particle-particle particle-mesh algorithm for electrolytes between charged dielectric interfaces.

Ionic distributions near charged interfaces control processes from colloidal aggregation to solvent flow in nanodevices. Such interfaces are often characterized by a jump in the permittivity, which gives rise to the surface polarization charge. This induced charge may significantly affect the ionic distributions so that efficient methods for modeling spatially varying dielectrics are needed. We formulate a method with O(N⁡log⁡N) scaling for electrolytes between charged planar interfaces with asymmetric dielectric contrasts. Our approach, which builds on earlier work, is based on combining image charges with the particle-particle particle-mesh algorithm and representing uniform surface charges via an electric field. This enables simulations of complex dielectric interactions that outperform most alternative methods in speed and accuracy. To make the method practically useful, we provide guidelines-based upon careful tests-for choosing optimal simulation parameters. Explicit expressions for the electrostatic forces are given to facilitate the implementation of our algorithm in standard molecular dynamics packages.

Tumor microenvironment characterization in triple-negative breast cancer identifies prognostic gene signature.

We aimed to elucidate the landscape of tumor microenvironment (TME) in triple-negative breast cancer (TNBC). Cohorts from Gene Expression Omnibus database (N = 107) and METABRIC (N = 299) were used as the training set and validation set, respectively. TME was evaluated via single-sample gene set enrichment analysis, and unsupervised clustering was used for cluster identification. Consequently, TNBC was classified into two distinct TME clusters (Cluster 1 and Cluster 2) according to predefined immune-related terms. Cluster 1 was characterized by low immune infiltration with poor prognosis; whereas, Cluster 2 was characterized by high immune infiltration with better survival probability. Further, Cluster 1 had larger tumor volumes, while Cluster 2 had smaller tumor volumes. Finally, a TME signature for prognosis stratification in TNBC was developed and validated. In summary, we comprehensively evaluated the TME of TNBC and constructed a TME signature that correlated with prognosis. Our results provide new insights for the immunotherapy of TNBC.

Harmonic surface mapping algorithm for molecular dynamics simulations of particle systems with planar dielectric interfaces.

We have developed an accurate and efficient method for molecular dynamics simulations of charged particles confined by planar dielectric interfaces. The algorithm combines the image-charge method for near field with the harmonic surface mapping, which converts the contribution of infinite far-field charges into a finite number of charges on an auxiliary spherical surface. We approximate the electrostatic potential of far-field charges via spherical harmonic expansion and determine the coefficients by fitting the Dirichlet-to-Neumann boundary condition, which only requires the potential within the simulation cell. Instead of performing the direct evaluation of spherical harmonic series expansion, we use Green's second identity to transform the series expansion into a spherical integral, which can be accurately represented by discrete charges on the sphere. Therefore, the fast multipole method can be readily employed to sum over all charges within and on the sphere, achieving truly linear O(N) complexity. Our algorithm can be applied to a broad range of charged complex fluids under dielectric confinement.

Dielectric Effects on Ion Transport in Polyelectrolyte Brushes.

Surface-grafted polyelectrolytes provide a versatile way to create functionalized interfaces and nanochannels with externally controllable properties. Understanding the behavior of ions within the brush-like assemblies is crucial for the further development of these devices. We demonstrate that the ion transport through the brushes is governed by the interplay of electrostatic ion-polymer binding and steric effects, leading to a mobility that depends nonmonotonically on grafting density. However, the ion-polymer binding can be modulated by the dielectric properties of the substrate. As a result, surface polarization suppresses ion mobility near insulating interfaces and enhances it near conducting interfaces, even causing a shift from nonmonotonic to monotonic variation with grafting density.