Non-homogeneous models of sequence evolution in the Bio++ suite of libraries and programs

Background  

Accurately modeling the sequence substitution process is required for the correct estimation of evolutionary parameters, be they phylogenetic relationships, substitution rates or ancestral states; it is also crucial to simulate realistic data sets. Such simulation procedures are needed to estimate the null-distribution of complex statistics, an approach referred to as parametric bootstrapping, and are also used to test the quality of phylogenetic reconstruction programs. It has often been observed that homologous sequences can vary widely in their nucleotide or amino-acid compositions, revealing that sequence evolution has changed importantly among lineages, and may therefore be most appropriately approached through non-homogeneous models. Several programs implementing such models have been developed, but they are limited in their possibilities: only a few particular models are available for likelihood optimization, and data sets cannot be easily generated using the resulting estimated parameters.     

Dodecyl Maltoside Protects Membrane Proteins In Vacuo

Molecular dynamics simulations have been used to characterize the effects of transfer from aqueous solution to a vacuum to inform our understanding of mass spectrometry of membrane-protein-detergent complexes. We compared two membrane protein architectures (an α-helical bundle versus a β-barrel) and two different detergent types (phosphocholines versus an alkyl sugar) with respect to protein stability and detergent packing. The β-barrel membrane protein remained stable as a protein-detergent complex in vacuum. Zwitterionic detergents formed conformationally destabilizing interactions with an α-helical membrane protein after detergent micelle inversion driven by dehydration in vacuum. In contrast, a nonionic alkyl sugar detergent resisted micelle inversion, maintaining the solution-phase conformation of the protein. This helps to explain the relative stability of membrane proteins in the presence of alkyl sugar detergents such as dodecyl maltoside.     

New perspectives of locomotor rehabilitation after stroke

The task-oriented approach incorporating treadmill walking for retraining gait early after stroke has contributed to promote locomotor recovery. To augment practice, training strategies such as mental practice and training in virtual environments are proposed. While the former offers more practice with less physical exertion, the latter allows safe practice in a variety of challenging environments. Work is under way to assess whether these new strategies can further enhance locomotor recovery.     

Darker female pigeons transmit more specific antibodies to their eggs than do paler ones

Melanin‐based coloration is widespread among vertebrates, yet the adaptive significance of such pigments remains elusive, particularly with regard to the link between melanin and immune‐mediated maternal effects. The aim of this study was to investigate whether melanin‐based coloration could signal the ability of mothers to mount a humoral response and to transfer maternal antibodies (Ab) to their young. We injected differently coloured (pale and dark) female feral pigeons (Columba livia) with Chlamydiae (a natural antigen) and Keyhole Limpet Haemocyanin (KLH, an artificial antigen), and found no significant difference in humoral response between differently coloured females. However, darker females transferred more Ab against Chlamydiae into their eggs than paler ones, despite similar circulating levels of Ab. In addition to this, melanin‐based coloration showed a high heritability value. This suggests that a genetically based coloured trait might be linked to the ability of females to transfer specific Ab against Chlamydiae (but not against KLH) to their offspring, independent of their ability to produce Ab. This suggests that transmission of maternal Ab is antigen dependent, and that melanin‐based coloration might signal female ability to transmit specific Ab against natural pathogens. © 2013 The Linnean Society of London     

Crystallographic Definition of the Epitope Promiscuity of the Broadly Neutralizing Anti-Human Immunodeficiency Virus Type 1 Antibody 2F5: Vaccine Design Implications

The quest to create a human immunodeficiency virus type 1 (HIV-1) vaccine capable of eliciting broadly neutralizing antibodies against Env has been challenging. Among other problems, one difficulty in creating a potent immunogen resides in the substantial overall sequence variability of the HIV envelope protein. The membrane-proximal region (MPER) of gp41 is a particularly conserved tryptophan-rich region spanning residues 659 to 683, which is recognized by three broadly neutralizing monoclonal antibodies (bnMAbs), 2F5, Z13, and 4E10. In this study, we first describe the variability of residues in the gp41 MPER and report on the invariant nature of 15 out of 25 amino acids comprising this region. Subsequently, we evaluate the ability of the bnMAb 2F5 to recognize 31 varying sequences of the gp41 MPER at a molecular level. In 19 cases, resulting crystal structures show the various MPER peptides bound to the 2F5 Fab′. A variety of amino acid substitutions outside the ⁶⁶⁴DKW⁶⁶⁶ core epitope are tolerated. However, changes at the ⁶⁶⁴DKW⁶⁶⁶ motif itself are restricted to those residues that preserve the aspartate''s negative charge, the hydrophobic alkyl-π stacking arrangement between the β-turn lysine and tryptophan, and the positive charge of the former. We also characterize a possible molecular mechanism of 2F5 escape by sequence variability at position 667, which is often observed in HIV-1 clade C isolates. Based on our results, we propose a somewhat more flexible molecular model of epitope recognition by bnMAb 2F5, which could guide future attempts at designing small-molecule MPER-like vaccines capable of eliciting 2F5-like antibodies.Eliciting broadly neutralizing antibodies (bnAbs) against primary isolates of human immunodeficiency virus type I (HIV-1) has been identified as a major milestone to attain in the quest for a vaccine in the fight against AIDS (12, 28). These antibodies would need to interact with HIV-1 envelope glycoproteins gp41 and/or gp120 (Env), target conserved regions and functional conformations of gp41/gp120 trimeric complexes, and prevent new HIV-1 fusion events with target cells (21, 57, 70, 71). Although a humoral response generating neutralizing antibodies against HIV-1 can be detected in HIV-1-positive individuals, the titers are often very low, and virus control is seldom achieved by these neutralizing antibodies (22, 51, 52, 66, 67). The difficulty in eliciting a broad and potent neutralizing antibody response against HIV-1 is thought to reside in the high degree of genetic diversity of the virus, in the heterogeneity of Env on the surface of HIV-1, and in the masking of functional regions by conformational covering, by an extensive glycan shield, or by the ability of some conserved domains to partition to the viral membrane (24, 25, 29, 30, 38, 39, 56, 68, 69). So far, vaccine trials using as immunogens mimics of Env in different conformations have primarily elicited antibodies with only limited neutralization potency across different HIV-1 clades although recent work has demonstrated more encouraging results (4, 12, 61).The use of conserved regions on gp41 and gp120 Env as targets for vaccine design has been mostly characterized by the very few anti-HIV-1 broadly neutralizing monoclonal antibodies (bnMAbs) that recognize them: the CD4 binding-site on gp120 (bnMAb b12), a CD4-induced gp120 coreceptor binding site (bnMAbs 17b and X5), a mannose cluster on the outer face of gp120 (bnMAb 2G12), and the membrane proximal external region (MPER) of gp41 (bnMAbs 2F5, Z13 and 4E10) (13, 29, 44, 58, 73). The gp41 MPER region is a particularly conserved part of Env that spans residues 659 to 683 (HXB2 numbering) (37, 75). Substitution and deletion studies have linked this unusually tryptophan-rich region to the fusion process of HIV-1, possibly involving a series of conformational changes (5, 37, 41, 49, 54, 74). Additionally, the gp41 MPER has been implicated in gp41 oligomerization, membrane leakage ability facilitating pore formation, and binding to the galactosyl ceramide receptor on epithelial cells for initial mucosal infection mediated by transcytosis (2, 3, 40, 53, 63, 64, 72). This wide array of roles for the gp41 MPER will put considerable pressure on sequence conservation, and any change will certainly lead to a high cost in viral fitness.Monoclonal antibody 2F5 is a broadly neutralizing monoclonal anti-HIV-1 antibody isolated from a panel of sera from naturally infected asymptomatic individuals. It reacts with a core gp41 MPER epitope spanning residues 662 to 668 with the linear sequence ELDKWAS (6, 11, 42, 62, 75). 2F5 immunoglobulin G binding studies and screening of phage display libraries demonstrated that the DKW core is essential for 2F5 recognition and binding (15, 36, 50). Crystal structures of 2F5 with peptides representing its core gp41 epitope reveal a β-turn conformation involving the central DKW residues, flanked by an extended conformation and a canonical α-helical turn for residues located at the N terminus and C terminus of the core, respectively (9, 27, 45, 47). In addition to binding to its primary epitope, evidence is accumulating that 2F5 also undergoes secondary interactions: multiple reports have demonstrated affinity of 2F5 for membrane components, possibly through its partly hydrophobic flexible elongated complementarity-determining region (CDR) H3 loop, and it has also been suggested that 2F5 might interact in a secondary manner with other regions of gp41 (1, 10, 23, 32, 33, 55). Altogether, even though the characteristics of 2F5 interaction with its linear MPER consensus epitope have been described extensively, a number of questions persist about the exact mechanism of 2F5 neutralization at a molecular level.One such ambiguous area of the neutralization mechanism of 2F5 is investigated in this study. Indeed, compared to bnMAb 4E10, 2F5 is the more potent neutralizing antibody although its breadth across different HIV-1 isolates is more limited (6, 35). In an attempt to shed light on the exact molecular requirements for 2F5 recognition of its primary gp41 MPER epitope, we performed structural studies of 2F5 Fab′ with a variety of peptides. The remarkable breadth of possible 2F5 interactions reveals a somewhat surprising promiscuity of the 2F5 binding site. Furthermore, we link our structural observations with the natural variation observed within the gp41 MPER and discuss possible routes of 2F5 escape from a molecular standpoint. Finally, our discovery of 2F5''s ability to tolerate a rather broad spectrum of amino acids in its binding, a spectrum that even includes nonnatural amino acids, opens the door to new ways to design small-molecule immunogens potentially capable of eliciting 2F5-like neutralizing antibodies.     

Costs and benefits of multiple resistance to insecticides for <Emphasis Type="Italic">Culex quinquefasciatus</Emphasis> mosquitoes

Background  

The evolutionary dynamics of xenobiotic resistance depends on how resistance mutations influence the fitness of their bearers, both in the presence and absence of xenobiotic selection pressure. In cases of multiple resistance, these dynamics will also depend on how individual resistance mutations interact with one another, and on the xenobiotics applied against them. We compared Culex quinquefasciatus mosquitoes harbouring two resistance alleles ace-1 ^R and Kdr ^R (conferring resistance to carbamate and pyrethroid insecticides, respectively) to mosquitoes bearing only one of the alleles, or neither allele. Comparisons were made in environments where both, only one, or neither type of insecticide was present.     

Metabolic Modeling of Dynamic Brain 13C NMR Multiplet Data: Concepts and Simulations with a Two-Compartment Neuronal-Glial Model

Metabolic modeling of dynamic ¹³C labeling curves during infusion of ¹³C-labeled substrates allows quantitative measurements of metabolic rates in vivo. However metabolic modeling studies performed in the brain to date have only modeled time courses of total isotopic enrichment at individual carbon positions (positional enrichments), not taking advantage of the additional dynamic ¹³C isotopomer information available from fine-structure multiplets in ¹³C spectra. Here we introduce a new ¹³C metabolic modeling approach using the concept of bonded cumulative isotopomers, or bonded cumomers. The direct relationship between bonded cumomers and ¹³C multiplets enables fitting of the dynamic multiplet data. The potential of this new approach is demonstrated using Monte-Carlo simulations with a brain two-compartment neuronal-glial model. The precision of positional and cumomer approaches are compared for two different metabolic models (with and without glutamine dilution) and for different infusion protocols ([1,6-¹³C₂]glucose, [1,2-¹³C₂]acetate, and double infusion [1,6-¹³C₂]glucose?+?[1,2-¹³C₂]acetate). In all cases, the bonded cumomer approach gives better precision than the positional approach. In addition, of the three different infusion protocols considered here, the double infusion protocol combined with dynamic bonded cumomer modeling appears the most robust for precise determination of all fluxes in the model. The concepts and simulations introduced in the present study set the foundation for taking full advantage of the available dynamic ¹³C multiplet data in metabolic modeling.