A family-based approach reveals the function of residues in the nuclear receptor ligand-binding domain

Literature studies, 3D structure data, and a series of sequence analysis techniques were combined to reveal important residues in the structure and function of the ligand-binding domain of nuclear hormone receptors. A structure-based multiple sequence alignment allowed for the seamless combination of data from many different studies on different receptors into one single functional model. It was recently shown that a combined analysis of sequence entropy and variability can divide residues in five classes; (1) the main function or active site, (2) support for the main function, (3) signal transduction, (4) modulator or ligand binding and (5) the rest. Mutation data extracted from the literature and intermolecular contacts observed in nuclear receptor structures were analyzed in view of this classification and showed that the main function or active site residues of the nuclear receptor ligand-binding domain are involved in cofactor recruitment. Furthermore, the sequence entropy-variability analysis identified the presence of signal transduction residues that are located between the ligand, cofactor and dimer sites, suggesting communication between these regulatory binding sites. Experimental and computational results agreed well for most residues for which mutation data and intermolecular contact data were available. This allows us to predict the role of the residues for which no functional data is available yet. This study illustrates the power of family-based approaches towards the analysis of protein function, and it points out the problems and possibilities presented by the massive amounts of data that are becoming available in the "omics era". The results shed light on the nuclear receptor family that is involved in processes ranging from cancer to infertility, and that is one of the more important targets in the pharmaceutical industry.     

Solution structure of the homodimeric core domain of Escherichia coli histidine kinase EnvZ.

Escherichia coli osmosensor EnvZ is a protein histidine kinase that plays a central role in osmoregulation, a cellular adaptation process involving the His-Asp phosphorelay signal transduction system. Dimerization of the transmembrane protein is essential for its autophosphorylation and phosphorelay signal transduction functions. Here we present the NMR-derived structure of the homodimeric core domain (residues 223-289) of EnvZ that includes His 243, the site of autophosphorylation and phosphate transfer reactions. The structure comprises a four-helix bundle formed by two identical helix-turn-helix subunits, revealing the molecular assembly of two active sites within the dimeric kinase.     

FIST: a sensory domain for diverse signal transduction pathways in prokaryotes and ubiquitin signaling in eukaryotes

MOTIVATION: Sensory domains that are conserved among Bacteria, Archaea and Eucarya are important detectors of common signals detected by living cells. Due to their high sequence divergence, sensory domains are difficult to identify. We systematically look for novel sensory domains using sensitive profile-based searches initiated with regions of signal transduction proteins where no known domains can be identified by current domain models. RESULTS: Using profile searches followed by multiple sequence alignment, structure prediction and domain architecture analysis, we have identified a novel sensory domain termed FIST, which is present in signal transduction proteins from Bacteria, Archaea and Eucarya. Chromosomal proximity of FIST-encoding genes to those coding for proteins involved in amino acid metabolism and transport suggest that FIST domains bind small ligands, such as amino acids.     

Structural determinants for EB1-mediated recruitment of APC and spectraplakins to the microtubule plus end

EB1 is a member of a conserved protein family that localizes to growing microtubule plus ends. EB1 proteins also recruit cell polarity and signaling molecules to microtubule tips. However, the mechanism by which EB1 recognizes cargo is unknown. Here, we have defined a repeat sequence in adenomatous polyposis coli (APC) that binds to EB1's COOH-terminal domain and identified a similar sequence in members of the microtubule actin cross-linking factor (MACF) family of spectraplakins. We show that MACFs directly bind EB1 and exhibit EB1-dependent plus end tracking in vivo. To understand how EB1 recognizes APC and MACFs, we solved the crystal structure of the EB1 COOH-terminal domain. The structure reveals a novel homodimeric fold comprised of a coiled coil and four-helix bundle motif. Mutational analysis reveals that the cargo binding site for MACFs maps to a cluster of conserved residues at the junction between the coiled coil and four-helix bundle. These results provide a structural understanding of how EB1 binds two regulators of microtubule-based cell polarity.     

The Nod-like receptor (NLR) family: a tale of similarities and differences

Innate immunity represents an important system with a variety of vital processes at the core of many diseases. In recent years, the central role of the Nod-like receptor (NLR) protein family became increasingly appreciated in innate immune responses. NLRs are classified as part of the signal transduction ATPases with numerous domains (STAND) clade within the AAA+ ATPase family. They typically feature an N-terminal effector domain, a central nucleotide-binding domain (NACHT) and a C-terminal ligand-binding region that is composed of several leucine-rich repeats (LRRs). NLRs are believed to initiate or regulate host defense pathways through formation of signaling platforms that subsequently trigger the activation of inflammatory caspases and NF-kB. Despite their fundamental role in orchestrating key pathways in innate immunity, their mode of action in molecular terms remains largely unknown. Here we present the first comprehensive sequence and structure modeling analysis of NLR proteins, revealing that NLRs possess a domain architecture similar to the apoptotic initiator protein Apaf-1. Apaf-1 performs its cellular function by the formation of a heptameric platform, dubbed apoptosome, ultimately triggering the controlled demise of the affected cell. The mechanism of apoptosome formation by Apaf-1 potentially offers insight into the activation mechanisms of NLR proteins. Multiple sequence alignment analysis and homology modeling revealed Apaf-1-like structural features in most members of the NLR family, suggesting a similar biochemical behaviour in catalytic activity and oligomerization. Evolutionary tree comparisons substantiate the conservation of characteristic functional regions within the NLR family and are in good agreement with domain distributions found in distinct NLRs. Importantly, the analysis of LRR domains reveals surprisingly low conservation levels among putative ligand-binding motifs. The same is true for the effector domains exhibiting distinct interfaces ensuring specific interactions with downstream target proteins. All together these factors suggest specific biological functions for individual NLRs.     

Molecular modeling of single polypeptide chain of calcium-binding protein p26olf from dimeric S100B(betabeta).

P26olf from olfactory tissue of frog, which may be involved in olfactory transduction or adaptation, is a Ca2+-binding protein with 217 amino acids. The p26olf molecule contains two homologous parts consisting of the N-terminal half with amino acids 1-109 and the C-terminal half with amino acids 110-217. Each half resembles S100 protein with about 100 amino acids and contains two helix-loop-helix Ca2+-binding structural motifs known as EF-hands: a normal EF-hand at the C-terminus and a pseudo EF-hand at the N-terminus. Multiple alignment of the two S100-like domains of p26olf with 18 S100 proteins indicated that the C-terminal putative EF-hand of each domain contains a four-residue insertion when compared with the typical EF-hand motifs in the S100 protein, while the N-terminal EF-hand is homologous to its pseudo EF-hand. We constructed a three-dimensional model of the p26olf molecule based on results of the multiple alignment and NMR structures of dimeric S100B(betabeta) in the Ca2+-free state. The predicted structure of the p26olf single polypeptide chain satisfactorily adopts a folding pattern remarkably similar to dimeric S100B(betabeta). Each domain of p26olf consists of a unicornate-type four-helix bundle and they interact with each other in an antiparallel manner forming an X-type four-helix bundle between the two domains. The two S100-like domains of p26olf are linked by a loop with no steric hindrance, suggesting that this loop might play an important role in the function of p26olf. The circular dichroism spectral data support the predicted structure of p26olf and indicate that Ca2+-dependent conformational changes occur. Since the C-terminal putative EF-hand of each domain fully keeps the helix-loop-helix motif having a longer Ca2+-binding loop, regardless of the four-residue insertion, we propose that it is a new, novel EF-hand, although it is unclear whether this EF-hand binds Ca2+. P26olf is a new member of the S100 protein family.     

DNA polymerase beta-like nucleotidyltransferase superfamily: identification of three new families, classification and evolutionary history

A detailed analysis of the polbeta superfamily of nucleotidyltransferases was performed using computer methods for iterative database search, multiple alignment, motif analysis and structural modeling. Three previously uncharacterized families of predicted nucleotidyltransferases are described. One of these new families includes small proteins found in all archaea and some bacteria that appear to consist of the minimal nucleotidyltransferase domain and may resemble the ancestral state of this superfamily. Another new family that is specifically related to eukaryotic polyA polymerases is typified by yeast Trf4p and Trf5p proteins that are involved in chromatin remodeling. The TRF family is represented by multiple members in all eukaryotes and may be involved in yet unknown nucleotide polymerization reactions required for maintenance of chromatin structure. Another new family of bacterial and archaeal nucleotidyltransferases is predicted to function in signal transduction since, in addition to the nucleotidyltransferase domain, these proteins contain ligand-binding domains. It is further shown that the catalytic domain of gamma proteobacterial adenylyl cyclases is homologous to the polbeta superfamily nucleotidyltransferases which emphasizes the general trend for the origin of signal-transducing enzymes from those involved in replication, repair and RNA processing. Classification of the polbeta superfamily into distinct families and examination of their phyletic distribution suggests that the evolution of this type of nucleotidyltransferases may have included bursts of rapid divergence linked to the emergence of new functions as well as a number of horizontal gene transfer events.     

Structural homology between rbs repressor and ribose binding protein implies functional similarity.

The deduced amino acid sequence of the rbs repressor, RbsR, of Escherichia coli is homologous over its C-terminal 272 residues to the entire sequence of the periplasmic ribose binding protein. RbsR is also homologous to a family of bacterial repressor proteins including LacI. This implies that the structure of the repressor consists of a two-domain binding protein portion attached to a DNA-binding domain having the four-helix structure of the LacI headpiece. The implications of these relationships to the mechanism of this class of repressors are discussed.     

Refining structural and functional predictions for secretasome components by comparative sequence analysis

Comparative sequence analysis of presenilins reveals the conserved transmembrane domain shared with leukocyte antigen CD47, possibly involved in signal transduction. Sensitive techniques of multiple sequence alignment extend the earlier observation of the aminopeptidase homology domain in nicastrin to suggest that this protein may be a catalytically active component of secretasome involved in proteolysis or co-proteolysis of presenilin or beta-amyloid.     

Rotational on-off switching of a hybrid membrane sensor kinase Tar-ArcB in Escherichia coli

Signal transduction in biological systems typically involves receptor proteins that possess an extracytosolic sensory domain connected to a cytosolic catalytic domain. Relatively little is known about the mechanism by which the signal is transmitted from the sensory site to the catalytic site. At least in the case of Tar (methyl-accepting chemotaxis protein for sensing aspartate) of Escherichia coli, vertical piston-like displacements of one transmembrane segment relative to the other within the monomer induced by ligand binding has been shown to modulate the catalytic activity of the cytosolic domain. The ArcB sensor kinase of E. coli is a transmembrane protein without a significant periplasmic domain. Here, we explore how the signal is conveyed to the catalytic site by analyzing the property of various Tar-ArcB hybrids. Our results suggest that, in contrast to the piston-like displacement that operates in Tar, the catalytic activity of ArcB is set by altering the orientation of the cytosolic domain of one monomer relative to the other in the homodimer. Thus, ArcB represents a distinct family of membrane receptor proteins whose catalytic activity is determined by rotational movements of the cytosolic domain.     

Aer on the inside looking out: paradigm for a PAS-HAMP role in sensing oxygen, redox and energy

Aer, the Escherichia coli aerotaxis (oxygen-sensing) receptor, is representative of a small class of receptors that face the cytoplasm in bacteria. Instead of sensing oxygen directly, Aer detects redox changes in the electron transport system or cytoplasm when the bacteria enter or leave a hypoxic microniche. As a result, Aer sensing also enables bacteria to avoid environments where carbon deficiency, unfavourable reduction potential or other insults would limit energy production. An FAD-binding PAS domain is the sensor for Aer and a HAMP domain interacts with the PAS domain to form an input-output module for signal transduction. By analogy to the first solution structure of an isolated HAMP domain from Archaeoglobus, Aer HAMP is proposed to fold into a four-helix bundle that rotates between a signal-on and signal-off conformation. Aer is the first protein in which a PAS-HAMP input-output module has been investigated. The structure and signal transduction mechanism of Aer is providing important insights into signalling by PAS and HAMP domains.     

Evolutionary Analysis for Functional Divergence of the Toll-Like Receptor Gene Family and Altered Functional Constraints

The Toll-like receptor (TLR) gene family consists of type 1 transmembrane receptors, which play essential roles in both innate immunity and adaptive immune response by ligand recognition and signal transduction. Using all available vertebrate TLR protein sequences, we inferred the phylogenetic tree and then characterized critical amino acid residues for functional divergence by detecting altered functional constraints after gene duplications. We found that the extracellular domain of TLR genes showed higher functional divergence than that of the cytoplasmic domain, particularly in the region between leucine-rich repeat (LRR) 10 and LRR 15 of TLR 4. Our finding supports the concept that sequence evolution in the extracellular domain may be responsible for the broad diversity of TLR ligand-binding affinity, providing a testable hypothesis for potential targets that could be verified by further experimentation.     

Different conformations of the kinase-on and kinase-off signaling states in the Aer HAMP domain

HAMP domains are sensory transduction modules that connect input and output domains in diverse signaling proteins from archaea, bacteria, and lower eukaryotes. Here, we employed in vivo disulfide cross-linking to explore the structure of the HAMP domain in the Escherichia coli aerotaxis receptor Aer. Using an Aer HAMP model based on the structure of Archaeoglobus fulgidus Af1503-HAMP, the closest residue pairs at the interface of the HAMP AS-1 and AS-2' helices were determined and then replaced with cysteines and cross-linked in vivo. Except for a unique discontinuity in AS-2, the data suggest that the Aer HAMP domain forms a parallel four-helix bundle that is similar to the structure of Af1503. The HAMP discontinuity was associated with a segment of AS-2 that was recently shown to interact with the Aer-PAS sensing domain. The four-helix HAMP bundle and its discontinuity were maintained in both the kinase-on and kinase-off states of Aer, although differences in the rates of disulfide formation also indicated the existence of different HAMP conformations in the kinase-on and kinase-off states. In particular, the kinase-on state was accompanied by significantly increased disulfide formation rates at the distal end of the HAMP four-helix bundle. This indicates that HAMP signaling may be associated with a tilting of the AS-1 and AS-2' helices, which may be the signal that is transmitted to the kinase control region of Aer.     

The ubiquitin domain superfold: structure-based sequence alignments and characterization of binding epitopes

Ubiquitin-like domains are present, apart from ubiquitin-like proteins themselves, in many multidomain proteins involved in different signal transduction processes. The sequence conservation for all ubiquitin superfold family members is rather poor, even between subfamily members, leading to mistakes in sequence alignments using conventional sequence alignment methods. However, a correct alignment is essential, especially for in silico methods that predict binding partners on the basis of sequence and structure. In this study, using 3D-structural information we have generated and manually corrected sequence alignments for proteins of the five ubiquitin superfold subfamilies. On the basis of this alignment, we suggest domains for which structural information will be useful to allow homology modelling. In addition, we have analysed the energetic and electrostatic properties of ubiquitin-like domains in complex with various functional binding proteins using the protein design algorithm FoldX. On the basis of an in silico alanine-scanning mutagenesis, we provide a detailed binding epitope mapping of the hotspots of the ubiquitin domain fold, involved in the interaction with different domains and proteins. Finally, we provide a consensus fingerprint sequence that identifies all sequences described to belong to the ubiquitin superfold family. It is possible that the method that we describe may be applied to other domain families sharing a similar fold but having low levels of sequence homology.     

Insights into eukaryotic multistep phosphorelay signal transduction revealed by the crystal structure of Ypd1p from Saccharomyces cerevisiae.

"Two-component" phosphorelay signal transduction systems constitute a potential target for antibacterial and antifungal agents, since they are found exclusively in prokaryotes and lower eukaryotes (yeast, fungi, slime mold, and plants) but not in mammalian organisms. Saccharomyces cerevisiae Ypd1p, a key intermediate in the osmosensing multistep phosphorelay signal transduction, catalyzes the phosphoryl group transfer between response regulators. Its 1.8 A structure, representing the first example of a eukaryotic phosphorelay protein, contains a four-helix bundle as in the HPt domain of Escherichia coli ArcB sensor kinase. However, Ypd1p has a 44-residue insertion between the last two helices of the helix bundle. The side-chain of His64, the site of phosphorylation, protrudes into the solvent. The structural resemblance between Ypd1p and ArcB HPt domain suggests that both prokaryotes and lower eukaryotes utilize the same basic protein fold for phosphorelay signal transduction. This study sheds light on the best characterized eukaryotic phosphorelay system.     

A model for the nucleotide-binding domains of ABC transporters based on the large domain of aspartate aminotransferase

ABC transporters are a large superfamily of integral membrane proteins involved in ATP-dependent transport across biological membranes. Members of this superfamily play roles in a number of phenomena of biomedical interest, including cystic fibrosis (CFTR) and multidrug resistance (P-glycoprotein, MRP). Most ABC transporters are predicted to consist of four domains, two membrane-spanning domains and two cytoplasmic domains. The latter contain conserved nucleotide-binding motifs. Attempts to determine the structure of ABC transporters and of their separate domains are in progress but have not yet been successful. To aid structure determination and possibly learn more about the domain boundaries, we set out to model nucleotide-binding domains (NBDs) of ABC transporters based on a known structure. Previous attempts to predict the 3D structure of NBDs were based solely on sequence similarity with known nucleotide-binding folds. We have analyzed the sequences of a number of nucleotide-binding domains with the algorithm THREADER, developed by D.T. Jones, and a possible fold was found in the structure of aspartate aminotransferase. We present a model for the N-terminal NBD of CFTR, based on the large domain of the A chain of aspartate aminotransferase. The model is refined using multiple sequence alignment, secondary structure prediction, and 3D-1D profiles. Our model seems to be in good agreement with known properties of nucleotide-binding domains and has some appealing characteristics compared with the previous models. Proteins 30:275–286, 1998. © 1998 Wiley-Liss, Inc.     

Three-dimensional structural model of the serine receptor ligand-binding domain.

Computer-based homology modeling techniques were used to construct a three-dimensional model of the Escherichia coli serine receptor ligand-binding domain based on the crystal structure of the Salmonella typhimurium aspartate receptor and the sequence homology between the two receptors. Residues that have been found in mutagenesis studies to be necessary for serine binding are located in a proposed serine-binding site. Several other mutations that affect swimming behavior require relatively small shifts in alpha-carbon positions in the model to give a minimized structure, suggesting that small changes in receptor conformation can affect the signaling state of the receptor.     

Sequence annotation of nuclear receptor ligand-binding domains by automated homology modeling

The quality of three-dimensional homology models derived from protein sequences provides an independent measure of the suitability of a protein sequence for a certain fold. We have used automated homology modeling and model assessment tools to identify putative nuclear hormone receptor ligand-binding domains in the genome of Caenorhabditis elegans. Our results indicate that the availability of multiple crystal structures is crucial to obtaining useful models in this receptor family. The majority of annotated mammalian nuclear hormone receptors could be assigned to a ligand-binding domain fold by using the best model derived from any of four template structures. This strategy also assigned the ligand-binding domain fold to a number of C.elegans. sequences without prior annotation. Interestingly, the retinoic acid receptor crystal structure contributed most to the number of sequences that could be assigned to a ligand-binding domain fold. Several causes for this can be suggested, including the high quality of this protein structure in terms of our assessment tools, similarity between the biological function or ligand of this receptor and the modeled genes and gene duplication in C.elegans.     

COG4849 is a novel family of nucleotidyltransferases

Only a minority of currently known protein families is characterized structurally. This makes homology-based structure modeling an essential instrument that can be viewed as the first approximation to experimental determination of protein structure. Using sequence similarity searches, we detected a distant similarity between a family of uncharacterized hypothetical proteins, COG4849, and the family of tRNA nucleotidyltransferases. The suggested remote homology between the N-terminal domain of COG4849 and the catalytic domain of tRNA nucleotidyltransferase was further supported by comparison of sequence profiles, methods for fold recognition and structure modeling. The combined multiple alignment of the two families reveals shared conservation of functionally important motifs and suggests the similarity in catalytic mechanisms of the performed reactions. Our results suggest that (i) the N-terminal domain of proteins from COG4849 shares structural similarity with the catalytic domain of tRNA nucleotidyltransferase, and (ii) this domain catalyzes the nucleotidyl transfer reaction involving two metal ions.