Amino acid pairing preferences in parallel beta-sheets in proteins

Statistical approaches have been applied to examine amino acid pairing preferences within parallel beta-sheets. The main chain hydrogen bonding pattern in parallel beta-sheets means that, for each residue pair, only one of the residues is involved in main chain hydrogen bonding with the strand containing the partner residue. We call this the hydrogen bonded (HB) residue and the partner residue the non-hydrogen bonded (nHB) residue, and differentiate between the favorability of a pair and that of its reverse pair, e.g. Asn(HB)-Thr(nHB)versus Thr(HB)-Asn(nHB). Significantly (p < or = 0.000001) favoured pairings were rationalised using stereochemical arguments. For instance, Asn(HB)-Thr(nHB) and Arg(HB)-Thr(nHB) were favoured pairs, where the residues adopted favoured chi1 rotamer positions that allowed side-chain interactions to occur. In contrast, Thr(HB)-Asn(nHB) and Thr(HB)-Arg(nHB) were not significantly favoured, and could only form side-chain interactions if the residues involved adopted less favourable chi1 conformations. The favourability of hydrophobic pairs e.g. Ile(HB)-Ile(nHB), Val(HB)-Val(nHB) and Leu(HB)-Ile(nHB) was explained by the residues adopting their most preferred chi1 and chi2 conformations, which enabled them to form nested arrangements. Cysteine-cysteine pairs are significantly favoured, although these do not form intrasheet disulphide bridges. Interactions between positively and negatively charged residues were asymmetrically preferred: those with the negatively charged residue at the HB position were more favoured. This trend was accounted for by the presence of general electrostatic interactions, which, based on analysis of distances between charged atoms, were likely to be stronger when the negatively charged residue is the HB partner. The Arg(HB)-Asp(nHB) interaction was an exception to this trend and its favorability was rationalised by the formation of specific side-chain interactions. This research provides rules that could be applied to protein structure prediction, comparative modelling and protein engineering and design. The methods used to analyse the pairing preferences are automated and detailed results are available (http://www.rubic.rdg.ac.uk/betapairprefsparallel/).

Amino acid pairing preferences in parallel beta-sheets in proteins

Statistical approaches have been applied to examine amino acid pairing preferences within parallel beta-sheets. The main chain hydrogen bonding pattern in parallel beta-sheets means that, for each residue pair, only one of the residues is involved in main chain hydrogen bonding with the strand containing the partner residue. We call this the hydrogen bonded (HB) residue and the partner residue the non-hydrogen bonded (nHB) residue, and differentiate between the favourability of a pair and that of its reverse pair, e.g. Asn(HB)-Thr(nHB) versus Thr(HB)-Asn(nHB). Significantly (p <= 0.000001) favoured pairings were rationalised using stereochemical arguments. For instance, Asn(HB)-Thr(nHB) and Arg(HB)-Thr(nHB) were favoured pairs, where the residues adopted favoured chi(1) rotamer positions that allowed side-chain interactions to occur. In contrast, Thr(HB)-Asn(nHB) and Thr(HB)-Arg(nHB) were not significantly favoured, and could only form side-chain interactions if the residues involved adopted less favourable chi(1) conformations. The favourability of hydrophobic pairs e.g. Ile(HB)-Ile(nHB), Val(HB)-Val(nHB) and Leu(HB)-Ile(nHB) was explained by the residues adopting their most preferred chi(1) and chi(2) conformations, which enabled them to form nested arrangements. Cysteine-cysteine pairs are significantly favoured, although these do not form intrasheet disulphide bridges. Interactions between positively and negatively charged residues were asymmetrically preferred: those with the negatively charged residue at the HB position were more favoured. This trend was accounted for by the presence of general electrostatic interactions, which, based on analysis of distances between charged atoms, were likely to be stronger when the negatively charged residue is the HB partner. The Arg(HB)-Asp(nHB) interaction was an exception to this trend and its favourability was rationalised by the formation of specific side-chain interactions. This research provides rules that could be applied to protein structure prediction, comparative modelling and protein engineering and design. The methods used to analyse the pairing preferences are automated and detailed results are available (http:// www.rubic.rdg.ac.uk/betapairprefsparallel/). (c) 2005 Elsevier Ltd. All rights reserved.

The slow road to the eukaryotic genome

The eukaryotic genome is a mosaic of eubacterial and archaeal genes in addition to those unique to itself. The mosaic may have arisen as the result of two prokaryotes merging their genomes, or from genes acquired from an endosymbiont of eubacterial origin. A third possibility is that the eukaryotic genome arose from successive events of lateral gene transfer over long periods of time. This theory does not exclude the endosymbiont, but questions whether it is necessary to explain the peculiar set of eukaryotic genes. We use phylogenetic studies and reconstructions of ancestral first appearances of genes on the prokaryotic phylogeny to assess evidence for the lateral gene transfer scenario. We find that phylogenies advanced to support fusion can also arise from a succession of lateral gene transfer events. Our reconstructions of ancestral first appearances of genes reveal that the various genes that make up the eukaryotic mosaic arose at different times and in diverse lineages on the prokaryotic tree, and were not available in a single lineage. Successive events of lateral gene transfer can explain the unusual mosaic structure of the eukaryotic genome, with its content linked to the immediate adaptive value of the genes its acquired. Progress in understanding eukaryotes may come from identifying ancestral features such as the eukaryotic splicesome that could explain why this lineage invaded, or created, the eukaryoticniche.

Phyloclimatic modeling: Combining phylogenetics and bioclimatic modeling

We investigate the impact of past climates on plant diversification by tracking the "footprint" of climate change on a phylogenetic tree. Diversity within the cosmopolitan carnivorous plant genus Drosera (Droseraceae) is focused within Mediterranean climate regions. We explore whether this diversity is temporally linked to Mediterranean-type climatic shifts of the mid-Miocene and whether climate preferences are conservative over phylogenetic timescales. Phyloclimatic modeling combines environmental niche (bioclimatic) modeling with phylogenetics in order to study evolutionary patterns in relation to climate change. We present the largest and most complete such example to date using Drosera. The bioclimatic models of extant species demonstrate clear phylogenetic patterns; this is particularly evident for the tuberous sundews from southwestern Australia (subgenus Ergaleium). We employ a method for establishing confidence intervals of node ages on a phylogeny using replicates from a Bayesian phylogenetic analysis. This chronogram shows that many clades, including subgenus Ergaleium and section Bryastrum, diversified during the establishment of the Mediterranean-type climate. Ancestral reconstructions of bioclimatic models demonstrate a pattern of preference for this climate type within these groups. Ancestral bioclimatic models are projected into palaeo-climate reconstructions for the time periods indicated by the chronogram. We present two such examples that each generate plausible estimates of ancestral lineage distribution, which are similar to their current distributions. This is the first study to attempt bioclimatic projections on evolutionary time scales. The sundews appear to have diversified in response to local climate development. Some groups are specialized for Mediterranean climates, others show wide-ranging generalism. This demonstrates that Phyloclimatic modeling could be repeated for other plant groups and is fundamental to the understanding of evolutionary responses to climate change.

A phyloclimatic study of Cyclamen

Background: The impact of global climate change on plant distribution, speciation and extinction is of current concern. Examining species climatic preferences via bioclimatic niche modelling is a key tool to study this impact. There is an established link between bioclimatic niche models and phylogenetic diversification. A next step is to examine future distribution predictions from a phylogenetic perspective. We present such a study using Cyclamen (Myrsinaceae), a group which demonstrates morphological and phenological adaptations to its seasonal Mediterranean-type climate. How will the predicted climate change affect future distribution of this popular genus of garden plants? Results: We demonstrate phylogenetic structure for some climatic characteristics, and show that most Cyclamen have distinct climatic niches, with the exception of several wide-ranging, geographically expansive, species. We reconstruct climate preferences for hypothetical ancestral Cyclamen. The ancestral Cyclamen lineage has a preference for the seasonal Mediterranean climate characteristic of dry summers and wet winters. Future bioclimatic niches, based on BIOCLIM and Maxent models, are examined with reference to a future climate scenario for the 2050s. Over the next 50 years we predict a northward shift in the area of climatic suitability, with many areas of current distribution becoming climatically unsuitable. The area of climatic suitability for every Cyclamen species is predicted to decrease. For many species, there may be no areas with a suitable climate regardless of dispersal ability, these species are considered to be at high risk of extinction. This risk is examined from a phylogenetic perspective. Conclusion: Examining bioclimatic niches from a phylogenetic perspective permits novel interpretations of these models. In particular, reconstruction of ancestral niches can provide testable hypothesis about the historical development of lineages. In the future we can expect a northwards shift in climatic suitability for the genus Cyclamen. If this proves to be the case then dispersal is the best chance of survival, which seems highly unlikely for ant-dispersed Cyclamen. Human-assisted establishment of Cyclamen species well outside their native ranges offers hope and could provide the only means of dispersal to potentially suitable future environments. Even without human intervention the phylogenetic perspective demonstrates that major lineages could survive climate change even if many species are lost.

Comparative expression analysis of transcription factor genes in the endostyle of invertebrate chordates

The endostyle of invertebrate chordates is a pharyngeal organ that is thought to be homologous with the follicular thyroid of vertebrates. Although thyroid-like features such as iodine-concentrating and peroxidase activities are located in the dorsolateral part of both ascidian and amphioxus endostyles, the structural organization and numbers of functional units are different. To estimate phylogenetic relationships of each functional zone with special reference to the evolution of the thyroid, we have investigated, in ascidian and amphioxus, the expression patterns of thyroid-related transcription factors such as TTF-2/MoxE4 and Pax2/5/8, as well as the forkhead transcription factors FoxQ1 and FoxA. Comparative gene expression analyses depicted an overall similarity between ascidians and amphioxus endostyles, while differences in expression patterns of these genes might be specifically related to the addition or elimination of a pair of glandular zones. Expressions of Ci-FoxE and BbFoxE4 suggest that the ancestral FoxE class might have been recruited for the formation of thyroid-like region in a possible common ancestor of chordates. Furthermore, coexpression of FoxE4, Pax2/5/8, and TPO in the dorsolateral part of both ascidian and amphioxus endostyles suggests that genetic basis of the thyroid function was already in place before the vertebrate lineage. (c) 2005 Wiley-Liss, Inc.

Limpets break Dollo's law

A new molecular phylogeny of the limpet molluscs (Calyptraeidae) reveals that coiled shells have independently re-evolved at least once in this family, which is a violation of Dollo's Law that complex ancestral states, once lost, are never reacquired. Reacquisition of the coiled ancestral state is remarkable in that uncoiled shells have been the most recent ancestral state for 20 million-100 million years. Adult coiling might have reevolved by the mechanism of prolonging the period during which genes for coiling are expressed in larvae. This and other developmental mechanisms could provide general routes for maintaining the potential to produce traits lost in distant ancestors.

The life history legacy of evolutionary body size change in carnivores

We estimate the body sizes of direct ancestors of extant carnivores, and examine selected aspects of life history as a function not only of species' current size, but also of recent changes in size. Carnivore species that have undergone marked recent evolutionary size change show life history characteristics typically associated with species closer to the ancestral body size. Thus, phyletic giants tend to mature earlier and have larger litters of smaller offspring at shorter intervals than do species of the same body size that are not phyletic giants. Phyletic dwarfs, by contrast, have slower life histories than nondwarf species of the same body size. We discuss two possible mechanisms for the legacy of recent size change: lag (in which life history variables cannot evolve as quickly as body size, leading to species having the 'wrong' life history for their body size) and body size optimization (in which life history and hence body size evolve in response to changes in energy availability); at present, we cannot distinguish between these alternatives. Our finding that recent body size changes help explain residual variation around life history allometries shows that a more dynamic view of character change enables comparative studies to make more precise predictions about species traits in the context of their evolutionary background.

Dispersal of NK homeobox gene clusters in amphioxus and humans

The Drosophila melanogaster genome has six physically clustered NK-related homeobox genes in just 180 kb. Here we show that the NK homeobox gene cluster was an ancient feature of bilaterian animal genomes, but has been secondarily split in chordate ancestry. The NK homeobox gene clusters of amphioxus and vertebrates are each split and dispersed at two equivalent intergenic positions. From the ancestral NK gene cluster, only the Tlx-Lbx and NK3-NK4 linkages have been retained in chordates. This evolutionary pattern is in marked contrast to the Hox and ParaHox gene clusters, which are compact in amphioxus and vertebrates, but have been disrupted in Drosophila.

Characterization of wheat puroindoline proteins

Puroindoline proteins were purified from selected UK-grown hexaploid wheats. Their identities were confirmed on the basis of capillary electrophoresis mobilities, relative molecular mass and N-terminal amino acid sequencing. Only one form of puroindoline-a protein was found in those varieties, regardless of endosperm texture. Three allelic forms of puroindoline-b protein were identified. Nucleotide sequencing of cDNA produced by RT-PCR of isolated mRNA indicated that these were the 'wild-type', found in soft wheats, puroindoline-b containing a Gly -> Ser amino acid substitution (position 46) and puroindoline-b containing a Trp -> Arg substitution (position 44). The latter two were found in hard wheats. Microheterogeneity, due to short extensions and/or truncations at the N-terminus and C-terminus, was detected for both puroindoline-a and puroindoline-b. The type of microheterogeneity observed was more consistent for puroindoline-a than for puroindoline-b, and may arise through slightly different post-translational processing pathways. A puroindoline-b allele corresponding to a Leu -> Pro substitution (position 60) was identified from the cDNA sequence of the hard variety Chablis, but no mature puroindoline-b protein was found in this or two other European varieties known to possess this puroindoline-b allele. Wheats possessing the puroindoline-b proteins with point mutations appeared to contain lower amounts of puroindoline protein. Such wheats have a hard endosperm texture, as do wheats from which puroindoline-a or puroindoline-b are absent. Our results suggest that point mutations in puroindoline-b genes may confer hard endosperm texture through accumulation of allelic forms of puroindoline-b proteins with altered functional properties and/or through lower amounts of puroindoline proteins.

Proteomic analysis of arginine adducts on glyoxal-modified ribonuclease

Accumulation of advanced glycation end-products (AGEs) on proteins is associated with the development of diabetic complications. Although the overall extent of modification of protein by AGEs is limited, localization of these modifications at a few critical sites might have a significant effect on protein structure and function. In the present study, we describe the sites of modification of RNase by glyoxal under physiological conditions. Arg(39) and Arg(85), which are closest to the active site of the enzyme, were identified as the primary sites of formation of the glyoxal-derived dihydroxyimidazolidine and hydroimidazolone adducts. Lower amounts of modification were detected at Arg(10), while Arg(33) appeared to be unmodified. We conclude that dihydroxyimidazolidine adducts are the primary products of modification of protein by glyoxal, that Arg(39) and Arg(85) are the primary sites of modification of RNase by glyoxal, and that modification of arginine residues during Maillard reactions of proteins is a highly selective process.

Mitochondrial and chloroplast DNA-based phylogeny of Pelargonium (Geraniaceae)

Overall phylogenetic relationships within the genus Pelargonium (Geraniaceae) were inferred based on DNA sequences from mitochondrial(mt)-encoded nad1 b/c exons and from chloroplast(cp)-encoded trnL (UAA) 5' exon-trnF (GAA) exon regions using two species of Geranium and Sarcocaulon vanderetiae as outgroups. The group II intron between nad1 exons b and c was found to be absent from the Pelargonium, Geranium, and Sarcocaulon sequences presented here as well as from Erodium, which is the first recorded loss of this intron in angiosperms. Separate phylogenetic analyses of the mtDNA and cpDNA data sets produced largely congruent topologies, indicating linkage between mitochondrial and chloroplast genome inheritance. Simultaneous analysis of the combined data sets yielded a well-resolved topology with high clade support exhibiting a basic split into small and large chromosome species, the first group containing two lineages and the latter three. One large chromosome lineage (x = 11) comprises species from sections Myrrhidium and Chorisma and is sister to a lineage comprising P. mutans (x = 11) and species from section Jenkinsonia (x = 9). Sister to these two lineages is a lineage comprising species from sections Ciconium (x = 9) and Subsucculentia (x = 10). Cladistic evaluation of this pattern suggests that x = 11 is the ancestral basic chromosome number for the genus.

Legume type and temperature effects on the toxicity of insecticide to the genus Callosobruchus (Coleoptera: Bruchidae)

Three bruchid pest species, Callosobruchus maculatus, Callosobruchus chinensis and Callosobruchus rhodesianus, were studied for their response to insecticide toxicity taking into account the separate and interactive effects of temperature and pre-adult food. The food types used were cowpea (Vigna unguiculata) and mungbean (Vigna radiata). Callosobruchus maculatus was the most tolerant to malathion and the least affected by temperature change while C. rhodesianus was the least tolerant. Over a 4 C range (23, 25, 27 C), there was generally a significant impact of temperature on the tolerance of the three species to the insecticide. The food type on which the insects developed influenced considerably the degree of insecticide tolerance. Callosobruchus maculatus and C. chinensis populations reared onmungbean had higher tolerance to malathion than their counterparts reared on cowpea, but the opposite was observed in C. rhodesianus populations. The food influence in this study suggested an ancestral cause or fitness cost depending on the species. The interaction of food-by-temperature had no significant effect on malathion toxicity to this genus. Correlation analysis showed C. chinensis to be relatively less sensitive to insecticide concentration over the range studied compared with the other two species.

Prevalence of mutations within the quinolone resistance-determining region of gyrA, gyrB, parC, and parE and association with antibiotic resistance in quinolone-resistant Salmonella enterica

Salmonella enterica isolates (n = 182) were examined for mutations in the quinolone resistance-determining region of gyrA, gyrB, parC, and parE. The frequency, location, and type of GyrA substitution varied with the serovar. Mutations were found in parC that encoded Thr57-Ser, Thr66-Ile, and Ser80-Arg substitutions. Mutations in the gyrB quinolone resistance-determining region were located at codon Tyr420-Cys or Arg437-Len. Novel mutations were also found in parE encoding Glu453-Gly, His461-Tyr, Ala498-Thr, Val512-Gly, and Ser518-Cys. Although it is counterintuitive, isolates with a mutation in both gyrA and parC were more susceptible to ciprofloxacin than were isolates with a mutation in gyrA alone.

Self-assembled arginine-coated peptide nanosheets in water

The surfactant-like peptide (Ala)6(Arg) is found to self-assemble into 3 nm-thick sheets in aqueous solution. Scanning transmission electron microscopy measurements of mass per unit area indicate a layer structure based on antiparallel dimers. At higher concentration the sheets wrap into unprecedented ultrathin helical ribbon and nanotube architectures.