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

Intrinsic Conformational Dynamics of Glycine and Alanine in Polarizable Molecular Dynamics Force Fields: Comparison to Spectroscopic Data

Molecular dynamics (MD) is a great tool for elucidating conformational dynamics of proteins and peptides in water at the atomistic level that often surpasses the level of detail available experimentally. Structure predictions, however, are limited by the accuracy of the underlying MD force field. This limitation is particularly stark in the case of intrinsically disordered peptides and proteins, which are characterized by solvent-accessible and disordered peptide regions and domains. Recent studies show that most additive MD force fields, including CHARMM36m, do not reproduce the intrinsic conformational distributions of guest amino acid residues x in cationic GxG peptides in water in line with experimental data. Positing that a lack of polarizability in additive MD force fields may be the culprit for the reported discrepancies, we here examine the conformational dynamics of guest glycine and alanine residues in cationic GxG peptides in water using two polarizable MD force fields, CHARMM Drude and AMOEBA. Our results indicate that while AMOEBA captures the experimental data better than CHARMM Drude, neither of the two polarizable force fields offers an improvement of the Ramachandran distributions of glycine and alanine residues in cationic GGG and GAG peptides, respectively, over CHARMM36m.

Conformational Manifold Sampled by Two Short Linear Motif Segments Probed by Circular Dichroism, Vibrational, and Nuclear Magnetic Resonance Spectroscopy

Disordered protein segments called short linear motifs (SLiM) serve as recognition sites for a variety of biological processes and act as targeting signals, modification, and ligand binding sites. While SLiMs do not adopt one of the known regular secondary structures, the conformational distribution might still reflect the structural propensities of their amino acid residues and possible interactions between them. In the past, conformational analyses of short peptides provided compelling evidence for the notion that individual residues are less conformationally flexible than locally expected for a random coil. Here, we combined various spectroscopies (NMR, IR, vibrational, and UV circular dichroism) to determine the Ramachandran plots of two SLiM motifs, i.e., GRRDSG and GRRTSG. They are two representatives of RxxS motifs that are capable of being phosphorylated by protein kinase A, an enzyme that plays a fundamental role in a variety of biological processes. Our results reveal that the nearest and non-nearest interactions between residues cause redistributions between polyproline II and β-strand basins while concomitantly stabilizing extended relative to turn-forming and helical structures. They also cause shifts in basin positions. With increasing temperature, β-strand populations become more populated at the expense of polyproline II. While molecular dynamics simulations with Amber ff14SB and CHARMM 36m force fields indicate residue-residue interactions, they do not account for the observed structural changes.

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.

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.