Characterization of MADS homeotic genes in the fern Ceratopteris richardii

The MADS genes encode a family of transcription factors, some of which control the identities of floral organs in flowering plants. To understand the role of MADS genes in the evolution of floral organs, five MADS genes (CMADS1, 2, 3, 4, and 6) were cloned from the fern Ceratopteris richardii, a nonflowering plant. A gene tree of partial amino acid sequences of seed plant and fern MADS genes showed that the fern genes form three subfamilies. All members of one of the fern MADS subfamilies have additional amino-terminal amino acids, which is a synapomorphic character of the AGAMOUS subfamily of the flowering plant MADS genes. Their structural similarity indicates a sister relationship between the two subfamilies. The temporal and spatial patterns of expression of the five fern MADS genes were assessed by Northern blot analyses and in situ hybridizations. CMADS1, 2, 3, and 4 are expressed similarly in the meristematic regions and primordia of sporophyte shoots and roots, as well as in reproductive structures, including sporophylls and sporangial initials, although the amount of expression in each tissue is different in each gene. CMADS6 is expressed in gametophytic tissues but not in sporophytic tissues. The lack of organ-specific expression of MADS genes in the reproductive structures of the fern sporophyte may indicate that the restriction of MADS gene expression to specific reproductive organs and the specialization of MADS gene functions as homeotic selector genes in the flowering plant lineage were important in floral organ evolution.

Crystal structures of eukaryotic translation initiation factor 5A from Methanococcus jannaschii at 1.8 Å resolution

Eukaryotic translation initiation factor 5A (eIF-5A) is a ubiquitous protein found in all eukaryotic cells. The protein is closely associated with cell proliferation in the G1–S stage of the cell cycle. Recent findings show that the eIF-5A proteins are highly expressed in tumor cells and act as a cofactor of the Rev protein in HIV-1-infected cells. The mature eIF is the only protein known to have the unusual amino acid hypusine, a post-translationally modified lysine. The crystal structure of eIF-5A from Methanococcus jannaschii (MJ eIF-5A) has been determined at 1.9 Å and 1.8 Å resolution in two crystal forms by using the multiple isomorphous replacement method and the multiwavelength anomalous diffraction method for the first crystal form and the molecular replacement method for the second crystal form. The structure consists of two folding domains, one of which is similar to the oligonucleotide-binding domain found in the prokaryotic cold shock protein and the translation initiation factor IF1 despite the absence of any significant sequence similarities. The 12 highly conserved amino acid residues found among eIF-5As include the hypusine site and form a long protruding loop at one end of the elongated molecule.

Structure-based assignment of the biochemical function of a hypothetical protein: A test case of structural genomics

Many small bacterial, archaebacterial, and eukaryotic genomes have been sequenced, and the larger eukaryotic genomes are predicted to be completely sequenced within the next decade. In all genomes sequenced to date, a large portion of these organisms’ predicted protein coding regions encode polypeptides of unknown biochemical, biophysical, and/or cellular functions. Three-dimensional structures of these proteins may suggest biochemical or biophysical functions. Here we report the crystal structure of one such protein, MJ0577, from a hyperthermophile, Methanococcus jannaschii, at 1.7-Å resolution. The structure contains a bound ATP, suggesting MJ0577 is an ATPase or an ATP-mediated molecular switch, which we confirm by biochemical experiments. Furthermore, the structure reveals different ATP binding motifs that are shared among many homologous hypothetical proteins in this family. This result indicates that structure-based assignment of molecular function is a viable approach for the large-scale biochemical assignment of proteins and for discovering new motifs, a basic premise of structural genomics.

Structural code for DNA recognition revealed in crystal structures of papillomavirus E2-DNA targets

Transcriptional regulation in papillomaviruses depends on sequence-specific binding of the regulatory protein E2 to several sites in the viral genome. Crystal structures of bovine papillomavirus E2 DNA targets reveal a conformational variant of B-DNA characterized by a roll-induced writhe and helical repeat of 10.5 bp per turn. A comparison between the free and the protein-bound DNA demonstrates that the intrinsic structure of the DNA regions contacted directly by the protein and the deformability of the DNA region that is not contacted by the protein are critical for sequence-specific protein/DNA recognition and hence for gene-regulatory signals in the viral system. We show that the selection of dinucleotide or longer segments with appropriate conformational characteristics, when positioned at correct intervals along the DNA helix, can constitute a structural code for DNA recognition by regulatory proteins. This structural code facilitates the formation of a complementary protein–DNA interface that can be further specified by hydrogen bonds and nonpolar interactions between the protein amino acids and the DNA bases.

Cloning and functional expression of a cDNA encoding a pheromone gland-specific acyl-CoA Δ11-desaturase of the cabbage looper moth, Trichoplusia ni

Desaturation of coenzyme-A esters of saturated fatty acids is a common feature of sex pheromone biosynthetic pathways in the Lepidoptera. The enzymes that catalyze this step share several biochemical properties with the ubiquitous acyl-CoA Δ9-desaturases of animals and fungi, suggesting a common ancestral origin. Unlike metabolic acyl-CoA Δ9-desaturases, pheromone desaturases have evolved unusual regio- and stereoselective activities that contribute to the remarkable diversity of chemical structures used as pheromones in this large taxonomic group. In this report, we describe the isolation of a cDNA encoding a pheromone gland desaturase from the cabbage looper moth, Trichoplusia ni, a species in which all unsaturated pheromone products are produced via a Δ11Z-desaturation mechanism. The largest ORF of the ≈1,250-bp cDNA encodes a 349-aa apoprotein (PDesat-Tn Δ11Z) with a predicted molecular mass of 40,240 Da. Its hydrophobicity profile is similar overall to those of rat and yeast Δ9-desaturases, suggesting conserved transmembrane topology. A 182-aa core domain delimited by conserved histidine-rich motifs implicated in iron-binding and catalysis has 72 and 58% similarity (including conservative substitutions) to acyl-CoA Δ9Z-desaturases of rat and yeast, respectively. Northern blot analysis revealed an ≈1,250-nt PDesat-Tn Δ11Z mRNA that is consistent with the spatial and temporal distribution of Δ11-desaturase enzyme activity. Genetic transformation of a desaturase-deficient strain of the yeast Saccharomyces cerevisiae with an expression plasmid encoding PDesat-Tn Δ11Z resulted in complementation of the strain’s fatty acid auxotrophy and the production of Δ11Z-unsaturated fatty acids.

Characterization of residual structure in the thermally denatured state of barnase by simulation and experiment: Description of the folding pathway

Residual structure in the denatured state of a protein may contain clues about the early events in folding. We have simulated by molecular dynamics the denatured state of barnase, which has been studied by NMR spectroscopy. An ensemble of 104 structures was generated after 2 ns of unfolding and following for a further 2 ns. The ensemble was heterogeneous, but there was nonrandom, residual structure with persistent interactions. Helical structure in the C-terminal portion of helix α1 (residues 13–17) and in helix α2 as well as a turn and nonnative hydrophobic clustering between β3 and β4 were observed, consistent with NMR data. In addition, there were tertiary contacts between residues in α1 and the C-terminal portion of the β-sheet. The simulated structures allow the rudimentary NMR data to be fleshed out. The consistency between simulation and experiment inspires confidence in the methods. A description of the folding pathway of barnase from the denatured to the native state can be constructed by combining the simulation with experimental data from φ value analysis and NMR.

Computer-based analysis of the binding steps in protein complex formation

Computer models were used to examine whether and under what conditions the multimeric protein complex is inhibited by high concentrations of one of its components—an effect analogous to the prozone phenomenon in precipitin tests. A series of idealized simple “ball-and-stick” structures representing small oligomeric complexes of protein molecules formed by reversible binding reactions were analyzed to determine the binding steps leading to each structure. The equilibrium state of each system was then determined over a range of starting concentrations and Kds and the steady-state concentration of structurally complete oligomer calculated for each situation. A strong inhibitory effect at high concentrations was shown by any protein molecule forming a bridge between two or more separable parts of the complex. By contrast, proteins linked to the outside of the complex by a single bond showed no inhibition whatsoever at any concentration. Nonbridging, multivalent proteins in the body of the complex could show an inhibitory effect or not depending on the structure of the complex and the strength of its bonds. On the basis of this study, we suggest that the prozone phenomenon will occur widely in living cells and that it could be a crucial factor in the regulation of protein complex formation.

Structural features of the kringle domain determine the intracellular degradation of under-γ-carboxylated prothrombin: Studies of chimeric rat/human prothrombin

Vitamin K antagonists such as warfarin inhibit the vitamin K-dependent γ-glutamyl carboxylation during protein processing and block the secretion of under-γ-carboxylated prothrombin (FII) in the rat but not in the human or bovine. Under-γ-carboxylated prothrombin is also secreted from warfarin-treated human (HepG2) cell cultures but is degraded in the endoplasmic reticulum in warfarin-treated rat (H-35) cell cultures. This differential response to warfarin has been shown to be determined by the structural difference in the proteins rather than by the origin of the cell line. When recombinant rat prothrombin (rFII) and human prothrombin (hFII) were expressed in a transformed human kidney cell line (HEK293), secretion of rFII but not hFII was drastically decreased in response to warfarin. To determine the structural signal required for this differential response, chimeric cDNAs with the propeptide/Gla domains, kringle domain, and serine protease domain exchanged between rFII and hFII were generated (FIIRHH and FIIHRR, FIIRRH and FIIHHR, FIIRHR and FIIHRH) and expressed in both warfarin-treated HEK293 cells and HepG2 cells. The presence of the hFII kringle domain changed the stability of rFII to that of hFII, and the rFII kringle domain changed the stability of hFII to that of rFII. The kringle domain therefore is critical in determining the metabolic fate of under-γ-carboxylated prothrombin precursors during processing. Prothrombin contains two kringle structures, and expression of additional rFII/hFII chimeras (FIIHrhH and FIIHhrH, FIIRrhR, and FIIRhrR) was used to determine that the first of the two kringles plays a more important role in the recognition process.

Crystal structure of the Holliday junction migration motor protein RuvB from Thermus thermophilus HB8

We report here the crystal structure of the RuvB motor protein from Thermus thermophilus HB8, which drives branch migration of the Holliday junction during homologous recombination. RuvB has a crescent-like architecture consisting of three consecutive domains, the first two of which are involved in ATP binding and hydrolysis. DNA is likely to interact with a large basic cleft, which encompasses the ATP-binding pocket and domain boundaries, whereas the junction-recognition protein RuvA may bind a flexible β-hairpin protruding from the N-terminal domain. The structures of two subunits, related by a noncrystallographic pseudo-2-fold axis, imply that conformational changes of motor protein coupled with ATP hydrolysis may reflect motility essential for its translocation around double-stranded DNA.

Mutagenesis of the regulatory domain of phenylalanine hydroxylase

The regulatory domain of phenylalanine hydroxylase (PAH, EC 1.14.16.1) consists of more than 100 amino acids at the N terminus, the removal of which significantly activates the enzyme. To study the regulatory properties controlled by the N terminus, a series of truncations and site-specific mutations were made in this region of rat PAH. These enzymes were expressed highly in Escherichia coli and purified through a pterin-conjugated Sepharose affinity column. The removal of the first 26 amino acids of the N terminus increased the activity by about 20-fold, but removal of the first 15 amino acids increased the activity by only 2-fold. Replacing serine-29 of rat PAH with cysteine from the same site of human PAH increased the activity by more than 4-fold. Mutation of serine to other amino acids with varying side chains: alanine, methionine, leucine, aspartic acid, asparagine, and arginine also resulted in significant activation, indicating a serine-specific inhibitory effect. But these site-specific mutants showed 30–40% lower activity when assayed with 6-methyl-5,6,7,8-tetrahydropterin. Stimulation of hydroxylase activity by preincubation of the enzyme with phenylalanine was inversely proportional to the activation state of all these mutants. Combined with recent crystal structures of PAH [Kobe, B. et al. (1999) Nat. Struct. Biol. 6, 442–448; and Erlandsen, H., Bjorgo, E., Flatmark, T. & Stevens, R. C. (2000) Biochemistry 39, 2208–2217], these data suggest that residues 16–26 have a controlling regulatory effect on the activity by interaction with the dihydroxypropyl side chain of (6R)-5,6,7,8-tetrahydrobiopterin. The serine/cysteine switch explains the difference in regulatory properties between human and rat PAH. The N terminus as a whole is important for maintaining rat PAH in an optimum catalytic conformation.

Rapid evolution of the DNA-binding site in LAGLIDADG homing endonucleases

Sequence analysis of chloroplast and mitochondrial large subunit rRNA genes from over 75 green algae disclosed 28 new group I intron-encoded proteins carrying a single LAGLIDADG motif. These putative homing endonucleases form four subfamilies of homologous enzymes, with the members of each subfamily being encoded by introns sharing the same insertion site. We showed that four divergent endonucleases from the I-CreI subfamily cleave the same DNA substrates. Mapping of the 66 amino acids that are conserved among the members of this subfamily on the 3-dimensional structure of I-CreI bound to its recognition sequence revealed that these residues participate in protein folding, homodimerization, DNA recognition and catalysis. Surprisingly, only seven of the 21 I-CreI amino acids interacting with DNA are conserved, suggesting that I-CreI and its homologs use different subsets of residues to recognize the same DNA sequence. Our sequence comparison of all 45 single-LAGLIDADG proteins identified so far suggests that these proteins share related structures and that there is a weak pressure in each subfamily to maintain identical protein–DNA contacts. The high sequence variability we observed in the DNA-binding site of homologous LAGLIDADG endonucleases provides insight into how these proteins evolve new DNA specificity.

Identification and properties of the crenarchaeal single-stranded DNA binding protein from Sulfolobus solfataricus

Single-stranded DNA binding proteins (SSBs) play central roles in cellular and viral processes involving the generation of single-stranded DNA. These include DNA replication, homologous recombination and DNA repair pathways. SSBs bind DNA using four ‘OB-fold’ (oligonucleotide/oligosaccharide binding fold) domains that can be organised in a variety of overall quaternary structures. Thus eubacterial SSBs are homotetrameric whilst the eucaryal RPA protein is a heterotrimer and euryarchaeal proteins vary significantly in their subunit compositions. We demonstrate that the crenarchaeal SSB protein is an abundant protein with a unique structural organisation, existing as a monomer in solution and multimerising on DNA binding. The protein binds single-stranded DNA distributively with a binding site size of ~5 nt per monomer. Sulfolobus SSB lacks the zinc finger motif found in the eucaryal and euryarchaeal proteins, possessing instead a flexible C-terminal tail, sensitive to trypsin digestion, that is not required for DNA binding. In comparison with Escherichia coli SSB, the tail may play a role in protein–protein interactions during DNA replication and repair.

Unwinding of nucleic acids by HCV NS3 helicase is sensitive to the structure of the duplex

Hepatitis C virus (HCV) helicase, non-structural protein 3 (NS3), is proposed to aid in HCV genome replication and is considered a target for inhibition of HCV. In order to investigate the substrate requirements for nucleic acid unwinding by NS3, substrates were prepared by annealing a 30mer oligonucleotide to a 15mer. The resulting 15 bp duplex contained a single-stranded DNA overhang of 15 nt referred to as the bound strand. Other substrates were prepared in which the 15mer DNA was replaced by a strand of peptide nucleic acid (PNA). The PNA–DNA substrate was unwound by NS3, but the observed rate of strand separation was at least 25-fold slower than for the equivalent DNA–DNA substrate. Binding of NS3 to the PNA–DNA substrate was similar to the DNA–DNA substrate, due to the fact that NS3 initially binds to the single-stranded overhang, which was identical in each substrate. A PNA–RNA substrate was not unwound by NS3 under similar conditions. In contrast, morpholino–DNA and phosphorothioate–DNA substrates were utilized as efficiently by NS3 as DNA–DNA substrates. These results indicate that the PNA–DNA and PNA–RNA heteroduplexes adopt structures that are unfavorable for unwinding by NS3, suggesting that the unwinding activity of NS3 is sensitive to the structure of the duplex.

Conserved stem–loop structures in the HIV-1 RNA region containing the A3 3′ splice site and its cis-regulatory element: possible involvement in RNA splicing

The HIV-1 transcript is alternatively spliced to over 30 different mRNAs. Whether RNA secondary structure can influence HIV-1 RNA alternative splicing has not previously been examined. Here we have determined the secondary structure of the HIV-1/BRU RNA segment, containing the alternative A3, A4a, A4b, A4c and A5 3′ splice sites. Site A3, required for tat mRNA production, is contained in the terminal loop of a stem–loop structure (SLS2), which is highly conserved in HIV-1 and related SIVcpz strains. The exon splicing silencer (ESS2) acting on site A3 is located in a long irregular stem–loop structure (SLS3). Two SLS3 domains were protected by nuclear components under splicing condition assays. One contains the A4c branch points and a putative SR protein binding site. The other one is adjacent to ESS2. Unexpectedly, only the 3′ A residue of ESS2 was protected. The suboptimal A3 polypyrimidine tract (PPT) is base paired. Using site-directed mutagenesis and transfection of a mini-HIV-1 cDNA into HeLa cells, we found that, in a wild-type PPT context, a mutation of the A3 downstream sequence that reinforced SLS2 stability decreased site A3 utilization. This was not the case with an optimized PPT. Hence, sequence and secondary structure of the PPT may cooperate in limiting site A3 utilization.

Detection of alkali metal ions in DNA crystals using state-of-the-art X-ray diffraction experiments

The observation of light metal ions in nucleic acids crystals is generally a fortuitous event. Sodium ions in particular are notoriously difficult to detect because their X-ray scattering contributions are virtually identical to those of water and Na+…O distances are only slightly shorter than strong hydrogen bonds between well-ordered water molecules. We demonstrate here that replacement of Na+ by K+, Rb+ or Cs+ and precise measurements of anomalous differences in intensities provide a particularly sensitive method for detecting alkali metal ion-binding sites in nucleic acid crystals. Not only can alkali metal ions be readily located in such structures, but the presence of Rb+ or Cs+ also allows structure determination by the single wavelength anomalous diffraction technique. Besides allowing identification of high occupancy binding sites, the combination of high resolution and anomalous diffraction data established here can also pinpoint binding sites that feature only partial occupancy. Conversely, high resolution of the data alone does not necessarily allow differentiation between water and partially ordered metal ions, as demonstrated with the crystal structure of a DNA duplex determined to a resolution of 0.6 Å.