Gapped BLAST and PSI-BLAST: a new generation of protein database search programs

Identification of a Frizzled-like cysteine rich domain in the extracellular region of developmental receptor tyrosine kinases

In Drosophila, members of the Frizzled family of tissue-polarity genes encode proteins that appear to function as cell-surface receptors for Wnts. The Frizzled genes belong to the seven transmembrane class of receptors (7TMR) and have on their extracellular region a cysteine-rich domain that has been implicated as the Wnt binding domain. This region has a characteristic spacing of ten cysteines, which has also been identified in FrzB (a secreted antagonist of Wnt signaling) and Smoothened (another 7TMR, which is involved in suppression of the hedgehog pathway). We have identified, using BLAST, sequence similarity between the cysteine-rich domain of Frizzled and several receptor tyrosine kinases, which have roles in development. These include the muscle-specific receptor tyrosine kinase (MuSK), the neuronal specific kinase (NSK2), and ROR1 and ROR2. At present, the ligands for these developmental tyrosine kinases are unknown. Our results suggest that Wnt-like ligands may bind to these developmental tyrosine kinases.

Cloning of human striatin cDNA (STRN), gene mapping to 2p22-p21, and preferential expression in brain

Rat striatin, a recently discovered calmodulin-binding protein belonging to the WD repeat family, is expressed in neurons, mostly in the striatum and motor and olfactory systems. Striatin is localized in the somato-dendritic compartment of neurons, mainly in the spines. It may play a role in dendritic Ca2+ signaling. Here we report the cloning and sequencing of human striatin cDNA (HGMW-approved symbol STRN), the localization of the gene to chromosome 2p22-p21, and its preferential expression in brain. The human cDNA sequence is predicted to encode a 780-amino-acid protein possessing eight WD repeats. Striatin is highly conserved between rat and human with 96% identity and 98% similarity at the amino acid level. Since the Caenorabditis elegans genome also contains a closely related striatin coding sequence, the function of striatin is likely to be well conserved.

Detecting distant relatives of mammalian LPS-binding and lipid transport proteins

In mammals, a family of four lipid binding proteins has been previously defined that includes two lipopolysaccharide binding proteins and two lipid transfer proteins. The first member of this family to have its three-dimensional structure determined is bactericidal/permeability-increasing protein (BPI). Using both the sequence and structure of BPI, along with recently developed sequence-sequence and sequence-structure similarity search methods, we have identified 13 distant members of the family in a diverse set of eukaryotes, including rat, chicken, Caenorhabditis elegans, and Biomphalaria galbrata. Although the sequence similarity between these 13 new members and any of the 4 original members of the BPI family is well below the "twilight zone," their high sequence-structure compatibility with BPI indicates they are likely to share its fold. These findings broaden the BPI family to include a member found in retina and brain, and suggest that a primitive member may have contained only one of the two similar domains of BPI.

Forum domain in Drosophila melanogaster cut locus possesses looped domains inside

We have studied the relationship between chromosomal forum domains and looped domains in the cut locus of Drosophila melanogaster . Forum domains were earlier detected by separation in pulsed-field gels of 50-150 kb chromosomal DNA fragments obtained after spontaneous non-random degradation of chromosomes. We have localized the boundary region where cleavage sites are scattered between two forum domains in the regulatory region of the cut locus. We have sequenced a 13 kb region spanning few kilobases from distal domain, the boundary region and part of the proximal forum domain where several scaffold associated regions (SARs) were observed. We conclude that forum domains and looped domains are physically different types of domains and belong to different levels of organization in eukaryotic chromosomes. The boundary region between the neighboring forum domains in the cut locus possesses the Doc element insertion and a micro-satellite stretch and thus might remind a small island of heterochromatin and correspond to so-called intercalary heterochromatin that is known to be located in the 7B1-2 band where the major part of the cut locus is reside.

Sirohaem sulfite reductase and other proteins encoded by genes at the dsr locus of Chromatium vinosum are involved in the oxidation of intracellular sulfur

The sequence of the dsr gene region of the phototrophic sulfur bacterium Chromatium vinosum D (DSMZ 180) was determined to clarify the in vivo role of 'reverse' sirohaem sulfite reductase. The dsrAB genes encoding dissimilatory sulfite reductase are part of a gene cluster, dsrABEFHCMK, that encodes four small, soluble proteins (DsrE, DsrF, DsrH and DsrC), a transmembrane protein (DsrM) with similarity to haem-b-binding polypeptides and a soluble protein (DsrK) resembling [4Fe-4S]-cluster-containing heterodisulfide reductase from methanogenic archaea. Northern hybridizations showed that expression of the dsr genes is increased by the presence of reduced sulfur compounds. The dsr genes are not only transcribed from a putative promoter upstream of dsrA but primary transcripts originating from (a) transcription start site(s) downstream of dsrB are also formed. Polar insertion mutations immediately upstream of dsrA, and in dsrB, dsrH and dsrM, led to an inability of the cells to oxidize intracellularly stored sulfur. The capability of the mutants to oxidize sulfide, thiosulfate and sulfite under photolithoautotrophic conditions was unaltered. Photoorganoheterotrophic growth was also unaffected. 'Reverse' sulfite reductase and DsrEFHCMK are, therefore, not essential for oxidation of sulfide or thiosulfate, but are obligatory for sulfur oxidation. These results, together with the finding that the sulfur globules of C. vinosum are located in the extracytoplasmic space whilst the dsr gene products appear to be either cytoplasmic or membrane-bound led to the proposal of new models for the pathway of sulfur oxidation in this phototrophic sulfur bacterium.

Purification and characterization of a thermostable class II fumarase from Thermus thermophilus

A thermostable fumarase was purified from a strain of Thermus thermophilus isolated from a Japanese hot spring. The maximum specific activity of the purified enzyme was 1740 units/mg at pH 8.0 and 85 degreesC. The enzyme was composed of four identical subunits with a molecular weight of 46,000 and displayed other enzymatic characteristics which are common to the class II fumarases. The thermal stability of the purified enzyme was remarkable, with over 80% of the activity remaining after a 24-h incubation at 90 degreesC. The enzyme was also resistant to chemical denaturants; 50% of the initial specific activity was detected in assay mixtures containing 0.8 M guanidine hydrochloride. The purified enzyme shared an extremely high sequence homology with Thermus aquaticus fumarase and Bacillus subtilis fumarase in the first 43 amino acid residues.

Secretion of elastinolytic enzymes and their propeptides by Pseudomonas aeruginosa

Elastase of Pseudomonas aeruginosa is synthesized as a preproenzyme. The signal sequence is cleaved off during transport across the inner membrane and, in the periplasm, proelastase is further processed. We demonstrate that the propeptide and the mature elastase are both secreted but that the propeptide is degraded extracellularly. In addition, reduction of the extracellular proteolytic activity led to the accumulation of unprocessed forms of LasA and LasD in the extracellular medium, which shows that these enzymes are secreted in association with their propeptides. Furthermore, a hitherto undefined protein with homology to a Streptomyces griseus aminopeptidase accumulated under these conditions.

Eucalypt NADP-dependent isocitrate dehydrogenase. cDNA cloning and expression in ectomycorrhizae

NADP-dependent isocitrate dehydrogenase (NADP-ICDH) activity is increased in roots of Eucalyptus globulus subsp. bicostata ex Maiden Kirkp. during colonization by the ectomycorrhizal fungus Pisolithus tinctorius Coker and Couch. To investigate the regulation of the enzyme expression, a cDNA (EgIcdh) encoding the NADP-ICDH was isolated from a cDNA library of E. globulus-P. tinctorius ectomycorrhizae. The putative polypeptide sequence of EgIcdh showed a high amino acid similarity with plant NADP-ICDHs. Because the deduced EgICDH protein lacks an amino-terminal targeting sequence and shows highest similarity to plant cytosolic ICDHs, it probably represents a cytoplasmic isoform. RNA analysis showed that the steady-state level of EgIcdh transcripts was enhanced nearly 2-fold in ectomycorrhizal roots compared with nonmycorrhizal roots. Increased accumulation of NADP-ICDH transcripts occurred as early as 2 d after contact and likely led to the observed increased enzyme activity. Indirect immunofluorescence microscopy indicated that NADP-ICDH was preferentially accumulated in the epidermis and stele parenchyma of nonmycorrhizal and ectomycorrhizal lateral roots. The putative role of cytosolic NADP-ICDH in ectomycorrhizae is discussed.

Membrane-bound lytic endotransglycosylase in Escherichia coli

The gene for a novel endotype membrane-bound lytic transglycosylase, emtA, was mapped at 26.7 min of the E. coli chromosome. EmtA is a lipoprotein with an apparent molecular mass of 22kDa. Overexpression of the emtA gene did not result in bacteriolysis in vivo, but the enzyme was shown to hydrolyze glycan strands isolated from murein by amidase treatment. The formation of tetra- and hexasaccharides, but no disaccharides, reflects the endospecificity of the enzyme. The products are characterized by the presence of 1,6-anhydromuramic acid, indicating a lytic transglycosylase reaction mechanism. EmtA may function as a formatting enzyme that trims the nascent murein strands produced by the murein synthesis machinery into proper sizes, or it may be involved in the formation of tightly controlled minor holes in the murein sacculus to facilitate the export of bulky compounds across the murein barrier.

Biosynthesis of 2-aceto-2-hydroxy acids: acetolactate synthases and acetohydroxyacid synthases

Two groups of enzymes are classified as acetolactate synthase (EC 4. 1.3.18). This review deals chiefly with the FAD-dependent, biosynthetic enzymes which readily catalyze the formation of acetohydroxybutyrate from pyruvate and 2-oxobutyrate, as well as of acetolactate from two molecules of pyruvate (the ALS/AHAS group). These enzymes are generally susceptible to inhibition by one or more of the branched-chain amino acids which are ultimate products of the acetohydroxyacids, as well as by several classes of herbicides (sulfonylureas, imidazolinones and others). Some ALS/AHASs also catalyze the (non-physiological) oxidative decarboxylation of pyruvate, leading to peracetic acid; the possible relationship of this process to oxygen toxicity is considered. The bacterial ALS/AHAS which have been well characterized consist of catalytic subunits (around 60 kDa) and smaller regulatory subunits in an alpha2beta2 structure. In the case of Escherichia coli isozyme III, assembly and dissociation of the holoenzyme has been studied. The quaternary structure of the eukaryotic enzymes is less clear and in plants and yeast only catalytic polypeptides (homologous to those of bacteria) have been clearly identified. The presence of regulatory polypeptides in these organisms cannot be ruled out, however, and genes which encode putative ALS/AHAS regulatory subunits have been identified in some cases. A consensus sequence can be constructed from the 21 sequences which have been shown experimentally to represent ALS/AHAS catalytic polypeptides. Many other sequences fit this consensus, but some genes identified as putative 'acetolactate synthase genes' are almost certainly not ALS/AHAS. The solution of the crystal structures of several thiamin diphosphate (ThDP)-dependent enzymes which are homologous to ALS/AHAS, together with the availability of many amino acid sequences for the latter enzymes, has made it possible for two laboratories to propose similar, reasonable models for a dimer of catalytic subunits of an ALS/AHAS. A number of characteristics of these enzymes can now be better understood on the basis of such models: the nature of the herbicide binding site, the structural role of FAD and the binding of ThDP-Mg2+. The models are also guides for experimental testing of ideas concerning structure-function relationships in these enzymes, e.g. the nature of the substrate recognition site. Among the important remaining questions is how the enzyme suppresses alternative reactions of the intrinsically reactive hydroxyethylThDP enamine formed by the decarboxylation of the first substrate molecule and specifically promotes its condensation with 2-oxobutyrate or pyruvate.

The PR domain of the Rb-binding zinc finger protein RIZ1 is a protein binding interface and is related to the SET domain functioning in chromatin-mediated gene expression

The PR domain, first noted as the PRDI-BF1-RIZ1 homologous region, defines a sub-class of zinc finger genes that appear to function as negative regulators of tumorigenesis. This family includes the MDS1-EVI1 gene inactivated in myeloid leukemia, the PRDI-BF1/BLIMP1 transcription repressor of c-myc involved in driving B-cell differentiation, and the RIZ gene, which encodes proteins capable of binding to the retinoblastoma tumor suppressor protein (Rb). The PR domain of MDS1-EVI1 is disrupted by translocations linked to myeloid leukemia, resulting in the activation of the PR-minus oncogenic product EVI1. Remarkably similar to MDS1-EVI1, RIZ gene also normally produces two protein products of different length, and the smaller protein RIZ2 lacks the PR domain of RIZ1 but is otherwise identical to RIZ1. These observations raise considerable interest to determine the function of PR. We show here that RIZ1 PR domain mediates protein-protein interaction. Recombinant fusion proteins of PR can bind to in vitro translated RIZ1 and RIZ2 proteins. The binding can be disrupted by amino acid substitutions at conserved residues of PR, suggesting that binding is specific. Of the three conserved exons of PR, the first two appear dispensable for binding, whereas the third exon is required. A region in the carboxyl terminus of RIZ proteins was mapped to be necessary and sufficient for PR binding. We also found that the PR domain shares significant sequence identity to the SET domain present in chromosomal proteins that function in modulating gene expression from yeast to mammals. Our data suggest that the PR domain is a derivative of SET domain and may function as protein binding interface in the regulation of chromatin-mediated gene expression.

Activation of extracellular signal-regulated kinase 2 by a novel Abl-binding protein, ST5

The human ST5 gene encodes three proteins with predicted molecular masses of 126, 82, and 70 kDa. These widely expressed proteins share a C-terminal region that bears significant sequence homology to a group of GDP/GTP exchange proteins for the Rab3 family of small GTP binding proteins. The N-terminal region of the largest ST5 protein, p126, contains two proline-rich sequences, PR1 and PR2, with consensus motifs similar to Src homology 3 (SH3) binding regions and to mitogen-activated protein kinase (MAPK) phosphorylation sites. Based on these properties, we sought to investigate the activity of ST5 proteins in signal transduction pathways. In vitro, p126 displayed preferential binding to c-Abl SH3, as compared with other SH3 domains. This interaction was mediated by the PR2 sequence. In vivo, expression of p126, but not p82 or p70, activated MAPK/ERK2 in response to EGF in COS-7 cells. Expression of c-Abl with p126 greatly enhanced this activity. Deletion of PR1 blocked the ability of p126 to activate ERK2. Deletion of PR2 did not affect the basal activity, but blocked the stimulatory effect of c-Abl. Whereas p82 expression had no effect on ERK2 activation by p126, p70 completely abrogated this activity. These observations suggest that ST5 can function as a signaling protein and can provide a link between c-Abl and ERK2.

The complete genome of Bacillus subtilis: from sequence annotation to data management and analysis

The completion of the entire 4.2-Mb genome sequence of the gram-positive bacterium Bacillus subtilis has been a milestone for biological studies on this model organism. This paper describes bioinformatics work related to this joint European and Japanese project: methods and strategies for gene annotation and detection of sequencing errors, using an integrated cooperative computer environment (Imagene); construction of a specialized database for data management and a WWW server for data retrieval (SubtiList); DNA sequence analysis, yielding striking results on oligonucleotide bias, repeated sequences, and codon usage, all landmarks of evolutionary events shaping the B. subtilis genome.

Derivation of the three-dimensional architecture of bacterial ribonuclease P RNAs from comparative sequence analysis

The secondary structure of bacterial RNase P RNA, a ribozyme responsible for the maturation of the 5' end of tRNAs, is well established on the basis of sequence comparison analysis. RNase P RNA secondary structures fall into two types, A and B, which share a common core formed by the assembly of two main folding domains, but differ in their peripheral elements.A revised alignment of 137 available sequences reveals new covariations allowing for the refinement of both types of secondary structures. Phylogenetic evidence is thus provided for the extension of stems P11, P14, P19, P10.1 and P15.1 through further canonical base-pairs or GAellipsisGA mismatches. These refinements led in turn to a new organization of the catalytic core, with coaxial stackings of helices P2 and P19 as well as P1 and P4. New inter-domain tertiary interactions involve loop L9 and helix P1 and loop L8 with helix P4. These features were incorporated into atomic-scale 3D models of RNase P RNA for representatives of each structural type, namely Escherichia coli and Bacillus subtilis. In each model, the juxtaposition of the core helices creates a cradle onto which the pre-tRNA substrate binds with most evolutionarily conserved residues converging towards the cleavage site. The inner cores of both types are stabilized similarly, albeit by different peripheral elements, emphasizing the modular and hierarchical organisation of the architecture of RNase P RNAs. Similarities are thus apparent between the type A modules, P16/P17/P6 and P13/P14, and their type B analogs, P5.1/P15.1 and P10. 1/P10.1a, respectively. Other noteworthy features of these models include compactness and good agreement with published crosslinking data.

Molecular cloning and expression of human cGMP-binding cGMP-specific phosphodiesterase (PDE5)

A human PDE5 cDNA has been isolated which contains an open reading frame encoding an 875 amino acid, 100,012 Da polypeptide, the expression of which yields a protein of the predicted size and is capable of hydrolyzing cGMP. The deduced amino acid sequence is very similar (95%) to that of bovine PDE5, and comprises a conserved cGMP-binding domain and catalytic domain. Northern analysis reveals a major and minor transcript of approximately 9 kb and approximately 8 kb respectively, thus indicating the existence of at least two splice variants, the major form being readily detected in bladder, colon, lung, pancreas, placenta, prostate, small intestine, and stomach.

Direct cloning and analysis of DNA sequences from a region of the Chinese hamster genome associated with aphidicolin-sensitive fragility

Fragile sites are reproducibly expressed and chemically induced decondensations on mitotic chromosomes observed under cytological conditions. They are classified both on the basis of the frequency with which they occur (rare and common) and in terms of the chemical agent used to induce expression in tissue culture cells. Aphidicolin-sensitive common fragile sites appear to be ubiquitous in humans and other mammals and have been considered as candidates of pathological importance. Recently DNA from FRA3B, the most highly expressed constitutive fragile site in the human genome, has been cloned although as yet the cause of the underlying fragility has not been identified. In this study we describe the isolation, using a direct cloning approach, of DNA from a region of the Chinese hamster genome associated with aphidicolin-inducible fragility. Cells of a human-hamster somatic cell hybrid were transfected with a pSV2HPRT vector while exposed to aphidicolin, an inhibitor of DNA polymerases alpha, delta and epsilon. FISH analysis of stable transfectant clones revealed that the ingoing plasmid DNA had preferentially integrated into fragile site-containing chromosomal bands. Plasmid rescue was used to recover DNA sequences flanking one such integration site in the hamster genome. We demonstrate by FISH analysis of metaphase cells induced with aphidicolin that the rescued DNA is from a region of fragility on Chinese hamster chromosome 2, distal to the DHFR locus. Analysis of the DNA sequences flanking the integration site revealed the overall A+T content of the 3,725 bp region sequenced to be 63.3%, with a highly [A].[T]-rich 156 bp region (86.5%) almost adjacent to the integration site. Computational analyses have identified strong homologies to Saccharomyces cerevisiae autonomous replicating sequences (ARS), polypyrimidine tracts, scaffold attachment site consensus sequences and a 24 bp consensus sequence highly conserved in eukaryotic replication origins, all of which appear to cluster around the [A].[T]-rich sequences. This domain also possesses structural characteristics which are common to both prokaryotic and eukaryotic origins of replications, in particular an unusually straight conformation of low thermal stability flanked either side by highly bent DNA segments. Further isolation and characterisation of DNA sequences from common fragile sites will facilitate studies into the underlying nature of these enigmatic regions of the mammalian genome, leading to a greater understanding of chromatin structure.

Rpp2, an essential protein subunit of nuclear RNase P, is required for processing of precursor tRNAs and 35S precursor rRNA in Saccharomyces cerevisiae

RPP2, an essential gene that encodes a 15.8-kDa protein subunit of nuclear RNase P, has been identified in the genome of Saccharomyces cerevisiae. Rpp2 was detected by sequence similarity with a human protein, Rpp20, which copurifies with human RNase P. Epitope-tagged Rpp2 can be found in association with both RNase P and RNase mitochondrial RNA processing in immunoprecipitates from crude extracts of cells. Depletion of Rpp2 protein in vivo causes accumulation of precursor tRNAs with unprocessed introns and 5' and 3' termini, and leads to defects in the processing of the 35S precursor rRNA. Rpp2-depleted cells are defective in processing of the 5.8S rRNA. Rpp2 immunoprecipitates cleave both yeast precursor tRNAs and precursor rRNAs accurately at the expected sites and contain the Rpp1 protein orthologue of the human scleroderma autoimmune antigen, Rpp30. These results demonstrate that Rpp2 is a protein subunit of nuclear RNase P that is functionally conserved in eukaryotes from yeast to humans.

Kir2.4: a novel K+ inward rectifier channel associated with motoneurons of cranial nerve nuclei

Members of the Kir2 subfamily of inwardly rectifying K+ channels characterized by their strong current rectification are widely expressed both in the periphery and in the CNS in mammals. We have cloned from rat brain a fourth subfamily member, designated Kir2.4 (IRK4), which shares 53-63% similarity to Kir2.1, Kir2.2, or Kir2.3 on the amino acid level. In situ hybridization analysis identifies Kir2.4 as the most restricted of all Kir subunits in the brain. Kir2. 4 transcripts are expressed predominantly in motoneurons of cranial nerve motor nuclei within the general somatic and special visceral motor cell column and thus are uniquely related to a functional system. Heterologous expression of Kir2.4 in Xenopus oocytes and mammalian cells gives rise to low-conductance channels (15 pS), with an affinity to the channel blockers Ba2+ (Ki = 390 microM) and Cs+ (Ki = 8.06 mM) 30-50-fold lower than in other Kir channels. Low Ba2+ sensitivity allows dissection of Kir2.4 currents from other Kir conductances in hypoglossal motoneurons (HMs) in rat brainstem slices. The finding that Ba2+-mediated block of Kir2.4 in HMs evokes tonic activity and increases the frequency of induced spike discharge indicates that Kir2.4 channels are of major importance in controlling excitability of motoneurons in situ.

Positional cloning of the gene for Nijmegen breakage syndrome

Nijmegen breakage syndrome (NBS), also known as ataxia-telangiectasia (AT) variant, is an autosomal recessive disorder characterized by microcephaly, growth retardation, severe combined immunodeficiency and a high incidence of lymphoid cancers. Cells from NBS patients display chromosome instability, hypersensitivity to ionizing radiation and abnormal cell-cycle regulation after irradiation, all of which are characteristics shared with AT. Recently, the NBS locus was mapped at 8q21 by two independent approaches, complementation studies and linkage analysis. Here, we report the positional cloning of the NBS gene, NBS1, from an 800-kb candidate region. The gene comprises 50 kb and encodes a protein of 754 amino acids. The amino-terminal region of the protein shows weak homology to the yeast XRS2, MEK1, CDS1 and SPK1 proteins. The gene is expressed at high levels in the testes, suggesting that it might be involved in meiotic recombination. We detected the same 5-bp deletion in 13 individuals, and conclude that it is likely to be a founder mutation.

OpaR, a homolog of Vibrio harveyi LuxR, controls opacity of Vibrio parahaemolyticus

Vibrio parahaemolyticus is an organism well adapted to communal life on surfaces. When grown on a surface or in a viscous layer, the bacterium induces a large gene system and differentiates to swarmer cells capable of movement over and colonization of surfaces. V. parahaemolyticus displays additional phenotypic versatility manifested as variable colony morphology, switching between translucent and opaque colony types. Although not itself luminescent, V. parahaemolyticus produces autoinducer molecules capable of inducing luminescence in Vibrio harveyi. To examine the role of quorum signaling in the lifestyles of V. parahaemolyticus, the functional homolog of the gene encoding the V. harveyi autoinducer-controlled transcriptional regulatory protein LuxR was cloned. Sequence analysis of the clone predicted an open reading frame with a deduced product 96% identical to LuxR. Introduction of the clone carrying the luxR-like locus into V. parahaemolyticus dramatically affected colony morphology, converting a translucent strain to an opaque one. When the coding sequence for the luxR homolog was placed under the control of the Ptac promoter, conversion to the opaque phenotype became inducible by isopropyl-beta-D-thiogalactopyranoside. Allelic disruption of the luxR-like gene on the chromosome of an opaque strain produced a translucent strain proficient in swarming ability. Primer extension mapping demonstrated opaR transcription in opaque but not translucent cell types. It is postulated that this gene, which has been named opaR, encodes a transcription factor controlling cell type. The underlying genetic basis for opaque-translucent variation may be the consequence of a genomic alteration detected in the opaR locus of opaque and translucent strains.

How far divergent evolution goes in proteins

In theory, mutations of protein sequences may eventually generate different functions as well as different structures. The observation of such records of protein evolution have been obscured by the dissipation of memory about the ancestors. In the past year, new advances in our understanding of divergent evolution were allowed by new protein structure determinations, including the ClpP proteases, steroid delta-isomerase, carboxypeptidase G2, the thrombin inhibitor triabin and the chloroplast Rieske protein. There is strong evidence for their distant homology with proteins of known structure despite significant functional or structural differences.

Beyond complete genomes: from sequence to structure and function

Computer analysis of complete prokaryotic genomes shows that microbial proteins are in general highly conserved--approximately 70% of them contain ancient conserved regions. This allows us to delineate families of orthologs across a wide phylogenetic range and, in many cases, predict protein functions with considerable precision. Sequence database searches using newly developed, sensitive algorithms result in the unification of such orthologous families into larger superfamilies sharing common sequence motifs. For many of these superfamilies, prediction of the structural fold and specific amino acid residues involved in enzymatic catalysis is possible. Taken together, sequence and structure comparisons provide a powerful methodology that can successfully complement traditional experimental approaches.

Protein fold irregularities that hinder sequence analysis

The detection of homologous protein sequences frequently provides useful predictions of function and structure. Methods for homology searching have continued to improve, such that very distant evolutionary relationships can now be detected. Little attention has been paid, however, to the problems of detecting homology when domains are inserted or permuted. Here we review recent occurrences of these phenomena and discuss methods that permit their detection.

Molecular characterization of Treponema pallidum mcp2, a putative chemotaxis protein gene

The nucleotide sequence of the Treponema pallidum mcp2 gene was determined. mcp2 encodes a 45.8-kDa protein whose deduced amino acid sequence has significant homology with the C-terminal region of bacterial methyl-accepting chemotaxis proteins (MCPs). The Mcp2 N terminus lacks the hydrophobic transmembrane regions present in most MCPs. An Mcp2 fusion protein was strongly reactive with antibody (HC23) to the highly conserved domain of MCPs and with rabbit syphilitic serum. Antibody HC23 reacted with six T. pallidum proteins, including a 45-kDa protein that may correspond to Mcp2. This protein was present in the aqueous phase from T. pallidum cells that were solubilized with Triton X-114 and phase partitioned.

Kir2.4: a novel K+ inward rectifier channel associated with motoneurons of cranial nerve nuclei

Members of the Kir2 subfamily of inwardly rectifying K+ channels characterized by their strong current rectification are widely expressed both in the periphery and in the CNS in mammals. We have cloned from rat brain a fourth subfamily member, designated Kir2.4 (IRK4), which shares 53-63% similarity to Kir2.1, Kir2.2, or Kir2.3 on the amino acid level. In situ hybridization analysis identifies Kir2.4 as the most restricted of all Kir subunits in the brain. Kir2. 4 transcripts are expressed predominantly in motoneurons of cranial nerve motor nuclei within the general somatic and special visceral motor cell column and thus are uniquely related to a functional system. Heterologous expression of Kir2.4 in Xenopus oocytes and mammalian cells gives rise to low-conductance channels (15 pS), with an affinity to the channel blockers Ba2+ (Ki = 390 microM) and Cs+ (Ki = 8.06 mM) 30-50-fold lower than in other Kir channels. Low Ba2+ sensitivity allows dissection of Kir2.4 currents from other Kir conductances in hypoglossal motoneurons (HMs) in rat brainstem slices. The finding that Ba2+-mediated block of Kir2.4 in HMs evokes tonic activity and increases the frequency of induced spike discharge indicates that Kir2.4 channels are of major importance in controlling excitability of motoneurons in situ.

Unification of protein families

Computational biology exploits the evolutionary connectivity between proteins and protein families to predict structural and functional properties of uncharacterized gene products. In the past year, conceptual and statistical refinements have substantially improved algorithms for the detection of remote homologues. In conjunction with the rapid growth of biological databases, the global organization of proteins into sequence families, functional families and structural families has become both pertinent and feasible.

Structural organization and sequence of CLN2, the defective gene in classical late infantile neuronal ceroid lipofuscinosis

Mutations in the CLN2 gene result in classical late infantile neuronal ceroid lipofuscinosis (LINCL), a fatal childhood neurodegenerative disease. In this report, we present the complete sequence of the human CLN2 gene and define its physical relationship with two other genes that have been previously mapped to chromosome 11p15. The CLN2 gene consists of 13 exons and 12 introns and spans 6.65 kb. By S1 mapping and primer extension, the 5'-terminus of the CLN2 mRNA was mapped to 32 nucleotides upstream of the proposed initiation codon. A number of other elements were found to be located in close proximity to CLN2, including the gene encoding transcription factor TAFII30, the gene encoding intregrin-linked kinase, and an approximately 914-bp fragment that is 82% identical to antithrombin III. In addition, an EST cDNA clone that is transcribed on the strand opposite to CLN2 and that overlaps a portion of the CLN2 gene was identified. Finally, a set of primer pairs are presented for the amplification of the coding sequences, putative promoter, and splice junctions of the CLN2 gene. Taken together, this information will facilitate the molecular analysis of and genetic testing for classical LINCL.

Current concepts of PLP and its role in the nervous system

Proteolipid protein (PLP) and its smaller isoform DM20 constitute the major myelin proteins of the CNS. Mutations of the X-linked Plp gene cause the heterogeneous syndromes of Pelizaeus-Merzbacher disease (PMD) and spastic paraplegia (SPG) in man and similar dysmyelinating disorders in a range of animal species. A variety of mutations including missense mutations, deletions, and duplications are responsible. Missense mutations cause a predicted alteration in primary structure of the encoded protein(s) and are generally associated with early onset of signs and generalised dysmyelination. The severity of the phenotype varies according to the particular codon involved and the influence of uncharacterised modifying genes. There is some evidence that the dysmyelination results from the altered protein acquiring a novel function deleterious to the oligodendrocyte's function. Transgenic mice carrying extra copies of the Plp gene provide a valid model of PMD/SPG due to gene duplication. Depending on the gene dosage, the phenotype can vary from early onset of severe and lethal dysmyelination through to a very late onset of a tract-specific demyelination and axonal degeneration. Mice with a null mutation of the Plp gene assemble and maintain normal amounts of myelin but develop a progressive axonopathy, again demonstrating tract specificity. The results indicate that the functions of PLP are far from clear. There is good evidence that it is involved in the formation of the intraperiod line of myelin, and the results from the knockout and transgenic mice suggest a role in the interaction of oligodendrocyte and axon.

Genomics and computational molecular biology

There has been a dramatic increase in the number of completely sequenced bacterial genomes during the past two years as a result of the efforts both of public genome agencies and the pharmaceutical industry. The availability of completely sequenced genomes permits more systematic analyses of genes, evolution and genome function than was otherwise possible. Using computational methods - which are used to identify genes and their functions including statistics, sequence similarity, motifs, profiles, protein folds and probabilistic models - it is possible to develop characteristic genome signatures, assign functions to genes, identify pathogenic genes, identify metabolic pathways, develop diagnostic probes and discover potential drug-binding sites. All of these directions are critical to understanding bacterial growth, pathogenicity and host-pathogen interactions.

Finding the genes in genomic DNA

Genome sequencing efforts will soon generate hundreds of millions of bases of human genomic DNA containing thousands of novel genes. In the past year, the accuracy of computational gene-finding methods has improved significantly, to the point where a reasonable approximation of the gene structures within an extended genomic region can often be predicted in advance of more detailed experimental studies.

Cloning genes encoding receptors related to chemoattractant receptors

We report the cloning of a novel human gene (GPR32) encoding a putative G-protein-coupled receptor (GPCR) of 356 amino acids and a related pseudogene psi GPR32. The deduced amino acid sequence of GPR32 shares 35-39% identity with members of the chemoattractant receptor family. psi GPR32 shares 93% nucleotide identity with GPR32. We identified a mouse EST encoding a putative GPCR (GPR33) of 309 amino acids. The deduced amino acid sequence of GPR33 shares 30-35% identity with members of the chemoattractant receptor family and 36% identity with the receptor encoded by GPR32. The human orthologue of GPR33 contains a single basepair substitution with respect to the mouse, resulting in the presence of an in-frame stop codon within the predicted second intracellular loop, demonstrating that it is a pseudogene. Through fluorescence in situ hybridization and physical mapping of YACs, both GPR32 and psi GPR32 were mapped to chromosomal 19, region q13.3, while psi GPR33 was mapped to chromosome 14q12.

Transposon tagging of the Defective embryo and meristems gene of tomato

The shoot and root apical meristems (SAMs and RAMs, respectively) of higher plants are mechanistically and structurally similar. This has led previously to the suggestion that the SAM and RAM represent modifications of a fundamentally homologous plan of organization. Despite recent interest in plant development, especially in the areas of meristem regulation, genes specifically required for the function of both the SAM and RAM have not yet been identified. Here, we report on a novel gene, Defective embryo and meristems (Dem), of tomato. This gene is required for the correct organization of shoot apical tissues of developing embryos, SAM development, and correct cell division patterns and meristem maintenance in roots. Dem was cloned using transposon tagging and shown to encode a novel protein of 72 kD with significant homology to YNV2, a protein of unknown function of Saccharomyces cerevisiae. Dem is expressed in root and shoot meristems and organ primordia but not in callus. The expression pattern of Dem mRNA in combination with the dem mutant phenotype suggests that Dem plays an important role within apical meristems.

Identification of a human member of the Ly-49 multigene family

Three classes of multigene family-encoded receptors enable NK cells to discriminate between polymorphic MHC class I molecules: Ly-49 homodimers, CD94/NKG2 heterodimers and the killer cell inhibitory receptors (KIR). Of these, CD94/NKG2 has been characterized in both rodents and humans. In contrast, Ly-49 family members have hitherto been found only in rodents, and KIR molecules only in the human. In this report, we describe a human cDNA, termed Ly-49L, that constitutes the first human member of the Ly-49 multi-gene family. Compared with rodent Ly-49 molecules, the Ly-49L sequence contains a premature stop codon and predicts a truncated protein that lacks the distal part of a C-terminal lectin domain. Evidence is presented that the premature stop codon results from incomplete excision of the intron between the first two lectin domain exons. Splice variants predicting a full-size Ly-49L protein were not detected. As demonstrated by Northern blot analysis, Ly-49L was transcribed by IL-2-activated NK cells, but not by freshly isolated B or T cells. PCR screening of a 22-clone yeast artificial chromosome contig localized the LY49L locus to the human NK gene complex on chromosome 12p12-p13. Southern blot analysis of genomic DNA showed a simple pattern with a full-length Ly-49L probe at low stringency hybridization conditions, suggesting that Ly-49L may be the only human member of the Ly-49 multigene family.

The Arcanobacterium (Actinomyces) pyogenes plasmid pAP1 is a member of the pIJ101/pJV1 family of rolling circle replication plasmids

The 2.4-kb plasmid pAP1 from Arcanobacterium (Actinomyces) pyogenes had sequence similarity within the putative replication protein and double-stranded origin with the pIJ101/pJV1 family of plasmids. pJGS84, a derivative of pAP1 containing a kanamycin resistance gene, was able to replicate in Escherichia coli and Corynebacterium pseudotuberculosis, as well as in A. pyogenes. Detection of single-stranded DNA intermediates of pJGS84 replication suggested that this plasmid replicates by the rolling circle mechanism.

Neurexophilins form a conserved family of neuropeptide-like glycoproteins

Neurexophilin was discovered as a neuronal glycoprotein that is copurified with neurexin Ialpha during affinity chromatography on immobilized alpha-latrotoxin (Petrenko et al., 1996). We have now investigated how neurexophilin interacts with neurexins, whether it is post-translationally processed by site-specific cleavage similar to neuropeptides, and whether related neuropeptide-like proteins are expressed in brain. Our data show that mammalian brains contain four genes for neurexophilins the products of which share a common structure composed of five domains: an N-terminal signal peptide, a variable N-terminal domain, a highly conserved central domain that is N-glycosylated, a short linker region, and a conserved C-terminal domain that is cysteine-rich. When expressed in pheochromocytoma (PC12) cells with a replication-deficient adenovirus, neurexophilin 1 was rapidly N-glycosylated and then slowly processed to a smaller mature form, probably by endoproteolytic cleavage. Similar expression experiments in other neuron-like cells and in fibroblastic cells revealed that N-glycosylation of neurexophilin 1 occurred in all cell types tested, whereas proteolytic processing was observed only in neuron-like cells. Finally, only recombinant neurexin Ialpha and IIIalpha but not neurexin Ibeta interacted with neurexophilin 1 and were preferentially bound to the processed mature form of neurexophilin. Together our data demonstrate that neurexophilins form a family of related glycoproteins that are proteolytically processed after synthesis and bind to alpha-neurexins. The structure and characteristics of neurexophilins indicate that they function as neuropeptides that may signal via alpha-neurexins.

Association of unconventional myosin MYO15 mutations with human nonsyndromic deafness DFNB3

DFNB3, a locus for nonsyndromic sensorineural recessive deafness, maps to a 3-centimorgan interval on human chromosome 17p11.2, a region that shows conserved synteny with mouse shaker-2. A human unconventional myosin gene, MYO15, was identified by combining functional and positional cloning approaches in searching for shaker-2 and DFNB3. MYO15 has at least 50 exons spanning 36 kilobases. Sequence analyses of these exons in affected individuals from three unrelated DFNB3 families revealed two missense mutations and one nonsense mutation that cosegregated with congenital recessive deafness.

SMART, a simple modular architecture research tool: identification of signaling domains

Accurate multiple alignments of 86 domains that occur in signaling proteins have been constructed and used to provide a Web-based tool (SMART: simple modular architecture research tool) that allows rapid identification and annotation of signaling domain sequences. The majority of signaling proteins are multidomain in character with a considerable variety of domain combinations known. Comparison with established databases showed that 25% of our domain set could not be deduced from SwissProt and 41% could not be annotated by Pfam. SMART is able to determine the modular architectures of single sequences or genomes; application to the entire yeast genome revealed that at least 6.7% of its genes contain one or more signaling domains, approximately 350 greater than previously annotated. The process of constructing SMART predicted (i) novel domain homologues in unexpected locations such as band 4.1-homologous domains in focal adhesion kinases; (ii) previously unknown domain families, including a citron-homology domain; (iii) putative functions of domain families after identification of additional family members, for example, a ubiquitin-binding role for ubiquitin-associated domains (UBA); (iv) cellular roles for proteins, such predicted DEATH domains in netrin receptors further implicating these molecules in axonal guidance; (v) signaling domains in known disease genes such as SPRY domains in both marenostrin/pyrin and Midline 1; (vi) domains in unexpected phylogenetic contexts such as diacylglycerol kinase homologues in yeast and bacteria; and (vii) likely protein misclassifications exemplified by a predicted pleckstrin homology domain in a Candida albicans protein, previously described as an integrin.

Highly specific protein sequence motifs for genome analysis

We present a method for discovering conserved sequence motifs from families of aligned protein sequences. The method has been implemented as a computer program called EMOTIF (http://motif. stanford.edu/emotif). Given an aligned set of protein sequences, EMOTIF generates a set of motifs with a wide range of specificities and sensitivities. EMOTIF also can generate motifs that describe possible subfamilies of a protein superfamily. A disjunction of such motifs often can represent the entire superfamily with high specificity and sensitivity. We have used EMOTIF to generate sets of motifs from all 7,000 protein alignments in the BLOCKS and PRINTS databases. The resulting database, called IDENTIFY (http://motif. stanford.edu/identify), contains more than 50,000 motifs. For each alignment, the database contains several motifs having a probability of matching a false positive that range from 10(-10) to 10(-5). Highly specific motifs are well suited for searching entire proteomes, while generating very few false predictions. IDENTIFY assigns biological functions to 25-30% of all proteins encoded by the Saccharomyces cerevisiae genome and by several bacterial genomes. In particular, IDENTIFY assigned functions to 172 of proteins of unknown function in the yeast genome.

Assessing sequence comparison methods with reliable structurally identified distant evolutionary relationships

Pairwise sequence comparison methods have been assessed using proteins whose relationships are known reliably from their structures and functions, as described in the SCOP database [Murzin, A. G., Brenner, S. E., Hubbard, T. & Chothia C. (1995) J. Mol. Biol. 247, 536-540]. The evaluation tested the programs BLAST [Altschul, S. F., Gish, W., Miller, W., Myers, E. W. & Lipman, D. J. (1990). J. Mol. Biol. 215, 403-410], WU-BLAST2 [Altschul, S. F. & Gish, W. (1996) Methods Enzymol. 266, 460-480], FASTA [Pearson, W. R. & Lipman, D. J. (1988) Proc. Natl. Acad. Sci. USA 85, 2444-2448], and SSEARCH [Smith, T. F. & Waterman, M. S. (1981) J. Mol. Biol. 147, 195-197] and their scoring schemes. The error rate of all algorithms is greatly reduced by using statistical scores to evaluate matches rather than percentage identity or raw scores. The E-value statistical scores of SSEARCH and FASTA are reliable: the number of false positives found in our tests agrees well with the scores reported. However, the P-values reported by BLAST and WU-BLAST2 exaggerate significance by orders of magnitude. SSEARCH, FASTA ktup = 1, and WU-BLAST2 perform best, and they are capable of detecting almost all relationships between proteins whose sequence identities are >30%. For more distantly related proteins, they do much less well; only one-half of the relationships between proteins with 20-30% identity are found. Because many homologs have low sequence similarity, most distant relationships cannot be detected by any pairwise comparison method; however, those which are identified may be used with confidence.

A unified statistical framework for sequence comparison and structure comparison

We present an approach for assessing the significance of sequence and structure comparisons by using nearly identical statistical formalisms for both sequence and structure. Doing so involves an all-vs.-all comparison of protein domains [taken here from the Structural Classification of Proteins (scop) database] and then fitting a simple distribution function to the observed scores. By using this distribution, we can attach a statistical significance to each comparison score in the form of a P value, the probability that a better score would occur by chance. As expected, we find that the scores for sequence matching follow an extreme-value distribution. The agreement, moreover, between the P values that we derive from this distribution and those reported by standard programs (e.g., BLAST and FASTA validates our approach. Structure comparison scores also follow an extreme-value distribution when the statistics are expressed in terms of a structural alignment score (essentially the sum of reciprocated distances between aligned atoms minus gap penalties). We find that the traditional metric of structural similarity, the rms deviation in atom positions after fitting aligned atoms, follows a different distribution of scores and does not perform as well as the structural alignment score. Comparison of the sequence and structure statistics for pairs of proteins known to be related distantly shows that structural comparison is able to detect approximately twice as many distant relationships as sequence comparison at the same error rate. The comparison also indicates that there are very few pairs with significant similarity in terms of sequence but not structure whereas many pairs have significant similarity in terms of structure but not sequence.

Measuring genome evolution

The determination of complete genome sequences provides us with an opportunity to describe and analyze evolution at the comprehensive level of genomes. Here we compare nine genomes with respect to their protein coding genes at two levels: (i) we compare genomes as "bags of genes" and measure the fraction of orthologs shared between genomes and (ii) we quantify correlations between genes with respect to their relative positions in genomes. Distances between the genomes are related to their divergence times, measured as the number of amino acid substitutions per site in a set of 34 orthologous genes that are shared among all the genomes compared. We establish a hierarchy of rates at which genomes have changed during evolution. Protein sequence identity is the most conserved, followed by the complement of genes within the genome. Next is the degree of conservation of the order of genes, whereas gene regulation appears to evolve at the highest rate. Finally, we show that some genomes are more highly organized than others: they show a higher degree of the clustering of genes that have orthologs in other genomes.

Impaired assembly of E1 decarboxylase of the branched-chain alpha-ketoacid dehydrogenase complex in type IA maple syrup urine disease

The E1 decarboxylase component of the human branched-chain ketoacid dehydrogenase complex comprises two E1alpha (45.5 kDa) and two E1beta (37.5 kDa) subunits forming an alpha2 beta2 tetramer. In patients with type IA maple syrup urine disease, the E1alpha subunit is affected, resulting in the loss of E1 and branched-chain ketoacid dehydrogenase catalytic activities. To study the effect of human E1alpha missense mutations on E1 subunit assembly, we have developed a pulse-chase labeling protocol based on efficient expression and assembly of human (His)6-E1alpha and untagged E1beta subunits in Escherichia coli in the presence of overexpressed chaperonins GroEL and GroES. Assembly of the two 35S-labeled E1 subunits was indicated by their co-extraction with Ni2+-nitrilotriacetic acid resin. The nine E1alpha maple syrup urine disease mutants studied showed aberrant kinetics of assembly with normal E1beta in the 2-h chase compared with the wild type and can be classified into four categories of normal (N222S-alpha and R220W-alpha), moderately slow (G245R-alpha), slow (G204S-alpha, A240P-alpha, F364C-alpha, Y368C-alpha, and Y393N-alpha), and no (T265R-alpha) assembly. Prolonged induction in E. coli grown in the YTGK medium or lowering of induction temperature from 37 to 28 degreesC (in the case of T265R-alpha), however, resulted in the production of mutant E1 proteins. Separation of purified E1 proteins by sucrose density gradient centrifugation showed that the wild-type E1 existed entirely as alpha2 beta2 tetramers. In contrast, a subset of E1alpha missense mutations caused the occurrence of exclusive alphabeta dimers (Y393N-alpha and F364C-alpha) or of both alpha2beta2 tetramers and lower molecular weight species (Y368C-alpha and T265R-alpha). Thermal denaturation at 50 degreesC indicated that mutant E1 proteins aggregated more rapidly than wild type (rate constant, 0.19 min-1), with the T265R-alpha mutant E1 most severely affected (rate constant, 4.45 min-1). The results establish that the human E1alpha mutations in the putative thiamine pyrophosphate-binding pocket that are studied, with the exception of G204S-alpha, have no effect on E1 subunit assembly. The T265R-alpha mutation adversely impacts both E1alpha folding and subunit interactions. The mutations involving the C-terminal aromatic residues impede both the kinetics of subunit assembly and the formation of the native alpha2 beta2 structure.

Neurexophilins form a conserved family of neuropeptide-like glycoproteins

Neurexophilin was discovered as a neuronal glycoprotein that is copurified with neurexin Ialpha during affinity chromatography on immobilized alpha-latrotoxin (Petrenko et al., 1996). We have now investigated how neurexophilin interacts with neurexins, whether it is post-translationally processed by site-specific cleavage similar to neuropeptides, and whether related neuropeptide-like proteins are expressed in brain. Our data show that mammalian brains contain four genes for neurexophilins the products of which share a common structure composed of five domains: an N-terminal signal peptide, a variable N-terminal domain, a highly conserved central domain that is N-glycosylated, a short linker region, and a conserved C-terminal domain that is cysteine-rich. When expressed in pheochromocytoma (PC12) cells with a replication-deficient adenovirus, neurexophilin 1 was rapidly N-glycosylated and then slowly processed to a smaller mature form, probably by endoproteolytic cleavage. Similar expression experiments in other neuron-like cells and in fibroblastic cells revealed that N-glycosylation of neurexophilin 1 occurred in all cell types tested, whereas proteolytic processing was observed only in neuron-like cells. Finally, only recombinant neurexin Ialpha and IIIalpha but not neurexin Ibeta interacted with neurexophilin 1 and were preferentially bound to the processed mature form of neurexophilin. Together our data demonstrate that neurexophilins form a family of related glycoproteins that are proteolytically processed after synthesis and bind to alpha-neurexins. The structure and characteristics of neurexophilins indicate that they function as neuropeptides that may signal via alpha-neurexins.

Neuropathy target esterase and a homologous Drosophila neurodegeneration-associated mutant protein contain a novel domain conserved from bacteria to man

The N-terminal amino acid sequences of proteolytic fragments of neuropathy target esterase (NTE), covalently labelled on its active-site serine by a biotinylated organophosphorus ester, were determined and used to deduce the location of this serine residue and to initiate cloning of its cDNA. A putative NTE clone, isolated from a human foetal brain cDNA library, encoded a 1327 residue polypeptide with no homology to any known serine esterases or proteases. The active-site serine of NTE (Ser-966) lay in the centre of a predicted hydrophobic helix within a 200-amino-acid C-terminal domain with marked similarity to conceptual proteins in bacteria, yeast and nematodes; these proteins may comprise a novel family of potential serine hydrolases. The Swiss Cheese protein which, when mutated, leads to widespread cell death in Drosophila brain [Kretzschmar, Hasan, Sharma, Heisenberg and Benzer (1997) J. Neurosci. 17, 7425-7432], was strikingly homologous to NTE, suggesting that genetically altered NTE may be involved in human neurodegenerative disease.

Function of Escherichia coli MsbA, an essential ABC family transporter, in lipid A and phospholipid biosynthesis

The Escherichia coli msbA gene, first identified as a multicopy suppressor of htrB mutations, has been proposed to transport nascent core-lipid A molecules across the inner membrane (Polissi, A., and Georgopoulos, C. (1996) Mol. Microbiol. 20, 1221-1233). msbA is an essential E. coli gene with high sequence similarity to mammalian Mdr proteins and certain types of bacterial ABC transporters. htrB is required for growth above 32 degreesC and encodes the lauroyltransferase that acts after Kdo addition during lipid A biosynthesis (Clementz, T., Bednarski, J., and Raetz, C. R. H. (1996) J. Biol. Chem. 271, 12095-12102). By using a quantitative new 32Pi labeling technique, we demonstrate that hexa-acylated species of lipid A predominate in the outer membranes of wild type E. coli labeled for several generations at 42 degreesC. In contrast, in htrB mutants shifted to 42 degreesC for 3 h, tetra-acylated lipid A species and glycerophospholipids accumulate in the inner membrane. Extra copies of the cloned msbA gene restore the ability of htrB mutants to grow at 42 degreesC, but they do not increase the extent of lipid A acylation. However, a significant fraction of the tetra-acylated lipid A species that accumulate in htrB mutants are transported to the outer membrane in the presence of extra copies of msbA. E. coli strains in which msbA synthesis is selectively shut off at 42 degreesC accumulate hexa-acylated lipid A and glycerophospholipids in their inner membranes. Our results support the view that MsbA plays a role in lipid A and possibly glycerophospholipid transport. The tetra-acylated lipid A precursors that accumulate in htrB mutants may not be transported as efficiently by MsbA as are penta- or hexa-acylated lipid A species.