Structural and dynamical analysis of the hydration of the Alzheimer's beta-amyloid peptide

Exploring Heterogeneous Dynamical Environment around an Ensemble of Aβ42 Peptide Monomer Conformations

Exploring the conformational properties of amyloid β (Aβ) peptides and the role of solvent (water) in guiding the dynamical environment at their interfaces is crucial for microscopic understanding of Aβ misfolding, which is involved in causing the most common neurodegenerative disorder, i.e., Alzheimer's disease. While numerous studies in the past have emphasized examining the conformational states of Aβ peptides, the role of water has not received much attention. Here, we have performed all-atom molecular dynamics simulations of several full-length Aβ₄₂ peptide monomers with different initial configurations. Our efforts are directed toward probing the origin of the heterogeneous dynamics of water around various segments of the Aβ peptide, identified as the two terminal segments (N-term and C-term) and the two hydrophobic segments (hp1 and hp2), along with the central turn region interconnecting hp1 and hp2. Our results revealed that water hydrating hp1, hp2, and turn (nonterminal segments) and C-term segments exhibit nonuniformly restricted translational as well as rotational motions. The degree of such restriction has been found to be correlated with the hydrogen bond relaxation time scales at the interface. Importantly, it is revealed that the water molecules around hp1 and, to some extent, around hp2, form relatively rigid hydration layers, compared to that around the other segments. Such rigid hydration layers arise due to relatively more solid-like caging motions resulting in relatively lesser hydration entropy. As hp1 and hp2 have been demonstrated to play a central role in Aβ aggregation, we believe that distinct water dynamics in the vicinity of these two segments, as outlined in this study, can provide vital information in understanding the early stages of the onset of the aggregation process of such peptides at higher concentration that can further aid toward advances in AD therapeutics.

How proteins modify water dynamics

Much of biology happens at the protein-water interface, so all dynamical processes in this region are of fundamental importance. Local structural fluctuations in the hydration layer can be probed by ¹⁷O magnetic relaxation dispersion (MRD), which, at high frequencies, measures the integral of a biaxial rotational time correlation function (TCF)-the integral rotational correlation time. Numerous ¹⁷O MRD studies have demonstrated that this correlation time, when averaged over the first hydration shell, is longer than in bulk water by a factor 3-5. This rotational perturbation factor (RPF) has been corroborated by molecular dynamics simulations, which can also reveal the underlying molecular mechanisms. Here, we address several outstanding problems in this area by analyzing an extensive set of molecular dynamics data, including four globular proteins and three water models. The vexed issue of polarity versus topography as the primary determinant of hydration water dynamics is resolved by establishing a protein-invariant exponential dependence of the RPF on a simple confinement index. We conclude that the previously observed correlation of the RPF with surface polarity is a secondary effect of the correlation between polarity and confinement. Water rotation interpolates between a perturbed but bulk-like collective mechanism at low confinement and an exchange-mediated orientational randomization (EMOR) mechanism at high confinement. The EMOR process, which accounts for about half of the RPF, was not recognized in previous simulation studies, where only the early part of the TCF was examined. Based on the analysis of the experimentally relevant TCF over its full time course, we compare simulated and measured RPFs, finding a 30% discrepancy attributable to force field imperfections. We also compute the full ¹⁷O MRD profile, including the low-frequency dispersion produced by buried water molecules. Computing a local RPF for each hydration shell, we find that the perturbation decays exponentially with a decay "length" of 0.3 shells and that the second and higher shells account for a mere 3% of the total perturbation measured by ¹⁷O MRD. The only long-range effect is a weak water alignment in the electric field produced by an electroneutral protein (not screened by counterions), but this effect is negligibly small for ¹⁷O MRD. By contrast, we find that the ¹⁷O TCF is significantly more sensitive to the important short-range perturbations than the other two TCFs examined here.

Effects of the dielectric properties of the ceramic-solvent interface on the binding of proteins to oxide ceramics: a non-local electrostatic approach

The rapid development of nanoscience and nanotechnology has raised many fundamental questions that significantly impede progress in these fields. In particular, understanding the physicochemical processes at the interface in aqueous solvents requires the development and application of efficient and accurate methods. In the present work we evaluate the electrostatic contribution to the energy of model protein-ceramic complex formation in an aqueous solvent. We apply a non-local (NL) electrostatic approach that accounts for the effects of the short-range structure of the solvent on the electrostatic interactions of the interfacial systems. In this approach the aqueous solvent is considered as a non-ionic liquid, with the rigid and strongly correlated dipoles of the water molecules. We have found that an ordered interfacial aqueous solvent layer at the protein- and ceramic-solvent interfaces reduces the charging energy of both the ceramic and the protein in the solvent, and significantly increases the electrostatic contribution to their association into a complex. This contribution in the presented NL approach was found to be significantly shifted with respect to the classical model at any dielectric constant value of the ceramics. This implies a significant increase of the adsorption energy in the protein-ceramic complex formation for any ceramic material. We show that for several biocompatible ceramics (for example HfO2, ZrO2, and Ta2O5) the above effect predicts electrostatically induced protein-ceramic complex formation. However, in the framework of the classical continuum electrostatic model (the aqueous solvent as a uniform dielectric medium with a high dielectric constant ∼80) the above ceramics cannot be considered as suitable for electrostatically induced complex formation. Our results also show that the protein-ceramic electrostatic interactions can be strong enough to compensate for the unfavorable desolvation effect in the process of protein-ceramic complex formation.

Non-Amyloid-β Component of Human α-Synuclein Oligomers Induces Formation of New Aβ Oligomers: Insight into the Mechanisms That Link Parkinson's and Alzheimer's Diseases

Parkinson's disease (PD) is characterized by the formation of Lewy bodies (LBs), of which their major component is the non-amyloid-β component (NAC) of α-synuclein (AS). Clinical studies have identified a link between PD and Alzheimer's disease (AD), but the question of why PD patients are at risk to develop various types of dementia, such as AD, is still elusive. In vivo studies have shown that Aβ can act as a seed for NAC/AS aggregation, promoting NAC/AS aggregation and thus contributing to the etiology of PD. However, the mechanisms by which NAC/AS oligomers interact with Aβ oligomers are still elusive. This work presents the interactions between NAC oligomers and Aβ oligomers at atomic resolution by applying extensive molecular dynamics simulations for an ensemble of cross-seeded NAC-Aβ(1-42) oligomers. The main conclusions of this study are as follows: first, the cross-seeded NAC-Aβ(1-42) oligomers represent polymorphic states, yet NAC oligomers prefer to interact with Aβ(1-42) oligomers to form double-layer over single-layer conformations due to electrostatic/hydrophobic interactions; second, among the single-layer conformations, the NAC oligomers induce formation of new β-strands in Aβ(1-42) oligomers, thus leading to new Aβ oligomer structures; and third, NAC oligomers stabilize the cross-β structure of Aβ oligomers, i.e., yielding compact Aβ fibril-like structures.

Molecular dynamics study of an insertion/duplication mutant of bacteriophage T4 lysozyme reveals the nature of α→β transition in full protein context

An α→β transition underlies the first step of disease causing amyloidogenesis in many proteins. In view of this, many studies have been carried out using peptide models to characterize these secondary structural transitions. In this paper we show that an insertion/duplication mutant 'L20' of bacteriophage T4 lysozyme (M. Sagermann, W. A. Baase and B. W. Matthews, Proc. Natl. Acad. Sci. U.S.A., 1999, 96, 6078) displays an α→β transition. We performed molecular dynamics (MD) simulation of L20, using the GROMACS package of programs and united atom GROMOS 53a6 force field for a time period of 600 ns at 300 K, in explicit water. Our MD simulation demonstrated that the transition occurs in a duplicated α-helical region inserted tandemly at the N-terminus of the 'parent' helix. We show that a C-terminal β-sheet anchors the parent helix while the loosely held N-terminal loop in the duplicate region is vulnerable to solvent attack and thus undergoes an α→β transition. Main chain-solvent interactions were seen to stabilize the observed β-structure. Thus L20 serves as a good protein model for characterization of α→β transition in a full length protein.

Investigating how peptide length and a pathogenic mutation modify the structural ensemble of amyloid beta monomer

The aggregation of amyloid beta (Aβ) peptides plays an important role in the development of Alzheimer's disease. Despite extensive effort, it has been difficult to characterize the secondary and tertiary structure of the Aβ monomer, the starting point for aggregation, due to its hydrophobicity and high aggregation propensity. Here, we employ extensive molecular dynamics simulations with atomistic protein and water models to determine structural ensembles for Aβ(42), Aβ(40), and Aβ(42)-E22K (the Italian mutant) monomers in solution. Sampling of a total of >700 microseconds in all-atom detail with explicit solvent enables us to observe the effects of peptide length and a pathogenic mutation on the disordered Aβ monomer structural ensemble. Aβ(42) and Aβ(40) have crudely similar characteristics but reducing the peptide length from 42 to 40 residues reduces β-hairpin formation near the C-terminus. The pathogenic Italian E22K mutation induces helix formation in the region of residues 20-24. This structural alteration may increase helix-helix interactions between monomers, resulting in altered mechanism and kinetics of Aβ oligomerization.

Effect of the ordered interfacial water layer in protein complex formation: A nonlocal electrostatic approach

Using a nonlocal electrostatic approach that incorporates the short-range structure of the contacting media, we evaluated the electrostatic contribution to the energy of the complex formation of two model proteins. In this study, we have demonstrated that the existence of an ordered interfacial water layer at the protein-solvent interface reduces the charging energy of the proteins in the aqueous solvent, and consequently increases the electrostatic contribution to the protein binding (change in free energy upon the complex formation of two proteins). This is in contrast with the finding of the continuum electrostatic model, which suggests that electrostatic interactions are not strong enough to compensate for the unfavorable desolvation effects.

Benchmarking implicit solvent folding simulations of the amyloid beta(10-35) fragment

A pathogenetic feature of Alzhemier disease is the aggregation of monomeric beta-amyloid proteins (Abeta) to form oligomers. Usually these oligomers of long peptides aggregate on time scales of microseconds or longer, making computational studies using atomistic molecular dynamics models prohibitively expensive and making it essential to develop computational models that are cheaper and at the same time faithful to physical features of the process. We benchmark the ability of our implicit solvent model to describe equilibrium and dynamic properties of monomeric Abeta(10-35) using all-atom Langevin dynamics (LD) simulations, since Alphabeta(10-35) is the only fragment whose monomeric properties have been measured. The accuracy of the implicit solvent model is tested by comparing its predictions with experiment and with those from a new explicit water MD simulation, (performed using CHARMM and the TIP3P water model) which is approximately 200 times slower than the implicit water simulations. The dependence on force field is investigated by running multiple trajectories for Alphabeta(10-35) using the CHARMM, OPLS-aal, and GS-AMBER94 force fields, whereas the convergence to equilibrium is tested for each force field by beginning separate trajectories from the native NMR structure, a completely stretched structure, and from unfolded initial structures. The NMR order parameter, S2, is computed for each trajectory and is compared with experimental data to assess the best choice for treating aggregates of Alphabeta. The computed order parameters vary significantly with force field. Explicit and implicit solvent simulations using the CHARMM force fields display excellent agreement with each other and once again support the accuracy of the implicit solvent model. Alphabeta(10-35) exhibits great flexibility, consistent with experiment data for the monomer in solution, while maintaining a general strand-loop-strand motif with a solvent-exposed hydrophobic patch that is believed to be important for aggregation. Finally, equilibration of the peptide structure requires an implicit solvent LD simulation as long as 30 ns.

Structures and free-energy landscapes of the wild type and mutants of the Abeta(21-30) peptide are determined by an interplay between intrapeptide electrostatic and hydrophobic interactions

The initial events in protein aggregation involve fluctuations that populate monomer conformations, which lead to oligomerization and fibril assembly. The highly populated structures, driven by a balance between hydrophobic and electrostatic interactions in the protease-resistant wild-type Abeta(21-30) peptide and mutants E22Q (Dutch), D23N (Iowa), and K28N, are analyzed using molecular dynamics simulations. Intrapeptide electrostatic interactions were connected to calculated pK(a) values that compare well with the experimental estimates. The pK(a) values of the titratable residues show that E22 and D23 side chains form salt bridges only infrequently with the K28 side chain. Contacts between E22-K28 are more probable in "dried" salt bridges, whereas D23-K28 contacts are more probable in solvated salt bridges. The strength of the intrapeptide hydrophobic interactions increases as D23N<WT<E22Q<K28A. Free-energy profiles and disconnectivity representation of the energy landscapes show that the monomer structures partition into four distinct basins. The hydrophobic interactions cluster the Abeta(21-30) peptide into two basins, differentiated by the relative position of the DVG(23-25) and GSN(25-27) fragments about the G25 residue. The E22Q mutation increases the population with intact VGSN turn compared to the wild-type (WT) peptide. The increase in the population of the structures in the aggregation-prone Basin I in E22Q, which occurs solely due to the difference in charge states between the Dutch mutant and the WT, gives a structural explanation of the somewhat larger aggregation rate in the mutant. The D23N mutation dramatically reduces the intrapeptide interactions. The K28A mutation increases the intrapeptide hydrophobic interactions that promote population of structures in Basin I and Basin II whose structures are characterized by hydrophobic interaction between V24 and K28 side chains but with well-separated ends of the backbone atoms in the VGSN turn. The intrapeptide electrostatic interactions in the WT and E22Q peptides roughen the free-energy surface compared to the K28A peptide. The D23N mutation has a flat free-energy surface, corresponding to an increased population of random coil-like structures with weak hydrophobic and electrostatic interactions. We propose that mutations or sequences that enhance the probability of occupying Basin I would promote aggregation of Abeta peptides.

Computational study on the structural diversity of amyloid Beta Peptide (abeta(10-35)) oligomers

We studied the oligomerization of Alzheimer amyloid beta peptide (Abeta) using a replica exchange molecular dynamics (REMD) simulation. The simulation was performed with Abeta(10-35) dimers, trimers, and tetramers. Extensive REMD simulations illustrated several possible oligomer conformations. As the size of the oligomer increased from a dimer to a tetramer, the number of possible configurations was reduced. We identified all the possible conformations for each oligomer and characterized their temperature dependence. It was found that the detailed structures of the oligomers, which may act as folding intermediates, are highly sensitive to the parameters of the simulation environment such as temperature and concentration. Structural diversities of Abeta oligomers suggest multiple pathways of the aggregation process.

Evaluation of the influence of the internal aqueous solvent structure on electrostatic interactions at the protein-solvent interface by nonlocal continuum electrostatic approach

The dielectric properties of the polar solvent on the protein-solvent interface at small intercharge distances are still poorly explored. To deconvolute this problem and to evaluate the pair-wise electrostatic interaction (PEI) energies of the point charges located at the protein-solvent interface we used a nonlocal (NL) electrostatic approach along with a static NL dielectric response function of water. The influence of the aqueous solvent microstructure (determined by a strong nonelectrostatic correlation effect between water dipoles within the orientational Debye polarization mode) on electrostatic interactions at the interface was studied in our work. It was shown that the PEI energies can be significantly higher than the energies evaluated by the classical (local) consideration, treating water molecules as belonging to the bulk solvent with a high dielectric constant. Our analysis points to the existence of a rather extended, effective low-dielectric interfacial water shell on the protein surface. The main dielectric properties of this shell (effective thickness together with distance- and orientation-dependent dielectric permittivity function) were evaluated. The dramatic role of this shell was demonstrated when estimating the protein association rate constants.

Diffusion and residence time of hydrogen peroxide and water in crowded protein environments

Reactive oxygen species (ROS) have important functions in cell signaling and, when present at overly high levels, may cause oxidation of important biological molecules. Kinetic models to study diffusion of ROS inside of mitochondria often assume dynamics similar to that in solution. However, it is well-known that separation of proteins in the cytosol or inside of mitochondria, where ROS are most predominant, can be smaller than 1 nm. Diffusion of small molecules can be better regarded as a percolation process. In this article, we report results of diffusivity and residence of water and hydrogen peroxide in the proximity of proteins. In carrying out this study, we found some issues with the conventional way of computing residence times by means of survival time correlation functions. The main problem is that particles remaining on the surface of a protein for long times and for which one has very poor statistics contribute significantly to the short time behavior of the survival time correlation function. We mathematically describe this problem and propose methodology to overcome it.

Synthesis, in vitro assay, and molecular modeling of new piperidine derivatives having dual inhibitory potency against acetylcholinesterase and Aβ1–42 aggregation for Alzheimer’s disease therapeutics

With the goal of developing Alzheimer's disease therapeutics, we have designed and synthesized new piperidine derivatives having dual action of acetylcholinesterase (AChE) and beta-amyloid peptide (Abeta) aggregation inhibition. For binding with the catalytic site of AChE, an ester with aromatic group was designed, and for the peripheral site, another aromatic group was considered. And for intercalating amyloid-beta oligomerization, long and linear conformation with a lipophilic group was considered. The synthetic methods employed for the structure with dual action depended on alcohols with an aromatic ring and the substituted benzoic acids, which are esterificated in the last step of the synthetic pathway. We screened these new derivatives through inhibition tests of acetylcholinesterase, butyrylcholinesterase (BChE), and Abeta(1-42) peptide aggregation, AChE-induced Abeta(1-42) aggregation. Our results displayed that compound 12 showed the best inhibitory potency and selectivity of AChE, and 29 showed the highest selectivity of BChE inhibition. Compounds 15 and 12 had inhibitory activities against Abeta(1-42) aggregation and AChE-induced Abeta aggregation. In the docking model, we confirmed that 4-chlorobenzene of 12 plays the parallel pi-pi stacking against the indole ring of Trp84 in the bottom gorge of AChE. Because the benzyhydryl moiety of 12 covered the peripheral site of AChE in a funnel-like shape, 12 showed good inhibitory potency against AChE and could inhibit AChE-induced Abeta(1-42) peptide aggregation.

Infrared Study of the Effect of Hydration on the Amide I Band and Aggregation Properties of Helical Peptides

The amide I' band of a polypeptide is sensitive not only to its secondary structure content but also to its environment. In this study we show how degrees of hydration affect the underlying spectral features of the amide I' band of two alanine-based helical peptides. This is achieved by solubilizing these peptides in the water pool of sodium bis(2-ethylhexyl)sulfosuccinate reverse micelles with different water contents or w0 values. In agreement with several earlier studies, our results show that the amide I' band arising from a group of dehydrated helical amides is centered at approximately 1650 cm-1, whereas hydration shifts this frequency toward lower wavenumbers. More importantly, temperature-dependent infrared studies further show that these helical peptides undergo a thermally induced conformational transition in reverse micelles of low w0 values (e.g., w0=6), resulting in soluble peptide aggregates rich in antiparallel beta-sheets. Interestingly, however, increasing w0 or water content leads to an increase in the onset temperature at which such beta-aggregates begin to form. Therefore, these results provide strong evidence suggesting that dehydration facilitates aggregate formation and that removal of water imposes a free energy barrier to peptide association and aggregation, a feature that has been suggested in recent simulation studies focusing on the mechanism of beta-amyloid formation.

The Structure of the Alzheimer Amyloid β 10-35 Peptide Probed through Replica-Exchange Molecular Dynamics Simulations in Explicit Solvent

The conformational states sampled by the Alzheimer amyloid beta (10-35) (Abeta 10-35) peptide were probed using replica-exchange molecular dynamics (REMD) simulations in explicit solvent. The Abeta 10-35 peptide is a fragment of the full-length Abeta 40/42 peptide that possesses many of the amyloidogenic properties of its full-length counterpart. Under physiological temperature and pressure, our simulations reveal that the Abeta 10-35 peptide does not possess a single unique folded state. Rather, this peptide exists as a mixture of collapsed globular states that remain in rapid dynamic equilibrium with each other. This conformational ensemble is dominated by random coil and bend structures with insignificant presence of an alpha-helical or beta-sheet structure. The 3D structure of Abeta 10-35 is seen to be defined by a salt bridge formed between the side-chains of K28 and D23. This salt bridge is also observed in Abeta fibrils and our simulations suggest that monomeric conformations of Abeta 10-35 contain pre-folded structural motifs that promote rapid aggregation of this peptide.

Statistical and molecular dynamics studies of buried waters in globular proteins

Buried solvent molecules are common in the core of globular proteins and contribute to structural stability. Folding necessitates the burial of polar backbone atoms in the protein core, whose hydrogen-bonding capacities should be satisfied on average. Whereas the residues in alpha-helices and beta-sheets form systematic main-chain hydrogen bonds, the residues in turns, coils and loops often contain polar atoms that fail to form intramolecular hydrogen bonds. The statistical analysis of 842 high resolution protein structures shows that well-resolved, internal water molecules preferentially reside near residues without alpha-helical and beta-sheet secondary structures. These buried waters most often form primary hydrogen bonds to main-chain atoms not involved in intramolecular hydrogen bonds, providing strong evidence that hydrating main-chain atoms is a key structural role of buried water molecules. Additionally, the average B-factor of protein atoms hydrogen-bonded to waters is smaller than that of protein atoms forming intramolecular hydrogen bonds, and the average B-factor of water molecules involved in primary hydrogen bonds with main-chain atoms is smaller than the average B-factor of water molecules involved in secondary hydrogen bonds to protein atoms that form concurrent intramolecular hydrogen bonds. To study the structural coupling between internal waters and buried polar atoms in detail we simulated the dynamics of wild-type FKBP12, in which a buried water, Wat137, forms one side-chain and multiple main-chain hydrogen bonds. We mutated E60, whose side-chain hydrogen bonds with Wat137, to Q, N, S or A, to modulate the multiplicity and geometry of hydrogen bonds to the water. Mutating E60 to a residue that is unable to form a hydrogen bond with Wat137 results in reorientation of the water molecule and leads to a structural readjustment of residues that are both near and distant to the water. We predict that the E60A mutation will result in a significantly reduced affinity of FKBP12 for its ligand FK506. The propensity of internal waters to hydrogen bond to buried polar atoms suggests that ordered water molecules may constitute fundamental structural components of proteins, particularly in regions where alpha-helical or beta-sheet secondary structure is not present.

Effect of single-point sequence alterations on the aggregation propensity of a model protein

Sequences of contemporary proteins are believed to have evolved through a process that optimized their overall fitness, including their resistance to deleterious aggregation. Biotechnological processing may expose therapeutic proteins to conditions that are much more conducive to aggregation than those encountered in a cellular environment. An important task of protein engineering is to identify alternative sequences that would protect proteins when processed at high concentrations without altering their native structure associated with specific biological function. Our computational studies exploit parallel tempering simulations of coarse-grained model proteins to demonstrate that isolated amino acid residue substitutions can result in significant changes in the aggregation resistance of the protein in a crowded environment while retaining protein structure in isolation. A thermodynamic analysis of protein clusters subject to competing processes of folding and association shows that moderate mutations can produce effects similar to those caused by changes in system conditions, including temperature, concentration, and solvent composition, that affect the aggregation propensity. The range of conditions where a protein can resist aggregation can therefore be tuned by sequence alterations, although the protein generally may retain its generic ability for aggregation.

Solvent participation in Serratia marcescens endonuclease complexes

The monomer and dimer of the bacterium Serratia marcescens endonuclease (SMnase) are each catalytically active and the two subunits of the dimer function independently of each other. Specific interfacial waters may play a role in stability, complex formation, and functionality. We performed molecular dynamics simulations of both a subunit of SMnase and its model built complex with DNA and analyzed the relation of the hydration sites to the catalytic mechanism. It was found that the binding of DNA has little influence on the global hydration properties of the protein, including occupancy and water residence time distributions. DNA and protein recognition in our model mainly involves direct contacts by hydrogen bond and hydrophobic interactions. Water-mediated contacts exist, but are less common. Three interior water clusters were identified for SMnase. One cluster around the active site in the monomer-DNA complex shows relatively strong interactions between hydration sites as well as between the sites and the biomolecules. The simulated cluster properties agreed well with experimental data. The magnesium ion shows ligand exchange. Although Mg2+ keeps six ligands during the entire simulation, upon the binding of DNA, Asn119 loses its coordination with Mg2+, while one nonbridging oxygen of the phosphate of a DNA residue and two oxygen atoms of the side chain of Glu127 become the ligands. Waters in a nearby cluster exchange and participate in the resolvation of groups in the presence of DNA. Water thus not only participates in the cleavage of DNA but also can stabilize the transition state and the leaving groups in our model.

High-resolution protein hydration NMR experiments: Probing how protein surfaces interact with water and other non-covalent ligands

High-resolution solution NMR experiments are extremely useful to characterize the location and the dynamics of hydrating water molecules at atomic resolution. However, these methods are severely limited by undesired incoherent transfer pathways such as those arising from exchange-relayed intra-molecular cross-relaxation. Here, we review several complementary exchange network editing methods that can be used in conjunction with other types of NMR hydration experiments such as magnetic relaxation dispersion and 1J(NC') measurements to circumvent these limitations. We also review several recent contributions illustrating how the original solution hydration NMR pulse sequence architecture has inspired new approaches to map other types of non-covalent interactions going well beyond the initial scope of hydration. Specifically, we will show how hydration NMR methods have evolved and have been adapted to binding site mapping, ligand screening, protein-peptide and peptide-lipid interaction profiling.

Spontaneous fibril formation by polyalanines; discontinuous molecular dynamics simulations

Fibrillary protein aggregates rich in beta-sheet structure have been implicated in the pathology of several neurodegenerative diseases. In this work, we investigate the formation of fibrils by performing discontinuous molecular dynamics simulations on systems containing 12 to 96 model Ac-KA(14)K-NH(2) peptides using our newly developed off-lattice, implicit-solvent, intermediate-resolution model, PRIME. We find that, at a low concentration, random-coil peptides assemble into alpha-helices at low temperatures. At intermediate concentrations, random-coil peptides assemble into alpha-helices at low temperatures and large beta-sheet structures at high temperatures. At high concentrations, the system forms beta-sheets over a wide range of temperatures. These assemble into fibrils above a critical temperature which decreases with concentration and exceeds the isolated peptide's folding temperature. At very high temperatures and all concentrations, the system is in a random-coil state. All of these results are in good qualitative agreement with those by Blondelle and co-workers on Ac-KA(14)K-NH(2) peptides. The fibrils observed in our simulations mimic the structural characteristics observed in experiments in terms of the number of sheets formed, the values of the intra- and intersheet separations, and the parallel peptide arrangement within each beta-sheet. Finally, we find that when the strength of the hydrophobic interaction between nonpolar side chains is high compared to the strength of hydrogen bonding, amorphous aggregates, rather than fibrillar aggregates, are formed.

Amyloid β-Peptide Oligomerization in Silico: Dimer and Trimer

Soluble oligomers of Alzheimer's amyloid beta protein (Abeta) may act as effectors of neurotoxicity in early stages of Alzheimer's disease. Detailed information about the structure of Abeta in atomistic level and the dynamics of assembly of monomeric Abeta into oligomeric structures is rather elusive. We have performed replica exchange molecular dynamics (REMD) simulations on the formation of the dimer and trimer of Abeta10-35 peptide. We have observed spontaneous formation of several basic structural units that may act as a template or an intermediate for further aggregation of Alzheimer's Abeta protein. Various conformers, including interlocking structures of experimentally known bend double beta strands, are identified.

Solvent and mutation effects on the nucleation of amyloid beta-protein folding

Experimental evidence suggests that the folding and aggregation of the amyloid beta-protein (Abeta) into oligomers is a key pathogenetic event in Alzheimer's disease. Inhibiting the pathologic folding and oligomerization of Abeta could be effective in the prevention and treatment of Alzheimer's disease. Here, using all-atom molecular dynamics simulations in explicit solvent, we probe the initial stages of folding of a decapeptide segment of Abeta, Abeta(21-30), shown experimentally to nucleate the folding process. In addition, we examine the folding of a homologous decapeptide containing an amino acid substitution linked to hereditary cerebral hemorrhage with amyloidosis-Dutch type, [Gln-22]Abeta(21-30). We find that: (i) when the decapeptide is in water, hydrophobic interactions and transient salt bridges between Lys-28 and either Glu-22 or Asp-23 are important in the formation of a loop in the Val-24-Lys-28 region of the wild-type decapeptide; (ii) in the presence of salt ions, salt bridges play a more prominent role in the stabilization of the loop; (iii) in water with a reduced density, the decapeptide forms a helix, indicating the sensitivity of folding to different aqueous environments; and (iv) the "Dutch" peptide in water, in contrast to the wild-type peptide, fails to form a long-lived Val-24-Lys-28 loop, suggesting that loop stability is a critical factor in determining whether Abeta folds into pathologic structures.

Vibrational energy relaxation in proteins

An overview of theories related to vibrational energy relaxation (VER) in proteins is presented. VER of a selected mode in cytochrome c is studied by using two theoretical approaches. One approach is the equilibrium simulation approach with quantum correction factors, and the other is the reduced model approach, which describes the protein as an ensemble of normal modes interacting through nonlinear coupling elements. Both methods result in similar estimates of the VER time (subpicoseconds) for a CD stretching mode in the protein at room temperature. The theoretical predictions are in accord with previous experimental data. A perspective on directions for the detailed study of time scales and mechanisms of VER in proteins is presented.

Kinetics of fibril formation by polyalanine peptides

Ordered beta-sheet complexes, termed amyloid fibrils, are the underlying structural components of the intra- and extracellular fibrillar protein deposits that are associated with a variety of human diseases, including Alzheimer's, Parkinson's, and the prion diseases. In this work, we investigated the kinetics of fibril formation using our newly developed off-lattice intermediate resolution model, PRIME. The model is simple enough to allow the treatment of large multichain systems while maintaining a fairly realistic description of protein dynamics without built-in bias toward any conformation when used in conjunction with constant temperature discontinuous molecular dynamics, a fast alternative to conventional molecular dynamics. Simulations were performed on systems containing 48-96 model Ac-KA14K-NH2 peptides. We found that fibril formation for polyalanines incorporate features that are characteristic of three models, the templated assembly, nucleated polymerization, and nucleated conformational conversion models, but that none of them gave a completely satisfactory description of the simulation kinetics. Fibril formation was nucleation-dependent, occurring after a lag time that decreased with increasing peptide concentration and increased with increasing temperature. Fibril formation appeared to be a conformational conversion process in which small amorphous aggregates --> beta-sheets --> ordered nucleus --> subsequent rapid growth of a small stable fibril or protofilament. Fibril growth in our simulations involved both beta-sheet elongation, in which the fibril grew by adding individual peptides to the end of each beta-sheet, and lateral addition, in which the fibril grew by adding already formed beta-sheets to its side. The initial rate of fibril formation increased with increasing concentration and decreased with increasing temperature.

Structure and function of amyloid in Alzheimer's disease

This review is focused on the structure and function of Alzheimer's amyloid deposits. Amyloid formation is a process in which normal well-folded cellular proteins undergo a self-assembly process that leads to the formation of large and ordered protein structures. Amyloid deposition, oligomerization, and higher order polymerization, and the structure adopted by these assemblies, as well as their functional relationship with cell biology are underscored. Numerous efforts have been directed to elucidate these issues and their relation with senile dementia. Significant advances made in the last decade in amyloid structure, dynamics and cell biology are summarized and discussed. The mechanism of amyloid neurotoxicity is discussed with emphasis on the Wnt signaling pathway. This review is focused on Alzheimer's amyloid fibrils in general and has been divided into two parts dealing with the structure and function of amyloid.

Phase diagrams describing fibrillization by polyalanine peptides

Amyloid fibrils are the structural components underlying the intra- and extracellular protein deposits that are associated with a variety of human diseases, including Alzheimer's, Parkinson's, and the prion diseases. In this work, we examine the thermodynamics of fibril formation using our newly-developed off-lattice intermediate-resolution protein model, PRIME. The model is simple enough to allow the treatment of large multichain systems while maintaining a fairly realistic description of protein dynamics when used in conjunction with constant-temperature discontinuous molecular dynamics, a fast alternative to conventional molecular dynamics. We conduct equilibrium simulations on systems containing 96 Ac-KA14K-NH2 peptides over a wide range of temperatures and peptide concentrations using the replica-exchange method. Based on measured values of the heat capacity, radius of gyration, and percentage of peptides that form the various structures, a phase diagram in the temperature-concentration plane is constructed delineating the regions where each structure is stable. There are four distinct single-phase regions: alpha-helices, fibrils, nonfibrillar beta-sheets, and random coils; and four two-phase regions: random coils/nonfibrillar beta-sheets, random coils/fibrils, fibrils/nonfibrillar beta-sheets, and alpha-helices/nonfibrillar beta-sheets. The alpha-helical region is at low temperature and low concentration. The nonfibrillar beta-sheet region is at intermediate temperatures and low concentrations and expands to higher temperatures as concentration is increased. The fibril region occurs at intermediate temperatures and intermediate concentrations and expands to lower as the peptide concentration is increased. The random-coil region is at high temperatures and all concentrations; this region shifts to higher temperatures as the concentration is increased.

Solvent Interactions and Protein Dynamics in Spin-labeled T4 Lysozyme

Aspects of T4 lysozyme dynamics and solvent interaction are investigated using atomically detailed Molecular Dynamics (MD) simulations. Two spin-labeled mutants of T4 lysozyme are analyzed (T4L-N40C and T4L-K48C), which have been found from electronic paramagnetic resonance (EPR) experiments to exhibit different mobilities at the site of spin probe attachment (N- and C-terminus of helix B, respectively). Similarities and differences in solvent distribution and diffusion around the spin label, as well as around exposed and buried residues within the protein, are discussed. The purpose is to capture possible strong interactions between the spin label (ring) and solvent molecules, which may affect EPR lineshapes. The effect of backbone motions on the water density profiles is also investigated. The focus is on the domain closure associated with the T4 lysozyme hinge-bending motion, which is analyzed by Essential Dynamics (ED). The N-terminus of helix B is found to be a "hinge" residue, which explains the high degree of flexibility and motional freedom at this site.

Protein-water interactions in ribonuclease A and angiogenin: A molecular dynamics study

It is known that water molecules play an important role in the biological functioning of proteins. The members of the ribonuclease A (RNase A) family of proteins, which are sequentially and structurally similar, are known to carry out the obligatory function of cleaving RNA and individually perform other diverse biological functions. Our focus is on elucidating whether the sequence and structural similarity lead to common hydration patterns, what the common hydration sites are and what the differences are. Extensive molecular dynamics simulations followed by a detailed analysis of protein-water interactions have been carried out on two members of the ribonuclease A superfamily-RNase A and angiogenin. The water residence times are analyzed and their relationship with the characteristic properties of the protein polar atoms, such as their accessible surface area and mean hydration, is studied. The capacity of the polar atoms to form hydrogen bonds with water molecules and participate in protein-water networks are investigated. The locations of such networks are identified for both proteins.

Elucidation of primary structure elements controlling early amyloid beta-protein oligomerization

Assembly of monomeric amyloid beta-protein (A beta) into oligomeric structures is an important pathogenetic feature of Alzheimer's disease. The oligomer size distributions of aggregate-free, low molecular weight A beta 40 and A beta 42 can be assessed quantitatively using the technique of photo-induced cross-linking of unmodified proteins. This approach revealed that low molecular weight A beta 40 is a mixture of monomer, dimer, trimer, and tetramer, in rapid equilibrium, whereas low molecular weight A beta 42 preferentially exists as pentamer/hexamer units (paranuclei), which self-associate to form larger oligomers. Here, photo-induced cross-linking of unmodified proteins was used to evaluate systematically the oligomerization of 34 physiologically relevant A beta alloforms, including those containing familial Alzheimer's disease-linked amino acid substitutions, naturally occurring N-terminal truncations, and modifications altering the charge, the hydrophobicity, or the conformation of the peptide. The most important structural feature controlling early oligomerization was the length of the C terminus. Specifically, the side-chain of residue 41 in A beta 42 was important both for effective formation of paranuclei and for self-association of paranuclei into larger oligomers. The side-chain of residue 42, and the C-terminal carboxyl group, affected paranucleus self-association. A beta 40 oligomerization was particularly sensitive to substitutions of Glu22 or Asp23 and to truncation of the N terminus, but not to substitutions of Phe19 or Ala21. A beta 42 oligomerization, in contrast, was largely unaffected by substitutions at positions 22 or 23 or by N-terminal truncations, but was affected significantly by substitutions of Phe19 or Ala21. These results reveal how specific regions and residues control A beta oligomerization and show that these controlling elements differ between A beta 40 and A beta 42.