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.

YPD™, PombePD™ and WormPD™: model organism volumes of the BioKnowledge™ Library, an integrated resource for protein information

The BioKnowledge Library is a relational database and web site (http://www.proteome.com) composed of protein-specific information collected from the scientific literature. Each Protein Report on the web site summarizes and displays published information about a single protein, including its biochemical function, role in the cell and in the whole organism, localization, mutant phenotype and genetic interactions, regulation, domains and motifs, interactions with other proteins and other relevant data. This report describes four species-specific volumes of the BioKnowledge Library, concerned with the model organisms Saccharo­myces cerevisiae (YPD), Schizosaccharomyces pombe (PombePD) and Caenorhabditis elegans (WormPD), and with the fungal pathogen Candida albicans (CalPD™). Protein Reports of each species are unified in format, easily searchable and extensively cross-referenced between species. The relevance of these comprehensively curated resources to analysis of proteins in other species is discussed, and is illustrated by a survey of model organism proteins that have similarity to human proteins involved in disease.

Structural Analysis and Molecular Model of a Self-Incompatibility RNase from Wild Tomato1

Self-incompatibility RNases (S-RNases) are an allelic series of style glycoproteins associated with rejection of self-pollen in solanaceous plants. The nucleotide sequences of S-RNase alleles from several genera have been determined, but the structure of the gene products has only been described for those from Nicotiana alata. We report on the N-glycan structures and the disulfide bonding of the S3-RNase from wild tomato (Lycopersicon peruvianum) and use this and other information to construct a model of this molecule. The S3-RNase has a single N-glycosylation site (Asn-28) to which one of three N-glycans is attached. S3-RNase has seven Cys residues; six are involved in disulfide linkages (Cys-16-Cys-21, Cys-46-Cys-91, and Cys-166-Cys-177), and one has a free thiol group (Cys-150). The disulfide-bonding pattern is consistent with that observed in RNase Rh, a related RNase for which radiographic-crystallographic information is available. A molecular model of the S3-RNase shows that four of the most variable regions of the S-RNases are clustered on one surface of the molecule. This is discussed in the context of recent experiments that set out to determine the regions of the S-RNase important for recognition during the self-incompatibility response.

Fine specificity of the antibody response to myelin basic protein in the central nervous system in multiple sclerosis: the minimal B-cell epitope and a model of its features.

T cells, B cells, and antibody are found in the white matter of the central nervous system in multiple sclerosis. The epitope center for the antibody response to human myelin basic protein (MBP) fits precisely the minimal epitope Pro85-Val-Val-His-Phe-Phe-Lys-Asn-Ile-Val-Thr-Pro96 for that reported for HLA DR2b (DRB1*1501)-restricted T cells that recognize MBP [Wucherpfenning, K.W., Sette, A., Southwood, S., Oseroff, C., Matsui, M., Strominger, J. & Hafler, D. (1994) J. Exp. Med. 179, 279-290], and overlaps with the reported DR2a-restricted epitope for T cells reactive to MBP [Martin, R., Howell, M. D., Jaraquemada, D., Furlage, M., Richert, J., Brostoff, S., Long, E. O., McFarlin, D. E. & McFarland, H. F. (1991) J. Exp. Med. 173, 19-24]. We describe a molecular model of this epitope.

Subunit-destabilizing mutations in Drosophila copper/zinc superoxide dismutase: neuropathology and a model of dimer dysequilibrium.

Mutations in Cu/Zn superoxide dismutase (SOD), a hallmark of familial amyotrophic lateral sclerosis (FALS) in humans, are shown here to confer striking neuropathology in Drosophila. Heterozygotes with one wild-type and one deleted SOD allele retain the expected 50% of normal activity for this dimeric enzyme. However, heterozygotes with one wild-type and one missense SOD allele show lesser SOD activities, ranging from 37% for a heterozygote carrying a missense mutation predicted from structural models to destabilize the dimer interface, to an average of 13% for several heterozygotes carrying missense mutations predicted to destabilize the subunit fold. Genetic and biochemical evidence suggests a model of dimer dysequilibrium whereby SOD activity in missense heterozygotes is reduced through entrapment of wild-type subunits into unstable or enzymatically inactive heterodimers. This dramatic impairment of the activity of wild-type subunits in vivo has implications for our understanding of FALS and for possible therapeutic strategies.

Both apolipoprotein E and A-I genes are present in a nonmammalian vertebrate and are highly expressed during embryonic development

Apolipoprotein E (apoE) is associated with several classes of plasma lipoproteins and mediates uptake of lipoproteins through its ability to interact with specific cell surface receptors. Besides its role in cardiovascular diseases, accumulating evidence has suggested that apoE could play a role in neurodegenerative diseases, such as Alzheimer disease. In vertebrates, apoA-I is the major protein of high-density lipoprotein. ApoA-I may play an important role in regulating the cholesterol content of peripheral tissues through the reverse cholesterol transport pathway. We have isolated cDNA clones that code for apoE and apoA-I from a zebrafish embryo library. Analysis of the deduced amino acid sequences showed the presence of a region enriched in basic amino acids in zebrafish apoE similar to the lipoprotein receptor-binding region of human apoE. We demonstrated by whole-mount in situ hybridization that apoE and apoA-I genes are highly expressed in the yolk syncytial layer, an extraembryonic structure implicated in embryonic and larval nutrition. ApoE transcripts were also observed in the deep cell layer during blastula stage, in numerous ectodermal derivatives after gastrulation, and after 3 days of development in a limited number of cells both in brain and in the eyes. Our data indicate that apoE can be found in a nonmammalian vertebrate and that the duplication events, from which apoE and apoA-I genes arose, occurred before the divergence of the tetrapod and teleost ancestors. Zebrafish can be used as a simple and useful model for studying the role of apolipoproteins in embryonic and larval nutrition and of apoE in brain morphogenesis and regeneration.

Cleavage of the antithrombin III binding site in heparin by heparinases and its implication in the generation of low molecular weight heparin

Heparin has been used as a clinical anticoagulant for more than 50 years, making it one of the most effective pharmacological agents known. Much of heparin's activity can be traced to its ability to bind antithrombin III (AT-III). Low molecular weight heparin (LMWH), derived from heparin by its controlled breakdown, maintains much of the antithrombotic activity of heparin without many of the serious side effects. The clinical significance of LMWH has highlighted the need to understand and develop chemical or enzymatic means to generate it. The primary enzymatic tools used for the production of LMWH are the heparinases from Flavobacterium heparinum, specifically heparinases I and II. Using pentasaccharide and hexasaccharide model compounds, we show that heparinases I and II, but not heparinase III, cleave the AT-III binding site, leaving only a partially intact site. Furthermore, we show herein that glucosamine 3-O sulfation at the reducing end of a glycosidic linkage imparts resistance to heparinase I, II, and III cleavage. Finally, we examine the biological and pharmacological consequences of a heparin oligosaccharide that contains only a partial AT-III binding site. We show that such an oligosaccharide lacks some of the functional attributes of heparin- and heparan sulfate-like glycosaminoglycans containing an intact AT-III site.

Loss of signaling through the G protein, Gz, results in abnormal platelet activation and altered responses to psychoactive drugs

Heterotrimeric G proteins mediate the earliest step in cell responses to external events by linking cell surface receptors to intracellular signaling pathways. Gz is a member of the Gi family of G proteins that is prominently expressed in platelets and brain. Here, we show that deletion of the α subunit of Gz in mice: (i) impairs platelet aggregation by preventing the inhibition of cAMP formation normally seen at physiologic concentrations of epinephrine, and (ii) causes the mice to be more resistant to fatal thromboembolism. Loss of Gzα also results in greatly exaggerated responses to cocaine, reduces the analgesic effects of morphine, and abolishes the effects of widely used antidepressant drugs that act as catecholamine reuptake inhibitors. These changes occur despite the presence of other Giα family members in the same cells and are not accompanied by detectable compensatory changes in the level of expression of other G protein subunits. Therefore, these results provide insights into receptor selectivity among G proteins and a model for understanding platelet function and the effects of psychoactive drugs.

Evolution by the birth-and-death process in multigene families of the vertebrate immune system

Concerted evolution is often invoked to explain the diversity and evolution of the multigene families of major histocompatibility complex (MHC) genes and immunoglobulin (Ig) genes. However, this hypothesis has been controversial because the member genes of these families from the same species are not necessarily more closely related to one another than to the genes from different species. To resolve this controversy, we conducted phylogenetic analyses of several multigene families of the MHC and Ig systems. The results show that the evolutionary pattern of these families is quite different from that of concerted evolution but is in agreement with the birth-and-death model of evolution in which new genes are created by repeated gene duplication and some duplicate genes are maintained in the genome for a long time but others are deleted or become nonfunctional by deleterious mutations. We found little evidence that interlocus gene conversion plays an important role in the evolution of MHC and Ig multigene families.

Binding and disruption of phospholipid bilayers by supramolecular RNA complexes

In an RNA world, RNAs would have regulated traffic through normally impermeable bilayer membranes. Using selection-amplification we previously found RNAs that bind stably and increase the ionic conductance of phospholipid membranes at high Mg2+ and Ca2+ concentrations. Now selection in reduced divalents yields RNAs that bind phosphatidylcholine liposomes under conditions closer to physiological. Such affinity for phospholipid membranes requires interactions between RNAs. In fact, we detected no functional monomeric membrane-binding RNAs. A membrane-active end-to-end heterotrimer consisting of 2 RNA 9 and 1 RNA 10 is defined by nucleotide protection, oligonucleotide competition, and mutant analysis. Oligomers of the heterotrimer bind stably, cause release of liposome-encapsulated solutes, and disrupt model black membranes. Individual RNA molecules do not show any of these activities. This novel mechanism of RNA binding to lipid membranes may not only regulate membrane permeability, but suggests that arrays of catalytic or structural RNAs on membranes are plausible. Finally, a selection met only by RNA complexes evokes new possibilities for selection-amplification itself.

Intracellular pH in adipocytes: effects of free fatty acid diffusion across the plasma membrane, lipolytic agonists, and insulin.

The main function of white adipose tissue is to store nutrient energy in the form of triglycerides. The mechanism by which free fatty acids (FFA) move into and out of the adipocyte has not been resolved. We show here that changes in intracellular pH (pH1) in adipocytes correlate with the movement of FFA across cellular membranes as predicted by the Kamp and Hamilton model of passive diffusion of FFA. Exposure of fat cells to lipolytic agents or external FFA results is a rapid intracellular acidification that is reversed by metabolism of the FFA or its removal by albumin. In contrast, insulin causes an alkalinization of the cell, consistent with its main function to promote esterification. Inhibition of Na+/H+ exchange in adipocytes does not prevent the changes in pHi caused by FFA, lipolytic agents, or insulin. A fatty acid dimer, which diffuses into the cell but is not metabolized, causes an irreversible acidification. Taken together, the data suggest that changes in pHi occur in adipocytes in response to the passive diffusion of un-ionized FFA (flip-flop) into and out of the cell and in response to their metabolism and production within the cell. These changes in pHi may, in turn, modulate hormonal signaling and metabolism with significant impact on cell function.

Shared fine specificity between T-cell receptors and an antibody recognizing a peptide/major histocompatibility class I complex.

Cytotoxic T cells recognize mosaic structures consisting of target peptides embedded within self-major histocompatibility complex (MHC) class I molecules. This structure has been described in great detail for several peptide-MHC complexes. In contrast, how T-cell receptors recognize peptide-MHC complexes have been less well characterized. We have used a complete set of singly substituted analogs of a mouse MHC class I, Kk-restricted peptide, influenza hemagglutinin (Ha)255-262, to address the binding specificity of this MHC molecule. Using the same peptide-MHC complexes we determined the fine specificity of two Ha255-262-specific, Kk-restricted T cells, and of a unique antibody, pSAN, specific for the same peptide-MHC complex. Independently, a model of the Ha255-262-Kk complex was generated through homology modeling and molecular mechanics refinement. The functional data and the model corroborated each other showing that peptide residues 1, 3, 4, 6, and 7 were exposed on the MHC surface and recognized by the T cells. Thus, the majority, and perhaps all, of the side chains of the non-primary anchor residues may be available for T-cell recognition, and contribute to the stringent specificity of T cells. A striking similarity between the specificity of the T cells and that of the pSAN antibody was found and most of the peptide residues, which could be recognized by the T cells, could also be recognized by the antibody.

Evolutionarily conserved Galphabetagamma binding surfaces support a model of the G protein-receptor complex.

The pivotal role of G proteins in sensory, hormonal, inflammatory, and proliferative responses has provoked intense interest in understanding how they interact with their receptors and effectors. Nonetheless, the locations of the receptors and effector binding sites remain poorly characterized, although nearly complete structures of the alphabetagamma heterotrimeric complex are available. Here we apply evolutionary trace (ET) analysis [Lichtarge, O., Bourne, H. R. & Cohen, F. E. (1996) J. Mol. Biol. 257, 342-358] to propose plausible locations for these sites. On each subunit, ET identifies evolutionarily selected surfaces composed of residues that do not vary within functional subgroups and that form spatial clusters. Four clusters correctly identify subunit interfaces, and additional clusters on Galpha point to likely receptor or effector binding sites. Our results implicate the conformationally variable region of Galpha in an effector binding role. Furthermore the range of predicted interactions between the receptor and Galphabetagamma, is sufficiently limited that we can build a low resolution and testable model of the receptor-G protein complex.

The kinesin walk: a dynamic model with elastically coupled heads.

Recently individual two-headed kinesin molecules have been studied in in vitro motility assays revealing a number of their peculiar transport properties. In this paper we propose a simple and robust model for the kinesin stepping process with elastically coupled Brownian heads that show all of these properties. The analytic and numerical treatment of our model results in a very good fit to the experimental data and practically has no free parameters. Changing the values of the parameters in the restricted range allowed by the related experimental estimates has almost no effect on the shape of the curves and results mainly in a variation of the zero load velocity that can be directly fitted to the measured data. In addition, the model is consistent with the measured pathway of the kinesin ATPase.

Dynamic friction and the origin of the complexity of earthquake sources.

We study a simple antiplane fault of finite length embedded in a homogeneous isotropic elastic solid to understand the origin of seismic source heterogeneity in the presence of nonlinear rate- and state-dependent friction. All the mechanical properties of the medium and friction are assumed homogeneous. Friction includes a characteristic length that is longer than the grid size so that our models have a well-defined continuum limit. Starting from a heterogeneous initial stress distribution, we apply a slowly increasing uniform stress load far from the fault and we simulate the seismicity for a few 1000 events. The style of seismicity produced by this model is determined by a control parameter associated with the degree of rate dependence of friction. For classical friction models with rate-independent friction, no complexity appears and seismicity is perfectly periodic. For weakly rate-dependent friction, large ruptures are still periodic, but small seismicity becomes increasingly nonstationary. When friction is highly rate-dependent, seismicity becomes nonperiodic and ruptures of all sizes occur inside the fault. Highly rate-dependent friction destabilizes the healing process producing premature healing of slip and partial stress drop. Partial stress drop produces large variations in the state of stress that in turn produce earthquakes of different sizes. Similar results have been found by other authors using the Burridge and Knopoff model. We conjecture that all models in which static stress drop is only a fraction of the dynamic stress drop produce stress heterogeneity.