Amino acid pairing at the N- and C-termini of helical segments in proteins

Structure-based evaluation of the envelope domain III-nonstructural protein 1 (EDIII-NS1) fusion as a dengue virus vaccine candidate

The lack of effective medicines or vaccines, combined with climate change and other environmental factors, annually subjects a significant proportion of the world's inhabitants to the risk of dengue virus (DENV) infection. These conditions increase the likelihood of exposure to mosquito-borne diseases such as dengue fever. Hence, many research approaches tend to develop efficient vaccine candidates against the dengue virus. Therefore, we used immunoinformatics and bioinformatics to design a construction for developing a candidate vaccine against dengue virus serotypes. In this study, the in silico structure, containing the non-structural protein 1 region (NS1) (consensus and epitope), the envelope domain III protein (EDIII) as the structural part of the virus construction, and the bc-loop of envelope domain II (EDII) as the neutralizing and protected epitope, were employed. We utilized in silico tools to enhance the immunogenicity and effectiveness of dengue virus vaccine candidates. Evaluations included refining and validating physicochemical characteristics, B and T-cell epitopes, homology modeling, and the three-dimensional structure to assess the designed vaccine's quality. In silico results for tertiary structure prediction and validation revealed high-quality modeling for all vaccine constructs. Additionally, the instructed model demonstrated stability throughout molecular dynamics simulation. The results of the immune simulation suggested that the titers of IgG and IgM could be raised to desirable values following injection into in vivo models. It can be concluded that the designed construct effectively induce humoral and cellular immunity and can be proposed as effective vaccine candidate against four dengue serotypes.Communicated by Ramaswamy H. Sarma.

Propensities of Some Amino Acid Pairings in α-Helices Vary with Length

The results of secondary structure prediction methods are widely used in applications in biotechnology and bioinformatics. However, the accuracy limit of these methods could be improved up to 92%. One approach to achieve this goal is to harvest information from the primary structure of the peptide. This study aims to contribute to this goal by investigating the variations in propensity of amino acid pairings to α-helices in globular proteins depending on helix length. (n):(n + 4) residue pairings were determined using a comprehensive peptide data set according to backbone hydrogen bond criterion which states that backbone hydrogen bond is the dominant driving force of protein folding. Helix length is limited to 13 to 26 residues. Findings of this study show that propensities of ALA:GLY and GLY:GLU pairings to α-helix in globular protein increase with increasing helix length but of ALA:ALA and ALA:VAL decrease. While the frequencies of ILE:ALA, LEU:ALA, LEU:GLN, LEU:GLU, LEU:LEU, MET:ILE and VAL:LEU pairings remain roughly constant with length, the 25 residue pairings have varying propensities in narrow helix lengths. The remaining pairings have no prominent propensity to α-helices.

Propensities of Amino Acid Pairings in Secondary Structure of Globular Proteins

A class of secondary structure prediction algorithms use the information from the statistics of the residue pairs found in secondary structural elements. Because the protein folding process is dominated by backbone hydrogen bonding, an approach based on backbone hydrogen-bonded residue pairings would improve the predicting capabilities of these class algorithms. The reliability of the prediction algorithms depends on the quality of the statistics, therefore, of the data set. In this study, it was aimed to determine the propensities of the backbone hydrogen-bonded residue pairings for secondary structural elements of α-helix and β-sheet in globular proteins using a new and comprehensive data set created from the peptides deposited in Worldwide Protein Data Bank. A master data set including 4882 globular peptide chains with resolution better than 2.5 Å, sequence identity smaller than 25% and length of no shorter than 100 residues were created. Separate data sub sets also were created for helix and sheet structures from master set and each sub set includes 4594 and 4483 chains, respectively. Backbone hydrogen-bonded residue pairings in helices and sheets were detected and the propensities of them were represented as odds ratios (observed/[random or expected]) in matrices. Propensities assigned by this study to the residue pairings in secondary structural elements (as helix, overall strands, parallel strands and antiparallel strands) differ from the previous studies by 19 to 34%. These dissimilarities are important and they would cause further improvements in secondary structure prediction algorithms.

SDRL: a sequence-dependent protein side-chain rotamer library

Since the introduction of the first protein side-chain rotamer library (RL) almost half a century ago, RLs have been components of many programs and algorithms in structural bioinformatics. Based on the dependence of side-chain dihedral angles on the local backbone, three types of RLs have been identified: backbone-independent, secondary-structure-dependent and backbone-dependent. In all previous studies, the effect of sequence specificity on side-chain conformational preferences was neglected. In the effort to develop a new class of RLs, we considered that the side-chain conformation of the central residue in each triplet on a protein backbone depends on the sequence of the triplet; therefore, we developed a sequence-dependent rotamer library (SDRL). To accomplish this, 400 possible triplet sequences for 18 natural amino acids as the central residue, which corresponds to 7200 triplet sequences in total, were considered. Searching the set of 11 546 selected PDB entries for the 7200 triplet sequences resulted in 2 364 541 instances occurring for 18 amino acids. Our results show that Leu and Val experience minimal impact from the adjacent residues in adopting side-chain conformations. Cys, Ile, Trp, His, Asp, Met, Glu, Gln, Arg and Lys, on the other hand, adopt their side-chain conformations mostly based on the adjacent residues on the backbone. The remaining residue types were moderately dependent on the adjacent residues. Using the new library, side-chain repacking algorithms can find preferred conformations of each residue more easily than with other backbone-independent RLs.

Mapping side chain interactions at protein helix termini

Background

Interactions that involve one or more amino acid side chains near the ends of protein helices stabilize helix termini and shape the geometry of the adjacent loops, making a substantial contribution to overall protein structure. Previous work has identified key helix-terminal motifs, such as Asx/ST N-caps, the capping box, and hydrophobic and electrostatic interactions, but important questions remain, including: 1) What loop backbone geometries are favoured by each motif? 2) To what extent are multi-amino acid motifs likely to represent genuine cooperative interactions? 3) Can new motifs be identified in a large, recent dataset using the latest bioinformatics tools?

Results

Three analytical tools are applied here to answer these questions. First, helix-terminal structures are partitioned by loop backbone geometry using a new 3D clustering algorithm. Next, Cascade Detection, a motif detection algorithm recently published by the author, is applied to each cluster to determine which sequence motifs are overrepresented in each geometry. Finally, the results for each motif are presented in a CapMap, a 3D conformational heatmap that displays the distribution of the motif’s overrepresentation across loop geometries, enabling the rapid isolation and characterization of the associated side chain interaction. This work identifies a library of geometry-specific side chain interactions that provides a new, detailed picture of loop structure near the helix terminus. Highlights include determinations of the favoured loop geometries for the Asx/ST N-cap motifs, capping boxes, “big” boxes, and other hydrophobic, electrostatic, H-bond, and pi stacking interactions, many of which have not been described before.

Conclusions

This work demonstrates that the combination of structural clustering and motif detection in the sequence space can efficiently identify side chain motifs and map them to the loop geometries which they support. Protein designers should find this study useful, because it identifies side chain interactions which are good candidates for inclusion in synthetic helix-terminal loops with specific desired geometries, since they are used in nature to support these geometries. The techniques described here can also be applied to map side chain interactions associated with other structural components of proteins such as beta and gamma turns.

Electronic supplementary material

The online version of this article (doi:10.1186/s12859-015-0671-4) contains supplementary material, which is available to authorized users.

Neighbor preferences of amino acids and context-dependent effects of amino acid substitutions in human, mouse, and dog

Amino acids show apparent propensities toward their neighbors. In addition to preferences of amino acids for their neighborhood context, amino acid substitutions are also considered to be context-dependent. However, context-dependence patterns of amino acid substitutions still remain poorly understood. Using relative entropy, we investigated the neighbor preferences of 20 amino acids and the context-dependent effects of amino acid substitutions with protein sequences in human, mouse, and dog. For 20 amino acids, the highest relative entropy was mostly observed at the nearest adjacent site of either N- or C-terminus except C and G. C showed the highest relative entropy at the third flanking site and periodic pattern was detected at G flanking sites. Furthermore, neighbor preference patterns of amino acids varied greatly in different secondary structures. We then comprehensively investigated the context-dependent effects of amino acid substitutions. Our results showed that nearly half of 380 substitution types were evidently context dependent, and the context-dependent patterns relied on protein secondary structures. Among 20 amino acids, P elicited the greatest effect on amino acid substitutions. The underlying mechanisms of context-dependent effects of amino acid substitutions were possibly mutation bias at a DNA level and natural selection. Our findings may improve secondary structure prediction algorithms and protein design; moreover, this study provided useful information to develop empirical models of protein evolution that consider dependence between residues.

Analyses of the general rule on residue pair frequencies in local amino acid sequences of soluble, ordered proteins

The amino acid sequences of soluble, ordered proteins with stable structures have evolved due to biological and physical requirements, thus distinguishing them from random sequences. Previous analyses have focused on extracting the features that frequently appear in protein substructures, such as α-helix and β-sheet, but the universal features of protein sequences have not been addressed. To clarify the differences between native protein sequences and random sequences, we analyzed 7368 soluble, ordered protein sequences, by inspecting the observed and expected occurrences of 400 amino acid pairs in local proximity, up to 10 residues along the sequence in comparison with their expected occurrence in random sequence. We found the trend that the hydrophobic residue pairs and the polar residue pairs are significantly decreased, whereas the pairs between a hydrophobic residue and a polar residue are increased. This trend was universally observed regardless of the secondary structure content but was not observed in protein sequences that include intrinsically disordered regions, indicating that it can be a general rule of protein foldability. The possible benefits of this rule are discussed from the viewpoints of protein aggregation and disorder, which are both caused by low-complexity regions of hydrophobic or polar residues.

The alphabet of intrinsic disorder: I. Act like a Pro: On the abundance and roles of proline residues in intrinsically disordered proteins

A significant fraction of every proteome is occupied by biologically active proteins that do not form unique three-dimensional structures. These intrinsically disordered proteins (IDPs) and IDP regions (IDPRs) have essential biological functions and are characterized by extensive structural plasticity. Such structural and functional behavior is encoded in the amino acid sequences of IDPs/IDPRs, which are enriched in disorder-promoting residues and depleted in order-promoting residues. In fact, amino acid residues can be arranged according to their disorder-promoting tendency to form an alphabet of intrinsic disorder that defines the structural complexity and diversity of IDPs/IDPRs. This review is the first in a series of publications dedicated to the roles that different amino acid residues play in defining the phenomenon of protein intrinsic disorder. We start with proline because data suggests that of the 20 common amino acid residues, this one is the most disorder-promoting.

Exploring ORFan domains in giant viruses: structure of mimivirus sulfhydryl oxidase R596

The mimivirus genome contains many genes that lack homologs in the sequence database and are thus known as ORFans. In addition, mimivirus genes that encode proteins belonging to known fold families are in some cases fused to domain-sized segments that cannot be classified. One such ORFan region is present in the mimivirus enzyme R596, a member of the Erv family of sulfhydryl oxidases. We determined the structure of a variant of full-length R596 and observed that the carboxy-terminal region of R596 assumes a folded, compact domain, demonstrating that these ORFan segments can be stable structural units. Moreover, the R596 ORFan domain fold is novel, hinting at the potential wealth of protein structural innovation yet to be discovered in large double-stranded DNA viruses. In the context of the R596 dimer, the ORFan domain contributes to formation of a broad cleft enriched with exposed aromatic groups and basic side chains, which may function in binding target proteins or localization of the enzyme within the virus factory or virions. Finally, we find evidence for an intermolecular dithiol/disulfide relay within the mimivirus R596 dimer, the first such extended, intersubunit redox-active site identified in a viral sulfhydryl oxidase.

Solution structure of the leader sequence of the patellamide precursor peptide, PatE1-34

The solution structure of the leader sequence of the patellamide precursor peptide was analysed by using CD and determined with NOE-restrained molecular dynamics calculations. This leader sequence is highly conserved in the precursor peptides of some other cyanobactins harbouring heterocycles, and is assumed to play a role in targeting the precursor peptide to the post-translational machinery. The sequence was observed to form an alpha-helix spanning residues 13-28 with a hydrophobic surface on one side of the helix. This hydrophobic surface is proposed to be the site of the initial binding with modifying enzymes.

Amino acid pair- and triplet-wise groupings in the interior of α-helical segments in proteins

A statistical approach has been applied to analyse primary structure patterns at inner positions of α-helices in proteins. A systematic survey was carried out in a recent sample of non-redundant proteins selected from the Protein Data Bank, which were used to analyse α-helix structures for amino acid pairing patterns. Only residues more than three positions apart from both termini of the α-helix were considered as inner. Amino acid pairings i, i+k (k=1, 2, 3, 4, 5), were analysed and the corresponding 20×20 matrices of relative global propensities were constructed. An analysis of (i, i+4, i+8) and (i, i+3, i+4) triplet patterns was also performed. These analysis yielded information on a series of amino acid patterns (pairings and triplets) showing either high or low preference for α-helical motifs and suggested a novel approach to protein alphabet reduction. In addition, it has been shown that the individual amino acid propensities are not enough to define the statistical distribution of these patterns. Global pair propensities also depend on the type of pattern, its composition and orientation in the protein sequence. The data presented should prove useful to obtain and refine useful predictive rules which can further the development and fine-tuning of protein structure prediction algorithms and tools.

Position-specific propensities of amino acids in the β-strand

Background

Despite the importance of β-strands as main building blocks in proteins, the propensity of amino acid in β-strands is not well-understood as it has been more difficult to determine experimentally compared to α-helices. Recent studies have shown that most of the amino acids have significantly high or low propensity towards both ends of β-strands. However, a comprehensive analysis of the sequence dependent amino acid propensities at positions between the ends of the β-strand has not been investigated.

Results

The propensities of the amino acids calculated from a large non-redundant database of proteins are found to be highly position-specific and vary continuously throughout the length of the β-strand. They follow an unexpected characteristic periodic pattern in inner positions with respect to the cap residues in both termini of β-strands; this periodic nature is markedly different from that of the α-helices with respect to the strength and pattern in periodicity. This periodicity is not only different for different amino acids but it also varies considerably for the amino acids belonging to the same physico-chemical group. Average hydrophobicity is also found to be periodic with respect to the positions from both termini of β-strands.

Conclusions

The results contradict the earlier perception of isotropic nature of amino acid propensities in the middle region of β-strands. These position-specific propensities should be of immense help in understanding the factors responsible for β-strand design and efficient prediction of β-strand structure in unknown proteins.

Residue patterning in helix interiorsThis paper is one of a selection of papers published in this special issue entitled “Canadian Society of Biochemistry, Molecular & Cellular Biology 52nd Annual Meeting — Protein Folding: Principles and Diseases” and has undergone the Journal's usual peer review process

The α-helix remains a focus of research because of its importance to protein folding and structure. Nevertheless, despite numerous empirical, computational, and theoretical studies, the fundamental structural properties governing their formation and stability are still unclear. We have examined the statistical occurrence of polar and apolar residue patterning in helical interiors in a large, non-redundant dataset, and compared these patterns with those found in other structural environments. While the familiar amphipathic distributions for both polar and apolar residues are evident, our analysis also finds significant differences between these distributions. Non-amphipathic signals can also be discerned within both distributions. Most interestingly, among various positional patterning, an analysis of immediate (i,i + 1) helical neighbours found: (i) clear neighbouring preferences, with high (low) occurrences of hydrophobics (hydrophilics) next to Gly, Pro, and short polar residues; (ii) high negative (positive) correlation between residue helical propensities and the degree of neighbouring hydrophobicity (hydrophilicity); and (iii) a preferred ordering among neighbours, implying an inherent helix directionality. Because (i,i + 1) helical pairs have limited side chain – side chain interactions, thermodynamic considerations cannot readily explain these observations, nor can evolutionary pressures that enhance tertiary interactions via amphipathicity, as this particular spacing does not segregate residues onto either the same or opposing helical faces. We suggest that the mechanism of helix formation may be partly responsible for these observations. In particular, the high negative correlation between residue helical propensities and neighbouring hydrophobicity suggests that hydrophobicity may play a more important role in helix formation than currently recognized.

Folding by numbers: primary sequence statistics and their use in studying protein folding

The exponential growth over the past several decades in the quantity of both primary sequence data available and the number of protein structures determined has provided a wealth of information describing the relationship between protein primary sequence and tertiary structure. This growing repository of data has served as a prime source for statistical analysis, where underlying relationships between patterns of amino acids and protein structure can be uncovered. Here, we survey the main statistical approaches that have been used for identifying patterns within protein sequences, and discuss sequence pattern research as it relates to both secondary and tertiary protein structure. Limitations to statistical analyses are discussed, and a context for their role within the field of protein folding is given. We conclude by describing a novel statistical study of residue patterning in β-strands, which finds that hydrophobic (i,i+2) pairing in β-strands occurs more often than expected at locations near strand termini. Interpretations involving β-sheet nucleation and growth are discussed.

Spatial structure and pH-dependent conformational diversity of dimeric transmembrane domain of the receptor tyrosine kinase EphA1

Eph receptors are found in a wide variety of cells in developing and mature tissues and represent the largest family of receptor tyrosine kinases, regulating cell shape, movements, and attachment. The receptor tyrosine kinases conduct biochemical signals across plasma membrane via lateral dimerization in which their transmembrane domains play an important role. Structural-dynamic properties of the homodimeric transmembrane domain of the EphA1 receptor were investigated with the aid of solution NMR in lipid bicelles and molecular dynamics in explicit lipid bilayer. EphA1 transmembrane segments associate in a right-handed parallel alpha-helical bundle, region (544-569)(2), through the N-terminal glycine zipper motif A(550)X(3)G(554)X(3)G(558). Under acidic conditions, the N terminus of the transmembrane helix is stabilized by an N-capping box formed by the uncharged carboxyl group of Glu(547), whereas its deprotonation results in a rearrangement of hydrogen bonds, fractional unfolding of the helix, and a realignment of the helix-helix packing with appearance of additional minor dimer conformation utilizing seemingly the C-terminal GG4-like dimerization motif A(560)X(3)G(564). This can be interpreted as the ability of the EphA1 receptor to adjust its response to ligand binding according to extracellular pH. The dependence of the pK(a) value of Glu(547) and the dimer conformational equilibrium on the lipid head charge suggests that both local environment and membrane surface potential can modulate dimerization and activation of the receptor. This makes the EphA1 receptor unique among the Eph family, implying its possible physiological role as an "extracellular pH sensor," and can have relevant physiological implications.

Acetylation of L12 increases interactions in the Escherichia coli ribosomal stalk complex

The ribosomal stalk complex in Escherichia coli consists of L10 and four copies of L7/L12, and is largely responsible for binding and recruiting translation factors. Structural characterisation of this stalk complex is difficult, primarily due to its dynamics. Here, we apply mass spectrometry to follow post-translational modifications and their effect on structural changes of the stalk proteins on intact ribosomes. Our results show that increased acetylation of L12 occurs during the stationary phase on ribosomes harvested from cells grown under optimal conditions. For cells grown in minimal medium, L12 acetylation and processing is altered, resulting in deficient removal of N-terminal methionine in approximately 50% of the L12 population, while processed L12 is almost 100% acetylated. Our results show also that N-acetylation of L12 correlates with an increased stability of the stalk complex in the gas phase. To investigate further the basis of this increased stability, we applied a solution phase hydrogen deuterium exchange protocol to compare the rate of deuterium incorporation in the proteins L9, L10, L11 and L12 as well as the acetylated form of L12 (L7), in situ on the ribosome. Results show that deuterium incorporation is consistently slower for L7 relative to L12 and for L10 when L7 is predominant. Our results imply a tightening of the interaction between L7 and L10 relative to that between L12 and L10. Since acetylation is predominant when cells are grown in minimal medium, we propose that these modifications form part of the cell's strategy to increase stability of the stalk complex under conditions of stress. More generally, our results demonstrate that it is possible to discern the influence of a 42 Da post-translational modification by mass spectrometry and to record subtle changes in hydrogen/deuterium exchange within the context of an intact 2.5 MDa particle.