Polyalanine α-helix microsolvation: assessing the energy of the peptide desolvation penalty with density functional theory

Conformational preference of dipeptide zwitterions in aqueous solvents

Proper description of solvent effects is challenging for theoretical methods, particularly if the solute is a zwitterion. Here, a series of theoretical procedures are used to determine the preferred solvated conformations of twelve hydrophobic dipeptides (Leu-Leu, Leu-Phe, Phe-Leu, Ile-Leu, Phe-Phe, Ala-Val, Val-Ala, Ala-Ile, Ile-Ala, Ile-Val, Val-Ile and Val-Val) in the zwitterionic state. First, the accuracy of density functional theory (DFT), combined with different implicit solvent models, for describing zwitterions in aqueous solvent is assessed by comparing the predicted against the experimental glycine tautomerization energy, i.e., the energetic difference between canonical and zwitterionic glycine in aqueous solvents. It is found that among the tested solvation schemes, the charge-asymmetric nonlocally determined local-electric solvation model (CANDLE) predicts an energetic difference in excellent agreement with the experimental value. Next, DFT-CANDLE is used to determine the most favorable solvated conformation for each of the investigated dipeptide zwitterions. The CANDLE-solvated structures are obtained by exploring the conformational space of each dipeptide zwitterion concatenating DFT calculations, in vacuum, with classical molecular dynamics simulations, in explicit solvents, and DFT calculations including explicit water molecules. It is found that the energetically most favorable conformations are similar to those of the dipeptide zwitterions in their respective crystal structures. Such structural agreement is indicative of the DFT-CANDLE accomplishment of the description of solvated zwitterions, and suggests that these biomolecules self-assemble as quasi-rigid objects.

Effect of strand register in the stability and reactivity of crystals from peptides forming amyloid fibrils

Amyloids are cytotoxic protein aggregates that deposit in human tissues, leading to several health disorders. Their aggregates can also exhibit catalytic properties, and they have been used as candidates for the development of functional biomaterials. Despite being polymorphic, amyloids often assemble as cross-β fibrils formed by in-register β sheet layers. Recent studies of some amyloidogenic protein segments revealed that they crystallize as antiparallel out-of-register β sheets. Such arrangement has been proposed to be responsible for the cytotoxicity in amyloid diseases, however, there is still no consensus on the molecular mechanism. Interestingly, two amyloidogenic peptide segments, NFGAILS and FGAILSS, arrange into out-of-register and in-register β sheets, respectively, even though they solely differ by one aminoacid residue at both termini. In this work, we used density functional theory (DFT) to address how the strand register contributes into the packing and molecular properties of the NFGAILS and FGAILSS crystals. Our results show that the out-of-register structure is substantially more stable, at 0 K, than the in-register one due to stronger inter-strand contacts. Based on an analysis of the electrostatic potential of the crystal slabs, it is suggested that the out-of-register may potentially interact with negatively charged groups, like those found in cell membranes. Moreover, calculated reactivity descriptors indicate a similar outcome, where only the out-of-register peptide exhibits intrinsic reactive surface sites at the exposed amine and carboxylic groups. It is therefore suggested that the out-of-register arrangement may indeed be crucial for amyloid cytotoxicity. The findings presented here could help to further our understanding of amyloid aggregation, function, and toxicity.

Origin of cooperativity in hydrogen bonding

Hydrogen bond cooperativity results from interactions between effective point dipoles which account for many-body contributions to cooperativity.