Softness and non-spherical shape define the phase behavior and the structural properties of lysozyme in aqueous solutions

Attenuation of pressure dips underneath piles of spherocylinders

The discrete element method (DEM) was used to simulate the piling of rod-like (elongated sphero-cylindrical) particles, mainly focusing on the effect of particle shape on the structural and force properties of the piles. In this work, rod-like particles of different aspect ratios were discharged on a flat surface to form wedge-shaped piles. The surface properties of the piles were characterized in terms of angle of repose and stress at the bottom of the piles. The results showed that the rise of the angle of repose became slower with the increase of particle aspect ratio. The pressure dip underneath the piles reached the maximum when the particle aspect ratio was around 1.6, beyond which the pressure dip phenomenon became attenuated. Both the pressure dip and the shear stress dip were quantitatively examined. The structure and forces inside the piles were further analyzed to understand the change in pressure dip, indicating that "bridging" or "arching" structures within the piles were the cause of the pressure dip.

Cluster Crystals Stabilized by Hydrophobic and Electrostatic Interactions

Cluster crystals are crystalline materials in which each site is occupied by multiple identical particles, atoms, colloids, or polymers. There are two classes of systems that make cluster crystals. One is composed of particles that interact via potentials that are bound at the origin and thus are able to penetrate each other. The other consists of non-interpenetrating particles whose interaction potential diverges at the origin. The goal of this work is to find which systems of the second class can make cluster crystals that are stable at room temperature. First, the general properties of the required potentials are established using an analytical model and Monte Carlo simulations. Next, we ask how such potentials can be constructed by combining hydrophobic attraction and electrostatic repulsion. A colloid model with a hard-sphere core and a repulsive wall is introduced to mimic the hydrophobic interaction. Charge is added to create long-range repulsion. A search in the parameter space of the colloid size, counterion type, and charge configuration uncovers several models for which effective colloid-colloid interaction, determined in explicit solvent as a potential of mean force, has the necessary shape. For the effective potential, cluster crystals are confirmed as low free-energy configurations in replica-exchange molecular dynamics simulations, which also generate the respective transition temperature. The model that exhibits a transition above room temperature is further studied in explicit solvent. Simulations on a 10 ns time scale show that crystalline conformations are stable below the target temperature but disintegrate rapidly above it, supporting the idea that hydrophobic and electrostatic interactions are sufficient to induce an assembly of cluster crystals. Finally, we discuss which physical systems are good candidates for experimental observations of cluster crystals.

Explicit-water theory for the salt-specific effects and Hofmeister series in protein solutions

Effects of addition of salts on stability of aqueous protein solutions are studied theoretically and the results are compared with experimental data. In our approach, all the interacting species, proteins, ions, and water molecules, are accounted for explicitly. Water molecules are modeled as hard spheres with four off-center attractive square-well sites. These sites serve to bind either another water or to solvate the ions or protein charges. The ions are represented as charged hard spheres, and decorated by attractive sites to allow solvation. Spherical proteins simultaneously possess positive and negative groups, represented by charged hard spheres, attached to the surface of the protein. The attractive square-well sites, mimicking the protein-protein van der Waals interaction, are located on the surface of the protein. To obtain numerical results, we utilized the energy route of Wertheim's associative mean spherical approximation. From measurable properties, we choose to calculate the second virial coefficient B2, which is closely related to the tendency of proteins to aggregate and eventually crystalize. Calculations are in agreement with experimental trends: (i) For low concentration of added salt, the alkali halide salts follow the inverse Hofmeister series. (ii) At higher concentration of added salt, the trend is reversed. (iii) When cations are varied, the salts follow the direct Hofmeister series. (iv) In contrast to the colloidal theories, our approach correctly predicts the non-monotonic behavior of B2 upon addition of salts. (v) With respect to anions, the theory predicts for the B2 values to follow different sequences below and above the iso-ionic point, as also confirmed experimentally. (vi) A semi-quantitative agreement between measured and calculated values for the second virial coefficient, as functions of pH of solution and added salt type and concentration, is obtained.

Protein folding, binding, and droplet formation in cell-like conditions

The many bystander macromolecules in the crowded cellular environments present both steric repulsion and weak attraction to proteins undergoing folding or binding and hence impact the thermodynamic and kinetic properties of these processes. The weak but nonrandom binding with bystander macromolecules may facilitate subcellular localization and biological function. Weak binding also leads to the emergence of a protein-rich droplet phase, which has been implicated in regulating a variety of cellular functions. All these important problems can now be addressed by realistic modeling of intermolecular interactions. Configurational sampling of concentrated protein solutions is an ongoing challenge.