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.

Replication enhancer elements within the open reading frame of tick-borne encephalitis virus and their evolution within the Flavivirus genus

We provide experimental evidence of a replication enhancer element (REE) within the capsid gene of tick-borne encephalitis virus (TBEV, genus Flavivirus). Thermodynamic and phylogenetic analyses predicted that the REE folds as a long stable stem–loop (designated SL6), conserved among all tick-borne flaviviruses (TBFV). Homologous sequences and potential base pairing were found in the corresponding regions of mosquito-borne flaviviruses, but not in more genetically distant flaviviruses. To investigate the role of SL6, nucleotide substitutions were introduced which changed a conserved hexanucleotide motif, the conformation of the terminal loop and the base-paired dsRNA stacking. Substitutions were made within a TBEV reverse genetic system and recovered mutants were compared for plaque morphology, single-step replication kinetics and cytopathic effect. The greatest phenotypic changes were observed in mutants with a destabilized stem. Point mutations in the conserved hexanucleotide motif of the terminal loop caused moderate virus attenuation. However, all mutants eventually reached the titre of wild-type virus late post-infection. Thus, although not essential for growth in tissue culture, the SL6 REE acts to up-regulate virus replication. We hypothesize that this modulatory role may be important for TBEV survival in nature, where the virus circulates by non-viraemic transmission between infected and non-infected ticks, during co-feeding on local rodents.

Synthesis and phase behaviour of a homologous series of polymethacrylate-based side-chain liquid crystal polymers

The homologous series of side chain liquid crystal polymers, the poly[x-(4-methoxyazobenzene- 40-oxy)alkyl methacrylate]s, has been prepared in which the length of the flexible alkyl spacer has been varied from 3 to 11 methylene units. All the polymers exhibit liquid crystalline behaviour. The propyl and butyl members show exclusively nematic behaviour. The pentyl, hexyl, octyl and decyl members show a nematic and a smectic A phase while the heptyl, nonyl and undecyl homologues exhibit only a smectic A phase. The smectic A phase has been studied using X-ray diffraction and assigned as a smectic A1 phase in which the side chains are fully overlapped and the backbones are confined to lie between the smectic layers. For the nonyl member an incommensurate smectic phase is observed. The dependence of the transition temperatures on the length of the flexible spacer is understood in terms of the average shapes of the side chains.

Computing the Cassels–Tate pairing on the 3-Selmer group of an elliptic curve

We extend the method of Cassels for computing the Cassels-Tate pairing on the 2-Selmer group of an elliptic curve, to the case of 3-Selmer groups. This requires significant modifications to both the local and global parts of the calculation. Our method is practical in sufficiently small examples, and can be used to improve the upper bound for the rank of an elliptic curve obtained by 3-descent.

Rapid word learning by 15-month-olds under tightly-controlled conditions

Infants (12 to 17 months) were taught 2 novel words for 2 images of novel objects, by pairing isolated auditory labels with to-be-associated images. Comprehension was tested using a preferential looking task in which the infant was presented with both images together with an isolated auditory label. The auditory label usually, but not always, matched one of the images. Infants looked preferentially at images that matched the auditory stimulus. The experiment controlled within-subjects for both side bias and preference for previously named items. Infants showed learning after 12 presentations of the new words. Evidence is presented that, in certain circumstances, the duration of longest look at a target may be a more robust measure of target preference than overall looking time. The experiment provides support for previous demonstrations of rapid word learning by pre-vocabulary spurt children, and offers some methodological improvements to the preferential looking task.

Soft topographic map for clustering and classification of bacteria

In this work a new method for clustering and building a topographic representation of a bacteria taxonomy is presented. The method is based on the analysis of stable parts of the genome, the so-called “housekeeping genes”. The proposed method generates topographic maps of the bacteria taxonomy, where relations among different type strains can be visually inspected and verified. Two well known DNA alignement algorithms are applied to the genomic sequences. Topographic maps are optimized to represent the similarity among the sequences according to their evolutionary distances. The experimental analysis is carried out on 147 type strains of the Gammaprotebacteria class by means of the 16S rRNA housekeeping gene. Complete sequences of the gene have been retrieved from the NCBI public database. In the experimental tests the maps show clusters of homologous type strains and present some singular cases potentially due to incorrect classification or erroneous annotations in the database.

Gene conversion and evolution of Xq28 duplicons involved in recurring inversions causing severe hemophilia A

Inversions breaking the 1041 bp int1h-1 or the 9.5-kb int22h-1 sequence of the F8 gene cause hemophilia A in 1/30,000 males. These inversions are due to homologous recombination between the above sequences and their inverted copies on the same DNA molecule, respectively, int1h-2 and int22h-2 or int22h-3. We find that (1) int1h and int22h duplicated more than 25 million years ago; (2) the identity of the copies (>99%) of these sequences in humans and other primates is due to gene conversion; (3) gene conversion is most frequent in the internal regions of int22h; (4) breakpoints of int22h-related inversions also tend to involve the internal regions of int22h; (5) sequence variations in a sample of human X chromosomes defined eight haplotypes of int22h-1 and 27 of int22h-2 plus int22h-3; (6) the latter two sequences, which lie, respectively, 500 and 600 kb telomeric to int22h-1 are five-fold more identical when in cis than when in trans, thus suggesting that gene conversion may be predominantly intrachromosomal; (7) int1h, int22h, and flanking sequences evolved at a rate of about 0.1% substitutions per million years during the divergence between humans and other primates, except for int1h during the human-chimpanzee divergence, when its rate of evolution was significantly lower. This is reminiscent of the slower evolution of palindrome arms in the male specific regions of the Y chromosome and we propose, as an explanation, that intrachromosomal gene conversion and cosegregation of the duplicated regions favors retention of the ancestral sequence and thus reduces the evolution rate.

Nuclear magnetic resonance structure shows that the severe acute respiratory syndrome coronavirus-unique domain contains a macrodomain fold

The nuclear magnetic resonance (NMR) structure of a central segment of the previously annotated severe acute respiratory syndrome (SARS)-unique domain (SUD-M, for "middle of the SARS-unique domain") in SARS coronavirus (SARS-CoV) nonstructural protein 3 (nsp3) has been determined. SUD-M(513-651) exhibits a macrodomain fold containing the nsp3 residues 528 to 648, and there is a flexibly extended N-terminal tail with the residues 513 to 527 and a C-terminal flexible tail of residues 649 to 651. As a follow-up to this initial result, we also solved the structure of a construct representing only the globular domain of residues 527 to 651 [SUD-M(527-651)]. NMR chemical shift perturbation experiments showed that SUD-M(527-651) binds single-stranded poly(A) and identified the contact area with this RNA on the protein surface, and electrophoretic mobility shift assays then confirmed that SUD-M has higher affinity for purine bases than for pyrimidine bases. In a further search for clues to the function, we found that SUD-M(527-651) has the closest three-dimensional structure homology with another domain of nsp3, the ADP-ribose-1 ''-phosphatase nsp3b, although the two proteins share only 5% sequence identity in the homologous sequence regions. SUD-M(527-651) also shows three-dimensional structure homology with several helicases and nucleoside triphosphate-binding proteins, but it does not contain the motifs of catalytic residues found in these structural homologues. The combined results from NMR screening of potential substrates and the structure-based homology studies now form a basis for more focused investigations on the role of the SARS-unique domain in viral infection.

EfeO-cupredoxins: major new members of the cupredoxin superfamily with roles in bacterial iron transport

The EfeUOB system of Escherichia coli is a tripartite, low pH, ferrous iron transporter. It resembles the high-affinity iron transporter (Ftr1p-Fet3p) of yeast in that EfeU is homologous to Ftr1p, an integral-membrane iron-permease. However, EfeUOB lacks an equivalent of the Fet3p component—the multicopper oxidase with three cupredoxin-like domains. EfeO and EfeB are periplasmic but their precise roles are unclear. EfeO consists primarily of a C-terminal peptidase-M75 domain with a conserved ‘HxxE’ motif potentially involved in metal binding. The smaller N-terminal domain (EfeO-N) is predicted to be cupredoxin (Cup) like, suggesting a previously unrecognised similarity between EfeO and Fet3p. Our structural modelling of the E. coli EfeO Cup domain identifies two potential metal-binding sites. Site I is predicted to bind Cu2+ using three conserved residues (C41 and 103, and E66) and M101. Of these, only one (C103) is conserved in classical cupredoxins where it also acts as a Cu ligand. Site II most probably binds Fe3+ and consists of four well conserved surface Glu residues. Phylogenetic analysis indicates that the EfeO-Cup domains form a novel Cup family, designated the ‘EfeO-Cup’ family. Structural modelling of two other representative EfeO-Cup domains indicates that different subfamilies employ distinct ligand sets at their proposed metal-binding sites. The ~100 efeO homologues in the bacterial sequence databases are all associated with various iron-transport related genes indicating a common role for EfeO-Cup proteins in iron transport, supporting a new copper-iron connection in biology.

NMR structure shows that the SARS-unique domain contains a macrodomain fold

The NMR structure of a central segment of the previously annotated "SARS-unique domain" (SUD-M; "middle of the SARS-unique domain") in the SARS coronavirus (SARS-CoV) non-structural protein 3 (nsp3) has been determined. SUD-M(513-651) exhibits a macrodomain fold containing the nsp3-residues 528-648, and there is a flexibly extended N-terminal tail with the residues 513-527 and a C-terminal flexible tail of residues 649-651. As a follow-up to this initial result, we also solved the structure of a construct representing only the globular domain of residues 527-651 [SUD-M(527-651)]. NMR chemical shift perturbation experiments showed that SUD-M(527-651) binds single-stranded poly-A and identified the contact area with this RNA on the protein surface, and electrophoretic mobility shift assays then confirmed that SUD-M has higher affinity for purine bases than for pyrimidine bases. In further search for clues to the function, we found that SUD-M(527-651) has the closest three-dimensional structure homology with another domain of nsp3, the ADP-ribose-1''-phosphatase nsp3b, although the two proteins share only 5% sequence identity in the homologous sequence regions. SUD-M(527-651) also shows 3D structure homology with several helicases and NTP-binding proteins, but it does not contain the motifs of catalytic residues found in these structural homologues. The combined results from NMR screening of potential substrates and the structure-based homology studies now form a basis for more focused investigations on the role of the SARS-unique domain in viral infection.

EfeUOB (YcdNOB) is a tripartite, acid-induced and CpxAR-regulated, low-pH Fe2+ transporter that is cryptic in Escherichia coli K-12 but functional in E-coli O157 : H7

Escherichia coli possesses iron transporters specific for either Fe2+ or Fe3+. Although Fe2+ is far more soluble than Fe3+, it rapidly oxidizes aerobically at pH >= 7. Thus, FeoAB, the major Fe2+ transporter of E. coli, operates anaerobically. However, Fe2+ remains stable aerobically under acidic conditions, although a low-pH Fe2+ importer has not been previously identified. Here we show that ycdNOB (efeUOB) specifies the first such transporter. efeUOB is repressed at high pH by CpxAR, and is Fe2+-Fur repressed. EfeU is homologous to the high-affinity iron permease, Ftr1p, of Saccharomyces cerevisiae and other fungi. EfeO is periplasmic with a cupredoxin N-terminal domain; EfeB is also periplasmic and is haem peroxidase-like. All three Efe proteins are required for Efe function. The efeU gene of E. coli K-12 is cryptic due to a frameshift mutation - repair of the single-base-pair deletion generates a functional EfeUOB system. In contrast, the efeUOB operon of the enterohaemorrhagic strain, O157:1147, lacks any frameshift and is functional. A 'wild-type' K-12 strain bearing a functional EfeUOB displays a major growth advantage under aerobic, low-pH, low-iron conditions when a competing metal is provided. Fe-55 transport assays confirm the ferrous iron specificity of EfeUOB.

Origin and evolution of 3'UTR of flaviviruses: long direct repeats as a basis for the formation of secondary structures and their significance for virus transmission

The 3' untranslated regions (3'UTRs) of flaviviruses are reviewed and analyzed in relation to short sequences conserved as direct repeats (DRs). Previously, alignments of the 3'UTRs have been constructed for three of the four recognized flavivirus groups, namely mosquito-borne, tick-borne, and nonclassified flaviviruses (MBFV, TBFV, and NCFV, respectively). This revealed (1) six long repeat sequences (LRSs) in the 3'UTR and open-reading frame (ORF) of the TBFV, (2) duplication of the 3'UTR of the NCFV by intramolecular recombination, and (3) the possibility of a common origin for all DRs within the MBFV. We have now extended this analysis and review it in the context of all previous published analyses. This has been achieved by constructing a robust alignment between all flaviviruses using the published DRs and secondary RNA structures as "anchors" to reveal additional homologies along the 3'UTR. This approach identified nucleotide regions within the MBFV, NKV (no-known vector viruses), and NCFV 3'UTRs that are homologous to different LRSs in the TBFV 3'UTR and ORF. The analysis revealed that some of the DRs and secondary RNA structures described individually within each flavivirus group share common evolutionary origins. The 3'UTR of flaviviruses, and possibly the ORF, therefore probably evolved through multiple duplication of an RNA domain, homologous to the LRS previously identified only in the TBFV. The short DRs in all virus groups appear to represent the evolutionary remnants of these domains rather than resulting from new duplications. The relevance of these flavivirus DRs to evolution, diversity, 3'UTR enhancer function, and virus transmission is reviewed.

Origin and evolution of a myxozoan worm

Buddenbrockia pluinatellae is an active, muscular, worm-shaped parasite of freshwater bryozoans. This rare and enigmatic animal has been assigned to the Myxozoa on the basis of 18S ribosomal DNA sequences and the presence of malacosporean spores. Here we report cloning of four homologous protein-coding genes from Buddenbrockia worms, the putatively conspecific sac-shaped parasite originally described as Tetracapsula bryozoides and the related sac-shaped parasite Tetracapsuloides bryosalmonae, the causative agent of proliferative kidney disease in salmonid fish. Analyses are consistent with the hypothesis that Buddenbrockia is indeed a malacosporean myxozoan, but do not provide support for conspecificity with either T. bryozoides or T. bryosalmonae. Implications for the evolution of worm-like body plans in the Myxozoa are discussed.

Ribonucleocapsid formation of severe acute respiratory syndrome coronavirus through molecular action of the N-terminal domain of N protein

Conserved among all coronaviruses are four structural proteins: the matrix (M), small envelope (E), and spike (S) proteins that are embedded in the viral membrane and the nucleocapsid phosphoprotein (N), which exists in a ribonucleoprotein complex in the lumen. The N-terminal domain of coronaviral N proteins (N-NTD) provides a scaffold for RNA binding, while the C-terminal domain (N-CTD) mainly acts as oligomerization modules during assembly. The C terminus of the N protein anchors it to the viral membrane by associating with M protein. We characterized the structures of N-NTD from severe acute respiratory syndrome coronavirus (SARS-CoV) in two crystal forms, at 1.17 A (monoclinic) and at 1.85 A (cubic), respectively, resolved by molecular replacement using the homologous avian infectious bronchitis virus (IBV) structure. Flexible loops in the solution structure of SARS-CoV N-NTD are now shown to be well ordered around the beta-sheet core. The functionally important positively charged beta-hairpin protrudes out of the core, is oriented similarly to that in the IBV N-NTD, and is involved in crystal packing in the monoclinic form. In the cubic form, the monomers form trimeric units that stack in a helical array. Comparison of crystal packing of SARS-CoV and IBV N-NTDs suggests a common mode of RNA recognition, but they probably associate differently in vivo during the formation of the ribonucleoprotein complex. Electrostatic potential distribution on the surface of homology models of related coronaviral N-NTDs suggests that they use different modes of both RNA recognition and oligomeric assembly, perhaps explaining why their nucleocapsids have different morphologies.