Effect of charged amino acid side chain length on lateral cross-strand interactions between carboxylate- and guanidinium-containing residues in a β-hairpin

Stability of collagen heterotrimer with same charge pattern and different charged residue identities

Salt bridges are important factors in maintaining the stability of proteins, and their contribution to protein folding has received much attention. Although the interaction energies, or stabilizing contributions, of individual salt bridges have been measured in various proteins, a systematic assessment of various types of salt bridges in a relatively uniform environment is still a valuable analysis. Here, we used a collagen heterotrimer as a host-guest platform to construct 48 heterotrimers with the same charge pattern. A variety of salt bridges were formed between the oppositely charged residues Lys, Arg, Asp, and Glu. The melting temperature (Tm) of the heterotrimers was measured with circular dichroism. The atomic structures of 10 salt bridges were shown in three x-ray crystals of heterotrimer. Molecular dynamics simulation based on the crystal structures indicated that strong, intermediate, and weak salt bridges have distinctive N-O distances. A linear regression model was used to predict the stability of heterotrimers with high accuracy (R2 = 0.93). We developed an online database to help readers understand how a salt bridge stabilizes collagen. This work will help us better understand the stabilizing mechanism of salt bridges in collagen folding and provide a new strategy to design collagen heterotrimers.

Diagonal Interactions between Glutamate and Arginine Analogs with Varying Side-Chain Lengths in a β-Hairpin

Cross-strand interactions are important for the stability of β-sheet structures. Accordingly, cross-strand diagonal interactions between glutamate and arginine analogs with varying side-chain lengths were studied in a series of β-hairpin peptides. The peptides were analyzed by homonuclear two-dimensional nuclear magnetic resonance methods. The fraction folded population and folding free energy of the peptides were derived from the chemical shift data. The fraction folded population trends could be rationalized using the strand propensity of the constituting residues, which was not the case for the peptides with lysine analogs, highlighting the difference between the arginine analogs and lysine analogs. Double-mutant cycle analysis was used to derive the diagonal ion-pairing interaction energetics. The most stabilizing diagonal cross-strand interaction was between the shortest residues (i.e., Asp2-Agp9), most likely due to the least side-chain conformational penalty for ion-pair formation. The diagonal interaction energetics in this study involving the arginine analogs appears to be consistent with and extend beyond our understanding of diagonal ion-pairing interactions involving lysine analogs. The results should be useful for designing β-strand-containing molecules to affect biological processes such as amyloid formation and protein-protein interactions.

The Effects of Charged Amino Acid Side-Chain Length on Diagonal Cross-Strand Interactions between Carboxylate- and Ammonium-Containing Residues in a β-Hairpin

The β-sheet is one of the common protein secondary structures, and the aberrant aggregation of β-sheets is implicated in various neurodegenerative diseases. Cross-strand interactions are an important determinant of β-sheet stability. Accordingly, both diagonal and lateral cross-strand interactions have been studied. Surprisingly, diagonal cross-strand ion-pairing interactions have yet to be investigated. Herein, we present a systematic study on the effects of charged amino acid side-chain length on a diagonal ion-pairing interaction between carboxylate- and ammonium-containing residues in a β-hairpin. To this end, 2D-NMR was used to investigate the conformation of the peptides. The fraction folded population and the folding free energy were derived from the chemical shift data. The fraction folded population for these peptides with potential diagonal ion pairs was mostly lower compared to the corresponding peptide with a potential lateral ion pair. The diagonal ion-pairing interaction energy was derived using double mutant cycle analysis. The Asp2-Dab9 (Asp: one methylene; Dab: two methylenes) interaction was the most stabilizing (−0.79 ± 0.14 kcal/mol), most likely representing an optimal balance between the entropic penalty to enable the ion-pairing interaction and the number of side-chain conformations that can accommodate the interaction. These results should be useful for designing β-sheet containing molecular entities for various applications.

Structural impact of thioamide incorporation into a β-hairpin

The thioamide is a naturally-occurring single atom substitution of the canonical amide bond. The exchange of oxygen to sulfur alters the amide's physical and chemical characteristics, thereby expanding its functionality. Incorporation of thioamides in prevalent secondary structures has demonstrated that they can either have stabilizing, destabilizing, or neutral effects. We performed a systematic investigation of the structural impact of thioamide incorporation in a β-hairpin scaffold with nuclear magnetic resonance (NMR). Thioamides as hydrogen bond donors did not increase the foldedness of the more stable "YKL" variant of this scaffold. In the less stable "HPT" variant of the scaffold, the thioamide could be stabilizing as a hydrogen bond donor and destabilizing as a hydrogen bond acceptor, but the extent of the perturbation depended upon the position of incorporation. To better understand these effects we performed structural modelling of the macrocyclic folded HPT variants. Finally, we compare the thioamide effects that we observe to previous studies of both side-chain and backbone perturbations to this β-hairpin scaffold to provide context for our observations.

Longer charged amino acids favor β-strand formation in hairpin peptides

Interactions between charged amino acids significantly influence the structure and function of proteins. The encoded charged amino acids Asp, Glu, Arg, and Lys have different number of hydrophobic methylenes linking the backbone to the charged functionality. It remains to be fully understood how does this difference in the number of methylenes affect protein structure stability. Protein secondary structures are the fundamental three-dimensional building blocks of protein structures. β-Sheet structures are particularly interesting, because these structures have been associated with a number of protein misfolding diseases. Herein, we report the effect of charged amino acid side chain length at two β-strand positions individually on the stability of a β-hairpin. The charged amino acids include side chains with a carboxylate, an ammonium, or a guanidinium group. The experimental peptides, fully folded reference peptides, and fully unfolded reference peptides were synthesized by solid phase peptide synthesis and analyzed by 2D NMR methods including TOCSY, DQF-COSY, and ROESY. Sequence specific assignments were performed for all peptides. The chemical shift data were used to derive the fraction folded population and the folding free energy for the experimental peptides. Results showed that the fraction folded population increased with increasing charged amino acid side chain length. These results should be useful for developing functional peptides that adopt the β-conformation.

Swapping the Positions in a Cross-Strand Lateral Ion-Pairing Interaction between Ammonium- and Carboxylate-Containing Residues in a β-Hairpin

Cross-strand lateral ion-pairing interactions are important for antiparallel β-sheet stability. Statistical studies suggested that swapping the position of cross-strand lateral residues should not significantly affect the interaction. Herein, we swapped the position of ammonium- and carboxylate-containing residues with different side-chain lengths in a cross-strand lateral ion-pairing interaction in a β-hairpin. The peptides were analyzed by 2D-NMR. The fraction folded population and folding free energy were derived from the chemical shift data. The ion-pairing interaction energy was derived using double mutant cycle analysis. The general trends for the fraction folded population and interaction energetics remained similar upon swapping the position of the interacting charged residues. The most stabilizing cross-strand interactions were between short residues, similar to the unswapped study. However, the fraction folded populations for most of the swapped peptides were higher compared to the corresponding unswapped peptides. Furthermore, subtle differences in the ion-pairing interaction energy upon swapping were observed, most likely due to the “unleveled” relative positioning of the interacting residues created by the inherent right-handed twist of the structure. These results should be useful for developing functional peptides that rely on lateral ion-pairing interactions across antiparallel β-strands.

Influence of Calcium Binding on Conformations and Motions of Anionic Polyamino Acids. Effect of Side Chain Length

Investigation of the effect of CaCl2 salt on conformations of two anionic poly(amino acids) with different side chain lengths, poly-(α-l glutamic acid) (PGA) and poly-(α-l aspartic acid) (PASA), was performed by atomistic molecular dynamics (MD) simulations. The simulations were performed using both unbiased MD and the Hamiltonian replica exchange (HRE) method. The results show that at low CaCl2 concentration adsorption of Ca2+ ions lead to a significant chain size reduction for both PGA and PASA. With the increase in concentration, the chains sizes partially recover due to electrostatic repulsion between the adsorbed Ca2+ ions. Here, the side chain length becomes important. Due to the longer side chain and its ability to distance the charged groups with adsorbed ions from both each other and the backbone, PGA remains longer in the collapsed state as the CaCl2 concentration is increased. The analysis of the distribution of the mineral ions suggests that both poly(amino acids) should induce the formation of mineral with the same structure of the crystal cell.

Effects of Arginine Deimination and Citrulline Side-Chain Length on Peptide Secondary Structure Formation

Post-translational modifications expand the chemical functionality of peptides and proteins beyond that originating from the encoded amino acids, but studies on the structural effects of these modifications have been limited. Arginine undergoes deimination to give citrulline (Cit), converting the positively charged guanidinium moiety into a neutral urea group. Herein, we report the effect of Arg deimination on secondary structure formation. To understand the reason for the number of methylene units in Cit, the effect of Cit side-chain length on secondary structure formation was also studied. Ala-based peptides and β-hairpin peptides were used to study α-helix and β-sheet formation, respectively. Peptides containing Cit analogues were prepared by an orthogonal protecting group strategy coupled with solid-phase carbamylation. The CD data for the Ala-based peptides were analyzed by using modified Lifson-Roig theory, showing that the helix propensity of Arg decreased upon deimination and that either shortening or lengthening Cit also decreased the helix propensity. The β-hairpin peptides were analyzed by NMR methods, showing minimal change in strand formation energetics upon Arg deimination. Altering the Cit side-chain length did not affect strand formation energetics either. These results should be useful for the preparation of urea-bearing systems and the design of peptides incorporating urea-bearing residues with varying side-chain length.

Domain Mapping of Heat Shock Protein 70 Reveals That Glutamic Acid 446 and Arginine 447 Are Critical for Regulating Superoxide Dismutase 2 Function

Stress-inducible heat shock protein 70 (hsp70) interacts with superoxide dismutase 2 (SOD2) in the cytosol after synthesis to transfer the enzyme to the mitochondria for subsequent activation. However, the structural basis for this interaction remains to be defined. To map the SOD2-binding site in hsp70, mutants of hsp70 were made and tested for their ability to bind SOD2. These studies showed that SOD2 binds in the amino acid 393-537 region of the chaperone. To map the hsp70-binding site in SOD2, we used a series of pulldown assays and showed that hsp70 binds to the amino-terminal domain of SOD2. To better define the binding site, we used a series of decoy peptides derived from the primary amino acid sequence in the SOD2-binding site in hsp70. This study shows that SOD2 specifically binds to hsp70 at ⁴⁴⁵GERAMT⁴⁵⁰ Small peptides containing GERAMT inhibited the transfer of SOD2 to the mitochondria and decreased SOD2 activity in vitro and in vivo To determine the amino acid residues in hsp70 that are critical for SOD2 interactions, we substituted each amino acid residue for alanine or more conservative residues, glutamine or asparagine, in the GERAMT-binding site. Substitutions of E446A/Q and R447A/Q inhibited the ability of the GERAMT peptide to bind SOD2 and preserved SOD2 function more than other substitutions. Together, these findings indicate that the GERAMT sequence is critical for hsp70-mediated regulation of SOD2 and that Glu⁴⁴⁶ and Arg⁴⁴⁷ cooperate with other amino acid residues in the GERAMT-binding site for proper chaperone-dependent regulation of SOD2 antioxidant function.

Coupling of the non-amyloid-component (NAC) domain and the KTK(E/Q)GV repeats stabilize the α-synuclein fibrils

The aggregates of α-synuclein (αS) are a major pathological hallmark of Parkinson's disease (PD) making their structure-function relationship important for rational drug design. Yet, the atomic structure of the αS aggregates is unavailable, making it difficult to understand the underlying aggregation mechanism. In this work, based on available experimental data, we examined plausible molecular structures of αS(20/30-110) fibrils for the first time by employing computational approaches. The optimized structure was used to investigate possible interactions with aggregation inhibitors. Our structural models characterize the essential properties of the five-layered fold of the αS fibril. The distribution of the β-strands and the topology of the five β-strands in the relatively stable models are in good agreement with experimental values. In particular, we find that the KTK(E/Q)GV repeat motifs significantly stabilize the αS fibrils. The charged residues within each repeat prefer exposure to the solvent in order to further stabilize the inter-layered interactions by salt-bridges. The organization of the repeat K(58)T(59)K(60)E(61)Q(62)V(63) between the β2 and β3 layers significantly affects the stability of the non-amyloid-component (NAC) domain. The coupling between the NAC domain and the KTKEGV repeats indicates that both regions can be potential binding sites for inhibitor design. The distinct binding modes of chemical agents that alter αS aggregation highlight the potential of our models in inhibitor design.