Search

The unique properties of suspensions containing both active (self-propelling) and passive matter, arising from the nonequilibrium nature of these systems, have been widely studied (e.g., enhanced diffusion, phase separation, and directed motion). Despite this, our understanding of the specific roles played by the relevant parameters of the constituent particles remains incomplete. For instance, to what extent are the velocity and density of swimmers qualitatively distinguishable when it comes to the resultant properties of the suspension as a whole, and when are they merely two different realizations of the same thing? Through the use of numerical simulations, containing both steric and hydrodynamic interactions, we investigate a new parameter, the effective energy density, and its ability to uniquely describe the dynamics and properties of a hybrid system of active and passive particles, including the rate of pair formation and the energy distribution amongst different constituent elements. This parameter depends on both the density and the swimming velocity of the active elements, unifying them into a single variable that surpasses the descriptive ability of either alone., Comment: 8 pages, 8 figures

Dissipative particle dynamics is a widely used mesoscale technique for the simulation of hydrodynamics (as well as immersed particles) utilizing coarse-grained molecular dynamics. While the method is capable of describing any fluid, the typical choice of the friction coefficient $\gamma$ and dissipative force cutoff $r_c$ yields an unacceptably low Schmidt number $Sc$ for the simulation of liquid water at standard temperature and pressure. There are a variety of ways to raise $Sc$, such as increasing $\gamma$ and $r_c$, but the relative cost of modifying each parameter (and the concomitant impact on numerical accuracy) has heretofore remained undetermined. We perform a detailed search over the parameter space, identifying the optimal strategy for the efficient and accuracy-preserving scaling of $Sc$, using both numerical simulations and theoretical predictions. The composite results recommend a parameter choice that leads to a speed improvement of a factor of three versus previously utilized strategies., Comment: Incorrect paper submitted

We implement the replica exchange molecular dynamics algorithm to study the interactions of a model peptide (WALP-16) with an explicitly represented DPPC membrane bilayer. We observe the spontaneous, unbiased insertion of WALP-16 into the DPPC bilayer and its folding into an a-helix with a trans-bilayer orientation. We observe that the insertion of the peptide into the DPPC bilayer precedes secondary structure formation. Although the peptide has some propensity to form a partially helical structure in the interfacial region of the DPPC/water system, this state is not a productive intermediate but rather an off-pathway trap for WALP-16 insertion. Equilibrium simulations show that the observed insertion/folding pathway mirrors the potential of mean force (PMF). Calculation of the enthalpic and entropic contributions to this PMF show that the surface bound conformation of WALP-16 is significantly lower in energy than other conformations, and that the insertion of WALP-16 into the bilayer without regular secondary structure is enthalpically unfavorable by 5-10 kcal/mol/residue. The observed insertion/folding pathway disagrees with the dominant conceptual model, which is that a surface bound helix is an obligatory intermediate for the insertion of a-helical peptides into lipid bilayers. In our simulations, the observed insertion/folding pathway is favored because of a large (> 100 kcal/mol) increase in system entropy that occurs when the unstructured WALP-16 peptide enters the lipid bilayer interior. The insertion/folding pathway that is lowest in free energy depends sensitively on the near cancellation of large enthalpic and entropic terms. This suggests that intrinsic membrane peptides may have a diversity of insertion/folding behaviors depending on the exact system of peptide and lipid under consideration., Comment: 30 pages, 5 figures

The interplay between structure-search of the native structure and desolvation in protein folding has been explored using a minimalist model. These results support a folding mechanism where most of the structural formation of the protein is achieved before water is expelled from the hydrophobic core. This view integrates water expulsion effects into the funnel energy landscape theory of protein folding. Comparisons to experimental results are shown for the SH3 protein. After the folding transition, a near-native intermediate with partially solvated hydrophobic core is found. This transition is followed by a final step that cooperatively squeezes out water molecules from the partially hydrated protein core., Comment: Proceedings of the National Academy of Science, 2002, Vol.99. 685-690

Recent developments in molecular theories and simulation of ions and polar molecules in water are reviewed. The hydration of imidazole and imidazolium solutes is used to exemplify the theoretical issues. The treatment of long-ranged electrostatic interactions in simulations is discussed extensively. It is argued that the Ewald approach is an easy way to get correct hydration free energies in the thermodynamic limit from molecular calculations; and that molecular simulations with Ewald interactions and periodic boundary conditions can also be more efficient than many common alternatives. The Ewald treatment permits a conclusive extrapolation to infinite system size. The picture that emerges from testing of simple models is that the most prominent failings of the simplest theories are associated with solvent proton conformations that lead to non-gaussian fluctuations of electrostatic potentials. Thus, the most favorable cases for the second-order perturbation theories are monoatomic positive ions. For polar and anionic solutes, continuum or gaussian theories are less accurate. The appreciation of the specific deficiencies of those simple models have led to new concepts, multistate gaussian and quasi-chemical theories, that address the cases for which the simpler theories fail. It is argued that, relative to direct dielectric continuum treatments, the quasi-chemical theories provide a better theoretical organization for the computational study of the electronic structure of solution species., Comment: invited J. Phys. Chem. feature article, 30 pages including 10 figures, minor stylistic revisions

Free energies of ionic solvation calculated from computer simulations exhibit a strong system size dependence. We perform a finite-size analysis based on a dielectric-continuum model with periodic boundary conditions. That analysis results in an estimate of the Born ion size. Remarkably, the finite-size correction applies to systems with only eight water molecules hydrating a sodium ion and results in an estimate of the Born radius of sodium that agrees with the experimental value., Comment: 2 EPS figures

Theories of solvation free energies often involve electrostatic potentials at the position of a solute charge. Simulation calculations that apply cutoffs and periodic boundary conditions based on molecular centers result in center-dependent contributions to electrostatic energies due to a systematic sorting of charges in radial shells. This sorting of charges induces a surface-charge density at the cutoff sphere or simulation-box boundary that depends on the choice of molecular centers. We identify a simple solution that gives correct, center-independent results, namely the radial integration of charge densities. Our conclusions are illustrated for a Lennard-Jones solute in water. The present results can affect the parameterization of force fields., Comment: 10 pages, 3 Encapsulated Postscript figures, LaTeX, uses RevTeX macro

An information theory model is used to construct a molecular explanation why hydrophobic solvation entropies measured in calorimetry of protein unfolding converge at a common temperature. The entropy convergence follows from the weak temperature dependence of occupancy fluctuations for molecular-scale volumes in water. The macroscopic expression of the contrasting entropic behavior between water and common organic solvents is the relative temperature insensitivity of the water isothermal compressibility. The information theory model provides a quantitative description of small molecule hydration and predicts a negative entropy at convergence. Interpretations of entropic contributions to protein folding should account for this result., Comment: Phys. Rev. Letts. (in press 1996), 3 pages, 3 figures

We study the chemical potential of water as a function of charge based on perturbation theory. By calculating the electrostatic-energy fluctuations of two states (fully charged and uncharged) we are able to determine accurate values for the dependence of the chemical potential on charge. We find identical results for the chemical-potential difference of fully charged and uncharged water from overlapping-histogram and acceptance-ratio methods and by smoothly connecting the curves of direct exponential averages. Our results agree with those of Rick and Berne (J. Am. Chem. Soc., 1994, 116, 3949) with respect to both the chemical-potential difference and its dependence on the charge coupling parameter. We observe significant deviations from simple Gaussian-fluctuation statistics. The dependence on the coupling parameter is not quadratic, as would be inferred from linear continuum models of electrostatics., Comment: 10 pages, 4 Postscript figures, uses revtex.sty

The hydration free energies of ions exhibit an approximately quadratic dependence on the ionic charge, as predicted by the Born model. We analyze this behavior using second-order perturbation theory. This provides effective methods to calculating free energies from equilibrium computer simulations. The average and the fluctuation of the electrostatic potential at charge sites appear as the first coefficients in a Taylor expansion of the free energy of charging. Combining the data from different charge states allows calculation of free-energy profiles as a function of the ionic charge. The first two Taylor coefficients of the free-energy profiles can be computed accurately from equi- librium simulations; but they are affected by a strong system-size dependence. We apply corrections for these finite-size effects by using Ewald lattice sum- mation and adding the self-interactions consistently. Results are presented for a model ion with methane-like Lennard-Jones parameters in SPC water. We find two very closely quadratic regimes with different parameters for positive and negative ions. We also studied the hydration free energy of potassium, calcium, fluoride, chloride, and bromide ions. We find negative ions to be solvated more strongly compared to positive ions of equal size, in agreement with experimen- tal data. We ascribe this preference of negative ions to their strong interac- tions with water hydrogens, which can penetrate the ionic van der Waals shell without direct energetic penalty in the models used. We also find a positive electrostatic potential at the center of uncharged Lennard-Jones particles in water, which favors negative ions. Regarding finite-system-size effects, we show that even using only 16 water molecules it is possible to calculate accu- rately the hydration free energy of sodium if self-interactions are considered., Comment: 17 pages, 6 Postscript figures, uses revtex.sty

Recent molecular simulation and integral equation results alkali-halide ion pair potentials-of-mean-force in water are discussed. Dielectric model calculations are implemented to check that these models produce that characteristic structure of contact and solvent-separated minima for oppositely charged ions in water under physiological thermodynamic conditions. Comparison of the dielectric model results with the most current molecular level information indicates that the dielectric model does not, however, provide an accurate description of these potentials-of-mean-force. We note that linear dielectric models correspond to modelistic implementations of second-order thermodynamic perturbation theory for the excess chemical potential of a distinguished solute molecule. Therefore, the molecular theory corresponding to the dielectric models is second-order thermodynamic perturbation theory for that excess chemical potential. The second-order, or fluctuation, term raises a technical computational issue of treatment of long-ranged interactions similar to the one which arises in calculation of the dielectric constant of the solvent. It is contended that the most important step for further development of dielectric models would be a separate assessment of the first-order perturbative term (equivalently the {\it potential at zero charge} ) which vanishes in the dielectric models but is generally nonzero. Parameterization of radii and molecular volumes should then be based of the second-order perturbative term alone. Illustrative initial calculations are presented and discussed., Comment: 37 pages and 8 figures. LA-UR-93-4205