Molecular dynamics study of water channels in natural and synthetic amyloid-ß fibrils.

We compared all-atom explicit solvent molecular dynamics simulations of three types of Aß(1-40) fibrils: brain-seeded fibrils (2M4J, with a threefold axial symmetry) and the other two, all-synthetic fibril polymorphs (2LMN and 2LMP, made under different fibrillization conditions). Fibril models were constructed using either a finite or an infinite number of layers made using periodic images. These studies yielded four conclusions. First, finite fibrils tend to unravel in a manner reminiscent of fibril dissolution, while infinite fibrils were more stable during simulations. Second, salt bridges in these fibrils remained stable in those fibrils that contained them initially, and those without salt bridges did not develop them over the time course of the simulations. Third, all fibrils tended to develop a "stagger" or register shift of ß-strands along the fibril axis. Fourth and most importantly, the brain-seeded, 2M4J, infinite fibrils allowed bidirectional transport of water in and out of the central longitudinal core of the fibril by rapidly developing gaps at the fibril vertices. 2LMP fibrils also showed this behavior, although to a lesser extent. The diffusion of water molecules in the fibril core region involved two dynamical states: a localized state and directed diffusion in the presence of obstacles. These observations provided support for the hypothesis that Aß fibrils could act as nanotubes. At least some Aß oligomers resembled fibrils structurally in having parallel, in-register ß-sheets and a sheet-turn-sheet motif. Thus, our findings could have implications for Aß cytotoxicity, which may occur through the ability of oligomers to form abnormal water and ion channels in cell membranes.

Ionic strength independence of charge distributions in solvation of biomolecules.

Electrostatic forces enormously impact the structure, interactions, and function of biomolecules. We perform all-atom molecular dynamics simulations for 5 proteins and 5 RNAs to determine the dependence on ionic strength of the ion and water charge distributions surrounding the biomolecules, as well as the contributions of ions to the electrostatic free energy of interaction between the biomolecule and the surrounding salt solution (for a total of 40 different biomolecule/solvent combinations). Although water provides the dominant contribution to the charge density distribution and to the electrostatic potential even in 1M NaCl solutions, the contributions of water molecules and of ions to the total electrostatic interaction free energy with the solvated biomolecule are comparable. The electrostatic biomolecule/solvent interaction energies and the total charge distribution exhibit a remarkable insensitivity to salt concentrations over a huge range of salt concentrations (20 mM to 1M NaCl). The electrostatic potentials near the biomolecule's surface obtained from the MD simulations differ markedly, as expected, from the potentials predicted by continuum dielectric models, even though the total electrostatic interaction free energies are within 11% of each other.

Long-time dynamics of Met-enkephalin: comparison of theory with Brownian dynamics simulations.

A recent theory for the long time dynamics of flexible chain molecules is applied for the first time to a peptide of biological importance, the neurotransmitter met-enkephalin. The dynamics of met-enkephalin is considerably more complicated than that of the previously studied glycine oligomers; met-enkephalin contains the interesting motions of phenyl groups and of side chains relative to the backbone, motions that are present in general flexible peptides. The theory extends the generalized Rouse (GR) model used to study the dynamics of polymers by providing a systematic procedure for including the contributions from the memory function matrices neglected in the GR theory. The new method describes the dynamics by time correlation functions instead of individual trajectories. These correlation functions are analytically expressed in terms of a set of equilibrium averages and the eigenvalues and eigenfunctions of the diffusion operator. The predictions of the theory are compared with Brownian dynamics (BD) simulations, so that both theory and simulation use identical potential functions and solvent models. The theory thus contains no adjustable parameters. Inclusion of the memory function contributions profoundly affects the dynamics. The theory produces very good agreement with the BD simulations for the global motions of met-enkephalin. It also correctly predicts the long-time relaxation rate for local motions.