Metaxin, a gene contiguous to both thrombospondin 3 and glucocerebrosidase, is required for embryonic development in the mouse: implications for Gaucher disease.

We have identified a murine gene, metaxin, that spans the 6-kb interval separating the glucocerebrosidase gene (GC) from the thrombospondin 3 gene on chromosome 3E3-F1. Metaxin and GC are transcribed convergently; their major polyadenylylation sites are only 431 bp apart. On the other hand, metaxin and the thrombospondin 3 gene are transcribed divergently and share a common promoter sequence. The cDNA for metaxin encodes a 317-aa protein, without either a signal sequence or consensus for N-linked glycosylation. Metaxin protein is expressed ubiquitously in tissues of the young adult mouse, but no close homologues have been found in the DNA or protein data bases. A targeted mutation (A-->G in exon 9) was introduced into GC by homologous recombination in embryonic stem cells to establish a mouse model for a mild form of Gaucher disease. A phosphoglycerate kinase-neomycin gene cassette was also inserted into the 3'-flanking region of GC as a selectable marker, at a site later identified as the terminal exon of metaxin. Mice homozygous for the combined mutations die early in gestation. Since the same amino acid mutation in humans is associated with mild type 1 Gaucher disease, we suggest that metaxin protein is likely to be essential for embryonic development in mice. Clearly, the contiguous gene organization at this locus limits targeting strategies for the production of murine models of Gaucher disease.

Three functional luciferase domains in a single polypeptide chain

We report a unique case of a gene containing three homologous and contiguous repeat sequences, each of which, after excision, cloning, and expression in Escherichia coli, is shown to code for a peptide catalyzing the same reaction as the native protein, Gonyaulax polyedra luciferase (Mr = 137). This enzyme, which catalyzes the light-emitting oxidation of a linear tetrapyrrole (dinoflagellate luciferin), exhibits no sequence similarities to other luciferases in databases. Sequence analysis also reveals an unusual evolutionary feature of this gene: synonymous substitutions are strongly constrained in the central regions of each of the repeated coding sequences.

A complete genome screen for genes predisposing to severe bipolar disorder in two Costa Rican pedigrees

Bipolar mood disorder (BP) is a debilitating syndrome characterized by episodes of mania and depression. We designed a multistage study to detect all major loci predisposing to severe BP (termed BP-I) in two pedigrees drawn from the Central Valley of Costa Rica, where the population is largely descended from a few founders in the 16th–18th centuries. We considered only individuals with BP-I as affected and screened the genome for linkage with 473 microsatellite markers. We used a model for linkage analysis that incorporated a high phenocopy rate and a conservative estimate of penetrance. Our goal in this study was not to establish definitive linkage but rather to detect all regions possibly harboring major genes for BP-I in these pedigrees. To facilitate this aim, we evaluated the degree to which markers that were informative in our data set provided coverage of each genome region; we estimate that at least 94% of the genome has been covered, at a predesignated threshold determined through prior linkage simulation analyses. We report here the results of our genome screen for BP-I loci and indicate several regions that merit further study, including segments in 18q, 18p, and 11p, in which suggestive lod scores were observed for two or more contiguous markers. Isolated lod scores that exceeded our thresholds in one or both families also occurred on chromosomes 1, 2, 3, 4, 5, 7, 13, 15, 16, and 17. Interesting regions highlighted in this genome screen will be followed up using linkage disequilibrium (LD) methods.

A widely expressed βIII spectrin associated with Golgi and cytoplasmic vesicles

Spectrin is an important structural component of the plasma membrane skeleton. Heretofore-unidentified isoforms of spectrin also associate with Golgi and other organelles. We have discovered another member of the β-spectrin gene family by homology searches of the GenBank databases and by 5′ rapid amplification of cDNA ends of human brain cDNAs. Collectively, 7,938 nucleotides of contiguous clones are predicted to encode a 271,294-Da protein, called βIII spectrin, with conserved actin-, protein 4.1-, and ankyrin-binding domains, membrane association domains 1 and 2, a spectrin dimer self-association site, and a pleckstrin-homology domain. βIII spectrin transcripts are concentrated in the brain and present in the kidneys, liver, and testes and the prostate, pituitary, adrenal, and salivary glands. All of the tested tissues contain major 9.0-kb and minor 11.3-kb transcripts. The human βIII spectrin gene (SPTBN2) maps to chromosome 11q13 and the mouse gene (Spnb3) maps to a syntenic region close to the centromere on chromosome 19. Indirect immunofluorescence studies of cultured cells using antisera specific to human βIII spectrin reveal a Golgi-associated and punctate cytoplasmic vesicle-like distribution, suggesting that βIII spectrin associates with intracellular organelles. This distribution overlaps that of several Golgi and vesicle markers, including mannosidase II, p58, trans-Golgi network (TGN)38, and β-COP and is distinct from the endoplasmic reticulum markers calnexin and Bip. Liver Golgi membranes and other vesicular compartment markers cosediment in vitro with βIII spectrin. βIII spectrin thus constitutes a major component of the Golgi and vesicular membrane skeletons.

Localization of the N terminus of hepatitis B virus capsid protein by peptide-based difference mapping from cryoelectron microscopy

Recently, cryoelectron microscopy of isolated macromolecular complexes has advanced to resolutions below 10 Å, enabling direct visualization of α-helical secondary structure. To help correlate such density maps with the amino acid sequences of the component proteins, we advocate peptide-based difference mapping, i.e., insertion of peptides, ≈10 residues long, at targeted points in the sequence and visualization of these peptides as bulk labels in cryoelectron microscopy-derived difference maps. As proof of principle, we have appended an extraneous octapeptide at the N terminus of hepatitis B virus capsid protein and determined its location on the capsid surface by difference imaging at 11 Å resolution. Hepatitis B virus capsids are icosahedral particles, ≈300 Å in diameter, made up of T-shaped dimers (subunit Mr, 16–21 kDa, depending on construct). The stems of the Ts protrude outward as spikes, whereas the crosspieces pack to form the contiguous shell. The two N termini per dimer reside on either side of the spike-stem, at the level at which it enters the shell. This location is consistent with formation of the known intramolecular disulfide bond between the cysteines at positions 61 and −7 (in the residual propeptide) in the “e-antigen” form of the capsid protein and has implications for why this clinically important antigen remains unassembled in vivo.

Synthetic genes for glycoprotein design and the elucidation of hydroxyproline-O-glycosylation codes

Design of hydroxyproline (Hyp)-rich glycoproteins (HRGPs) offers an approach for the structural and functional analysis of these wall components, which are broadly implicated in plant growth and development. HRGPs consist of multiple small repetitive “glycomodules” extensively O-glycosylated through the Hyp residues. The patterns of Hyp-O-glycosylation are putatively coded by the primary sequence as described by the Hyp contiguity hypothesis, which predicts contiguous Hyp residues to be attachment sites of small arabinooligosaccharides (1–5 Ara residues/Hyp); while clustered, noncontiguous Hyp residues are sites of arabinogalactan polysaccharide attachment. As a test, we designed two simple HRGPs as fusion proteins with green fluorescent protein. The first was a repetitive Ser-Hyp motif that encoded only clustered noncontiguous Hyp residues, predicted polysaccharide addition sites. The resulting glycoprotein had arabinogalactan polysaccharide O-linked to all Hyp residues. The second construct, based on the consensus sequence of a gum arabic HRGP, contained both arabinogalactan and arabinooligosaccharide addition sites and, as predicted, gave a product that contained both saccharide types. These results identify an O-glycosylation code of plants.

Transgenic rats reveal functional conservation of regulatory controls between the Fugu isotocin and rat oxytocin genes

We have asked whether comparative genome analysis and rat transgenesis can be used to identify functional regulatory domains in the gene locus encoding the hypothalamic neuropeptides oxytocin (OT) and vasopressin. Isotocin (IT) and vasotocin (VT) are the teleost homologues of these genes. A contiguous stretch of 46 kb spanning the Fugu IT-VT locus has been sequenced, and nine putative genes were found. Unlike the OT and vasopressin genes, which are closely linked in the mammalian genome in a tail-to-tail orientation, Fugu IT and VT genes are linked head to tail and are separated by five genes. When a cosmid containing the Fugu IT-VT locus was introduced into the rat genome, we found that the Fugu IT gene was specifically expressed in rat hypothalamic oxytocinergic neurons and mimicked the response of the endogenous OT gene to an osmotic stimulus. These data show that cis-acting elements and trans-acting factors mediating the cell-specific and physiological regulation of the OT and IT genes are conserved between mammals and fish. The combination of Fugu genome analysis and transgenesis in a mammal is a powerful tool for identifying and analyzing conserved vertebrate regulatory elements.

In vitro reconstruction of the Aspergillus (= Emericella) nidulans genome

A physical map of the 31-megabase Aspergillus nidulans genome is reported, in which 94% of 5,134 cosmids are assigned to 49 contiguous segments. The physical map is the result of a two-way ordering process, in which clones and probes were ordered simultaneously on a binary DNA/DNA hybridization matrix. Compression by elimination of redundant clones resulted in a minimal map, which is a chromosome walk. Repetitive DNA is nonrandomly dispersed in the A. nidulans genome, reminiscent of heterochromatic banding patterns of higher eukaryotes. We hypothesize gene clusters may arise by horizontal transfer and spread by transposition to explain the nonrandom pattern of repeats along chromosomes.

Ancient mtDNA sequences in the human nuclear genome: A potential source of errors in identifying pathogenic mutations

Nuclear-localized mtDNA pseudogenes might explain a recent report describing a heteroplasmic mtDNA molecule containing five linked missense mutations dispersed over the contiguous mtDNA CO1 and CO2 genes in Alzheimer’s disease (AD) patients. To test this hypothesis, we have used the PCR primers utilized in the original report to amplify CO1 and CO2 sequences from two independent ρ° (mtDNA-less) cell lines. CO1 and CO2 sequences amplified from both of the ρ° cells, demonstrating that these sequences are also present in the human nuclear DNA. The nuclear pseudogene CO1 and CO2 sequences were then tested for each of the five “AD” missense mutations by restriction endonuclease site variant assays. All five mutations were found in the nuclear CO1 and CO2 PCR products from ρ° cells, but none were found in the PCR products obtained from cells with normal mtDNA. Moreover, when the overlapping nuclear CO1 and CO2 PCR products were cloned and sequenced, all five missense mutations were found, as well as a linked synonymous mutation. Unlike the findings in the original report, an additional 32 base substitutions were found, including two in adjacent tRNAs and a two base pair deletion in the CO2 gene. Phylogenetic analysis of the nuclear CO1 and CO2 sequences revealed that they diverged from modern human mtDNAs early in hominid evolution about 770,000 years before present. These data would be consistent with the interpretation that the missense mutations proposed to cause AD may be the product of ancient mtDNA variants preserved as nuclear pseudogenes.

A Discrete Stage of Baculovirus GP64-mediated Membrane Fusion

Viral fusion protein trimers can play a critical role in limiting lipids in membrane fusion. Because the trimeric oligomer of many viral fusion proteins is often stabilized by hydrophobic 4-3 heptad repeats, higher-order oligomers might be stabilized by similar sequences. There is a hydrophobic 4-3 heptad repeat contiguous to a putative oligomerization domain of Autographa californica multicapsid nucleopolyhedrovirus envelope glycoprotein GP64. We performed mutagenesis and peptide inhibition studies to determine if this sequence might play a role in catalysis of membrane fusion. First, leucine-to-alanine mutants within and flanking the amino terminus of the hydrophobic 4-3 heptad repeat motif that oligomerize into trimers and traffic to insect Sf9 cell surfaces were identified. These mutants retained their wild-type conformation at neutral pH and changed conformation in acidic conditions, as judged by the reactivity of a conformationally sensitive mAb. These mutants, however, were defective for membrane fusion. Second, a peptide encoding the portion flanking the GP64 hydrophobic 4-3 heptad repeat was synthesized. Adding peptide led to inhibition of membrane fusion, which occurred only when the peptide was present during low pH application. The presence of peptide during low pH application did not prevent low pH–induced conformational changes, as determined by the loss of a conformationally sensitive epitope. In control experiments, a peptide of identical composition but different sequence, or a peptide encoding a portion of the Ebola GP heptad motif, had no effect on GP64-mediated fusion. Furthermore, when the hemagglutinin (X31 strain) fusion protein of influenza was functionally expressed in Sf9 cells, no effect on hemagglutinin-mediated fusion was observed, suggesting that the peptide does not exert nonspecific effects on other fusion proteins or cell membranes. Collectively, these studies suggest that the specific peptide sequences of GP64 that are adjacent to and include portions of the hydrophobic 4-3 heptad repeat play a dynamic role in membrane fusion at a stage that is downstream of the initiation of protein conformational changes but upstream of lipid mixing.

A general method for determining helix packing in membrane proteins in situ: Helices I and II are close to helix VII in the lactose permease of Escherichia coli

It was previously shown that coexpression of the lactose permease of Escherichia coli in two contiguous fragments leads to functional complementation. We demonstrate here that site-directed thiol crosslinking of coexpressed permease fragments can be used to determine helix proximity in situ without the necessity of purifying the permease. After coexpression of the six N-terminal (N6) and six C-terminal (C6) transmembrane helices, each with a single Cys residue, crosslinking was carried out in native membranes and assessed by the mobility of anti-C-terminal-reactive polypeptides on immunoblots. A Cys residue at position 242 or 245 (helix VII) forms a disulfide with a Cys residue at either position 28 or 29 (helix I), but not with a Cys residue at position 27, which is on the opposite face of helix I, thereby indicating that the face of helix I containing Pro-28 and Phe-29 is close to helix VII. Similarly, a Cys residue at position 242 or 245 (helix VII) forms a disulfide with a Cys residue at either position 52 or 53 (helix II), but not with a Cys residue at position 54. Furthermore, low-efficiency crosslinking is observed between a Cys residue at position 52 or 53 and a Cys residue at position 361 (helix XI). The results indicate that helix VII lies in close proximity to both helices I and II and that helix II is also close to helix XI. The method should be applicable to a number of different polytopic membrane proteins.

A regulator of G protein signaling interaction surface linked to effector specificity

Proteins of the regulator of G protein signaling (RGS) family accelerate GTP hydrolysis by the α subunits (Gα) of G proteins, leading to rapid recovery of signaling cascades. Many different RGS proteins can accelerate GTP hydrolysis by an individual Gα, and GTP hydrolysis rates of different Gαs can be enhanced by the same RGS protein. Consequently, the mechanisms for specificity in RGS regulation and the residues involved remain unclear. Using the evolutionary trace (ET) method, we have identified a cluster of residues in the RGS domain that includes the RGS-Gα binding interface and extends to include additional functionally important residues on the surface. One of these is within helix α3, two are in α5, and three are in the loop connecting α5 and α6. A cluster of surface residues on Gα previously identified by ET, and composed predominantly of residues from the switch III region and helix α3, is spatially contiguous with the ET-identified residues in the RGS domain. This cluster includes residues proposed to interact with the γ subunit of Gtα's effector, cGMP phosphodiesterase (PDEγ). The proximity of these clusters suggests that they form part of an interface between the effector and the RGS-Gα complex. Sequence variations in these residues correlate with PDEγ effects on GTPase acceleration. Because ET identifies residues important for all members of a protein family, these residues likely form a general site for regulation of G protein-coupled signaling cascades, possibly by means of effector interactions.

Designed protein tetramer zipped together with a hydrophobic Alzheimer homology: A structural clue to amyloid assembly

Limited solubility and precipitation of amyloidogenic sequences such as the Alzheimer peptide (β-AP) are major obstacles to a molecular understanding of protein fibrillation and deposition processes. Here we have circumvented the solubility problem by stepwise engineering a β-AP homology into a soluble scaffold, the monomeric protein S6. The S6 construct with the highest β-AP homology crystallizes as a tetramer that is linked by the β-AP residues forming intermolecular antiparallel β-sheets. This construct also shows increased coil aggregation during refolding, and a 14-mer peptide encompassing the engineered sequence forms fibrils. Mutational analysis shows that intermolecular association is linked to the overall hydrophobicity of the sticky sequence and implies the existence of “structural gatekeepers” in the wild-type protein, that is, charged side chains that prevent aggregation by interrupting contiguous stretches of hydrophobic residues in the primary sequence.

Molecular characterization of the phenylacetic acid catabolic pathway in Pseudomonas putida U: The phenylacetyl-CoA catabolon

Fourteen different genes included in a DNA fragment of 18 kb are involved in the aerobic degradation of phenylacetic acid by Pseudomonas putida U. This catabolic pathway appears to be organized in three contiguous operons that contain the following functional units: (i) a transport system, (ii) a phenylacetic acid activating enzyme, (iii) a ring-hydroxylation complex, (iv) a ring-opening protein, (v) a β-oxidation-like system, and (vi) two regulatory genes. This pathway constitutes the common part (core) of a complex functional unit (catabolon) integrated by several routes that catalyze the transformation of structurally related molecules into a common intermediate (phenylacetyl-CoA).

Recent improvements in prediction of protein structure by global optimization of a potential energy function

Recent improvements of a hierarchical ab initio or de novo approach for predicting both α and β structures of proteins are described. The united-residue energy function used in this procedure includes multibody interactions from a cumulant expansion of the free energy of polypeptide chains, with their relative weights determined by Z-score optimization. The critical initial stage of the hierarchical procedure involves a search of conformational space by the conformational space annealing (CSA) method, followed by optimization of an all-atom model. The procedure was assessed in a recent blind test of protein structure prediction (CASP4). The resulting lowest-energy structures of the target proteins (ranging in size from 70 to 244 residues) agreed with the experimental structures in many respects. The entire experimental structure of a cyclic α-helical protein of 70 residues was predicted to within 4.3 Å α-carbon (Cα) rms deviation (rmsd) whereas, for other α-helical proteins, fragments of roughly 60 residues were predicted to within 6.0 Å Cα rmsd. Whereas β structures can now be predicted with the new procedure, the success rate for α/β- and β-proteins is lower than that for α-proteins at present. For the β portions of α/β structures, the Cα rmsd's are less than 6.0 Å for contiguous fragments of 30–40 residues; for one target, three fragments (of length 10, 23, and 28 residues, respectively) formed a compact part of the tertiary structure with a Cα rmsd less than 6.0 Å. Overall, these results constitute an important step toward the ab initio prediction of protein structure solely from the amino acid sequence.