Determining suitable lego-structures to estimate stability of larger peptide nanostructures using computational methods

Structural Water Stabilizes Protein Motifs in Liquid Protein Phase: The Folding Mechanism of Short β-Sheets Coupled to Phase Transition

Macromolecular associates, such as membraneless organelles or lipid-protein assemblies, provide a hydrophobic environment, i.e., a liquid protein phase (LP), where folding preferences can be drastically altered. LP as well as the associated phase change from water (W) is an intriguing phenomenon related to numerous biological processes and also possesses potential in nanotechnological applications. However, the energetic effects of a hydrophobic yet water-containing environment on protein folding are poorly understood. Here, we focus on small β-sheets, the key motifs of proteins, undergoing structural changes in liquid–liquid phase separation (LLPS) and also model the mechanism of energy-coupled unfolding, e.g., in proteases, during W → LP transition. Due to the importance of the accurate description for hydrogen bonding patterns, the employed models were studied by using quantum mechanical calculations. The results demonstrate that unfolding is energetically less favored in LP by ~0.3–0.5 kcal·mol⁻¹ per residue in which the difference further increased by the presence of explicit structural water molecules, where the folded state was preferred by ~1.2–2.3 kcal·mol⁻¹ per residue relative to that in W. Energetics at the LP/W interfaces was also addressed by theoretical isodesmic reactions. While the models predict folded state preference in LP, the unfolding from LP to W renders the process highly favorable since the unfolded end state has >1 kcal·mol⁻¹ per residue excess stabilization.

Protein Stability and Unfolding Following Glycine Radical Formation

Glycine (Gly) residues are particularly susceptible to hydrogen abstraction; which results in the formation of the capto-dative stabilized C_α-centered Gly radical (GLR) on the protein backbone. We examined the effect of GLR formation on the structure of the Trp cage; tryptophan zipper; and the villin headpiece; three fast-folding and stable miniproteins; using all-atom (OPLS-AA) molecular dynamics simulations. Radicalization changes the conformation of the GLR residue and affects both neighboring residues but did not affect the stability of the Trp zipper. The stability of helices away from the radical center in villin were also affected by radicalization; and GLR in place of Gly15 caused the Trp cage to unfold within 1 µs. These results provide new evidence on the destabilizing effects of protein oxidation by reactive oxygen species.

A theoretical comparison of self-assembling alpha- and beta-peptide nanostructures: toward design of beta-barrel frameworks

Self-assembling peptide-based nanotubes are among the most investigated bioactive compounds as a result of their numerous potential applications as novel biomaterials. To support rational bottom-up design of such artificial nanosystems, here we investigate structural and energetic properties of various sheet-derived nanotubes. We carried out high level quantum chemical calculations on large models, composed of up to 32 amino acids, and characterized structures from extended beta-sheets to the molecular framework of beta-barrel proteins. Surprisingly, enzyme-resistant nonnatural beta-peptides have an affinity to form nanotubes that is remarkably higher than that of natural alpha-peptides. We analyzed the stability of both systems depending on (i) parallel or antiparallel orientation, (ii) the number of peptide strands, and (iii) the formed hydrogen bond pattern. Applicability is outlined by investigating guest molecules in the tubes. It is hoped that the structural and energetic data presented here will be effectively used in the design of novel peptide nanosystems.

Thermodynamic Functions of Conformational Changes: Conformational Network of Glycine Diamide Folding, Entropy Lowering, and Informational Accumulation

A refined grid of a conformational potential energy surface (PES) and a conformational entropy surface for glycine diamide was generated by ab initio molecular computations. The possible network of reaction paths was recognized in terms of the linear combinations of internal coordinates corresponding to conrotatory and disrotatory modes of motions. Such a Woodward-Hoffmann-like path selection principle was detected for the folding of this peptide from extended to some virtually cyclic structure. It seemed reasonable to assume that this principle (or its generalized form) might be applicable to protein folding. A reaction path network was projected on the potential energy, and a continuous entropy surface was constructed under the condition of reduced dimensionality. The low entropy of the folded conformation indicated an information accumulation between 326% and 1414% with respect to the fully extended or unfolded structure. It is found that the location of existing and 'latent' critical points on the surface is revealed by the extrema and inflection points of the entropy curve.

Prediction of Folding Preference of 10 kDa Silk-like Proteins Using a Lego Approach and ab Initio Calculations

Because of their great flexibility and strength resistance, both spider silks and silkworm silks are of increasing scientific and commercial interest. Despite numerous spectroscopic and theoretical studies, several structural properties at the atomic level have yet to be identified. The present theoretical investigation focuses on these issues by studying three silk-like model peptides: (AG)(64), [(AG)(4)EG](16), and [(AG)(4)PEG](16), using a Lego-type approach to construct these polypeptides. On the basis of these examples it is shown that thermoneutral isodesmic reactions and ab initio calculations provide a capable method to investigate structural properties of repetitive polypeptides. The most probable overall fold schema of these molecules with respect to the type of embedded hairpin structures were determined at the ab initio level of theory (RHF/6-311++G(d,p)//RHF/3-21G). Further on, analysis is carried out on the possible hairpin and turn regions and on their effect on the global fold. In the case of the (AG)(64) model peptide, the optimal beta-sheet/turn ratio was also determined, which provided good support for experimental observations. In addition, lateral shearing of a hairpin "folding unit" was investigated at the quantum chemical level to explain the mechanical properties of spider silk. The unique mechanical characteristics of silk bio-compounds are now investigated at the atomic level.