A revised model of the active site of alternative oxidase

The plant mitochondrial protein alternative oxidase catalyses dioxygen dependent ubiquinol oxidation to yield ubiquinone and water. A structure of this protein has previously been proposed based on an assumed structural homology to the di-iron carboxylate family of proteins. However, these authors suggested the protein has a very different topology than the known structures of di-iron carboxylate proteins. We have re-examined this model and based on comparison of recent sequences and structural data on di-iron carboxylate proteins we present a new model of the alternative oxidase which allows prediction of active site residues and a possible membrane binding motif.     

Crystal structure of a BCL-W domain-swapped dimer: implications for the function of BCL-2 family proteins

The prosurvival and proapoptotic proteins of the BCL-2 family share a similar three-dimensional fold despite their opposing functions. However, many biochemical studies highlight the requirement for conformational changes for the functioning of both types of proteins, although structural data to support such changes remain elusive. Here, we describe the X-ray structure of dimeric BCL-W that reveals a major conformational change involving helices α3 and α4 hinging away from the core of the protein. Biochemical and functional studies reveal that the α4-α5 hinge region is required for dimerization of BCL-W, and functioning of both pro- and antiapoptotic BCL-2 proteins. Hence, this structure reveals a conformational flexibility not seen in previous BCL-2 protein structures and provides insights into how these regulators of apoptosis can change conformation to exert their function.     

ECOD: An Evolutionary Classification of Protein Domains

Understanding the evolution of a protein, including both close and distant relationships, often reveals insight into its structure and function. Fast and easy access to such up-to-date information facilitates research. We have developed a hierarchical evolutionary classification of all proteins with experimentally determined spatial structures, and presented it as an interactive and updatable online database. ECOD (Evolutionary Classification of protein Domains) is distinct from other structural classifications in that it groups domains primarily by evolutionary relationships (homology), rather than topology (or “fold”). This distinction highlights cases of homology between domains of differing topology to aid in understanding of protein structure evolution. ECOD uniquely emphasizes distantly related homologs that are difficult to detect, and thus catalogs the largest number of evolutionary links among structural domain classifications. Placing distant homologs together underscores the ancestral similarities of these proteins and draws attention to the most important regions of sequence and structure, as well as conserved functional sites. ECOD also recognizes closer sequence-based relationships between protein domains. Currently, approximately 100,000 protein structures are classified in ECOD into 9,000 sequence families clustered into close to 2,000 evolutionary groups. The classification is assisted by an automated pipeline that quickly and consistently classifies weekly releases of PDB structures and allows for continual updates. This synchronization with PDB uniquely distinguishes ECOD among all protein classifications. Finally, we present several case studies of homologous proteins not recorded in other classifications, illustrating the potential of how ECOD can be used to further biological and evolutionary studies.     

Effects of point mutations on protein structure are nonexponentially distributed

Inspection of structure changes in proteins borne by altering their sequences brings understanding of physics, functioning and evolution of existing proteins, and helps engineer modified ones. On single amino acid substitutions, the most frequent mutation type, shifts in backbone conformation are typically small, raising doubts if and how such minor modifications could drive evolutionary divergence. Here, we report that the distribution of magnitudes of structure change on such substitutions is heavy-tailed--whereas protein structures are robust to most substitutions, changes much larger than average occur with raised odds compared to what would be expected for exponential distribution with the same mean. This nonexponential behavior allows for reconciling the apparent contradiction between the observed conservation of protein structures and the substantial evolutionary plasticity implied in their diversity. The presence of the heavy tail in the distribution promotes structure divergence, facilitating exploration of new functionality, and conformations within folds, as well as exploration of structure space for new folds.     

Structural and functional analysis of the middle segment of hsp90: implications for ATP hydrolysis and client protein and cochaperone interactions

Activation of client proteins by the Hsp90 molecular chaperone is dependent on binding and hydrolysis of ATP, which drives a molecular clamp via transient dimerization of the N-terminal domains. The crystal structure of the middle segment of yeast Hsp90 reveals considerable evolutionary divergence from the equivalent regions of other GHKL protein family members such as MutL and GyrB, including an additional domain of new fold. Using the known structure of the N-terminal nucleotide binding domain, a model for the Hsp90 dimer has been constructed. From this structure, residues implicated in the ATPase-coupled conformational cycle and in interactions with client proteins and the activating cochaperone Aha1 have been identified, and their roles functionally characterized in vitro and in vivo.     

Crystal structure of E. coli Hsp100 ClpB nucleotide-binding domain 1 (NBD1) and mechanistic studies on ClpB ATPase activity

E. coli Hsp100 ClpB was recently identified as a critical part in a multi-chaperone system to play important roles in protein folding, protein transport and degradation in cell physiology. ClpB contains two nucleotide-binding domains (NBD1 and NBD2) within their primary sequences. NBD1 and NBD2 of ClpB can be classified as members of the large ATPase family known as ATPases associated with various cellular activities (AAA). To investigate how ClpB performs its ATPase activities for its chaperone activity, we have determined the crystal structure of ClpB nucleotide-binding domain 1 (NBD1) by MAD method to 1.80 A resolution. The NBD1 monomer structure contains one domain that comprises 11 alpha-helices and six beta-strands. When compared with the typical AAA structures, the crystal structure of ClpB NBD1 reveals a novel AAA topology with six-stranded beta-sheet as its core. The N-terminal portion of NBD1 structure has an extra beta-strand flanked by two extra alpha-helices that are not present in other AAA structures. Moreover, the NBD1 structure does not have a C-terminal helical domain as other AAA proteins do. No nucleotide molecule is bound with ClpB NBD1 in the crystal structure probably due to lack of the C-terminal helix domain in the structure. Isothermal titration calorimetry (ITC) studies of ClpB NBD1 and other ClpB deletion mutations showed that either ClpB NBD1 or NBD2 alone does not bind to nucleotides. However, ClpB NBD2 combined with ClpB C-terminal fragment can interact with one ADP or ATP molecule. ITC data also indicated that full-length ClpB could bind two ADP molecules or one ATP analogue ATPgammaS molecule. Further ATPase activity studies of ClpB and ClpB deletion mutants showed that only wild-type ClpB have ATPase activity. None of ClpB NBD1 domain, NBD2 domain and NBD2 with C-terminal fragment has detectable ATPase activities. On the basis of our structural and mutagenesis data, we proposed a "see-saw" model to illustrate the mechanisms by which ClpB performs its ATPase activities for chaperone functions.     

Structural diversity in the Mycobacteria DUF3349 superfamily

A protein superfamily with a “Domain of Unknown Function,”, DUF3349 (PF11829), is present predominately in Mycobacterium and Rhodococcus bacterial species suggesting that these proteins may have a biological function unique to these bacteria. We previously reported the inaugural structure of a DUF3349 superfamily member, Mycobacterium tuberculosis Rv0543c. Here, we report the structures determined for three additional DUF3349 proteins: Mycobacterium smegmatis MSMEG_1063 and MSMEG_1066 and Mycobacterium abscessus MAB_3403c. Like Rv0543c, the NMR solution structure of MSMEG_1063 revealed a monomeric five α‐helix bundle with a similar overall topology. Conversely, the crystal structure of MSMEG_1066 revealed a five α‐helix protein with a strikingly different topology and a tetrameric quaternary structure that was confirmed by size exclusion chromatography. The NMR solution structure of a fourth member of the DUF3349 superfamily, MAB_3403c, with 18 residues missing at the N‐terminus, revealed a monomeric α‐helical protein with a folding topology similar to the three C‐terminal helices in the protomer of the MSMEG_1066 tetramer. These structures, together with a GREMLIN‐based bioinformatics analysis of the DUF3349 primary amino acid sequences, suggest two subfamilies within the DUF3349 family. The division of the DUF3349 into two distinct subfamilies would have been lost if structure solution had stopped with the first structure in the DUF3349 family, highlighting the insights generated by solving multiple structures within a protein superfamily. Future studies will determine if the structural diversity at the tertiary and quaternary levels in the DUF3349 protein superfamily have functional roles in Mycobacteria and Rhodococcus species with potential implications for structure‐based drug discovery.     

Sequence similarity as a predictor of the transmembrane topology of membrane-intrinsic subunits of bacterial respiratory chain enzymes

Integral membrane proteins usually have a predominantly alpha-helical secondary structure in which transmembrane segments are connected by membrane-extrinsic loops. Although a number of membrane protein structures have been reported in recent years, in most cases transmembrane topologies are initially predicted using a variety of theoretical techniques, including hydropathy analyses and the "positive inside" rule. We have explored the use of plots of the distribution of sequence similarity within families of membrane proteins comprising homeomorphic domains as a new method for the prediction/verification of the orientation of transmembrane topology models within certain families of multimeric respiratory chain enzymes. Within such proteins, analyses of sequence similarity can: i) identify heme and/or quinol binding sites; ii) identify potential electron-transfer conduits to/from prosthetic groups; and iii) locate regions defining potential subunit-subunit interactions. We mined emerging bioinformatic data for sequences of 11 families of membrane-intrinsic proteins that are part of multimeric respiratory chain complexes that also have membrane-extrinsic subunits. The sequences of each family were then aligned and the resultant alignments converted into a graphical format recording an empirical measure of the sequence similarity plotted versus residue position. In each case, this plot was compared to the predicted transmembrane topology. With one exception, there is a strong correlation between the existence     

Genomic structure of the human ezrin gene

The ERM proteins, ezrin, radixin, and moesin, act as linkers between the plasma membrane and actin cytoskeleton. They are involved in a variety of cellular functions, such as cell adhesion, migration, and the organization of cell surface structures, and are highly homologous, both in protein sequence and in functional activity, with merlin/schwannomin, a neurofibromatosis-2-associated tumor-suppressor protein. We report here the genomic structure and intron junction sequences of the human ezrin gene. Ezrin consists of 13 exons and spans approximately 24 kb genomic DNA. The coding parts of the exons range in size from 12 bp to 275 bp and the introns from 182 bp to 7 kb. The genomic structures of ezrin and moesin are highly conserved, suggesting their recent divergence. Radiation hybrid mapping has refined the location of ezrin to the interval between D6S442 and D6S281. Received: 1 June 1998 / Accepted: 25 August 1998     

Emergence of scale-free distribution in protein–protein interaction networks based on random selection of interacting domain pairs

Recent analyses of biological and artificial networks have revealed a common network architecture, called scale-free topology. The origin of the scale-free topology has been explained by using growth and preferential attachment mechanisms. In a cell, proteins are the most important carriers of function, and are composed of domains as elemental units responsible for the physical interaction between protein pairs. Here, we propose a model for protein–protein interaction networks that reveals the emergence of two possible topologies. We show that depending on the number of randomly selected interacting domain pairs, the connectivity distribution follows either a scale-free distribution, even in the absence of the preferential attachment, or a normal distribution. This new approach only requires an evolutionary model of proteins (nodes) but not for the interactions (edges). The edges are added by means of random interaction of domain pairs. As a result, this model offers a new mechanistic explanation for understanding complex networks with a direct biological interpretation because only protein structures and their functions evolved through genetic modifications of amino acid sequences. These findings are supported by numerical simulations as well as experimental data.     

Structure-based inference of molecular functions of proteins of unknown function from Berkeley Structural Genomics Center

Advances in sequence genomics have resulted in an accumulation of a huge number of protein sequences derived from genome sequences. However, the functions of a large portion of them cannot be inferred based on the current methods of sequence homology detection to proteins of known functions. Three-dimensional structure can have an important impact in providing inference of molecular function (physical and chemical function) of a protein of unknown function. Structural genomics centers worldwide have been determining many 3-D structures of the proteins of unknown functions, and possible molecular functions of them have been inferred based on their structures. Combined with bioinformatics and enzymatic assay tools, the successful acceleration of the process of protein structure determination through high throughput pipelines enables the rapid functional annotation of a large fraction of hypothetical proteins. We present a brief summary of the process we used at the Berkeley Structural Genomics Center to infer molecular functions of proteins of unknown function.     

Do proteins learn to evolve? The Hopfield network as a basis for the understanding of protein evolution

Correlations between amino-acid residues can be observed in sets of aligned protein sequences, and the analysis of their statistical and evolutionary significance and distribution has been thoroughly investigated. In this paper, we present a model based on such covariations in protein sequences in which the pairs of residues that have mutual influence combine to produce a system analogous to a Hopfield neural network. The emergent properties of such a network, such as soft failure and the connection between network architecture and stored memory, have close parallels in known proteins. This model suggests that an explanation for observed characters of proteins such as the diminution of function by substitutions distant from the active site, the existence of protein folds (superfolds) that can perform several functions based on one architecture, and structural and functional resilience to destabilizing substitutions might derive from their inherent network-like structure. This model may also provide a basis for mapping the relationship between structure, function and evolutionary history of a protein family, and thus be a powerful tool for rational engineering.     

Bioinformatics Analysis Reveals Abundant Short Alpha-Helices as a Common Structural Feature of Oomycete RxLR Effector Proteins

RxLR effectors represent one of the largest and most diverse effector families in oomycete plant pathogens. These effectors have attracted enormous attention since they can be delivered inside the plant cell and manipulates host immunity. With the exceptions of a signal peptide and the following RxLR-dEER and C-terminal W/Y/L motifs identified from the sequences themselves, nearly no functional domains have been found. Recently, protein structures of several RxLRs were revealed to comprise alpha-helical bundle repeats. However, approximately half of all RxLRs lack obvious W/Y/L motifs, which are associated with helical structures. In this study, secondary structure prediction of the putative RxLR proteins was performed. We found that the C-terminus of the majority of these RxLR proteins, irrespective of the presence of W/Y/L motifs, contains abundant short alpha-helices. Since a large-scale experimental determination of protein structures has been difficult to date, results of the current study extend our understanding on the oomycete RxLR effectors in protein secondary structures from individual members to the entire family. Moreover, we identified less alpha-helix-rich proteins from secretomes of several oomycete and fungal organisms in which RxLRs have not been identified, providing additional evidence that these organisms are unlikely to harbor RxLR-like proteins. Therefore, these results provide additional information that will aid further studies on the evolution and functional mechanisms of RxLR effectors.     

Evolutionary Dynamics on Protein Bi-stability Landscapes can Potentially Resolve Adaptive Conflicts

Experimental studies have shown that some proteins exist in two alternative native-state conformations. It has been proposed that such bi-stable proteins can potentially function as evolutionary bridges at the interface between two neutral networks of protein sequences that fold uniquely into the two different native conformations. Under adaptive conflict scenarios, bi-stable proteins may be of particular advantage if they simultaneously provide two beneficial biological functions. However, computational models that simulate protein structure evolution do not yet recognize the importance of bi-stability. Here we use a biophysical model to analyze sequence space to identify bi-stable or multi-stable proteins with two or more equally stable native-state structures. The inclusion of such proteins enhances phenotype connectivity between neutral networks in sequence space. Consideration of the sequence space neighborhood of bridge proteins revealed that bi-stability decreases gradually with each mutation that takes the sequence further away from an exactly bi-stable protein. With relaxed selection pressures, we found that bi-stable proteins in our model are highly successful under simulated adaptive conflict. Inspired by these model predictions, we developed a method to identify real proteins in the PDB with bridge-like properties, and have verified a clear bi-stability gradient for a series of mutants studied by Alexander et al. (Proc Nat Acad Sci USA 2009, 106:21149–21154) that connect two sequences that fold uniquely into two different native structures via a bridge-like intermediate mutant sequence. Based on these findings, new testable predictions for future studies on protein bi-stability and evolution are discussed.     

Structure, tissue distribution and genomic organization of the murine RRM-type RNA binding proteins TIA-1 and TIAR.

TIA-1 and TIAR are RNA binding proteins of the RNA recognition motif (RRM)/ribonucleoprotein (RNP) family that have been implicated as effectors of apoptotic cell death. We report the structures of murine TIA-1 and TIAR (mTIA-1 and mTIAR) deduced from cDNA cloning, the mRNA and protein tissue distribution of mTIA-1 and mTIAR, and the exon-intron structures of the mTIA-1 and mTIAR genes. Both mTIA-1 and mTIAR are comprised of three approximately 100 amino acid N-terminal RRM domains and a approximately 90 amino acid C-terminal auxiliary domain. This subfamily of RRM proteins is evolutionarily well conserved; mTIA-1 and mTIAR are 80% similar to each other, and 96 and 99% similar to hTIA-1 and hTIAR, respectively. The overall exon-intron structures of the mTIA-1 and mTIAR genes are also similar to each other, as well as to the human TIA-1 gene structure. While Northern blot analysis reveals that mTIA-1 and mTIAR mRNAs have a broad tissue distribution, mTIA-1 and mTIAR proteins are predominantly expressed in brain, testis and spleen. At least two isoforms of both mTIA-1 and mTIAR are generated by alternative splicing. Murine TIA-1 isoforms including or lacking the exon 5 encoded sequences are expressed at a ratio of approximately 1:1, whereas mTIAR isoforms including or lacking the 5'-end of exon 3 sequences are expressed in a approximately 1:6 ratio. Molecular characterization of murine TIA-1 and TIAR RNA binding proteins provides the basis for a genetic analysis of the functional roles of these proteins during mammalian development.     

Crystal structure of the human supernatant protein factor

Supernatant protein factor (SPF) promotes the epoxidation of squalene catalyzed by microsomes. Several studies suggest its in vivo role in the cholesterol biosynthetic pathway by a yet unknown mechanism. SPF belongs to a family of lipid binding proteins called CRAL_TRIO, which include yeast phosphatidylinositol transfer protein Sec14 and tocopherol transfer protein TTP. The crystal structure of human SPF at a resolution of 1.9 A reveals a two domain topology. The N-terminal 275 residues form a Sec14-like domain, while the C-terminal 115 residues consist of an eight-stranded jelly-roll barrel similar to that found in many viral protein structures. The ligand binding cavity has a peculiar horseshoe-like shape. Contrary to the Sec14 crystal structure, the lipid-exchange loop is in a closed conformation, suggesting a mechanism for lipid exchange.