Clustering in complex fluids

Molecular Dynamics Investigation of Giant Clustering in Small-Molecule Solutions: The Case of Aqueous PEHA

The importance of the formation of giant clusters in solution, in nature and industry, is increasingly recognized. However, relatively little attention has been paid to the formation of giant clusters in solutions of small, relatively soluble but nonamphiphilic molecules. In this work, we present a general methodology based on molecular dynamics that can be used to investigate such systems. As a case study, we focus on the formation of apparently stable clusters of pentaethylenehexamine (PEHA) in water. These clusters have been used as templates for the construction of bioinspired silica nanoparticles. To better understand clustering in this system, we study the effect of PEHA protonation state (neutral, +1, and +2) and counterion type (chloride or acetate) on PEHA clustering in dilute aqueous solutions (200 and 400 mM) using large-scale classical molecular dynamics. We find that large stable clusters are formed by singly charged PEHA with chloride or acetate as the counterion, although it is not clear for the case with acetate whether bulk phase separation, that might lead to precipitation, would eventually occur. Large clusters also appear to be stable for doubly charged PEHA with acetate, the less soluble counterion. We attribute this behavior to a form of complex coacervation, observed here for relatively small and highly soluble molecules (PEHA + counterion) rather than the large polyions usually found to form such coacervates. We discuss whether this behavior might also be described by an effective SALR (short-range attraction, long-range repulsion) interaction. This work might help future studies of additives for the design of novel bioinspired templated nanomaterials and of giant clustering in small-molecule solutions more generally.

Molecular Dynamics Investigation of Clustering in Aqueous Glycine Solutions

Recent experiments with undersaturated aqueous glycine solutions have repeatedly exhibited the presence of giant liquid-like clusters or nanodroplets around 100 nm in diameter. These nanodroplets re-appear even after careful efforts for their removal and purification of the glycine solution. The composition of these clusters is not clear, although it has been suggested that they are mainly composed of glycine, a small and very soluble amino acid. To gain insights into this phenomenon, we study the aggregation of glycine in aqueous solutions at concentrations below the experimental solubility limit using large-scale molecular dynamics simulations under ambient conditions. Three protonation states of glycine (zwitterion = GLZ, anion = GLA, and cation = GLC) are simulated using molecular force fields based on the 1.14*CM1A partial charge scheme, which incorporates the OPLS all-atom force field and TIP3P water. When initiated from dispersed states, we find that giant clusters do not form in our simulations unless salt impurities are present. Moreover, if simulations are initiated from giant cluster states, we find that they tend to dissolve in the absence of salt impurities. Therefore, the simulation results provide little support for the possibility that the giant clusters seen in experiments are composed purely of glycine (and water). Considering that strenuous efforts are made in experiments to remove impurities such as salt, we propose that the giant clusters observed might instead result from the aggregation of reaction products of aqueous glycine, such as diketopiperazine or other oligoglycines which may be difficult to separate from glycine using conventional methods, or their co-aggregation with glycine.

Equilibrium mechanisms of self-limiting assembly

Self-assembly is a ubiquitous process in synthetic and biological systems, broadly defined as the spontaneous organization of multiple subunits (e.g. macromolecules, particles) into ordered multi-unit structures. The vast majority of equilibrium assembly processes give rise to two states: one consisting of dispersed disassociated subunits, and the other, a bulk-condensed state of unlimited size. This review focuses on the more specialized class of self-limiting assembly, which describes equilibrium assembly processes resulting in finite-size structures. These systems pose a generic and basic question, how do thermodynamic processes involving non-covalent interactions between identical subunits "measure" and select the size of assembled structures? In this review, we begin with an introduction to the basic statistical mechanical framework for assembly thermodynamics, and use this to highlight the key physical ingredients that ensure equilibrium assembly will terminate at finite dimensions. Then, we introduce examples of self-limiting assembly systems, and classify them within this framework based on two broad categories: self-closing assemblies and open-boundary assemblies. These include well-known cases in biology and synthetic soft matter - micellization of amphiphiles and shell/tubule formation of tapered subunits - as well as less widely known classes of assemblies, such as short-range attractive/long-range repulsive systems and geometrically-frustrated assemblies. For each of these self-limiting mechanisms, we describe the physical mechanisms that select equilibrium assembly size, as well as potential limitations of finite-size selection. Finally, we discuss alternative mechanisms for finite-size assemblies, and draw contrasts with the size-control that these can achieve relative to self-limitation in equilibrium, single-species assemblies.

Structure and kinetics in colloidal films with competing interactions

Using computer simulation we explore how two-dimensional systems of colloids with a combination of short-range attractive and long-range repulsive interactions generate complex structures and kinetics. Cooperative effects mean the attractive potential, despite being very short-ranged compared to the repulsion, can have significant effects on large-scale structure. By considering the number of particles occupying a notional "repulsion zone" defined by the repulsion length scale, we classify different characteristic structural regimes in which the combination of attraction and repulsion leads to different structural and kinetic outcomes, such as compact clustering, chain labyrinths, and coexisting clusters and chains. In some regimes small changes in repulsion range and/or area fraction can change timescales of structural evolution by many orders of magnitude.

Nanosecond Tracer Diffusion as a Probe of the Solution Structure and Molecular Mobility of Protein Assemblies: The Case of Ovalbumin

Protein diffusion is not only an important process ensuring biological function but can also be used as a probe to obtain information on structural properties of protein assemblies in liquid solutions. Here, we explore the oligomerization state of ovalbumin at high protein concentrations by means of its short-time self-diffusion. We employ high-resolution incoherent quasielastic neutron scattering to access the self-diffusion on nanosecond timescales, on which interparticle contacts are not altered. Our results indicate that ovalbumin in aqueous (D₂O) solutions occurs in increasingly large assemblies of its monomeric subunits with rising protein concentration. It changes from nearly monomeric toward dimeric and ultimately larger than tetrameric complexes. Simultaneously, we access information on the internal molecular mobility of ovalbumin on the nanometer length scale and compare it with results obtained for bovine serum albumin, immunoglobulin, and β-lactoglobulin.

The extended law of corresponding states when attractions meet repulsions

Short-range attractive colloids show well-defined phase behaviour in the absence of repulsions, and highly intriguing equilibrium gelation in the presence of long-range repulsions. We present the state diagram of short-range attractive colloids with repulsions that range from fully screened to intermediately ranged, i.e. longer-ranged than the attractions, but shorter ranged than the colloid size. We demonstrate that although the macroscopic phase behaviour does not change perceptibly, there is a dramatic increase of inhomogeneities once the repulsions become longer-ranged than the attractions. The interaction potentials are characterized with small angle neutron scattering, and used to renormalize the state diagram with the minimum in the interaction potential, min[U(r)], and with the reduced second virial coefficient, B2*. We find that the extended law of corresponding states captures the onset of phase separation for shorter ranged repulsions, but fails for longer ranged repulsions. Instead, for a given model of U(r), the transition from visually homogeneous fluid to phase separation and/or gelation can be rescaled with min[U(r)] over the full range of repulsions. Finally, we suggest a generic state diagram to describe the effect of repulsions on short-range attractive systems.

Crowding-Controlled Cluster Size in Concentrated Aqueous Protein Solutions: Structure, Self- and Collective Diffusion

We investigate the concentration-controlled formation of clusters in β-lactoglobulin (BLG) protein solutions combining structural and dynamical scattering techniques. The static structure factor from small-angle X-ray scattering as well as de-Gennes narrowing in the nanosecond diffusion function D(q) from neutron spin echo spectroscopy support a picture of cluster formation. Using neutron backscattering spectroscopy, a monotonous increase of the average hydrodynamic cluster radius is monitored over a broad protein concentration range, corresponding to oligomeric structures of BLG ranging from the native dimers up to roughly four dimers. The results suggest that BLG forms compact clusters that are static on the observation time scale of several nanoseconds. The presented analysis provides a general framework to access the structure and dynamics of macromolecular assemblies in solution.

Structural transitions and hysteresis in clump- and stripe-forming systems under dynamic compression

Using numerical simulations, we study the dynamical evolution of particles interacting via competing long-range repulsion and short-range attraction in two dimensions. The particles are compressed using a time-dependent quasi-one dimensional trough potential that controls the local density, causing the system to undergo a series of structural phase transitions from a low density clump lattice to stripes, voids, and a high density uniform state. The compression proceeds via slow elastic motion that is interrupted with avalanche-like bursts of activity as the system collapses to progressively higher densities via plastic rearrangements. The plastic events vary in magnitude from small rearrangements of particles, including the formation of quadrupole-like defects, to large-scale vorticity and structural phase transitions. In the dense uniform phase, the system compresses through row reduction transitions mediated by a disorder-order process. We characterize the rearrangement events by measuring changes in the potential energy, the fraction of sixfold coordinated particles, the local density, and the velocity distribution. At high confinements, we find power law scaling of the velocity distribution during row reduction transitions. We observe hysteresis under a reversal of the compression when relatively few plastic rearrangements occur. The decompressing system exhibits distinct phase morphologies, and the phase transitions occur at lower compression forces as the system expands compared to when it is compressed.

Role of clusters in nonclassical nucleation and growth of protein crystals

The development of multistep nucleation theory has spurred on experimentalists to find intermediate metastable states that are relevant to the solidification pathway of the molecule under interest. A great deal of studies focused on characterizing the so-called "precritical clusters" that may arise in the precipitation process. However, in macromolecular systems, the role that these clusters might play in the nucleation process and in the second stage of the precipitation process, i.e., growth, remains to a great extent unknown. Therefore, using biological macromolecules as a model system, we have studied the mesoscopic intermediate, the solid end state, and the relationship that exists between them. We present experimental evidence that these clusters are liquid-like and stable with respect to the parent liquid and metastable compared with the emerging crystalline phase. The presence of these clusters in the bulk liquid is associated with a nonclassical mechanism of crystal growth and can trigger a self-purifying cascade of impurity-poisoned crystal surfaces. These observations demonstrate that there exists a nontrivial connection between the growth of the macroscopic crystalline phase and the mesoscopic intermediate which should not be ignored. On the other hand, our experimental data also show that clusters existing in protein solutions can significantly increase the nucleation rate and therefore play a relevant role in the nucleation process.