Monte Carlo Simulations of Complexation between Weak Polyelectrolytes and a Charged Nanoparticle. Influence of Polyelectrolyte Chain Length and Concentration

Unusual Aspects of Charge Regulation in Flexible Weak Polyelectrolytes

This article reviews the state of the art of the studies on charge regulation (CR) effects in flexible weak polyelectrolytes (FWPE). The characteristic of FWPE is the strong coupling of ionization and conformational degrees of freedom. After introducing the necessary fundamental concepts, some unconventional aspects of the the physical chemistry of FWPE are discussed. These aspects are: (i) the extension of statistical mechanics techniques to include ionization equilibria and, in particular, the use of the recently proposed Site Binding-Rotational Isomeric State (SBRIS) model, which allows the calculation of ionization and conformational properties on the same foot; (ii) the recent progresses in the inclusion of proton equilibria in computer simulations; (iii) the possibility of mechanically induced CR in the stretching of FWPE; (iv) the non-trivial adsorption of FWPE on ionized surfaces with the same charge sign as the PE (the so-called "wrong side" of the isoelectric point); (v) the influence of macromolecular crowding on CR.

Adsorption of flexible proteins in the 'wrong side' of the isoelectric point: Casein macropeptide as a model system

We analyze the conditions of the adsorption of a flexible peptide onto a charged substrate in the 'wrong side' of the isoelectric point (WSIP), i.e. when surface and peptide charges have the same sign. As a model system, we focus on the casein macropeptide (CMP), both in the aglycosylated (aCMP) and fully glycosydated (gCMP) forms. We model the substrate as a uniformly charged plane while CMP is treated as a bead-and-spring model including electrostatic interactions, excluded volume effects and acid/base equilibria. Adsorption coverage, aminoacid charges and concentration profiles are computed by means of Monte Carlo simulations at fixed pH and salt concentration. We conclude that for different reasons the CMP can be adsorbed to both positively and negatively charged surfaces in the WSIP. For negatively charged surfaces, WSIP adsorption is due to the patchy distribution of charges: the peptide is attached to the surface by the positively charged end of the chain, while the repulsion of the surface for the negatively charged tail is screened by the small ions of the added salt. This effect increases with salt concentration. Conversely, a positively charged substrate induces strong charge regulation of the peptide: the acidic groups are deprotonated, and the peptide becomes negatively charged. This effect is stronger at low salt concentrations and it is more intense for gCMP than for aCMP, due to the presence of the additional sialic groups in gCMP.

Self-Assembly of Miktoarm Star Polyelectrolytes in Solutions with Various Ionic Strengths

We studied the self-assembly of miktoarm star polyelectrolytes with different numbers of arms in solutions with various ionic strengths using coarse-grained molecular dynamic simulations. Spherical micelles are obtained for star polyelectrolytes with fewer arms, whereas wormlike clusters are obtained for star polyelectrolytes with more arms at a low ionic strength environment, with hydrophilic arms showing a stretched conformation. The number of clusters shows an overall decreasing tendency with increasing the number of arms in star polyelectrolytes due to strong electrostatic coupling between polycations and polyanions. The formation of wormlike clusters follows an overall stepwise pathway with an intermittent association-dissociation process for star polyelectrolytes with weak electrostatic coupling. These computational results can provide relevant physical insights to understand the self-assembly mechanism of star polyelectrolytes in solvents with various ionic strengths and to design star polyelectrolytes with functional groups that can fine-tune self-assembled structures for specific applications.

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.

Inducing pH control over the critical micelle concentration of zwitterionic surfactants via polyacids adsorption: Effect of chain length and structure

HYPOTHESIS

The critical concentration above which micelles form from zwitterionic surfactant solutions and their thermodynamic stability is affected by the interaction with weak Brønsted polyacid chains (A_n) via the formation of charged hydrogen bonds between the latter and anionic moieties.

EXPERIMENTS

The interaction between zwitterionic micelles and polyacids capable of forming hydrogen bonds, and its dependence on the environmental pH and polymer structure, has been studied with constant-pH simulations and a restricted primitive model for all electrolytes.

FINDINGS

At low pH, the formation of polyacid/micelle complexes is witnessed independently of the polymer size or structure, so that the concentration above which micelles form is substantially decreased compared to polyacid-free cases. Upon rising pH, polymer desorption takes place within a narrow range of pH values, its location markedly depending on the size and structure of polyacids, and on the relative disposition between headgroup charged moieties. Thus, the desorption onset for long linear polyacids (A₆₀) interacting with sulphobetainic headgroups is roughly two pH units higher than for six decameric chains (6A₁₀) adsorbed onto micelles bearing phosphorylcholinic headgroups. This effect, together with the preferential desorption of chain ends at intermediate pH, may be exploited for drug delivery purposes or building advanced metamaterials.

On Concept of Hybrid in Colloid Sciences

The concepts hybrid and hybridization are common in many scientific fields, as in the taxonomic parts of botany and zoology, in modern genetic, and in the quantum–mechanical theory of atomic–molecular orbitals, which are of foremost relevance in most aspects of modern chemistry. Years later, scientists applied the concept hybrid to colloids, if the particles’ domains are endowed with functionalities differing each from the other in nature and/or composition. For such denomination to be fully valid, the domains belonging to a given hybrid must be recognizable each from another in terms of some intrinsic features. Thus, the concept applies to particles where a given domain has its own physical state, functionality, or composition. Literature examples in this regard are many. Different domains that are present in hybrid colloids self-organize, self-sustain, and self-help, according to the constraints dictated by kinetic and/or thermodynamic stability rules. Covalent, or non-covalent, bonds ensure the formation of such entities, retaining the properties of a given family, in addition to those of the other, and, sometimes, new ones. The real meaning of this behavior is the same as in zoology; mules are pertinent examples, since they retain some features of their own parents (i.e., horses and donkeys) but also exhibit completely new ones, such as the loss of fertility. In colloid sciences, the concept hybrid refers to composites with cores of a given chemical type and surfaces covered by moieties differing in nature, or physical state. This is the result of a mimicry resembling the ones met in a lot of biological systems and foods, too. Many combinations may occur. Silica nanoparticles on which polymers/biopolymers are surface-bound (irrespective of whether binding is covalent or not) are pertinent examples. Here, efforts are made to render clear the concept, which is at the basis of many applications in the biomedical field, and not only. After a historical background and on some features of the species taking part to the formation of hybrids, we report on selected cases met in modern formulations of mixed, and sometimes multifunctional, colloid entities.

On the distribution of hydrophilic polyelectrolytes and their counterions around zwitterionic micelles: the possible impact on the charge density in solution

Despite their charge neutrality, micelles composed of surfactants with zwitterionic headgroups selectively accumulate anions at their hydrophobic core/solution interphase due to electrostatic interactions if headgroup positive moieties are the innermost. This tendency may be markedly enhanced if polyions substitute simple ions. To investigate this possibility, solutions composed of zwitterionic micelles and hydrophilic polyanions have been investigated with Monte Carlo simulations representing the studied systems via primitive electrolyte models. Structural and energetic properties are obtained to highlight the impact of connecting simple ions into polyions on the interactions between electrolytes and micelles. Despite the latter, polyanions conserve their conformational properties. A marked increase in the concentration of charged species inside the micellar corona is, instead, found when polyions are present independently of their charge sign or the headgroup structure. Thus, polyelectrolytes act as "shuttle" for all charged species, with the potential of increasing reactions rates involving the latter due to mass effects. Besides, results for the polyions/micelles mixing free energy and Helmholtz energy profiles indicate that the critical micelle concentration is impacted minimally by hydrophilic polyelectrolytes, an outcome agreeing with experiments. This finding is entirely due to weak enthalpic effects while mixing hydrophilic polyions and micelles. A strong reduction in the screening of the micelle negative charge, acquired following the adsorption of anions in the corona and due to counterions layering just outside it (the so called "chameleon effect"), is forecasted when polyanions substitute monovalent anions.

Influence of macromolecular crowding on the charge regulation of intrinsically disordered proteins

In this work we study the coupling between ionization and conformational properties of two IDPs, histatin-5 and β-amyloid 42, in the presence of neutral and charged crowders. The latter is modeled to resemble bovine serum albumin (BSA). With this aim, semi-grand canonical Monte Carlo simulations are performed, so that the IDP charge is a dynamic property, undergoing protonation/deprotonation processes. Both ionization properties (global and specific amino acid charge and binding capacitance) and radius of gyration are analyzed in a large range of pH values and salt concentrations. Without crowder agents, the titration curve of histatin-5, a polycation, is salt-dependent while that of β-amyloid 42, a polyampholyte, is almost unaffected. The salt concentration is found to be particularly relevant at pH values where the protein binding capacitance (directly linked with charge fluctuation) is larger. Upon addition of neutral crowders, charge regulation is observed in histatin-5, while for β-amyloid 42 this effect is very small. The main mechanism for charge regulation is found to be the effective increase in the ionic strength due to the excluded volume. In the presence of charged crowders, a significant increase in the charge of both IDPs is observed in almost all the pH range. In this case, the IDP charge is altered not only by the increase in the effective ionic strength but also by its direct electrostatic interaction with the charged crowders.

Role of pKA in Charge Regulation and Conformation of Various Peptide Sequences

Peptides containing amino acids with ionisable side chains represent a typical example of weak ampholytes, that is, molecules with multiple titratable acid and base groups, which generally exhibit charge regulating properties upon changes in pH. Charged groups on an ampholyte interact electrostatically with each other, and their interaction is coupled to conformation of the (macro)molecule, resulting in a complex feedback loop. Their charge-regulating properties are primarily determined by the pKA of individual ionisable side-chains, modulated by electrostatic interactions between the charged groups. The latter is determined by the amino acid sequence in the peptide chain. In our previous work we introduced a simple coarse-grained model of a flexible peptide. We validated it against experiments, demonstrating its ability to quantitatively predict charge on various peptides in a broad range of pH. In the current work, we investigated two types of peptide sequences: diblock and alternating, each of them consisting of an equal number of amino acids with acid and base side-chains. We showed that changing the sequence while keeping the same overall composition has a profound effect on the conformation, whereas it practically does not affect total charge on the peptide. Nevertheless, the sequence significantly affects the charge state of individual groups, showing that the zero net effect on the total charge is a consequence of unexpected cancellation of effects. Furthermore, we investigated how the difference between the pKA of acid and base side chains affects the charge and conformation of the peptide, showing that it is possible to tune the charge-regulating properties by following simple guiding principles based on the pKA and on the amino acid sequence. Our current results provide a theoretical basis for understanding of the complex coupling between the ionisation and conformation in flexible polyampholytes, including synthetic polymers, biomimetic materials and biological molecules, such as intrinsically disordered proteins, whose function can be regulated by changes in the pH.

Lipid-Nucleic Acid Complexes: Physicochemical Aspects and Prospects for Cancer Treatment

Cancer is an extremely complex disease, typically caused by mutations in cancer-critical genes. By delivering therapeutic nucleic acids (NAs) to patients, gene therapy offers the possibility to supplement, repair or silence such faulty genes or to stimulate their immune system to fight the disease. While the challenges of gene therapy for cancer are significant, the latter approach (a type of immunotherapy) starts showing promising results in early-stage clinical trials. One important advantage of NA-based cancer therapies over synthetic drugs and protein treatments is the prospect of a more universal approach to designing therapies. Designing NAs with different sequences, for different targets, can be achieved by using the same technologies. This versatility and scalability of NA drug design and production on demand open the way for more efficient, affordable and personalized cancer treatments in the future. However, the delivery of exogenous therapeutic NAs into the patients’ targeted cells is also challenging. Membrane-type lipids exhibiting permanent or transient cationic character have been shown to associate with NAs (anionic), forming nanosized lipid-NA complexes. These complexes form a wide variety of nanostructures, depending on the global formulation composition and properties of the lipids and NAs. Importantly, these different lipid-NA nanostructures interact with cells via different mechanisms and their therapeutic potential can be optimized to promising levels in vitro. The complexes are also highly customizable in terms of surface charge and functionalization to allow a wide range of targeting and smart-release properties. Most importantly, these synthetic particles offer possibilities for scaling-up and affordability for the population at large. Hence, the versatility and scalability of these particles seem ideal to accommodate the versatility that NA therapies offer. While in vivo efficiency of lipid-NA complexes is still poor in most cases, the advances achieved in the last three decades are significant and very recently a lipid-based gene therapy medicine was approved for the first time (for treatment of hereditary transthyretin amyloidosis). Although the path to achieve efficient NA-delivery in cancer therapy is still long and tenuous, these advances set a new hope for more treatments in the future. In this review, we attempt to cover the most important biophysical and physicochemical aspects of non-viral lipid-based gene therapy formulations, with a perspective on future cancer treatments in mind.

Theoretical Modeling of Chemical Equilibrium in Weak Polyelectrolyte Layers on Curved Nanosystems

Surface functionalization with end-tethered weak polyelectrolytes (PE) is a versatile way to modify and control surface properties, given their ability to alter their degree of charge depending on external cues like pH and salt concentration. Weak PEs find usage in a wide range of applications, from colloidal stabilization, lubrication, adhesion, wetting to biomedical applications such as drug delivery and theranostics applications. They are also ubiquitous in many biological systems. Here, we present an overview of some of the main theoretical methods that we consider key in the field of weak PE at interfaces. Several applications involving engineered nanoparticles, synthetic and biological nanopores, as well as biological macromolecules are discussed to illustrate the salient features of systems involving weak PE near an interface or under (nano)confinement. The key feature is that by confining weak PEs near an interface the degree of charge is different from what would be expected in solution. This is the result of the strong coupling between structural organization of weak PE and its chemical state. The responsiveness of engineered and biological nanomaterials comprising weak PE combined with an adequate level of modeling can provide the keys to a rational design of smart nanosystems.

Critical adsorption of multiple polyelectrolytes onto a nanosphere: splitting the adsorption-desorption transition boundary

Employing extensive Monte Carlo computer simulations, we investigate in detail the properties of multichain adsorption of charged flexible polyelectrolytes (PEs) onto oppositely charged spherical nanoparticles (SNPs). We quantify the conditions of critical adsorption-the phase-separation curve between the adsorbed and desorbed states of the PEs-as a function of the SNP surface-charge density and the concentration of added salt. We study the degree of fluctuations of the PE-SNP electrostatic binding energy, which we use to quantify the emergence of the phase subtransitions, including a series of partially adsorbed PE configurations. We demonstrate how the phase-separation adsorption-desorption boundary shifts and splits into multiple subtransitions at low-salt conditions, thereby generalizing and extending the results for critical adsorption of a single PE onto the SNP. The current findings are relevant for finite concentrations of PEs around the attracting SNP, such as the conditions for PE adsorption onto globular proteins carrying opposite electric charges.

Interface Counterion Localization Induces a Switch between Tight and Loose Configurations of Knotted Weak Polyacid Rings despite Intermonomer Coulomb Repulsions

Stochastic simulations have been used to investigate the conformational behavior of knotted weak polyacid rings as a function of pH. Different from the commonly expected ionization-repulsion-expansion scheme upon increasing pH, theoretical results suggest a nonmonotonic behavior of the gyration radius R_g². Polyelectrolyte recontraction at high ionization is induced by the weakening of Coulomb repulsion due to counterions (CIs) localizing at the interphase between the polymer and solvent, and the more marked it appears, the more complex is the knot topology. Compared with strong polyelectrolytic species of identical ionization, weak polyacids present tighter knots due to their ability to localize neutral monomers inside the tangled part. Increasing the solvent Bjerrum length enhances CIs localization, lowering the pH at which polyacids start decreasing their average size. A similar effect is also obtained by increasing the amount of "localizable" cations by adding salts.

Monte Carlo study of the effects of macroion charge distribution on the ionization and adsorption of weak polyelectrolytes and concurrent counterion release

HYPOTHESIS

Adsorption of weak polyelectrolytes onto charged nanoparticles, and concurrent effects such as spatial partitioning of ions may be influenced by details of the polyelectrolyte structure (linear or star-like) and size, by the mobility of the nanoparticle surface charge, or the valence of the nanoparticle counterions.

EXPERIMENTS

Ionization and complexation of weak polyelectrolytes on spherical macroions with monovalent and divalent countrions has been studied with constant-pH Monte Carlo titrations and primitive electrolyte models for linear and star-like polymers capable, also, of forming charged hydrogen bonds. Nanoparticles surface charge has been represented either as a single colloid-centered total charge (CCTC) or as surface-tethered mobile monovalent spherical charges (SMMSC).

FINDINGS

Differences in the average number of adsorbed polyelectrolyte arms and their average charge, and in the relative amount of macroion counterions (m-CI's) released upon polymer adsorption are found between CCTC and SMMSC nanoparticles. The amount of the counterions released also depends on the polymer structure. As CCTC adsorbs a lower number of star-like species arms, the degree of condensation of polymer counterions (p-CI's) onto the polyelectrolyte is also substantially higher for the CCTC colloid, with a concurrent decrease of the osmotic coefficient values.

Nanoparticle-Hydrogel Composites: From Molecular Interactions to Macroscopic Behavior

Hydrogels are materials used in a variety of applications, ranging from tissue engineering to drug delivery. The incorporation of nanoparticles to yield composite hydrogels has gained substantial momentum over the years since these afford tailor-making and extend material mechanical properties far beyond those achievable through molecular design of the network component. Here, we review different procedures that have been used to integrate nanoparticles into hydrogels; the types of interactions acting between polymers and nanoparticles; and how these underpin the improved mechanical and optical properties of the gels, including the self-healing ability of these composite gels, as well as serving as the basis for future development. In a less explored approach, hydrogels have been used as dispersants of nanomaterials, allowing a larger exposure of the surface of the nanomaterial and thus a better performance in catalytic and sensor applications. Furthermore, the reporting capacity of integrated nanoparticles in hydrogels to assess hydrogel properties, such as equilibrium swelling and elasticity, is highlighted.

Simulations of ionization equilibria in weak polyelectrolyte solutions and gels

This article recapitulates the state of the art regarding simulations of ionization equilibria of weak polyelectrolyte solutions and gels. We start out by reviewing the essential thermodynamics of ionization and show how the weak polyelectrolyte ionization differs from the ionization of simple weak acids and bases. Next, we describe simulation methods for ionization reactions, focusing on two methods: the constant-pH ensemble and the reaction ensemble. After discussing the advantages and limitations of both methods, we review the existing simulation literature. We discuss coarse-grained simulations of weak polyelectrolytes with respect to ionization equilibria, conformational properties, and the effects of salt, both in good and poor solvent conditions. This is followed by a discussion of branched star-like weak polyelectrolytes and weak polyelectrolyte gels. At the end we touch upon the interactions of weak polyelectrolytes with other polymers, surfaces, nanoparticles and proteins. Although proteins are an important class of weak polyelectrolytes, we explicitly exclude simulations of protein ionization equilibria, unless they involve protein-polyelectrolyte interactions. Finally, we try to identify gaps and open problems in the existing simulation literature, and propose challenges for future development.

Explicit Ion Effects on the Charge and Conformation of Weak Polyelectrolytes

The titration behavior of weak polyelectrolytes is of high importance, due to their uses in new technologies including nanofiltration and drug delivery applications. A comprehensive picture of polyelectrolyte titration under relevant conditions is currently lacking, due to the complexity of systems involved in the process. One must contend with the inherent structural and solvation properties of the polymer, the presence of counterions, and local chemical equilibria enforced by background salt concentration and solution acidity. Moreover, for these cases, the systems of interest have locally high concentrations of monomers, induced by polymer connectivity or confinement, and thus deviate from ideal titration behavior. This work furthers knowledge in this limit utilizing hybrid Monte Carlo⁻Molecular Dynamics simulations to investigate the influence of salt concentration, pK a , pH, and counterion valence in determining the coil-to-globule transition of poorly solvated weak polyelectrolytes. We characterize this transition at a range of experimentally relevant salt concentrations and explicitly examine the role multivalent salts play in determining polyelectrolyte ionization behavior and conformations. These simulations serve as an essential starting point in understanding the complexation between weak polyelectrolytes and ion rejection of self-assembled copolymer membranes.

Role of Associative Charging in the Entropy-Energy Balance of Polyelectrolyte Complexes

Polyelectrolytes may be classified into two primary categories (strong and weak) depending on how their charge state responds to the local environment. Both of these find use in many applications, including drug delivery, gene therapy, layer-by-layer films, and fabrication of ion filtration membranes. The mechanism of polyelectrolyte complexation is, however, still not completely understood, though experimental investigations suggest that entropy gain due to release of counterions is the key driving force for strong polyelectrolyte complexation. Here we perform a comprehensive thermodynamic investigation through coarse-grained molecular simulations permitting us to calculate the free energy of complex formation. Importantly, our expanded-ensemble methods permit the explicit separation of energetic and entropic contributions to the free energy. Our investigations indicate that entropic contributions indeed dominate the free energy of complex formation for strong polyelectrolytes, but are less important than energetic contributions when weak electrostatic coupling or weak polyelectrolytes are present. Our results provide a new view of the free energy of polyelectrolyte complex formation driven by polymer association, which should also arise in systems with large charge spacings or bulky counterions, both of which act to weaken ion-polymer binding.

pH-Dependent Polyelectrolyte Bridging of Charged Nanoparticles

Systems comprised of polyelectrolytes and charged nanoparticles are of great technological interest, being common components in formulations among other uses. The colloidal stability of formulations is an important issue, and thus a lot of effort has been made to study the interactions of individual components in these systems. Here, the complexation and adsorption of an annealed (pH-dependent) polyelectrolyte to two spherical nanoparticles has been studied using coarse-grained Monte Carlo simulations. This has been done mainly by varying the solution pH and separation distance (concentration) between the nanoparticles. The polyelectrolyte charge distribution is seen to vary with nanoparticle separation distance, and its ability to bridge both nanoparticles changes with pH. The flexible polyelectrolyte creates compact, multilink bridges at short nanoparticle separation distances and evolves to a stretched single-link bridge at longer distances, where a larger fraction of the polyelectrolyte wraps around the nanoparticles. The annealed polyelectrolyte is also compared with a quenched polyelectrolyte of similar fixed fractional charge. Here, a difference is found in the adsorption ability at low pH/ionization due to the ability of the annealed polyelectrolytes to concentrate charges in the vicinity of the nanoparticle. At intermediate polyelectrolyte charge fractions and with increasing nanoparticle separation distances, the annealed system is able to link nanoparticles at larger distances as compared to the quenched, in good agreement with experimental observations. The results in this work contribute to the understanding of the effect of annealed polyelectrolytes and pH variations in the phase behavior of polyelectrolyte-nanoparticle systems, potentially aiding in the design and optimization of pH-responsive systems.

Hydrogel-Based Drug Delivery Nanosystems for the Treatment of Brain Tumors

Chemotherapy is commonly associated with limited effectiveness and unwanted side effects in normal cells and tissues, due to the lack of specificity of therapeutic agents to cancer cells when systemically administered. In brain tumors, the existence of both physiological barriers that protect tumor cells and complex resistance mechanisms to anticancer drugs are additional obstacles that hamper a successful course of chemotherapy, thus resulting in high treatment failure rates. Several potential surrogate therapies have been developed so far. In this context, hydrogel-based systems incorporating nanostructured drug delivery systems (DDS) and hydrogel nanoparticles, also denoted nanogels, have arisen as a more effective and safer strategy than conventional chemotherapeutic regimens. The former, as a local delivery approach, have the ability to confine the release of anticancer drugs near tumor cells over a long period of time, without compromising healthy cells and tissues. Yet, the latter may be systemically administered and provide both loading and targeting properties in their own framework, thus identifying and efficiently killing tumor cells. Overall, this review focuses on the application of hydrogel matrices containing nanostructured DDS and hydrogel nanoparticles as potential and promising strategies for the treatment and diagnosis of glioblastoma and other types of brain cancer. Some aspects pertaining to computational studies are finally addressed.

Monte Carlo simulations of weak polyelectrolyte microgels: pH-dependence of conformation and ionization

In this study, we investigated the effect of pH on single weak acidic polyelectrolyte microgels under salt-free conditions with (i) varying microgel concentration, (ii) varying content of acidic groups and (iii) different crosslinking densities using Monte Carlo simulations under explicit consideration of the protonation/deprotonation reaction. We assessed both global properties, such as the degree of ionization, the degree of swelling and the counterion distribution, and local properties such as the radial network ionization profile and the ionization along the polymer chains as a function of pH. We found a pronounced suppression of the pH-dependent ionization of the microgels, as compared to the ideal titration behavior and a shift of the titration curve to a higher pH originating in the proximity of acidic groups in the microgel. In contrast to macroscopic gels, counterions can leave the microgel, resulting in an effective charge of the network, which hinders the ionization. A decreasing microgel concentration leads to an increased effective charge of the microgel and a more pronounced shift of the titration curve. The number of acidic groups showed only a weak effect on the ionization behavior of the microgels. For two different microgels with different crosslinking densities, similar scaling of the gel size was observed. A distinct transition from an uncharged and unswollen to a highly charged and expanded polymer network was observed for all investigated microgels. The degree of swelling mainly depends on the degree of ionization. An inhomogeneous distribution of the degree of ionization along the radial profile of the microgel was found.