Visual cavity analysis in molecular simulations

Molecular surfaces provide a useful mean for analyzing interactions between biomolecules; such as identification and characterization of ligand binding sites to a host macromolecule. We present a novel technique, which extracts potential binding sites, represented by cavities, and characterize them by 3D graphs and by amino acids. The binding sites are extracted using an implicit function sampling and graph algorithms. We propose an advanced cavity exploration technique based on the graph parameters and associated amino acids. Additionally, we interactively visualize the graphs in the context of the molecular surface. We apply our method to the analysis of MD simulations of Proteinase 3, where we verify the previously described cavities and suggest a new potential cavity to be studied

Somatostatin receptors 2 and 5 are preferentially expressed in proliferating endothelium

Angiogenesis is characterised by activation, migration and proliferation of endothelial cells and is central to the pathology of cancer, cardiovascular disease and chronic inflammation. Somatostatin is an inhibitory polypeptide that acts through five receptors (sst 1, 2, 3, 4, 5). Sst has previously been reported in endothelium, but their role remains obscure. Here, we report the expression of sst in human umbilical vein endothelial cells (HUVECs) in vitro, during proliferation and quiescence. A protocol for culturing proliferating and quiescent HUVECs was established, and verified by analysing cell cycle distribution in propidium-iodide-stained samples using flow cytometry. Sst mRNA was then quantified in nine proliferating and quiescent HUVEC lines using quantitative reverse transcriptase–polymerase chain reaction. Sst 2 and 5 were preferentially expressed in proliferating HUVECs. All samples were negative for sst 4. Sst 1 and 3 expression and cell cycle progression were unrelated. Immunostaining for sst 2 and 5 showed positivity in proliferating but not quiescent cells, confirming sst 2 and 5 protein expression. Inhibition of proliferating cells with somatostatin analogues Octreotide and SOM230, which have sst 5 activity, was found (Octreotide 10−10–10−6 M: 48.5–70.2% inhibition; SOM230 10−9–10−6 M: 44.9–65.4% inhibition) in a dose-dependent manner, suggesting that sst 5 may have functional activity in proliferation. Dynamic changes in sst 2 and 5 expression during the cell cycle and the inhibition of proliferation with specific analogues suggest that these receptors may have a role in angiogenesis

Using a Human Challenge Model of Infection to Measure Vaccine Efficacy: A Randomised, Controlled Trial Comparing the Typhoid Vaccines M01ZH09 with Placebo and Ty21a

Background Typhoid persists as a major cause of global morbidity. While several licensed vaccines to prevent typhoid are available, they are of only moderate efficacy and unsuitable for use in children less than two years of age. Development of new efficacious vaccines is complicated by the human host-restriction of Salmonella enterica serovar Typhi (S. Typhi) and lack of clear correlates of protection. In this study, we aimed to evaluate the protective efficacy of a single dose of the oral vaccine candidate, M01ZH09, in susceptible volunteers by direct typhoid challenge. Methods and Findings We performed a randomised, double-blind, placebo-controlled trial in healthy adult participants at a single centre in Oxford (UK). Participants were allocated to receive one dose of double-blinded M01ZH09 or placebo or 3-doses of open-label Ty21a. Twenty-eight days after vaccination, participants were challenged with 104CFU S. Typhi Quailes strain. The efficacy of M01ZH09 compared with placebo (primary outcome) was assessed as the percentage of participants reaching pre-defined endpoints constituting typhoid diagnosis (fever and/or bacteraemia) during the 14 days after challenge. Ninety-nine participants were randomised to receive M01ZH09 (n = 33), placebo (n = 33) or 3-doses of Ty21a (n = 33). After challenge, typhoid was diagnosed in 18/31 (58.1% [95% CI 39.1 to 75.5]) M01ZH09, 20/30 (66.7% [47.2 to 87.2]) placebo, and 13/30 (43.3% [25.5 to 62.6]) Ty21a vaccine recipients. Vaccine efficacy (VE) for one dose of M01ZH09 was 13% [95% CI -29 to 41] and 35% [-5 to 60] for 3-doses of Ty21a. Retrospective multivariable analyses demonstrated that pre-existing anti-Vi antibody significantly reduced susceptibility to infection after challenge; a 1 log increase in anti-Vi IgG resulting in a 71% decrease in the hazard ratio of typhoid diagnosis ([95% CI 30 to 88%], p = 0.006) during the 14 day challenge period. Limitations to the study included the requirement to limit the challenge period prior to treatment to 2 weeks, the intensity of the study procedures and the high challenge dose used resulting in a stringent model. Conclusions Despite successfully demonstrating the use of a human challenge study to directly evaluate vaccine efficacy, a single-dose M01ZH09 failed to demonstrate significant protection after challenge with virulent Salmonella Typhi in this model. Anti-Vi antibody detected prior to vaccination played a major role in outcome after challenge

Brain Expressed microRNAs Implicated in Schizophrenia Etiology

BACKGROUND: Protein encoding genes have long been the major targets for research in schizophrenia genetics. However, with the identification of regulatory microRNAs (miRNAs) as important in brain development and function, miRNAs genes have emerged as candidates for schizophrenia-associated genetic factors. Indeed, the growing understanding of the regulatory properties and pleiotropic effects that miRNA have on molecular and cellular mechanisms, suggests that alterations in the interactions between miRNAs and their mRNA targets may contribute to phenotypic variation. METHODOLOGY/PRINCIPAL FINDINGS: We have studied the association between schizophrenia and genetic variants of miRNA genes associated with brain-expression using a case-control study design on three Scandinavian samples. Eighteen known SNPs within or near brain-expressed miRNAs in three samples (Danish, Swedish and Norwegian: 420/163/257 schizophrenia patients and 1006/177/293 control subjects), were analyzed. Subsequently, joint analysis of the three samples was performed on SNPs showing marginal association. Two SNPs rs17578796 and rs1700 in hsa-mir-206 (mir-206) and hsa-mit-198 (mir-198) showed nominal significant allelic association to schizophrenia in the Danish and Norwegian sample respectively (P = 0.0021 & p = 0.038), of which only rs17578796 was significant in the joint sample. In-silico analysis revealed that 8 of the 15 genes predicted to be regulated by both mir-206 and mir-198, are transcriptional targets or interaction partners of the JUN, ATF2 and TAF1 connected in a tight network. JUN and two of the miRNA targets (CCND2 and PTPN1) in the network have previously been associated with schizophrenia. CONCLUSIONS/SIGNIFICANCE: We found nominal association between brain-expressed miRNAs and schizophrenia for rs17578796 and rs1700 located in mir-206 and mir-198 respectively. These two miRNAs have a surprising large number (15) of targets in common, eight of which are also connected by the same transcription factors

The Comparative Osteology of the Petrotympanic Complex (Ear Region) of Extant Baleen Whales (Cetacea: Mysticeti)

Anatomical comparisons of the ear region of baleen whales (Mysticeti) are provided through detailed osteological descriptions and high-resolution photographs of the petrotympanic complex (tympanic bulla and petrosal bone) of all extant species of mysticete cetaceans. Salient morphological features are illustrated and identified, including overall shape of the bulla, size of the conical process of the bulla, morphology of the promontorium, and the size and shape of the anterior process of the petrosal. We place our comparative osteological observations into a phylogenetic context in order to initiate an exploration into petrotympanic evolution within Mysticeti.The morphology of the petrotympanic complex is diagnostic for individual species of baleen whale (e.g., sigmoid and conical processes positioned at midline of bulla in Balaenoptera musculus; confluence of fenestra cochleae and perilymphatic foramen in Eschrichtius robustus), and several mysticete clades are united by derived characteristics. Balaenids and neobalaenids share derived features of the bulla, such as a rhomboid shape and a reduced anterior lobe (swelling) in ventral aspect, and eschrichtiids share derived morphologies of the petrosal with balaenopterids, including loss of a medial promontory groove and dorsomedial elongation of the promontorium. Monophyly of Balaenoidea (Balaenidae and Neobalaenidae) and Balaenopteroidea (Balaenopteridae and Eschrichtiidae) was recovered in phylogenetic analyses utilizing data exclusively from the petrotympanic complex.This study fills a major gap in our knowledge of the complex structures of the mysticete petrotympanic complex, which is an important anatomical region for the interpretation of the evolutionary history of mammals. In addition, we introduce a novel body of phylogenetically informative characters from the ear region of mysticetes. Our detailed anatomical descriptions, illustrations, and comparisons provide valuable data for current and future studies on the phylogenetic relationships, evolution, and auditory physiology of mysticetes and other cetaceans throughout Earth's history

From Structure Prediction to Genomic Screens for Novel Non-Coding RNAs

Non-coding RNAs (ncRNAs) are receiving more and more attention not only as an abundant class of genes, but also as regulatory structural elements (some located in mRNAs). A key feature of RNA function is its structure. Computational methods were developed early for folding and prediction of RNA structure with the aim of assisting in functional analysis. With the discovery of more and more ncRNAs, it has become clear that a large fraction of these are highly structured. Interestingly, a large part of the structure is comprised of regular Watson-Crick and GU wobble base pairs. This and the increased amount of available genomes have made it possible to employ structure-based methods for genomic screens. The field has moved from folding prediction of single sequences to computational screens for ncRNAs in genomic sequence using the RNA structure as the main characteristic feature. Whereas early methods focused on energy-directed folding of single sequences, comparative analysis based on structure preserving changes of base pairs has been efficient in improving accuracy, and today this constitutes a key component in genomic screens. Here, we cover the basic principles of RNA folding and touch upon some of the concepts in current methods that have been applied in genomic screens for de novo RNA structures in searches for novel ncRNA genes and regulatory RNA structure on mRNAs. We discuss the strengths and weaknesses of the different strategies and how they can complement each other

The study of atmospheric ice-nucleating particles via microfluidically generated droplets

Ice-nucleating particles (INPs) play a significant role in the climate and hydrological cycle by triggering ice formation in supercooled clouds, thereby causing precipitation and affecting cloud lifetimes and their radiative properties. However, despite their importance, INP often comprise only 1 in 10³–10⁶ ambient particles, making it difficult to ascertain and predict their type, source, and concentration. The typical techniques for quantifying INP concentrations tend to be highly labour-intensive, suffer from poor time resolution, or are limited in sensitivity to low concentrations. Here, we present the application of microfluidic devices to the study of atmospheric INPs via the simple and rapid production of monodisperse droplets and their subsequent freezing on a cold stage. This device offers the potential for the testing of INP concentrations in aqueous samples with high sensitivity and high counting statistics. Various INPs were tested for validation of the platform, including mineral dust and biological species, with results compared to literature values. We also describe a methodology for sampling atmospheric aerosol in a manner that minimises sampling biases and which is compatible with the microfluidic device. We present results for INP concentrations in air sampled during two field campaigns: (1) from a rural location in the UK and (2) during the UK’s annual Bonfire Night festival. These initial results will provide a route for deployment of the microfluidic platform for the study and quantification of INPs in upcoming field campaigns around the globe, while providing a benchmark for future lab-on-a-chip-based INP studies

Rickettsia Phylogenomics: Unwinding the Intricacies of Obligate Intracellular Life

BACKGROUND: Completed genome sequences are rapidly increasing for Rickettsia, obligate intracellular alpha-proteobacteria responsible for various human diseases, including epidemic typhus and Rocky Mountain spotted fever. In light of phylogeny, the establishment of orthologous groups (OGs) of open reading frames (ORFs) will distinguish the core rickettsial genes and other group specific genes (class 1 OGs or C1OGs) from those distributed indiscriminately throughout the rickettsial tree (class 2 OG or C2OGs). METHODOLOGY/PRINCIPAL FINDINGS: We present 1823 representative (no gene duplications) and 259 non-representative (at least one gene duplication) rickettsial OGs. While the highly reductive (approximately 1.2 MB) Rickettsia genomes range in predicted ORFs from 872 to 1512, a core of 752 OGs was identified, depicting the essential Rickettsia genes. Unsurprisingly, this core lacks many metabolic genes, reflecting the dependence on host resources for growth and survival. Additionally, we bolster our recent reclassification of Rickettsia by identifying OGs that define the AG (ancestral group), TG (typhus group), TRG (transitional group), and SFG (spotted fever group) rickettsiae. OGs for insect-associated species, tick-associated species and species that harbor plasmids were also predicted. Through superimposition of all OGs over robust phylogeny estimation, we discern between C1OGs and C2OGs, the latter depicting genes either decaying from the conserved C1OGs or acquired laterally. Finally, scrutiny of non-representative OGs revealed high levels of split genes versus gene duplications, with both phenomena confounding gene orthology assignment. Interestingly, non-representative OGs, as well as OGs comprised of several gene families typically involved in microbial pathogenicity and/or the acquisition of virulence factors, fall predominantly within C2OG distributions. CONCLUSION/SIGNIFICANCE: Collectively, we determined the relative conservation and distribution of 14354 predicted ORFs from 10 rickettsial genomes across robust phylogeny estimation. The data, available at PATRIC (PathoSystems Resource Integration Center), provide novel information for unwinding the intricacies associated with Rickettsia pathogenesis, expanding the range of potential diagnostic, vaccine and therapeutic targets

Discovery of barley miRNAs through deep sequencing of short reads

Background: MicroRNAs are important components of the regulatory network of biological systems and thousands have been discovered in both animals and plants. Systematic investigations performed in species with sequenced genomes such as Arabidopsis, rice, poplar and Brachypodium have provided insights into the evolutionary relationships of this class of small RNAs among plants. However, miRNAs from barley, one of the most important cereal crops, remain unknown. Results: We performed a large scale study of barley miRNAs through deep sequencing of small RNAs extracted from leaves of two barley cultivars. By using the presence of miRNA precursor sequences in related genomes as one of a number of supporting criteria, we identified up to 100 miRNAs in barley. Of these only 56 have orthologs in wheat, rice or Brachypodium that are known to be expressed, while up to 44 appear to be specifically expressed in barley. Conclusions: Our study, the first large scale investigation of small RNAs in barley, has identified up to 100 miRNAs. We demonstrate that reliable identification of miRNAs via deep sequencing in a species whose genome has not been sequenced requires a more careful analysis of sequencing errors than is commonly performed. We devised a read filtering procedure for dealing with errors. In addition, we found that the use of a large dataset of almost 35 million reads permits the use of read abundance distributions along putative precursor sequences as a practical tool for isolating miRNAs in a large background of reads originating from other non-coding and coding RNAs. This study therefore provides a generic approach for discovering novel miRNAs where no genome sequence is available.Andreas W Schreiber, Bu-Jun Shi, Chun-Yuan Huang, Peter Langridge, Ute Bauman