Intrinsic Conformational Dynamics of Alanine in Water/Ethanol Mixtures: An Experiment-Driven Molecular Dynamics Study

How do salt and lipids affect conformational dynamics of Aβ42 monomers in water?

It is well established that amyloid β-protein (Aβ) self-assembly is involved in triggering of Alzheimer's disease. On the other hand, evidence of physiological function of Aβ interacting with lipids has only begun to emerge. Details of Aβ-lipid interactions, which may underlie physiological and pathological activities of Aβ, are not well understood. Here, the effects of salt and 1,2-dimyristoyl-sn-glycero-3-phosphocholine (DMPC) lipids on conformational dynamics of Aβ42 monomer in water are examined by all-atom molecular dynamics (MD). We acquired six sets of 250 ns long MD trajectories for each of the three lipid concentrations (0, 27, and 109 mM) in the absence and presence of 150 mM salt. Ten replica trajectories per set are used to enhance sampling of Aβ42 conformational space. We show that salt facilitates long-range tertiary contacts in Aβ42, resulting in more compact Aβ42 conformations. By contrast, addition of lipids results in lipid-concentration dependent Aβ42 unfolding concomitant with enhanced stability of the turn in the A21-A30 region. At the high lipid concentration, salt enables the N-terminal region of Aβ42 to form long-range tertiary contacts and interact with lipids, which results in formation of a parallel β-strand. Aβ42 forms stable lipid-protein complexes whereby the protein is adhered to the lipid cluster rather than embedded into it. We propose that the inability of Aβ42 monomer to get embedded into the lipid cluster may be important for facilitating repair of leaks in the blood-brain barrier without penetrating and damaging cellular membranes.

Dynamical Anomaly of Aqueous Amphiphilic Solutions: Connection to Solution H-Bond Fluctuation Dynamics?

We have investigated the possible connection between "dynamical anomaly" observed in time-resolved fluorescence measurements of reactive and nonreactive solute-centered relaxation dynamics in aqueous binary mixtures of different amphiphiles and the solution intra- and interspecies H-bond fluctuation dynamics. Earlier studies have connected the anomalous thermodynamic properties of binary mixtures at very low amphiphile concentrations to the structural distortion of water. This is termed as "structural anomaly." Interestingly, the abrupt changes in the composition-dependent average rates of solute relaxation dynamics occur at amphiphile mole fractions approximately twice as large as those where structural anomalies appear. We have investigated this anomalous solution dynamical aspect by considering (water + tertiary butanol) as a model system and performed molecular dynamics simulations at several tertiary butanol (TBA) concentrations covering the extremely dilute to the moderately concentrated regimes. The "dynamical anomaly" has been followed via monitoring the composition dependence of the intra- and interspecies H-bond fluctuations and reorientational relaxations of TBA and water molecules. Solution structural aspects have been followed via examining the tetrahedral order parameter, radial and spatial distribution functions, numbers of H bonds per water and TBA molecules, and the respective populations participating in H-bond formation. Our simulations reveal abrupt changes in the H-bond fluctuations and reorientational dynamics and tetrahedral order parameter at amphiphile concentrations differing approximately by a factor of 2 and corroborates well with the steady-state and the time-resolved spectroscopic measurements. This work therefore explains, following a uniform and cogent manner, both the experimentally observed structural and dynamical anomalies in microscopic terms.

Do molecular dynamics force fields accurately model Ramachandran distributions of amino acid residues in water?

Molecular dynamics (MD) is a powerful tool for studying intrinsically disordered proteins, however, its reliability depends on the accuracy of the force field. We assess Amber ff19SB, Amber ff14SB, OPLS-AA/M, and CHARMM36m with respect to their capacity to capture intrinsic conformational dynamics of 14 guest residues x (=G, A, L, V, I, F, Y, D^P, E^P, R, C, N, S, T) in GxG peptides in water. The MD-derived Ramachandran distribution of each guest residue is used to calculate 5 J-coupling constants and amide I' band profiles to facilitate a comparison to spectroscopic data through reduced χ² functions. We show that the Gaussian model, optimized to best fit the experimental data, outperforms all MD force fields by an order of magnitude. The weaknesses of the MD force fields are: (i) insufficient variability of the polyproline II (pPII) population among the guest residues; (ii) oversampling of antiparallel at the expense of transitional β-strand region; (iii) inadequate sampling of turn-forming conformations for ionizable and polar residues; and (iv) insufficient guest residue-specificity of the Ramachandran distributions. Whereas Amber ff19SB performs worse than the other three force fields with respect to χ² values, it accounts for residue-specific pPII content better than the other three force fields. Additional testing of residue-specific RSFF1 and Amber ff14SB combined with TIP4P/2005 on six guest residues x (=A, I, F, D^P, R, S) reveals that residue specificity derived from protein coil libraries or an improved water model alone do not result in significantly lower χ² values.

Soluble State of Villin Headpiece Protein as a Tool in the Assessment of MD Force Fields

Protein self-assembly plays an important role in cellular processes. Whereas molecular dynamics (MD) represents a powerful tool in studying assembly mechanisms, its predictions depend on the accuracy of underlying force fields, which are known to overly promote protein assembly. We here examine villin headpiece domain, HP36, which remains soluble at concentrations amenable to MD studies. The experimental characterization of soluble HP36 at concentrations of 0.05 to 1 mM reveals concentration-independent 90% monomeric and 10% dimeric populations. Extensive all-atom MD simulations at two protein concentrations, 0.9 and 8.5 mM, probe the HP36 dimer population, stability, and kinetics of dimer formation within two MD force fields, Amber ff14SB and CHARMM36m. MD results demonstrate that whereas CHARMM36m captures experimental HP36 monomer populations at the lower concentration, both force fields overly promote HP36 association at the higher concentration. Moreover, contacts stabilizing HP36 dimers are force-field-dependent. CHARMM36m produces consistently higher HP36 monomer populations, lower association rates, and weaker dependence of these quantities on the protein concentration than Amber ff14SB. Nonetheless, the highest monomer populations and dissociation constants are observed when the TIP3P water model in Amber ff14SB is replaced by TIP4P/2005, showcasing the critical role of the water model in addressing the protein solubility problem in MD.