In the absence of experimental assignments, the empirical charge/hydropathy correlation for the prediction of natively unfolded protein sequences (Uversky, V. N.; Gillespie, J. R.; Fink, A. L. Proteins: Struct., Funct., Genet. 2000, 41, 415-427) provides perhaps the most intuitive description of gross polypeptide conformation. The success of this correlation rests on an essential chain length independence of the boundary line between expanded and compact conformations, conversely stabilized by highly charged/weakly hydrophobic residues or weakly charged/highly hydrophobic residues, respectively. We present extensive simulation results for coarse-grained polypeptides over a wide range of sequence hydrophobicities, charges, and lengths. A coil-to-globule transition in sequence composition space analogous to the charge/hydropathy correlation is observed. A near sequence length independent stability boundary is only found when counterions for the charged peptides are explicitly included, as a result of counterion condensation stabilization of repulsive electrostatic interactions on the globule surface. The observed counterion adsorption is shown to be in quantitative agreement with theoretical condensation predictions. We argue that alternate functionalities, beyond charge and hydrophobicity, empirically known to correlate with conformational disorder can be incorporated into our minimalist polypeptide model to study the interplay between independent predictors of unfolded sequences.

Following the conclusions of an information theory analysis that hydrophobic hydration is dictated by the equation of state of liquid water, we perform simulations of ten different water models to examine the correlation between the fidelity of each model to the experimental density of liquid water and the accuracy of its description of methane hydration. We find that the three- and five-point water models provide an inferior description of both the liquid density and methane solubility compared to the four-point water models. Of the four-point water models, TIP4P/2005 provides the best description of both the aqueous equation-of-state and methane hydration thermodynamics. When the optimized potentials for liquid simulation united-atom description for methane is used, we find that while the entropy and heat capacity of methane hydration are in excellent agreement with experiment, the chemical potential and enthalpy are systematically shifted upwards. We subsequently reoptimize the methane interaction to accurately reproduce the experimental solubilities as a function of temperature by accounting for missing attractive interactions.

Although hot, cold, and high pressure denaturation are well characterized, the possibility of negative pressure unfolding has received much less attention. Proteins under negative pressure, however, are important in applications such as medical ultrasound, and the survival of biopoloymers in the xylem and adjacent parenchyma cells of vascular plants. In addition, negative pressure unfolding is fundamentally important in obtaining a complete understanding of protein stability and naturally complements previous studies of high pressure denaturation. We use extensive replica-exchange molecular dynamics (REMD) simulations and thermodynamic analysis to obtain folding/unfolding equilibrium phase diagrams for the miniprotein trp-cage (α-structure, 20-residue), the GB1 β-hairpin (β-structure, 16-residue), and the AK16 peptide (α-helix, 16-residue). Although the trp-cage is destabilized by negative pressure, the GB1 β-hairpin and AK16 peptide are stabilized by this condition.

The fluid phase diagram of trimer particles composed of one central attractive bead and two repulsive beads was determined as a function of simple geometric parameters using flat-histogram Monte Carlo methods. A variety of self-assembled structures were obtained including spherical micelle-like clusters, elongated clusters, and densely packed cylinders, depending on both the state conditions and shape of the trimer. Advanced simulation techniques were employed to determine transitions between self-assembled structures and macroscopic phases using thermodynamic and structural definitions. Simple changes in particle geometry yield dramatic changes in phase behavior, ranging from macroscopic fluid phase separation to molecular-scale self-assembly. In special cases, both self-assembled, elongated clusters and bulk fluid phase separation occur simultaneously. Our work suggests that tuning particle shape and interactions can yield superstructures with controlled architecture.

The Free Energy and Advanced Sampling Simulation Toolkit (FEASST) is a free, open-source, modular program to conduct molecular and particle-based simulations with Metropolis, Wang-Landau, and Transition-Matrix Monte Carlo methods. FEASST is implemented in C++ and may be imported as a module within Python 2 or 3. This document describes the initial public release version 1.0.

We systematically investigate the assembly of binary multi-flavored colloidal mixtures in two dimensions. In these mixtures all pairwise interactions between species may be tuned independently. This introduces an additional degree of freedom over more traditional binary mixtures with fixed mixing rules, which is anticipated to open new avenues for directed self-assembly. At present, colloidal self-assembly into non-trivial lattices tends to require either high pressures for isotropically interacting particles, or the introduction of directionally anisotropic interactions. Here we demonstrate tunable assembly into a plethora of structures which requires neither of these conditions. We develop a minimal model that defines a three-dimensional phase space containing one dimension for each pairwise interaction, then employ various computational techniques to map out regions of this phase space in which the system self-assembles into these different morphologies. We then present a mean-field model that is capable of reproducing these results for size-symmetric mixtures, which reveals how to target different structures by tuning pairwise interactions, solution stoichiometry, or both. Concerning particle size asymmetry, we find that domains in this model’s phase space, corresponding to different morphologies, tend to undergo a continuous “rotation” whose magnitude is proportional to the size asymmetry. Such continuity enables one to estimate the relative stability of different lattices for arbitrary size asymmetries. Owing to its simplicity and accuracy, we expect this model to serve as a valuable design tool for engineering binary colloidal crystals from multi-flavored components.

Asymmetric autocatalysis and inhibition have been proposed as key processes in the spontaneous emergence of chiral symmetry breaking in a prebiotic world. An elementary lattice model is formulated to simulate the kinetics of chiral symmetry breaking via autocatalysis and inhibition in a mixture of prochiral reactants, chiral products, and inert solvent. Starting from a chirally unbiased initial state, spontaneous symmetry breaking occurs in spite of equal a priori probability for creating either product enantiomer, and the coupled reaction-diffusion processes subsequently amplify the random early-stage symmetry breaking. The processes of reaction and diffusion are kinetically intertwined in a way leading to competition in the appearance of enantiomeric excess. An effective transition temperature can be identified below which spontaneous symmetry breaking appears. In the absence of inhibition, reactions are predominantly autocatalytic under both reaction control (fast diffusion, slow reaction) or diffusion control (fast reaction, slow diffusion) conditions. In the presence of inhibition, simulations with different system sizes converge to the same transition temperature under reaction control conditions, and in this limit the reactions are predominantly nonautocatalytic.

Trimers with one attractive bead and two repulsive beads, similar to recently synthesized trimer patchy colloids, were simulated with flat-histogram Monte Carlo methods to obtain the stable self-assembled structures for different shapes and interaction potentials. Extended corresponding states principle was successfully applied to self-assembling systems in order to approximately collapse the results for models with the same shape, but different interaction range. This helps us directly compare simulation results with previous experiment, and good agreement was found between the two. In addition, a variety of self-assembled structures were observed by varying the trimer geometry, including spherical clusters, elongated clusters, monolayers, and spherical shells. In conclusion, our results help to compare simulations and experiments, via extended corresponding states, and we predict the formation of self-assembled structures for trimer shapes that have not been experimentally synthesized.