Direct observation of chaperone-induced changes in a protein folding pathway.

How chaperone interactions affect protein folding pathways is a central problem in biology. With the use of optical tweezers and all-atom molecular dynamics simulations, we studied the effect of chaperone SecB on the folding and unfolding pathways of maltose binding protein (MBP) at the single-molecule level. In the absence of SecB, we find that the MBP polypeptide first collapses into a molten globulelike compacted state and then folds into a stable core structure onto which several alpha helices are finally wrapped. Interactions with SecB completely prevent stable tertiary contacts in the core structure but have no detectable effect on the folding of the external alpha helices. It appears that SecB only binds to the extended or molten globulelike structure and retains MBP in this latter state. Thus during MBP translocation, no energy is required to disrupt stable tertiary interactions.

Effect of membrane environment on proton permeation through gramicidin A channels.

Multistate empirical valence bond simulations were employed to study proton transport through gramicidin A channels embedded in two different lipid bilayers, glycerol 1-monooleate (GMO) and diphytanolphosphatidylcholine (DiPhPC). Free energy barriers to proton permeation were derived using a new internal reaction coordinate describing the proton permeation process. The large quantitative and qualitative differences between the two systems are discussed in terms of local bilayer structures, ordering of interfacial water, and channel flexibility in the two environments.

Mechanisms of passive ion permeation through lipid bilayers: insights from simulations.

Multistate empirical valence bond and classical molecular dynamics simulations were used to explore mechanisms for passive ion leakage through a dimyristoyl phosphatidylcholine lipid bilayer. In accordance with a previous study on proton leakage (Biophys. J. 2005, 88, 3095), it was found that the permeation mechanism must be a highly concerted one, in which ion, solvent, and membrane coordinates are coupled. The presence of the ion itself significantly alters the response of those coordinates, suggesting that simulations of transmembrane water structures without explicit inclusion of the ionic solute are insufficient for elucidating transition mechanisms. The properties of H(+), Na(+), OH(-), and bare water molecules in the membrane interior were compared, both by biased sampling techniques and by constructing complete and unbiased transition paths. It was found that the anomalous difference in leakage rates between protons and other cations can be largely explained by charge delocalization effects rather than the usual kinetic picture (Grotthuss hopping of the proton). Permeability differences between anions and cations through phosphatidylcholine bilayers are correlated with suppression of favorable membrane breathing modes by cations.

Flexible simple point-charge water model with improved liquid-state properties.

In order to introduce flexibility into the simple point-charge (SPC) water model, the impact of the intramolecular degrees of freedom on liquid properties was systematically studied in this work as a function of many possible parameter sets. It was found that the diffusion constant is extremely sensitive to the equilibrium bond length and that this effect is mainly due to the strength of intermolecular hydrogen bonds. The static dielectric constant was found to be very sensitive to the equilibrium bond angle via the distribution of intermolecular angles in the liquid: A larger bond angle will increase the angle formed by two molecular dipoles, which is particularly significant for the first solvation shell. This result is in agreement with the work of Hochtl et al. [J. Chem. Phys. 109, 4927 (1998)]. A new flexible simple point-charge water model was derived by optimizing bulk diffusion and dielectric constants to the experimental values via the equilibrium bond length and angle. Due to the large sensitivities, the parametrization only slightly perturbs the molecular geometry of the base SPC model. Extensive comparisons of thermodynamic, structural, and kinetic properties indicate that the new model is much improved over the standard SPC model and its overall performance is comparable to or even better than the extended SPC model.

A coarse-grained model for double-helix molecules in solution: spontaneous helix formation and equilibrium properties.

A new reductionist coarse-grained model is presented for double-helix molecules in solution. As with such models for lipid bilayers and micelles, the level of description is both particulate and mesoscopic. The particulate (bead-and-spring) nature of the model makes for a simple implementation in standard molecular dynamics simulation codes and allows for investigation of thermomechanic properties without preimposing any (form of) response function. The mesoscopic level of description--where groups of atoms are condensed into coarse-grained beads--causes long-range interactions to be effectively screened, which greatly enhances the efficiency and scalability of simulations. Without imposing local or global order parameters, a linear initial configuration of the model molecule spontaneously assembles into a double helix due to the interplay between three contributions: hydrophobic/hydrophilic interactions between base pairs, backbone, and solvent; phosphate-phosphate repulsion along the backbone; and favorable base-pair stacking energy. We present results for the process of helix formation as well as for the equilibrium properties of the final state, and investigate how both depend on the input parameters. The current model holds promise for two routes of investigation: First, within a limited set of generic parameters, the effect of local (atomic-scale) perturbations on overall helical properties can be systematically studied. Second, since the efficiency allows for a direct simulation of both small and large (>100 base pairs) systems, the model presents a testground for systematic coarse-graining methods.

Protons may leak through pure lipid bilayers via a concerted mechanism.

Protons are known to permeate pure lipid bilayers at a rate that is anomalous compared to those of other small monovalent cations. The prevailing mechanism via which they cross the membrane is still unclear, and it is unknown how to probe the mechanism directly by experiment. One of the more popular theories assumes the formation of membrane-spanning single-file water wires providing a matrix along which the protons can "hop" over the barrier. However, free energy calculations on such structures (without the presence of an excess proton) suggest that this mechanism alone cannot account for the observed permeation rates. We use the multistate empirical valence bond method to directly study water structures surrounding a (delocalized) excess proton on its way through the membrane. We find that membrane-spanning networks, rather than single-file chains, are formed around the proton. We also find that such structures are considerably stabilized in the presence of the proton, with lifetimes of several hundreds of picoseconds. The observed structures are suggestive of a new, concerted, mechanism and provide some direction for further investigation.

A new perspective on the coarse-grained dynamics of fluids.

A computational methodology is presented that is designed to model, at a coarse-grained level, the mesoscale dynamics of fluids and potentially other forms of soft matter. Within a molecular dynamics simulation, "ghost" particles of a specific size, corresponding to the fundamental length-scale of coarse-graining, are used as micro-probes designed to respond to local mesoscale fluid flows and stress gradients. A subsequent coarse-grained model is then developed that incorporates both the coarse-grained mesoscale dynamics and isothermal compressibility of the original microscopic system. The method is applied to water and methanol. A contrast with dissipative particle dynamics (DPD) is also presented.