The effect of sequence on the conformational stability of a model heteropolymer in explicit water

We investigate the properties of a two-dimensional lattice heteropolymer model for a protein in which water is explicitly represented. The model protein distinguishes between hydrophobic and polar monomers through the effect of the hydrophobic monomers on the entropy and enthalpy of the hydrogen bonding of solvation shell water molecules. As experimentally observed, model heteropolymer sequences fold into stable native states characterized by a hydrophobic core to avoid unfavorable interactions with the solvent. These native states undergo cold, pressure, and thermal denaturation into distinct configurations for each type of unfolding transition. However, the heteropolymer sequence is an important element, since not all sequences will fold into stable native states at positive pressures. Simulation of a large collection of sequences indicates that these fall into two general groups, those exhibiting highly stable native structures and those that do not. Statistical analysis of important patterns in sequences shows a strong tendency for observing long blocks of hydrophobic or polar monomers in the most stable sequences. Statistical analysis also shows that alternation of hydrophobic and polar monomers appears infrequently among the most stable sequences. These observations are not absolute design rules and, in practice, these are not sufficient to rationally design very stable heteropolymers. We also study the effect of mutations on improving the stability of the model proteins, and demonstrate that it is possible to obtain a very stable heteropolymer from directed evolution of an initially unstable heteropolymer.     

Triangle Distribution and Equation of State for Classical Rigid Disks

The triangle distribution function f ⁽³⁾ for three mutual near neighbors in the plane describes basic aspects of short-range order and statistical thermodynamics in two-dimensional many-particle systems. This paper examines prospects for constructing a self-consistent calculation for the rigid-disk-system f ⁽³⁾. We present several identities obeyed by f ⁽³⁾. A rudimentary closure suggested by scaled-particle theory is introduced. In conjunction with three of the basic identities, this closure leads to an unique f ⁽³⁾ over the entire density range. The pressure equation of state exhibits qualitatively correct behaviors in both the low-density and the close-packed limits, but no intervening phase transition appears. We discuss extensions to improved disk closures, and to the three-dimensional rigid-sphere system.     

Integrating diffusion maps with umbrella sampling: application to alanine dipeptide

Nonlinear dimensionality reduction techniques can be applied to molecular simulation trajectories to systematically extract a small number of variables with which to parametrize the important dynamical motions of the system. For molecular systems exhibiting free energy barriers exceeding a few k(B)T, inadequate sampling of the barrier regions between stable or metastable basins can lead to a poor global characterization of the free energy landscape. We present an adaptation of a nonlinear dimensionality reduction technique known as the diffusion map that extends its applicability to biased umbrella sampling simulation trajectories in which restraining potentials are employed to drive the system into high free energy regions and improve sampling of phase space. We then propose a bootstrapped approach to iteratively discover good low-dimensional parametrizations by interleaving successive rounds of umbrella sampling and diffusion mapping, and we illustrate the technique through a study of alanine dipeptide in explicit solvent.     

Waterlike glass polyamorphism in a monoatomic isotropic Jagla model

We perform discrete-event molecular dynamics simulations of a system of particles interacting with a spherically-symmetric (isotropic) two-scale Jagla pair potential characterized by a hard inner core, a linear repulsion at intermediate separations, and a weak attractive interaction at larger separations. This model system has been extensively studied due to its ability to reproduce many thermodynamic, dynamic, and structural anomalies of liquid water. The model is also interesting because: (i) it is very simple, being composed of isotropically interacting particles, (ii) it exhibits polyamorphism in the liquid phase, and (iii) its slow crystallization kinetics facilitate the study of glassy states. There is interest in the degree to which the known polyamorphism in glassy water may have parallels in liquid water. Motivated by parallels between the properties of the Jagla potential and those of water in the liquid state, we study the metastable phase diagram in the glass state. Specifically, we perform the computational analog of the protocols followed in the experimental studies of glassy water. We find that the Jagla potential calculations reproduce three key experimental features of glassy water: (i) the crystal-to-high-density amorphous solid (HDA) transformation upon isothermal compression, (ii) the low-density amorphous solid (LDA)-to-HDA transformation upon isothermal compression, and (iii) the HDA-to-very-high-density amorphous solid (VHDA) transformation upon isobaric annealing at high pressure. In addition, the HDA-to-LDA transformation upon isobaric heating, observed in water experiments, can only be reproduced in the Jagla model if a free surface is introduced in the simulation box. The HDA configurations obtained in cases (i) and (ii) are structurally indistinguishable, suggesting that both processes result in the same glass. With the present parametrization, the evolution of density with pressure or temperature is remarkably similar to the corresponding experimental measurements on water. Our simulations also suggest that the Jagla potential may reproduce features of the HDA-VHDA transformations observed in glassy water upon compression and decompression. Snapshots of the system during the HDA-VHDA and HDA-LDA transformations reveal a clear segregation between LDA and HDA but not between HDA and VHDA, consistent with the possibility that LDA and HDA are separated by a first order transformation as found experimentally, whereas HDA and VHDA are not. Our results demonstrate that a system of particles with simple isotropic pair interactions, a Jagla potential with two characteristic length scales, can present polyamorphism in the glass state as well as reproducing many of the distinguishing properties of liquid water. While most isotropic pair potential models crystallize readily on simulation time scales at the low temperatures investigated here, the Jagla potential is an exception, and is therefore a promising model system for the study of glass phenomenology.     

Generating inherent structures of liquids: comparison of local minimization algorithms

Properties of minima obtained by two different local minimization algorithms, the conjugate gradient and the limited-memory Broyden-Fletcher-Goldfarb-Shanno, are compared for a liquid whose particles interact via a modified Lennard-Jones pair interaction. The average properties of inherent structures obtained by the two minimization procedures are shown to agree within statistical error, though for a given starting configuration the two algorithms may not quench to identical minima.