Enhanced protein crystallization around the metastable critical point

We report on a computer-simulation study of homogeneous crystal nucleation in a model for globular proteins. We find that the presence of a metastable vapour-liquid critical point drastically changes the pathway for the formation of a critical nucleus. But what is more important, the large density fluctuations near the critical point also lowers the free-energy barrier to nucleation and hence increases the nucleation rate. As␣the location of the vapour-liquid critical point can be controlled by changing the solvent conditions, our simulation results suggest a guided approach to protein crystallization. Received: 4 June 1998/Accepted: 3 September 1998 / Published online: 10 December 1998     

Free-energy landscape of nucleation with an intermediate metastable phase studied using capillarity approximation

Capillarity approximation is used to study the free-energy landscape of nucleation when an intermediate metastable phase exists. The critical nucleus that corresponds to the saddle point of the free-energy landscape as well as the whole free-energy landscape can be studied using this capillarity approximation, and various scenarios of nucleation and growth can be elucidated. In this study, we consider a model in which a stable solid phase nucleates within a metastable vapor phase when an intermediate metastable liquid phase exists. We predict that a composite critical nucleus that consists of a solid core and a liquid wetting layer as well as pure liquid and pure solid critical nuclei can exist depending not only on the supersaturation of the liquid phase relative to that of the vapor phase but also on the wetting behavior of the liquid surrounding the solid. The existence of liquid critical nucleus indicates that the phase transformation from metastable vapor to stable solid occurs via the intermediate metastable liquid phase, which is quite similar to the scenario of nucleation observed in proteins and colloidal systems. By studying the minimum-free-energy path on the free-energy landscape, we can study the evolution of the composition of solid and liquid within nuclei which is not limited to the critical nucleus.     

Steady-state nucleation rate and flux of composite nucleus at saddle point

The steady-state nucleation rate and flux of composite nucleus at the saddle point is studied by extending the theory of binary nucleation. The Fokker-Planck equation that describes the nucleation flux is derived using the Master equation for the growth of the composite nucleus, which consists of the core of the final stable phase surrounded by a wetting layer of the intermediate metastable phase nucleated from a metastable parent phase recently evaluated by Iwamatsu [J. Chem. Phys. 134, 164508 (2011)]. The Fokker-Planck equation is similar to that used in the theory of binary nucleation, but the non-diagonal elements exist in the reaction rate matrix. First, the general solution for the steady-state nucleation rate and the direction of nucleation flux is derived. Next, this information is then used to study the nucleation of composite nucleus at the saddle point. The dependence of steady-state nucleation rate as well as the direction of nucleation flux on the reaction rate in addition to the free-energy surface is studied using a model free-energy surface. The direction of nucleation current deviates from the steepest-descent direction of the free-energy surface. The results show the importance of two reaction rate constants: one from the metastable environment to the intermediate metastable phase and the other from the metastable intermediate phase to the stable new phase. On the other hand, the gradient of the potential Φ or the Kramers crossover function (the commitment or splitting probability) is relatively insensitive to reaction rates or free-energy surface.     

A_a型Patchy胶体粒子的相态结构

根据热力学微扰理论研究了Aa型Patchy胶体的相态结构,考察了Patch之间的缔合强度及Patch数目对体系相态结构的调控机制.利用相平衡原理给出了Patchy胶体的流体、玻璃态固体和面心立方晶体之间转变及溶胶-凝胶转变的相图,讨论了玻璃态和晶态固相的成核机制、临界现象和相变问题.研究结果表明,Patchy胶体粒子之间的缔合作用和Patch数目可以显著地调控体系的三相点、临界温度和临界密度等特征.在高温条件下,Patchy胶体以一次成核方式结晶;而在低温条件下则以两步成核方式逐步成核结晶,中间经过非晶态的玻璃态固相作为过渡.说明Patchy粒子之间的缔合作用对其相态结构具有决定性影响,因而成为调控体系聚集态结构的重要因素.     

Pathways of cluster growth and kinetic slowing down in a model of short-range attractive colloids

We present the results of an extensive 3D Brownian dynamics simulation of the self-assembly of colloidal particles for a short-range attractive model that is quenched below its metastable critical point. In particular, results are obtained in the small-volume-fraction, low-temperature region in which we find so-called sticky beads that diffuse around the system, without reaching a final large cluster on the timescale of our simulation. For larger volume fractions in this low-temperature regime, a gel forms as the result of kinetically slowed down spinodal decomposition, as shown earlier for other short-range attractive models (Foffi, G.; De Michele, C.; Sciortino, F.; Tartaglia, P. Phys. Rev. Lett. 2005, 94, 078301. Zaccarelli, E. J. Phys.: Condens. Matter 2007, 19, 323101). We also show that for quenches below the critical point but above the intersection of the binodal with the glass line, two-step crystallization takes place. For sufficiently small volume fractions, the first step is the nucleation of dense fluid drops, followed by the second step of crystallization within these drops, as first proposed for a model of protein crystallization for quenches just above the metastable critical point (ten Wolde, P. R.; Frenkel, D. Science 1997, 277, 1975). For larger values of the volume fraction, the initial step is spinodal decomposition that leads to the formation of an interconnected network of low- and high-density fluids. The second step is crystallization that takes place within the dense fluid phase.     

Analytic calculation of phase diagrams for solutions containing colloids or globular proteins

Second-order Barker–Henderson perturbation theory gives phase diagrams for colloid and protein solutions that include stable and metastable fluid–fluid, solid–fluid, and solid–solid phases. The potential of mean force is described by a hard-sphere interacting with a Yukawa potential. Calculations for different ranges of attraction show that, as expected, fluid–fluid coexistence becomes metastable when the potential becomes short-ranged. For a very short-ranged Yukawa potential, the phase diagram shows isostructural solid–solid equilibria with a critical point. To test more simplified models, phase diagrams from second-order Barker–Henderson perturbation theory are compared with those from the random-phase approximation for the fluid phase and the van der Waals theory for the solid phase; this comparison shows significantly different phase diagrams. Moreover, with a potential of mean force with primary and secondary minima, calculations using second-order perturbation theory identify conditions where colloidal and protein solutions can present two fluid–fluid regions, each with a critical point; however, the higher-density fluid–fluid region is likely to be metastable. The analytic calculations described here may be useful for interpretation of experimental phase diagrams and for guiding design of separation processes.     

Stresses inside critical nuclei

The usual derivation of classical nucleation theory is inappropriate for crystal nucleation. In particular, it leads to a seriously flawed estimate of the pressure inside a critical nucleus. This has consequences for the prediction of possible metastable phases during the nucleation process. In this paper, we reanalyze the theory for crystal nucleation based on the thermodynamics of small crystals suspended in a liquid, due to Mullins (J. Chem. Phys. 1984, 81, 1436). As an illustration of the difference between the classical picture and the present approach, we consider a numerical study of crystal nucleation in binary mixtures of hard spherical colloids with a size ratio of 1:10. The stable crystal phase of this system can be either dense or expanded. We find that, in the vicinity of the solid-solid critical point where the crystallites are highly compressible, small crystal nuclei are less dense than large nuclei. This phenomenon cannot be accounted for by either classical nucleation theory or by the Gibbsian droplet model.     

Phase field theory of interfaces and crystal nucleation in a eutectic system of fcc structure: II. Nucleation in the metastable liquid immiscibility region

In the second part of our paper, we address crystal nucleation in the metastable liquid miscibility region of eutectic systems that is always present, though experimentally often inaccessible. While this situation resembles the one seen in single component crystal nucleation in the presence of a metastable vapor-liquid critical point addressed in previous works, it is more complex because of the fact that here two crystal phases of significantly different compositions may nucleate. Accordingly, at a fixed temperature below the critical point, six different types of nuclei may form: two liquid-liquid nuclei: two solid-liquid nuclei; and two types of composite nuclei, in which the crystalline core has a liquid "skirt," whose composition falls in between the compositions of the solid and the initial liquid phases, in addition to nuclei with concentric alternating composition shells of prohibitively high free energy. We discuss crystalline phase selection via exploring/identifying the possible pathways for crystal nucleation.     

On the coupling between slow diffusion transport and barrier crossing in nucleation

We model the coupling between slow diffusion transport and nucleation using the diffusion equation, an Ostwald-Freundlich boundary condition, and a mass balance linking nucleus size to flux across the nucleus-solution interface. The model retains some characteristics of the classical nucleation theory because of the common theoretical foundations behind classical nucleation theory and the Ostwald-Freundlich equation. For example, the classically critical-sized nucleus in the uniform supersaturated concentration field is an unstable equilibrium point. However, the model also shows that certain types of concentration profiles can drive a classically pre-critical nucleus over the nucleation barrier. We identify the separatrix as a function of both nucleus size and characteristics of the local concentration field. Our analysis may be useful for understanding the effects of local concentration fluctuations and especially for understanding the role of dense precursor particles in driving two-step nucleation processes. Our analysis may also provide a starting point for further statistical field theory analyses of local concentration fluctuations and their effects on nucleation rates.     

Phase behavior of self-associating fluids with weaker dispersion interactions between bonded particles

In this study, we explore the global phase behavior of a simple model for self-associating fluids where association reduces the strength of the dispersion interactions between bonded particles. Recent research shows that this type of behavior likely explains the thermodynamic properties of strongly polar fluids and certain micellar solutions. Based on Wertheim's theory of associating liquids [M. S. Wertheim, J. Stat. Phys. 42, 459 (1986); 42, 477 (1986)], our model takes into account the effect that dissimilar particle interactions have on the equilibrium constant for self-association in the system. We find that weaker interactions between bonded molecules tend to favor the dissociation of chains at any temperature and density. This effect stabilizes a monomeric liquid phase at high densities, enriching the global phase behavior of the system. In particular, for systems in which the energy of mixing between bonded and unbonded species is positive, we find a triple point involving a vapor, a dense phase of chain aggregates, and a monomeric liquid. Phase coexistence between the vapor and the monomeric fluid is always more stable at temperatures above the triple point, but a highly associated fluid may exist as a metastable phase under these conditions. The presence of this metastable phase may explain the characteristic nucleation behavior of the liquid phase in strongly dipolar fluids.     

Nucleation of bulk phases in the HCl/H2O system

We report experimental results on the low-temperature uptake of HCl on H(2)O ice (ice). HCl was deposited on the surface at greater than monolayer amounts at 85 K, and the ice substrate was heated. The temperature dependence of the HCl vapor pressure from this phase was measured from 110 to 150 K, with the nucleation of a bulk hydrate phase observed at 150 K. Measurements were conducted in a closed system by simultaneous application of gas phase mass spectrometry and surface spectroscopy to characterize vapor/solid equilibrium and the nucleation of bulk hydrate phases. Combining the nucleation data reported here with data we reported previously (180 to 200 K) and data from two other laboratories (165 and 170 K), the thermodynamic boundaries for the nucleation of both the metastable bulk solution and bulk hydrate phases subsequent to monolayer adsorption of HCl have been determined. The nucleation of the metastable bulk solution phase occurs promptly at monolayer coverage at the ice/liquid coexistence boundary on the binary bulk phase diagram. The nucleation of the bulk hexahydrate occurs from this metastable solution along a locus of points defining a state of constant solution free energy. This measured free energy is -51.2 +/- 0.9 kJ/mol. Finally, the temperature dependence of the HCl vapor pressure from the low-temperature phase is reported here for the first time and is consistent with that of the metastable solution predicted by this thermodynamic model of uptake, extending the range of validity of this model of adsorption followed by bulk solution and hydrate nucleation to a lower bound in temperature of 110 K.     

Crystal nucleation for a model of globular proteins

A continuum model of globular proteins proposed by Talanquer and Oxtoby [J. Chem. Phys. 109, 223 (1998)] is investigated numerically, with particular emphasis on the region near the metastable fluid-fluid coexistence curve. Classical nucleation theory is shown to be invalid not only in the vicinity of the metastable critical point but also close to the liquidus line. An approximate analytic solution is also presented for the shape and properties of the nucleating crystal droplet.     

A finite-size scaling study of a model of globular proteins

Grand canonical Monte Carlo simulations are used to explore the metastable fluid-fluid coexistence curve of the modified Lennard-Jones model of globular proteins of ten Wolde and Frenkel [Science, 277, 1975 (1997)]. Using both mixed-field finite-size scaling and histogram-reweighting methods, the joint distribution of density and energy fluctuations is analyzed at coexistence to accurately determine the critical-point parameters. The subcritical coexistence region is explored using the recently developed hyper parallel tempering Monte Carlo simulation method along with histogram reweighting to obtain the density distributions. The phase diagram for the metastable fluid-fluid coexistence curve is calculated in close proximity to the critical point, a region previously unattained by simulations.     

Thermodynamics of diamond nucleation on the nanoscale

To have a clear insight into the diamond nucleation upon the hydrothermal synthesis and the reduction of carbide (HSRC), we performed the thermodynamic approach on the nanoscale to elucidate the diamond nucleation taking place in HSRC supercritical-fluid systems taking into account the capillary effect of the nanosized curvature of the diamond critical nuclei, based on the carbon thermodynamic equilibrium phase diagram. These theoretical analyses showed that the nanosize-induced interior pressure of diamond nuclei could drive the metastable phase region of the diamond nucleation in HSRC into the new stable phase region of diamond in the carbon phase diagram. Accordingly, the diamond nucleation is preferable to the graphite phase formation in the competing growth between diamond and graphite upon HSRC. Meanwhile, we predicted that 400 MPa should be the threshold pressure for the diamond synthesis by HSRC in the metastable phase region of diamond, based on the proposed thermodynamic nucleation on the nanoscale.     

Solution crystallisation via a submerged liquid-liquid phase boundary: oiling out

In many situations the process of crystallisation from solution is known to occur via metastable crystalline states (polymorphs or solvates). Here we present what we believe to be a novel example of small molecule crystallisation in which the initial separation of a solute rich liquid phase precedes the crystallisation event. We believe this occurs because a submerged liquid-liquid phase boundary is accessible within the metastable zone of the crystal nucleation process.     

Nucleation in a simple model for protein solutions with anisotropic interactions

A lattice analog of density functional theory is used to explore the structural and thermodynamic properties of critical nuclei in mixtures of particles with attractive anisotropic interactions. Protein molecules are assumed to occupy the sites on a regular cubic lattice, with effective directional interactions that mimic hydrogen bonding and the solvation forces induced by water. Interaction parameters are chosen to qualitatively reproduce the phase behavior of protein solutions. Our model predicts that critical nuclei of the solidlike phase have nonspherical shapes, and that their specific geometry depends on the nature of the anisotropic interactions. Molecules tend to align in distinctive ways in the core and in the interfacial region of these critical clusters, and the width and structure of the interface are highly affected by the presence of a metastable fluid-fluid critical point. Close to the critical region, the height of the barrier to nucleation is strongly reduced; this effect is enhanced by increasing the anisotropy of the intermolecular interactions. Unlike systems with short-range isotropic interactions, nucleation in our model is initiated by highly ordered clusters in which the order-disorder transition is confined to the interfacial region.     

Nucleation on cylindrical plates: sharp transitions and double barriers

We apply methods of density-functional theory in statistical mechanics to study the properties of droplets and bubbles formed on a single cylindrical plate or between two such disks immersed in a metastable fluid. Our approach allows us to analyze the properties of different types of aggregates and investigate the effect of disk size, disk separation, and solid-fluid interactions on the dynamics of a liquid-vapor phase transition. The finite size of disks induces nucleation phenomena that are not observed in the cases of either a planar wall or a slit pore. Heterogeneous nucleation on a single disk is characterized by the existence of two distinct types of critical nuclei that control the phase-transition dynamics at different supersaturations. Asymmetric droplets or bubbles formed on one side of the disk are the preferred nucleation path at high supersaturations. However, these types of aggregates become unstable close to the binodal, where they abruptly collapse into nuclei that engulf the cylindrical plates. Droplet or bubble nucleation in between two disks may occur through a free-energy barrier with one or two maxima depending on the value of the system parameters and the supersaturation. Metastable droplets or bubbles corresponding to local minima of the free energy are observed forming between two plates only after density fluctuations in the system achieve a critical size. These types of aggregates only exist for cylindrical plates larger than a minimum size given a fixed distance between the disks. The stability of these droplets and bubbles decreases when the plates are separated.     

Nucleation of self-associating fluids: free versus activated association

We use the density functional theory of statistical mechanics in a square gradient approximation to analyze the structure, size, and work of formation of critical nuclei in self-associating fluids where association reduces the strength of the interactions between bonded particles. This effect is expected in systems of strongly dipolar particles that associate into chains. In this work we analyze the nucleation behavior of two types of self-associating fluids: a system comprised of particles that can freely associate, and a system in which the association process involves a thermally activated initiation step. For the first case, we explore the properties of critical nuclei in fluids that exhibit a metastable critical point between a vapor phase and a highly associated liquid phase. In fluids where the association dynamics involves an initiation step, we investigate the nucleation behavior in the vicinity of the polymerization transition. In both cases critical nuclei undergo a structural transition that shares many of the features of the coil-globule transition reported in Monte Carlo simulations of strongly dipolar Stockmayer fluids. Our results suggest that the sharp structural transition observed in these simulations is evidence of the existence of a second-order or nearly second-order association transition in these model fluids.     

Lattice model of equilibrium polymerization. V. Scattering properties and the width of the critical regime for phase separation

Dynamic clustering associated with self-assembly in many complex fluids can qualitatively alter the shape of phase boundaries and produce large changes in the scale of critical fluctuations that are difficult to comprehend within the existing framework of theories of critical phenomena for nonassociating fluids. In order to elucidate the scattering and critical properties of associating fluids, we consider several models of equilibrium polymerization that describe widely occurring types of associating fluids at equilibrium and that exhibit the well defined cluster geometry of linear polymer chains. Specifically, a Flory-Huggins-type lattice theory is used, in conjunction with the random phase approximation, to compute the correlation length amplitude xi(o) and the Ginzburg number Gi corresponding, respectively, to the scale of composition fluctuations and to a parameter characterizing the temperature range over which Ising critical behavior is exhibited. Our calculations indicate that upon increasing the interparticle association energy, the polymer chains become increasingly long in the vicinity of the critical point, leading naturally to a more asymmetric phase boundary. This increase in the average degree of polymerization implies, in turn, a larger xi(o) and a drastically reduced width of the critical region (as measured by Gi). We thus obtain insight into the common appearance of asymmetric phase boundaries in a wide range of "complex" fluids and into the observation of apparent mean field critical behavior even rather close to the critical point.