Measurement of genetic structure within populations using Moran's spatial autocorrelation statistics.

Spatial structure of genetic variation within populations, an important interacting influence on evolutionary and ecological processes, can be analyzed in detail by using spatial autocorrelation statistics. This paper characterizes the statistical properties of spatial autocorrelation statistics in this context and develops estimators of gene dispersal based on data on standing patterns of genetic variation. Large numbers of Monte Carlo simulations and a wide variety of sampling strategies are utilized. The results show that spatial autocorrelation statistics are highly predictable and informative. Thus, strong hypothesis tests for neutral theory can be formulated. Most strikingly, robust estimators of gene dispersal can be obtained with practical sample sizes. Details about optimal sampling strategies are also described.

Random mutagenesis of Thermus aquaticus DNA polymerase I: concordance of immutable sites in vivo with the crystal structure.

Expression of Thermus aquaticus (Taq) DNA polymerase I (pol I) in Escherichia, coli complements the growth defect caused by a temperature-sensitive mutation in the host pol I. We replaced the nucleotide sequence encoding amino acids 659-671 of the O-helix of Taq DNA pol I, corresponding to the substrate binding site, with an oligonucleotide containing random nucleotides. Functional Taq pol I mutants were selected based on colony formation at the nonpermissive temperature. By using a library with 9% random substitutions at each of 39 positions, we identified 61 active Taq pol I mutants, each of which contained from one to four amino acid substitutions. Some amino acids, such as alanine-661 and threonine-664, were tolerant of several or even many diverse replacements. In contrast, no replacements or only conservative replacements were identified at arginine-659, lysine-663, and tyrosine-671. By using a library with totally random nucleotides at five different codons (arginine-659, arginine-660, lysine-663, phenylalanine-667, and glycine-668), we confirmed that arginine-659 and lysine-663 were immutable, and observed that only tyrosine substituted for phenylalanine-667. The two immutable residues and the two residues that tolerate only highly conservative replacements lie on the side of O-helix facing the incoming deoxynucleoside triphosphate, as determined by x-ray analysis. Thus, we offer a new approach to assess concordance of the active conformation of an enzyme, as interpreted from the crystal structure, with the active conformation inferred from in vivo function.

Structure, tissue distribution, and chromosomal localization of the prepronociceptin gene.

Nociceptin (orphanin FQ), the newly discovered natural agonist of opioid receptor-like (ORL1) receptor, is a neuropeptide that is endowed with pronociceptive activity in vivo. Nociceptin is derived from a larger precursor, prepronociceptin (PPNOC), whose human, mouse, and rat genes we have now isolated. The PPNOC gene is highly conserved in the three species and displays organizational features that are strikingly similar to those of the genes of preproenkephalin, preprodynorphin, and preproopiomelanocortin, the precursors to endogenous opioid peptides, suggesting the four genes belong to the same family-i.e., have a common evolutionary origin. The PPNOC gene encodes a single copy of nociceptin as well as of other peptides whose sequence is strictly conserved across murine and human species; hence it is likely to be neurophysiologically significant. Northern blot analysis shows that the PPNOC gene is predominantly transcribed in the central nervous system (brain and spinal cord) and, albeit weakly, in the ovary, the sole peripheral organ expressing the gene. By using a radiation hybrid cell line panel, the PPNOC gene was mapped to the short arm of human chromosome 8 (8p21), between sequence-tagged site markers WI-5833 and WI-1172, in close proximity of the locus encoding the neurofilament light chain NEFL. Analysis of yeast artificial chromosome clones belonging to the WC8.4 contig covering the 8p21 region did not allow to detect the presence of the gene on these yeast artificial chromosomes, suggesting a gap in the coverage within this contig.

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.

Highly stable expression of a foreign gene from rabies virus vectors.

A reverse genetics approach was applied to generate a chimeric nonsegmented negative strand RNA virus, rabies virus (RV) of the Rhabdoviridae family, that expresses a foreign protein. DNA constructs containing the entire open reading frame of the bacterial chloramphenicol acetyltransferase (CAT) gene and an upstream RV cistron border sequence were inserted either into the nontranslated pseudogene region of a full-length cDNA copy of the RV genome or exchanged with the pseudogene region. After intracellular T7 RNA polymerase-driven expression of full-length antigenome RNA transcripts and RV nucleoprotein, phosphoprotein and polymerase from transfected plasmids, RVs transcribing novel monocistronic mRNAs and expressing CAT at high levels, were recovered. The chimeric viruses possessed the growth characteristics of standard RV and were genetically stable upon serial cell culture passages. CAT activity was still observed in cell cultures infected with viruses passaged for more than 25 times. Based on the unprecedented stability of the chimeric RNA genomes, which is most likely due to the structure of the rhabdoviral ribonucleoprotein complex, we predict the successful future use of recombinant rhabdovirus vectors for displaying foreign antigens or delivering therapeutic genes.

High yield conversion of doxorubicin to 2-pyrrolinodoxorubicin, an analog 500-1000 times more potent: structure-activity relationship of daunosamine-modified derivatives of doxorubicin.

A convenient, high yield conversion of doxorubicin to 3'-deamino-3'-(2''-pyrroline-1''-yl)doxorubicin is described. This daunosamine-modified analog of doxorubicin is 500-1000 times more active in vitro than doxorubicin. The conversion is effected by using a 30-fold excess of 4-iodobutyraldehyde in anhydrous dimethylformamide. The yield is higher than 85%. A homolog of this compound, 3'-deamino-3'-(1'',3''-tetrahydropyridine-1''-yl)doxorubicin, was also synthesized by using 5-iodovaleraldehyde. In this homolog, the daunosamine nitrogen is incorporated into a six- instead of a five-membered ring. This analog was 30-50 times less active than its counterpart with a five-membered ring. A similar structure-activity relationship was found when 3'-deamino-3'-(3''-pyrrolidone-1''-yl)doxorubicin (containing a five-membered ring) and 3'-deamino-3'-(3''-piperidone-1''-yl)doxorubicin (with a six-membered ring) were tested in vitro, the former being 5 times more potent than the latter. To further elucidate structure-activity relationships, 3'-deamino-3'-(pyrrolidine-1''-yl)doxorubicin, 3'-deamino-3'-(isoindoline-2''-yl)doxorubicin, 3'-deamino-3'-(2''-methyl-2''-pyrroline-1''-yl)doxorubicin, and 3'-deamino-3'-(3''-pyrroline-1''-yl)doxorubicin were also synthesized and tested. All the analogs were prepared by using reactive halogen compounds for incorporating the daunosamine nitrogen of doxorubicin into a five- or six-membered ring. These highly active antineoplastic agents can be used for incorporation into targeted cytotoxic analogs of luteinizing hormone-releasing hormone intended for cancer therapy.

The wing of the enhancer-binding domain of Mu phage transposase is flexible and is essential for efficient transposition.

A tetramer of the Mu transposase (MuA) pairs the recombination sites, cleaves the donor DNA, and joins these ends to a target DNA by strand transfer. Juxtaposition of the recombination sites is accomplished by the assembly of a stable synaptic complex of MuA protein and Mu DNA. This initial critical step is facilitated by the transient binding of the N-terminal domain of MuA to an enhancer DNA element within the Mu genome (called the internal activation sequence, IAS). Recently we solved the three-dimensional solution structure of the enhancer-binding domain of Mu phage transposase (residues 1-76, MuA76) and proposed a model for its interaction with the IAS element. Site-directed mutagenesis coupled with an in vitro transposition assay has been used to assess the validity of the model. We have identified five residues on the surface of MuA that are crucial for stable synaptic complex formation but dispensable for subsequent events in transposition. These mutations are located in the loop (wing) structure and recognition helix of the MuA76 domain of the transposase and do not seriously perturb the structure of the domain. Furthermore, in order to understand the dynamic behavior of the MuA76 domain prior to stable synaptic complex formation, we have measured heteronuclear 15N relaxation rates for the unbound MuA76 domain. In the DNA free state the backbone atoms of the helix-turn-helix motif are generally immobilized whereas the residues in the wing are highly flexible on the pico- to nanosecond time scale. Together these studies define the surface of MuA required for enhancement of transposition in vitro and suggest that a flexible loop in the MuA protein required for DNA recognition may become structurally ordered only upon DNA binding.

Primary structure and functional expression of a cGMP-gated potassium channel.

Cyclic nucleotides modulate potassium (K) channel activity in many cells and are thought to act indirectly by inducing channel protein phosphorylation. Herein we report the isolation from rabbit of a gene encoding a K channel (Kcn1) that is specifically activated by cGMP and not by cAMP. Analysis of the deduced amino acid sequence (725 amino acids) indicates that, in addition to a core region that is highly homologous to Shaker K channels, Kcn1 also contains a cysteine-rich region similar to that of ligand-gated ion channels and a cyclic nucleotide-binding region. Northern blot analysis detects gene expression in kidney, aorta, and brain. Kcn1 represents a class of K channels that may be specifically regulated by cGMP and could play an important role in mediating the effects of substances, such as nitric oxide, that increase intracellular cGMP.

Lyophilization-induced reversible changes in the secondary structure of proteins.

Changes in the secondary structure of some dozen different proteins upon lyophilization of their aqueous solutions have been investigated by means of Fourier-transform infrared spectroscopy in the amide III band region. Dehydration markedly (but reversibly) alters the secondary structure of all the proteins studied, as revealed by both the quantitative analysis of the second derivative spectra and the Gaussian curve fitting of the original infrared spectra. Lyophilization substantially increases the beta-sheet content and lowers the alpha-helix content of all proteins. In all but one case, proteins become more ordered upon lyophilization.

Receptive field structure in the visual cortex: does selective stimulation induce plasticity?

Sensory areas of adult cerebral cortex can reorganize in response to long-term alterations in patterns of afferent signals. This long-term plasticity is thought to play a crucial role in recovery from injury and in some forms of learning. However, the degree to which sensory representations in primary cortical areas depend on short-term (i.e., minute to minute) stimulus variations remains unclear. A traditional view is that each neuron in the mature cortex has a fixed receptive field structure. An alternative view, with fundamentally different implications for understanding cortical function, is that each cell's receptive field is highly malleable, changing according to the recent history of the sensory environment. Consistent with the latter view, it has been reported that selective stimulation of regions surrounding the receptive field induces a dramatic short-term increase in receptive field size for neurons in the visual cortex [Pettet, M. W. & Gilbert, C. D. (1992) Proc. Natl. Acad. Sci. USA 89, 8366-8370]. In contrast, we report here that there is no change in either the size or the internal structure of the receptive field following several minutes of surround stimulation. However, for some cells, overall responsiveness increases. These results suggest that dynamic alterations of receptive field structure do not underlie short-term plasticity in the mature primary visual cortex. However, some degree of short-term adaptability could be mediated by changes in responsiveness.

Molecular cloning and sequence analysis of expansins--a highly conserved, multigene family of proteins that mediate cell wall extension in plants.

Expansins are unusual proteins discovered by virtue of their ability to mediate cell wall extension in plants. We identified cDNA clones for two cucumber expansins on the basis of peptide sequences of proteins purified from cucumber hypocotyls. The expansin cDNAs encode related proteins with signal peptides predicted to direct protein secretion to the cell wall. Northern blot analysis showed moderate transcript abundance in the growing region of the hypocotyl and no detectable transcripts in the nongrowing region. Rice and Arabidopsis expansin cDNAs were identified from collections of anonymous cDNAs (expressed sequence tags). Sequence comparisons indicate at least four distinct expansin cDNAs in rice and at least six in Arabidopsis. Expansins are highly conserved in size and sequence (60-87% amino acid sequence identity and 75-95% similarity between any pairwise comparison), and phylogenetic trees indicate that this multigene family formed before the evolutionary divergence of monocotyledons and dicotyledons. Sequence and motif analyses show no similarities to known functional domains that might account for expansin action on wall extension. A series of highly conserved tryptophans may function in expansin binding to cellulose or other glycans. The high conservation of this multigene family indicates that the mechanism by which expansins promote wall extensin tolerates little variation in protein structure.

Is water structure around hydrophobic groups clathrate-like?

The term "clathrate structure" is quantified for solvation of nonpolar groups by enumerating hydrogen-bonded ring sizes both in the solvation shell and through the shell-bulk interface and comparing it to a bulk control using the ST4 water model. For clathrate-like structure to be evident, the distributions along the hydrophobic surface are expected to be dominated by pentagons, with significant depletion of hexagons and larger polygons. While the distribution in this region is indeed distinguished by a large number of pentagons, there are significant contributions from hexagons and larger rings as well. Calculated polygon distributions through the shell-bulk interface indicate that when water structure is highly cooperative along the hydrophobic surface, hydrogen-bonded pathways leading back into bulk are then reduced. These results are qualitatively consistent with the observation that hydrophobicity is proportional to the nonpolar solute surface area.

Three-dimensional structure of a cysteine-rich repeat from the low-density lipoprotein receptor.

The low-density lipoprotein (LDL) receptor plays a central role in mammalian cholesterol metabolism, clearing lipoproteins which bear apolipoproteins E and B-100 from plasma. Mutations in this molecule are associated with familial hypercholesterolemia, a condition which leads to an elevated plasma cholesterol concentration and accelerated atherosclerosis. The N-terminal segment of the LDL receptor contains a heptad of cysteine-rich repeats that bind the lipoproteins. Similar repeats are present in related receptors, including the very low-density lipoprotein receptor and the LDL receptor-related protein/alpha 2-macroglobulin receptor, and in proteins which are functionally unrelated, such as the C9 component of complement. The first repeat of the human LDL receptor has been expressed in Escherichia coli as a glutathione S-transferase fusion protein, and the cleaved and purified receptor module has been shown to fold to a single, fully oxidized form that is recognized by the monoclonal antibody IgG-C7 in the presence of calcium ions. The three-dimensional structure of this module has been determined by two-dimensional NMR spectroscopy and shown to consist of a beta-hairpin structure, followed by a series of beta turns. Many of the side chains of the acidic residues, including the highly conserved Ser-Asp-Glu triad, are clustered on one face of the module. To our knowledge, this structure has not previously been described in any other protein and may represent a structural paradigm both for the other modules in the LDL receptor and for the homologous domains of several other proteins. Calcium ions had only minor effects on the CD spectrum and no effect on the 1H NMR spectrum of the repeat, suggesting that they induce no significant conformational change.

Structure and stability of a second molten globule intermediate in the apomyoglobin folding pathway.

Apomyoglobin folding proceeds through a molten globule intermediate (low-salt form; I1) that has been characterized by equilibrium (pH 4) and kinetic (pH 6) folding experiments. Of the eight alpha-helices in myoglobin, three (A, G, and H) are structured in I1, while the rest appear to be unfolded. Here we report on the structure and stability of a second intermediate, the trichloroacetate form of the molten globule intermediate (I2), which is induced either from the acid-unfolded protein or from I1 by > or = 5 mM sodium trichloroacetate. Circular dichroism measurements monitoring urea- and acid-induced unfolding indicate that I2 is more highly structured and more stable than I1. Although I2 exhibits properties closer to those of the native protein, one-dimensional NMR spectra show that it maintains the lack of fixed side-chain structure that is the hallmark of a molten globule. Amide proton exchange and 1H-15N two-dimensional NMR experiments are used to identify the source of the extra helicity observed in I2. The results reveal that the existing A, G, and H helices present in I1 have become more stable in I2 and that a fourth helix--the B helix--has been incorporated into the molten globule. Available evidence is consistent with I2 being an on-pathway intermediate. The data support the view that apomyoglobin folds in a sequential fashion through a single pathway populated by intermediates of increasing structure and stability.

Folding of a nascent polypeptide chain in vitro: cooperative formation of structure in a protein module.

We have prepared a family of peptide fragments of the 64-residue chymotrypsin inhibitor 2, corresponding to its progressive elongation from the N terminus. The growing polypeptide chain has little tendency to form stable structure until it is largely synthesized, and what structures are formed are nonnative and lack, in particular, the native secondary structural elements of alpha-helix and beta-sheet. These elements then develop as sufficient tertiary interactions are made in the nearly full-length chain. The growth of structure in the small module is highly cooperative and does not result from the hierarchical accretion of substructures.