Complex evolution of orthologous and paralogous decarboxylase genes

The decarboxylases are involved in neurotransmitter synthesis in animals, and in pathways of secondary metabolism in plants. Different decarboxylase proteins are characterized for their different substrate specificities, but are encoded by homologous genes. We study, within a maximum-likelihood framework, the evolutionary relationships among dopa decarboxylase (Ddc), histidine decarboxylase (Hdc) and alpha-methyldopa hypersensitive (amd) in animals, and tryptophan decarboxylase (Wdc) and tyrosine decarboxylase (Ydc) in plants. The evolutionary rates are heterogeneous. There are differences between paralogous genes in the same lineages: 4.13 x 10(-10) nucleotide substitutions per site per year in mammalian Ddc vs. 1.95 in Hdc; between orthologous genes in different lineages, 7.62 in dipteran Ddc vs. 4.13 in mammalian Ddc; and very large temporal variations in some lineages, from 3.7 up to 54.9 in the Drosophila Ddc lineage. Our results are inconsistent with the molecular clock hypothesis.     

Unexpected foraminiferal diversity revealed by small-subunit rDNA analysis of Antarctic sediment

Studies of benthic Foraminifera typically rely on the morphological identification of dried specimens. This approach can introduce sampling bias against small, delicate, or morphologically ambiguous forms. To overcome this limitation, we extracted total DNA from sediment followed by PCR using group- and species-specific primers. Phylogenetic analyses revealed that approximately ninety percent of the PCR products represented previously undescribed sequence types that group with undersampled members of the allogromiid Foraminifera. We also used a modification of this technique to track individual species in sediment fractions too fine for normal morphological identification, and to confirm species placement of morphologically ambiguous foraminiferans. We were able to identify the DNA of several large foraminiferal species in fine fractions in a seasonally-dependent manner, indicating that in some seasons the majority of the standing stock of these species exists as gametes/juveniles. The approach outlined here represents a powerful strategy for exploring the total diversity of benthic foraminiferal communities.     

The comet assay for DNA damage and repair

The comet assay (single-cell gel electrophoresis) is a simple method for measuring deoxyribonucleic acid (DNA) strand breaks in eukaryotic cells. Cells embedded in agarose on a microscope slide are lysed with detergent and high salt to form nucleoids containing supercoiled loops of DNA linked to the nuclear matrix. Electrophoresis at high pH results in structures resembling comets, observed by fluorescence microscopy; the intensity of the comet tail relative to the head reflects the number of DNA breaks. The likely basis for this is that loops containing a break lose their supercoiling and become free to extend toward the anode. The assay has applications in testing novel chemicals for genotoxicity, monitoring environmental contamination with genotoxins, human biomonitoring and molecular epidemiology, and fundamental research in DNA damage and repair. The sensitivity and specificity of the assay are greatly enhanced if the nucleoids are incubated with bacterial repair endonucleases that recognize specific kinds of damage in the DNA and convert lesions to DNA breaks, increasing the amount of DNA in the comet tail. DNA repair can be monitored by incubating cells after treatment with damaging agent and measuring the damage remaining at intervals. Alternatively, the repair activity in a cell extract can be measured by incubating it with nucleoids containing specific damage.     

N-acetylaspartate in the vertebrate brain: metabolism and function

N-Acetyl-l-aspartate (NAA) is an amino acid that is present in the vertebrate brain. Its concentration is one of the highest of all free amino acids and, although NAA is synthesized and stored primarily in neurons, it cannot be hydrolyzed in these cells. Furthermore, neuronal NAA is dynamic and turns over more than once each day by virtue of its continuous efflux, in a regulated intercompartmental cycling via extracellular fluids, between neurons and a second compartment in oligodendrocytes. The metabolism of NAA, between its anabolic compartment in neurons and its catabolic compartment in oligodendrocytes, and its possible physiological role in the brain has been the subject of much speculation. There are two human inborn errors in metabolism of NAA. One is Canavan disease (CD), in which there is a buildup of NAA (hyperacetylaspartia) and associated spongiform leukodystrophy, caused by a lack of aspartoacylase activity. The other is a singular human case of lack of NAA (hypoacetylaspartia), where the enzyme that synthesizes NAA is apparently absent. There are two animal models currently available for studies of CD. One is a rat with a natural deletion of the catabolic enzyme, and the other a gene knockout mouse. In addition to the presence of NAA in neurons, its prominence in ¹H nuclear magnetic resonance spectroscopic studies has led to its wide use in diagnostic human medicine as both an indicator of brain pathology and of disease progression in a variety of CNS diseases. In this review, various hypotheses regarding the metabolism of NAA and its possible role in the CNS are evaluated. Based on this analysis, it is concluded that although NAA may have several functions in the CNS, an important role of the NAA intercompartmental system is osmoregulatory, and in this role it may be the primary mechanism for the removal of intracellular water, against a water gradient, from myelinated neurons.     

A path from primary protein sequence to ligand recognition

A novel method to organize protein structural information based solely on sequence is presented. The method clusters proteins into families that correlate with the three-dimensional protein structure and the conformation of the bound ligands. This procedure was applied to nicotinamide adenine dinucleotide [NAD(P)]-utilizing enzymes to identify a total of 94 sequence families, 53 of which are structurally characterized. Each of the structurally characterized proteins within a sequence family correlates to a single protein fold and to a common bound conformation of NAD(P). A wide range of structural folds is identified that recognize NAD(P), including Rossmann folds and beta/alpha barrels. The defined sequence families can be used to identify the type and prevalence of NAD(P)-utilizing enzymes in the proteomes of sequenced organisms. The proteome of Mycobacterium tuberculosis was mined to generate a proteome-wide profile of NAD(P)-utilizing enzymes coded by this organism. This enzyme family comprises approximately 6% of the open reading frames, with the largest subgroup being the Rossmann fold, short-chain dehydrogenases. The preponderance of short-chain dehydrogenases correlates strongly with the phenotype of M. tuberculosis, which is characterized as having one of the most complex prokaryotic cell walls.     

Evolutionary analysis reveals collective properties and specificity in the C-type lectin and lectin-like domain superfamily

Members of the C-type lectin/C-type lectin-like domain (CTL/CTLD) superfamily share a common fold and are involved in a variety of functions, such as generalized defense mechanisms against foreign agents, discrimination between healthy and pathogen-infected cells, and endocytosis and blood coagulation. In this work we used ConSurf, a computer program recently developed in our lab, to perform an evolutionary analysis of this superfamily in order to further identify characteristics of all or part of its members. Given a set of homologous proteins in the form of multiple sequence alignment (MSA) and an inferred phylogenetic tree, ConSurf calculates the conservation score in every alignment position, taking into account the relationships between the sequences and the physicochemical similarity between the amino acids. The scores are then color-coded onto the three-dimensional structure of one of the homologous proteins. We provide here and at http://ashtoret.tau.ac.il/ approximately sharon a detailed analysis of the conservation pattern obtained for the entire superfamily and for two subgroups of proteins: (a) 21 CTLs and (b) 11 heterodimeric CTLD toxins. We show that, in general, proteins of the superfamily have one face that is constructed mostly of conserved residues and another that is not, and we suggest that the former face is involved in binding to other proteins or domains. In the CTLs examined we detected a region of highly conserved residues, corresponding to the known calcium- and carbohydrate-binding site of the family, which is not conserved throughout the entire superfamily, and in the CTLD toxins we found a patch of highly conserved residues, corresponding to the known dimerization region of these proteins. Our analysis also detected patches of conserved residues with yet unknown function(s).     

Free-energy landscape of a chameleon sequence in explicit water and its inherent alpha/beta bifacial property

A sequence in yeast MATalpha2/MCM1/DNA complex that folds into alpha-helix or beta-hairpin depending on the surroundings has been known as "chameleon" sequence. We obtained the free-energy landscape of this sequence by using a generalized-ensemble method, multicanonical molecular dynamics simulation, to sample the conformational space. The system was expressed with an all-atom model in explicit water, and the initial conformation for the simulation was a random one. The free-energy landscape demonstrated that this sequence inherently has an ability to form either alpha or beta structure: The conformational distribution in the landscape consisted of two alpha-helical clusters with different packing patterns of hydrophobic residues, and four beta-hairpin clusters with different strand-strand interaction patterns. Narrow pathways connecting the clusters were found, and analysis on the pathways showed that a compact structure formed at the N-terminal root of the chameleon sequence controls the cluster-cluster transitions. The free-energy landscape indicates that a small conditional change induces alpha-beta transitions. Additional unfolding simulations done with replacing amino acids showed that the chameleon sequence has an advantage to form an alpha-helix. Current study may be useful to understand the mechanism of diseases resulting from abnormal chain folding, such as amyloid disease.     

Mapping Hsp47 binding site(s) using CNBr peptides derived from type I and type II collagen

As a crucial molecular chaperone in collagen biosynthesis, Hsp47 interacts with the nascent form as well as the mature triple-helical form of procollagen. The location(s) of Hsp47 binding sites on the collagen molecule are, as yet, unknown. We have examined the substrate specificity of Hsp47 in vitro using well-characterized CNBr peptide fragments of type I and type II collagen along with radiolabeled, recombinant Hsp47. Interaction of these peptides with Hsp47 bound to collagen-coated microtiter wells showed several binding sites for Hsp47 along the length of the alpha1 and alpha2 chains of type I collagen and the alpha1 chain of type II collagen, with the N-terminal regions showing the strongest affinities. The latter observation was also supported by the results of a ligand-blot assay. Except for two peptides in the alpha2(I) chain, peptides that showed substantial binding to Hsp47 did so in their triple-helical and not random-coil form. Unlike earlier studies that used peptide models for collagen, the results obtained here on fragments of type I and type II collagen identify, for the first time, binding of Hsp47 to specific regions of the collagen molecule. They also point to additional structural requirements for Hsp47 binding besides the known preference for third-position Arg residues and the triple-helical conformation.     

Genetic structure in a montane ranid frog: restricted gene flow and nuclear-mitochondrial discordance

There is substantial debate over the criteria that should be used to group populations of a species into distinct units for conservation (e.g. evolutionarily significant units, management units, distinct population segments). However, in practice molecular genetic differentiation is often the only or main criterion used to identify such units. Most genetic studies attempting to define conservation units in animals use a single molecular marker, most often mitochondrial, and use samples from a limited number of populations throughout the species' range. Although there are many benefits to using mtDNA, certain features can cause it to show patterns of differentiation among populations that do not reflect the history of differentiation at the nuclear genome where loci controlling traits of adaptive significance presumably occur. Here we illustrate an example of such mitochondrial-nuclear discordance in a ranid frog, and show how using mtDNA or nuclear loci alone could have led to very different conservation recommendations. We also found very high genetic differentiation among populations on a local scale, and discuss the conservation implications of our results.