Loop folds in proteins and evolutionary conservation of folding nuclei

We show that loops of close contacts involving hydrophobic residues are important in protein folding. Contrary to Berezovsky Berezovsky and Trifonov (J Biomol Struct Dyn 20, 5-6, 2002) the loops important in protein folding usually are much larger in size than 23-31 residues, being instead comparable to the size of the protein for single domain proteins. Additionally what is important are not single loop contacts, but a highly interconnected network of such loop contacts, which provides extra stability to a protein fold and which leads to their conservation in evolution.     

Cloning of a full-length complementary DNA for an Artemia salina glycine-rich protein. Structural relationship with RNA binding proteins

Overlapping cDNAs have been isolated containing all the coding sequences for Artemia salina protein GRP33, a glycine-rich protein (16.6 mol % glycine), with a molecular weight of 32,992. GRP33 is closely related to HD40, the major protein component of Artemia heterogeneous nuclear ribonucleoprotein particles, and shares certain characteristics with other RNA binding proteins. The C-terminal region (123 amino acids) contains 39 glycine residues. This region has multiple arginine residues flanked by glycines, resembling the glycine-dimethylarginine clusters present in other RNA binding proteins. Secondary structure predictions for the protein reveal two distinct domains: a hydrophilic C-terminal domain with an extended conformation and a larger N-terminal domain with a number of alpha-helices and beta-sheets.     

A computational method for the design of nested proteins by loop‐directed domain insertion

The computational design of novel nested proteins—in which the primary structure of one protein domain (insert) is flanked by the primary structure segments of another (parent)—would enable the generation of multifunctional proteins. Here we present a new algorithm, called Loop‐Directed Domain Insertion (LooDo), implemented within the Rosetta software suite, for the purpose of designing nested protein domain combinations connected by flexible linker regions. Conformational space for the insert domain is sampled using large libraries of linker fragments for linker‐to‐parent domain superimposition followed by insert‐to‐linker superimposition. The relative positioning of the two domains (treated as rigid bodies) is sampled efficiently by a grid‐based, mutual placement compatibility search. The conformations of the loop residues, and the identities of loop as well as interface residues, are simultaneously optimized using a generalized kinematic loop closure algorithm and Rosetta EnzymeDesign, respectively, to minimize interface energy. The algorithm was found to consistently sample near‐native conformations and interface sequences for a benchmark set of structurally similar but functionally divergent domain‐inserted enzymes from the α/β hydrolase superfamily, and discriminates well between native and nonnative conformations and sequences, although loop conformations tended to deviate from the native conformations. Furthermore, in cross‐domain placement tests, native insert‐parent domain combinations were ranked as the best‐scoring structures compared to nonnative domain combinations. This algorithm should be broadly applicable to the design of multi‐domain protein complexes with any combination of inserted or tandem domain connections.     

Characteristics and prediction of domain linker sequences in multi-domain proteins

To facilitate swift structural characterizations, structural genomic/proteomic projects need to divide large multi-domain proteins into structural domains and to determine their structures separately. Thus, the assignment of structural domains based solely on sequence information, especially on the physico-chemical properties of the amino acid sequences, could be very helpful for such projects. In this study, we examined the characteristics of domain linker sequences, which are loop sequences connecting two structural domains. To this end, we prepared a set of 101 non-redundant multi-domain protein sequences with known structures, and performed an analysis of the linker sequences. The analysis revealed that the frequencies of five (Pro, Gly, Asp, Asn, Lys) amino acid residues differed significantly between the linker and non-linker loop sequences. Moreover, we observed a similar deviation for the residue pair frequencies between the two types of loop sequences. Finally, we describe an automated method, based on the above analysis, to detect loops that have high probabilities of being domain linkers in a protein sequence.     

Helix packing in polytopic membrane proteins: role of glycine in transmembrane helix association.

The nature and distribution of amino acids in the helix interfaces of four polytopic membrane proteins (cytochrome c oxidase, bacteriorhodopsin, the photosynthetic reaction center of Rhodobacter sphaeroides, and the potassium channel of Streptomyces lividans) are studied to address the role of glycine in transmembrane helix packing. In contrast to soluble proteins where glycine is a noted helix breaker, the backbone dihedral angles of glycine in transmembrane helices largely fall in the standard alpha-helical region of a Ramachandran plot. An analysis of helix packing reveals that glycine residues in the transmembrane region of these proteins are predominantly oriented toward helix-helix interfaces and have a high occurrence at helix crossing points. Moreover, packing voids are generally not formed at the position of glycine in folded protein structures. This suggests that transmembrane glycine residues mediate helix-helix interactions in polytopic membrane proteins in a fashion similar to that seen in oligomers of membrane proteins with single membrane-spanning helices. The picture that emerges is one where glycine residues serve as molecular notches for orienting multiple helices in a folded protein complex.     

Investigating by CD the molecular mechanism of elasticity of elastomeric proteins

Elastomeric proteins are widespread in the animal kingdom, and their main function is to confer elasticity and resilience to organs and tissues. Besides common functional properties, elastomeric proteins share a common sequence design. They are usually constituted by repetitive sequences with a high content of glycine residues. From a conformational point of view, all the elastomeric proteins since now analyzed show a dynamic equilibria between folded (mainly beta-turns) and extended (polyproline II and beta-strands) conformations that could be at the origin of the high entropy of the relaxed state. As a matter of fact, elastin, lamprin, abductin, as well as the PEVK domain of titin share the same conformational ensemble, thus pointing to a common molecular mechanism as the origin of elasticity. CD spectroscopy represents the proper spectroscopic technique to be used overall because of its particular sensitivity to the presence of PPII structure. Its use in the molecular studies of elastin, abductin, and lamprin as well as the recently analyzed protein resilin will be presented.     

Relationship between structure and amino acid sequence of strongly twisted and coiled β-hairpins in globular proteins

β-Hairpins are widespread in proteins, and it is possible to find them both within β-sheets and separately. In this work, a comparative analysis of amino acid sequences of β-strands within strongly twisted β-hairpins from different structural protein subclasses has been conducted. Strongly twisted and coiled β-hairpin generates in the space a right double helix out of β-strands that are connected by a loop region (connections). The frequencies of amino acid residues on the internal (concave) and external (convex) surfaces of strongly twisted β-hairpins have been determined (220 β-hairpins from nonhomologous proteins were studied). The concave surface of these β-hairpins is mainly generated by hydrophobic residues, while the convex surface by hydrophilic residues; accordingly, the alternation of hydrophobic internal and hydrophilic external residues is observed in their amino acid sequences. Amino acid residues of glycine and alanine (especially in places of the largest twisting of the strands) were anomalously frequently found in internal positions of strongly twisted and coiled β-hairpins. It was established that internal positions never contain the proline residues, while external positions in the twisting region contain them in a relatively large amount. It was demonstrated that at least one amino acid residue in α_L- or ε-conformation is required for generation of relatively short (up to 7 amino acid residues) connection. As a rule, these positions are occupied by glycines. Thus, not only the alternation of hydrophobic and hydrophilic amino acid residues, but also the presence of one or two glycine residues in the connection region and the excess of glycines and alanines in the places of the largest strand twisting on the concave surface, as well as the presence of prolines on the convex surface, are required to generate a strongly twisted and coiled β-hairpin.     

Structures of N-terminally acetylated proteins

Primary structures of 250 characterized proteins with N-terminally acetylated residues were correlated with residue distributions and other data. Excluding multiple forms derived from characterized species variants, the structures represent 105 different types of acetylated proteins. Results of comparisons extend previous suggestions based on fewer structures and define relationships further. The N-terminal residue that is acetylated is of a limited type and is frequently a small residue, with a heavy over-representation of serine and alanine. However, the occurrence of methionine at the acetylated position is also high, whereas that of glycine is less frequent than previously estimated. Lysine is over-represented in the N-terminal region, as is aspartic and glutamic acids at a few positions close to the acetylated N-terminus (especially the adjacent position). Finally, distributions of branched-chain residues in the N-terminal region of acetylated proteins are altered in relation to those of proteins in general, isoleucine is over-represented, and leucine and valine are under-represented. The results suggest that alpha-amino-acetylated proteins have special residues in N-terminally non-hydrophobic structures. Data are compatible with a protective function for acetylation but do not exclude further role(s) in processing or other special functions.     

Identification of a glycine motif required for packing in EmrE, a multidrug transporter from Escherichia coli

Glycine residues may play functional and structural roles in membrane proteins. In this work we studied the role of glycine residues in EmrE, a small multidrug transporter from Escherichia coli. EmrE extrudes various drugs across the plasma membrane in exchange with protons and, as a result, confers resistance against their toxic effects. Each of 12 glycine residues was replaced by site-directed mutagenesis. Four of the 12 glycine residues in EmrE are evolutionary conserved within the small multidrug resistance family of multidrug transporters. Our analysis reveals that only two (Gly-67 and Gly-97) of these four highly conserved residues are essential for transporter activity. Moreover, two glycine positions that are less conserved, Gly-17 and Gly-90, demonstrate also a nil phenotype when substituted. Our present results identifying Gly-17 and Gly-67 as irreplaceable reinforce the importance of previously defined functional clusters. Two essential glycine residues, Gly-90 and Gly-97, form a protein motif in which glycine residues are separated by six other residues (GG7). Upon substitution of glycine in these positions, the protein ability to form dimers is impaired as evaluated by cross-linking and pull-down experiments.     

Molecular basis for substrate discrimination by glycine transporters

Glycine is an inhibitory neurotransmitter in the spinal cord and brain stem, where it acts on strychnine-sensitive glycine receptors, and is also an excitatory neurotransmitter throughout the brain and spinal cord, where it acts on the N-methyl-d-aspartate family of receptors. There are two Na(+)/Cl(-)-dependent glycine transporters, GLYT1 and GLYT2, which control extracellular glycine concentrations and these transporters show differences in substrate selectivity and blocker sensitivity. A bacterial Na(+)-dependent leucine transporter (LeuT(Aa)) has recently been crystallized and its structure determined. When the amino acid residues within the leucine binding site of LeuT(Aa) are aligned with residues of the two glycine transporters there are a number of identical residues and also some key differences. In this report, we demonstrate that the LeuT(Aa) structure represents a good working model of the Na(+)/Cl(-)-dependent neurotransmitters and that differences in substrate selectivity can be attributed to a single difference of a glycine residue in transmembrane domain 6 of GLYT1 for a serine residue at the corresponding position of GLYT2.     

SPXX, a frequent sequence motif in gene regulatory proteins

A new DNA-binding unit, composed of four amino acid residues and common in gene regulatory proteins, is proposed. The occurrences of the sequences Ser-Pro-X-X (SPXX) and Thr-Pro-X-X (TPXX) in gene regulatory proteins are compared with those in general proteins. These sequences are found more frequently in gene regulatory proteins including homoeotic gene products, segmentation gene products, steroid hormone receptors and certain oncogene products, than they are in DNA-binding proteins that are not directly involved in gene regulation, such as the core histones, or in general proteins. It is therefore suggested that these sequences contribute to DNA-binding in a manner important for gene regulation. Amino acid residues characteristic of the types of proteins are found as the variable residues X: basic residues, Lys and Arg, in histones, H1 and sea urchin spermatogenous H2B; Tyr in RNA polymerase II; and Ser, Thr, Ala, Leu and Pro in other gene regulatory proteins S(T)PXX sequences are located on either side of other DNA-recognizing units such as Zn fingers, helix-turn-helices, and cores of histones. The structure of a S(T)PXX sequence is presumed to be a beta-turn I stabilized by two hydrogen bonds, and its potential mode of DNA-binding is discussed.     

Structural basis for diversity of the EF-hand calcium-binding proteins

The calcium binding proteins of the EF-hand super-family are involved in the regulation of all aspects of cell function. These proteins exhibit a great diversity of composition, structure, Ca2+-binding and target interaction properties. Here, our current understanding of the Ca2+-binding mechanism is assessed. The structures of the EF-hand motifs containing 11-14 amino acid residues in the Ca2+-binding loop are analyzed within the framework of the recently proposed two-step Ca2+-binding mechanism. A hypothesis is put forward that in all EF-hand proteins the Ca2+-binding and the resultant conformational responses are governed by the central structure connecting the Ca2+-binding loops in the two-EF-hand domain. This structure, named EFbeta-scaffold, defines the position of the bound Ca2+, and coordinates the function of the N-terminal (variable and flexible) with the C-terminal (invariable and rigid) parts of the Ca2+-binding loop. It is proposed that the nature of the first ligand of the Ca2+-binding loop is an important determinant of the conformational change. Additional factors, including the interhelical contacts, the length, structure and flexibility of the linker connecting the EF-hand motifs, and the overall energy balance provide the fine-tuning of the Ca2+-induced conformational change in the EF-hand proteins.     

A mutation in Escherichia coli tRNA nucleotidyltransferase that affects only AMP incorporation is in a sequence often associated with nucleotide-binding proteins

Escherichia coli strain 5C15 contains a mutation in the cca gene that decreases AMP incorporation by tRNA nucleotidyltransferase while leaving CMP incorporation unaffected. Earlier studies of the purified mutant enzyme suggested that the mutation was localized to the AMP-incorporating site. In order to analyze this mutation in more detail, the cca gene from strain 5C15 was cloned into plasmid pUC8. Analysis of tRNA nucleotidyltransferase activity in extracts of a strain transformed with this plasmid demonstrated an elevated level of CMP incorporation, but low AMP incorporation, as expected from the properties of the original mutant. Sequence analysis of the mutant cca gene revealed only a single G to A point mutation leading to a glycine to aspartic acid substitution at position 70 of the peptide chain. The amino acid change was localized to one of two Gly-X-Gly-X-X-Gly sequences present in the protein. This sequence has been identified previously near the nucleotide-binding domain of various proteins, but it has not been noted in enzymes that incorporate nucleotide residues. However, other sequences often associated with ATP-binding domains are not found in tRNA nucleotidyltransferase. The implications of these findings for our understanding of nucleotide-binding domains are discussed.     

Sequence-structure signals of 3D domain swapping in proteins

Three-dimensional domain swapping occurs when two or more identical proteins exchange identical parts of their structure to generate an oligomeric unit. It affects proteins with diverse sequences and structures, and is expected to play important roles in evolution, functional regulation and even conformational diseases. Here, we search for traces of domain swapping in the protein sequence, by means of algorithms that predict the structure and stability of proteins using database-derived potentials. Regions whose sequences are not optimal with regard to the stability of the native structure, or showing marked intrinsic preferences for non-native conformations in absence of tertiary interactions are detected in most domain-swapping proteins. These regions are often located in areas crucial in the swapping process and are likely to influence it on a kinetic or thermodynamic level. In addition, cation-pi interactions are frequently observed to zip up the edges of the interface between intertwined chains or to involve hinge loop residues, thereby modulating stability. We end by proposing a set of mutations altering the swapping propensities, whose experimental characterization would contribute to refine our in silico derived hypotheses.     

Pattern recognition of sequence similarities in globular proteins by Fourier analysis: A novel approach to molecular evolution

A new algorithm is introduced for analyzing gene-duplication-independent (orthologous) and gene-duplication-dependent amino acid sequence similarities between proteins of different species. It is based on the calculation of an autocorrelation function D(x) as a Fourier series analogous to that used in crystal analysis by x-ray diffraction. The primary structure of the protein is decomposed into "homopolypeptide-defective sequences" containing identical or similar amino acid residues and vacancies corresponding to the missing amino acid residues. The Fourier transforms F(h) simulating the diffraction patterns of defective linear gratings corresponding to the defective homopolypeptide sequences are calculated. The squared F(h) values are then used as coefficients of Fourier series corresponding to the autocorrelation functions D(x). A peak of D(x) corresponds to a vector of length x, which is the distance between two identical amino acid residues. It is pointed out that optical diffraction methods, instead of computer methods, would also be useful. It is shown through a number of examples that this method allows satisfactory pattern recognition of homologies and internal duplications of an initial segment of the polypeptide chain. In the latter case the value of the above method may be seen from the fact that it detects repeated duplications in proteins such as spinach ferredoxin and myoglobin, for which other methods had either failed or given inconclusive results. The above approach appears most promising for studies of molecular evolution and structure-sequence correlations.     

Identification of a novel mammalian annexin. cDNA cloning, sequence analysis, and ubiquitous expression of the annexin XI gene.

Annexins (or lipocortins) are a family of at least 10 structurally related calcium- and phospholipid-binding proteins. Each protein consists of a conserved core domain having four (or eight) repeats of a segment approximately 70 amino acids in length and a nonconserved, usually short, amino-terminal domain. To date, amino acid sequences for eight distinct mammalian annexins have been predicted from cDNAs. This report describes an additional member of this family, bovine annexin XI, identified by cDNA cloning and sequence analysis. The 503-amino acid deduced protein consists of a core domain of four annexin repeats and a long amino-terminal domain rich in glycine, proline, and tyrosine. This novel annexin gene is expressed in a wide variety of tissues and isolated cells in culture.     

(1)H/(15)N heteronuclear NMR spectroscopy shows four dynamic domains for phospholamban reconstituted in dodecylphosphocholine micelles

We report the backbone dynamics of monomeric phospholamban in dodecylphosphocholine micelles using (1)H/(15)N heteronuclear NMR spectroscopy. Phospholamban is a 52-amino acid membrane protein that regulates Ca-ATPase in cardiac muscle. Phospholamban comprises three structural domains: a transmembrane domain from residues 22 to 52, a connecting loop from 17 to 21, and a cytoplasmic domain from 1 to 16 that is organized in an "L"-shaped structure where the transmembrane and the cytoplasmic domain form an angle of approximately 80 degrees (Zamoon et al., 2003; Mascioni et al., 2002). T(1), T(2), and (1)H/(15)N nuclear Overhauser effect values measured for the amide backbone resonances were interpreted using the model-free approach of Lipari and Szabo. The results point to the existence of four dynamic domains, revealing the overall plasticity of the cytoplasmic helix, the flexible loop, and part of the transmembrane domain (residues 22-30). In addition, using Carr-Purcell-Meiboom-Gill-based experiments, we have characterized phospholamban dynamics in the micros-ms timescale. We found that the majority of the residues in the cytoplasmic domain, the flexible loop, and the first ten residues of the transmembrane domain undergo dynamics in the micros-ms range, whereas minimal dynamics were detected for the transmembrane domain. Hydrogen/deuterium exchange factors measured at different temperatures support the existence of slow motion in both the loop and the cytoplasmic helix. We propose that these dynamic properties are critical factors in the biomolecular recognition of phospholamban by Ca-ATPase and other interacting proteins such as protein kinase A and protein phosphatase 1.     

Sequence and structural homology among membrane-associated domains of CFTR and certain transporter proteins

A sequence comparison of the two membrane-associated (MA) domains of the cystic fibrosis transmembrane conductance regulator (CFTR), multidrug resistance transporter (MDR), and -factor pheromone export system (STE6) proteins, each of which are believed to contain a total of 12 transmembrane (TM) segments, reveals significant amino acid homology and length conservation in the loop regions that connect individual TM sequences. Similar structural homology is observed between these proteins, hemolysin B (HLYB) and the major histocompatibility-linked peptide transporter, HAM1, the latter two which contain a single MA domain composed of six TM segments. In addition, there are specific sequences that are conserved within the TM segments of the five different membrane proteins. This observation suggests that the folding topologies of the MA domains of MDR, STE6, and CFTR in the plasma membrane are likely to be very similar. The sequence analysis also reveals that there are three characteristic motifs (a pair of aromatic residues, LTLXXXXXXP and GXXL) that are conserved in MDR, STE6, HLYB, HAM1, but not in CFTR. We propose that although CFTR may be evolutionarily related to these other membrane proteins, it belongs to a separate subclass.     

Structural and dynamic differences between normal and transforming N-ras gene products: a 31P and isotope-edited 1H NMR study

[15N]Glycine was biosynthetically incorporated into normal cellular N-ras p21 and a position 12 transforming mutant, in order to produce p21 proteins containing several site-specific NMR probes at or near activating positions in the guanine nucleotide binding domain. We have previously assigned all five glycine resonances located in loops directly involved in binding of guanosine diphosphate in the wild-type p21 protein [Campbell-Burk, S., Papastavros, M. Z., McCormick, F., & Redfield, A. G. (1989) Proc. Natl. Acad. Sci. U.S.A 86, 817-820]. In this report, the corresponding glycine resonances in the p21 mutant have been assigned, and spectral differences between normal and mutant p21-guanosine diphosphate (p21.GDP) complexes have been investigated. Our combined 1H[15N] and 31P NMR results show that substitution of aspartate for glycine-12 produces perturbations in the phosphoryl binding domain, near the point of the mutation. Although many of the remaining glycines were unaffected, spectral differences were also observed outside the GDP binding domain. Two of the five active-site glycines in wild-type p21.GDP have very slow amide proton exchange rates with water (kappa less than 2.8 x 10(-5) s-1). The active-site glycines are located in solvent-exposed loops, so their apparent solvent inaccessibility may result from strong hydrogen bond formation between glycine amide protons and bound guanine diphosphate and/or other nearby groups in p21.