Multiscale simulations of polyzwitterions in aqueous bulk solutions and brush array configurations

Photo-responsive anti-fouling polyzwitterionic brushes: a mesoscopic simulation

The antifouling effects of a toothbrush-shaped photo-responsive polyzwitterionic membrane were studied via dissipative particle dynamics simulations in this work. The results reveal that the membrane modified by spiropyran methacrylate brushes displays photo-switchable and antifouling capability due to the photo-induced ring-opening reaction. Namely, surface morphology and hydrophilicity change in response to visible or UV light irradiation, which can be observed visually by protein adsorption and desorption. Further study indicates that: (1) brush-modification density can influence the structure and properties of the membrane. With low modification density, systems cannot establish an intact selective layer, which hinders the antifouling ability; as the modification density increases, the intact selective layer can be formed, which is conducive to the expression of photo-responsiveness and antifouling capability. (2) Factors of toothbrush-hair length and grafting ratio can influence the establishment of a light-responsive surface: as the grafting ratio and toothbrush-hair length increase, the light-responsive surface is gradually formed, meanwhile, the antifouling ability can be continuously reinforced under UV light irradiation. (3) As the brushes switch into a zwitterionic merocyanine state under UV exposure, the selective layer swelling becomes stronger than that with a hydrophobic spiropyran state under visible exposure. This is owing to the enhanced interaction between zwitterionic brushes and water, which is the root of the antifouling effect. The present work is expected to provide some guidelines for the design and development of novel antifouling membrane surfaces.

Learning Force Field Parameters from Differentiable Particle-Field Molecular Dynamics

We develop ∂-HylleraasMD (∂-HyMD), a fully end-to-end differentiable molecular dynamics software based on the Hamiltonian hybrid particle-field formalism, and use it to establish a protocol for automated optimization of force field parameters. ∂-HyMD is templated on the recently released HylleraaasMD software, while using the JAX autodiff framework as the main engine for the differentiable dynamics. ∂-HyMD exploits an embarrassingly parallel optimization algorithm by spawning independent simulations, whose trajectories are simultaneously processed by reverse mode automatic differentiation to calculate the gradient of the loss function, which is in turn used for iterative optimization of the force-field parameters. We show that parallel organization facilitates the convergence of the minimization procedure, avoiding the known memory and numerical stability issues of differentiable molecular dynamics approaches. We showcase the effectiveness of our implementation by producing a library of force field parameters for standard phospholipids, with either zwitterionic or anionic heads and with saturated or unsaturated tails. Compared to the all-atom reference, the force field obtained by ∂-HyMD yields better density profiles than the parameters derived from previously utilized gradient-free optimization procedures. Moreover, ∂-HyMD models can predict with good accuracy properties not included in the learning objective, such as lateral pressure profiles, and are transferable to other systems, including triglycerides.

Current status, challenges and prospects of antifouling materials for oncology applications

Targeted therapy has become crucial to modern translational science, offering a remedy to conventional drug delivery challenges. Conventional drug delivery systems encountered challenges related to solubility, prolonged release, and inadequate drug penetration at the target region, such as a tumor. Several formulations, such as liposomes, polymers, and dendrimers, have been successful in advancing to clinical trials with the goal of improving the drug's pharmacokinetics and biodistribution. Various stealth coatings, including hydrophilic polymers such as PEG, chitosan, and polyacrylamides, can form a protective layer over nanoparticles, preventing aggregation, opsonization, and immune system detection. As a result, they are classified under the Generally Recognized as Safe (GRAS) category. Serum, a biological sample, has a complex composition. Non-specific adsorption of chemicals onto an electrode can lead to fouling, impacting the sensitivity and accuracy of focused diagnostics and therapies. Various anti-fouling materials and procedures have been developed to minimize the impact of fouling on specific diagnoses and therapies, leading to significant advancements in recent decades. This study provides a detailed analysis of current methodologies using surface modifications that leverage the antifouling properties of polymers, peptides, proteins, and cell membranes for advanced targeted diagnostics and therapy in cancer treatment. In conclusion, we examine the significant obstacles encountered by present technologies and the possible avenues for future study and development.

Zwitterion-Modified Nanogel Responding to Temperature and Ionic Strength: A Dissipative Particle Dynamics Simulation

The self-assembly and stimuli-responsive properties of nanogel poly(n-isopropylacrylamide) (p(NIPAm)) and zwitterion-modified nanogel poly(n-isopropylacrylamide-co-sulfobetainemethacrylate) (p(NIPAm-co-SBMA)) were explored by dissipative particle dynamics simulations. Simulation results reveal that for both types of nanogel, it is beneficial to form spherical nanogels at polymer concentrations of 5-10%. When the chain length (L) elongates from 10 to 40, the sizes of the nanogels enlarge. As for the p(NIPAm) nanogel, it shows thermoresponsiveness; when it switches to the hydrophilic state, the nanogel swells, and vice versa. The zwitterion-modified nanogel p(NIPAm-co-SBMA) possesses thermoresponsiveness and ionic strength responsiveness concurrently. At 293 K, both hydrophilic p(NIPAm) and superhydrophilic polysulfobetaine methacrylate (pSBMA) could appear on the outer surface of the nanogel; however, at 318 K, superhydrophilic pSBMA is on the outer surface to cover the hydrophobic p(NIPAm) core. As the temperature rises, the nanogel shrinks and remains antifouling all through. The salt-responsive property can be reflected by the nanogel size; the volumes of the nanogels in saline systems are larger than those in salt-free systems as the ionic condition inhibits the shrinkage of the zwitterionic pSBMA. This work exhibits the temperature-responsive and salt-responsive behavior of zwitterion-modified-pNIPAm nanogels at the molecular level and provides guidance in antifouling nanogel design.

Leaflet Tensions Control the Spatio-Temporal Remodeling of Lipid Bilayers and Nanovesicles

Biological and biomimetic membranes are based on lipid bilayers, which consist of two monolayers or leaflets. To avoid bilayer edges, which form when the hydrophobic core of such a bilayer is exposed to the surrounding aqueous solution, a single bilayer closes up into a unilamellar vesicle, thereby separating an interior from an exterior aqueous compartment. Synthetic nanovesicles with a size below 100 nanometers, traditionally called small unilamellar vesicles, have emerged as potent platforms for the delivery of drugs and vaccines. Cellular nanovesicles of a similar size are released from almost every type of living cell. The nanovesicle morphology has been studied by electron microscopy methods but these methods are limited to a single snapshot of each vesicle. Here, we review recent results of molecular dynamics simulations, by which one can monitor and elucidate the spatio-temporal remodeling of individual bilayers and nanovesicles. We emphasize the new concept of leaflet tensions, which control the bilayers' stability and instability, the transition rates of lipid flip-flops between the two leaflets, the shape transformations of nanovesicles, the engulfment and endocytosis of condensate droplets and rigid nanoparticles, as well as nanovesicle adhesion and fusion. To actually compute the leaflet tensions, one has to determine the bilayer's midsurface, which represents the average position of the interface between the two leaflets. Two particularly useful methods to determine this midsurface are based on the density profile of the hydrophobic lipid chains and on the molecular volumes.

Globular Proteins and Where to Find Them within a Polymer Brush-A Case Study

Protein adsorption by polymerized surfaces is an interdisciplinary topic that has been approached in many ways, leading to a plethora of theoretical, numerical and experimental insight. There is a wide variety of models trying to accurately capture the essence of adsorption and its effect on the conformations of proteins and polymers. However, atomistic simulations are case-specific and computationally demanding. Here, we explore universal aspects of the dynamics of protein adsorption through a coarse-grained (CG) model, that allows us to explore the effects of various design parameters. To this end, we adopt the hydrophobic-polar (HP) model for proteins, place them uniformly at the upper bound of a CG polymer brush whose multibead-spring chains are tethered to a solid implicit wall. We find that the most crucial factor affecting the adsorption efficiency appears to be the polymer grafting density, while the size of the protein and its hydrophobicity ratio come also into play. We discuss the roles of ligands and attractive tethering surfaces to the primary adsorption as well as secondary and ternary adsorption in the presence of attractive (towards the hydrophilic part of the protein) beads along varying spots of the backbone of the polymer chains. The percentage and rate of adsorption, density profiles and the shapes of the proteins, alongside with the respective potential of mean force are recorded to compare the various scenarios during protein adsorption.

Prediction of zwitterion hydration and ion association properties using machine learning

Molecular dynamics simulations were performed to study the hydration and ion association properties of a library of zwitterionic molecules with varying charged moieties and spacer chemistries in pure water and with Na⁺ and Cl^- ions. The structure and dynamics of associations were calculated using the radial distribution and residence time correlation function. Resulting association properties are used as target variables for a machine learning model, with cheminformatic descriptors of the molecule subunits used as descriptors. Prediction of hydration properties revealed that steric and hydrogen bonding descriptors were of greatest importance and there was influence from the cationic moiety on the anionic moiety hydration properties. Ion association properties prediction performed poorly, which is attributed to the role of hydration layers in ion association dynamics. This study is the first to quantitatively describe the influence of subunit chemistry on hydration and ion association properties of zwitterions. These quantitative descriptions supplement prior studies of zwitterion association and previously described design principles.

Conformations and dynamic behaviors of confined wormlike chains in a pressure-driven flow

The conformations and dynamic behaviors of wormlike chains confined by a slit in a pressure-driven flow were investigated using dissipative particle dynamics method. The wormlike chains exhibit varying conformations due to the varying shear stresses across the slit. The wormlike chain solution can be well described by the power-law fluid, and the power-law index decreases with the increase in chain rigidity. We also presented that the wormlike chain undergoes tumbling motion in the vicinity of the wall in the presence of pressure-driven flow. We also found that the wormlike chains can migrate both away from the wall and slightly away from the slit center, and the migration away from the slit center increases as the chain rigidity is increased because of hydrodynamic interactions induced in a more rigid wormlike chain.

Modeling Polyzwitterion-Based Drug Delivery Platforms: A Perspective of the Current State-of-the-Art and Beyond

Drug delivery platforms are anticipated to have biocompatible and bioinert surfaces. PEGylation of drug carriers is the most approved method since it improves water solubility and colloid stability and decreases the drug vehicles' interactions with blood components. Although this approach extends their biocompatibility, biorecognition mechanisms prevent them from biodistribution and thus efficient drug transfer. Recent studies have shown (poly)zwitterions to be alternatives for PEG with superior biocompatibility. (Poly)zwitterions are super hydrophilic, mainly stimuli-responsive, easy to functionalize and they display an extremely low protein adsorption and long biodistribution time. These unique characteristics make them already promising candidates as drug delivery carriers. Furthermore, since they have highly dense charged groups with opposite signs, (poly)zwitterions are intensely hydrated under physiological conditions. This exceptional hydration potential makes them ideal for the design of therapeutic vehicles with antifouling capability, i.e., preventing undesired sorption of biologics from the human body in the drug delivery vehicle. Therefore, (poly)zwitterionic materials have been broadly applied in stimuli-responsive "intelligent" drug delivery systems as well as tumor-targeting carriers because of their excellent biocompatibility, low cytotoxicity, insignificant immunogenicity, high stability, and long circulation time. To tailor (poly)zwitterionic drug vehicles, an interpretation of the structural and stimuli-responsive behavior of this type of polymer is essential. To this end, a direct study of molecular-level interactions, orientations, configurations, and physicochemical properties of (poly)zwitterions is required, which can be achieved via molecular modeling, which has become an influential tool for discovering new materials and understanding diverse material phenomena. As the essential bridge between science and engineering, molecular simulations enable the fundamental understanding of the encapsulation and release behavior of intelligent drug-loaded (poly)zwitterion nanoparticles and can help us to systematically design their next generations. When combined with experiments, modeling can make quantitative predictions. This perspective article aims to illustrate key recent developments in (poly)zwitterion-based drug delivery systems. We summarize how to use predictive multiscale molecular modeling techniques to successfully boost the development of intelligent multifunctional (poly)zwitterions-based systems.