Sequence conservation in families whose members have little or no sequence similarity: the four-helical cytokines and cytochromes

Lessons on protein structure from interleukin-4: All disulfides are not created equal

Interleukin-4 (IL-4) is a hematopoietic cytokine composed by a four-helix bundle stabilized by an antiparallel beta-sheet and three disulfide bonds: Cys3-Cys127, Cys24-Cys65, and Cys46-Cys99. IL-4 is involved in several immune responses associated to infection, allergy, autoimmunity, and cancer. Besides its physiological relevance, IL-4 is often used as a "model" for protein design and engineering. Hence, to understand the role of each disulfide in the structure and dynamics of IL-4, we carried out several spectroscopic analyses (circular dichroism [CD], fluorescence, nuclear magnetic resonance [NMR]), and molecular dynamics (MD) simulations on wild-type IL-4 and four IL-4 disulfide mutants. All disulfide mutants showed loss of structure, altered interhelical angles, and looser core packings, showing that all disulfides are relevant for maintaining the overall fold and stability of the four-helix bundle motif, even at very low pH. In the absence of the disulfide connecting both protein termini Cys3-Cys127, C3T-IL4 showed a less packed protein core, loss of secondary structure (~9%) and fast motions on the sub-nanosecond time scale (lower S2 order parameters and larger τc correlation time), especially at the two protein termini, loops, beginning of helix A and end of helix D. In the absence of Cys24-Cys65, C24T-IL4 presented shorter alpha-helices (14% loss in helical content), altered interhelical angles, less propensity to form the small anti-parallel beta-sheet and increased dynamics. Simultaneously deprived of two disulfides (Cys3-Cys127 and Cys24-Cys65), IL-4 formed a partially folded "molten globule" with high 8-anilino-1-naphtalenesulphonic acid-binding affinity and considerable loss of secondary structure (~50%decrease), as shown by the far UV-CD, NMR, and MD data.

Functional equivalence of germ plasm organizers

The proteins Oskar (Osk) in Drosophila and Bucky ball (Buc) in zebrafish act as germ plasm organizers. Both proteins recapitulate germ plasm activities but seem to be unique to their animal groups. Here, we discover that Osk and Buc show similar activities during germ cell specification. Drosophila Osk induces additional PGCs in zebrafish. Surprisingly, Osk and Buc do not show homologous protein motifs that would explain their related function. Nonetheless, we detect that both proteins contain stretches of intrinsically disordered regions (IDRs), which seem to be involved in protein aggregation. IDRs are known to rapidly change their sequence during evolution, which might obscure biochemical interaction motifs. Indeed, we show that Buc binds to the known Oskar interactors Vasa protein and nanos mRNA indicating conserved biochemical activities. These data provide a molecular framework for two proteins with unrelated sequence but with equivalent function to assemble a conserved core-complex nucleating germ plasm.

Multicellular organisms use gametes for their propagation. Gametes are formed from germ cells, which are specified during embryogenesis in some animals by the inheritance of RNP granules known as germ plasm. Transplantation of germ plasm induces extra germ cells, whereas germ plasm ablation leads to the loss of gametes and sterility. Therefore, germ plasm is key for germ cell formation and reproduction. However, the molecular mechanisms of germ cell specification by germ plasm in the vertebrate embryo remain an unsolved question. Proteins, which assemble the germ plasm, are known as germ plasm organizers. Here, we show that the two germ plasm organizers Oskar from the fly and Bucky ball from the fish show similar functions by using a cross species approach. Both are intrinsically disordered proteins, which rapidly changed their sequence during evolution. Moreover, both proteins still interact with conserved components of the germ cell specification pathway. These data might provide a first example of two proteins with the same biological role, but distinct sequence.

Using THz Spectroscopy, Evolutionary Network Analysis Methods, and MD Simulation to Map the Evolution of Allosteric Communication Pathways in c-Type Lysozymes

It is now widely accepted that protein function is intimately tied with the navigation of energy landscapes. In this framework, a protein sequence is not described by a distinct structure but rather by an ensemble of conformations. And it is through this ensemble that evolution is able to modify a protein's function by altering its landscape. Hence, the evolution of protein functions involves selective pressures that adjust the sampling of the conformational states. In this work, we focus on elucidating the evolutionary pathway that shaped the function of individual proteins that make-up the mammalian c-type lysozyme subfamily. Using both experimental and computational methods, we map out specific intermolecular interactions that direct the sampling of conformational states and accordingly, also underlie shifts in the landscape that are directly connected with the formation of novel protein functions. By contrasting three representative proteins in the family we identify molecular mechanisms that are associated with the selectivity of enhanced antimicrobial properties and consequently, divergent protein function. Namely, we link the extent of localized fluctuations involving the loop separating helices A and B with shifts in the equilibrium of the ensemble of conformational states that mediate interdomain coupling and concurrently moderate substrate binding affinity. This work reveals unique insights into the molecular level mechanisms that promote the progression of interactions that connect the immune response to infection with the nutritional properties of lactation, while also providing a deeper understanding about how evolving energy landscapes may define present-day protein function.

Structural Insights into the Neutralization Properties of the Fully Human, Anti-interferon Monoclonal Antibody Sifalimumab

We report the three-dimensional structure of human interferon α-2A (IFN-α2A) bound to the Fab fragment of a therapeutic monoclonal antibody (sifalimumab; IgG1/κ). The structure of the corresponding complex was solved at a resolution of 3.0 Å using molecular replacement and constitutes the first reported structure of a human type I IFN bound to a therapeutic antibody. This study revealed the major contribution made by the first complementarity-determining region in each of sifalimumab light and heavy chains. These data also provided the molecular basis for sifalimumab mechanism of action. We propose that its interferon-neutralizing properties are the result of direct competition for IFN-α2A binding to the IFN receptor subunit 1 (IFNAR1) and do not involve inhibiting IFN-α2A binding to the IFN receptor subunit 2 (IFNAR2).

Using THz time-scale infrared spectroscopy to examine the role of collective, thermal fluctuations in the formation of myoglobin allosteric communication pathways and ligand specificity

In this investigation we use THz time-scale spectroscopy to conduct an initial set of studies on myoglobin with the aim of providing further insight into the global, collective thermal fluctuations in the protein that have been hypothesized to play a prominent role in the dynamic formation of transient ligand channels as well as in shaping the molecular level basis for ligand discrimination. Using the two ligands O2 and CO, we have determined that the perturbation from the heme-ligand complex has a strong influence on the characteristics of the myoglobin collective dynamics that are excited upon binding. Further, the differences detected in the collective protein motions in Mb-O2 compared with those in Mb-CO appear to be intimately tied with the pathways of long-range allosteric communication in the protein, which ultimately determine the trajectories selected by the respective ligands on the path to and from the heme-binding cavity.

Pierced Lasso Bundles are a new class of knot-like motifs

A four-helix bundle is a well-characterized motif often used as a target for designed pharmaceutical therapeutics and nutritional supplements. Recently, we discovered a new structural complexity within this motif created by a disulphide bridge in the long-chain helical bundle cytokine leptin. When oxidized, leptin contains a disulphide bridge creating a covalent-loop through which part of the polypeptide chain is threaded (as seen in knotted proteins). We explored whether other proteins contain a similar intriguing knot-like structure as in leptin and discovered 11 structurally homologous proteins in the PDB. We call this new helical family class the Pierced Lasso Bundle (PLB) and the knot-like threaded structural motif a Pierced Lasso (PL). In the current study, we use structure-based simulation to investigate the threading/folding mechanisms for all the PLBs along with three unthreaded homologs as the covalent loop (or lasso) in leptin is important in folding dynamics and activity. We find that the presence of a small covalent loop leads to a mechanism where structural elements slipknot to thread through the covalent loop. Larger loops use a piercing mechanism where the free terminal plugs through the covalent loop. Remarkably, the position of the loop as well as its size influences the native state dynamics, which can impact receptor binding and biological activity. This previously unrecognized complexity of knot-like proteins within the helical bundle family comprises a completely new class within the knot family, and the hidden complexity we unraveled in the PLBs is expected to be found in other protein structures outside the four-helix bundles. The insights gained here provide critical new elements for future investigation of this emerging class of proteins, where function and the energetic landscape can be controlled by hidden topology, and should be take into account in ab initio predictions of newly identified protein targets.

An HIV-1 envelope glycoprotein trimer with an embedded IL-21 domain activates human B cells

Broadly neutralizing antibodies (bNAbs) that target the HIV-1 envelope glycoproteins (Env) can prevent virus acquisition, but several Env properties limit its ability to induce an antibody response that is of sufficient quantity and quality. The immunogenicity of Env can be increased by fusion to co-stimulatory molecules and here we describe novel soluble Env trimers with embedded interleukin-4 (IL-4) or interleukin-21 (IL-21) domains, designed to activate B cells that recognize Env. In particular, the chimeric Env_IL-21 molecule activated B cells efficiently and induced the differentiation of antibody secreting plasmablast-like cells. We studied whether we could increase the activity of the embedded IL-21 by designing a chimeric IL-21/IL-4 (ChimIL-21/4) molecule and by introducing amino acid substitutions in the receptor binding domain of IL-21 that were predicted to enhance its binding. In addition, we incorporated IL-21 into a cleavable Env trimer and found that insertion of IL-21 did not impair Env cleavage, while Env cleavage did not impair IL-21 activity. These studies should guide the further design of chimeric proteins and Env_IL-21 may prove useful in improving antibody responses against HIV-1.

Molecular cloning and genomic characterization of novel leptin-like genes in salmonids provide new insight into the evolution of the Leptin gene family

In the current study we describe the identification of novel leptin B homologous gene/s in the four salmonid species Atlantic salmon (Salmo salar), rainbow trout (Oncorhynchus mykiss), brown trout (Salmo trutta) and Arctic charr (Salvelinus alpinus). Homology modeling of Salmo salar (Ss) LepB1/B2 suggests that the protein satisfies parameters as long-chain four helical cytokine family and that the basic structural pattern of the protein follows that of human leptin (Zhang et al., 1997). Importantly, the docking studies suggested the SsLepB has binding affinity to the AA residues that identify the leptin binding and FNIII domains of the SsLep receptor (Rønnestad et al., 2010). Phylogenetic analyses support that LepB paralogs have most probably originated by 4R whole genome duplication (WGD) before speciation of the salmonid lineages. LepB1 and LepB2 genes are both present in the two closest relatives, the Atlantic salmon and the brown trout, while rainbow trout and charr have only preserved the long LepB1 variant in their genome. We have defined the sites of SsLepB mRNA expression at key life stages in Atlantic salmon and found that SsLepB1 and SsLepB2, although to different extent, were expressed in redundant and mostly complementary fashion in brain and gills throughout the lifecycle, suggesting that this pair of paralogs is likely undergoing early stages of subfunctionalization. Furthermore, we have quantified the expression profiles of SsLepB genes and of other two recently duplicated salmon leptins (SsLepA1, SsLepA2) during early development and show evidence that in fish, as in mammals and amphibians, leptin could play important roles in growth and development. This study provides an essential groundwork to further elucidate structural and functional evolution of this important hormone in salmonids as well as in other teleosts.

Tissue-protective cytokines: structure and evolution

Cytokines are important mediators of host defense and immunity, and were first identified for their role in immunity to infections. It was then found that some of them are pathogenic mediators in inflammatory diseases and much of the emphasis is now on pro-inflammatory cytokines, also in consideration of the fact that TNF inhibitors became effective drugs in chronic inflammatory diseases. The recent studies on the tissue-protective activities of erythropoietin (EPO) led to the term "tissue-protective cytokine." We discuss here how tissue-protective actions might be common to other cytokines, particularly those of the 4-alpha helical structural superfamily.

Structural basis for the dual recognition of helical cytokines IL-34 and CSF-1 by CSF-1R

Lacking any discernible sequence similarity, interleukin-34 (IL-34) and colony stimulating factor 1 (CSF-1) signal through a common receptor CSF-1R on cells of mononuclear phagocyte lineage. Here, the crystal structure of dimeric IL-34 reveals a helical cytokine fold homologous to CSF-1, and we further show that the complex architecture of IL-34 bound to the N-terminal immunoglobulin domains of CSF-1R is similar to the CSF-1/CSF-1R assembly. However, unique conformational adaptations in the receptor domain geometry and intermolecular interface explain the cross-reactivity of CSF-1R for two such distantly related ligands. The docking adaptations of the IL-34 and CSF-1 quaternary complexes, when compared to the stem cell factor assembly, draw a common evolutionary theme for transmembrane signaling. In addition, the structure of IL-34 engaged by a Fab fragment reveals the mechanism of a neutralizing antibody that can help deconvolute IL-34 from CSF-1 biology, with implications for therapeutic intervention in diseases with myeloid pathogenic mechanisms.

Comprehensively surveying structure and function of RING domains from Drosophila melanogaster

Using a complete set of RING domains from Drosophila melanogaster, all the solved RING domains and cocrystal structures of RING-containing ubiquitin-ligases (RING-E3) and ubiquitin-conjugating enzyme (E2) pairs, we analyzed RING domains structures from their primary to quarternary structures. The results showed that: i) putative orthologs of RING domains between Drosophila melanogaster and the human largely occur (118/139, 84.9%); ii) of the 118 orthologous pairs from Drosophila melanogaster and the human, 117 pairs (117/118, 99.2%) were found to retain entirely uniform domain architectures, only Iap2/Diap2 experienced evolutionary expansion of domain architecture; iii) 4 evolutionary structurally conserved regions (SCRs) are responsible for homologous folding of RING domains at the superfamily level; iv) besides the conserved Cys/His chelating zinc ions, 6 equivalent residues (4 hydrophobic and 2 polar residues) in the SCRs possess good-consensus and conservation- these 4 SCRs function in the structural positioning of 6 equivalent residues as determinants for RING-E3 catalysis; v) members of these RING proteins located nucleus, multiple subcellular compartments, membrane protein and mitochondrion are respectively 42 (42/139, 30.2%), 71 (71/139, 51.1%), 22 (22/139, 15.8%) and 4 (4/139, 2.9%); vi) CG15104 (Topors) and CG1134 (Mul1) in C3HC4, and CG3929 (Deltex) in C3H2C3 seem to display broader E2s binding profiles than other RING-E3s; vii) analyzing intermolecular interfaces of E2/RING-E3 complexes indicate that residues directly interacting with E2s are all from the SCRs in RING domains. Of the 6 residues, 2 hydrophobic ones contribute to constructing the conserved hydrophobic core, while the 2 hydrophobic and 2 polar residues directly participate in E2/RING-E3 interactions. Based on sequence and structural data, SCRs, conserved equivalent residues and features of intermolecular interfaces were extracted, highlighting the presence of a nucleus for RING domain fold and formation of catalytic core in which related residues and regions exhibit preferential evolutionary conservation.

HIV protein sequence hotspots for crosstalk with host hub proteins

HIV proteins target host hub proteins for transient binding interactions. The presence of viral proteins in the infected cell results in out-competition of host proteins in their interaction with hub proteins, drastically affecting cell physiology. Functional genomics and interactome datasets can be used to quantify the sequence hotspots on the HIV proteome mediating interactions with host hub proteins. In this study, we used the HIV and human interactome databases to identify HIV targeted host hub proteins and their host binding partners (H2). We developed a high throughput computational procedure utilizing motif discovery algorithms on sets of protein sequences, including sequences of HIV and H2 proteins. We identified as HIV sequence hotspots those linear motifs that are highly conserved on HIV sequences and at the same time have a statistically enriched presence on the sequences of H2 proteins. The HIV protein motifs discovered in this study are expressed by subsets of H2 host proteins potentially outcompeted by HIV proteins. A large subset of these motifs is involved in cleavage, nuclear localization, phosphorylation, and transcription factor binding events. Many such motifs are clustered on an HIV sequence in the form of hotspots. The sequential positions of these hotspots are consistent with the curated literature on phenotype altering residue mutations, as well as with existing binding site data. The hotspot map produced in this study is the first global portrayal of HIV motifs involved in altering the host protein network at highly connected hub nodes.

PeptideMine--a webserver for the design of peptides for protein-peptide binding studies derived from protein-protein interactomes

Background

Signal transduction events often involve transient, yet specific, interactions between structurally conserved protein domains and polypeptide sequences in target proteins. The identification and validation of these associating domains is crucial to understand signal transduction pathways that modulate different cellular or developmental processes. Bioinformatics strategies to extract and integrate information from diverse sources have been shown to facilitate the experimental design to understand complex biological events. These methods, primarily based on information from high-throughput experiments, have also led to the identification of new connections thus providing hypothetical models for cellular events. Such models, in turn, provide a framework for directing experimental efforts for validating the predicted molecular rationale for complex cellular processes. In this context, it is envisaged that the rational design of peptides for protein-peptide binding studies could substantially facilitate the experimental strategies to evaluate a predicted interaction. This rational design procedure involves the integration of protein-protein interaction data, gene ontology, physico-chemical calculations, domain-domain interaction data and information on functional sites or critical residues.

Results

Here we describe an integrated approach called "PeptideMine" for the identification of peptides based on specific functional patterns present in the sequence of an interacting protein. This approach based on sequence searches in the interacting sequence space has been developed into a webserver, which can be used for the identification and analysis of peptides, peptide homologues or functional patterns from the interacting sequence space of a protein. To further facilitate experimental validation, the PeptideMine webserver also provides a list of physico-chemical parameters corresponding to the peptide to determine the feasibility of using the peptide for in vitro biochemical or biophysical studies.

Conclusions

The strategy described here involves the integration of data and tools to identify potential interacting partners for a protein and design criteria for peptides based on desired biochemical properties. Alongside the search for interacting protein sequences using three different search programs, the server also provides the biochemical characteristics of candidate peptides to prune peptide sequences based on features that are most suited for a given experiment. The PeptideMine server is available at the URL: http://caps.ncbs.res.in/peptidemine

Molecular dissection of human interleukin-31-mediated signal transduction through site-directed mutagenesis

Interleukin (IL)-31 is a recently described cytokine, preferentially produced by T helper 2 lymphocytes and associated with skin diseases, such as atopic dermatitis. IL-31 is a member of the four alpha-helix bundle cytokine family and is related to the IL-6 subgroup. Its heterodimeric membrane receptor is composed of the gp130-like receptor (GPL) subunit associated to the oncostatin M receptor subunit. We identified critical amino acids implicated in the ligand receptor interaction by computational analysis combined with site-directed mutagenesis. Six IL-31 residues selected for their putative involvement in cytokine receptor contact sites were alanine-substituted, and the corresponding proteins were expressed in mammalian and bacterial systems. Biochemical, membrane binding, cell signaling, and cell proliferation analyses showed that mutation E44A, E106A, or H110A abolished IL-31 binding to GPL and the subsequent signaling events. A second ligand receptor-binding site involved Lys(134), with alanine substitution leading to a protein that still binds GPL, but is unable to recruit the second receptor subunit and the subsequent signaling pathways. The results indicate that IL-31 recognizes its receptor complex through two different binding sites, and we propose a three-dimensional model for IL-31.

Finding of residues crucial for supersecondary structure formation

This work evaluates the hypothesis that proteins with an identical supersecondary structure (SSS) share a unique set of residues--SSS-determining residues--even though they may belong to different protein families and have very low sequence similarities. This hypothesis was tested on two groups of sandwich-like proteins (SPs). Proteins in each group have an identical SSS, but their sequence similarity is below the "twilight zone." To find the SSS-determining residues specific to each group, a unique structure-based algorithm of multiple sequences alignment was developed. The units of alignment are individual strands and loops rather than whole sequences. The algorithm is based on the alignment of residues that form hydrogen bonds between corresponding strands. Structure-based alignment revealed that 30-35% of the positions in the sequences in each group of proteins are "conserved positions" occupied either by hydrophobic-only or hydrophilic-only residues. Moreover, each group of SPs is characterized by a unique set of SSS-determining residues found at the conserved positions. The set of SSS-determining residues has very high sensitivity and specificity for identifying proteins with a corresponding SSS: It is an "amino acid tag" that brands a sequence as having a particular SSS. Thus, the sets of SSS-determining residues can be used to classify proteins and to predict the SSS of a query amino acid sequence.

Why there is more to protein evolution than protein function: splicing, nucleosomes and dual-coding sequence

There is considerable variation in the rate at which different proteins evolve. Why is this? Classically, it has been considered that the density of functionally important sites must predict rates of protein evolution. Likewise, amino acid choice is usually assumed to reflect optimal protein function. In the present article, we briefly review evidence suggesting that this protein function-centred view is too simplistic. In particular, we concentrate on how selection acting during the protein's production history can also affect protein evolutionary rates and amino acid choice. Exploring the role of selection at the DNA and RNA level, we specifically address how the need (i) to specify exonic splice enhancer motifs in pre-mRNA, and (ii) to ensure nucleosome positioning on DNA have an impact on amino acid choice and rates of evolution. For both, we review evidence that sequence affected by more than one coding demand is particularly constrained. Strikingly, in mammals, splicing-related constraints are quantitatively as important as expression parameters in predicting rates of protein evolution. These results indicate that there is substantially more to protein evolution than protein functional constraints.

Folding energy landscape of cytochrome cb562

Cytochrome cb(562) is a variant of an Escherichia coli four-helix bundle b-type heme protein in which the porphyrin prosthetic group is covalently ligated to the polypeptide near the terminus of helix 4. Studies from other laboratories have shown that the apoprotein folds rapidly without the formation of intermediates, whereas the holoprotein loses heme before native structure can be attained. Time-resolved fluorescence energy transfer (TRFET) measurements of cytochrome cb(562) refolding triggered using an ultrafast continuous-flow mixer (150 micros dead time) reveal that heme attachment to the polypeptide does not interfere with rapid formation of the native structure. Analyses of the TRFET data produce distributions of Trp-59-heme distances in the protein before, during, and after refolding. Characterization of the moments and time evolution of these distributions provides compelling evidence for a refolding mechanism that does not involve significant populations of intermediates. These observations suggest that the cytochrome b(562) folding energy landscape is minimally frustrated and able to tolerate the introduction of substantial perturbations (i.e., the heme prosthetic group) without the formation of deep misfolded traps.

The role of genomic data in the discovery, annotation and evolutionary interpretation of the interferon-lambda family

BACKGROUND

Type-I interferons, type-II interferons, and the IL-10 family are helical cytokines with similar three-dimensional folds. However, their homologous relationship is difficult to detect on the basis of sequence alone. We have previously described the discovery of the human type-III interferons (IFN lambda-1, -2, -3 or IL-29, IL-28A, IL-28B), which required a combination of manual and computational techniques applied to predicted protein sequences.

PRINCIPAL FINDINGS

Here we describe how the use of gene structure analysis and comparative genomics enabled a more extensive understanding of these genes early in the discovery process. More recently, additional mammalian genome sequences have shown that there are between one and potentially nine copies of interferon lambda genes in each genome, and that several species have single exon versions of the interferon lambda gene.

SIGNIFICANCE

The variable number of single exon type-I interferons in mammals, along with recently identified genes in zebrafish homologous to interferons allows a story of interferon evolution to be proposed. This model suggests that the gene duplications and single exon retrotransposons of mammalian type-III interferons are positively selected for within a genome. These characteristics are also shared with the fish interferons and could be responsible for the generation of the IL10 family and also the single exon type-I interferons.

"Hot cores" in proteins: comparative analysis of the apolar contact area in structures from hyper/thermophilic and mesophilic organisms

Background

A wide variety of stabilizing factors have been invoked so far to elucidate the structural basis of protein thermostability. These include, amongst the others, a higher number of ion-pairs interactions and hydrogen bonds, together with a better packing of hydrophobic residues. It has been frequently observed that packing of hydrophobic side chains is improved in hyperthermophilic proteins, when compared to their mesophilic counterparts. In this work, protein crystal structures from hyper/thermophilic organisms and their mesophilic homologs have been compared, in order to quantify the difference of apolar contact area and to assess the role played by the hydrophobic contacts in the stabilization of the protein core, at high temperatures.

Results

The construction of two datasets was carried out so as to satisfy several restrictive criteria, such as minimum redundancy, resolution and R-value thresholds and lack of any structural defect in the collected structures. This approach allowed to quantify with relatively high precision the apolar contact area between interacting residues, reducing the uncertainty due to the position of atoms in the crystal structures, the redundancy of data and the size of the dataset. To identify the common core regions of these proteins, the study was focused on segments that conserve a similar main chain conformation in the structures analyzed, excluding the intervening regions whose structure differs markedly. The results indicated that hyperthermophilic proteins underwent a significant increase of the hydrophobic contact area contributed by those residues composing the alpha-helices of the structurally conserved regions.

Conclusion

This study indicates the decreased flexibility of alpha-helices in proteins core as a major factor contributing to the enhanced termostability of a number of hyperthermophilic proteins. This effect, in turn, may be due to an increased number of buried methyl groups in the protein core and/or a better packing of alpha-helices with the rest of the structure, caused by the presence of hydrophobic beta-branched side chains.

The Existence of Multiple Conformers of Interleukin-21 Directs Engineering of a Superpotent Analogue

The high resolution three-dimensional structure of human interleukin (hIL)-21 has been resolved by heteronuclear NMR spectroscopy. Overall, the hIL-21 structure is dominated by a well defined central four-helical bundle, arranged in an up-up-down-down topology, as observed for other cytokines. A segment of the hIL-21 molecule that includes the third helical segment, helix C, is observed to exist in two distinct and interchangeable states. In one conformer, the helix C segment is presented in a regular, alpha-helical conformation, whereas in the other conformer, this segment is largely disordered. A structure-based sequence alignment of hIL-21 with receptor complexes of the related cytokines, interleukin-2 and -4, implied that this particular segment is involved in receptor binding. An hIL-21 analog was designed to stabilize the region around helix C through the introduction of a segment grafted from hIL-4. This novel hIL-21 analog was demonstrated to exhibit a 10-fold increase in potency in a cellular assay.

The selection of acceptable protein mutations

We have determined the general constraints that govern sequence divergence in proteins that retain entirely, or very largely, the same structure and function. To do this we collected data from three different groups of orthologous sequences: those found in humans and mice, in humans and chickens, and in Escherichia coli and Salmonella enterica. In total, these organisms have 21,738 suitable pairs of orthologs, and these contain nearly 2 million mutations. The three groups differ greatly in the taxa from which they come and/or in the time that separates them from their last common ancestor. Nevertheless, the results we obtain from the three different groups are strikingly similar. For each group, the orthologous sequence pairs were assigned to six different divergence categories on the basis of their sequence identities. For categories with the same divergence, common accepted mutations have similar frequencies and rank orders in the three groups. With divergence, the width of the range of common mutations grows in the same manner in each group. We examined the distribution of mutations in protein structures. With increasing divergence, mutations increase at different rates in the buried, intermediate, and exposed regions of protein structures in a manner that explains the exponential relationship between the divergence of structure and sequence. This work implies that commonly allowed mutations are selected by a set of general constraints that are well defined and whose nature varies with divergence.

Massive sequence perturbation of the Raf ras binding domain reveals relationships between sequence conservation, secondary structure propensity, hydrophobic core organization and stability

The contributions of specific residues to the delicate balance between function, stability and folding rates could be determined, in part by [corrected] comparing the sequences of structures having identical folds, but insignificant sequence homology. Recently, we have devised an experimental strategy to thoroughly explore residue substitutions consistent with a specific class of structure. Using this approach, the amino acids tolerated at virtually all residues of the c-Raf/Raf1 ras binding domain (Raf RBD), an exemplar of the common beta-grasp ubiquitin-like topology, were obtained and used to define the sequence determinants of this fold. Herein, we present analyses suggesting that more subtle sequence selection pressure, including propensity for secondary structure, the hydrophobic core organization and charge distribution are imposed on the Raf RBD sequence. Secondly, using the Gibbs free energies (DeltaG(F-U)) obtained for 51 mutants of Raf RBD, we demonstrate a strong correlation between amino acid conservation and the destabilization induced by truncating mutants. In addition, four mutants are shown to significantly stabilize Raf RBD native structure. Two of these mutations, including the well-studied R89L, are known to severely compromise binding affinity for ras. Another stabilized mutant consisted of a deletion of amino acid residues E104-K106. This deletion naturally occurs in the homologues a-Raf and b-Raf and could indicate functional divergence. Finally, the combination of mutations affecting five of 78 residues of Raf RBD results in stabilization of the structure by approximately 12 kJ mol(-1) (DeltaG(F-U) is -22 and -34 kJ mol(-1) for wt and mutant, respectively). The sequence perturbation approach combined with sequence/structure analysis of the ubiquitin-like fold provide a basis for the identification of sequence-specific requirements for function, stability and folding rate of the Raf RBD and structural analogues, highlighting the utility of conservation profiles as predictive tools of structural organization.

Pattern similarity study of functional sites in protein sequences: lysozymes and cystatins

BACKGROUND

Although it is generally agreed that topography is more conserved than sequences, proteins sharing the same fold can have different functions, while there are protein families with low sequence similarity. An alternative method for profile analysis of characteristic conserved positions of the motifs within the 3D structures may be needed for functional annotation of protein sequences. Using the approach of quantitative structure-activity relationships (QSAR), we have proposed a new algorithm for postulating functional mechanisms on the basis of pattern similarity and average of property values of side-chains in segments within sequences. This approach was used to search for functional sites of proteins belonging to the lysozyme and cystatin families.

RESULTS

Hydrophobicity and beta-turn propensity of reference segments with 3-7 residues were used for the homology similarity search (HSS) for active sites. Hydrogen bonding was used as the side-chain property for searching the binding sites of lysozymes. The profiles of similarity constants and average values of these parameters as functions of their positions in the sequences could identify both active and substrate binding sites of the lysozyme of Streptomyces coelicolor, which has been reported as a new fold enzyme (Cellosyl). The same approach was successfully applied to cystatins, especially for postulating the mechanisms of amyloidosis of human cystatin C as well as human lysozyme.

CONCLUSION

Pattern similarity and average index values of structure-related properties of side chains in short segments of three residues or longer were, for the first time, successfully applied for predicting functional sites in sequences. This new approach may be applicable to studying functional sites in un-annotated proteins, for which complete 3D structures are not yet available.

S4: structure-based sequence alignments of SCOP superfamilies

S4 is an automatically generated database of multiple structure-based sequence alignments of protein superfamilies in the SCOP database. All structural domains that do not share more than 40% sequence identity as defined by the ASTRAL compendium of protein structures are included. The alignments are constructed using pairwise structural alignments to generate residue equivalences that are then integrated into multiple alignments using sequence alignment tools. We describe the database and give examples showing how the automatically generated S4 alignments compare favourably to hand-crafted alignments. Available at: http://compbio.mds.qmw.ac.uk/S4.html.

Recognition of the Helical Cytokine Fold

Cytokines are an important class of proteins responsible for intercellular communication. The helical cytokines have a four-helix bundle fold, and they have remained largely intractable for sequence alignment methods due to their high evolutionary divergence. This paper presents a method that is specifically designed to recognize the helical cytokine fold in preprotein sequences such as full-length cDNA translations or transcripts predicted by gene finding methods. A protein fold is modeled by multiple sequence profiles, each representing a structurally conserved region. Nonstructural profiles are used to represent additional signals found in preprotein sequences. Profiles are connected by loop regions, each of a specified minimum and maximum length. A model for the helical cytokines is created by progressively improving a placement of four amphipathic helices onto training sequences. The sensitivity and specificity of the method are evaluated by a cross-validation procedure, demonstrating that cytokines with no intrafamily sequence similarity can be recognized. The method has been successfully used for the discovery of several new helical cytokines in the human genome.

Evolutionarily conserved regions and hydrophobic contacts at the superfamily level: The case of the fold-type I, pyridoxal-5'-phosphate-dependent enzymes

The wealth of biological information provided by structural and genomic projects opens new prospects of understanding life and evolution at the molecular level. In this work, it is shown how computational approaches can be exploited to pinpoint protein structural features that remain invariant upon long evolutionary periods in the fold-type I, PLP-dependent enzymes. A nonredundant set of 23 superposed crystallographic structures belonging to this superfamily was built. Members of this family typically display high-structural conservation despite low-sequence identity. For each structure, a multiple-sequence alignment of orthologous sequences was obtained, and the 23 alignments were merged using the structural information to obtain a comprehensive multiple alignment of 921 sequences of fold-type I enzymes. The structurally conserved regions (SCRs), the evolutionarily conserved residues, and the conserved hydrophobic contacts (CHCs) were extracted from this data set, using both sequence and structural information. The results of this study identified a structural pattern of hydrophobic contacts shared by all of the superfamily members of fold-type I enzymes and involved in native interactions. This profile highlights the presence of a nucleus for this fold, in which residues participating in the most conserved native interactions exhibit preferential evolutionary conservation, that correlates significantly (r = 0.70) with the extent of mean hydrophobic contact value of their apolar fraction.

The linked conservation of structure and function in a family of high diversity: the monomeric cupredoxins

The monomeric cupredoxins are a highly divergent family of copper binding electron transport proteins that function in photosynthesis and respiration. To determine how function and structure are conserved in the context of large sequence differences, we have carried out a detailed analysis of the cupredoxins of known structure and their sequence homologs. The common structure of the cupredoxins is formed by a sandwich of two beta sheets which support a copper binding site. The structure of the deeply buried core is intimately coupled to the binding site on the surface of the protein; in each protein the conserved regions form one continuous substructure that extends from the surface active site and through the center of the molecule. Residues around the active site are conserved for functional reasons, while those deeper in the structure will be conserved for structural reasons. Together the two sets support each other.

Structural studies of the Nudix hydrolase DR1025 from Deinococcus radiodurans and its ligand complexes

We have determined the crystal structure, at 1.4A, of the Nudix hydrolase DR1025 from the extremely radiation resistant bacterium Deinococcus radiodurans. The protein forms an intertwined homodimer by exchanging N-terminal segments between chains. We have identified additional conserved elements of the Nudix fold, including the metal-binding motif, a kinked beta-strand characterized by a proline two positions upstream of the Nudix consensus sequence, and participation of the N-terminal extension in the formation of the substrate-binding pocket. Crystal structures were also solved of DR1025 crystallized in the presence of magnesium and either a GTP analog or Ap(4)A (both at 1.6A resolution). In the Ap(4)A co-crystal, the electron density indicated that the product of asymmetric hydrolysis, ATP, was bound to the enzyme. The GTP analog bound structure showed that GTP was bound almost identically as ATP. Neither nucleoside triphosphate was further cleaved.