Structural repeats and evolution of tobacco mosaic virus coat protein and RNA

The three-dimensional structure of the tobacco mosaic virus (TMV) coat protein disk suggests a possible pathway for the early evolution of the virus self-assembly mechanism.The coat protein contains a 2-fold repeated structural pattern in the folding of both its four alpha helices (A,B,C,D), which run alternately forward and back along the radius of the disk, and the four-stranded antiparallel pleated sheet which links these helices to the hydrophobic girdle at the outer rim of the disk. Helices A and B can be approximately superposed on C and D by a screw rotation about a molecular pseudo-dyad axis which lies nearly parallel to the plane of the protein disk. This operation relates 29 pairs of α-carbon positions with a root-mean-square deviation of 1.77 Å. A second pseudo-dyad in the pleated-sheet region relates 14 more atom pairs with a deviation of 2.32 Å and forms a distorted continuation of the relationship between the helices. The helix dyad also relates repeated pairs of functionally important amino acids which take part in intersubunit contacts.We have analysed these structural repeats and tested their significance by comparing them with repeats in other “helix quartet” proteins, cytochrome b⁵ and the hemerythrins, as well as with an irregular helix cluster in thermolysin. TMV is noticeably more repetitive than the others, including hemerythrin which is thought to have evolved by gene duplication.We propose that the primitive TMV coat protein was a dimeric structure of two smaller units paired about a 2-fold axis. Each unit was a pair of helices, linked at the inner radius of the virus rod by a short bend, where the RNA binding site formed, and connected at the outer radius by two short strands of beta sheet. A tandem gene duplication joined the two units and formed the present helix quartet. The flexible loop which now runs into the centre of the virus and connects helix C to helix D developed later. The assembly origin RNA may have evolved from part of the coat protein RNA which codes for this loop.     

Predominant end-products of prophage Mu DNA transposition during the lytic cycle are replicon fusions

The structure of l-arabinose-binding protein (M_r 33, 100), an essential component of the osmotic shock-sensitive, high-affinity l-arabinose transport system in Escherichia coli, has been determined at 2.4 Å resolution. The phases were solved by the method of multiple isomorphous replacement, using four derivatives, p-chloromercuribenzenesulfonate and CdI₂ (data to 2.4 Å resolution), and p-chloromercurinitrophenol and (NH₄)₂PtCl₄ to 3.5 Å resolution. A final mean figure of merit of 0.65 was obtained for 9628 reflections.With the aid of the amino acid sequence determined by Hogg &; Hermodson (1977), a complete model of the protein molecule has been determined using initially an optical comparator. The entire model was subsequently examined in detail using a computer graphic system.The protein molecule is ellipsoidal (axial ratio of 2:1), and consists of two globular domains (designated P and Q). Each domain is made from two separate polypeptide chain segments. Despite the discontinuity in the folding, the arrangements of the secondary structure in the two domains are very similar. Both domains contain a six-stranded parallel β-sheet (with the exception of the sixth anti-parallel strand in the Q domain) flanked by two α-helices on either side. The packing topology is α/β. A C-terminal helix is shared by both domains.The two domains show significant conformational similarity but lack sequence homology. A comparison of the two domains revealed that of the 139 α-carbons in the P domain and 152 in the Q domain, 92 were found to be equivalent with a root-mean-square distance of 2.6 Å.The cleft formed by the packing of the two domains is predominantly lined with hydrophilic residues. The sugar-binding site is located in this cleft.     

Genomic organization of the human multidrug resistance (MDR1) gene and origin of P-glycoproteins

The MDR1 gene, responsible for multidrug resistance in human cells, encodes a broad specificity efflux pump (P-glycoprotein). P-glycoprotein consists of two similar halves, each half including a hydrophobic transmembrane region and a nucleotide-binding domain. On the basis of sequence homology between the N-terminal and C-terminal halves of P-glycoprotein, we have previously suggested that this gene arose by duplication of a primordial gene. We have now determined the complete intron/exon structure of the MDR1 gene by direct sequencing of cosmid clones and enzymatic amplification of genomic DNA segments. The MDR1 gene includes 28 introns, 26 of which interrupt the protein-coding sequence. Although both halves of the protein-coding sequence are composed of approximately the same number of exons, only two intron pairs, both within the nucleotide-binding domains, are located at conserved positions in the two halves of the protein. The other introns occur at different locations in the two halves of the protein and in most cases interrupt the coding sequence at different positions relative to the open reading frame. These results suggest that the P-glycoprotein arose by fusion of genes for two related but independently evolved proteins rather than by internal duplication.     

Human hexokinase: sequences of amino- and carboxyl-terminal halves are homologous

cDNA clones encoding human hexokinase have been isolated from an adult kidney library. Analysis of this 917 amino acid protein (Mr = 102,519) indicates that the sequences of the NH2- and COOH-terminal halves, corresponding to the regulatory and catalytic domains, respectively, are homologous; and that eukaryotic hexokinases evolved by duplication of a gene encoding a protein of 450 amino acids. The COOH-terminal half of the protein created by this gene duplication retained the glucose binding site and glucose phosphorylating activity while the substrate binding sites of the NH2-terminal half evolved into a new allosteric effector site.     

Three-fold structural pattern in the soybean trypsin inhibitor (Kunitz)

The molecule contains three very similar irregular Y-shaped lobes of antiparallel twisted β-sheet, which are grouped symmetrically round a central axis and linked by hydrogen bonds to form a six-stranded barrel. Each lobe can be superposed on either neighbour by a rotation of approximately 120 °. Of the 160 residues seen in the X-ray electron density map, 101 may be superposed onto other residues within a root-mean-square distance of 2.1 Å. The bond which reacts with trypsin lies on a loop between the first two lobes. It is suggested that the protein evolved from a primitive symmetrical trimer of identical subunits by tandem gene triplication.     

Evolution of the ATP-binding-cassette transmembrane transporters of vertebrates

The ATP-binding-cassette transmembrane transporters (ABC transporters) known from vertebrates belong to four major subfamilies: (1) the P- glycoproteins (Pgp); (2) the cystic fibrosis transmembrane conductance regulators (CFTR); (3) the Tap proteins encoded with the major histocompatibility complex of mammals; and (4) the peroxisomal membrane proteins. Both Pgp and CFTR have a structure suggesting a past internal gene duplication; a phylogenetic analysis indicated that these duplications occurred independently, while an independent tandem gene duplication occurred in the case of the Tap family. Both the Pgp and Tap proteins show evidence of relationship to bacterial ABC transporters lacking internal duplication, and both are significantly more closely related to the HlyB and MsbA families of transporters from purple bacteria than they are to ABC transporters from nonpurple bacteria. The simplest hypothesis to explain this observation is that eukaryotic Pgp and Tap genes are descended from a mitochondrial gene or genes that were subsequently translocated to the nuclear genome. The Pgp genes of eukaryotes are characterized by a remarkable degree of convergent evolution between the ATP-binding cassettes of their N- terminal and C-terminal halves, whereas no such convergence is seen between the two halves of CFTR genes or between the duplicated Tap genes. Exon 13 of the CFTR gene, which encodes a putative regulatory domain not found in other ABC transporters apart from CFTR, showed high levels of both synonymous and nonsynonymous difference in comparisons among different mammalian species, suggesting that this region is a mutational hot spot.      

Gene duplication,genome duplication,and the functional diversification of vertebrate globins

The functional diversification of the vertebrate globin gene superfamily provides an especially vivid illustration of the role of gene duplication and whole-genome duplication in promoting evolutionary innovation. For example, key globin proteins that evolved specialized functions in various aspects of oxidative metabolism and oxygen signaling pathways (hemoglobin [Hb], myoglobin [Mb], and cytoglobin [Cygb]) trace their origins to two whole-genome duplication events in the stem lineage of vertebrates. The retention of the proto-Hb and Mb genes in the ancestor of jawed vertebrates permitted a physiological division of labor between the oxygen-carrier function of Hb and the oxygen-storage function of Mb. In the Hb gene lineage, a subsequent tandem gene duplication gave rise to the proto α- and β-globin genes, which permitted the formation of multimeric Hbs composed of unlike subunits (α₂β₂). The evolution of this heteromeric quaternary structure was central to the emergence of Hb as a specialized oxygen-transport protein because it provided a mechanism for cooperative oxygen-binding and allosteric regulatory control. Subsequent rounds of duplication and divergence have produced diverse repertoires of α- and β-like globin genes that are ontogenetically regulated such that functionally distinct Hb isoforms are expressed during different stages of prenatal development and postnatal life. In the ancestor of jawless fishes, the proto Mb and Hb genes appear to have been secondarily lost, and the Cygb homolog evolved a specialized respiratory function in blood-oxygen transport. Phylogenetic and comparative genomic analyses of the vertebrate globin gene superfamily have revealed numerous instances in which paralogous globins have convergently evolved similar expression patterns and/or similar functional specializations in different organismal lineages.     

The mechanism of vault opening from the high resolution structure of the N‐terminal repeats of MVP

Vaults are ubiquitous ribonucleoprotein complexes involved in a diversity of cellular processes, including multidrug resistance, transport mechanisms and signal transmission. The vault particle shows a barrel‐shaped structure organized in two identical moieties, each consisting of 39 copies of the major vault protein MVP. Earlier data indicated that vault halves can dissociate at acidic pH. The crystal structure of the vault particle solved at 8 Å resolution, together with the 2.1‐Å structure of the seven N‐terminal domains (R1–R7) of MVP, reveal the interactions governing vault association and provide an explanation for a reversible dissociation induced by low pH. The structural comparison with the recently published 3.5 Å model shows major discrepancies, both in the main chain tracing and in the side chain assignment of the two terminal domains R1 and R2.     

α/β Barrel evolution and the modular assembly of enzymes: Emerging trends in the flavin oxidase/dehydrogenase family

α/β barrels have an ill-defined origin. Evidence exists which favours their divergent evolution from a common ancestral barrel and convergent evolution to a stable fold. However, recent sequence and structural information for the flavin oxidase/dehydrogenase family of barrel enzymes indicate that sub-families of α/β barrels have evolved divergently. The modular fusion of barrel domains with core structures from other gene families has also contributed to the evolution of related but catalytically distinct enzyme molecules within each sub-family of the flavin oxidases/dehydrogenases. An analysis of the structures and sequences of the flavin oxidases/dehydrogenases has now enabled studies focusing on the evolutionary origins and modular assembly of this important family of proteins to be initiated.     

Comparison of the folding of 2-Keto-3-deoxy-6-phosphogluconate aldolase,triosephosphate isomerase and pyruvate kinase: Implications in molecular evolution

The 8-fold α/β barrel conformation of 2-keto-3-deoxy-6-phosphogluconate aldolase has been compared to that of triosephosphate isomerase and the A-domain of pyruvate kinase. There are eight supersecondary structure units (α/β) in each of these proteins, and the comparisons were carried out in orientations corresponding to each of the possible congruences, i.e. first to first, first to second,… of the supersecondary structure units. The comparison of the C_α structure of the main chain folding of the three enzymes indicated about 150 equivalences with rootmean-square differences of about 3.1 Å, with no orientational preference, including the aldolase with itself. In addition, there is no sequence homology between the aldolase and the isomerase, and no indication of gene duplication in the former. The lack of orientational preference among the three enzymes suggests convergence to a fold of exceptional stability. However, all three enzymes activate a CH bond adjacent to a carbonyl, and their active sites correspond to the f strand, F helix region of the α/β barrel, thus contradicting the foregoing and suggesting divergent evolution from a common precursor. Other and similar arguments are also presented for and against convergent evolution of these three strikingly similar enzymes.     

Concerted evolution of exons and introns in the MHC-linked tenascin-X gene of mammals.

Two modes of evolution of repeated domains in proteins have been described: (1) a conservative mode, whereby individual domains are conserved across gene duplication and speciation events, and (2) a concerted mode, whereby repeat domains become homogenized within a gene, presumably by intragenic partial duplication and/or gene conversion. The evolution of repeated EGF-like and fibronection-type-III-like (Fn-III) domains in the vertebrate extracellular matrix proteins tenascin-X (TNX) and tenascin-C (TNC) was studied by comparisons between human and mouse orthologs and between the paralogous TNC and TNX genes. The EGF-like repeats have largely been homogenized within each gene by concerted evolution since the duplication of the two genes but have been conserved since the divergence of rodents and primates. The Fn-III domains of TNC have likewise mainly evolved in a conservative fashion since the divergence of rodents and primates. In contrast, the Fn-III repeats of TNX fall into three distinct categories with regard to mode of evolution: (1) The three C-terminal repeats have been conserved since before duplication of the TNX and TNC genes. (2) Certain other repeats have been homogenized within each gene since gene duplication but have been conserved since the divergence of rodents and primates. (3) Still other repeats have evolved in a concerted fashion in rodent and primate lineages since their divergence. Remarkably, certain introns adjacent to the exons encoding these concertedly evolving Fn-III repeats have themselves evolved in a concerted fashion. This is the first known example of concerted evolution of repeated introns within a protein-coding gene.     

Structure of trichosanthin at 1.88 Å resolution

Trichosanthin (TCS) is one of the single chain ribosome-inactivating proteins (RIPs). The crystals of the orthorhombic form of trichosanthin have been obtained from a citrate buffer (pH 5.4) with KC1 as the precipitant. The crystal belongs to the space group P2₁2₁2₁ with a = 38.31, b = 76.22, c = 79.21 Å. The structure was solved by molecular replacement method and refined using the programs XPLOR and PROLSQ to an R-factor of 0.191 for the reflections within the 6–1.88 Å resolution range. The bond length and bond angle in the protein molecule have root-mean-square deviations from ideal value of 0.013 Å and 3.3°, respectively. The refined model includes 247 residues and 197 water molecules. The TCS molecule consists of two structural domains. The large domain contains six α-helices, a six stranded sheet, and an antiparallel β-sheet. The small domain has a largest α-helix, which shows a distinct bend. The possible active site of the molecule located on the cleft between two domains was proposed. In the active site Arg-163 and Glu-160, Glu-189 and Arg-122 form two ion pairs, Glu-189 and Gln-156 are hydrogen bonded to each other. Three water molecules are bonded to the residues in the active site region. The structures of TCS molecule and ricin A-chain (RTA) superimpose quite well, showing that the structures of the two protein molecules are homologous. Comparison of the structures of the TCS molecule in this orthorhombic crystal with that in the monoclinic crystal indicates that there are no essential differences of the structures between the two protein crystals. © 1994 Wiley-Liss, Inc.     

Gene duplication in the evolution of the two complementing domains of gram-negative bacterial tetracycline efflux proteins

The resistance of Gram- bacteria to the broad-spectrum antibiotic tetracycline (Tc) results from energy-dependent drug efflux mediated by the tet gene product, the cytoplasmic membrane Tet protein. Amino acid (aa) sequences deduced from total tet nucleotide sequences of three different resistance determinants (classes A, B and C) indicate that the protein products [Tet(A), Tet(B), and Tet(C)] share a common ancestor. Hydropathic analysis of Tet sequences predicts twelve transmembrane segments in each protein, with six occurring in each half of the molecule. More importantly, the linear distributions of these segments in the N- and C-terminal halves are nearly identical, suggesting that the two halves of each Tet protein are related by a process of tandem gene duplication and divergence. Indeed, a variable but significant conservation of sequence was detected among the N- and C-terminal halves for all possible comparisons of the three proteins. Such conservation was not observed within other prokaryotic integral membrane proteins or when other prokaryotic proteins were compared to Tet halves. Similarity, both in sequence and in predicted transmembrane structural organization, strongly suggests that a common ancestor of Tet(A), Tet(B), and Tet(C) arose by duplication of a gene reading frame specifying a transmembrane protein of approximately 200 aa residues. The two halves of Tet proteins correspond to the two domains, alpha and beta, which have distinct, complementary roles in Tc efflux. Nevertheless, selective constraints to function in the cytoplasmic membrane have apparently led to maintenance of similar patterns of secondary structural organization in these complementary domains.     

The crystal structure of a TIR domain from Arabidopsis thaliana reveals a conserved helical region unique to plants

Plants use a highly evolved immune system to exhibit defense response against microbial infections. The plant TIR domain, together with the nucleotide‐binding (NB) domain and/or a LRR region, forms a type of molecule, named resistance (R) proteins, that interact with microbial effector proteins and elicit hypersensitive responses against infection. Here, we report the first crystal structure of a plant TIR domain from Arabidopsis thaliana (AtTIR) solved at a resolution of 2.0 Å. The structure consists of five β‐strands forming a parallel β‐sheet at the core of the protein. The β‐strands are connected by a series of α‐helices and the overall fold mimics closely that of other mammalian and bacterial TIR domains. However, the region of the αD‐helix reveals significant differences when compared with other TIR structures, especially the αD3‐helix that corresponds to an insertion only present in plant TIR domains. Available mutagenesis data suggest that several conserved and exposed residues in this region are involved in the plant TIR signaling function.     

tRNA Creation by Hairpin Duplication

Many studies have suggested that the modern cloverleaf structure of tRNA may have arisen through duplication of a primordial hairpin, but the timing of this duplication event has been unclear. Here we measure the level of sequence identity between the two halves of each of a large sample of tRNAs and compare this level to that of chimeric tRNAs constructed either within or between groups defined by phylogeny and/or specificity. We find that actual tRNAs have significantly more matches between the two halves than do random sequences that can form the tRNA structure, but there is no difference in the average level of matching between the two halves of an individual tRNA and the average level of matching between the two halves of the chimeric tRNAs in any of the sets we constructed. These results support the hypothesis that the modern tRNA cloverleaf arose from a single hairpin duplication prior to the divergence of modern tRNA specificities and the three domains of life. [Reviewing Editor: Dr. Niles Lehman]     

NADP-Dependent enzymes. II: Evolution of the mono- and dinucleotide binding domains

Nicotinamide adenine dinucleotides [NAD and NADP with both referred to as NAD(P)] are among the more diffuse redox cofactors. Despite their stereochemical similarity where the only difference is a phosphomonoester on the ribose near the adenine of NADP, they show different biochemical reactivities with NAD behaving as an oxidant and NADP as a reductant. NAD(P)-dependent enzymes generally share a common open α/β fold with few exceptions only recently structurally characterized. This study of the molecular evolution of the NAD(P) binding domains, possible given the large number of known molecular structures, addresses two main questions: 1) can a common fold exist in different biological systems (divergent evolution) and 2) does a relationship exist among similar biological systems that display different folds (convergent evolution)? Both the structures of mono- and dinucleotide binding domains have been classified by cluster analysis based on the similarity evaluated by their main chain Cα superposition. Moreover, the cofactor conformations and the stereochemical characteristics of their pockets have also been classified by analogous methods on the basis of the published tertiary structures. Two primary results appear: 1) the classification of the mononucleotide binding domains is different from that of the dinucleotide binding folds and 2) both divergent and convergent evolutionary pathways can be hypothesized, the latter less frequently observed and less pronounced but nevertheless evident. The generally accepted hypothesis that dinucleotide binding domains have evolved by gene duplication of primordial genes coding for the smaller mononucleotide binding domains is acceptable but the two halves of the resulting dinucleotide binding domains are evolutionarly uncorrelated. The NH₂-terminal mononucleotide binding domain is less variable than the COOH-terminal half, probably because it involves the binding of the ADP moiety of NAD(P) invariant in all examined systems. There is evidence to postulate that evolutionary pathways for NAD(P)-dependent enzymes are both divergent and convergent. In fact, nearly all combinations of similarity/dissimilarity in overall fold, cofactor conformation, and cofactor binding pocket structural characteristics for each enzyme pair examined are possible. The NAD(P)-dependent enzymes apparently provide a canonical example of an evolutionary principle that “anything goes.” © 1997 Wiley-Liss Inc.     

Comparative evolution of photosynthetic genes in response to polyploid and nonpolyploid duplication

The likelihood of duplicate gene retention following polyploidy varies by functional properties (e.g. gene ontologies or protein family domains), but little is known about the effects of whole-genome duplication on gene networks related by a common physiological process. Here, we examined the effects of both polyploid and nonpolyploid duplications on genes encoding the major functional groups of photosynthesis (photosystem I, photosystem II, the light-harvesting complex, and the Calvin cycle) in the cultivated soybean (Glycine max), which has experienced two rounds of whole-genome duplication. Photosystem gene families exhibit retention patterns consistent with dosage sensitivity (preferential retention of polyploid duplicates and elimination of nonpolyploid duplicates), whereas Calvin cycle and light-harvesting complex gene families do not. We observed similar patterns in barrel medic (Medicago truncatula), which shared the older genome duplication with soybean but has evolved independently for approximately 50 million years, and in Arabidopsis (Arabidopsis thaliana), which experienced two nested polyploidy events independent from the legume duplications. In both soybean and Arabidopsis, Calvin cycle gene duplicates exhibit a greater capacity for functional differentiation than do duplicates within the photosystems, which likely explains the greater retention of ancient, nonpolyploid duplicates and larger average gene family size for the Calvin cycle relative to the photosystems.     

The 1.8 Å Resolution Structure of Hevamine,a Plant Chitinase/Lysozyme,and Analysis of the Conserved Sequence and Structure Motifs of Glycosyl Hydrolase Family 18

The three-dimensional structure of hevamine, a plant enzyme with chitinase and lysozyme activity, has been refined at 1.8 Å resolution to an R-factor of 14.9% and a freeR-factor of 19.6%. The final model consists of all 273 amino acid residues and 206 ordered water molecules. Two non-prolinecis-peptides were identified, involving Phe32 and Trp255, both of which are implicated in substrate binding.Other glycosyl hydrolase family 18 proteins with known three-dimen sional structure are bacterial chitinase A, endo-β-N-acetylglucosaminidase F₁, endo-β-N-acetylglucosaminidase H, and the two plant proteins concanavalin B and narbonin, which have no known enzymatic activity. All these structures contain a (βα)₈barrel fold, with the two family 18 consensus regions roughly corresponding to the third and fourth barrel strands. This confirms the grouping of these proteins into family 18, which was only based on weak and local sequence similarity. The substrate specificity of the enzymes is determined by the loops following the barrel strands that form the substrate binding site. All enzymes have an aspartic acid and a glutamic acid residue in positions identical with Asp 125 and the catalytic Glu127 of hevamine. The lack of chitinase activity of concanavalin B and narbonin can be explained by the absence of one of these carboxylate groups, and by differences in the loops that form the substrate-binding cleft in hevamine.     

X-ray analyses of aspartic proteinases. The three-dimensional structure at 2.1 A resolution of endothiapepsin

The molecular structure of endothiapepsin (EC 3.4.23.6), the aspartic proteinase from Endothia parasitica, has been refined to a crystallographic R-factor of 0.178 at 2.1 A resolution. The positions of 2389 protein non-hydrogen atoms have been determined and the present model contains 333 solvent molecules. The structure is bilobal, consisting of two predominantly beta-sheet domains that are related by an approximate 2-fold axis. Of approximately 170 residues, 65 are topologically equivalent when one lobe is superimposed on the other. Twenty beta-strands are arranged as five beta-sheets and are connected by regions involving 29 turns and four helices. A central sheet involves three antiparallel strands from each lobe organized around the dyad axis. Each lobe contains a further local dyad that passes through two sheets arranged as a sandwich and relates two equivalent motifs of four antiparallel strands (a, b, c, d) followed by a helix or an irregular helical region. Sheets 1N and 1C, each contain two interpenetrating psi structures contributed by strands c,d,d' and c',d',d, which are related by the intralobe dyad. A further sheet, 2N or 2C, is formed from two extended beta-hairpins from strands b,c and b',c' that fold above the sheets 1N and 1C, respectively, and are hydrogen-bonded around the local intralobe dyad. Asp32 and Asp215 are related by the interlobe dyad and form an intricate hydrogen-bonded network with the neighbouring residues and comprise the most symmetrical part of the structure. The side-chains of the active site aspartate residues are held coplanar and the nearby main chain makes a "fireman's grip" hydrogen-bonding network. Residues 74 to 83 from strands a'N and b'N in the N-terminal lobe form a beta-hairpin loop with high thermal parameters. This "flap" projects over the active site cleft and shields the active site from the solvent region. Shells of water molecules are found on the surface of the protein molecule and large solvent channels are observed within the crystal. There are only three regions of intermolecular contacts and the crystal packing is stabilized by many solvent molecules forming a network of hydrogen bonds. The three-dimensional structure of endothiapepsin is found to be similar to two other fungal aspartic proteinases, penicillopepsin and rhizopuspepsin. Even though sequence identities of endothiapepsin with rhizopuspepsin and penicillopepsin are only 41% and 51%, respectively, a superposition of the three-dimensional structures of these three enzymes shows that 237 residues (72%) are within a root-mean-square distance of 1.0 A.