Characterizations of diverse residue clusters in protein three-dimensional structures.

We present new methods for identifying and analyzing statistically significant residue clusters that occur in three-dimensional (3D) protein structures. Residue clusters of different kinds occur in many contexts. They often feature the active site (e.g., in substrate binding), the interface between polypeptide units of protein complexes, regions of protein-protein and protein-nucleic acid interactions, or regions of metal ion coordination. The methods are illustrated with 3D clusters centering on four themes. (i) Acidic or histidine-acidic clusters associated with metal ions. (ii) Cysteine clusters including coordination of metals such as zinc or iron-sulfur structures, cysteine knots prominent in growth factors, multiple sets of buried disulfide pairings that putatively nucleate the hydrophobic core, or cysteine clusters of mostly exposed disulfide bridges. (iii) Iron-sulfur proteins and charge clusters. (iv) 3D environments of multiple histidine residues. Study of diverse 3D residue clusters offers a new perspective on protein structure and function. The algorithms can aid in rapid identification of distinctive sites, suggest correlations among protein structures, and serve as a tool in the analysis of new structures.

Large-scale protein structure modeling of the Saccharomyces cerevisiae genome

The function of a protein generally is determined by its three-dimensional (3D) structure. Thus, it would be useful to know the 3D structure of the thousands of protein sequences that are emerging from the many genome projects. To this end, fold assignment, comparative protein structure modeling, and model evaluation were automated completely. As an illustration, the method was applied to the proteins in the Saccharomyces cerevisiae (baker’s yeast) genome. It resulted in all-atom 3D models for substantial segments of 1,071 (17%) of the yeast proteins, only 40 of which have had their 3D structure determined experimentally. Of the 1,071 modeled yeast proteins, 236 were related clearly to a protein of known structure for the first time; 41 of these previously have not been characterized at all.

The RESID Database of protein structure modifications and the NRL-3D Sequence–Structure Database

The RESID Database is a comprehensive collection of annotations and structures for protein post-translational modifications including N-terminal, C-terminal and peptide chain cross-link modifications. The RESID Database includes systematic and frequently observed alternate names, Chemical Abstracts Service registry numbers, atomic formulas and weights, enzyme activities, taxonomic range, keywords, literature citations with database cross-references, structural diagrams and molecular models. The NRL-3D Sequence–Structure Database is derived from the three-dimensional structure of proteins deposited with the Research Collaboratory for Structural Bioinformatics Protein Data Bank. The NRL-3D Database includes standardized and frequently observed alternate names, sources, keywords, literature citations, experimental conditions and searchable sequences from model coordinates. These databases are freely accessible through the National Cancer Institute–Frederick Advanced Biomedical Computing Center at these web sites: http://www.ncifcrf.gov/RESID, http://www.ncifcrf.gov/ NRL-3D; or at these National Biomedical Research Foundation Protein Information Resource web sites: http://pir.georgetown.edu/pirwww/dbinfo/resid.html, http://pir.georgetown.edu/pirwww/dbinfo/nrl3d.html

Photoactive yellow protein: A structural prototype for the three-dimensional fold of the PAS domain superfamily

PAS domains are found in diverse proteins throughout all three kingdoms of life, where they apparently function in sensing and signal transduction. Although a wealth of useful sequence and functional information has become recently available, these data have not been integrated into a three-dimensional (3D) framework. The very early evolutionary development and diverse functions of PAS domains have made sequence analysis and modeling of this protein superfamily challenging. Limited sequence similarities between the ∼50-residue PAS repeats and one region of the bacterial blue-light photosensor photoactive yellow protein (PYP), for which ground-state and light-activated crystallographic structures have been determined to high resolution, originally were identified in sequence searches using consensus sequence probes from PAS-containing proteins. Here, we found that by changing a few residues particular to PYP function, the modified PYP sequence probe also could select PAS protein sequences. By mapping a typical ∼150-residue PAS domain sequence onto the entire crystallographic structure of PYP, we show that the PAS sequence similarities and differences are consistent with a shared 3D fold (the PAS/PYP module) with obvious potential for a ligand-binding cavity. Thus, PYP appears to prototypically exhibit all the major structural and functional features characteristic of the PAS domain superfamily: the shared PAS/PYP modular domain fold of ∼125–150 residues, a sensor function often linked to ligand or cofactor (chromophore) binding, and signal transduction capability governed by heterodimeric assembly (to the downstream partner of PYP). This 3D PAS/PYP module provides a structural model to guide experimental testing of hypotheses regarding ligand-binding, dimerization, and signal transduction.

The design, computer modeling, solution structure, and biological evaluation of synthetic analogs of bryostatin 1

The bryostatins are a unique family of emerging cancer chemotherapeutic candidates isolated from marine bryozoa. Although the biochemical basis for their therapeutic activity is not known, these macrolactones exhibit high affinities for protein kinase C (PKC) isozymes, compete for the phorbol ester binding site on PKC, and stimulate kinase activity in vitro and in vivo. Unlike the phorbol esters, they are not first-stage tumor promoters. The design, computer modeling, NMR solution structure, PKC binding, and functional assays of a unique class of synthetic bryostatin analogs are described. These analogs (7b, 7c, and 8) retain the putative recognition domain of the bryostatins but are simplified through deletions and modifications in the C4-C14 spacer domain. Computer modeling of an analog prototype (7a) indicates that it exists preferentially in two distinct conformational classes, one in close agreement with the crystal structure of bryostatin 1. The solution structure of synthetic analog 7c was determined by NMR spectroscopy and found to be very similar to the previously reported structures of bryostatins 1 and 10. Analogs 7b, 7c, and 8 bound strongly to PKC isozymes with Ki = 297, 3.4, and 8.3 nM, respectively. Control 7d, like the corresponding bryostatin derivative, exhibited weak PKC affinity, as did the derivative, 9, lacking the spacer domain. Like bryostatin, acetal 7c exhibited significant levels of in vitro growth inhibitory activity (1.8–170 ng/ml) against several human cancer cell lines, providing an important step toward the development of simplified, synthetically accessible analogs of the bryostatins.

Molecular cloning and characterization of an invertebrate cellular retinoic acid binding protein

We have cloned a cDNA and gene from the tobacco hornworm, Manduca sexta, which is related to the vertebrate cellular retinoic acid binding proteins (CRABPs). CRABPs are members of the superfamily of lipid binding proteins (LBPs) and are thought to mediate the effects of retinoic acid (RA) on morphogenesis, differentiation, and homeostasis. This discovery of a Manduca sexta CRABP (msCRABP) demonstrates the presence of a CRABP in invertebrates. Compared with bovine/murine CRABP I, the deduced amino acid sequence of msCRABP is 71% homologous overall and 88% homologous for the ligand binding pocket. The genomic organization of msCRABP is conserved with other CRABP family members and the larger LBP superfamily. Importantly, the promoter region contains a motif that resembles an RA response element characteristic of the promoter region of most CRABPs analyzed. Three-dimensional molecular modeling based on postulated structural homology with bovine/murine CRABP I shows msCRABP has a ligand binding pocket that can accommodate RA. The existence of an invertebrate CRABP has significant evolutionary implications, suggesting CRABPs appeared during the evolution of the LBP superfamily well before vertebrate/invertebrate divergence, instead of much later in evolution in selected vertebrates.

Protein thermostability above 100°C: A key role for ionic interactions

The discovery of hyperthermophilic microorganisms and the analysis of hyperthermostable enzymes has established the fact that multisubunit enzymes can survive for prolonged periods at temperatures above 100°C. We have carried out homology-based modeling and direct structure comparison on the hexameric glutamate dehydrogenases from the hyperthermophiles Pyrococcus furiosus and Thermococcus litoralis whose optimal growth temperatures are 100°C and 88°C, respectively, to determine key stabilizing features. These enzymes, which are 87% homologous, differ 16-fold in thermal stability at 104°C. We observed that an intersubunit ion-pair network was substantially reduced in the less stable enzyme from T. litoralis, and two residues were then altered to restore these interactions. The single mutations both had adverse effects on the thermostability of the protein. However, with both mutations in place, we observed a fourfold improvement of stability at 104°C over the wild-type enzyme. The catalytic properties of the enzymes were unaffected by the mutations. These results suggest that extensive ion-pair networks may provide a general strategy for manipulating enzyme thermostability of multisubunit enzymes. However, this study emphasizes the importance of the exact local environment of a residue in determining its effects on stability.

Genetic definition of a protein-splicing domain: Functional mini-inteins support structure predictions and a model for intein evolution

Inteins are protein-splicing elements, most of which contain conserved sequence blocks that define a family of homing endonucleases. Like group I introns that encode such endonucleases, inteins are mobile genetic elements. Recent crystallography and computer modeling studies suggest that inteins consist of two structural domains that correspond to the endonuclease and the protein-splicing elements. To determine whether the bipartite structure of inteins is mirrored by the functional independence of the protein-splicing domain, the entire endonuclease component was deleted from the Mycobacterium tuberculosis recA intein. Guided by computer modeling studies, and taking advantage of genetic systems designed to monitor intein function, the 440-aa Mtu recA intein was reduced to a functional mini-intein of 137 aa. The accuracy of splicing of several mini-inteins was verified. This work not only substantiates structure predictions for intein function but also supports the hypothesis that, like group I introns, mobile inteins arose by an endonuclease gene invading a sequence encoding a small, functional splicing element.

RNA-controlled polymorphism in the in vivo assembly of 180-subunit and 120-subunit virions from a single capsid protein

Repeated, specific interactions between capsid protein (CP) subunits direct virus capsid assembly and exemplify regulated protein–protein interactions. The results presented here reveal a striking in vivo switch in CP assembly. Using cryoelectron microscopy, three-dimensional image reconstruction, and molecular modeling, we show that brome mosaic virus (BMV) CP can assemble in vivo two remarkably distinct capsids that selectively package BMV-derived RNAs in the absence of BMV RNA replication: a 180-subunit capsid indistinguishable from virions produced in natural infections and a previously unobserved BMV capsid type with 120 subunits arranged as 60 CP dimers. Each such dimer contains two CPs in distinct, nonequivalent environments, in contrast to the quasi-equivalent CP environments throughout the 180-subunit capsid. This 120-subunit capsid utilizes most of the CP interactions of the 180-subunit capsid plus nonequivalent CP–CP interactions. Thus, the CP of BMV, and perhaps other viruses, can encode CP–CP interactions that are not apparent from mature virions and may function in assembly or disassembly. Shared structural features suggest that the 120- and 180-subunit capsids share assembly steps and that a common pentamer of CP dimers may be an important assembly intermediate. The ability of a single CP to switch between distinct capsids by means of alternate interactions also implies reduced evolutionary barriers between different capsid structures. The in vivo switch between alternate BMV capsids is controlled by the RNA packaged: a natural BMV genomic RNA was packaged in 180-subunit capsids, whereas an engineered mRNA containing only the BMV CP gene was packaged in 120-subunit capsids. RNA features can thus direct the assembly of a ribonucleoprotein complex between alternate structural pathways.

Key residues revealed in a major conformational epitope of the U1-70K protein

Epitopes depending on three-dimensional folding of proteins have during recent years been acknowledged to be main targets for many autoantibodies. However, a detailed resolution of conformation-dependent epitopes has to date not been achieved in spite of its importance for understanding the complex interaction between an autoantigen and the immune system. In analysis of immunodominant epitopes of the U1-70K protein, the major autoantigen recognized by human ribonucleoprotein (RNP)-positive sera, we have used diversely mutated recombinant Drosophila melanogaster 70K proteins as antigens in assays for human anti-RNP antibodies. Thus, the contribution of individual amino acids to antigenicity could be assayed with the overall structure of the major antigenic domain preserved, and analysis of how antigenicity can be reconstituted rather than obliterated was enabled. Our results reveal that amino acid residue 125 is situated at a crucial position for recognition by human anti-RNP autoantibodies and that flanking residues at positions 119–126 also appear to be of utmost importance for recognition. These results are discussed in relation to structural models of RNA-binding domains, and tertiary structure modeling indicates that the residues 119–126 are situated at easily accessible positions in the end of an α-helix in the RNA binding region. This study identifies a major conformation-dependent epitope of the U1-70K protein and demonstrates the significance of individual amino acids in conformational epitopes. Using this model, we believe it will be possible to analyze other immunodominant regions in which protein conformation has a strong impact.

A functional genetic screen identifies regions at the C-terminal tail and death-domain of death-associated protein kinase that are critical for its proapoptotic activity

Death-associated protein kinase (DAP-kinase) is a Ca+2/calmodulin-regulated serine/threonine kinase with a multidomain structure that participates in apoptosis induced by a variety of signals. To identify regions in this protein that are critical for its proapoptotic activity, we performed a genetic screen on the basis of functional selection of short DAP-kinase-derived fragments that could protect cells from apoptosis by acting in a dominant-negative manner. We expressed a library of randomly fragmented DAP-kinase cDNA in HeLa cells and treated these cells with IFN-γ to induce apoptosis. Functional cDNA fragments were recovered from cells that survived the selection, and those in the sense orientation were examined further in a secondary screen for their ability to protect cells from DAP-kinase-dependent tumor necrosis factor-α-induced apoptosis. We isolated four biologically active peptides that mapped to the ankyrin repeats, the “linker” region, the death domain, and the C-terminal tail of DAP-kinase. Molecular modeling of the complete death domain provided a structural basis for the function of the death-domain-derived fragment by suggesting that the protective fragment constitutes a distinct substructure. The last fragment, spanning the C-terminal serine-rich tail, defined a new regulatory region. Ectopic expression of the tail peptide (17 amino acids) inhibited the function of DAP-kinase, whereas removal of this region from the complete protein caused enhancement of the killing activity, indicating that the C-terminal tail normally plays a negative regulatory role. Altogether, this unbiased screen highlighted functionally important regions in the protein and revealed an additional level of regulation of DAP-kinase apoptotic function that does not affect the catalytic activity.

Determination of the binding constants of the centromere protein Cbf1 to all 16 centromere DNAs of Saccharomyces cerevisiae

Cbf1p is a Saccharomyces cerevisiae chromatin protein belonging to the basic region helix–loop–helix leucine zipper (bHLHzip) family of DNA binding proteins. Cbf1p binds to a conserved element in the 5′-flanking region of methionine biosynthetic genes and to centromere DNA element I (CDEI) of S.cerevisiae centromeric DNA. We have determined the apparent equilibrium dissociation constants of Cbf1p binding to all 16 CDEI DNAs in gel retardation assays. Binding constants of full-length Cbf1p vary between 1.7 and 3.8 nM. However, the dissociation constants of a Cbf1p deletion variant that has been shown to be fully sufficient for Cbf1p function in vivo vary in a range between 3.2 and 12 nM. In addition, native polyacrylamide gel electrophoresis revealed distinct changes in the 3D structure of the Cbf1p/CEN complexes. We also show that the previously reported DNA binding stimulation activity of the centromere protein p64 functions on both the Cbf1 full-length protein and a deletion variant containing only the bHLHzip domain of Cbf1p. Our results suggest that centromeric DNA outside the consensus CDEI sequence and interaction of Cbf1p with adjacent centromere proteins contribute to the complex formation between Cbf1p and CEN DNA.

Site-specific mutations in the mature form of human IL-18 with enhanced biological activity and decreased neutralization by IL-18 binding protein

IL-18 can be considered a proinflammatory cytokine mediating disease as well as an immunostimulatory cytokine that is important for host defense against infection and cancer. The high-affinity, constitutively expressed, and circulating IL-18 binding protein (IL-18BP), which competes with cell surface receptors for IL-18 and neutralizes IL-18 activity, may act as a natural antiinflammatory as well as immunosuppressive molecule. In the present studies, the IL-18 precursor caspase-1 cleavage site was changed to a factor Xa site, and, after expression in Escherichia coli, mature IL-18 was generated by factor Xa cleavage. Mature IL-18 generated by factor Xa cleavage was fully active. Single point mutations in the mature IL-18 peptide were made, and the biological activities of the wild-type (WT) IL-18 were compared with those of the mutants. Mutants E42A and K89A exhibited 2-fold increased activity compared with WT IL-18. A double mutant, E42A plus K89A, exhibited 4-fold greater activity. Unexpectedly, IL-18BP failed to neutralize the double mutant E42A plus K89A compared with WT IL-18. The K89A mutant was intermediate in being neutralized by IL-18BP, whereas neutralization of the E42A mutant was comparable to that in the WT IL-18. The identification of E42 and K89 in the mature IL-18 peptide is consistent with previous modeling studies of IL-18 binding to IL-18BP and explains the unusually high affinity of IL-18BP for IL-18.

X-ray structure of the human hyperplastic discs protein: An ortholog of the C-terminal domain of poly(A)-binding protein

The poly(A)-binding protein (PABP) recognizes the 3′ mRNA poly(A) tail and plays an essential role in eukaryotic translation initiation and mRNA stabilization/degradation. PABP is a modular protein, with four N-terminal RNA-binding domains and an extensive C terminus. The C-terminal region of PABP is essential for normal growth in yeast and has been implicated in mediating PABP homo-oligomerization and protein–protein interactions. A small, proteolytically stable, highly conserved domain has been identified within this C-terminal segment. Remarkably, this domain is also present in the hyperplastic discs protein (HYD) family of ubiquitin ligases. To better understand the function of this conserved region, an x-ray structure of the PABP-like segment of the human HYD protein has been determined at 1.04-Å resolution. The conserved domain adopts a novel fold resembling a right-handed supercoil of four α-helices. Sequence profile searches and comparative protein structure modeling identified a small ORF from the Arabidopsis thaliana genome that encodes a structurally similar but distantly related PABP/HYD domain. Phylogenetic analysis of the experimentally determined (HYD) and homology modeled (PABP) protein surfaces revealed a conserved feature that may be responsible for binding to a PABP interacting protein, Paip1, and other shared interaction partners.

Role of Rab3 GDP/GTP Exchange Protein in Synaptic Vesicle Trafficking at the Mouse Neuromuscular Junction

The Rab3 small G protein family consists of four members, Rab3A, -3B, -3C, and -3D. Of these members, Rab3A regulates Ca2+-dependent neurotransmitter release. These small G proteins are activated by Rab3 GDP/GTP exchange protein (Rab3 GEP). To determine the function of Rab3 GEP during neurotransmitter release, we have knocked out Rab3 GEP in mice. Rab3 GEP−/− mice developed normally but died immediately after birth. Embryos at E18.5 showed no evoked action potentials of the diaphragm and gastrocnemius muscles in response to electrical stimulation of the phrenic and sciatic nerves, respectively. In contrast, axonal conduction of the spinal cord and the phrenic nerve was not impaired. Total numbers of synaptic vesicles, especially those docked at the presynaptic plasma membrane, were reduced at the neuromuscular junction ∼10-fold compared with controls, whereas postsynaptic structures and functions appeared normal. Thus, Rab3 GEP is essential for neurotransmitter release and probably for formation and trafficking of the synaptic vesicles.