Search

Molecular dynamics was used to compute radical and orientational distribution functions of water molecules around the different molecules: argon, methane and neopentane. It was found that the orientational distribution of water molecules on a narrow shell proximate to the neopentane molecule was very similar to that near a paraffin wall, in contrast to argon and methane.

Heterogeneous ice growth exhibits a maximum in freezing rate arising from the competition between kinetics and the thermodynamic driving force between the solid and liquid states. Here, we use molecular dynamics simulations to elucidate the atomistic details of this competition, focusing on water properties in the interfacial region along the secondary prismatic direction. The crystal growth velocity is maximized when the efficiency of converting interfacial water molecules to ice, collectively known as the attachment kinetics, is greatest. We find water molecules that contact the intermediate ice layer in concave regions along the atomistically roughened surface are more likely to freeze directly. An increased roughening of the solid surface at large undercoolings consequently plays an important limiting role in the rate of ice growth, as water molecules are unable to integrate into increasingly deeper surface pockets. These results provide insight into the molecular mechanisms for self-assembly of solid phases that are important in many biological and atmospheric processes.

Longstanding mechanistic questions about the role of protecting osmolyte trimethylamine N-oxide (TMAO) that favors protein folding and the denaturing osmolyte urea are addressed by studying their effects on the folding of uncharged polymer chains. Using atomistic molecular dynamics simulations, we show that 1 M TMAO and 7 M urea solutions act dramatically differently on these model polymer chains. Their behaviors are sensitive to the strength of the attractive dispersion interactions of the chain with its environment: when these dispersion interactions are sufficiently strong, TMAO suppresses the formation of extended conformations of the hydrophobic polymer as compared to water while urea promotes the formation of extended conformations. Similar trends are observed experimentally for real protein systems. Quite surprisingly, we find that both protecting and denaturing osmolytes strongly interact with the polymer, seemingly in contrast with existing explanations of the osmolyte effect on proteins. We show that what really matters for a protective osmolyte is its effective depletion as the polymer conformation changes, which leads to a negative change in the preferential binding coefficient. For TMAO, there is a much more favorable free energy of insertion of a single osmolyte near collapsed conformations of the polymer than near extended conformations. By contrast, urea is preferentially stabilized next to the extended conformation and thus has a denaturing effect.

We study the kinetics of assembly of two plates of varying hydrophobicity, including cases where drying occurs and water strongly solvates the plate surfaces. The potential of mean force and molecular-scale hydrodynamics are computed from molecular dynamics simulations in explicit solvent as a function of particle separation. In agreement with our recent work on nanospheres [J. Phys. Chem. B 2012,116, 378-389], regions of high friction are found to be engendered by large and slow solvent fluctuations. These slow fluctuations can be due to either drying or confinement. The mean first passage times for assembly are computed by means of molecular dynamics simulations in explicit solvent and by Brownian dynamics simulations along the reaction path. Brownian dynamics makes use of the potential of mean force and hydrodynamic profile that we determined. Surprisingly, we find reasonable agreement between full-scale molecular dynamics and Brownian dynamics, despite the role of slow solvent relaxation in the assembly process. We found that molecular-scale hydrodynamic interactions are essential in describing the kinetics of assembly.

Solvent plays an important role in the relative motion of nanoscopic bodies, and the study of such phenomena can help elucidate the mechanism of hydrophobic assembly, as well as the influence of solvent-mediated effects on in vivo motion in crowded cellular environments. Here we study important aspects of this problem within the framework of Brownian dynamics. We compute the free energy surface that the Brownian particles experience and their hydrodynamic interactions from molecular dynamics simulations in explicit solvent. We find that molecular scale effects dominate at short distances, thus giving rise to deviations from the predictions of continuum hydrodynamic theory. Drying phenomena, solvent layering, and fluctuations engender distinct signatures of the molecular scale. The rate of assembly in the diffusion-controlled limit is found to decrease from molecular scale hydrodynamic interactions, in opposition to the free energy driving force for hydrophobic assembly, and act to reinforce the influence of the free energy surface on the association of more hydrophilic bodies.

A small change in the Hamiltonian scaling in Replica Exchange with Solute Tempering (REST) is found to improve its sampling efficiency greatly, especially for the sampling of aqueous protein solutions in which there are large-scale solute conformation changes. Like the original REST (REST1), the new version (which we call REST2) also bypasses the poor scaling with system size of the standard Temperature Replica Exchange Method (TREM), reducing the number of replicas (parallel processes) from what must be used in TREM. This reduction is accomplished by deforming the Hamiltonian function for each replica in such a way that the acceptance probability for the exchange of replica configurations does not depend on the number of explicit water molecules in the system. For proof of concept, REST2 is compared with TREM and with REST1 for the folding of the trpcage and β-hairpin in water. The comparisons confirm that REST2 greatly reduces the number of CPUs required by regular replica exchange and greatly increases the sampling efficiency over REST1. This method reduces the CPU time required for calculating thermodynamic averages and for the ab initio folding of proteins in explicit water.

The binding affinity between a probe hydrophobic particle and model hydrophobic plates with different charge (or dipole) densities in water was investigated through molecular dynamics simulation free-energy perturbation calculations. We observed a reduced binding affinity when the plates are charged, in agreement with previous findings. With increased charge density, the plates can change from "hydrophobic like" (pulling the particle into the interplate region) to "hydrophilic like" (ejecting the particle out of the interplate region), demonstrating the competition between hydrophobic and electrostatic interactions. The reduction of the binding affinity is quadratically dependent on the magnitude of the charge for symmetric systems, but linear and cubic terms also make a contribution for asymmetric systems. Statistical perturbation theory explains these results and shows when and why implicit solvent models fail.

We studied by molecular dynamics simulations the temperature dependence of hydrophobic association and drying transition of large-scale solutes. Similar to the behavior of small solutes, we found the association process to be characterized by a large negative heat capacity change. The origin of this large change in heat capacity is the high fragility of hydrogen bonds between water molecules at the interface with hydrophobic solutes; an increase in temperature breaks more hydrogen bonds at the interface than in the bulk. With increasing temperature, both entropy and enthalpy changes for association strongly decrease, while the change in free energy weakly varies, exhibiting a small minimum at high temperatures. At around T=Ts=360 K, the change in entropy is zero, a behavior similar to the solvation of small nonpolar solutes. Unexpectedly, we find that at Ts, there is still a substantial orientational ordering of the interfacial water molecules relative to the bulk. Nevertheless, at this point, the change in entropy vanishes due to a compensating contribution of translational entropy. Thus, at Ts, there is rotational order and translational disorder of the interfacial water relative to bulk water. In addition, we studied the temperature dependence of the drying-wetting transition. By calculating the contact angle of water on the hydrophobic surface at different temperatures, we compared the critical distance observed in the simulations with the critical distance predicted by macroscopic theory. Although the deviations of the predicted from the observed values are very small (8-23%), there seems to be an increase in the deviations with an increase in temperature. We suggest that these deviations emerge due to increased fluctuations, characterizing finite systems, as the temperature increases.

We apply the recently developed replica exchange with solute tempering (REST) to three large solvated peptide systems: an alpha-helix, a beta-hairpin, and a TrpCage, with these peptides defined as the "central group". We find that our original implementation of REST is not always more efficient than the replica exchange method (REM). Specifically, we find that exchanges between folded (F) and unfolded (U) conformations with vastly different structural energies are greatly reduced by the nonappearance of the water self-interaction energy in the replica exchange acceptance probabilities. REST, however, is expected to remain useful for a large class of systems for which the energy gap between the two states is not large, such as weakly bound protein-ligand complexes. Alternatively, a shell of water molecules can be incorporated into the central group, as discussed in the original paper.

Parallel tempering (or the replica exchange method (REM)) is a powerful method for speeding up the sampling of conformational states of systems with rough energy landscapes, like proteins, where stable conformational states can be separated by large energy barriers. The usual implementation of the REM is performed on local computer clusters (or parallel processors) where the different replicas must be run synchronously. Here, we present serial replica exchange (SREM), a method that is equivalent to the standard REM in terms of efficiency yet runs asynchronously on a distributed network of computers. A second advantage is the method's greatly enhanced fault tolerance, which enables the study of biological systems on worldwide distributed computing environments, such as Folding@Home. For proof of concept, we apply the SREM to a single alanine dipeptide molecule in explicit water. We show that the SREM reproduces the thermodynamic and structural properties determined by the REM.

The effect of salts on the solvent-induced interactions between hydrophobic particles dispersed in explicit aqueous solution is investigated as a function of the salt's ionic charge density by molecular dynamics simulations. We demonstrate that aggregates of the hydrophobic particles can be formed or dissolved in response to changes in the charge density of the ions. Ions with high charge density increase the propensity of the hydrophobic particles to aggregate. This corresponds to stronger hydrophobic interactions and a decrease in the solubility (salting-out) of the hydrophobic particles. Ions with low charge density can either increase or decrease the propensity for aggregation depending on whether the concentration of the salt is low or high, respectively. At low concentrations of low charge density ions, the aggregate forms a "micelle-like" structure in which the ions are preferentially adsorbed at the surface of the aggregate. These "micelle-like" structures can be soluble in water so that the electrolyte can both increase the solubility and increase aggregation at the same time. We also find, that at the concentration of the hydrophobic particles studied (approximately 0.75 m), the aggregation process resembles a first-order transition in finite systems.

A hydrophobic aided replica exchange method (HAREM) is introduced to accelerate the simulation of all-atom protein folding in explicit solvent. This method is based on exaggerating the hydrophobic effect of various protein amino acids in water by attenuating the protein-water attractive interactions (mimicking the Chaperon effect) while leaving other interactions among protein atoms and water molecules unchanged. The method is applied to a small representative protein, the alpha-helix 3K(I), and it is found that the HAREM method successfully folds the protein within 4 ns, while the regular replica exchange method does not fold the same protein within 5 ns, even with many more replicas.

The fluctuating elastic boundary (FEB) model for molecular dynamics has recently been developed and validated through simulations of liquid argon. In the FEB model, a flexible boundary which consists of particles connected by springs is used to confine the solvated system, thereby eliminating the need for periodic boundary conditions. In this study, we extend this model to the simulation of bulk water and solvated alanine dipeptide. Both the confining potential and boundary particle interaction functions are modified to preserve the structural integrity of the boundary and prevent the leakage of the solute-solvent system through the boundary. A broad spectrum of structural and dynamic properties of liquid water are computed and compared with those obtained from conventional periodic boundary condition simulations. The applicability of the model to biomolecular simulations is investigated through the analysis of conformational population distribution of solvated alanine dipeptide. In most cases we find remarkable agreement between the two simulation approaches.

Molecular dynamics simulations are performed to study the dynamics of interfacial water confined in the interdomain region of a two-domain protein, BphC enzyme. The results show that near the protein surface the water diffusion constant is much smaller and the water-water hydrogen bond lifetime is much longer than that in bulk. The diffusion constant and hydrogen bond lifetime can vary by a factor of as much as 2 in going from the region near the hydrophobic domain surface to the bulk. Water molecules in the first solvation shell persist for a much longer time near local concave sites than near convex sites. Also, the water layer survival correlation time shows that on average water molecules near the extended hydrophilic surfaces have longer residence times than those near hydrophobic surfaces. These results indicate that local surface curvature and hydrophobicity have a significant influence on water dynamics.

To isolate the effects of the inclusion of polarizability in the force field model on the structure and dynamics of the solvating water in differing electrostatic environments of proteins, we present the results of molecular dynamics simulations of the bovine pancreatic trypsin inhibitor (BPTI) in water with force fields that explicitly include polarization for both the protein and the water. We use three model potentials for water and two model potentials for the protein. Two of the water models and one of the protein models are polarizable. A total of six systems were simulated representing all combinations of these polarizable and nonpolarizable protein and water force fields. We find that all six systems behave in a similar manner in regions of the protein that are weakly electrostatic (either hydrophobic or weakly hydrophilic). However, in the vicinity of regions of the protein with relatively strong electrostatic fields (near positively or negatively charged residues), we observe that the water structure and dynamics are dependent on both the model of the protein and the model of the water. We find that a large part of the dynamical dependence can be described by small changes in the local environments of each region that limit the local density of non-hydrogen-bonded waters, precisely the water molecules that facilitate the dynamical relaxation of the water-water hydrogen bonds. We introduce a simple method for rescaling for this effect. When this is done, we are able to effectively isolate the influence of polarizability on the dynamics. We find that the solvating water's relaxation is most affected when both the protein and the water models are polarizable. However, when only one model (or neither) is polarizable, the relaxation is similar regardless of the models used.

Hydrogen-bond (H-bond) dynamics in the air-water interface is studied by molecular dynamics simulations. The analysis reveals that the dynamics of breaking and forming hydrogen bonds in the air-water interface is faster than that in bulk water for the polarizable water models. This is in contrast to the results found on a protein surface. We show that the difference stems from more rapid translational diffusion in the interface. When the effect of pair diffusion is eliminated, the hydrogen-bond dynamics in the interface is observed to be slower than that in the bulk. This occurs because the number of water molecules adjacent to a hydrogen-bonded pair and available to accept or donate a hydrogen bond is smaller in the interface than in the bulk. The comparison between polarizable water models and fixed-charge models highlights the potential importance of the polarization effect in the water-vapor interface.

A new approach is developed to study the dynamics of the localized process in solutions and other condensed phase systems. The approach employs a fluctuating elastic boundary (FEB) model which encloses the simulated system in an elastic bag that mimics the effects of the bulk solvent. This alleviates the need for periodic boundary conditions and allows for a reduction in the number of solvent molecules that need to be included in the simulation. The boundary bag is modeled as a mesh of quasi-particles connected by elastic bonds. The FEB model allows for volume and density fluctuations characteristic of the bulk system, and the shape of the boundary fluctuates during the course of the simulation to adapt to the configuration fluctuations of the explicit solute-solvent system inside. The method is applied to the simulation of a Lennard-Jones model of liquid argon. Various structural and dynamical quantities are computed and compared with those obtained from conventional periodic boundary simulations. The agreement between the two is excellent in most cases, thus validating the viability of the FEB method.