Thomas M. Truskett

Excess-entropy-based anomalies for a waterlike fluid

Many thermodynamic and dynamic properties of water display unusual behavior at low enough temperatures. In a recent study, Yan et al. [Phys. Rev. Lett. 95, 130604 (2005)] identified a spherically symmetric two-scale potential that displays many of the same anomalous properties as water. More specifically, for select parametrizations of the potential, one finds that the regions where isothermal compression anomalously (i) decreases the fluids structural order, (ii) increases its translational self-diffusivity, and (iii) increases its entropy form nested domes in the temperature-density plane. These property relationships are similar to those found for more realistic models of water. In this work, the authors provide evidence that suggests that the anomalous regions specified above can all be linked through knowledge of the excess entropy. Specifically, the authors show how entropy scaling relationships developed by Rosenfeld [Phys. Rev. A 15, 2545 (1977)] can be used to describe the region of diffusivity anomalies and to predict the state conditions for which anomalous viscosity and thermal conductivity behavior might be found.

Layering and Position-Dependent Diffusive Dynamics of Confined Fluids

We study the diffusive dynamics of a hard-sphere fluid confined between parallel smooth hard walls. The position-dependent diffusion coefficient normal to the walls is larger in regions of high local packing density. High density regions also have the largest available volume, consistent with the fast local diffusivity. Indeed, local and global diffusivities as a function of the Widom insertion probability approximately collapse onto a master curve. Parallel and average normal diffusivities are strongly coupled at high densities and deviate from bulk fluid behavior.

Thermodynamics predicts how confinement modifies the dynamics of the equilibrium hard-sphere fluid.

We study how confining the equilibrium hard-sphere fluid to restrictive one- and two-dimensional channels with smooth interacting walls modifies its structure, dynamics, and entropy using molecular dynamics and transition-matrix Monte Carlo simulations. Although confinement strongly affects local structuring, the relationships between self-diffusivity, excess entropy, and average fluid density are, to an excellent approximation, independent of channel width or particle-wall interactions. Thus, thermodynamics can be used to predict how confinement impacts dynamics.

Relationship between thermodynamics and dynamics of supercooled liquids

Diffusivity, a measure for how rapidly a fluid self-mixes, shows an intimate, but seemingly fragmented, connection to thermodynamics. On one hand, the “configurational” contribution to entropy (related to the number of mechanically stable configurations that fluid molecules can adopt) has long been considered key for predicting supercooled liquid dynamics near the glass transition. On the other hand, the excess entropy (relative to ideal gas) provides a robust scaling for the diffusivity of fluids above the freezing point. Here we provide, to our knowledge, the first evidence that excess entropy also captures how supercooling a fluid modifies its diffusivity, suggesting that dynamics, from ideal gas to glass, is related to a single, standard thermodynamic quantity.

Concentrated dispersions of equilibrium protein nanoclusters that reversibly dissociate into active monomers

Stabilizing proteins at high concentration is of broad interest in drug delivery, for treatment of cancer and many other diseases. Herein, we create highly concentrated antibody dispersions (up to 260 mg/mL) comprising dense equilibrium nanoclusters of protein (monoclonal antibody 1B7, polyclonal sheep immunoglobulin G, and bovine serum albumin) molecules which, upon dilution in vitro or administration in vivo, remain conformationally stable and biologically active. The extremely concentrated environment within the nanoclusters (∼700 mg/mL) provides conformational stability to the protein through a novel self-crowding mechanism, as shown by computer simulation, while the primarily repulsive nanocluster interactions result in colloidally stable, transparent dispersions. The nanoclusters are formed by adding trehalose as a cosolute which strengthens the short-ranged attraction between protein molecules. The protein cluster diameter was reversibly tuned from 50 to 300 nm by balancing short-ranged attraction against long-ranged electrostatic repulsion of weakly charged protein at a pH near the isoelectric point. This behavior is described semiquantitatively with a free energy model which includes the fractal dimension of the clusters. Upon dilution of the dispersion in vitro, the clusters rapidly dissociated into fully active protein monomers as shown with biophysical analysis (SEC, DLS, CD, and SDS-PAGE) and sensitive biological assays. Since the concept of forming nanoclusters by tuning colloid interactions is shown to be general, it is likely applicable to a variety of biological therapeutics, mitigating the need to engineer protein stability through amino acid modification. In vivo subcutaneous injection into mice results in indistinguishable pharmacokinetics versus a standard antibody solution. Stable protein dispersions with low viscosities may potentially enable patient self-administration by subcutaneous injection of antibody therapeutics being discovered and developed.

Effects of lengthscales and attractions on the collapse of hydrophobic polymers in water

We present results from extensive molecular dynamics simulations of collapse transitions of hydrophobic polymers in explicit water focused on understanding effects of lengthscale of the hydrophobic surface and of attractive interactions on folding. Hydrophobic polymers display parabolic, protein-like, temperature-dependent free energy of unfolding. Folded states of small attractive polymers are marginally stable at 300 K and can be unfolded by heating or cooling. Increasing the lengthscale or decreasing the polymer–water attractions stabilizes folded states significantly, the former dominated by the hydration contribution. That hydration contribution can be described by the surface tension model, ΔG = γ(T)ΔA, where the surface tension, γ, is lengthscale-dependent and decreases monotonically with temperature. The resulting variation of the hydration entropy with polymer lengthscale is consistent with theoretical predictions of Huang and Chandler [Huang DM, Chandler D (2000) Proc Natl Acad Sci USA 97:8324–8327] that explain the blurring of entropy convergence observed in protein folding thermodynamics. Analysis of water structure shows that the polymer–water hydrophobic interface is soft and weakly dewetted, and is characterized by enhanced interfacial density fluctuations. Formation of this interface, which induces polymer folding, is strongly opposed by enthalpy and favored by entropy, similar to the vapor–liquid interface.

Anomalous structure and dynamics of the Gaussian-core fluid.

It is known that there are thermodynamic states for which the Gaussian-core fluid displays anomalous properties such as expansion upon isobaric cooling (density anomaly) and increased single-particle mobility upon isothermal compression (self-diffusivity anomaly). Here, we investigate how temperature and density affect its short-range translational structural order, as characterized by the two-body excess entropy. We find that there is a wide range of conditions for which the short-range translational order of the Gaussian-core fluid decreases upon isothermal compression (structural order anomaly). As we show, the origin of the structural anomaly is qualitatively similar to that of other anomalous fluids (e.g., water or colloids with short-range attractions) and is connected to how compression affects static correlations at different length scales. Interestingly, we find that the self-diffusivity of the Gaussian-core fluid obeys a scaling relationship with the two-body excess entropy that is very similar to the one observed for a variety of simple liquids. One consequence of this relationship is that the state points for which structural, self-diffusivity, and density anomalies of the Gaussian-core fluid occur appear as cascading regions on the temperature-density plane; a phenomenon observed earlier for models of waterlike fluids. There are, however, key differences between the anomalies of Gaussian-core and waterlike fluids, and we discuss how those can be qualitatively understood by considering the respective interparticle potentials of these models. Finally, we note that the self-diffusivity of the Gaussian-core fluid obeys different scaling laws depending on whether the two-body or total excess entropy is considered. This finding, which deserves more comprehensive future study, appears to underscore the significance of higher-body correlations for the behavior of fluids with bounded interactions.

Structural Ensemble of an Intrinsically Disordered Polypeptide

Intrinsically disordered proteins (IDPs), which play key roles in cell signaling and regulation, do not display specific tertiary structure when isolated in solution. Instead, they dynamically explore an ensemble of unfolded configurations, adopting more stable, ordered structures only after binding to their ligands. Whether ligands induce IDP structural changes upon binding or simply bind to pre-existing conformers that are populated within the IDPs structural ensemble is not well understood. Molecular simulations can provide information with the spatiotemporal resolution necessary to resolve these issues. Here, we report on the conformational ensemble of a 15-residue wild-type p53 fragment from the TAD domain and its mutant (TAD-P27L) obtained by replica exchange molecular dynamics simulation using an optimized (fully atomistic, explicit solvent) protein model and the experimental validation of the simulation results. We use a clustering method based on structural similarity to identify conformer states populated by the peptides in solution from the simulated ensemble. We show that p53 populates solution structures that strongly resemble the ligand (MDM2)-bound structure, but at the same time, the conformational free-energy landscape is relatively flat in the absence of the ligand.

Equilibrium gold nanoclusters quenched with biodegradable polymers

Although sub-100 nm nanoclusters of metal nanoparticles are of interest in many fields including biomedical imaging, sensors, and catalysis, it has been challenging to control their morphologies and chemical properties. Herein, a new concept is presented to assemble equilibrium Au nanoclusters of controlled size by tuning the colloidal interactions with a polymeric stabilizer, PLA(1k)-b-PEG(10k)-b-PLA(1k). The nanoclusters form upon mixing a dispersion of ~5 nm Au nanospheres with a polymer solution followed by partial solvent evaporation. A weakly adsorbed polymer quenches the equilibrium nanocluster size and provides steric stabilization. Nanocluster size is tuned from ~20 to ~40 nm by experimentally varying the final Au nanoparticle concentration and the polymer/Au ratio, along with the charge on the initial Au nanoparticle surface. Upon biodegradation of the quencher, the nanoclusters reversibly and fully dissociate to individual ~5 nm primary particles. Equilibrium cluster size is predicted semiquantitatively with a free energy model that balances short-ranged depletion and van der Waals attractions with longer-ranged electrostatic repulsion, as a function of the Au and polymer concentrations. The close spacings of the Au nanoparticles in the clusters produce strong NIR extinction over a broad range of wavelengths from 650 to 900 nm, which is of practical interest in biomedical imaging.

Synergistic formation and stabilization of oil-in-water emulsions by a weakly interacting mixture of zwitterionic surfactant and silica nanoparticles.

Oil-in-water emulsions were formed and stabilized at low amphiphile concentrations by combining hydrophilic nanoparticles (NPs) (i.e., bare colloidal silica) with a weakly interacting zwitterionic surfactant, caprylamidopropyl betaine, to generate a high hydrophilic-lipophilic balance. The weak interaction of the NPs with surfactant was quantified with contact angle measurements. Emulsions were characterized by static light scattering to determine the droplet size distributions, optical photography to quantify phase separation due to creaming, and both optical and electron microscopy to determine emulsion microstructure. The NPs and surfactant acted synergistically to produce finer emulsions with a greater stability to coalescence relative to the behavior with either NPs or surfactant alone. As a consequence of the weak adsorption of the highly hydrophilic surfactant on the anionic NPs along with the high critical micelle concentration, an unusually large surfactant concentration was available to adsorb at the oil-water interface and lower the interfacial tension. The synergy for emulsion formation and stabilization for the two amphiphiles was even greater in the case of a high-salinity synthetic seawater aqueous phase. Here, higher NP adsorption at the oil-water interface was caused by electrostatic screening of interactions between (1) NPs and the anionic oil-water interface and (2) between the NPs. This greater adsorption as well as partial flocculation of the NPs provided a more efficient barrier to droplet coalescence.