Nonequilibrium response spectroscopy of voltage-sensitive ion channel gating.

We describe a new electrophysiological technique called nonequilibrium response spectroscopy, which involves application of rapidly fluctuating (as high as 14 kHz) large-amplitude voltage clamp waveforms to ion channels. As a consequence of the irreversible (in the sense of Carnot) exchange of energy between the fluctuating field and the channel protein, the gating response is exquisitely sensitive to features of the kinetics that are difficult or impossible to adequately resolve by means of traditional stepped potential protocols. Here we focus on the application of dichotomous (telegraph) noise voltage fluctuations, a broadband Markovian colored noise that fluctuates between two values. Because Markov kinetic models of channel gating can be embedded within higher-dimensional Markov models that take into account the effects of the voltage fluctuations, many features of the response of the channels can be calculated algebraically. This makes dichotomous noise and its generalizations uniquely suitable for model selection and kinetic analysis. Although we describe its application to macroscopic ionic current measurements, the nonequilibrium response method can also be applied to gating and single channel current recording techniques. We show how data from the human cardiac isoform (hH1a) of the Na+ channel expressed in mammalian cells can be acquired and analyzed, and how these data reveal hidden aspects of the molecular kinetics that are not revealed by conventional methods.     

Molecular cloning and analysis of zebrafish voltage-gated sodium channel beta subunit genes: implications for the evolution of electrical signaling in vertebrates

Background  

Action potential generation in excitable cells such as myocytes and neurons critically depends on voltage-gated sodium channels. In mammals, sodium channels exist as macromolecular complexes that include a pore-forming alpha subunit and 1 or more modulatory beta subunits. Although alpha subunit genes have been cloned from diverse metazoans including flies, jellyfish, and humans, beta subunits have not previously been identified in any non-mammalian species. To gain further insight into the evolution of electrical signaling in vertebrates, we investigated beta subunit genes in the teleost Danio rerio (zebrafish).     

Molecular evolution and modern human origins

While molecular evolutionists may be fascinated by the features and history of a particular gene or DNA segment, evolutionary anthropologists are often more interested in the activities and history of groups of people. We may want to know, for instance, when and where humans have migrated, how much exchange between groups has taken place, and how population sizes have changed. Population genetic theory provides the hope that through analyses of genetic data we will gain insight into the history of populations. Genetic data from extant human populations are now accruing at a remarkable rate. We might, therefore, expect to have answers in hand. There remains, however, a wide gap between the available theory and data; too often we fail to draw firm conclusions because our interpretation of analytic results requires that we make myriad assumptions about our data. In any one instance, these assumptions might include estimates of mutation rate, mutational mechanism, population sizes, the role that natural selection has played, and the rate of migration among groups. Often these assumptions are implicit, invisible to most. How, then, are we to make any progress? © 1998 Wiley-Liss, Inc.     

Molecular evolution of ubiquitin genes

Ubiquitin is a singular protein with multiple functions. It is probably the most slowly evolving protein known, is encoded by genes with a unique structure, and provides an intriguing case study for various aspects of molecular evolution. In particular, the multiple ubiquitin-coding repeats which have been characterized in man, yeast and a slime mould graphically illustrate the dynamics of concerted evolution, but cast doubts on the effectiveness of this process for unlinked arrays in this repeat family.     

Molecular phylogenies of plastid origins and algal evolution

Summary An overview of recent molecular analyses regarding origins of plastids in algal lineages is presented. Since different phylogenetic analyses can yield contradictory views of algal plastid origins, we have examined the effect of two distance measurement methods and two distance matrix tree-building methods upon topologies for the ribulose-1,5-bisphosphate carboxylase/oxygenase large subunit nucleotide sequence data set. These results are contrasted to those from bootstrap parsimony analysis of nucleotide sequence data subsets. It is shown that the phylogenetic information contained within nucleotide sequences for the chloroplast-encoded gene for the large subunit of ribulose-1,5-bisphosphate carboxylase/oxygenase, integral to photosynthesis, indicates an independent origin for this plastid gene in different plant taxa. This finding is contrasted to contrary results derived from 16S rRNA sequences. Possible explanations for discrepancies observed for these two different molecules are put forth. Other molecular sequence data which address questions of early plant evolution and the eubacterial origins of algal organelles are discussed. Offprint requests to: W. Martin     

Molecular evolution of the murine tspy genes

Molecular aspects of murine evolution were studied by sequencing, and subsequently comparing, introns of the Y-chromosomal tspy genes from Apodemus agrarius, A. sylvaticus, A. flavicollis, Mus platythrix (subgenus Pyromys), M. booduga (subgenus Leggada), and from species of the subgenus Mus, including M. cervicolor, M. macedonicus and M. spretus. Estimates of nucleotide substitution rates in these lineages were in perfect agreement with phylogenetic data previously published by She et al. (1990), Catzeflis et al. (1992; 1993), and Lyon et al. (1996). The only exception was provided by a comparatively late divergence of M. spretus and M. macedonicus. Our data also suggest that M. booduga diverged from the subgenus Mus about 3 Myr ago.     

Molecular evolution of nitrate reductase genes

To understand the evolutionary mechanisms and relationships of nitrate reductases (NRs), the nucleotide sequences encoding 19 nitrate reductase (NR) genes from 16 species of fungi, algae, and higher plants were analyzed. The NR genes examined show substantial sequence similarity, particularly within functional domains, and large variations in GC content at the third codon position and intron number. The intron positions were different between the fungi and plants, but conserved within these groups. The overall and nonsynonymous substitution rates among fungi, algae, and higher plants were estimated to be 4.33 × 10⁻¹⁰ and 3.29 × 10⁻¹⁰ substitutions per site per year. The three functional domains of NR genes evolved at about one-third of the rate of the N-terminal and the two hinge regions connecting the functional domains. Relative rate tests suggested that the nonsynonymous substitution rates were constant among different lineages, while the overall nucleotide substitution rates varied between some lineages. The phylogenetic trees based on NR genes correspond well with the phylogeny of the organisms determined from systematics and other molecular studies. Based on the nonsynonymous substitution rate, the divergence time of monocots and dicots was estimated to be about 340 Myr when the fungi–plant or algae–higher plant divergence times were used as reference points and 191 Myr when the rice–barley divergence time was used as a reference point. These two estimates are consistent with other estimates of divergence times based on these reference points. The lack of consistency between these two values appears to be due to the uncertainty of the reference times. Received: 10 April 1995 / Accepted: 10 September 1995     

Molecular origins of asymmetric proton conduction in the influenza M2 channel

The M2 proton channel of influenza A is embedded into the viral envelope and allows acidification of the virion when the external pH is lowered. In contrast, no outward proton conductance is observed when the internal pH is lowered, although outward current is observed at positive voltage. Residues Trp41 and Asp44 are known to play a role in preventing pH-driven outward conductance, but the mechanism for this is unclear. We investigate this issue using classical molecular dynamics simulations with periodic proton hops. When all key His37 residues are neutral, inward proton movement is much more facile than outward movement if the His are allowed to shuttle the proton. The preference for inward movement increases further as the charge on the His37 increases. Analysis of the trajectories reveals three factors accounting for this asymmetry. First, in the outward direction, Asp44 traps the hydronium by strong electrostatic interactions. Secondly, Asp44 and Trp41 orient the hydronium with the protons pointing inward, hampering outward Grotthus hopping. As a result, the effective barrier is lower in the inward direction. Trp41 adds to the barrier by weakly H-bonding to potential H⁺ acceptors. Finally, for charged His, the H₃O⁺ in the inner vestibule tends to get trapped at lipid-lined fenestrations of the cone-shaped channel. Simulations qualitatively reproduce the experimentally observed higher outward conductance of mutants. The ability of positive voltage, unlike proton gradient, to induce an outward current appears to arise from its ability to bias H₃O⁺ and the waters around it toward more H-outward orientations.     

Molecular considerations in the evolution of bacterial genes

Summary Synonymous and nonsynonymous substitution rates at the loci encoding glyceraldehyde-3-phosphate dehydrogenase (gap) and outer membrane protein 3A (ompA) were examined in 12 species of enteric bacteria. By examining homologous sequences in species of varying degrees of relatedness and of known phylogenetic relationships, we analyzed the patterns of synonymous and nonsynonymous substitutions within and among these genes. Although both loci accumulate synonymous substitutions at reduced rates due to codon usage bias, portions of thegap andompA reading frames show significant deviation in synonymous substitution rates not attributable to local codon bias. A paucity of synonymous substitutions in portions of theompA gene may reflect selection for a novel mRNA secondary structure. In addition, these studies allow comparisons of homologous protein-coding sequences (gap) in plants, animals, and bacteria, revealing differences in evolutionary constraints on this glycolytic enzyme in these lineages.     

Molecular regulation of cardiac ryanodine receptor ion channel

The release of Ca(2+) ions from intracellular stores is a key step in a wide variety of cellular functions. In striated muscle, the release of Ca(2+) from the sarcoplasmic reticulum (SR) leads to muscle contraction. Ca(2+) release occurs through large, high-conductance Ca(2+) release channels, also known as ryanodine receptors (RyRs) because they bind the plant alkaloid ryanodine with high affinity and specificity. The RyRs are isolated as 30S protein complexes comprised of four 560 kDa RyR2 subunits and four 12 kDa FK506 binding protein (FKBP12) subunits. Multiple endogenous effector molecules and posttranslational modifications regulate the RyRs. This review focuses on current research toward understanding the control of the isolated cardiac Ca(2+) release channel/ryanodine receptor (RyR2) by Ca(2+), calmodulin, thiol oxidation/reduction and nitrosylation, and protein phosphorylation.     

Molecular genetics of the VDAC ion channel: Structural model and sequence analysis

The voltage-dependent anion-selective channel of the outer mitochondrial membrane provides a unique system in which to study the molecular basis of voltage gating of ion flow. We have cloned and sequenced acDNA coding for this protein in yeast. From the derived amino acid sequence, we have generated a preliminary model for the secondary structure of the protein which suggests that the protein forms a -barrel type structure. Comparison of the VDAC amino acid sequence with that of the bacterial porins has indicated that the two classes of molecules appear to be unrelated.     

Molecular dynamics simulation of the antiamoebin ion channel: linking structure and conductance

Molecular-dynamics simulations were carried out to ascertain which of the potential multimeric forms of the transmembrane peptaibol channel, antiamoebin, is consistent with its measured conductance. Estimates of the conductance obtained through counting ions that cross the channel and by solving the Nernst-Planck equation yield consistent results, indicating that the motion of ions inside the channel can be satisfactorily described as diffusive. The calculated conductance of octameric channels is markedly higher than the conductance measured in single channel recordings, whereas the tetramer appears to be nonconducting. The conductance of the hexamer was estimated to be 115 ± 34 pS and 74 ± 20 pS, at 150 mV and 75 mV, respectively, in satisfactory agreement with the value of 90 pS measured at 75 mV. On this basis, we propose that the antiamoebin channel consists of six monomers. Its pore is large enough to accommodate K⁺ and Cl⁻ with their first solvation shells intact. The free energy barrier encountered by K⁺ is only 2.2 kcal/mol whereas Cl⁻ encounters a substantially higher barrier of nearly 5 kcal/mol. This difference makes the channel selective for cations. Ion crossing events are shown to be uncorrelated and follow Poisson statistics.     

Molecular evolution of cycloidea-like genes in Fabaceae

The cycloidea (CYC) gene controls floral symmetry in snapdragon (Antirrhinum majus). We investigated the evolution of CYC-like genes in some species of legumes that have zygomorphic flowers. Two to four CYC-like genes were isolated from a single species. The results of NJ and ML analyses indicate that CYC-like genes in legumes group into two monophyletic clades; one group consists of eight CYC-like genes (Clade 1) and the other contains three CYC-like genes and TB1 of maize (Clade 2). These phylogenetic trees and the Shimodaira–Hasegawa test suggest that Clade 1 is a sister of the original CYC group (Clade 3). Moreover, the result of the GeneTree analysis showed that the CYC-like genes experienced repeated duplication events during the evolution of legumes. We herein speculate as to the role of CYC-like genes in legumes and discuss the evolutionary processes that these genes have undergone. Current address (Jun Yokoyama and Masayuki Maki): Division of Ecology and Evolutionary Biology, Graduate School of Life Sciences, Tohoku University, Aoba, Sendai 980-8578, Japan     

Molecular evolution of genes in avian genomes

Background

Obtaining a draft genome sequence of the zebra finch (Taeniopygia guttata), the second bird genome to be sequenced, provides the necessary resource for whole-genome comparative analysis of gene sequence evolution in a non-mammalian vertebrate lineage. To analyze basic molecular evolutionary processes during avian evolution, and to contrast these with the situation in mammals, we aligned the protein-coding sequences of 8,384 1:1 orthologs of chicken, zebra finch, a lizard and three mammalian species.

Results

We found clear differences in the substitution rate at fourfold degenerate sites, being lowest in the ancestral bird lineage, intermediate in the chicken lineage and highest in the zebra finch lineage, possibly reflecting differences in generation time. We identified positively selected and/or rapidly evolving genes in avian lineages and found an over-representation of several functional classes, including anion transporter activity, calcium ion binding, cell adhesion and microtubule cytoskeleton.

Conclusions

Focusing specifically on genes of neurological interest and genes differentially expressed in the unique vocal control nuclei of the songbird brain, we find a number of positively selected genes, including synaptic receptors. We found no evidence that selection for beneficial alleles is more efficient in regions of high recombination; in fact, there was a weak yet significant negative correlation between ω and recombination rate, which is in the direction predicted by the Hill-Robertson effect if slightly deleterious mutations contribute to protein evolution. These findings set the stage for studies of functional genetics of avian genes.     

Molecular evolution of Drosophila odorant receptor genes

A total of 752 odorant receptor (Or) genes, including pseudogenes, were identified in 11 Drosophila species and named after their orthologs in Drosophila melanogaster. The 813 Or genes, including 61 from D. melanogaster, were classified into 59 orthologous groups that are well supported by gene phylogeny. By reconciling with the gene family phylogeny, we estimated the number of gene duplication/loss events and intron gain/loss events in the species phylogeny. We found that these events are particularly frequent in Drosophila grimshawi, Drosophila willistoni, and obscura group. More than half of the duplicated genes stay as tandem arrays, whose size range from 2 to 8. These genes vary in sequence and some likely underwent positive selection, indicating that the gene duplication was important for flies to acquire new olfactory functions. We hypothesize that Or genes conferred the basic olfactory repertoire to ancestral flies before the speciation of the Drosophila and Sophophora subgenera about 40 Mya. This repertoire has been largely maintained in the current species, whereas lineage-specific gene duplication seems to have led to additional specialization in some species in response to specific ecological conditions.     

Molecular evolution of human visual pigment genes

By comparing the published DNA sequences for (a) the genes encoding the human visual color pigments (red, green, and blue) with (b) the genes encoding human, bovine, and Drosophila rhodopsins, a phylogenetic tree for the mammalian pigment genes has been constructed. This evolutionary tree shows that the common ancestor of the visual color pigment genes diverged first from that of the rhodopsin genes; then the common ancestor of the red and green pigment genes and the ancestor of the blue pigment gene diverged; and finally the red and green pigment genes diverged from each other much more recently. Nucleotide substitutions in the rhodopsin genes are best explained by the neutral theory of molecular evolution. However, important functional adaptations seem to have occurred twice during the evolution of the color pigment genes in humans: first, to the common ancestor of the three color pigment genes after its divergence from the ancestor of the rhodopsin gene and, second, to the ancestor of the red pigment gene after its divergence from that of the green pigment gene.