A conformational switch driven by phosphorylation regulates the activity of the evolutionarily conserved SNARE Ykt6.

Ykt6 is a soluble N-ethylmaleimide sensitive factor activating protein receptor (SNARE) critically involved in diverse vesicular fusion pathways. While most SNAREs rely on transmembrane domains for their activity, Ykt6 dynamically cycles between the cytosol and membrane-bound compartments where it is active. The mechanism that regulates these transitions and allows Ykt6 to achieve specificity toward vesicular pathways is unknown. Using a Parkinson's disease (PD) model, we found that Ykt6 is phosphorylated at an evolutionarily conserved site which is regulated by Ca2+ signaling. Through a multidisciplinary approach, we show that phosphorylation triggers a conformational change that allows Ykt6 to switch from a closed cytosolic to an open membrane-bound form. In the phosphorylated open form, the spectrum of protein interactions changes, leading to defects in both the secretory and autophagy pathways, enhancing toxicity in PD models. Our studies reveal a mechanism by which Ykt6 conformation and activity are regulated with potential implications for PD.

Computational 3D histological phenotyping of whole zebrafish by X-ray histotomography.

Organismal phenotypes frequently involve multiple organ systems. Histology is a powerful way to detect cellular and tissue phenotypes, but is largely descriptive and subjective. To determine how synchrotron-based X-ray micro-tomography (micro-CT) can yield 3-dimensional whole-organism images suitable for quantitative histological phenotyping, we scanned whole zebrafish, a small vertebrate model with diverse tissues, at ~1 micron voxel resolutions. Micro-CT optimized for cellular characterization (histotomography) allows brain nuclei to be computationally segmented and assigned to brain regions, and cell shapes and volumes to be computed for motor neurons and red blood cells. Striking individual phenotypic variation was apparent from color maps of computed densities of brain nuclei. Unlike histology, the histotomography also allows the study of 3-dimensional structures of millimeter scale that cross multiple tissue planes. We expect the computational and visual insights into 3D cell and tissue architecture provided by histotomography to be useful for reference atlases, hypothesis generation, comprehensive organismal screens, and diagnostics.

Rigid Embedding of Fixed and Stained, Whole, Millimeter-Scale Specimens for Section-free 3D Histology by Micro-Computed Tomography.

For over a hundred years, the histological study of tissues has been the gold standard for medical diagnosis because histology allows all cell types in every tissue to be identified and characterized. Our laboratory is actively working to make technological advances in X-ray micro-computed tomography (micro-CT) that will bring the diagnostic power of histology to the study of full tissue volumes at cellular resolution (i.e., an X-ray Histo-tomography modality). Toward this end, we have made targeted improvements to the sample preparation pipeline. One key optimization, and the focus of the present work, is a straightforward method for rigid embedding of fixed and stained millimeter-scale samples. Many of the published methods for sample immobilization and correlative micro-CT imaging rely on placing the samples in paraffin wax, agarose, or liquids such as alcohol. Our approach extends this work with custom procedures and the design of a 3-dimensional printable apparatus to embed the samples in an acrylic resin directly into polyimide tubing, which is relatively transparent to X-rays. Herein, sample preparation procedures are described for the samples from 0.5 to 10 mm in diameter, which would be suitable for whole zebrafish larvae and juveniles, or other animals and tissue samples of similar dimensions. As proof of concept, we have embedded the specimens from Danio, Drosophila, Daphnia, and a mouse embryo; representative images from 3-dimensional scans for three of these samples are shown. Importantly, our methodology leads to multiple benefits including rigid immobilization, long-term preservation of laboriously-created resources, and the ability to re-interrogate samples.

The S6 gate in regulatory Kv6 subunits restricts heteromeric K+ channel stoichiometry.

The Shaker-like family of voltage-gated K+ channels comprises four functionally independent gene subfamilies, Shaker (Kv1), Shab (Kv2), Shaw (Kv3), and Shal (Kv4), each of which regulates distinct aspects of neuronal excitability. Subfamily-specific assembly of tetrameric channels is mediated by the N-terminal T1 domain and segregates Kv1-4, allowing multiple channel types to function independently in the same cell. Typical Shaker-like Kv subunits can form functional channels as homotetramers, but a group of mammalian Kv2-related genes (Kv5.1, Kv6s, Kv8s, and Kv9s) encodes subunits that have a "silent" or "regulatory" phenotype characterized by T1 self-incompatibility. These channels are unable to form homotetramers, but instead heteromerize with Kv2.1 or Kv2.2 to diversify the functional properties of these delayed rectifiers. While T1 self-incompatibility predicts that these heterotetramers could contain up to two regulatory (R) subunits, experiments show a predominance of 3:1R stoichiometry in which heteromeric channels contain a single regulatory subunit. Substitution of the self-compatible Kv2.1 T1 domain into the regulatory subunit Kv6.4 does not alter the stoichiometry of Kv2.1:Kv6.4 heteromers. Here, to identify other channel structures that might be responsible for favoring the 3:1R stoichiometry, we compare the sequences of mammalian regulatory subunits to independently evolved regulatory subunits from cnidarians. The most widespread feature of regulatory subunits is the presence of atypical substitutions in the highly conserved consensus sequence of the intracellular S6 activation gate of the pore. We show that two amino acid substitutions in the S6 gate of the regulatory subunit Kv6.4 restrict the functional stoichiometry of Kv2.1:Kv6.4 to 3:1R by limiting the formation and function of 2:2R heteromers. We propose a two-step model for the evolution of the asymmetric 3:1R stoichiometry, which begins with evolution of self-incompatibility to establish the regulatory phenotype, followed by drift of the activation gate consensus sequence under relaxed selection to limit stoichiometry to 3:1R.

Comparative analysis of fixation and embedding techniques for optimized histological preparation of zebrafish.

In recognition of the importance of zebrafish as a model organism for studying human disease, we have created zebrafish content for a web-based reference atlas of microanatomy for comparing histology and histopathology between model systems and with humans (http://bio-atlas.psu.edu). Fixation, decalcification, embedding, and sectioning of zebrafish were optimized to maximize section quality. A comparison of protocols involving six fixatives showed that 10% Neutral Buffered Formalin at 21°C for 24h yielded excellent results. Sectioning of juveniles and adults requires bone decalcification; EDTA at 0.35M produced effective decalcification in 21-day-old juveniles through adults (≥~3Months). To improve section plane consistency in sets of larvae, we have developed new array casting molds based on the outside contours of larvae derived from 3D microCT images. Tissue discontinuity in sections, a common barrier to creating quality sections of zebrafish, was minimized by processing and embedding the formalin-fixed zebrafish tissues in plasticized forms of paraffin wax, and by periodic hydration of the block surface in ice water between sets of sections. Optimal H&E (Hematoxylin and Eosin) staining was achieved through refinement of standard protocols. High quality slide scans produced from glass histology slides were digitally processed to maximize image quality, and experimental replicates posted as full slides as part of this publication. Modifications to tissue processing are still needed to eliminate the need for block surface hydration. The further addition of slide collections from other model systems and 3D tools for visualizing tissue architecture would greatly increase the utility of the digital atlas.

FKBP12 contributes to α-synuclein toxicity by regulating the calcineurin-dependent phosphoproteome.

Calcineurin is an essential Ca2+-dependent phosphatase. Increased calcineurin activity is associated with α-synuclein (α-syn) toxicity, a protein implicated in Parkinson's Disease (PD) and other neurodegenerative diseases. Calcineurin can be inhibited with Tacrolimus through the recruitment and inhibition of the 12-kDa cis-trans proline isomerase FK506-binding protein (FKBP12). Whether calcineurin/FKBP12 represents a native physiologically relevant assembly that occurs in the absence of pharmacological perturbation has remained elusive. We leveraged α-syn as a model to interrogate whether FKBP12 plays a role in regulating calcineurin activity in the absence of Tacrolimus. We show that FKBP12 profoundly affects the calcineurin-dependent phosphoproteome, promoting the dephosphorylation of a subset of proteins that contributes to α-syn toxicity. Using a rat model of PD, partial elimination of the functional interaction between FKBP12 and calcineurin, with low doses of the Food and Drug Administration (FDA)-approved compound Tacrolimus, blocks calcineurin's activity toward those proteins and protects against the toxic hallmarks of α-syn pathology. Thus, FKBP12 can endogenously regulate calcineurin activity with therapeutic implications for the treatment of PD.

Synchrotron microCT imaging of soft tissue in juvenile zebrafish reveals retinotectal projections.

Biomedical research and clinical diagnosis would benefit greatly from full volume determinations of anatomical phenotype. Comprehensive tools for morphological phenotyping are central for the emerging field of phenomics, which requires high-throughput, systematic, accurate, and reproducible data collection from organisms affected by genetic, disease, or environmental variables. Theoretically, complete anatomical phenotyping requires the assessment of every cell type in the whole organism, but this ideal is presently untenable due to the lack of an unbiased 3D imaging method that allows histopathological assessment of any cell type despite optical opacity. Histopathology, the current clinical standard for diagnostic phenotyping, involves the microscopic study of tissue sections to assess qualitative aspects of tissue architecture, disease mechanisms, and physiological state. However, quantitative features of tissue architecture such as cellular composition and cell counting in tissue volumes can only be approximated due to characteristics of tissue sectioning, including incomplete sampling and the constraints of 2D imaging of 5 micron thick tissue slabs. We have used a small, vertebrate organism, the zebrafish, to test the potential of microCT for systematic macroscopic and microscopic morphological phenotyping. While cell resolution is routinely achieved using methods such as light sheet fluorescence microscopy and optical tomography, these methods do not provide the pancellular perspective characteristic of histology, and are constrained by the limited penetration of visible light through pigmented and opaque specimens, as characterizes zebrafish juveniles. Here, we provide an example of neuroanatomy that can be studied by microCT of stained soft tissue at 1.43 micron isotropic voxel resolution. We conclude that synchrotron microCT is a form of 3D imaging that may potentially be adopted towards more reproducible, large-scale, morphological phenotyping of optically opaque tissues. Further development of soft tissue microCT, visualization and quantitative tools will enhance its utility.

Bilaterian Giant Ankyrins Have a Common Evolutionary Origin and Play a Conserved Role in Patterning the Axon Initial Segment.

In vertebrate neurons, the axon initial segment (AIS) is specialized for action potential initiation. It is organized by a giant 480 Kd variant of ankyrin G (AnkG) that serves as an anchor for ion channels and is required for a plasma membrane diffusion barrier that excludes somatodendritic proteins from the axon. An unusually long exon required to encode this 480Kd variant is thought to have been inserted only recently during vertebrate evolution, so the giant ankyrin-based AIS scaffold has been viewed as a vertebrate adaptation for fast, precise signaling. We re-examined AIS evolution through phylogenomic analysis of ankyrins and by testing the role of ankyrins in proximal axon organization in a model multipolar Drosophila neuron (ddaE). We find giant isoforms of ankyrin in all major bilaterian phyla, and present evidence in favor of a single common origin for giant ankyrins and the corresponding long exon in a bilaterian ancestor. This finding raises the question of whether giant ankyrin isoforms play a conserved role in AIS organization throughout the Bilateria. We examined this possibility by looking for conserved ankyrin-dependent AIS features in Drosophila ddaE neurons via live imaging. We found that ddaE neurons have an axonal diffusion barrier proximal to the cell body that requires a giant isoform of the neuronal ankyrin Ank2. Furthermore, the potassium channel shal concentrates in the proximal axon in an Ank2-dependent manner. Our results indicate that the giant ankyrin-based cytoskeleton of the AIS may have evolved prior to the radiation of extant bilaterian lineages, much earlier than previously thought.

Novel Molecular Interactions of Acylcarnitines and Fatty Acids with Myoglobin.

Previous research has indicated that long-chain fatty acids can bind myoglobin (Mb) in an oxygen-dependent manner. This suggests that oxy-Mb may play an important role in fuel delivery in Mb-rich muscle fibers (e.g. type I fibers and cardiomyocytes), and raises the possibility that Mb also serves as an acylcarnitine-binding protein. We report for the first time the putative interaction and affinity characteristics for different chain lengths of both fatty acids and acylcarnitines with oxy-Mb using molecular dynamic simulations and isothermal titration calorimetry experiments. We found that short- to medium-chain fatty acids or acylcarnitines (ranging from C2:0 to C10:0) fail to achieve a stable conformation with oxy-Mb. Furthermore, our results indicate that C12:0 is the minimum chain length essential for stable binding of either fatty acids or acylcarnitines with oxy-Mb. Importantly, the empirical lipid binding studies were consistent with structural modeling. These results reveal that: (i) the lipid binding affinity for oxy-Mb increases as the chain length increases (i.e. C12:0 to C18:1), (ii) the binding affinities of acylcarnitines are higher when compared with their respective fatty acid counterparts, and (iii) both fatty acids and acylcarnitines bind to oxy-Mb in 1:1 stoichiometry. Taken together, our results support a model in which oxy-Mb is a novel regulator of long-chain acylcarnitine and fatty acid pools in Mb-rich tissues. This has important implications for physiological fuel management during exercise, and relevance to pathophysiological conditions (e.g. fatty acid oxidation disorders and cardiac ischemia) where long-chain acylcarnitine accumulation is evident.

Identification of Respiratory Syncytial Virus Nonstructural Protein 2 Residues Essential for Exploitation of the Host Ubiquitin System and Inhibition of Innate Immune Responses.

UNLABELLED: Respiratory syncytial virus (RSV) is a leading cause of lower respiratory tract infection in young children worldwide. The RSV nonstructural protein 2 (NS2) is a multifunctional protein that primarily acts to antagonize the innate immune system by targeting STAT2 for proteasomal degradation. We investigated the structural determinants of NS2 important for interaction with the host ubiquitin system to degrade STAT2 during infection. We found that NS2 expression enhances ubiquitination of host proteins. Bioinformatics analysis provided a platform for identification of specific residues that limit NS2-induced ubiquitination. Combinations of multiple mutations displayed an additive effect on reducing NS2-induced ubiquitination. Using a reverse genetics system, we generated recombinant RSV (rRSV) containing NS2 ubiquitin mutations, which maintained their effect on ubiquitin expression during infection. Interestingly, STAT2 degradation activity was ablated in the NS2 ubiquitin mutant rRSV. In addition, NS2 ubiquitin mutations decreased rRSV replication, indicating a correlation between NS2's ubiquitin function and antagonism of innate immune signaling to enhance viral replication. Our approach of targeting NS2 residues required for NS2 inhibition of immune responses provides a mechanism for attenuating RSV for vaccine development. IMPORTANCE: RSV has been circulating globally for more than 60 years, causing severe respiratory disease in pediatric, elderly, and immunocompromised populations. Production of a safe, effective vaccine against RSV is a public health priority. The NS2 protein is an effective target for prevention and treatment of RSV due to its antagonistic activity against the innate immune system. However, NS2-deleted RSV vaccine candidates rendered RSV overattenuated or poorly immunogenic. Alternatively, we can modify essential NS2 structural features to marginally limit viral growth while maintaining immune responses, providing the necessary balance between antigenicity and safety required for an effective vaccine. We coupled bioinformatics analysis with reverse genetics to introduce mutations into RSV's negative-sense genome. In this way we constructed rRSV NS2 ubiquitin mutants that limited NS2's ability to antagonize the innate immune system, thereby attenuating rRSV growth and increasing innate immune responses.

Functional Characterization of Cnidarian HCN Channels Points to an Early Evolution of Ih.

HCN channels play a unique role in bilaterian physiology as the only hyperpolarization-gated cation channels. Their voltage-gating is regulated by cyclic nucleotides and phosphatidylinositol 4,5-bisphosphate (PIP2). Activation of HCN channels provides the depolarizing current in response to hyperpolarization that is critical for intrinsic rhythmicity in neurons and the sinoatrial node. Additionally, HCN channels regulate dendritic excitability in a wide variety of neurons. Little is known about the early functional evolution of HCN channels, but the presence of HCN sequences in basal metazoan phyla and choanoflagellates, a protozoan sister group to the metazoans, indicate that the gene family predates metazoan emergence. We functionally characterized two HCN channel orthologs from Nematostella vectensis (Cnidaria, Anthozoa) to determine which properties of HCN channels were established prior to the emergence of bilaterians. We find Nematostella HCN channels share all the major functional features of bilaterian HCNs, including reversed voltage-dependence, activation by cAMP and PIP2, and block by extracellular Cs+. Thus bilaterian-like HCN channels were already present in the common parahoxozoan ancestor of bilaterians and cnidarians, at a time when the functional diversity of voltage-gated K+ channels was rapidly expanding. NvHCN1 and NvHCN2 are expressed broadly in planulae and in both the endoderm and ectoderm of juvenile polyps.

Bimodal regulation of an Elk subfamily K+ channel by phosphatidylinositol 4,5-bisphosphate.

Phosphatidylinositol 4,5-bisphosphate (PIP2) regulates Shaker K+ channels and voltage-gated Ca2+ channels in a bimodal fashion by inhibiting voltage activation while stabilizing open channels. Bimodal regulation is conserved in hyperpolarization-activated cyclic nucleotide-gated (HCN) channels, but voltage activation is enhanced while the open channel state is destabilized. The proposed sites of PIP2 regulation in these channels include the voltage-sensor domain (VSD) and conserved regions of the proximal cytoplasmic C terminus. Relatively little is known about PIP2 regulation of Ether-á-go-go (EAG) channels, a metazoan-specific family of K+ channels that includes three gene subfamilies, Eag (Kv10), Erg (Kv11), and Elk (Kv12). We examined PIP2 regulation of the Elk subfamily potassium channel human Elk1 to determine whether bimodal regulation is conserved within the EAG K+ channel family. Open-state stabilization by PIP2 has been observed in human Erg1, but the proposed site of regulation in the distal C terminus is not conserved among EAG family channels. We show that PIP2 strongly inhibits voltage activation of Elk1 but also stabilizes the open state. This stabilization produces slow deactivation and a mode shift in voltage gating after activation. However, removal of PIP2 has the net effect of enhancing Elk1 activation. R347 in the linker between the VSD and pore (S4-S5 linker) and R479 near the S6 activation gate are required for PIP2 to inhibit voltage activation. The ability of PIP2 to stabilize the open state also requires these residues, suggesting an overlap in sites central to the opposing effects of PIP2 on channel gating. Open-state stabilization in Elk1 requires the N-terminal eag domain (PAS domain + Cap), and PIP2-dependent stabilization is enhanced by a conserved basic residue (K5) in the Cap. Our data shows that PIP2 can bimodally regulate voltage gating in EAG family channels, as has been proposed for Shaker and HCN channels. PIP2 regulation appears fundamentally different for Elk and KCNQ channels, suggesting that, although both channel types can regulate action potential threshold in neurons, they are not functionally redundant.

Molecular dynamic simulations reveal the structural determinants of Fatty Acid binding to oxy-myoglobin.

The mechanism(s) by which fatty acids are sequestered and transported in muscle have not been fully elucidated. A potential key player in this process is the protein myoglobin (Mb). Indeed, there is a catalogue of empirical evidence supporting direct interaction of globins with fatty acid metabolites; however, the binding pocket and regulation of the interaction remains to be established. In this study, we employed a computational strategy to elucidate the structural determinants of fatty acids (palmitic & oleic acid) binding to Mb. Sequence analysis and docking simulations with a horse (Equus caballus) structural Mb reference reveals a fatty acid-binding site in the hydrophobic cleft near the heme region in Mb. Both palmitic acid and oleic acid attain a "U" shaped structure similar to their conformation in pockets of other fatty acid-binding proteins. Specifically, we found that the carboxyl head group of palmitic acid coordinates with the amino group of Lys45, whereas the carboxyl group of oleic acid coordinates with both the amino groups of Lys45 and Lys63. The alkyl tails of both fatty acids are supported by surrounding hydrophobic residues Leu29, Leu32, Phe33, Phe43, Phe46, Val67, Val68 and Ile107. In the saturated palmitic acid, the hydrophobic tail moves freely and occasionally penetrates deeper inside the hydrophobic cleft, making additional contacts with Val28, Leu69, Leu72 and Ile111. Our simulations reveal a dynamic and stable binding pocket in which the oxygen molecule and heme group in Mb are required for additional hydrophobic interactions. Taken together, these findings support a mechanism in which Mb acts as a muscle transporter for fatty acid when it is in the oxygenated state and releases fatty acid when Mb converts to deoxygenated state.

The BRCA1 tumor suppressor binds to inositol 1,4,5-trisphosphate receptors to stimulate apoptotic calcium release.

The inositol 1,4,5-trisphosphate receptor (IP3R) is a ubiquitously expressed endoplasmic reticulum (ER)-resident calcium channel. Calcium release mediated by IP3Rs influences many signaling pathways, including those regulating apoptosis. IP3R activity is regulated by protein-protein interactions, including binding to proto-oncogenes and tumor suppressors to regulate cell death. Here we show that the IP3R binds to the tumor suppressor BRCA1. BRCA1 binding directly sensitizes the IP3R to its ligand, IP3. BRCA1 is recruited to the ER during apoptosis in an IP3R-dependent manner, and, in addition, a pool of BRCA1 protein is constitutively associated with the ER under non-apoptotic conditions. This is likely mediated by a novel lipid binding activity of the first BRCA1 C terminus domain of BRCA1. These findings provide a mechanistic explanation by which BRCA1 can act as a proapoptotic protein.

Major diversification of voltage-gated K+ channels occurred in ancestral parahoxozoans.

We examined the origins and functional evolution of the Shaker and KCNQ families of voltage-gated K(+) channels to better understand how neuronal excitability evolved. In bilaterians, the Shaker family consists of four functionally distinct gene families (Shaker, Shab, Shal, and Shaw) that share a subunit structure consisting of a voltage-gated K(+) channel motif coupled to a cytoplasmic domain that mediates subfamily-exclusive assembly (T1). We traced the origin of this unique Shaker subunit structure to a common ancestor of ctenophores and parahoxozoans (cnidarians, bilaterians, and placozoans). Thus, the Shaker family is metazoan specific but is likely to have evolved in a basal metazoan. Phylogenetic analysis suggested that the Shaker subfamily could predate the divergence of ctenophores and parahoxozoans, but that the Shab, Shal, and Shaw subfamilies are parahoxozoan specific. In support of this, putative ctenophore Shaker subfamily channel subunits coassembled with cnidarian and mouse Shaker subunits, but not with cnidarian Shab, Shal, or Shaw subunits. The KCNQ family, which has a distinct subunit structure, also appears solely within the parahoxozoan lineage. Functional analysis indicated that the characteristic properties of Shaker, Shab, Shal, Shaw, and KCNQ currents evolved before the divergence of cnidarians and bilaterians. These results show that a major diversification of voltage-gated K(+) channels occurred in ancestral parahoxozoans and imply that many fundamental mechanisms for the regulation of action potential propagation evolved at this time. Our results further suggest that there are likely to be substantial differences in the regulation of neuronal excitability between ctenophores and parahoxozoans.

TRPV channel-mediated calcium transients in nociceptor neurons are dispensable for avoidance behaviour.

Animals need to sense and react to potentially dangerous environments. TRP ion channels participate in nociception, presumably via Ca(2+) influx, in most animal species. However, the relationship between ion permeation and animals' nocifensive behaviour is unknown. Here we use an invertebrate animal model with relevance for mammalian pain. We analyse the putative selectivity filter of OSM-9, a TRPV channel, in osmotic avoidance behaviour of Caenorhabditis elegans. Using mutagenized OSM-9 expressed in the head nociceptor neuron, ASH, we study nocifensive behaviour and Ca(2+) influx. Within the selectivity filter, M(601)-F(609), Y604G strongly reduces avoidance behaviour and eliminates Ca(2+) transients. Y604F also abolishes Ca(2+) transients in ASH, while sustaining avoidance behaviour, yet it disrupts behavioral plasticity. Homology modelling of the OSM-9 pore suggests that Y(604) may assume a scaffolding role. Thus, aromatic residues in the OSM-9 selectivity filter are critical for pain behaviour and ion permeation. These findings have relevance for understanding evolutionary roots of mammalian nociception.

Functional evolution of Erg potassium channel gating reveals an ancient origin for IKr.

Mammalian Ether-a-go-go related gene (Erg) family voltage-gated K(+) channels possess an unusual gating phenotype that specializes them for a role in delayed repolarization. Mammalian Erg currents rectify during depolarization due to rapid, voltage-dependent inactivation, but rebound during repolarization due to a combination of rapid recovery from inactivation and slow deactivation. This is exemplified by the mammalian Erg1 channel, which is responsible for IKr, a current that repolarizes cardiac action potential plateaus. The Drosophila Erg channel does not inactivate and closes rapidly upon repolarization. The dramatically different properties observed in mammalian and Drosophila Erg homologs bring into question the evolutionary origins of distinct Erg K(+) channel functions. Erg channels are highly conserved in eumetazoans and first evolved in a common ancestor of the placozoans, cnidarians, and bilaterians. To address the ancestral function of Erg channels, we identified and characterized Erg channel paralogs in the sea anemone Nematostella vectensis. N. vectensis Erg1 (NvErg1) is highly conserved with respect to bilaterian homologs and shares the IKr-like gating phenotype with mammalian Erg channels. Thus, the IKr phenotype predates the divergence of cnidarians and bilaterians. NvErg4 and Caenorhabditis elegans Erg (unc-103) share the divergent Drosophila Erg gating phenotype. Phylogenetic and sequence analysis surprisingly indicates that this alternate gating phenotype arose independently in protosomes and cnidarians. Conversion from an ancestral IKr-like gating phenotype to a Drosophila Erg-like phenotype correlates with loss of the cytoplasmic Ether-a-go-go domain. This domain is required for slow deactivation in mammalian Erg1 channels, and thus its loss may partially explain the change in gating phenotype.

Reevaluation of the evolutionary events within recA/RAD51 phylogeny.

BACKGROUND: The recA/RAD51 gene family encodes a diverse set of recombinase proteins that affect homologous recombination, DNA-repair, and genome stability. The recA gene family is expressed across all three domains of life - Eubacteria, Archaea, and Eukaryotes - and even in some viruses. To date, efforts to resolve the deep evolutionary origins of this ancient protein family have been hindered by the high sequence divergence between paralogous groups (i.e. ~30% average pairwise identity). RESULTS: Through large taxon sampling and the use of a phylogenetic algorithm designed for inferring evolutionary events in highly divergent paralogs, we obtained a robust, parsimonious and more refined phylogenetic history of the recA/RAD51 superfamily. CONCLUSIONS: In summary, our model for the evolution of recA/RAD51 family provides a better understanding of the ancient origin of recA proteins and the multiple events that lead to the diversification of recA homologs in eukaryotes, including the discovery of additional RAD51 sub-families.

PHYRN: a robust method for phylogenetic analysis of highly divergent sequences.

Both multiple sequence alignment and phylogenetic analysis are problematic in the "twilight zone" of sequence similarity (≤ 25% amino acid identity). Herein we explore the accuracy of phylogenetic inference at extreme sequence divergence using a variety of simulated data sets. We evaluate four leading multiple sequence alignment (MSA) methods (MAFFT, T-COFFEE, CLUSTAL, and MUSCLE) and six commonly used programs of tree estimation (Distance-based: Neighbor-Joining; Character-based: PhyML, RAxML, GARLI, Maximum Parsimony, and Bayesian) against a novel MSA-independent method (PHYRN) described here. Strikingly, at "midnight zone" genetic distances (~7% pairwise identity and 4.0 gaps per position), PHYRN returns high-resolution phylogenies that outperform traditional approaches. We reason this is due to PHRYN's capability to amplify informative positions, even at the most extreme levels of sequence divergence. We also assess the applicability of the PHYRN algorithm for inferring deep evolutionary relationships in the divergent DANGER protein superfamily, for which PHYRN infers a more robust tree compared to MSA-based approaches. Taken together, these results demonstrate that PHYRN represents a powerful mechanism for mapping uncharted frontiers in highly divergent protein sequence data sets.

Conserved GXXXG- and S/T-like motifs in the transmembrane domains of NS4B protein are required for hepatitis C virus replication.

Hepatitis C virus (HCV) nonstructural protein 4B (NS4B) is an integral membrane protein, which plays an important role in the organization and function of the HCV replication complex (RC). Although much is understood about its amphipathic N-terminal and C-terminal domains, we know very little about the role of the transmembrane domains (TMDs) in NS4B function. We hypothesized that in addition to anchoring NS4B into host membranes, the TMDs are engaged in intra- and intermolecular interactions required for NS4B structure/function. To test this hypothesis, we have engineered a chimeric JFH1 genome containing the Con1 NS4B TMD region. The resulting virus titers were greatly reduced from those of JFH1, and further analysis indicated a defect in genome replication. We have mapped this incompatibility to NS4B TMD1 and TMD2 sequences, and we have defined putative TMD dimerization motifs (GXXXG in TMD2 and TMD3; the S/T cluster in TMD1) as key structural/functional determinants. Mutations in each of the putative motifs led to significant decreases in JFH1 replication. Like most of the NS4B chimeras, mutant proteins had no negative impact on NS4B membrane association. However, some mutations led to disruption of NS4B foci, implying that the TMDs play a role in HCV RC formation. Further examination indicated that the loss of NS4B foci correlates with the destabilization of NS4B protein. Finally, we have identified an adaptive mutation in the NS4B TMD2 sequence that has compensatory effects on JFH1 chimera replication. Taken together, these data underscore the functional importance of NS4B TMDs in the HCV life cycle.