Crystal structure of D-amino acid oxidase: a case of active site mirror-image convergent evolution with flavocytochrome b2.

D-amino acid oxidase is the prototype of the FAD-dependent oxidases. It catalyses the oxidation of D-amino acids to the corresponding alpha-ketoacids. The reducing equivalents are transferred to molecular oxygen with production of hydrogen peroxide. We have solved the crystal structure of the complex of D-amino acid oxidase with benzoate, a competitive inhibitor of the substrate, by single isomorphous replacement and eightfold averaging. Each monomer is formed by two domains with an overall topology similar to that of p-hydroxybenzoate hydroxylase. The benzoate molecule lays parallel to the flavin ring and is held in position by a salt bridge with Arg-283. Analysis of the active site shows that no side chains are properly positioned to act as the postulated base required for the catalytic carboanion mechanism. On the contrary, the benzoate binding mode suggests a direct transfer of the substrate alpha-hydrogen to the flavin during the enzyme reductive half-reaction.The active site Of D-amino acid oxidase exhibits a striking similarity with that of flavocytochrome b2, a structurally unrelated FMN-dependent flavoenzyme. The active site groups (if these two enzymes are in fact superimposable once the mirror-image of the flavocytochrome b2 active site is generated with respect to the flavin plane. Therefore, the catalytic sites of D-amino acid oxidase and flavocytochrome b2 appear to have converged to a highly similar but enantiomeric architecture in order to catalvze similar reactions (oxidation of alpha-amino acids or alpha-hydroxy acids), although with opposite stereochemistry.

Molecular cloning of violaxanthin de-epoxidase from romaine lettuce and expression in Escherichia coli.

Plants need to avoid or dissipate excess light energy to protect photosystem II (PSII) from photoinhibitory damage. Higher plants have a conserved system that dissipates excess energy as heat in the light-harvesting complexes of PSII that depends on the transthylakoid delta pH and violaxanthin de-epoxidase (VDE) activity. To our knowledge, we report the first cloning of a cDNA encoding VDE and expression of functional enzyme in Escherichia coli. VDE is nuclear encoded and has a transit peptide with characteristic features of other lumen-localized proteins. The cDNA encodes a putative polypeptide of 473 aa with a calculated molecular mass of 54,447 Da. Cleavage of the transit peptide results in a mature putative polypeptide of 348 aa with a calculated molecular mass of 39,929 Da, close to the apparent mass of the purified enzyme (43 kDa). The protein has three interesting domains including (i) a cysteine-rich region, (ii) a lipocalin signature, and (iii) a highly charged region. The E. coli expressed enzyme de-epoxidizes violaxanthin sequentially to antheraxanthin and zeaxanthin, and is inhibited by dithiothreitol, similar to VDE purified from chloroplasts. This confirms that the cDNA encodes an authentic VDE of a higher plant and is unequivocal evidence that the same enzyme catalyzes the two-step mono de-epoxidation reaction. The cloning of VDE opens new opportunities for examining the function and evolution of the xanthophyll cycle, and possibly enhancing light-stress tolerance of plants.

A new high molecular weight immunoglobulin class from the carcharhine shark: implications for the properties of the primordial immunoglobulin.

All immunoglobulins and T-cell receptors throughout phylogeny share regions of highly conserved amino acid sequence. To identify possible primitive immunoglobulins and immunoglobulin-like molecules, we utilized 3' RACE (rapid amplification of cDNA ends) and a highly conserved constant region consensus amino acid sequence to isolate a new immunoglobulin class from the sandbar shark Carcharhinus plumbeus. The immunoglobulin, termed IgW, in its secreted form consists of 782 amino acids and is expressed in both the thymus and the spleen. The molecule overall most closely resembles mu chains of the skate and human and a new putative antigen binding molecule isolated from the nurse shark (NAR). The full-length IgW chain has a variable region resembling human and shark heavy-chain (VH) sequences and a novel joining segment containing the WGXGT motif characteristic of H chains. However, unlike any other H-chain-type molecule, it contains six constant (C) domains. The first C domain contains the cysteine residue characteristic of C mu1 that would allow dimerization with a light (L) chain. The fourth and sixth domains also contain comparable cysteines that would enable dimerization with other H chains or homodimerization. Comparison of the sequences of IgW V and C domains shows homology greater than that found in comparisons among VH and C mu or VL, or CL thereby suggesting that IgW may retain features of the primordial immunoglobulin in evolution.

The molecular biology of apoptosis.

All multicellular organisms have mechanisms for killing their own cells, and use physiological cell death for defence, development, homeostasis, and aging. Apoptosis is a morphologically recognizable form of cell death that is implemented by a mechanism that has been conserved throughout evolution from nematode to man. Thus homologs of the genes that implement cell death in nematodes also do so in mammals, but in mammals the process is considerably more complex, involving multiple isoforms of the components of the cell death machinery. In some circumstances this allows independent regulation of pathways that converge upon a common end point. A molecular understanding of this mechanism may allow design of therapies that either enhance or block cell death at will.

Population dynamics of flaviviruses revealed by molecular phylogenies.

The phylogeny of 123 complete envelope gene sequences was reconstructed in order to understand the evolution of tick- and mosquito-borne flaviviruses. An analysis of phylogenetic tree structure reveals a continual and asymmetric branching process in the tick-borne flaviviruses, compared with an explosive radiation in the last 200 years in viruses transmitted by mosquitoes. The distinction between these two viral groups probably reflects differences in modes of dispersal, propagation, and changes in the size of host populations. The most serious implication of this work is that growing human populations are being exposed to an expanding range of increasingly diverse viral strains.

ndhF sequence evolution and the major clades in the sunflower family.

An extensive sequence comparison of the chloroplast ndhF gene from all major clades of the largest flowering plant family (Asteraceae) shows that this gene provides approximately 3 times more phylogenetic information than rbcL. This is because it is substantially longer and evolves twice as fast. The 5' region (1380 bp) of ndhF is very different from the 3' region (855 bp) and is similar to rbcL in both the rate and the pattern of sequence change. The 3' region is more A+T-rich, has higher levels of nonsynonymous base substitution, and shows greater transversion bias at all codon positions. These differences probably reflect different functional constraints on the 5' and 3' regions of ndhF. The two patterns of base substitutions of ndhF are particularly advantageous for phylogenetic reconstruction because the conserved and variable segments can be used for older and recent groups, respectively. Phylogenetic analyses of 94 ndhF sequences provided much better resolution of relationships than previous molecular and morphological phylogenies of the Asteraceae. The ndhF tree identified five major clades: (i) the Calyceraceae is the sister family of Asteraceae; (ii) the Barnadesioideae is monophyletic and is the sister group to the rest of the family; (iii) the Cichorioideae and its two basal tribes Mutisieae and Cardueae are paraphyletic; (iv) four tribes of Cichorioideae (Lactuceae, Arctoteae, Liabeae, and Vernonieae) form a monophyletic group, and these are the sister clade of the Asteroideae; and (v) the Asteroideae is monophyletic and includes three major clades.

Certain observations concerning the effects of epistasis on complex traits and the evolution of genomes.

Thesis (Ph.D.)--University of Washington, 2016-05

Functional recycling of C2 domains throughout evolution: A comparative study of synaptotagmin, protein kinase C and phospholipase C by sequence, structural and modelling approaches

The C2 domain is one of the most frequent and widely distributed calcium-binding motifs. Its structure comprises an eight-stranded beta-sandwich with two structural types as if the result of a circular permutation. Combining sequence, structural and modelling information, we have explored, at different levels of granularity, the functional characteristics of several families of C2 domains. At the coarsest level,the similarity correlates with key structural determinants of the C2 domain fold and, at the finest level, with the domain architecture of the proteins containing them, highlighting the functional diversity between the various subfamilies. The functional diversity appears as different conserved surface patches throughout this common fold. In some cases, these patches are related to substrate-binding sites whereas in others they correspond to interfaces of presumably permanent interaction between other domains within the same polypeptide chain. For those related to substrate-binding sites, the predictions overlap with biochemical data in addition to providing some novel observations. For those acting as protein-protein interfaces' our modelling analysis suggests that slight variations between families are a result of not only complementary adaptations in the interfaces involved but also different domain architecture. In the light of the sequence and structural genomic projects, the work presented here shows that modelling approaches along with careful sub-typing of protein families will be a powerful combination for a broader coverage in proteomics. (C) 2003 Elsevier Ltd. All rights reserved.

Molecular phylogeny and biogeography of the dung beetle genus Temnoplectron Westwood (Scarabaeidae: Scarabaeinae) from Australia’s wet tropics

The landscape of the Australian Wet Tropics can be described as islands of montane rainforest Surrounded by warmer or more xeric habitats. Historical glaciation cycles have caused expansion and contraction of these rainforest islands leading to consistent patterns of genetic divergence within species of vertebrates. To explore whether this dynamic history has promoted speciation in endemic and diverse groups Of insects, we used a combination of mtDNA sequencing and morphological characters to estimate relationships and the tempo of divergence among Australian representatives of the dung beetle genus Temnoplectron. This phylogenetic hypothesis shares a number of well-supported clades with a previously published phylogenetic hypothesis based on morphological data. though statistical support for several nodes is weak. Sister species relationships well-supported in both tree topologies. and a tree obtained by combining the two data sets. suggest that speciation has mostly been allopatric. We identify a number of speciation barriers, which coincide with phylogeographic breaks found in vertebrate species. Large sequence divergences between species emphasize that speciation events are ancient (pre-Pleistocene). The flightless, rainforest species appear to have speciated rapidly. but also in the distant past. (C) 2003 Elsevier Inc. All rights reserved.

Selection pressure-driven evolution of the Epstein-Barr virus-encoded oncogene LMP1 in virus isolates from southeast Asia

The geographically constrained distribution of Epstein-Barr virus (EBV)-associated nasopharyngeal carcinoma (NPC) in southeast Asian populations suggests that both viral and host genetics may influence disease risk. Although susceptibility loci have been mapped within the human genome, the role of viral genetics in the focal distribution of NPC remains an enigma. Here we report a molecular phylogenetic analysis of an NPC-associated viral oncogene, LMP1, in a large panel of EBV isolates from southeast Asia and from Papua New Guinea, Africa, and Australia, regions of the world where NPC is and is not endemic, respectively. This analysis revealed that LMP1 sequences show a distinct geographic structure, indicating that the southeast Asian isolates have evolved as a lineage distinct from those of Papua New Guinea, African, and Australian isolates. Furthermore, a likelihood ratio test revealed that the C termini of the LMP1 sequences of the southeast Asian lineage are under significant positive selection pressure, particularly at some sites within the C-terminal activator regions. We also present evidence that although the N terminus and transmembrane region of LMP1 have undergone recombination, the C-terminal region of the gene has evolved without any history of recombination. Based on these observations, we speculate that selection pressure may be driving the LMP1 sequences in virus isolates from southeast Asia towards a more malignant phenotype, thereby influencing the endemic distribution of NPC in this region.

Evolution of polydnaviruses as insect immune suppressors

Polydnaviruses (PDVs) are endogenous particles that are used by some endoparasitic hymenoptera to disrupt host immunity and development. Recent analyses of encapsidated PDV genes have increased the number of known PDV gene families, which are often closely related to insect genes. Several PDV proteins inactivate host haemocytes by damaging their actin cytoskeleton. These proteins share no significant sequence homology and occur in polyphyletic PDV genera, possibly indicating that convergent evolution has produced functionally similar immune-suppressive molecules causing a haemocyte phenotype characterised by damaged cytoskeleton and inactivation. These phenomena provide further insights into the immune-suppressive activity of PDVs and raise interesting questions about PDV evolution, a topic that has puzzled researchers ever since the discovery of PDVs.

Evolution of the regular zone of histone H1 in Fabaceae plants

An analysis of the historic H1 subtype, H1-1, in eight legumes belonging to four genera of the tribe Vicieae (Pisum, Lathyrus, Lens, and Vicia), revealed an extended region consisting of the tandemly repeated AKPAAK motifs. We named this region the Regular zone (RZ). The AKPAAK motifs are organized into two blocks separated by a short (two or six amino acids) intervening sequence (IS). The distal block contains six AKPAAK motifs, while the number of repeats in the proximal block varies from six in V. faba to seven in the other species. In V. hirsuta, the first two repeated units of the proximal block are octapeptides AKAKPAAK. The apparent rate of synonymous substitutions in the blocks of RZ is much higher than in the rest of the gene. This can be explained by repeat shuffling within each block. In the C-domain of the orthologous H1 subtype froth Medicago truncatula (tribe Trifolieae), a region corresponding to the RZ of Vicieae species was found. It also consists of two blocks of AKPAAK motifs (four and three repeats in the proximal and distal blocks, respectively). These blocks are separated by a 20-amino acid IS. The first 20 amino acids of the Medicago RZ are not part of AKPAAK repeats. We hypothesise that the RZ has most probably evolved as a result of an expansion of AKPAAK repeats from two separate sites in the C-domain. This process started tens of millions of years ago and was most likely directed by positive selection.

Evolution of resistance to sulfadoxine-pyrimethamine in Plasmodium falciparum

The development of resistance to sulfadoxine-pyrimethamine by Plasmodium parasites is a major problem for the effective treatment of malaria, especially P. falciparum malaria. Although the molecular basis for parasite resistance is known, the factors promoting the development and transmission of these resistant parasites are less clear. This paper reports the results of a quantitative comparison of factors previously hypothesized as important for the development of drug resistance, drug dosage, time of treatment, and drug elimination half-life, with an in-host dynamics model of P. falciparum malaria in a malaria-naive host. The results indicate that the development of drug resistance can be categorized into three stages. The first is the selection of existing parasites with genetic mutations in the dihydrofolate reductase or dihydropteroate synthetase gene. This selection is driven by the long half-life of the sulfadoxine-pyrimethamine combination. The second stage involves the selection of parasites with allelic types of higher resistance within the host during an infection. The timing of treatment relative to initiation of a specific anti-P. falciparum EMP1 immune response is an important factor during this stage, as is the treatment dosage. During the third stage, clinical treatment failure becomes prevalent as the parasites develop sufficient resistance mutations to survive therapeutic doses of the drug combination. Therefore, the model output reaffirms the importance of correct treatment of confirmed malaria cases in slowing the development of parasite resistance to sulfadoxine-pyrimethamine.

Design strategies for controlling the molecular weight and rate using reversible addition-fragmentation chain transfer mediated living radical polymerization

Living radical polymerization has allowed complex polymer architectures to be synthesized in bulk, solution, and water. The most versatile of these techniques is reversible addition-fragmentation chain transfer (RAFT), which allows a wide range of functional and nonfunctional polymers to be made with predictable molecular weight distributions (MWDs), ranging from very narrow to quite broad. The great complexity of the RAFT mechanism and how the kinetic parameters affect the rate of polymerization and MWD are not obvious. Therefore, the aim of this article is to provide useful insights into the important kinetic parameters that control the rate of polymerization and the evolution of the MWD with conversion. We discuss how a change in the chain-transfer constant can affect the evolution of the MWD. It is shown how we can, in principle, use only one RAFT agent to obtain a poly-mer with any MWD. Retardation and inhibition are discussed in terms of (1) the leaving R group reactivity and (2) the intermediate radical termination model versus the slow fragmentation model. (c) 2005 Wiley Periodicals, Inc.

Molecular analysis reveals tighter social regulation of immigration in patrilocal populations than in matrilocal populations

Human social organization can deeply affect levels of genetic diversity. This fact implies that genetic information can be used to study social structures, which is the basis of ethnogenetics. Recently, methods have been developed to extract this information from genetic data gathered from subdivided populations that have gone through recent spatial expansions, which is typical of most human populations. Here, we perform a Bayesian analysis of mitochondrial and Y chromosome diversity in three matrilocal and three patrilocal groups from northern Thailand to infer the number of males and females arriving in these populations each generation and to estimate the age of their range expansion. We find that the number of male immigrants is 8 times smaller in patrilocal populations than in matrilocal populations, whereas women move 2.5 times more in patrilocal populations than in matrilocal populations. In addition to providing genetic quantification of sex-specific dispersal rates in human populations, we show that although men and women are exchanged at a similar rate between matrilocal populations, there are far fewer men than women moving into patrilocal populations. This finding is compatible with the hypothesis that men are strictly controlling male immigration and promoting female immigration in patrilocal populations and that immigration is much less regulated in matrilocal populations.