Space trajectories optimization using variable-chromosome-length genetic algorithms

The problem of optimal design of a multi-gravity-assist space trajectories, with free number of deep space maneuvers (MGADSM) poses multi-modal cost functions. In the general form of the problem, the number of design variables is solution dependent. To handle global optimization problems where the number of design variables varies from one solution to another, two novel genetic-based techniques are introduced: hidden genes genetic algorithm (HGGA) and dynamic-size multiple population genetic algorithm (DSMPGA). In HGGA, a fixed length for the design variables is assigned for all solutions. Independent variables of each solution are divided into effective and ineffective (hidden) genes. Hidden genes are excluded in cost function evaluations. Full-length solutions undergo standard genetic operations. In DSMPGA, sub-populations of fixed size design spaces are randomly initialized. Standard genetic operations are carried out for a stage of generations. A new population is then created by reproduction from all members based on their relative fitness. The resulting sub-populations have different sizes from their initial sizes. The process repeats, leading to increasing the size of sub-populations of more fit solutions. Both techniques are applied to several MGADSM problems. They have the capability to determine the number of swing-bys, the planets to swing by, launch and arrival dates, and the number of deep space maneuvers as well as their locations, magnitudes, and directions in an optimal sense. The results show that solutions obtained using the developed tools match known solutions for complex case studies. The HGGA is also used to obtain the asteroids sequence and the mission structure in the global trajectory optimization competition (GTOC) problem. As an application of GA optimization to Earth orbits, the problem of visiting a set of ground sites within a constrained time frame is solved. The J2 perturbation and zonal coverage are considered to design repeated Sun-synchronous orbits. Finally, a new set of orbits, the repeated shadow track orbits (RSTO), is introduced. The orbit parameters are optimized such that the shadow of a spacecraft on the Earth visits the same locations periodically every desired number of days.

Variable impact on mortality of AIDS-defining events diagnosed during combination antiretroviral therapy: not all AIDS-defining conditions are created equal

BACKGROUND: The extent to which mortality differs following individual acquired immunodeficiency syndrome (AIDS)-defining events (ADEs) has not been assessed among patients initiating combination antiretroviral therapy. METHODS: We analyzed data from 31,620 patients with no prior ADEs who started combination antiretroviral therapy. Cox proportional hazards models were used to estimate mortality hazard ratios for each ADE that occurred in >50 patients, after stratification by cohort and adjustment for sex, HIV transmission group, number of antiretroviral drugs initiated, regimen, age, date of starting combination antiretroviral therapy, and CD4+ cell count and HIV RNA load at initiation of combination antiretroviral therapy. ADEs that occurred in <50 patients were grouped together to form a "rare ADEs" category. RESULTS: During a median follow-up period of 43 months (interquartile range, 19-70 months), 2880 ADEs were diagnosed in 2262 patients; 1146 patients died. The most common ADEs were esophageal candidiasis (in 360 patients), Pneumocystis jiroveci pneumonia (320 patients), and Kaposi sarcoma (308 patients). The greatest mortality hazard ratio was associated with non-Hodgkin's lymphoma (hazard ratio, 17.59; 95% confidence interval, 13.84-22.35) and progressive multifocal leukoencephalopathy (hazard ratio, 10.0; 95% confidence interval, 6.70-14.92). Three groups of ADEs were identified on the basis of the ranked hazard ratios with bootstrapped confidence intervals: severe (non-Hodgkin's lymphoma and progressive multifocal leukoencephalopathy [hazard ratio, 7.26; 95% confidence interval, 5.55-9.48]), moderate (cryptococcosis, cerebral toxoplasmosis, AIDS dementia complex, disseminated Mycobacterium avium complex, and rare ADEs [hazard ratio, 2.35; 95% confidence interval, 1.76-3.13]), and mild (all other ADEs [hazard ratio, 1.47; 95% confidence interval, 1.08-2.00]). CONCLUSIONS: In the combination antiretroviral therapy era, mortality rates subsequent to an ADE depend on the specific diagnosis. The proposed classification of ADEs may be useful in clinical end point trials, prognostic studies, and patient management.

Expression Analysis of the Theileria parva Subtelomere-Encoded Variable Secreted Protein Gene Family

Background The intracellular protozoan parasite Theileria parva transforms bovine lymphocytes inducing uncontrolled proliferation. Proteins released from the parasite are assumed to contribute to phenotypic changes of the host cell and parasite persistence. With 85 members, genes encoding subtelomeric variable secreted proteins (SVSPs) form the largest gene family in T. parva. The majority of SVSPs contain predicted signal peptides, suggesting secretion into the host cell cytoplasm. Methodology/Principal Findings We analysed SVSP expression in T. parva-transformed cell lines established in vitro by infection of T or B lymphocytes with cloned T. parva parasites. Microarray and quantitative real-time PCR analysis revealed mRNA expression for a wide range of SVSP genes. The pattern of mRNA expression was largely defined by the parasite genotype and not by host background or cell type, and found to be relatively stable in vitro over a period of two months. Interestingly, immunofluorescence analysis carried out on cell lines established from a cloned parasite showed that expression of a single SVSP encoded by TP03_0882 is limited to only a small percentage of parasites. Epitope-tagged TP03_0882 expressed in mammalian cells was found to translocate into the nucleus, a process that could be attributed to two different nuclear localisation signals. Conclusions Our analysis reveals a complex pattern of Theileria SVSP mRNA expression, which depends on the parasite genotype. Whereas in cell lines established from a cloned parasite transcripts can be found corresponding to a wide range of SVSP genes, only a minority of parasites appear to express a particular SVSP protein. The fact that a number of SVSPs contain functional nuclear localisation signals suggests that proteins released from the parasite could contribute to phenotypic changes of the host cell. This initial characterisation will facilitate future studies on the regulation of SVSP gene expression and the potential biological role of these enigmatic proteins.

Diversification in the South American Pampas: the genetic and morphological variation of the widespread Petunia axillaris complex (Solanaceae)

Understanding the spatiotemporal distribution of genetic variation and the ways in which this distribution is connected to the ecological context of natural populations is fundamental for understanding the nature and mode of intraspecific and, ultimately, interspecific differentiation. The Petunia axillaris complex is endemic to the grasslands of southern South America and includes three subspecies: P.a.axillaris, P.a.parodii and P.a.subandina. These subspecies are traditionally delimited based on both geography and floral morphology, although the latter is highly variable. Here, we determined the patterns of genetic (nuclear and cpDNA), morphological and ecological (bioclimatic) variation of a large number of P.axillaris populations and found that they are mostly coincident with subspecies delimitation. The nuclear data suggest that the subspecies are likely independent evolutionary units, and their morphological differences may be associated with local adaptations to diverse climatic and/or edaphic conditions and population isolation. The demographic dynamics over time estimated by skyline plot analyses showed different patterns for each subspecies in the last 100000years, which is compatible with a divergence time between 35000 and 107000years ago between P.a.axillaris and P.a.parodii, as estimated with the IMa program. Coalescent simulation tests using Approximate Bayesian Computation do not support previous suggestions of extensive gene flow between P.a.axillaris and P.a.parodii in their contact zone.

Crystal chemistry and hydrogen bonding of rustumite Ca10(Si2O7)2(SiO4)(OH)2Cl2 with variable OH, Cl, F

Three samples of the skarn mineral rustumite Ca10(Si2O7)2(SiO4)(OH)2Cl2, space group C2/c, a ≈7.6, b ≈ 18.5, c ≈ 15.5 Å, β ≈ 104°, with variable OH, Cl, F content were investigated by electron microprobe, single-crystal X-ray structure refinements, and Raman spectroscopy. “Rust1LCl” is a low chlorine rustumite Ca10(Si2O7)2(SiO4)(OH1.88F0.12)(Cl1.28,OH0.72) from skarns associated with the Rize batholith near Ikizedere, Turkey. “Rust2F” is a F-bearing rustumite Ca10(Si2O7)2(SiO4)(OH1.13F0.87) (Cl1 96OH0.04) from xenoliths in ignimbrites of the Upper Chegem Caldera, Northern Caucasus, Russia. “Rust3LClF” represents a low-Cl, F-bearing rustumite Ca10(Si2O7)2(SiO4)0.87(H4O4)0.13(OH1.01F0.99) (Cl1.00 OH1.00) from altered merwinite skarns of the Birkhin massif, Baikal Lake area, Eastern Siberia, Russia. Rustumite from Birkhin massif is characterized by a significant hydrogarnet-like or fluorine substitution at the apices of the orthosilicate group, leading to specific atomic displacements. The crystal structures including hydrogen positions have been refined from single-crystal X-ray data to R1 = 0.0205 (Rust1_LCl), R1 = 0.0295 (Rust2_F), and R1 = 0.0243 (Rust3_LCl_F), respectively. Depletion in Cl and replacement by OH is associated with smaller unit-cell dimensions. The substitution of OH by F leads to shorter hydrogen bonds O-H⋯F instead of O-H⋯OH. Raman spectra for all samples have been measured and confirm slight strengthening of the hydrogen bonds with uptake of F.This study discusses the complex crystal chemistry of the skarn mineral rustumite and may provide a wider understanding of the chemical reactions related to contact metamorphism of limestones.

Transcriptional regulation of the Borrelia burgdorferi antigenically variable VlsE surface protein.

The Lyme disease agent Borrelia burgdorferi can persistently infect humans and other animals despite host active immune responses. This is facilitated, in part, by the vls locus, a complex system consisting of the vlsE expression site and an adjacent set of 11 to 15 silent vls cassettes. Segments of nonexpressed cassettes recombine with the vlsE region during infection of mammalian hosts, resulting in combinatorial antigenic variation of the VlsE outer surface protein. We now demonstrate that synthesis of VlsE is regulated during the natural mammal-tick infectious cycle, being activated in mammals but repressed during tick colonization. Examination of cultured B. burgdorferi cells indicated that the spirochete controls vlsE transcription levels in response to environmental cues. Analysis of PvlsE::gfp fusions in B. burgdorferi indicated that VlsE production is controlled at the level of transcriptional initiation, and regions of 5' DNA involved in the regulation were identified. Electrophoretic mobility shift assays detected qualitative and quantitative changes in patterns of protein-DNA complexes formed between the vlsE promoter and cytoplasmic proteins, suggesting the involvement of DNA-binding proteins in the regulation of vlsE, with at least one protein acting as a transcriptional activator.

Variable effects of the conserved RNA hairpin element upon 3' end processing of histone pre-mRNA in vitro.

We have studied the requirements for efficient histone-specific RNA 3' processing in nuclear extract from mammalian tissue culture cells. Processing is strongly impaired by mutations in the pre-mRNA spacer element that reduce the base-pairing potential with U7 RNA. Moreover, by exchanging the hairpin and spacer elements of two differently processed H4 genes, we find that this difference is exclusively due to the spacer element. Finally, processing is inhibited by the addition of competitor RNAs, if these contain a wild-type spacer sequence, but not if their spacer element is mutated. Conversely, the importance of the hairpin for histone RNA 3' processing is highly variable: A hairpin mutant of the H4-12 gene is processed with almost wild-type efficiency in extract from K21 mouse mastocytoma cells but is strongly affected in HeLa cell extract, whereas an identical hairpin mutant of the H4-1 gene is affected in both extracts. The hairpin defect of H4-12-specific RNA in HeLa cells can be overcome by a compensatory mutation that increases the base complementarity to U7 snRNA. Very similar results were also obtained in RNA competition experiments: processing of H4-12-specific RNA can be competed by RNA carrying a wild-type hairpin element in extract from HeLa, but not K21 cells, whereas processing of H4-1-specific RNA can be competed in both extracts. With two additional histone genes we obtained results that were in one case intermediate and in the other similar to those obtained with H4-1. These results suggest that hairpin binding factor(s) can cooperatively support the ability of U7 snRNPs to form an active processing complex, but is(are) not directly involved in the processing mechanism.

Simple-Offer vs. Complex-Offer Auctions in Deregulated Electricity Markets

In my recent experimental research of wholesale electricity auctions, I discovered that the complex structure of the offers leaves a lot of room for strategic behavior, which consequently leads to anti- competitive and inefficient outcomes in the market. A specific feature of these complex-offer auctions is that the sellers submit not only the quantities and the minimum prices at which they are willing to sell, but also the start-up fees that are designed to reimburse the fixed start-up costs of the generation plants. In this paper, using the experimental method I compare the performance of two complex-offer auctions (COAs) against the performance of a simple-offer auction (SOA), in which the sellers have to recover all their generation costs --- fixed and variable ---through a uniform market-clearing price. I find that the SOA significantly reduces consumer prices and lowers price volatility. It mitigates anti-competitive effects that are present in the COAs and achieves allocative efficiency more quickly.

An empirical evaluation of the Random Forests classifier models for variable selection in a large-scale lung cancer case-control study

Random Forests™ is reported to be one of the most accurate classification algorithms in complex data analysis. It shows excellent performance even when most predictors are noisy and the number of variables is much larger than the number of observations. In this thesis Random Forests was applied to a large-scale lung cancer case-control study. A novel way of automatically selecting prognostic factors was proposed. Also, synthetic positive control was used to validate Random Forests method. Throughout this study we showed that Random Forests can deal with large number of weak input variables without overfitting. It can account for non-additive interactions between these input variables. Random Forests can also be used for variable selection without being adversely affected by collinearities. ^ Random Forests can deal with the large-scale data sets without rigorous data preprocessing. It has robust variable importance ranking measure. Proposed is a novel variable selection method in context of Random Forests that uses the data noise level as the cut-off value to determine the subset of the important predictors. This new approach enhanced the ability of the Random Forests algorithm to automatically identify important predictors for complex data. The cut-off value can also be adjusted based on the results of the synthetic positive control experiments. ^ When the data set had high variables to observations ratio, Random Forests complemented the established logistic regression. This study suggested that Random Forests is recommended for such high dimensionality data. One can use Random Forests to select the important variables and then use logistic regression or Random Forests itself to estimate the effect size of the predictors and to classify new observations. ^ We also found that the mean decrease of accuracy is a more reliable variable ranking measurement than mean decrease of Gini. ^

Geochemistry of the lower sheeted dike complex of ODP Holes 137-504B and 140-504B

Sixty-three samples representing 379 m of sheeted dikes from Deep Sea Drilling Project/Ocean Drilling Program Site 504B have been analyzed for major and selected trace elements by X-ray fluorescence. The samples range from microcrystalline aphyric basalts to moderately phyric (2%-10% phenocrysts) diabase that are typically multiply saturated with plagioclase, olivine, and clinopyroxene, in order of relative abundance. All analyzed samples are classified as Group D compositions with moderate to slightly elevated compatible elements (MgÆ-value = 0.65% ± 0.03%; Al2O3 = 15.5% ± 0.8%; CaO = 13.0% ± 0.3%; Ni = 114 ± 29 ppm), and unusually depleted levels of moderate to highly incompatible elements (Nb < 1 ppm; Zr = 44 ± 7 ppm; Rb < 0.5 ppm; Ba ~ 1 ppm; P2O5 = 0.07% ± 0.02%). These compositions are consistent with a multistage melting of a normal ocean ridge basaltic mantle source followed by extensive fractionation of olivine, plagioclase, and clinopyroxene. Leg 140 aphyric to sparsely phyric (0%-2% phenocrysts) basalts and diabases are compositionally indistinguishable from similarly phyric samples at higher levels in the hole. An examination of the entire crustal section, from the overlying volcanics through the sheeted dikes observed in Leg 140, reveals no significant trends indicating the enrichment or depletion of Costa Rica Rift Zone source magmas over time. Similarly, significant trends toward increased or decreased differentiation cannot be identified, although compositional patterns reflecting variable amounts of phenocryst addition are apparent at various depths. Below ? 1700 mbsf to the bottom of the Leg 140 section, there is a broadly systematic pattern of Zn depletion with depth, the result of high-temperature hydrothermal leaching. This zone of depletion is thought to be a significant source of Zn for the hydrothermal fluids depositing metal sulfides at ridge-crest hydrothermal vents and the sulfide-mineralization zone, located in the transition between pillow lavas and sheeted dikes. Localized zones of intense alteration (60%-95% recrystallization) are present on a centimeter to meter scale in many lithologic units. Within these zones, normally immobile elements Ti, Zr, Y, and rare-earth elements are strongly depleted compared with "fresher" samples centimeters away. The extent of compositional variability of these elements tends to obscure primary igneous trends if the highly altered samples are not identified or removed. At levels up to 40% (or possibly 60%) recrystallization, Ti, Zr, and Y retain their primary signatures. Although the mechanisms are unclear, it is possible that these intense alteration zones are a source of Y and rare-earth elements for the typically rare-earth-element-enriched hydrothermal vent fluids of mid-ocean ridges.

Estudio numérico de flujo turbulento cargado con partículas sólidas a través de canales y tuberías de sección variable

La presente Tesis Doctoral tiene como objetivo el estudio de flujo turbulento cargado con partículas sólidas a través de canales y tuberías de sección constante usando un enfoque Euleriano-Lagrangiano. El campo de flujo de la fase de transporte (aire) se resuelve usando simulación de grandes escalas (LES), implementada en un programa de volúmenes finitos mientras que las ecuaciones gobernantes de la fase dispersa son resueltas por medio de un algoritmo de seguimiento Lagrangiano de partículas que ha sido desarrollado y acoplado al programa que resuelve el flujo. Se estudia de manera sistemática y progresiva la interacción fluido→partícula (one-way coupling), a través de diferentes configuraciones geométricas en coordenadas cartesianas (canales de sección constante y variable) y en coordenadas cilíndricas (tuberías de sección constante y sección variable) abarcando diferentes números de Reynolds y diferentes tamaños de partículas; todos los resultados obtenidos han sido comparados con datos publicados previamente. El estudio de flujo multifásico a través de, tuberías de sección variable, ha sido abordada en otras investigaciones mayoritariamente de forma experimental o mediante simulación usando modelos de turbulencia menos complejos y no mediante LES. El patrón de flujo que se verifica en una tubería con expansión es muy complejo y dicha configuración geométrica se halla en múltiples aplicaciones industriales que involucran el transporte de partículas sólidas, por ello es de gran interés su estudio. Como hecho innovador, en esta tesis no solo se resuelven las estadísticas de velocidad del fluido y las partículas en tuberías con diferentes tamaños de expansión y diferentes regímenes de flujo sino que se caracteriza, usando diversas formulaciones del número de Stokes y el parámetro de arrastre, el ingreso y acumulación de partículas dentro de la zona de recirculación, obteniéndose resultados coincidentes con datos experimentales. ABSTRACT The objective of this Thesis research is to study the turbulent flow laden with solid particles through channels and pipes with using Eulerian-Lagrangian approach. The flow field of the transport phase (air ) is solved using large eddy simulation ( LES ) implemented in a program of finite volume while the governing equations of the dispersed phase are resolved by means of a particle Lagrangian tracking algorithm which was developed and coupled to principal program flow solver . We studied systematically and progressively the fluid interaction → particle ( one- way coupling ) , through different geometric configurations in Cartesian coordinates ( channel with constant and variable section) and in cylindrical coordinates ( pipes with constant section and variable section ) covering different Reynolds numbers and different particle sizes, all results have been compared with previously published data . The study of multiphase flow through, pipes with variable section has been addressed in other investigations predominantly experimentally or by simulation using less complex models and no turbulence by LES. The flow pattern is verified in a pipe expansion is very complex and this geometry is found in many industrial applications involving the transport of solid particles, so it is of great interest to study. As an innovator fact , in this Thesis not only finds fluid velocity statistics and particles with different sizes of pipe expansion and different flow regimes but characterized, using various formulations of the Stokes number and the drag parameter are resolved, the entry and accumulation of particles within the recirculation zone , matching results obtained with experimental data.

Molecular characteristics of fibroblast growth factor–fibroblast growth factor receptor–heparin-like glycosaminoglycan complex

Fibroblast growth factor (FGF) family plays key roles in development, wound healing, and angiogenesis. Understanding of the molecular nature of interactions of FGFs with their receptors (FGFRs) has been seriously limited by the absence of structural information on FGFR or FGF–FGFR complex. In this study, based on an exhaustive analysis of the primary sequences of the FGF family, we determined that the residues that constitute the primary receptor-binding site of FGF-2 are conserved throughout the FGF family, whereas those of the secondary receptor binding site of FGF-2 are not. We propose that the FGF–FGFR interaction mediated by the ‘conserved’ primary site interactions is likely to be similar if not identical for the entire FGF family, whereas the ‘variable’ secondary sites, on both FGF as well as FGFR mediates specificity of a given FGF to a given FGFR isoform. Furthermore, as the pro-inflammatory cytokine interleukin 1 (IL-1) and FGF-2 share the same structural scaffold, we find that the spatial orientation of the primary receptor-binding site of FGF-2 coincides structurally with the IL-1β receptor-binding site when the two molecules are superimposed. The structural similarities between the IL-1 and the FGF system provided a framework to elucidate molecular principles of FGF–FGFR interactions. In the FGF–FGFR model proposed here, the two domains of a single FGFR wrap around a single FGF-2 molecule such that one domain of FGFR binds to the primary receptor-binding site of the FGF molecule, while the second domain of the same FGFR binds to the secondary receptor-binding site of the same FGF molecule. Finally, the proposed model is able to accommodate not only heparin-like glycosaminoglycan (HLGAG) interactions with FGF and FGFR but also FGF dimerization or oligomerization mediated by HLGAG.

Regulation of T cell receptor (TCR) β gene expression by CD3 complex signaling in immature thymocytes: Implications for TCRβ allelic exclusion

During αβ thymocyte development, clonotype-independent CD3 complexes are expressed at the cell surface before the pre-T cell receptor (TCR). Signaling through clonotype-independent CD3 complexes is required for expression of rearranged TCRβ genes. On expression of a TCRβ polypeptide chain, the pre-TCR is assembled, and TCRβ locus allelic exclusion is established. We investigated the putative contribution of clonotype-independent CD3 complex signaling to TCRβ locus allelic exclusion in mice single-deficient or double-deficient for CD3ζ/η and/or p56lck. These mice display defects in the expression of endogenous TCRβ genes in immature thymocytes, proportional to the severity of CD3 complex malfunction. Exclusion of endogenous TCRβ VDJ (variable, diversity, joining) rearrangements by a functional TCRβ transgene was severely compromised in the single-deficient and double-deficient mutant mice. In contrast to wild-type mice, most of the CD25+ double-negative (DN) thymocytes of the mutant mice failed to express the TCRβ transgene, suggesting defective expression of the TCRβ transgene similar to endogenous TCRβ genes. In the mutant mice, a proportion of CD25+ DN thymocytes that failed to express the transgene expressed endogenous TCRβ polypeptide chains. Many double-positive cells of the mutant mice coexpressed endogenous and transgenic TCRβ chains or more than one endogenous TCRβ chain. The data suggest that signaling through clonotype-independent CD3 complexes may contribute to allelic exclusion of the TCRβ locus by inducing the expression of rearranged TCRβ genes in CD25+ DN thymocytes.

A family of phase-variable restriction enzymes with differing specificities generated by high-frequency gene rearrangements

The hsd genes of Mycoplasma pulmonis encode restriction and modification enzymes exhibiting a high degree of sequence similarity to the type I enzymes of enteric bacteria. The S subunits of type I systems dictate the DNA sequence specificity of the holoenzyme and are required for both the restriction and the modification reactions. The M. pulmonis chromosome has two hsd loci, both of which contain two hsdS genes each and are complex, site-specific DNA inversion systems. Embedded within the coding region of each hsdS gene are a minimum of three sites at which DNA inversions occur to generate extensive amino acid sequence variations in the predicted S subunits. We show that the polymorphic hsdS genes produced by gene rearrangement encode a family of functional S subunits with differing DNA sequence specificities. In addition to creating polymorphisms in hsdS sequences, DNA inversions regulate the phase-variable production of restriction activity because the other genes required for restriction activity (hsdR and hsdM) are expressed only from loci that are oriented appropriately in the chromosome relative to the hsd promoter. These data cast doubt on the prevailing paradigms that restriction systems are either selfish or function to confer protection from invasion by foreign DNA.

Conserved structural features of the synaptic fusion complex: SNARE proteins reclassified as Q- and R-SNAREs

SNARE [soluble NSF (N-ethylmaleimide-sensitive fusion protein) attachment protein receptor] proteins are essential for membrane fusion and are conserved from yeast to humans. Sequence alignments of the most conserved regions were mapped onto the recently solved crystal structure of the heterotrimeric synaptic fusion complex. The association of the four α-helices in the synaptic fusion complex structure produces highly conserved layers of interacting amino acid side chains in the center of the four-helix bundle. Mutations in these layers reduce complex stability and cause defects in membrane traffic even in distantly related SNAREs. When syntaxin-4 is modeled into the synaptic fusion complex as a replacement of syntaxin-1A, no major steric clashes arise and the most variable amino acids localize to the outer surface of the complex. We conclude that the main structural features of the neuronal complex are highly conserved during evolution. On the basis of these features we have reclassified SNARE proteins into Q-SNAREs and R-SNAREs, and we propose that fusion-competent SNARE complexes generally consist of four-helix bundles composed of three Q-SNAREs and one R-SNARE.