Molecular dynamics simulations for microscopic behavior of water molecules in the vicinity of zwitterionic self-assembled monolayers

Zwitterionic Biomaterials

The term "zwitterionic polymers" refers to polymers that bear a pair of oppositely charged groups in their repeating units. When these oppositely charged groups are equally distributed at the molecular level, the molecules exhibit an overall neutral charge with a strong hydration effect via ionic solvation. The strong hydration effect constitutes the foundation of a series of exceptional properties of zwitterionic materials, including resistance to protein adsorption, lubrication at interfaces, promotion of protein stabilities, antifreezing in solutions, etc. As a result, zwitterionic materials have drawn great attention in biomedical and engineering applications in recent years. In this review, we give a comprehensive and panoramic overview of zwitterionic materials, covering the fundamentals of hydration and nonfouling behaviors, different types of zwitterionic surfaces and polymers, and their biomedical applications.

Advanced density-based methods for the characterization of materials, binding events, and kinetics

Investigations of the densities of chemicals and materials bring valuable insights into the fundamental understanding of matter and processes. Recently, advanced density-based methods have been developed with wide measurement ranges (i.e. 0-23 g cm^-3), high resolutions (i.e. 10^-6 g cm^-3), compatibility with different types of samples and the requirement of extremely low volumes of sample (as low as a single cell). Certain methods, such as magnetic levitation, are inexpensive, portable and user-friendly. Advanced density-based methods are, therefore, beneficially used to obtain absolute density values, composition of mixtures, characteristics of binding events, and kinetics of chemical and biological processes. Herein, the principles and applications of magnetic levitation, acoustic levitation, electrodynamic balance, aqueous multiphase systems, and suspended microchannel resonators for materials science are discussed.

Interactions between Acrylate/Methacrylate Biomaterials and Organic Foulants Evaluated by Molecular Dynamics Simulations of Simplified Binary Mixtures

Improving hydrophilicity is a key factor for enhancing the biocompatibility of polymer surfaces. Nevertheless, previous studies have reported that poly(2-methoxyethyl acrylate) (PMEA) surfaces demonstrate markedly better biocompatibility than more hydrophilic poly(2-hydroxyethyl methacrylate) (PHEMA) surfaces. In this work, the origins of the excellent biocompatibility of the PMEA surface are investigated using molecular dynamics (MD) simulations of simplified binary mixtures of acrylate/methacrylate trimers and organic solvents, with n-nonane, 1,5-pentanediol, or 1-octanol serving as the probe organic foulants. The interactions between the acrylate/methacrylate trimers and solvent molecules were evaluated by calculating the radial distribution function (RDF), with the resulting curves indicating that the 2-methoxyethyl acrylate (MEA) trimer has a lower affinity for n-nonane molecules than the 2-hydroxyethyl methacrylate (HEMA) trimer. This result agrees with the experimental consensus that the biocompatibility of PMEA surfaces is better than that of PHEMA surfaces, supporting the hypothesis that the affinity between an acrylate/methacrylate trimer and a foulant molecule in a simplified binary mixture is a significant factor in determining a surface's antifouling properties. The RDF curves obtained for the other two solvent systems exhibited behavior that further highlighted the advantages of the PMEA surfaces as biocompatible polymers. In addition, the validity of employing the second virial coefficient (B₂) as a predictor of antifouling properties was explored. The order of the B₂ values of different binary mixtures indicated that the MEA trimers have the lowest affinities with n-nonane molecules, which confirms that although PMEA is more hydrophobic than PHEMA, it exhibits better biocompatibility. This analysis demonstrates that the MEA's weaker miscibility with nonpolar foulants contributes to the excellent biocompatibility of PMEA. Thus, B₂ is a promising criterion for assessing the miscibility between acrylate/methacrylate materials and nonpolar organic foulants, which indicates the potential for predicting the antifouling properties of acrylate/methacrylate polymer materials by evaluating the value of B₂.

Effects of Side-Chain Spacing and Length on Hydration States of Poly(2-methoxyethyl acrylate) Analogues: A Molecular Dynamics Study

Hydration states of polymers are known to directly influence the adsorption of biomolecules. Particularly, intermediate water (IW) has been found able to prevent protein adsorption. Experimental studies have examined the IW content and nonthrombogenicity of poly(2-methoxyethyl acrylate) analogues with different side-chain spacings and lengths, which are HPx (x is the number of backbone carbons in a monomer) and PMCyA (y is the number of carbons in-between ester and ether oxygens of the side-chain) series, respectively. HPx was reported to possess more IW content but lower nonthrombogenicity compared to PMCyA with analogous composition. To understand the reason for the conflict, molecular dynamics simulations were conducted to elucidate the difference in the properties between the HPx and PMCyA. Simulation results showed that the presence of more methylene groups in the side chain more effectively prohibits water penetration in the polymer than those in the polymer backbone, causing a lower IW content in the PMCyA. At a high water content, the methoxy oxygen in the shorter side chain of the HPx cannot effectively bind water compared to that in the PMCyA side chain. HPx side chains may have more room to contact with molecules other than water (e.g., proteins), causing experimentally less nonthrombogenicity of HPx than that of PMCyA. In summary, theoretical simulations successfully explained the difference in the effects of side-chain spacing and length in atomistic scale.

Molecular Dynamics Study on the Water Mobility and Side-Chain Flexibility of Hydrated Poly(ω-methoxyalkyl acrylate)s

Intermediate water (IW) is known to play an important role in the antifouling property of biocompatible polymers. However, how IW prevents protein adsorption is still unclear. To understand the role of IW in the antifouling mechanism, molecular dynamics simulation was used to investigate the dynamic properties of water and side-chains for hydrated poly(ω-methoxyalkyl acrylate)s (PMCxA, where x indicates the number of methylene carbons) with x = 1-6 and poly(n-butyl acrylate) (PBA) in this study. Since the polymers uptake more water than their equilibrium water content (EWC) at the polymer/water interface, we analyzed the hydrated polymers at a water content higher than that of EWC. It was found that the water molecules interacting with one polymer oxygen atom (BW1), of which most are IW molecules, in PMC2A exhibit the lowest mobility, while those in PBA and PMC1A show a higher mobility. The result was consistent with the expectation that the biocompatible polymer with a long-resident hydration layer possesses good antifouling property. Through the detailed analysis of side-chain binding with three different types of BW1 molecules, we found that the amount of side-chains simultaneously interacting with two BW1 molecules, which exhibit the highest flexibility among the three kinds of side-chains, is the lowest for PMC1A. The high mobility of BW1 is thus suggested as the main factor for the poor protein adsorption resistance of PMC1A even though it possesses enough IW content and relatively flexible side-chains. Contrarily, a maximum amount of side-chains simultaneously interacting with two BW1 molecules was found in the hydrated PMC3A. The moderate side-chain length of PMC3A allows side-chains to simultaneously interact with two BW1 molecules and minimizes the hydrophobic part attractively interacting with a protein at the polymer/water interface. The unique structure of PMC3A may be the reason causing the best protein adsorption resistance among the PMCxAs.

Molecular Understanding of Ion Effect on Polyzwitterion Conformation in an Aqueous Environment

Polyzwitterions (PZs) are promising materials for the antifouling in reverse osmosis and nanofiltration membrane technology for water treatment. Fundamental understanding of the structure and molecular interactions involving zwitterions is crucial to the optimal design of antifouling in membrane separation. Here we employ the umbrella sampling and molecular dynamics simulations to investigate molecular interactions between sulfobetaine/carboxybetaine zwitterions and different metal ions (Na⁺, K⁺, and Ca²⁺) in an aqueous solution. The simulation results show that these ions can form stable or metastable contact ionic/solvent-shared-ionic pairs with zwitterions. Simulations at different grafting densities of PZ brush arrays reveal complex competitive association mechanisms, which are attributed to nonbonded electrostatic and van der Waals interactions among zwitterions, water molecules, and different metal ions in an aqueous environment. While the high-grafting density of the PZ brush array leads to a strong branch association between different zwitterions in water, this association is decreased at intermediate- and low-grafting densities due to strong zwitterion-water interactions. More importantly, adding ions into water at intermediate- and low-grafting densities further breaks down the zwitterion branch association, resulting in a randomly oriented and dispersed branch configuration with significant swelling of the polymers. The degree of swelling depends on the type of ions, which further changes the surface electrostatic potential of PZ coatings.

Interaction Mechanisms of Zwitterions with Opposite Dipoles in Aqueous Solutions

Zwitterionic groups have been widely used in antibiofouling surfaces to resist nonspecific adsorption of proteins and other biomolecules. The interactions among zwitterionic groups have attracted considerable attention in bioengineering, whereas the understanding of their nanomechanical mechanism still remains limited. In this work, the interaction mechanisms between two zwitterionic groups with opposite dipoles, i.e., phosphorylcholine (PC) and sulfobetaine (SB), have been investigated via direct force measurements using an atomic force microscope (AFM) and dynamic adsorption tests using the quartz crystal microbalance with dissipation monitoring technique (QCM-D) in aqueous solutions. The AFM force measurements show that the adhesive forces between contacted zwitterionic surfaces during separation in both symmetric and asymmetric configurations were close, mainly due to the enforced alignment of opposing dipole pairs via complementary orientations under confinement. The solution salinity and pH had almost negligible influence on the adhesion measured during surface separation. The QCM-D adsorption tests of PC-headed lipid on PC and SB surfaces showed some degree of adsorption of lipid molecules on the SB surface, whereas not on the PC surface. The different adsorption behaviors indicate that because the outermost negatively charged sulfonic group on the SB faced the aqueous solution, this configuration could facilitate it to form an attractive electrostatic interaction with the PC head of lipid molecules in the solution. This work shows that in addition to hydration and steric interactions, the zwitterionic dipole-induced interactions play an important role in the adhesion and antifouling behaviors of the zwitterionic molecules and surfaces. The improved fundamental understanding provides useful insights into the development of new functional materials and coatings with antifouling applications.

Hydration States of Cholinium Phosphate-Type Ionic Liquids as a Function of Water Content

We investigated the hydration states of cholinium phosphate-type ionic liquids (ILs) in relation to ion structure, focusing on the influence of the hydroxyl group of the cation and the alkyl chain length of the anion. Water activity measurements provided information on the macroscopic hydration states of the hydrated ILs, while NMR measurements and molecular dynamics simulations clearly showed the microscopic interactions and coordination of the water molecules. The hydrogen bonding networks in these ILs were influenced by the anion structure and water content, and the mobility of water molecules was influenced by the number of hydroxyl groups in the cation and anion.

Molecular Simulations of the Hydration Behavior of a Zwitterion Brush Array and Its Antifouling Property in an Aqueous Environment

We carried out umbrella sampling and molecular dynamics (MD) simulations to investigate molecular interactions between sulfobetaine zwitterions or between sulfobetaine brushes in different media. Simulation results show that it is more energetically favorable for the two sulfobetaine zwitterions or brushes to be fully hydrated in aqueous solutions than in vacuum where strong ion pairs are formed. Structural properties of the hydrated sulfobetaine brush array and its antifouling behavior against a foulant gel are subsequently studied through steered MD simulations. We find that sulfobetaine brush arrays with different grafting densities have different structures and antifouling mechanisms. At a comparably higher grafting density, the sulfobetaine brush array exhibits a more organized structure which can hold a tightly bound hydration water layer at the interface. Compression of this hydration layer results in a strong repulsive force. However, at a comparably lower grafting density, the brush array exhibits a randomly oriented structure in which the antifouling of the brush array is through the deformation of the sulfobetaine branches.