Improving sequence-based fold recognition by using 3D model quality assessment

Motivation: The ability of a simple method (MODCHECK) to determine the sequence–structure compatibility of a set of structural models generated by fold recognition is tested in a thorough benchmark analysis. Four Model Quality Assessment Programs (MQAPs) were tested on 188 targets from the latest LiveBench-9 automated structure evaluation experiment. We systematically test and evaluate whether the MQAP methods can successfully detect native-likemodels. Results: We show that compared with the other three methods tested MODCHECK is the most reliable method for consistently performing the best top model selection and for ranking the models. In addition, we show that the choice of model similarity score used to assess a model's similarity to the experimental structure can influence the overall performance of these tools. Although these MQAP methods fail to improve the model selection performance for methods that already incorporate protein three dimension (3D) structural information, an improvement is observed for methods that are purely sequence-based, including the best profile–profile methods. This suggests that even the best sequence-based fold recognition methods can still be improved by taking into account the 3D structural information.

Rapid protein domain assignment from amino acid sequence using predicted secondary structure

The elucidation of the domain content of a given protein sequence in the absence of determined structure or significant sequence homology to known domains is an important problem in structural biology. Here we address how successfully the delineation of continuous domains can be accomplished in the absence of sequence homology using simple baseline methods, an existing prediction algorithm (Domain Guess by Size), and a newly developed method (DomSSEA). The study was undertaken with a view to measuring the usefulness of these prediction methods in terms of their application to fully automatic domain assignment. Thus, the sensitivity of each domain assignment method was measured by calculating the number of correctly assigned top scoring predictions. We have implemented a new continuous domain identification method using the alignment of predicted secondary structures of target sequences against observed secondary structures of chains with known domain boundaries as assigned by Class Architecture Topology Homology (CATH). Taking top predictions only, the success rate of the method in correctly assigning domain number to the representative chain set is 73.3%. The top prediction for domain number and location of domain boundaries was correct for 24% of the multidomain set (±20 residues). These results have been put into context in relation to the results obtained from the other prediction methods assessed

Profiling dehydrin gene sequence and physiological parameters in drought tolerant and susceptible spring wheat cultivars.

Physiological and yield traits such as stomatal conductance (mmol m-2s-1), Leaf relative water content (RWC %) and grain yield per plant were studied in a separate experiment. Results revealed that five out of sixteen cultivars viz. Anmol, Moomal, Sarsabz, Bhitai and Pavan, appeared to be relatively more drought tolerant. Based on morphophysiological results, studies were continued to look at these cultivars for drought tolerance at molecular level. Initially, four well recognized primers for dehydrin genes (DHNs) responsible for drought induction in T. durum L., T. aestivum L. and O. sativa L. were used for profiling gene sequence of sixteen wheat cultivars. The primers amplified the DHN genes variably like Primer WDHN13 (T. aestivum L.) amplified the DHN gene in only seven cultivars whereas primer TdDHN15 (T. durum L.) amplified all the sixteen cultivars with even different DNA banding patterns some showing second weaker DNA bands. Third primer TdDHN16 (T. durum L.) has shown entirely different PCR amplification prototype, specially showing two strong DNA bands while fourth primer RAB16C (O. sativa L.) failed to amplify DHN gene in any of the cultivars. Examination of DNA sequences revealed several interesting features. First, it identified the two exon/one intron structure of this gene (complete sequences were not shown), a feature not previously described in the two database cDNA sequences available from T. aestivum L. (gi|21850). Secondly, the analysis identified several single nucleotide polymorphisms (SNPs), positions in gene sequence. Although complete gene sequence was not obtained for all the cultivars, yet there were a total of 38 variable positions in exonic (coding region) sequence, from a total gene length of 453 nucleotides. Matrix of SNP shows these 37 positions with individual sequence at positions given for each of the 14 cultivars (sequence of two cultivars was not obtained) included in this analysis. It demonstrated a considerable diversity for this gene with only three cultivars i.e. TJ-83, Marvi and TD-1 being similar to the consensus sequence. All other cultivars showed a unique combination of SNPs. In order to prove a functional link between these polymorphisms and drought tolerance in wheat, it would be necessary to conduct a more detailed study involving directed mutation of this gene and DHN gene expression.

Complete sequence and molecular epidemiology of IncK epidemic plasmid encoding bla(CTX-M-14)

Antimicrobial drug resistance is a global challenge for the 21st century with the emergence of resistant bacterial strains worldwide. Transferable resistance to beta-lactam antimicrobial drugs, mediated by production of extended-spectrum beta-lactamases (ESBLs), is of particular concern. In 2004, an ESBL-carrying IncK plasmid (pCT) was isolated from cattle in the United Kingdom. The sequence was a 93,629-bp plasmid encoding a single antimicrobial drug resistance gene, bla(CTX-M-14). From this information, PCRs identifying novel features of pCT were designed and applied to isolates from several countries, showing that the plasmid has disseminated worldwide in bacteria from humans and animals. Complete DNA sequences can be used as a platform to develop rapid epidemiologic tools to identify and trace the spread of plasmids in clinically relevant pathogens, thus facilitating a better understanding of their distribution and ability to transfer between bacteria of humans and animals.

Sequence and phylogenetic analysis of viper venom serine proteases

Snakebites are a major neglected tropical disease responsible for as many as 95000 deaths every year worldwide. Viper venom serine proteases disrupt haemostasis of prey and victims by affecting various stages of the blood coagulation system. A better understanding of their sequence, structure, function and phylogenetic relationships will improve the knowledge on the pathological conditions and aid in the development of novel therapeutics for treating snakebites. A large dataset for all available viper venom serine proteases was developed and analysed to study various features of these enzymes. Despite the large number of venom serine protease sequences available, only a small proportion of these have been functionally characterised. Although, they share some of the common features such as a C-terminal extension, GWG motif and disulphide linkages, they vary widely between each other in features such as isoelectric points, potential N-glycosylation sites and functional characteristics. Some of the serine proteases contain substitutions for one or more of the critical residues in catalytic triad or primary specificity pockets. Phylogenetic analysis clustered all the sequences in three major groups. The sequences with substitutions in catalytic triad or specificity pocket clustered together in separate groups. Our study provides the most complete information on viper venom serine proteases to date and improves the current knowledge on the sequence, structure, function and phylogenetic relationships of these enzymes. This collective analysis of venom serine proteases will help in understanding the complexity of envenomation and potential therapeutic avenues.

Multi-target out-of-sequence data association: tracking using graphical models

When performing data fusion, one often measures where targets were and then wishes to deduce where targets currently are. There has been recent research on the processing of such out-of-sequence data. This research has culminated in the development of a number of algorithms for solving the associated tracking problem. This paper reviews these different approaches in a common Bayesian framework and proposes an architecture that orthogonalises the data association and out-of-sequence problems such that any combination of solutions to these two problems can be used together. The emphasis is not on advocating one approach over another on the basis of computational expense, but rather on understanding the relationships among the algorithms so that any approximations made are explicit. Results for a multi-sensor scenario involving out-of-sequence data association are used to illustrate the utility of this approach in a specific context.

Multi-target out-of-sequence data association

In data fusion systems, one often encounters measurements of past target locations and then wishes to deduce where the targets are currently located. Recent research on the processing of such out-of-sequence data has culminated in the development of a number of algorithms for solving the associated tracking problem. This paper reviews these different approaches in a common Bayesian framework and proposes an architecture that orthogonalises the data association and out-of-sequence problems such that any combination of solutions to these two problems can be used together. The emphasis is not on advocating one approach over another on the basis of computational expense, but rather on understanding the relationships between the algorithms so that any approximations made are explicit.

Reversible thermal transition of polydiacetylene based on KTTKS collagen sequence

Here we explore the physico-chemical properties of a peptide amphiphile obtained by chemical conjugation of the collagenstimulating peptide KTTKS with 10,12-pentacosadiynoic acid which photopolymerizes as a stable and extended polydiacetylene. We investigate the self-assembly of this new polymer and rationalize its peculiar behavior in terms of a thermal conformational transition. Surprisingly, this polymer shows a thermal transition associated with a non-cooperative increase in b-sheet content at high temperature.

Conformation and self-association of peptide amphiphiles based on the KTTKS collagen sequence

Studying peptide amphiphiles (PAs), we investigate the influence of alkyl chain length on the aggregation behavior of the collagen-derived peptide KTTKS with applications ranging from antiwrinkle cosmetic creams to potential uses in regenerative medicine. We have studied synthetic peptides amphiphiles C14− KTTKS (myristoyl Lys-Thr-Thr-Lys-Ser) and C18−KTTKS(stearoyl-Lys-Thr Thr-Lys-Ser) to investigate in detail their physicochemical properties. It is presumed that the hydrophobic chain in these self-assembling peptide amphiphiles enhances peptide permeation across the skin compared to KTTKS alone. Subsequently Cn−KTTKS should act as a prodrug and release the peptide by enzymatic cleavage. Our results should be useful in the further development of molecules with collagen-stimulating activity.

Distribution, gene sequence and expression in vivo of the plasmid encoded fimbrial antigen of Salmonella serotype Enteritidis

The pefA gene which encoded the serotype associated plasmid (SAP) mediated fimbrial major subunit antigen of Salmonella enterica serotype Typhimurium shared genetic identity with 128 of 706 salmonella isolates as demonstrated by dot (colony) hybridization. Seventy-seven of 113 isolates of Typhimurium and individual isolates of serotypes Bovis-morbificans, Cholerae-suis and Enteritidis phage type 9b hybridized pefA strongly, whereas 48 isolates of Enteritidis hybridized pefA weakly and one Enteritidis isolate of phage type 14b failed to hybridize. Individual isolates of 294 serotypes and 247 individual isolates of serotype Dublin did not hybridize pefA. Southern hybridization of plasmids extracted from Enteritidis demonstrated that the pefA gene probe hybridized strongly an atypical SAP of 80 kb in size harboured by one Enteritidis isolate of phage-type 9b, whereas the typical SAP of 58 kb in size harboured by 48 Enteritidis isolates hybridized weakly. One Enteritidis isolate of phage type 14b which failed to hybridize pefA in dot (colony) hybridization experiments was demonstrated to be plasmid free. A cosmid library of Enteritidis phage type 4 expressed in Escherichia coli K12 was screened by hybridization for the presence of pef sequences. Recombinant clones which were deduced to harbour the entire pef operon elaborated a PEF-like fimbrial structure at the cell surface. The PEF-like fimbrial antigen was purified from one cosmid clone and used in western blot experiments with sera from chickens infected with Enteritidis phage-type 4. Seroconversion to the fimbrial antigen was observed which indicated that the Enteritidis PEF-like fimbrial structure was expressed at some stage during infection. Nucleotide sequence analysis demonstrated that the pefA alleles of Typhimurium and Enteritidis phage-type 4 shared 76% DNA nucleotide and 82% deduced amino acid sequence identity.

Sequence analysis and distribution of an IS3-like insertion element isolated from Salmonella enteritidis

The nucleotide sequence of a 3 kb region immediately upstream of the sef operon operon of Salmonella enteritidis was determined. A 1230 base pair insertion sequence which shared sequence identity (> 75%) with members of the IS3 family was revealed. This element, designated IS1230, had almost identical (90% identity) terminal inverted repeats to Escherichia coli IS3 but unlike other IS3-like sequences lacked the two characteristic open reading frames which encode the putative transposase. S. enteritidis possessed only one copy of this insertion sequence although Southern hybridisation analysis of restriction digests of genomic DNA revealed another fragment located in a region different from the sef operon which hybridised weakly which suggested the presence of an IS1230 homologue. The distribution of IS1230 and IS1230-like elements was shown to be widespread amongst salmonellas and the patterns of restriction fragments which hybridised differed significantly between Salmonella serotypes and it is suggested that IS1230 has potential for development as a differential diagnostic tool.

Sequence heterogeneity of an mpb70 gene analogue in Mycobacterium kansasii

The protein antigen MPB70 is a major component of culture supernatants of Mycobacterium bovis and Is an active ingredient of bovine PPD used for skin-testing cattle for tuberculosis. we have shown that Mycobacterium kansasii possesses a similar gene that cross-reacts in a PCR test for M. bovis. Single strand conformational polymorphism analysis, and the DNA sequence of the PCR product, shows differences between M. kansasii strains, supporting the suggestion that M. kansasii is not a homogeneous species.

Sequence analysis and distribution in Salmonella enterica serovars of IS3-like elements

The genome of Salmonella enterica serovar Enteritidis was shown to possess three IS3-like insertion elements, designated IS1230A, B and C, and each was cloned and their respective deoxynucleotide sequences determined. Mutations in elements IS1230A and B resulted in frameshifts in the open reading frames that encoded a putative transposase to be inactive. IS1230C was truncated at nucleotide 774 relative to IS1230B and therefore did not possess the 3' terminal inverted repeat. The three IS1230 derivatives were closely related to each other based on nucleotide sequence similarity. IS1230A was located adjacent to the sef operon encoding SEF14 fimbriae located at minute 97 of the genome of S. Enteritidis. IS1230B was located adjacent to the umuDC operon at minute 42.5 on the genome, itself located near to one terminus of an 815-kb genome inversion of S. Enteritidis relative to S. Typhimurium. IS1230C was located next to attB, the bacteriophage P22 attachment site, and proB, encoding gamma-glutamyl phosphate reductase. A truncated 3' remnant of IS1230, designated IS1230T, was identified in a clinical isolate of S. Typhimurium DT193 strain 2391. This element was located next to attB adjacent to which were bacteriophage P22-like sequences. Southern hybridisation of total genomic DNA from eighteen phage types of S. Enteritidis and eighteen definitive types of S. Typhimurium showed similar, if not identical, restriction fragment profiles in the respective serovars when probed with IS1230A.

Phylogenetic relationships withinPelargonium sect. Peristera (Geraniaceae) inferred from nrDNA and cpDNA sequence comparisons

Phylogenetic analysis of nrDNA ITS and trnL (UAA) 5′ exon-trnF (GAA) chloroplast DNA sequences from 17 species ofPelargonium sect.Peristera, together with nine putative outgroups, suggests paraphyly for the section and a close relationship between the highly disjunct South African and Australian species of sect.Peristera. Representatives fromPelargonium sectt.Reniformia, Ligularia s. l. andIsopetalum (the St. Helena endemicP. cotyledonis) appear to be nested within thePeristera clade. The close relationship between the South African and AustralianPeristera is interpreted as being caused by long-range dispersal to Australia, probably as recent as the late Pliocene.

On the use of a time sequence of surface pressures in four-dimensional data assimilation

The possibility of using a time sequence of surface pressure observations in four-dimensional data assimilation is being investigated. It is shown that a linear multilevel quasi-geostrophic model can be updated successfully with surface data alone, provided the number of time levels are at least as many as the number of vertical levels. It is further demonstrated that current statistical analysis procedures are very inefficient to assimilate surface observations, and it is shown by numerical experiments that the vertical interpolation must be carried out using the structure of the most dominating baroclinic mode in order to obtain a satisfactory updating. Different possible ways towards finding a practical solution are being discussed.