Combining Scattering Experiments and Colloid Theory to Characterize Charge Effects in Concentrated Antibody Solutions.

Charges and their contribution to protein-protein interactions are essential for the key structural and dynamic properties of monoclonal antibody (mAb) solutions. In fact, they influence the apparent molecular weight, the static structure factor, the collective diffusion coefficient, or the relative viscosity, and their concentration dependence. Further, charges play an important role in the colloidal stability of mAbs. There exist standard experimental tools to characterize mAb net charges, such as the measurement of the electrophoretic mobility, the second virial coefficient, or the diffusion interaction parameter. However, the resulting values are difficult to directly relate to the actual overall net charge of the antibody and to theoretical predictions based on its known molecular structure. Here, we report the results of a systematic investigation of the solution properties of a charged IgG1 mAb as a function of concentration and ionic strength using a combination of electrophoretic measurements, static and dynamic light scattering, small-angle X-ray scattering, and tracer particle-based microrheology. We analyze and interpret the experimental results using established colloid theory and coarse-grained computer simulations. We discuss the potential and limits of colloidal models for the description of the interaction effects of charged mAbs, in particular pointing out the importance of incorporating shape and charge anisotropy when attempting to predict structural and dynamic solution properties at high concentrations.

Softness matters: effects of compression on the behavior of adsorbed microgels at interfaces.

Deformable colloids and macromolecules adsorb at interfaces as they decrease the interfacial energy between the two media. The deformability, or softness, of these particles plays a pivotal role in the properties of the interface. In this study, we employ a comprehensive in situ approach, combining neutron reflectometry with molecular dynamics simulations, to thoroughly examine the profound influence of softness on the structure of microgel Langmuir monolayers under compression. Lateral compression of both hard and soft microgel particle monolayers induces substantial structural alterations, leading to an amplified protrusion of the microgels into the aqueous phase. However, a critical distinction emerges: hard microgels are pushed away from the interface, in stark contrast to the soft ones, which remain firmly anchored to it. Concurrently, on the air-exposed side of the monolayer, lateral compression induces a flattening of the surface of the hard monolayer. This phenomenon is not observed for the soft particles as the monolayer is already extremely flat even in the absence of compression. These findings significantly advance our understanding of the key role of softness on both the equilibrium phase behavior of the monolayer and its effect when soft colloids are used as stabilizers of responsive interfaces and emulsions.

Two-step deswelling in the Volume Phase Transition of thermoresponsive microgels.

Thermoresponsive microgels are one of the most investigated types of soft colloids, thanks to their ability to undergo a Volume Phase Transition (VPT) close to ambient temperature. However, this fundamental phenomenon still lacks a detailed microscopic understanding, particularly regarding the presence and the role of charges in the deswelling process. This is particularly important for the widely used poly(N-isopropylacrylamide)-based microgels, where the constituent monomers are neutral but charged groups arise due to the initiator molecules used in the synthesis. Here, we address this point combining experiments with state-of-the-art simulations to show that the microgel collapse does not happen in a homogeneous fashion, but through a two-step mechanism, entirely attributable to electrostatic effects. The signature of this phenomenon is the emergence of a minimum in the ratio between gyration and hydrodynamic radii at the VPT. Thanks to simulations of microgels with different cross-linker concentrations, charge contents, and charge distributions, we provide evidence that peripheral charges arising from the synthesis are responsible for this behavior and we further build a universal master curve able to predict the two-step deswelling. Our results have direct relevance on fundamental soft condensed matter science and on applications where microgels are involved, ranging from materials to biomedical technologies.

Interfacial Fluid Rheology of Soft Particles.

In situ interfacial rheology and numerical simulations are used to investigate microgel monolayers in a wide range of packing fractions, ζ_{2D}. The heterogeneous particle compressibility determines two flow regimes characterized by distinct master curves. To mimic the microgel architecture and reproduce experiments, an interaction potential combining a soft shoulder with the Hertzian model is introduced. In contrast to bulk conditions, the elastic moduli vary nonmonotonically with ζ_{2D} at the interface, confirming long-sought predictions of reentrant behavior for Hertzian-like systems.

Using Cluster Theory to Calculate the Experimental Structure Factors of Antibody Solutions.

Monoclonal antibody solutions are set to become a major therapeutic tool in the years to come, capable of targeting various diseases by clever design of their antigen binding site. However, the formulation of stable solutions suitable for patient self-administration typically presents challenges, as a result of the increase in viscosity that often occurs at high concentrations. Here, we establish a link between the microscopic molecular details and the resulting properties of an antibody solution through the characterization of clusters, which arise in the presence of self-associating antibodies. In particular, we find that experimental small-angle X-ray scattering data can be interpreted by means of analytical models previously exploited for the study of polymeric and colloidal objects, based on the presence of such clusters. The latter are determined by theoretical calculations and supported by computer simulations of a coarse-grained minimal model, in which antibodies are treated as Y-shaped colloidal molecules and attractive domains are designed as patches. Using the theoretically predicted cluster size distributions, we are able to describe the experimental structure factors over a wide range of concentration and salt conditions. We thus provide microscopic evidence for the well-established fact that the concentration-dependent increase in viscosity is originated by the presence of clusters. Our findings bring new insights on the self-assembly of monoclonal antibodies, which can be exploited for guiding the formulation of stable and effective antibody solutions.

Correction: Concentration and temperature dependent interactions and state diagram of dispersions of copolymer microgels.

Correction for 'Concentration and temperature dependent interactions and state diagram of dispersions of copolymer microgels' by José Ruiz-Franco et al., Soft Matter, 2023, 19, 3614-3628, https://doi.org/10.1039/D3SM00120B.

Concentration and temperature dependent interactions and state diagram of dispersions of copolymer microgels.

We investigate by means of small angle neutron scattering experiments and numerical simulations the interactions and inter-particle arrangements of concentrated dispersions of copolymer poly(N-isopropylacrylamide)-poly(ethylene glycol methyl ether methacrylate) (PNIPAM-PEGMA) microgels across the volume phase transition (VPT). The scattering data of moderately concentrated dispersions are accurately modeled at all temperatures by using a star polymer form factor and static structure factors calculated from the effective potential obtained from simulations. Interestingly, for temperatures below the VPT temperature (VPTT), the radius of gyration and blob size of the particles significantly decrease with increasing the effective packing fraction in the non-overlapping regime. This is attributed to the presence of charges in the system associated with the use of an ionic initiator in the synthesis. Simulations using the experimentally corroborated interaction potential are used to explore the state diagram in a wide range of effective packing fractions. Below and slightly above the VPTT, the system undergoes an arrest transition mainly driven by the soft repulsion between the particles. Only well above the VPTT the system is found to phase separate before arresting. Our results highlight the versatility and potential of copolymer PNIPAM-PEGMA microgels to explore different kinds of arrested states balancing attraction and repulsion by changing temperature and packing fraction.

Auxetic polymer networks: The role of crosslinking, density, and disorder.

Low-crosslinked polymer networks have recently been found to behave auxetically when subjected to small tensions, that is, their Poisson's ratio ν becomes negative. In addition, for specific state points, numerical simulations revealed that diamond-like networks reach the limit of mechanical stability, exhibiting values of ν = -1, a condition that we define as hyper-auxeticity. This behavior is interesting per se for its consequences in materials science but is also appealing for fundamental physics because the mechanical instability is accompanied by evidence of criticality. In this work, we deepen our understanding of this phenomenon by performing a large set of equilibrium and stress-strain simulations in combination with phenomenological elasticity theory. The two approaches are found to be in good agreement, confirming the above results. We also extend our investigations to disordered polymer networks and find that the hyper-auxetic behavior also holds in this case, still manifesting a similar critical-like behavior as in the diamond one. Finally, we highlight the role of the number density, which is found to be a relevant control parameter determining the elastic properties of the system. The validity of the results under disordered conditions paves the way for an experimental investigation of this phenomenon in real systems, such as hydrogels.

Structure and elasticity of model disordered, polydisperse, and defect-free polymer networks.

The elasticity of disordered and polydisperse polymer networks is a fundamental problem of soft matter physics that is still open. Here, we self-assemble polymer networks via simulations of a mixture of bivalent and tri- or tetravalent patchy particles, which result in an exponential strand length distribution analogous to that of experimental randomly cross-linked systems. After assembly, the network connectivity and topology are frozen and the resulting system is characterized. We find that the fractal structure of the network depends on the number density at which the assembly has been carried out, but that systems with the same mean valence and same assembly density have the same structural properties. Moreover, we compute the long-time limit of the mean-squared displacement, also known as the (squared) localization length, of the cross-links and of the middle monomers of the strands, showing that the dynamics of long strands is well described by the tube model. Finally, we find a relation connecting these two localization lengths at high density and connect the cross-link localization length to the shear modulus of the system.

Universality class of the motility-induced critical point in large scale off-lattice simulations of active particles.

We perform large-scale computer simulations of an off-lattice two-dimensional model of active particles undergoing a motility-induced phase separation (MIPS) to investigate the system's critical behaviour close to the critical point of the MIPS curve. By sampling steady-state configurations for large system sizes and performing finite size scaling analysis we provide exhaustive evidence that the critical behaviour of this active system belongs to the Ising universality class. In addition to the scaling observables that are also typical of passive systems, we study the critical behaviour of the kinetic temperature difference between the two active phases. This quantity, which is always zero in equilibrium, displays instead a critical behavior in the active system which is well described by the same exponent of the order parameter in agreement with mean-field theory.

Molecular insights on poly(N-isopropylacrylamide) coil-to-globule transition induced by pressure.

By using extensive all-atom molecular dynamics simulations of an atactic linear polymer chain, we provide microscopic insights into poly(N-isopropylacrylamide) (PNIPAM) coil-to-globule transition addressing the roles played by both temperature and pressure. We detect a coil-to-globule transition up to large pressures, showing a reentrant behavior of the critical temperature with increasing pressure in agreement with experimental observations. Furthermore, again confirming the experimental findings, we report the existence at high pressures of a new kind of globular state. It is characterized by a more structured hydration shell that is closer to PNIPAM hydrophobic domains, as compared to the globular state observed at atmospheric pressure. Our results highlight that temperature and pressure induce a PNIPAM coil-to-globule transition through different molecular mechanisms, opening the way for a systematic use of both thermodynamic variables to tune the location of the transition and the properties of the associated swollen/collapsed states.

Effective potentials induced by mixtures of patchy and hard co-solutes.

The addition of co-solutes to colloidal suspensions is often employed to induce tunable depletion interactions. In this work, we investigate effective colloidal interactions arising from binary co-solute mixtures of hard spheres and patchy particles. By changing the relative concentration of the two species, we show that the resulting effective potential Veff continuously changes from the one obtained for a single-component hard sphere co-solute to that mediated by the single-component patchy particle co-solute. Interestingly, we find that, independent of the relative concentration of the two components, the resulting Veff is additive, i.e., it is well-described by the linear combination of the effective interactions mediated by respective pure co-solutes. However, a breakdown of the additivity occurs when the co-solute mixture is close to the onset of a demixing transition. These results represent a step forward in understanding and predicting colloidal behavior in complex and crowded environments and for exploiting this knowledge to design targeted colloidal superstructures.

Dynamical properties of different models of elastic polymer rings: Confirming the link between deformation and fragility.

We report extensive numerical simulations of different models of 2D polymer rings with internal elasticity. We monitor the dynamical behavior of the rings as a function of the packing fraction to address the effects of particle deformation on the collective response of the system. In particular, we compare three different models: (i) a recently investigated model [N. Gnan and E. Zaccarelli, Nat. Phys. 15, 683 (2019)] where an inner Hertzian field providing the internal elasticity acts on the monomers of the ring, (ii) the same model where the effect of such a field on the center of mass is balanced by opposite forces, and (iii) a semi-flexible model where an angular potential between adjacent monomers induces strong particle deformations. By analyzing the dynamics of the three models, we find that in all cases, there exists a direct link between the system fragility and particle asphericity. Among the three, only the first model displays anomalous dynamics in the form of a super-diffusive behavior of the mean-squared displacement and of a compressed exponential relaxation of the density auto-correlation function. We show that this is due to the combination of internal elasticity and the out-of-equilibrium force self-generated by each ring, both of which are necessary ingredients to induce such a peculiar behavior often observed in experiments of colloidal gels. These findings reinforce the role of particle deformation, connected to internal elasticity, in driving the dynamical response of dense soft particles.

Glass and Jamming Rheology in Soft Particles Made of PNIPAM and Polyacrylic Acid.

The phase behaviour of soft colloids has attracted great attention due to the large variety of new phenomenologies emerging from their ability to pack at very high volume fractions. Here we report rheological measurements on interpenetrated polymer network microgels composed of poly(N-isopropylacrylamide) (PNIPAM) and polyacrylic acid (PAAc) at fixed PAAc content as a function of weight concentration. We found three different rheological regimes characteristic of three different states: a Newtonian shear-thinning fluid, an attractive glass characterized by a yield stress, and a jamming state. We discuss the possible molecular mechanisms driving the formation of these states.

Crowding in the Eye Lens: Modeling the Multisubunit Protein ß-Crystallin with a Colloidal Approach.

We present a multiscale characterization of aqueous solutions of the bovine eye lens protein ßH crystallin from dilute conditions up to dynamical arrest, combining dynamic light scattering, small-angle x-ray scattering, tracer-based microrheology, and neutron spin echo spectroscopy. We obtain a comprehensive explanation of the observed experimental signatures from a model of polydisperse hard spheres with additional weak attraction. In particular, the model predictions quantitatively describe the multiscale dynamical results from microscopic nanometer cage diffusion over mesoscopic micrometer gradient diffusion up to macroscopic viscosity. Based on a comparative discussion with results from other crystallin proteins, we suggest an interesting common pathway for dynamical arrest in all crystallin proteins, with potential implications for the understanding of crowding effects in the eye lens.

A Colloid Approach to Self-Assembling Antibodies.

Concentrated solutions of monoclonal antibodies have attracted considerable attention due to their importance in pharmaceutical formulations; yet, their tendency to aggregate and the resulting high viscosity pose considerable problems. Here we tackle this problem by a soft condensed matter physics approach, which combines a variety of experimental measurements with a patchy colloid model, amenable of analytical solution. We thus report results of structural antibodies and dynamic properties obtained through scattering methods and microrheological experiments. We model the data using a colloid-inspired approach, explicitly taking into account both the anisotropic shape of the molecule and its charge distribution. Our simple patchy model is able to disentangle self-assembly and intermolecular interactions and to quantitatively describe the concentration-dependence of the osmotic compressibility, collective diffusion coefficient, and zero shear viscosity. Our results offer new insights on the key problem of antibody formulations, providing a theoretical and experimental framework for a quantitative assessment of the effects of additional excipients or chemical modifications and a prediction of the resulting viscosity.

Numerical modelling of non-ionic microgels: an overview.

Microgels are complex macromolecules. These colloid-sized polymer networks possess internal degrees of freedom and, depending on the polymer(s) they are made of, can acquire a responsiveness to variations of the environment (temperature, pH, salt concentration, etc.). Besides being valuable for many practical applications, microgels are also extremely important to tackle fundamental physics problems. As a result, these last years have seen a rapid development of protocols for the synthesis of microgels, and more and more research has been devoted to the investigation of their bulk properties. However, from a numerical standpoint the picture is more fragmented, as the inherently multi-scale nature of microgels, whose bulk behaviour crucially depends on the microscopic details, cannot be handled at a single level of coarse-graining. Here we present an overview of the methods and models that have been proposed to describe non-ionic microgels at different length-scales, from the atomistic to the single-particle level. We especially focus on monomer-resolved models, as these have the right level of details to capture the most important properties of microgels, responsiveness and softness. We suggest that these microscopic descriptions, if realistic enough, can be employed as starting points to develop the more coarse-grained representations required to investigate the behaviour of bulk suspensions.

Numerical insights on ionic microgels: structure and swelling behaviour.

Recent progress has been made in the numerical modelling of neutral microgel particles with a realistic, disordered structure. In this work we extend this approach to the case of co-polymerised microgels where a thermoresponsive polymer is mixed with acidic groups. We compare the cases where counterions directly interact with microgel charges or are modelled implicitly through a Debye-Hückel description. We do so by performing extensive numerical simulations of single microgels across the volume phase transition (VPT) varying the temperature and the fraction of charged monomers. We find that the presence of charges considerably alters the microgel structure, quantified by the monomer density profiles and by the form factors of the microgels, particularly close to the VPT. We observe significant deviations between the implicit and explicit models, with the latter comparing more favourably to available experiments. In particular, we observe a shift of the VPT temperature to larger values as the amount of charged monomers increases. We also find that below the VPT the microgel-counterion complex is almost neutral, while it develops a net charge above the VPT. Interestingly, under these conditions the collapsed microgel still retains a large amount of counterions inside its structure. Since these interesting features cannot be captured by the implicit model, our results show that it is crucial to explicitly include the counterions in order to realistically model ionic thermoresponsive microgels.

Rheological investigation of gels formed by competing interactions: A numerical study.

A transition from solid-like to liquid-like behavior occurs when colloidal gels are subjected to a prolonged exposure to a steady shear. This phenomenon, which is characterized by a yielding point, is found to be strongly dependent on the packing fraction. However, it is not yet known how the effective inter-particle potential affects this transition. To this aim, we present a numerical investigation of the rheology of equilibrium gels in which a short-range depletion is complemented by a long-range electrostatic interaction. We observe a single yielding event in the stress-strain curve, occurring at a fixed strain. The stress overshoot is found to follow a power-law dependence on the Péclet number, with an exponent larger than that found in depletion gels, suggesting that its value may depend systematically on the underlying colloid-colloid interactions. We also establish a mapping between equilibrium states and steady states under shear, which allows us to identify the structural modifications induced by the presence of the shear. Remarkably, we find that steady states corresponding to the same Péclet number, obtained by different combinations of shear rate and solvent viscosity, show identical structural and rheological properties. Our results highlight the importance of understanding the coupling between colloidal interactions, solvent effects, and flow to be able to describe the microscopic organization of colloidal particles under shear.