Microsolvation Effect on the Twist of β-sheets

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.

Crocin Inhibits the Fibrillation of Human α-synuclein and Disassembles Mature Fibrils: Experimental Findings and Mechanistic Insights from Molecular Dynamics Simulation

The aggregation of human alpha-synuclein (hαS) is pivotally implicated in the development of most types of synucleinopathies. Molecules that can inhibit or reverse the aggregation process of amyloidogenic proteins have potential therapeutic value. The anti-aggregating activity of multiple carotenoid compounds has been reported over the past decades against a growing list of amyloidogenic polypeptides. Here, we aimed to determine whether crocin, the main carotenoid glycoside component of saffron, would inhibit hαS aggregation or could disassemble its preformed fibrils. By employing a series of biochemical and biophysical techniques, crocin was exhibited to inhibit hαS fibrillation in a dose-dependent fashion by stabilizing very early aggregation intermediates in off-pathway non-toxic conformations with little β-sheet content. We also observed that crocin at high concentrations could efficiently destabilize mature fibrils and disassemble them into seeding-incompetent intermediates by altering their β-sheet conformation and reshaping their structure. Our atomistic molecular dynamics (MD) simulations demonstrated that crocin molecules bind to both the non amyloid-β component (NAC) region and C-terminal domain of hαS. These interactions could thereby stabilize the autoinhibitory conformation of the protein and prevent it from adopting aggregation-prone structures. MD simulations further suggested that ligand molecules prefer to reside longitudinally along the fibril axis onto the edges of the inter-protofilament interface where they establish hydrogen and hydrophobic bonds with steric zipper stabilizing residues. These interactions turned out to destabilize hαS fibrils by altering the interstrand twist angles, increasing the rigidity of the fibril core, and elevating its radius of gyration. Our findings suggest the potential pharmaceutical implication of crocin in synucleinopathies.

PSIQUE: Protein Secondary Structure Identification on the Basis of Quaternions and Electronic Structure Calculations

The secondary structure is important in protein structure analysis, classification, and modeling. We have developed a novel method for secondary structure assignment, termed PSIQUE, based on the potential energy surface (PES) of polyalanine obtained using an infinitely long chain model and density functional theory calculations. First, uniform protein segments are determined in terms of a difference of quaternions between neighboring amino acids along the protein backbone. Then, the identification of the secondary structure motifs is carried out based on the minima found in the PES. PSIQUE shows good agreement with other secondary structure assignment methods. However, it provides better discrimination of subtle secondary structures (e.g., helix types) and termini and produces more uniform segments while also accounting for local distortions. Overall, PSIQUE provides a precise and reliable assignment of secondary structures, so it should be helpful for the detailed characterization of the protein structure.

Origin of cooperativity in hydrogen bonding

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

The folding of acetyl(Ala)28NH2 and acetyl(Ala)40NH2 extended strand peptides into antiparallel β-sheets. A density functional theory study of β-sheets with β-turns

We report ONIOM calculations using B3LYP/D95** and AM1 on β-sheet formation from acetyl(Ala)(N)NH(2) (N = 28 or 40). The sheets contain from one to four β-turns for N = 28 and up to six for N = 40. We have obtained four types of geometrically optimized structures. All contain only β-turns. They differ from each other in the types of β-turns formed. The unsolvated sheets containing two turns are most stable. Aqueous solvation (using the SM5.2 and CPCM methods) reduces the stabilities of the folded structures compared to the extended strands.

Aqueous solvation of polyalanine α-helices with specific water molecules and with the CPCM and SM5.2 aqueous continuum models using density functional theory

We present density functional theory (DFT) calculations at the X3LYP/D95(d,p) level on the solvation of polyalanine α-helices in water. The study includes the effects of discrete water molecules and the CPCM and AMSOL SM5.2 solvent continuum model both separately and in combination. We find that individual water molecules cooperatively hydrogen-bond to both the C- and N-termini of the helix, which results in increases in the dipole moment of the helix/water complex to more than the vector sum of their individual dipole moments. These waters are found to be more stable than in bulk solvent. On the other hand, individual water molecules that interact with the backbone lower the dipole moment of the helix/water complex to below that of the helix itself. Small clusters of waters at the termini increase the dipole moments of the helix/water aggregates, but the effect diminishes as more waters are added. We discuss the somewhat complex behavior of the helix with the discrete waters in the continuum models.