Molecular Simulation of the Binding of Amyloid Beta to Apolipoprotein A-I in High-Density Lipoproteins

Disrupted clearance of amyloid beta (Aβ) from the brain enhances its aggregation and formation of amyloid plaques in Alzheimer's disease. The most abundant protein constituent of circulating high-density lipoprotein (HDL) particles, apoA-I, readily crosses the blood-brain barrier from periphery circulation, exhibits low-micromolar binding affinity for soluble, neurotoxic forms of Aβ, and modulates Aβ aggregation and toxicity in vitro. Its highly conserved N-terminal sequence, 42LNLKLLD48 ('LN'), has been proposed as a binding region for Aβ. However, high-resolution structural characterisation of the mechanism of HDL-Aβ interaction is very difficult to attain. Molecular dynamics simulations were conducted to investigate for the first time the interaction of Aβ and the 'LN' segment of apoA-I. Favourable binding of Aβ by HDLs was found to be driven by hydrophobic and hydrogen-bonding interactions predominantly between the 'LN' segment of apoA-I and Aβ. Preferential binding of Aβ may proceed in small, protein-rich HDLs whereby solvent-exposed hydrophobic 'LN' segments of apoA-I interact specifically with Aβ, stabilising it on the HDL surface in a possibly non-amyloidogenic conformation, facilitating effective Aβ clearance. These findings rationalise the potentially therapeutic role of HDLs in reducing Aβ aggregation and toxicity, and of peptide mimics of the apoA-I interacting region in blocking Aβ aggregation.

Characterization of the Conformations of Amyloid Beta 42 in Solution That May Mediate Its Initial Hydrophobic Aggregation

Intrinsically disordered peptides, such as amyloid β42 (Aβ42), lack a well-defined structure in solution. Aβ42 can undergo abnormal aggregation and amyloidogenesis in the brain, forming fibrillar plaques, a hallmark of Alzheimer's disease. The insoluble fibrillar forms of Aβ42 exhibit well-defined, cross β-sheet structures at the molecular level and are less toxic than the soluble, intermediate disordered oligomeric forms. However, the mechanism of initial interaction of monomers and subsequent oligomerization is not well understood. The structural disorder of Aβ42 adds to the challenges of determining the structural properties of its monomers, making it difficult to understand the underlying molecular mechanism of pathogenic aggregation. Certain regions of Aβ42 are known to exhibit helical propensity in different physiological conditions. NMR spectroscopy has shown that the Aβ42 monomer at lower pH can adopt an α-helical conformation and as the pH is increased, the peptide switches to β-sheet conformation and aggregation occurs. CD spectroscopy studies of aggregation have shown the presence of an initial spike in the amount of α-helical content at the start of aggregation. Such an increase in α-helical content suggests a mechanism wherein the peptide can expose critical non-polar residues for interaction, leading to hydrophobic aggregation with other interacting peptides. We have used molecular dynamics simulations to characterize in detail the conformational landscape of monomeric Aβ42 in solution to identify molecular properties that may mediate the early stages of oligomerization. We hypothesized that conformations with α-helical structure have a higher probability of initiating aggregation because they increase the hydrophobicity of the peptide. Although random coil conformations were found to be the most dominant, as expected, α-helical conformations are thermodynamically accessible, more so than β-sheet conformations. Importantly, for the first time α-helical conformations are observed to increase the exposure of aromatic and hydrophobic residues to the aqueous solvent, favoring their hydrophobically driven interaction with other monomers to initiate aggregation. These findings constitute a first step toward characterizing the mechanism of formation of disordered, low-order oligomers of Aβ42.

Association and Dissociation of an Aqueous Amphiphile at Elevated Temperatures

The hydrophobic interaction is often thought to increase with increasing temperature. Although there is good experimental evidence for decreased aqueous solubility and increased clustering of both nonpolar and amphiphilic molecules as temperature is increased, the detailed nature of the changes in intermolecular interactions with temperature remain unknown. By use of isotope substitution neutron scattering difference measurements on a 0.04 mole fraction solution of tert-butanol in water as the solute clustering passes through a temperature maximum, the changes in local intermolecular structures are examined. Although, as expected, the solute molecules cluster through increased contact between their nonpolar head groups with the exclusion of water, the detailed geometry of the mutual interactions changes as temperature increases. As the clustering breaks up with further temperature increase, the local structures formed do not mirror those that were found in the low-temperature dispersed system: the disassembly process is not the reverse of assembly. The clusters formed by the solute head groups are reminiscent of structures that are found in systems of spherical molecules, modulated by the additional constraint of near-maximal hydrogen bonding between the polar tails of the alcohol and the solvent water. Although the overall temperature behavior is qualitatively what would be expected of a hydrophobically driven system, the way the system resolves the competing interactions and their different temperature dependencies is complex, suggesting it could be misleading to think of the aggregation of aqueous amphiphiles solely in terms of a hydrophobic driving force.

Low-temperature and high-pressure induced swelling of a hydrophobic polymer-chain in aqueous solution

We report molecular dynamics simulations of a hydrophobic polymer-chain in aqueous solution between 260 K and 420 K at pressures of 1 bar, 3000 bar, and 4500 bar. The simulations reveal a hydrophobically collapsed structure at low pressures and high temperatures. At 3000 bar and about 260 K and at 4500 bar and about 260 K, however, an abrupt transition to a swelled state is observed. The transition is driven by a smaller volume and a remarkably strong lower enthalpy of the swelled state, indicating a steep positive slope of the corresponding transition line. The swelling is strongly stabilized by the energetically favorable state of water in the polymer's hydrophobic first hydration shell at low temperatures. This finding is consistent with the observation of a positive heat capacity of hydrophobic solvation. Moreover, the slope and location of the estimated swelling transition line for the collapsed hydrophobic chain coincides remarkably well with the cold denaturation transition of proteins.

Experimental configurational landscapes in aqueous solutions

Structures and interactions between molecules in solution are modulated by the solvent. Changes in solvent conditions can lead to structural changes and transitions such as the assembly processes seen in micelle formation and protein folding. In the case of even quite complex liquid systems, we can now explore experimentally the configurational energy landscapes that underlie these processes. Using an aqueous solution of an amphiphile as an example, the structural transitions induced by changes in temperature, concentration and added salt are examined at the molecular level, and some critical regions of the landscape identified. Moreover, the potentials of mean force that quantitatively describe the solvent-modulated interaction between molecules in solution can now be experimentally accessed.

Heat capacity effects associated with the hydrophobic hydration and interaction of simple solutes: A detailed structural and energetical analysis based on molecular dynamics simulations

We examine the SPCE [H. J. C. Berendsen et al., J. Chem. Phys. 91, 6269 (1987)] and TIP5P [M. W. Mahoney and W. L. Jorgensen, J. Chem. Phys 112, 8910 (2000)] water models using a temperature series of molecular-dynamics simulations in order to study heat-capacity effects associated with the hydrophobic hydration and interaction of xenon particles. The temperature interval between 275 and 375 K along the 0.1-MPa isobar is studied. For all investigated models and state points we calculate the excess chemical potential for xenon employing the Widom particle insertion technique. The solvation enthalpy and excess heat capacity is obtained from the temperature dependence of the chemical potentials and, alternatively, directly by Ewald summation, as well as a reaction field based method. All three methods provide consistent results. In addition, the reaction field technique allows a separation of the solvation enthalpy into solute/solvent and solvent/solvent parts. We find that the solvent/solvent contribution to the excess heat capacity is dominating, being about one order of magnitude larger than the solute/solvent part. This observation is attributed to the enlarged heat capacity of the water molecules in the hydration shell. A detailed spatial analysis of the heat capacity of the water molecules around a pair of xenon particles at different separations reveals that even more enhanced heat capacity of the water located in the bisector plane between two adjacent xenon atoms is responsible for the maximum of the heat capacity found for the desolvation barrier distance, recently reported by Shimizu and Chan [J. Am. Chem. Soc. 123, 2083 (2001)]. The about 60% enlarged heat capacity of water in the concave part of the joint xenon-xenon hydration shell is the result of a counterplay of strengthened hydrogen bonds and an enhanced breaking of hydrogen bonds with increasing temperature. Differences between the two models with respect to the heat capacity in the xenon-xenon contact state are attributed to the different water model bulk heat capacities, and to the different spatial extension of the structure effect introduced by the hydrophobic particles. Similarities between the different states of water in the joint xenon-xenon hydration shell and the properties of stretched water are discussed.

Temperature dependence of the hydrophobic hydration and interaction of simple solutes: An examination of five popular water models

We examine five different popular rigid water models (SPC, SPCE, TIP3P, TIP4P, and TIP5P) using molecular dynamics simulations in order to investigate the hydrophobic hydration and interaction of apolar Lennard-Jones solutes as a function of temperature in the range between 275 and 375 K along the 0.1 MPa isobar. For all investigated models and state points we calculate the excess chemical potential for the noble gases and methane employing the Widom particle insertion technique. All water models exhibit too small hydration entropies, but show a clear hierarchy. TIP3P shows poorest agreement with experiment, whereas TIP5P is closest to the experimental data at lower temperatures and SPCE is closest at higher temperatures. As a first approximation, this behavior can be rationalized as a temperature shift with respect to the solvation behavior found in real water. A rescaling procedure inspired by the information theory model of Hummer et al. [Chem. Phys. 258, 349 (2000)] suggests that the different solubility curves for the different models and real water can be largely explained on the basis of the different density curves at constant pressure. In addition, the models that give a good representation of the water structure at ambient conditions (TIP5P, SPCE, and TIP4P) show considerably better agreement with the experimental data than the ones which exhibit less structured O-O correlation functions (SPC and TIP3P). In the second part of the paper we calculate the hydrophobic interaction between xenon particles directly from a series of 60 ns simulation runs. We find that the temperature dependence of the association is to a large extent related to the strength of the solvation entropy. Nevertheless, differences between the models seem to require a more detailed molecular picture. The TIP5P model shows by far the strongest temperature dependence. The suggested density rescaling is also applied to the chemical potential in the xenon-xenon contact-pair configuration, indicating the presence of a temperature where the hydrophobic interaction turns into purely repulsive. The predicted association for xenon in real water suggests the presence of a strong variation with temperature, comparable to the behavior found for TIP5P water. Comparing different water models and experimental data we conclude that a proper description of density effects is an important requirement for a water model to account correctly for the correct description of the hydrophobic effects. A water model exhibiting a density maximum at the correct temperature is desirable.

Molecular theory of hydrophobic effects: "She is too mean to have her name repeated."

This paper reviews the molecular theory of hydrophobic effects relevant to biomolecular structure and assembly in aqueous solution. Recent progress has resulted in simple, validated molecular statistical thermodynamic theories and clarification of confusing theories of decades ago. Current work is resolving effects of wider variations of thermodynamic state, e.g., pressure denaturation of soluble proteins, and more exotic questions such as effects of surface chemistry in treating stability of macromolecular structures in aqueous solution.

Solvation thermodynamics of cyclohexane

The solvation thermodynamics of cyclohexane in pure liquid phase and in water is analyzed by means of the theoretical approach developed by Lee. The sum of the work of cavity creation and the dispersive solute–solvent interaction energy reproduces well the experimental Gibbs energy values over the whole temperature range 5–100°C. This implies that the purely structural solvent reorganization is an exactly compensating process in both liquids. The dispersive solute–solvent interaction energy is larger in magnitude in cyclohexane than in water, whereas the work of cavity creation is larger in water than in cyclohexane. Therefore, both terms contrast the transfer of cyclohexane from pure liquid phase to water, determining its hydrophobicity. This mechanism qualitatively corresponds to that operative in the case of benzene.Key words: hydrophobicity, cavity creation, excluded volume, dispersive interactions, enthalpy–entropy compensation.

Surface geochemistry of the clay minerals

Clay minerals are layer type aluminosilicates that figure in terrestrial biogeochemical cycles, in the buffering capacity of the oceans, and in the containment of toxic waste materials. They are also used as lubricants in petroleum extraction and as industrial catalysts for the synthesis of many organic compounds. These applications derive fundamentally from the colloidal size and permanent structural charge of clay mineral particles, which endow them with significant surface reactivity. Unraveling the surface geochemistry of hydrated clay minerals is an abiding, if difficult, topic in earth sciences research. Recent experimental and computational studies that take advantage of new methodologies and basic insights derived from the study of concentrated ionic solutions have begun to clarify the structure of electrical double layers formed on hydrated clay mineral surfaces, particularly those in the interlayer region of swelling 2:1 layer type clay minerals. One emerging trend is that the coordination of interlayer cations with water molecules and clay mineral surface oxygens is governed largely by cation size and charge, similarly to a concentrated ionic solution, but the location of structural charge within a clay layer and the existence of hydrophobic patches on its surface provide important modulations. The larger the interlayer cation, the greater the influence of clay mineral structure and hydrophobicity on the configurations of adsorbed water molecules. This picture extends readily to hydrophobic molecules adsorbed within an interlayer region, with important implications for clay-hydrocarbon interactions and the design of catalysts for organic synthesis.

Hydration and conformational equilibria of simple hydrophobic and amphiphilic solutes

We consider whether the continuum model of hydration optimized to reproduce vacuum-to-water transfer free energies simultaneously describes the hydration free energy contributions to conformational equilibria of the same solutes in water. To this end, transfer and conformational free energies of idealized hydrophobic and amphiphilic solutes in water are calculated from explicit water simulations and compared to continuum model predictions. As benchmark hydrophobic solutes, we examine the hydration of linear alkanes from methane through hexane. Amphiphilic solutes were created by adding a charge of +/-1e to a terminal methyl group of butane. We find that phenomenological continuum parameters fit to transfer free energies are significantly different from those fit to conformational free energies of our model solutes. This difference is attributed to continuum model parameters that depend on solute conformation in water, and leads to effective values for the free energy/surface area coefficient and Born radii that best describe conformational equilibrium. In light of these results, we believe that continuum models of hydration optimized to fit transfer free energies do not accurately capture the balance between hydrophobic and electrostatic contributions that determines the solute conformational state in aqueous solution.