On cooperative effects and aggregation of GNNQQNY and NNQQNY peptides

Cooperativity and Frustration Effects (or Lack Thereof) in Polarizable and Non-polarizable Force Fields

Understanding cooperativity and frustration is crucial for studying biological processes such as molecular recognition and protein aggregation. Force fields have been extensively utilized to explore cooperativity in the formation of protein secondary structures and self-assembled systems. Multiple studies have demonstrated that polarizable force fields provide more accurate descriptions of this phenomenon compared to fixed-charge pairwise nonpolarizable force fields, thanks to the incorporation of polarization effects. In this study, we assess the performance of the AMOEBA polarizable force field and the AMBER and OPLS nonpolarizable pairwise force fields in capturing positive and negative cooperativity recently explored in neutral and charged molecular clusters using density functional theory. Our findings show that polarizable and nonpolarizable force fields qualitatively reproduce the relative cooperativity observed in electron structure calculations. However, AMBER and OPLS fail to describe absolute cooperativity. In contrast, AMOEBA accounts for the absolute cooperativity by considering interactions beyond pairwise interactions. According to the energy decomposition analysis, it is observed that the electrostatic interactions calculated with the AMBER and OPLS force fields seem to play an important and counterintuitive role in reproducing the adiabatic interaction energies calculated with density functional theory. However, it is important to note that these force fields, due to their nature, do not explicitly incorporate many-body effects, which limits their ability to accurately describe cooperativity. On the other hand, frustration in polarizable and nonpolarizable force fields is caused by changes in bond stretching and angle bending terms of the building blocks when they are forming a complex.

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.

Simulations of Amyloid-Forming Peptides in the Crystal State

There still is little treatment available for amyloid diseases, despite their significant impact on individuals and the social and economic implications for society. One reason for this is that the physical nature of amyloid formation is not understood sufficiently well. Therefore, fundamental research at the molecular level remains necessary to support the development of therapeutics. A few structures of short peptides from amyloid-forming proteins have been determined. These can in principle be used as scaffolds for designing aggregation inhibitors. Attempts to this end have often used the tools of computational chemistry, in particular molecular simulation. However, few simulation studies of these peptides in the crystal state have been presented so far. Hence, to validate the capability of common force fields (AMBER19SB, CHARMM36m, and OPLS-AA/M) to yield insight into the dynamics and structural stability of amyloid peptide aggregates, we have performed molecular dynamics simulations of twelve different peptide crystals at two different temperatures. From the simulations, we evaluate the hydrogen bonding patterns, the isotropic B-factors, the change in energy, the Ramachandran plots, and the unit cell parameters and compare the results with the crystal structures. Most crystals are stable in the simulations but for all force fields there is at least one that deviates from the experimental crystal, suggesting more work is needed on these models.

Understanding the unusual stiffness of hydrophobic dipeptide crystals

The hydrophobic diphenylalanine peptide crystal is known to be unusually stiff, with an experimental Young's modulus in the range of 19-27 GPa. Here it is shown by means of density functional theory calculations that phenylalanine-leucine, leucine-phenylalanine, alanine-valine, valine-alanine and valine-valine hydrophobic dipeptide crystals are also unusually stiff, with Young's moduli in the range of 19.7-33.3 GPa. To further our understanding of the origin of that unusual stiffness, a linear correlation is established between Young's modulus and the strength and orientation of the hydrogen bond network developed along the crystals, showing that stiffness in these materials is primarily dictated by hydrogen bonding.

Homogeneous Molecular Systems are Positively Cooperative, but Charged Molecular Systems are Negatively Cooperative

Molecular systems bound together through noncovalent interactions are ubiquitous in nature, many of which are involved in essential life processes, yet little is known about the principles governing their structure, stability, and function. Cooperativity as one of the intrinsic properties in these systems plays a key role. In this work, on the basis of our recent quantification scheme of the cooperativity effect, we present a general pattern to identify which systems are positively cooperative and which are negatively cooperative. We show that cooperativity in homogeneous molecular systems is positive, but cooperativity in charged molecular systems is negative. We also employ analytical tools from energetics and information perspectives to appreciate the origin of the cooperativity effect. We find that positive cooperativity is dominated by the exchange-correlation interaction and steric effect, whereas negative cooperativity is governed by the electrostatic interaction. Our results should have strong implications for better understanding molecular recognition, protein folding, signal transduction, allosteric regulation, and other processes.

Quantification and origin of cooperativity: insights from density functional reactivity theory

Cooperativity is a widely used chemical concept whose existence is ubiquitous in chemical and biological systems but whose quantification is still controversial and origin much less appreciated. In this work, using the interaction energy of a molecular system, which is composed of multiple copies of a building block, we propose a quantitative measurement to evaluate the cooperativity effect. This quantification approach is then applied to six molecular systems, i.e., water cluster, argon cluster, protonated water cluster, zinc atom cluster, water cluster on top of a graphene sheet, and alpha helix of glycine amino acids, each with up to 20 copies of the building block. Cooperativity is seen in all these systems. Both positive and negative cooperativity effects are observed. Employing the two energy partition schemes in density functional theory and the information-theoretic quantities such as Shannon entropy, Fisher information, information gain, etc., we then examine the origin of the cooperativity effect for these systems. Strong linear correlations between the cooperativity measure and some of these theoretical quantities have been unveiled. With these correlations, we are able to quantitatively account for their origin of cooperativity. Our results show that the interactions governing the existence and validity of the cooperativity effect are complicated. An opposite mechanism in enthalpy-entropy compensation for positive and negative cooperativity has been unveiled. These results should provide new insights and understandings from a different viewpoint about the nature and origin of cooperativity to appreciate this vastly important chemical concept.

Origin of cooperativity in hydrogen bonding

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