Charge regulation mechanism in end-tethered weak polyampholytes

A sustainable bioremediation of vanadium from marine environment and value-addition using potential thraustochytrids

Rising concerns about global environmental degradation underscore the pressing need for effective solutions to combat heavy metal pollution. Industries such as semiconductor and steel production discharge vanadium into marine ecosystems, posing significant risks to both marine life and human health. The current study investigates efficacy of utilizing marine thraustochytrid for efficient vanadium removal outcompeting other microbial sources. By optimizing pH and temperature conditions during harvesting, achieved a remarkable 50.80 % enhancement in vanadium removal efficiency, from 19.31 to 29.12 mg/L. Furthermore, chelating agents EDTA and citric acid supplementation demonstrated promising enhancements, reaching up to 31.21 and 32.59 mg/L, respectively. Notably, vanadium-treated biomass supplemented with citric acid exhibited maximum enhancement in lipid content, from 58.47 to 75.34 %, indicating thraustochytrid's potential for biofuel production. This study presents a sustainable approach for industrial-scale vanadium bioremediation, aligning with Sustainable Development Goals focused on dual benefits of environmental protection and renewable energy.

On the Validity of Constant pH Simulations

Constant pH (cpH) simulations are now a standard tool for investigating charge regulation in coarse-grained models of polyelectrolytes and colloidal systems. Originally developed for studying solutions with implicit ions, extending this method to systems with explicit ions or solvents presents several challenges. Ensuring proper charge neutrality within the simulation cell requires performing titration moves in sync with the insertion or deletion of ions, a crucial aspect often overlooked in the literature. Contrary to the prevailing views, cpH simulations are inherently grand-canonical, meaning that the controlled pH is that of the reservoir. The presence of the Donnan potential between the implicit reservoir and the simulation cell introduces significant differences between titration curves calculated for open and closed systems; the pH of an isolated (closed) system is different from the pH of the reservoir for the same protonation state of the polyelectrolyte. To underscore this point, in this paper, we will compare the titration curves calculated using the usual cpH algorithm with those from the exact canonical simulation algorithm. In the latter case, titration moves adhere to the correct detailed balance condition, and pH is calculated using the recently introduced surface Widom insertion algorithm. Our findings reveal a very significant difference between the titration isotherms obtained using the standard cpH algorithm and the canonical titration algorithm, emphasizing the importance of using the correct simulation approach when studying charge regulation of polyelectrolytes, proteins, and colloidal particles.

Induced phase transformation in ionizable colloidal nanoparticles

Acid-base equilibria directly influence the functionality and behavior of particles in a system. Due to the ionizing effects of acid-base functional groups, particles will undergo charge exchange. The degree of ionization and their intermolecular and electrostatic interactions are controlled by varying the pH and salt concentration of the solution in a system. Although the pH can be tuned in experiments, it is hard to model this effect using simulations or theoretical approaches. This is due to the difficulty in treating charge regulation and capturing the cooperative effects in a colloidal suspension with Coulombic interaction. In this work, we analyze a suspension of ionizable colloidal particles via molecular dynamics (MD) simulations, along with Monte Carlo simulations for charge regulation (MC-CR) and derive a phase diagram of the system as a function of pH. It is observed that as pH increases, particles functionalized with acid groups change their arrangement from face-centered cubic (FCC) packing to a disordered state. We attribute these transitions to an increase in the degree of charge polydispersity arising from an increase in pH. Our work shows that charge regulation leads to amorphous solids in colloids when the mean nanoparticle charge is sufficiently high.

Aggregation-Induced Asymmetric Charge States of Amino Acids in Supramolecular Nanofibers

Electrostatic interactions contribute critically to the kinetic pathways and thermodynamic outcomes of peptide self-assembly involving one or more than one charged amino acids. While it is well understood in protein folding that those amino acids with acidic/basic side chains could shift their pKas when placed in a hydrophobic microenvironment, to what extent aggregation of monomeric peptide units from the bulk solution could alter their charged status and how this change in pKa values would reciprocally impact their assembly outcomes. Here, we design and analyze two solution systems containing peptide amphiphiles with hydrocarbon chains of different lengths to determine the factor of deprotonation on assembly. Our results suggest that models of supramolecular nanofibers with uniformly distributed, fully charged amino acids are oversimplified. We demonstrate, with molecular dynamics simulations, and validate with experimental results that asymmetric, different protonation states of the peptides lead to distinct nanostructures after self-assembly. The results give estimates on the electrostatic interactions in peptide amphiphiles required for their self-assembly and shed light on modeling molecular assembly systems containing charged amino acids.

A new method for reactive constant pH simulations

We present a simulation method that allows us to calculate the titration curves for systems undergoing protonation/deprotonation reactions-such as charged colloidal suspensions with acidic/basic surface groups, polyelectrolytes, polyampholytes, and proteins. The new approach allows us to simultaneously obtain titration curves both for systems in contact with salt and acid reservoir (semi-grand canonical ensemble) and for isolated suspensions (canonical ensemble). To treat the electrostatic interactions, we present a new method based on Ewald summation-which accounts for the existence of both Bethe and Donnan potentials within the simulation cell. We show that the Donnan potential dramatically affects the pH of a suspension. Counterintuitively, we find that in suspensions with a large volume fraction of nanoparticles and low ionic strength, the number of deprotonated groups can be 100% larger in an isolated system, compared to a system connected to a reservoir by a semi-permeable membrane-both systems being at exactly the same pH.

Tailoring Interactions of Random Copolymer Polyelectrolyte Complexes to Remove Nanoplastic Contaminants from Water

We investigate the usage of polyelectrolyte complex materials for water remediation purposes, specifically their ability to remove nanoplastics from water, on which there is currently little to no prior research. We demonstrate that oppositely charged random copolymers are effective at quantitatively removing nanoplastic contamination from aqueous solution. The mechanisms underlying this remediation ability are explored through computational simulations, with corroborating quartz crystal microbalance adsorption experiments. We find that hydrophobic nanostructures and interactions likely play an important role.

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.

Reactive Monte Carlo simulations for charge regulation of colloidal particles

We use a reactive Monte Carlo simulation method and the primitive model of electrolyte to study acid-base equilibrium that controls charge regulation in colloidal systems. The simulations are performed in a semi-grand canonical ensemble in which colloidal suspension is in contact with a reservoir of salt and strong acid. The interior of colloidal particles is modeled as a low dielectric medium, different from the surrounding water. The effective colloidal charge is calculated for different numbers of surface acidic groups, pH, salt concentrations, and types of electrolyte. In the case of potassium chloride, the titration curves are compared with the experimental measurements obtained using potentiometric titration. A good agreement is found between simulations and experiments. In the case of lithium chloride, the specific ionic adsorption is taken into account through the partial dehydration of lithium ion.

Transfer Matrix Model of pH Effects in Polymeric Complex Coacervation

Oppositely charged polyelectrolytes can undergo an associative phase separation, in a process known as polymeric complex coacervation. This phenomenon is driven by the electrostatic attraction between polyanion and polycation species, leading to the formation of a polymer-dense coacervate phase and a coexisting polymer-dilute supernatant phase. There has been significant recent interest in the physical origin and features of coacervation; yet notably, experiments often use weak polyelectrolytes the charge state of which depends on solution pH, while theoretical or computational efforts typically assume strong polyelectrolytes that remain fully charged. There have been only a few efforts to address this limitation, and thus there has been little exploration of how pH can affect complex coacervation. In this paper, we modify a transfer matrix theory of coacervation to account for acid-base equilibria, taking advantage of its ability to directly account for some local ion correlations that will affect monomer charging. We show that coacervation can stabilize the charged state of a weak polyelectrolyte via the proximity of oppositely charged monomers, and can lead to asymmetric phase diagrams where the positively and negatively charged polyelectrolytes exhibit different behaviors near the pK_a of either chain. Specifically, there is a partitioning of one of the salt species to a coacervate to maintain electroneutrality when one of the polyelectrolytes is only partially charged. This results in the depletion of the same salt species in the supernatant, and overall can suppress phase separation. We also demonstrate that, when one of the species is only partially charged, mixtures that are off-stoichiometric in volume fraction but stoichiometric in charge exhibit the greatest propensity to form coacervate phases.