Sequence of the 1918 pandemic influenza virus nonstructural gene (NS) segment and characterization of recombinant viruses bearing the 1918 NS genes

The influenza A virus pandemic of 1918���1919 resulted in an estimated 20–40 million deaths worldwide. The hemagglutinin and neuraminidase sequences of the 1918 virus were previously determined. We here report the sequence of the A/Brevig Mission/1/18 (H1N1) virus nonstructural (NS) segment encoding two proteins, NS1 and nuclear export protein. Phylogenetically, these genes appear to be close to the common ancestor of subsequent human and classical swine strain NS genes. Recently, the influenza A virus NS1 protein was shown to be a type I IFN antagonist that plays an important role in viral pathogenesis. By using the recently developed technique of generating influenza A viruses entirely from cloned cDNAs, the hypothesis that the 1918 virus NS1 gene played a role in virulence was tested in a mouse model. In a BSL3+ laboratory, viruses were generated that possessed either the 1918 NS1 gene alone or the entire 1918 NS segment in a background of influenza A/WSN/33 (H1N1), a mouse-adapted virus derived from a human influenza strain first isolated in 1933. These 1918 NS viruses replicated well in tissue culture but were attenuated in mice as compared with the isogenic control viruses. This attenuation in mice may be related to the human origin of the 1918 NS1 gene. These results suggest that interaction of the NS1 protein with host-cell factors plays a significant role in viral pathogenesis.

Mutagenesis of the Glucose-1-Phosphate-Binding Site of Potato Tuber ADP-Glucose Pyrophosphorylase1

Lysine (Lys)-195 in the homotetrameric ADP-glucose pyrophosphorylase (ADPGlc PPase) from Escherichia coli was shown previously to be involved in the binding of the substrate glucose-1-phosphate (Glc-1-P). This residue is highly conserved in the ADPGlc PPase family. Site-directed mutagenesis was used to investigate the function of this conserved Lys residue in the large and small subunits of the heterotetrameric potato (Solanum tuberosum) tuber enzyme. The apparent affinity for Glc-1-P of the wild-type enzyme decreased 135- to 550-fold by changing Lys-198 of the small subunit to arginine, alanine, or glutamic acid, suggesting that both the charge and the size of this residue influence Glc-1-P binding. These mutations had little effect on the kinetic constants for the other substrates (ATP and Mg2+ or ADP-Glc and inorganic phosphate), activator (3-phosphoglycerate), inhibitor (inorganic phosphate), or on the thermal stability. Mutagenesis of the corresponding Lys (Lys-213) in the large subunit had no effect on the apparent affinity for Glc-1-P by substitution with arginine, alanine, or glutamic acid. A double mutant, SK198RLK213R, was also obtained that had a 100-fold reduction of the apparent affinity for Glc-1-P. The data indicate that Lys-198 in the small subunit is directly involved in the binding of Glc-1-P, whereas they appear to exclude a direct role of Lys-213 in the large subunit in the interaction with this substrate.

Molecular characterization of the mycobacterial heparin-binding hemagglutinin, a mycobacterial adhesin

Although it generally is accepted that the interaction of Mycobacterium tuberculosis with alveolar macrophages is a key step in the pathogenesis of tuberculosis, interactions with other cell types, especially epithelial cells, also may be important. In this study we describe the molecular characterization of a mycobacterial heparin-binding hemagglutinin (HBHA), a protein that functions as an adhesin for epithelial cells. The structural gene was cloned from M. tuberculosis and bacillus Calmette–Guérin, and the sequence was found to be identical between the two species. The calculated Mr was smaller than the observed Mr when analyzed by SDS/PAGE. This difference can be attributed to the Lys/Pro-rich repeats that occur at the C-terminal end of the protein and to a putative carbohydrate moiety. Glycosylation of HBHA appears to protect the protein from proteolytic degradation, which results in the removal of the C-terminal Lys/Pro-rich region responsible for binding of HBHA to sulfated carbohydrates. Evidence suggests that glycosylation is also important for HBHA-mediated hemagglutination and for certain immunologic properties of the protein. Finally, the absence of a signal peptide in the coding region of HBHA raises the possibility that this protein is not secreted via the general secretion pathway.

A glutamate residue in the catalytic center of the yeast chorismate mutase restricts enzyme activity to acidic conditions

Chorismate mutase acts at the first branchpoint of aromatic amino acid biosynthesis and catalyzes the conversion of chorismate to prephenate. Comparison of the x-ray structures of allosteric chorismate mutase from the yeast Saccharomyces cerevisiae with Escherichia coli chorismate mutase/prephenate dehydratase suggested conserved active sites between both enzymes. We have replaced all critical amino acid residues, Arg-16, Arg-157, Lys-168, Glu-198, Thr-242, and Glu-246, of yeast chorismate mutase by aliphatic amino acid residues. The resulting enzymes exhibit the necessity of these residues for catalytic function and provide evidence of their localization at the active site. Unlike some bacterial enzymes, yeast chorismate mutase has highest activity at acidic pH values. Replacement of Glu-246 in the yeast chorismate mutase by glutamine changes the pH optimum for activity of the enzyme from a narrow to a broad pH range. These data suggest that Glu-246 in the catalytic center must be protonated for maximum catalysis and restricts optimal activity of the enzyme to low pH.

Activation of the c-Jun N-terminal kinase pathway by a novel protein kinase related to human germinal center kinase

The c-Jun N-terminal kinase (JNK), or stress-activated protein kinase plays a crucial role in cellular responses stimulated by environmental stress and proinflammatory cytokines. However, the mechanisms that lead to the activation of the JNK pathway have not been elucidated. We have isolated a cDNA encoding a novel protein kinase that has significant sequence similarities to human germinal center kinase (GCK) and human hematopoietic progenitor kinase 1. The novel GCK-like kinase (GLK) has a nucleotide sequence that encodes an ORF of 885 amino acids with 11 kinase subdomains. Endogenous GLK could be activated by UV radiation and proinflammatory cytokine tumor necrosis factor α. When transiently expressed in 293 cells, GLK specifically activated the JNK, but not the p42/44MAPK/extracellular signal-regulated kinase or p38 kinase signaling pathways. Interestingly, deletion of amino acids 353–835 in the putative C-terminal regulatory region, or mutation of Lys-35 in the putative ATP-binding domain, markedly reduced the ability of GLK to activate JNK. This result indicates that both kinase activity and the C-terminal region of GLK are required for maximal activation of JNK. Furthermore, GLK-induced JNK activation could be inhibited by a dominant-negative mutant of mitogen-activated protein kinase kinase kinase 1 (MEKK1) or mitogen-activated protein kinase kinase 4/SAPK/ERK kinase 1 (SEK1), suggesting that GLK may function upstream of MEKK1 in the JNK signaling pathway.

Yeast chorismate mutase in the R state: Simulations of the active site

The isomerization of chorismate to prephenate by chorismate mutase in the biosynthetic pathway that forms Tyr and Phe involves C5—O (ether) bond cleavage and C1—C9 bond formation in a Claisen rearrangement. Development of negative charge on the ether oxygen, stabilized by Lys-168 and Glu-246, is inferred from the structure of a complex with a transition state analogue (TSA) and from the pH-rate profile of the enzyme and the E246Q mutant. These studies imply a protonated Glu-246 well above pH 7. Here, several 500-ps molecular dynamics simulations test the stability of enzyme–TSA complexes by using a solvated system with stochastic boundary conditions. The simulated systems are (i) protonated Glu-246 (stable), (ii) deprotonated Glu-246 (unstable), (iii) deprotonated Glu-246 plus one H2O between Glu-246 and the ether oxygen (unstable), (iv) the E246Q mutant (stable), and (v) addition of OH− between protonated Glu-246 and the ether oxygen. In (v), a local conformational change of Lys-168 displaced the OH− into the solvent region, suggesting a possible rate-determining step that precedes the catalytic step. In a 500-ps simulation of the enzyme complexed with the reactant chorismate or the product prephenate, no water molecule remained near the oxygen of the ligand. Calculations using the linearized Poisson–Boltzmann equation show that the effective pKa of Glu-246 is shifted from 5.8 to 8.1 as the negative charge on the ether oxygen of the TSA is changed from −0.56 electron to −0.9 electron. Altogether, these results support retention of a proton on Glu-246 to high pH and the absence of a water molecule in the catalytic steps.

Dual role of the nuclear factor of activated T cells insert region in DNA recognition and cooperative contacts to activator protein 1

The transcription factors nuclear factor of activated T cells (NFAT) and activator protein 1 (AP-1) coordinately regulate cytokine gene expression in activated T-cells by binding to closely juxtaposed sites in cytokine promoters. The structural basis for cooperative binding of NFAT and AP-1 to these sites, and indeed for the cooperative binding of transcription factors to composite regulatory elements in general, is not well understood. Mutagenesis studies have identified a segment of AP-1, which lies at the junction of its DNA-binding and dimerization domains (basic region and leucine zipper, respectively), as being essential for protein–protein interactions with NFAT in the ternary NFAT/AP-1/DNA complex. In a model of the ternary complex, the segment of NFAT nearest AP-1 is the Rel insert region (RIR), a feature that is notable for its hypervariability in size and in sequence amongst members of the Rel transcription factor family. Here we have used mutational analysis to study the role of the NFAT RIR in binding to DNA and AP-1. Parallel yeast one-hybrid screening assays in combination with alanine-scanning mutagenesis led to the identification of four amino acid residues in the RIR of NFAT2 (also known as NFATC1 or NFATc) that are essential for cooperativity with AP-1 (Ile-544, Glu-545, Thr-551, and Ile-553), and three residues that are involved in interactions with DNA (Lys-538, Arg-540, and Asn-541). These results were confirmed and extended through in vitro binding assays. We thus conclude that the NFAT RIR plays an essential dual role in DNA recognition and cooperative binding to AP-1 family transcription factors.

Phosphorylation of Tyrosine 992, 1068, and 1086 Is Required for Conformational Change of the Human Epidermal Growth Factor Receptor C-Terminal Tail

We reported previously that a conformation-specific antibody, Ab P2, to a 16-amino acid peptide (Glu-Gly-Tyr-Lys-Lys-Lys-Tyr-Gln-Gln-Val-Asp-Glu-Glu-Phe-Leu-Arg) of the cytoplasmic domain of the β-type platelet-derived growth factor receptor also recognizes the epidermal growth factor (EGF) receptor. Although the antibody is not directed to phosphotyrosine, it recognizes in immunoprecipitation the activated and hence phosphorylated form of both receptors. In P2 peptide, there are two tripeptide sequences, Asp-Glu-Glu and Tyr-Gln-Gln, that are also present in the EGF receptor. Our present studies using either EGF receptor C-terminal deletion mutants or point mutations (Tyr→Phe) and our previous studies on antibody inhibition by P2-derived peptides suggest that Gln-Gln in combination with Asp-Glu-Glu forms a high-affinity complex with Ab P2 and that such complex formation is dependent on tyrosine phosphorylation. Of the five phosphate acceptor sites in the EGF receptor, clustered in the extreme C-terminal tail, phosphorylation of three tyrosine residues (992, 1068, and 1086) located between Asp-Glu-Glu and Gln-Gln is necessary for Ab P2 binding. In contrast, the acceptor sites Tyr 1173 and 1148 play no role in the conformation change. Asp-Glu-Glu and Gln-Gln are located 169 amino acids apart, and it is highly likely that the interactions among three negatively charged phosphotyrosine residues in the receptor C terminus may result in the bending of the peptide chain in such a way that these two peptides come close to each other to form an antibody-binding site. Such a possibility is also supported by our finding that receptor dephosphorylation results in complete loss of Ab P2–binding activity. In conclusion, we have identified a domain within the cytoplasmic part of the EGF receptor whose conformation is altered by receptor phosphorylation; furthermore, we have identified the tyrosine residues that positively regulate this conformation.

A Specific Point Mutant at Position 1 of the Influenza Hemagglutinin Fusion Peptide Displays a Hemifusion Phenotype

We showed previously that substitution of the first residue of the influenza hemagglutinin (HA) fusion peptide Gly1 with Glu abolishes fusion activity. In the present study we asked whether this striking phenotype was due to the charge or side-chain volume of the substituted Glu. To do this we generated and characterized six mutants with substitutions at position 1: Gly1 to Ala, Ser, Val, Glu, Gln, or Lys. We found the following. All mutants were expressed at the cell surface, could be cleaved from the precursor (HA0) to the fusion permissive form (HA1-S-S-HA2), bound antibodies against the major antigenic site, bound red blood cells, and changed conformation at low pH. Only Gly, Ala, and Ser supported lipid mixing during fusion with red blood cells. Only Gly and Ala supported content mixing. Ser HA, therefore, displayed a hemifusion phenotype. The hemifusion phenotype of Ser HA was confirmed by electrophysiological studies. Our findings indicate that the first residue of the HA fusion peptide must be small (e.g., Gly, Ala, or Ser) to promote lipid mixing and must be small and apolar (e.g., Gly or Ala) to support both lipid and content mixing. The finding that Val HA displays no fusion activity underscores the idea that hydrophobicity is not the sole factor dictating fusion peptide function. The surprising finding that Ser HA displays hemifusion suggests that the HA ectodomain functions not only in the first stage of fusion, lipid mixing, but also, either directly or indirectly, in the second stage of fusion, content mixing.

Studies on the role of the hydrophobic domain of Ost4p in interactions with other subunits of yeast oligosaccharyl transferase

In the yeast, Saccharomyces cerevisiae, oligosaccharyl transferase (OT), which catalyzes the transfer of dolichol-linked oligosaccharide chains to nascent polypeptides in the endoplasmic reticulum, consists of nine nonidentical membrane protein subunits. Genetic and biochemical evidence indicated these nine proteins exist in three subcomplexes. Three of the OT subunits (Ost4p, Ost3p, and Stt3p) have been proposed to exist in one subcomplex. To investigate the interaction of these three membrane proteins, initially we carried out a mutational analysis of Ost4p, which is an extraordinarily small membrane protein containing only 36 amino acid residues. This analysis indicated that when single amino acid residues in a region close to the luminal face of the putative transmembrane domain of Ost4p were changed into an ionizable amino acid such as Lys or Asp, growth at 37°C and OT activity measured in vitro were impaired. In addition, using immunoprecipitation techniques and Western blot analysis, we found that with these mutations the interaction between Ost4p, Ost3p, and Stt3p was disrupted. Introduction of Lys or Asp residues at other positions in the putative transmembrane domain or at the N or C terminus of Ost4p had no effect on disrupting subunit interactions or impairing the activity of OT. These findings suggest that a localized region of the putative transmembrane domain of Ost4p mediates in stabilization of the interaction with the two other OT subunits (Ost3p and Stt3p) in a subcomplex in the endoplasmic reticulum membrane.

A meanfield approach to the thermodynamics of a protein-solvent system with application to the oligomerization of the tumor suppressor p53

The thermodynamic stability and oligomerization status of the tumor suppressor p53 tetramerization domain have been studied experimentally and theoretically. A series of hydrophilic mutations at Met-340 and Leu-344 of human p53 were designed to disrupt the hydrophobic dimer–dimer interface of the tetrameric oligomerization domain of p53 (residues 325–355). Meanfield calculations of the free energy of the solvated mutants as a function of interdimer distance were compared with experimental data on the thermal stability and oligomeric state (tetramer, dimer, or equilibrium mixture of both) of each mutant. The calculations predicted a decreasing stability and oligomeric state for the following amino acids at residue 340: Met (tetramer) > Ser Asp, His, Gln, > Glu, Lys (dimer), whereas the experimental results showed the following order: Met (tetramer) > Ser > Gln > His, Lys > Asp, Glu (dimers). For residue 344, the calculated trend was Leu (tetramer) > Ala > Arg, Gln, Lys (dimer), and the experimental trend was Leu (tetramer) > Ala, Arg, Gln, Lys (dimer). The discrepancy for the lysine side chain at residue 340 is attributed to the dual nature of lysine, both hydrophobic and charged. The incorrect prediction of stability of the mutant with Asp at residue 340 is attributed to the fact that within the meanfield approach, we use the wild-type backbone configuration for all mutants, but low melting temperatures suggest a softening of the α-helices at the dimer–dimer interface. Overall, this initial application of meanfield theory toward a protein-solvent system is encouraging for the application of the theoretical model to more complex systems.

Catalytic mechanism of the adenylyl and guanylyl cyclases: Modeling and mutational analysis

The adenylyl and guanylyl cyclases catalyze the formation of 3′,5′-cyclic adenosine or guanosine monophosphate from the corresponding nucleoside 5′-triphosphate. The guanylyl cyclases, the mammalian adenylyl cyclases, and their microbial homologues function as pairs of homologous catalytic domains. The crystal structure of the rat type II adenylyl cyclase C2 catalytic domain was used to model by homology a mammalian adenylyl cyclase C1-C2 domain pair, a homodimeric adenylyl cyclase of Dictyostelium discoideum, a heterodimeric soluble guanylyl cyclase, and a homodimeric membrane guanylyl cyclase. Mg2+ATP or Mg2+GTP were docked into the active sites based on known stereochemical constraints on their conformation. The models are consistent with the activities of seven active-site mutants. Asp-310 and Glu-432 of type I adenylyl cyclase coordinate a Mg2+ ion. The D310S and D310A mutants have 10-fold reduced Vmax and altered [Mg2+] dependence. The NTP purine moieties bind in mostly hydrophobic pockets. Specificity is conferred by a Lys and an Asp in adenylyl cyclase, and a Glu, an Arg, and a Cys in guanylyl cyclase. The models predict that an Asp from one domain is a general base in the reaction, and that the transition state is stabilized by a conserved Asn-Arg pair on the other domain.

Digging for dead genes: an analysis of the characteristics of the pseudogene population in the Caenorhabditis elegans genome

Pseudogenes are non-functioning copies of genes in genomic DNA, which may either result from reverse transcription from an mRNA transcript (processed pseudogenes) or from gene duplication and subsequent disablement (non-processed pseudogenes). As pseudogenes are apparently ‘dead’, they usually have a variety of obvious disablements (e.g., insertions, deletions, frameshifts and truncations) relative to their functioning homologs. We have derived an initial estimate of the size, distribution and characteristics of the pseudogene population in the Caenorhabditis elegans genome, performing a survey in ‘molecular archaeology’. Corresponding to the 18 576 annotated proteins in the worm (i.e., in Wormpep18), we have found an estimated total of 2168 pseudogenes, about one for every eight genes. Few of these appear to be processed. Details of our pseudogene assignments are available from http://bioinfo.mbb.yale.edu/genome/worm/pseudogene. The population of pseudogenes differs significantly from that of genes in a number of respects: (i) pseudogenes are distributed unevenly across the genome relative to genes, with a disproportionate number on chromosome IV; (ii) the density of pseudogenes is higher on the arms of the chromosomes; (iii) the amino acid composition of pseudogenes is midway between that of genes and (translations of) random intergenic DNA, with enrichment of Phe, Ile, Leu and Lys, and depletion of Asp, Ala, Glu and Gly relative to the worm proteome; and (iv) the most common protein folds and families differ somewhat between genes and pseudogenes—whereas the most common fold found in the worm proteome is the immunoglobulin fold and the most common ‘pseudofold’ is the C-type lectin. In addition, the size of a gene family bears little overall relationship to the size of its corresponding pseudogene complement, indicating a highly dynamic genome. There are in fact a number of families associated with large populations of pseudogenes. For example, one family of seven-transmembrane receptors (represented by gene B0334.7) has one pseudogene for every four genes, and another uncharacterized family (represented by gene B0403.1) is approximately two-thirds pseudogenic. Furthermore, over a hundred apparent pseudogenic fragments do not have any obvious homologs in the worm.

On the role of the K-proton transfer pathway in cytochrome c oxidase

Cytochrome c oxidase is a membrane-bound enzyme that catalyzes the four-electron reduction of oxygen to water. This highly exergonic reaction drives proton pumping across the membrane. One of the key questions associated with the function of cytochrome c oxidase is how the transfer of electrons and protons is coupled and how proton transfer is controlled by the enzyme. In this study we focus on the function of one of the proton transfer pathways of the R. sphaeroides enzyme, the so-called K-proton transfer pathway (containing a highly conserved Lys(I-362) residue), leading from the protein surface to the catalytic site. We have investigated the kinetics of the reaction of the reduced enzyme with oxygen in mutants of the enzyme in which a residue [Ser(I-299)] near the entry point of the pathway was modified with the use of site-directed mutagenesis. The results show that during the initial steps of oxygen reduction, electron transfer to the catalytic site (to form the “peroxy” state, Pr) requires charge compensation through the proton pathway, but no proton uptake from the bulk solution. The charge compensation is proposed to involve a movement of the K(I-362) side chain toward the binuclear center. Thus, in contrast to what has been assumed previously, the results indicate that the K-pathway is used during oxygen reduction and that K(I-362) is charged at pH ≈ 7.5. The movement of the Lys is proposed to regulate proton transfer by “shutting off” the protonic connectivity through the K-pathway after initiation of the O2 reduction chemistry. This “shutoff” prevents a short-circuit of the proton-pumping machinery of the enzyme during the subsequent reaction steps.

Copper-catalyzed oxidation of the recombinant SHa(29–231) prion protein

Metal-catalyzed oxidation may result in structural damage to proteins and has been implicated in aging and disease, including neurological disorders such as Alzheimer's disease and amyotrophic lateral sclerosis. The selective modification of specific amino acid residues with high metal ion affinity leads to subtle structural changes that are not easy to detect but may have dramatic consequences on physical and functional properties of the oxidized protein molecules. PrP contains a histidine-rich octarepeat domain that binds copper. Because copper-binding histidine residues are particularly prone to metal-catalyzed oxidation, we investigated the effect of this reaction on the recombinant prion protein SHaPrP(29–231). Using Cu2+/ascorbate, we oxidized SHaPrP(29–231) in vitro. Oxidation was demonstrated by liquid chromatography/mass spectrometry, which showed the appearance of protein species of higher mass, including increases in multiples of 16, characteristic of oxygen incorporation. Digestion studies using Lys C indicate that the 29–101 region, which includes the histidine-containing octarepeats, is particularly affected by oxidation. Oxidation was time- and copper concentration-dependent and was evident with copper concentrations as low as 1 μM. Concomitant with oxidation, SHaPrP(29–231) suffered aggregation and precipitation, which was nearly complete after 15 min, when the prion protein was incubated at 37°C with a 6-fold molar excess of Cu2+. These findings indicate that PrP, a copper-binding protein, may be particularly susceptible to metal-catalyzed oxidation and that oxidation triggers an extensive structural transition leading to aggregation.