Discordant and chameleon sequences: their distribution and implications for amyloidogenicity

Identification of ambiguous encoding in protein secondary structure is paramount to develop an understanding of key protein segments underlying amyloid diseases. We investigate two types of structurally ambivalent peptides, which were hypothesized in the literature as indicators of amyloidogenic proteins: discordant α-helices and chameleon sequences. Chameleon sequences are peptides discovered experimentally in different secondary-structure types. Discordant α-helices are α-helical stretches with strong β-strand propensity or prediction. To assess the distribution of these features in known protein structures, and their potential role in amyloidogenesis, we analyzed the occurrence of discordant α-helices and chameleon sequences in nonredundant sets of protein domains (n = 4263) and amyloidogenic proteins extracted from the literature (n = 77). Discordant α-helices were identified if discordance was observed between known secondary structures and secondary-structure predictions from the GOR-IV and PSIPRED algorithms. Chameleon sequences were extracted by searching for identical sequence words in α-helices and β-strands. We defined frustrated chameleons and very frustrated chameleons based on varying degrees of total β propensity ≥α propensity. To our knowledge, this is the first study to discern statistical relationships between discordance, chameleons, and amyloidogenicity. We observed varying enrichment levels for some categories of discordant and chameleon sequences in amyloidogenic sequences. Chameleon sequences are also significantly enriched in proteins that have discordant helices, indicating a clear link between both phenomena. We identified the first set of discordant-chameleonic protein segments we predict may be involved in amyloidosis. We present a detailed analysis of discordant and chameleons segments in the family of one of the amyloidogenic proteins, the Prion Protein.     

Foldability and chameleon propensity of fold‐switching protein sequences

It has been shown recently by Porter & Looger that a significant number of proteins exist that can form more than one stable fold. This note examines the sequences of these fold‐switching proteins by (a) calculating their foldability scores recently introduced by the present author and (b) comparing the propensity of chameleon sequences in fold switchers and in non fold switchers. It has been found that the average foldability score of the fold switchers indicates weaker foldability. As for the propensity of chameleon sequences of length 5 to 7 it was found, somewhat surprisingly, that there is only a very small difference between the fold switchers and the non fold switchers. Furthermore, when looking at amino acid propensities, for several amino acids there was even an opposing trend in the deviation of their propensities from the overall amino acid propensities.     

Effects of chameleon dispense-to-plunge speed on particle concentration,complex formation,and final resolution: A case study using the Neisseria gonorrhoeae ribonucleotide reductase inactive complex

Ribonucleotide reductase (RNR) is an essential enzyme that converts ribonucleotides to deoxyribonucleotides and is a promising antibiotic target, but few RNRs have been structurally characterized. We present the use of the chameleon, a commercially-available piezoelectric cryogenic electron microscopy plunger, to address complex denaturation in the Neisseria gonorrhoeae class Ia RNR. Here, we characterize the extent of denaturation of the ring-shaped complex following grid preparation using a traditional plunger and using a chameleon with varying dispense-to-plunge times. We also characterize how dispense-to-plunge time influences the amount of protein sample required for grid preparation and preferred orientation of the sample. We demonstrate that the fastest dispense-to-plunge time of 54 ms is sufficient for generation of a data set that produces a high quality structure, and that a traditional plunging technique or slow chameleon dispense-to-plunge times generate data sets limited in resolution by complex denaturation. The 4.3 Å resolution structure of Neisseria gonorrhoeae class Ia RNR in the inactive α4β4 oligomeric state solved using the chameleon with a fast dispense-to-plunge time yields molecular information regarding similarities and differences to the well studied Escherichia coli class Ia RNR α4β4 ring.     

Computer programs for the analysis of nucleotide sequences (MALK)

A system for the computer analysis of nucleic acid and protein sequences ("Helix") is described. Format of the DNA sequences is EMBL--compatible and may be easily commented with the help of convenient menus. "Helix" has also following possibilities: an effective alignment of gele reading data and formation of the final sequence; simple making of recombined molecules "in calcular"; calculations of nucleotide and dinucleotide distribution along the sequence; looking for coding frames; calculations percentage of codons and amino acids in coding frames; searching for direct and inverted repeats; sequences alignment; protein secondary structure prediction; restriction mapping; DNA--protein translation. "Helix" also contain programs for RNA-structure prediction, looking for homologies throughover the EMAL bank, choosing optimal sequence for probes and searching promoters. All the programs are written at FORTRAN-77 and automatically translated into FORTRAN-4. "Helix" require only 64 kbite.     

IMMUNOHISTOCHEMICAL LOCALIZATION OF GABA IN CHAMELEON RETINA (CHAMALEO CHAMALEO)

We used a policlonal antiserum against GABA and demonstated GABA-immunoreactivity (GABA-IR) in several populations of amacrine cells in the inner nuclear layer (INL), and other cells in the inner plexiform layer (IPL) of the central and peripheral retina of the chameleon. Horizontal cells do not contain GABA-IR and the chameleon retina is therefore an exception among non-mammals. GABA-IR was not seen in cell bodies in the position of photoreceptor, bipolar and interplexiform cells suggesting that GABA is not involved in synaptic transmission in the outer plexiform layer of chameleon retina.     

Grafting a short chameleon sequence from αB crystallin into a β-sheet scaffold protein

Many protein and peptide sequences are self-assembled into β-sheet-rich fibrous structures called amyloids. Their atomic details provide insights into fundamental knowledge related to amyloid diseases. To study the detailed structure of the amyloid, we have developed a model system that mimics the self-assembling process of the amyloid within a water-soluble protein, termed peptide self-assembly mimic (PSAM). PSAM enables capturing of a peptide sequence within a water-soluble protein, thus making structural and energetics-related studies possible. In this work, we extend our PSAM approach to a naturally occurring chameleon sequence from αB crystallin. We chose “Val–Leu–Gly–Asp–Val (VLGDV)”, a five amino-acid sequence, which forms a β-turn in the native structure and a β-barrel in the amyloid oligomer cylindrin, as a grafting sequence to the PSAM scaffold. The crystal structure revealed that the sequence grafting induced β-sheet bending at the grafted site. We further investigated the role of the central glycine residue and found that its role in the β-sheet bending is dependent on the neighboring residues. The ability of PSAM to observe the structural alterations induced by the grafted sequence provides an opportunity to evaluate the structural impact of a sequence from the peptide self-assembly.     

Analysis of protein chameleon sequence characteristics

Conversion of local structural state of a protein from an α-helix to a β-strand is usually associated with a major change in the tertiary structure. Similar changes were observed during the self assembly of amyloidogenic proteins to form fibrils, which are implicated in severe diseases conditions, e.g., Alzheimer disease. Studies have emphasized that certain protein sequence fragments known as chameleon sequences do not have a strong preference for either helical or the extended conformations. Surprisingly, the information on the local sequence neighborhood can be used to predict their secondary at a high accuracy level. Here we report a large scale-analysis of chameleon sequences to estimate their propensities to be associated with different local structural states such as α -helices, β-strands and coils. With the help of the propensity information derived from the amino acid composition, we underline their complexity, as more than one quarter of them prefers coil state over to the regular secondary structures. About half of them show preference for both α-helix and β-sheet conformations and either of these two states is favored by the rest.     

ChSeq: A database of chameleon sequences

Chameleon sequences (ChSeqs) refer to sequence strings of identical amino acids that can adopt different conformations in protein structures. Researchers have detected and studied ChSeqs to understand the interplay between local and global interactions in protein structure formation. The different secondary structures adopted by one ChSeq challenge sequence‐based secondary structure predictors. With increasing numbers of available Protein Data Bank structures, we here identify a large set of ChSeqs ranging from 6 to 10 residues in length. The homologous ChSeqs discovered highlight the structural plasticity involved in biological function. When compared with previous studies, the set of unrelated ChSeqs found represents an about 20‐fold increase in the number of detected sequences, as well as an increase in the longest ChSeq length from 8 to 10 residues. We applied secondary structure predictors on our ChSeqs and found that methods based on a sequence profile outperformed methods based on a single sequence. For the unrelated ChSeqs, the evolutionary information provided by the sequence profile typically allows successful prediction of the prevailing secondary structure adopted in each protein family. Our dataset will facilitate future studies of ChSeqs, as well as interpretations of the interplay between local and nonlocal interactions. A user‐friendly web interface for this ChSeq database is available at prodata.swmed.edu/chseq .     

Certain heptapeptide and large sequences representing an entire helix, strand or coil conformation in proteins are associated as chameleon sequences

Helices, strands and coils in proteins of known three-dimensional structure, corresponding to heptapeptide and large sequences (‘probe’ peptides), were scanned against peptide sequences of variable length, comprising seven or more residues that correspond to a different conformation (‘target’ peptides) in protein crystal structures available from the Protein Data Bank (PDB). Where the ‘probe’ and ‘target’ peptide sequences exactly match, they correspond to ‘chameleon’ sequences in protein structures. We observed ∼548 heptapeptide and large chameleon sequences that included peptides in the coil conformation from 53,794 PDB files that were analyzed. However, after excluding several chameleon peptides based on the quality of protein structure data, redundancy and peptides associated with cloning artifacts, such as, histidine-tags, we observed only ten chameleon peptides in structurally different proteins and the maximum length comprised seven amino acid residues. Our analysis suggests that the quality of protein structure data is important for identifying possibly, the ‘true chameleons’ in PDB. Majority of the chameleon sequences correspond to an entire strand in one protein that is observed as part of helix sequence in another protein. The heptapeptide chameleons are characterized with a high propensity of alanine, leucine and valine amino acid residues. The total hydropathy values range between −11.2 and 22.9, the difference in solvent accessibility between 2.0 Å² and 373 Å² units and the difference in total number of residue neighbor contacts between 0 and 7 residues. Our work identifies for the first time heptapeptide and large sequences that correspond to a single complete helix, strand or coil, which adopt entirely different secondary structures in another protein.     

Nonrandom co-occurrence patterns of rainforest chameleons

Null‐model analysis of co‐occurrence patterns is a powerful tool to identify ‘structure’ in community ecology data sets. We evaluated the community structure of chameleons in rainforest regions of Nigeria and Cameroon using available data in the literature, including peer‐reviewed articles and unpublished environmental reports to industries. We performed Monte Carlo simulations (5000 iterations, using the sequential swap algorithm) under several model assumptions to derive co‐occurrence patterns among species. Food and spatial (habitat) segregation patterns in both lowland rainforest and montane forest were investigated. We subjected four indices of co‐occurrence patterns (C‐ratio, number of checkerboard species pairs, number of species combinations, and V‐score) to randomization procedures. Overall, the chameleon communities do not show random organization, but instead exhibit precise deterministic patterns. In lowland rainforest, chameleon communities are assembled deterministically along the food niche resource axis, but not along the habitat niche resource axis. The opposite holds for chameleon communities in montane rainforest. We predict that these patterns can be generalized to other regions of tropical Africa, thus helping to determine the general structure of chameleon communities in tropical African forests.     

Free-energy landscape of a chameleon sequence in explicit water and its inherent alpha/beta bifacial property

A sequence in yeast MATalpha2/MCM1/DNA complex that folds into alpha-helix or beta-hairpin depending on the surroundings has been known as "chameleon" sequence. We obtained the free-energy landscape of this sequence by using a generalized-ensemble method, multicanonical molecular dynamics simulation, to sample the conformational space. The system was expressed with an all-atom model in explicit water, and the initial conformation for the simulation was a random one. The free-energy landscape demonstrated that this sequence inherently has an ability to form either alpha or beta structure: The conformational distribution in the landscape consisted of two alpha-helical clusters with different packing patterns of hydrophobic residues, and four beta-hairpin clusters with different strand-strand interaction patterns. Narrow pathways connecting the clusters were found, and analysis on the pathways showed that a compact structure formed at the N-terminal root of the chameleon sequence controls the cluster-cluster transitions. The free-energy landscape indicates that a small conditional change induces alpha-beta transitions. Additional unfolding simulations done with replacing amino acids showed that the chameleon sequence has an advantage to form an alpha-helix. Current study may be useful to understand the mechanism of diseases resulting from abnormal chain folding, such as amyloid disease.     

Patterns of morphological variation and correlates of habitat use in Chameleons

The lizard family Chameleonidae is one of the most distinctive taxa of all vertebrates. Nonetheless, despite great intrafamilial diversity, little research has been conducted on morphological variation among chameleons. As a first step in this direction, we took morphological measurements on the head, limbs, spines, and tail of 56 species. Our goals were to investigate whether morphological variation correlates with differences in ecology and to examine whether correlations exist among different aspects of morphology. Based on existing information, species were classified either as arboreal or terrestrial, the latter referring to species that are known to use the ground on a regular basis. This study confirms that considerable intrafamilial variation exists among chameleons and that these traits appear to be evolutionarily quite labile. Once the effects of size are removed, functionally related traits tend to covary; however, few correlations are observed between non-functionally related traits. Many differences in the lengths of the limbs and head elements were detected between terrestrial and arboreal species, but the functional and selective significance of these differences is not clear. Further research on chameleon behaviour and ecology is required to understand the factors contributing to chameleon morphological diversity. © 2002 The Linnean Society of London, Biological Journal of the Linnean Society , 2002, 76 , 91–103.     

Conformational contagion in a protein: structural properties of a chameleon sequence

Certain sequences, known as chameleon sequences, take both alpha- and beta-conformations in natural proteins. We demonstrate that a wild chameleon sequence fused to the C-terminal alpha-helix or beta-sheet in foreign stable proteins from hyperthermophiles forms the same conformation as the host secondary structure. However, no secondary structural formation is observed when the sequence is attached to the outside of the secondary structure. These results indicate that this sequence inherently possesses an ability to make either alpha- or beta-conformation, depending on the sequentially neighboring secondary structure if little other nonlocal interaction occurs. Thus, chameleon sequences take on a satellite state through contagion by the power of a secondary structure. We propose this "conformational contagion" as a new nonlocal determinant factor in protein structure and misfolding related to protein conformational diseases.     

The solution structure of a chimeric LEKTI domain reveals a chameleon sequence

The conversion of an alpha-helical to a beta-strand conformation and the presence of chameleon sequences are fascinating from the perspective that such structural features are implicated in the induction of amyloid-related fatal diseases. In this study, we have determined the solution structure of a chimeric domain (Dom1PI) from the multidomain Kazal-type serine proteinase inhibitor LEKTI using multidimensional NMR spectroscopy. This chimeric protein was constructed to investigate the reasons for differences in the folds of the homologous LEKTI domains 1 and 6 [Lauber, T., et al. (2003) J. Mol. Biol. 328, 205-219]. In Dom1PI, two adjacent phenylalanine residues (F28 and F29) of domain 1 were substituted with proline and isoleucine, respectively, as found in the corresponding P4' and P5' positions of domain 6. The three-dimensional structure of Dom1PI is significantly different from the structure of domain 1 and closely resembles the structure of domain 6, despite the sequence being identical to that of domain 1 except for the two substituted phenylalanine residues and being only 31% identical to the sequence of domain 6. The mutation converted a short 3(10)-helix into an extended loop conformation and parts of the long COOH-terminal alpha-helix of domain 1 into a beta-hairpin structure. The latter conformational change occurs in a sequence stretch distinct from the region containing the substituted residues. Therefore, this switch from an alpha-helical structure to a beta-hairpin structure indicates a chameleon sequence of seven residues. We conclude that the secondary structure of Dom1PI is determined not only by the local protein sequence but also by nonlocal interactions.     

Characterization of the near native conformational states of the SAM domain of Ste11 protein by NMR spectroscopy

The sterile alpha motif or SAM domain is one of the most frequently present protein interaction modules with diverse functional attributions. SAM domain of the Ste11 protein of budding yeast plays important roles in mitogen‐activated protein kinase cascades. In the current study, urea‐induced, at subdenaturing concentrations, structural, and dynamical changes in the Ste11 SAM domain have been investigated by nuclear magnetic resonance spectroscopy. Our study revealed that a number of residues from Helix 1 and Helix 5 of the Ste11 SAM domain display plausible alternate conformational states and largest chemical shift perturbations at low urea concentrations. Amide proton (H/D) exchange experiments indicated that Helix 1, loop, and Helix 5 become more susceptible to solvent exchange with increased concentrations of urea. Notably, Helix 1 and Helix 5 are directly involved in binding interactions of the Ste11 SAM domain. Our data further demonstrate that the existence of alternate conformational states around the regions involved in dimeric interactions in native or near native conditions. Proteins 2014; 82:2957–2969. © 2014 Wiley Periodicals, Inc.     

Evolutionary origins of viviparity in Chamaeleonidae

Historically, an understanding of viviparity and its evolution in Old World chameleons (Chamaeleonidae) has lagged behind that of other squamate families. Not only is reproductive information scarce or entirely absent for most chameleon species, but the literature reveals no consensus as to the frequency and ecological circumstances under which chameleon viviparity evolved. We integrated information on reproductive modes for nearly all chameleon species with recently published family-scale phylogenetic and ecological analyses to clarify aspects of reproductive evolution in chameleons. Ancestral-trait reconstructions, after accounting for phylogenetic uncertainty, indicated that viviparity has arisen a minimum of three times in Chamaeleonidae, with each origin of live birth in closed-canopy forests. Our maximum-likelihood optimization therefore did not support the previous hypotheses of one, two or four origins of viviparity in the family. Past claims that arboreality would not allow for evolution of viviparity were also not supported, nor was a recent suggestion that viviparity has reverted to oviparity. However, cold climates of high latitudes and elevations may have selected for viviparity in arboreal chameleons. While peritoneal pigmentation may facilitate viviparity, its role as an exaptation rather than an adaptation remains equivocal without data from a wider range of chameleon species. Based on a comprehensive review of reproductive modes throughout the family, our study has resolved the number of origins of viviparity in Chamaeleonidae and provided evidence that live birth evolved under arboreal conditions on three separate occasions in this enigmatic squamate group. This study also reveals the value of using phylogenetic analysis in a manner that is robust to uncertainty (rather than simple correlational approaches) when the goal is to reconstruct evolutionary sequences and selective pressures.     

Conservation of metal-coordinating residues

As a result of rapid advances in genome sequencing, the pace of discovery of new protein sequences has surpassed that of structure and function determination by orders of magnitude. This is also true for metal-binding proteins, that is, proteins that bind one or more metal atoms necessary for their biological function. While metal binding site geometry and composition have been extensively studied, no large scale investigation of metal-coordinating residue conservation has been pursued so far. In pursuing this analysis, we were able to corroborate anecdotal evidence that certain residues are preferred to others for binding to certain metals. The conservation of most metal-coordinating residues is correlated with residue preference in a statistically significant manner. Additionally, we also established a statistically significant difference in conservation between metal-coordinating and noncoordinating residues. These results could be useful for providing better insight to functional importance of metal-coordinating residues, possibly aiding metal binding site prediction and design, metal-protein complex structure prediction, drug discovery, as well as model fitting to electron-density maps produced by X-ray crystallography.     

Disordered C-terminal domain of tyrosyl-tRNA synthetase: secondary structure prediction.

The C-terminal domain (residues 320-419) of tyrosyl-tRNA synthetase (TyrRS) from Bacillus stearothermophilus is disordered in the crystal structure and involved in the binding of the anticodon arm of tRNA(Tyr). The sequences of 11 TyrRSs of prokaryotic or mitochondrial origins were aligned and the alignment showed the existence of conserved residues in the sequences of the C-terminal domains. A consensus could be deduced from the application of five programs of secondary structure prediction to the 11 sequences of the query set. These results suggested that the sequences of the C-terminal domains determined a precise and conserved secondary structure. They predicted that the C-terminal domain would have a mixed fold (alpha/beta or alpha+beta), with the alpha-helices in the first half of the sequence and the beta-strands mainly in its second half. Several programs of fold recognition from sequence alone, by threading onto known structures, were applied but none of them identified a type of fold that would be common to the different sequences of the query set. Therefore, the fold of the C-terminal, anticodon binding domain might be novel.     

The Cancer Mutation D83V Induces an α-Helix to β-Strand Conformation Switch in MEF2B

MEF2B is a major target of somatic mutations in non-Hodgkin lymphoma. Most of these mutations are non-synonymous substitutions of surface residues in the MADS-box/MEF2 domain. Among them, D83V is the most frequent mutation found in tumor cells. The link between this hotspot mutation and cancer is not well understood. Here we show that the D83V mutation induces a dramatic α-helix to β-strand switch in the MEF2 domain. Located in an α-helix region rich in β-branched residues, the D83V mutation not only removes the extensive helix stabilization interactions but also introduces an additional β-branched residue that further shifts the conformation equilibrium from α-helix to β-strand. Cross-database analyses of cancer mutations and chameleon sequences revealed a number of well-known cancer targets harboring β-strand favoring mutations in chameleon α-helices, suggesting a commonality of such conformational switch in certain cancers and a new factor to consider when stratifying the rapidly expanding cancer mutation data.     

Mind the gaps: evidence of bias in estimates of multiple sequence alignments

Multiple sequence alignment (MSA) is a crucial first step in the analysis of genomic and proteomic data. Commonly occurring sequence features, such as deletions and insertions, are known to affect the accuracy of MSA programs, but the extent to which alignment accuracy is affected by the positions of insertions and deletions has not been examined independently of other sources of sequence variation. We assessed the performance of 6 popular MSA programs (ClustalW, DIALIGN-T, MAFFT, MUSCLE, PROBCONS, and T-COFFEE) and one experimental program, PRANK, on amino acid sequences that differed only by short regions of deleted residues. The analysis showed that the absence of residues often led to an incorrect placement of gaps in the alignments, even though the sequences were otherwise identical. In data sets containing sequences with partially overlapping deletions, most MSA programs preferentially aligned the gaps vertically at the expense of incorrectly aligning residues in the flanking regions. Of the programs assessed, only DIALIGN-T was able to place overlapping gaps correctly relative to one another, but this was usually context dependent and was observed only in some of the data sets. In data sets containing sequences with non-overlapping deletions, both DIALIGN-T and MAFFT (G-INS-I) were able to align gaps with near-perfect accuracy, but only MAFFT produced the correct alignment consistently. The same was true for data sets that comprised isoforms of alternatively spliced gene products: both DIALIGN-T and MAFFT produced highly accurate alignments, with MAFFT being the more consistent of the 2 programs. Other programs, notably T-COFFEE and ClustalW, were less accurate. For all data sets, alignments produced by different MSA programs differed markedly, indicating that reliance on a single MSA program may give misleading results. It is therefore advisable to use more than one MSA program when dealing with sequences that may contain deletions or insertions, particularly for high-throughput and pipeline applications where manual refinement of each alignment is not practicable.