Multiplicative random walk Metropolis-Hastings on the real line

In this article we propose multiplication based random walk Metropolis Hastings (MH) algorithm on the real line. We call it the random dive MH (RDMH) algorithm. This algorithm, even if simple to apply, was not studied earlier in Markov chain Monte Carlo literature. The associated kernel is shown to have standard properties like irreducibility, aperiodicity and Harris recurrence under some mild assumptions. These ensure basic convergence (ergodicity) of the kernel. Further the kernel is shown to be geometric ergodic for a large class of target densities on $\mathbb{R}$. This class even contains realistic target densities for which random walk or Langevin MH are not geometrically ergodic. Three simulation studies are given to demonstrate the mixing property and superiority of RDMH to standard MH algorithms on real line. A share-price return data is also analyzed and the results are compared with those available in the literature

Four small puzzles that Rosetta doesn't solve

A complete macromolecule modeling package must be able to solve the simplest structure prediction problems. Despite recent successes in high resolution structure modeling and design, the Rosetta software suite fares poorly on deceptively small protein and RNA puzzles, some as small as four residues. To illustrate these problems, this manuscript presents extensive Rosetta results for four well-defined test cases: the 20-residue mini-protein Trp cage, an even smaller disulfide-stabilized conotoxin, the reactive loop of a serine protease inhibitor, and a UUCG RNA tetraloop. In contrast to previous Rosetta studies, several lines of evidence indicate that conformational sampling is not the major bottleneck in modeling these small systems. Instead, approximations and omissions in the Rosetta all-atom energy function currently preclude discriminating experimentally observed conformations from de novo models at atomic resolution. These molecular "puzzles" should serve as useful model systems for developers wishing to make foundational improvements to this powerful modeling suite.Comment: Published in PLoS One as a manuscript for the RosettaCon 2010 Special Collectio

Work-related musculoskeletal disorders : A survey of physical therapists in Izmir-Turkey

BACKGROUND: This study was planned to collect data about causes, prevalence and responses to work-related musculoskeletal disorders reported by physiotherapists employed in Izmir, Turkey. METHOD: A two-page survey with closed ended questions was used as the data collected method. This survey was distributed to 205 physiotherapists working in Izmir, Turkey, and 120 physiotherapists answered. Questions included occupational history of physiotherapists and musculoskeletal symptoms, special areas, tasks, job-related risk factors, injury prevention strategies, and responses to injury. RESULTS: Eighty-five percent of the physiotherapists have had a musculoskeletal injury once or more in their lifetime. Injuries have been occurred mostly in low back (26 %), hand-wrist (18 %), shoulders (14 %) and neck (12 %). The highest risk factor in causing the injury was transferring the patient at 15%. Sixty-nine percent of physiotherapists visited a physician for their injury and sixty-seven percent of the respondents indicated that they had not limited their patient contact time as a result to their injury CONCLUSIONS: According to the results of this study, the rate of musculoskeletal disorders in physiotherapists in Izmir-Turkey has been found to be high due to their profession. Respondents felt that a change in work habits was required in order to decrease the risk of another injury

Integrating transposable elements in the 3D genome

Chromosome organisation is increasingly recognised as an essential component of genome regulation, cell fate and cell health. Within the realm of transposable elements (TEs) however, the spatial information of how genomes are folded is still only rarely integrated in experimental studies or accounted for in modelling. Whilst polymer physics is recognised as an important tool to understand the mechanisms of genome folding, in this commentary we discuss its potential applicability to aspects of TE biology. Based on recent works on the relationship between genome organisation and TE integration, we argue that existing polymer models may be extended to create a predictive framework for the study of TE integration patterns. We suggest that these models may offer orthogonal and generic insights into the integration profiles (or "topography") of TEs across organisms. In addition, we provide simple polymer physics arguments and preliminary molecular dynamics simulations of TEs inserting into heterogeneously flexible polymers. By considering this simple model, we show how polymer folding and local flexibility may generically affect TE integration patterns. The preliminary discussion reported in this commentary is aimed to lay the foundations for a large-scale analysis of TE integration dynamics and topography as a function of the three-dimensional host genome

Real-Time Imaging and Quantification of Amyloid-β Peptide Aggregates by Novel Quantum-Dot Nanoprobes

Background: Protein aggregation plays a major role in the pathogenesis of neurodegenerative disorders, such as Alzheimer’s disease. However, direct real-time imaging of protein aggregation, including oligomerization and fibrillization, has never been achieved. Here we demonstrate the preparation of fluorescent semiconductor nanocrystal (quantum dot; QD)-labeled amyloid-b peptide (QDAb) and its advanced applications. Methodology/Principal Findings: The QDAb construct retained Ab oligomer-forming ability, and the sizes of these oligomers could be estimated from the relative fluorescence intensities of the imaged spots. Both QDAb coaggregation with intact Ab42 and insertion into fibrils were detected by fluorescence microscopy. The coaggregation process was observed by real-time 3D imaging using slit-scanning confocal microscopy, which showed a typical sigmoid curve with 1.5 h in the lag-time and 12 h until saturation. Inhibition of coaggregation using an anti-Ab antibody can be observed as 3D images on a microscopic scale. Microglia ingested monomeric QDAb more significantly than oligomeric QDAb, and the ingested QDAb was mainly accumulated in the lysosome. Conclusions/Significance: These data demonstrate that QDAb is a novel nanoprobe for studying Ab oligomerization an

Atomic-accuracy prediction of protein loop structures through an RNA-inspired ansatz

Consistently predicting biopolymer structure at atomic resolution from sequence alone remains a difficult problem, even for small sub-segments of large proteins. Such loop prediction challenges, which arise frequently in comparative modeling and protein design, can become intractable as loop lengths exceed 10 residues and if surrounding side-chain conformations are erased. This article introduces a modeling strategy based on a 'stepwise ansatz', recently developed for RNA modeling, which posits that any realistic all-atom molecular conformation can be built up by residue-by-residue stepwise enumeration. When harnessed to a dynamic-programming-like recursion in the Rosetta framework, the resulting stepwise assembly (SWA) protocol enables enumerative sampling of a 12 residue loop at a significant but achievable cost of thousands of CPU-hours. In a previously established benchmark, SWA recovers crystallographic conformations with sub-Angstrom accuracy for 19 of 20 loops, compared to 14 of 20 by KIC modeling with a comparable expenditure of computational power. Furthermore, SWA gives high accuracy results on an additional set of 15 loops highlighted in the biological literature for their irregularity or unusual length. Successes include cis-Pro touch turns, loops that pass through tunnels of other side-chains, and loops of lengths up to 24 residues. Remaining problem cases are traced to inaccuracies in the Rosetta all-atom energy function. In five additional blind tests, SWA achieves sub-Angstrom accuracy models, including the first such success in a protein/RNA binding interface, the YbxF/kink-turn interaction in the fourth RNA-puzzle competition. These results establish all-atom enumeration as a systematic approach to protein structure that can leverage high performance computing and physically realistic energy functions to more consistently achieve atomic resolution.Comment: Identity of four-loop blind test protein and parts of figures 5 have been omitted in this preprint to ensure confidentiality of the protein structure prior to its public releas

The Role of the BMP Signaling Antagonist Noggin in the Development of Prostate Cancer Osteolytic Bone Metastasis

Members of the BMP and Wnt protein families play a relevant role in physiologic and pathologic bone turnover. Extracellular antagonists are crucial for the modulation of their activity. Lack of expression of the BMP antagonist noggin by osteoinductive, carcinoma-derived cell lines is a determinant of the osteoblast response induced by their bone metastases. In contrast, osteolytic, carcinoma-derived cell lines express noggin constitutively. We hypothesized that cancer cell-derived noggin may contribute to the pathogenesis of osteolytic bone metastasis of solid cancers by repressing bone formation. Intra-osseous xenografts of PC-3 prostate cancer cells induced osteolytic lesions characterized not only by enhanced osteoclast-mediated bone resorption, but also by decreased osteoblast-mediated bone formation. Therefore, in this model, uncoupling of the bone remodeling process contributes to osteolysis. Bone formation was preserved in the osteolytic lesions induced by noggin-silenced PC-3 cells, suggesting that cancer cell-derived noggin interferes with physiologic bone coupling. Furthermore, intra-osseous tumor growth of noggin-silenced PC-3 cells was limited, most probably as a result of the persisting osteoblast activity. This investigation provides new evidence for a model of osteolytic bone metastasis where constitutive secretion of noggin by cancer cells mediates inhibition of bone formation, thereby preventing repair of osteolytic lesions generated by an excess of osteoclast-mediated bone resorption. Therefore, noggin suppression may be a novel strategy for the treatment of osteolytic bone metastases

DNaseI Hypersensitivity and Ultraconservation Reveal Novel, Interdependent Long-Range Enhancers at the Complex Pax6 Cis-Regulatory Region

The PAX6 gene plays a crucial role in development of the eye, brain, olfactory system and endocrine pancreas. Consistent with its pleiotropic role the gene exhibits a complex developmental expression pattern which is subject to strict spatial, temporal and quantitative regulation. Control of expression depends on a large array of cis-elements residing in an extended genomic domain around the coding region of the gene. The minimal essential region required for proper regulation of this complex locus has been defined through analysis of human aniridia-associated breakpoints and YAC transgenic rescue studies of the mouse smalleye mutant. We have carried out a systematic DNase I hypersensitive site (HS) analysis across 200 kb of this critical region of mouse chromosome 2E3 to identify putative regulatory elements. Mapping the identified HSs onto a percent identity plot (PIP) shows many HSs correspond to recognisable genomic features such as evolutionarily conserved sequences, CpG islands and retrotransposon derived repeats. We then focussed on a region previously shown to contain essential long range cis-regulatory information, the Pax6 downstream regulatory region (DRR), allowing comparison of mouse HS data with previous human HS data for this region. Reporter transgenic mice for two of the HS sites, HS5 and HS6, show that they function as tissue specific regulatory elements. In addition we have characterised enhancer activity of an ultra-conserved cis-regulatory region located near Pax6, termed E60. All three cis-elements exhibit multiple spatio-temporal activities in the embryo that overlap between themselves and other elements in the locus. Using a deletion set of YAC reporter transgenic mice we demonstrate functional interdependence of the elements. Finally, we use the HS6 enhancer as a marker for the migration of precerebellar neuro-epithelium cells to the hindbrain precerebellar nuclei along the posterior and anterior extramural streams allowing visualisation of migratory defects in both pathways in Pax6(Sey/Sey) mice

The JCMT Plane Survey: early results from the l = 30 degree field

We present early results from the JCMT Plane Survey (JPS), which has surveyed the northern inner Galactic plane between longitudes l=7 and l=63 degrees in the 850-{\mu}m continuum with SCUBA-2, as part of the James Clerk Maxwell Telescope Legacy Survey programme. Data from the l=30 degree survey region, which contains the massive star-forming regions W43 and G29.96, are analysed after approximately 40% of the observations had been completed. The pixel-to-pixel noise is found to be 19 mJy/beam, after a smooth over the beam area, and the projected equivalent noise levels in the final survey are expected to be around 10 mJy/beam. An initial extraction of compact sources was performed using the FellWalker method resulting in the detection of 1029 sources above a 5-{\sigma} surface-brightness threshold. The completeness limits in these data are estimated to be around 0.2 Jy/beam (peak flux density) and 0.8 Jy (integrated flux density) and are therefore probably already dominated by source confusion in this relatively crowded section of the survey. The flux densities of extracted compact sources are consistent with those of matching detections in the shallower ATLASGAL survey. We analyse the virial and evolutionary state of the detected clumps in the W43 star-forming complex and find that they appear younger than the Galactic-plane average

Rationalization and Design of the Complementarity Determining Region Sequences in an Antibody-Antigen Recognition Interface

Protein-protein interactions are critical determinants in biological systems. Engineered proteins binding to specific areas on protein surfaces could lead to therapeutics or diagnostics for treating diseases in humans. But designing epitope-specific protein-protein interactions with computational atomistic interaction free energy remains a difficult challenge. Here we show that, with the antibody-VEGF (vascular endothelial growth factor) interaction as a model system, the experimentally observed amino acid preferences in the antibody-antigen interface can be rationalized with 3-dimensional distributions of interacting atoms derived from the database of protein structures. Machine learning models established on the rationalization can be generalized to design amino acid preferences in antibody-antigen interfaces, for which the experimental validations are tractable with current high throughput synthetic antibody display technologies. Leave-one-out cross validation on the benchmark system yielded the accuracy, precision, recall (sensitivity) and specificity of the overall binary predictions to be 0.69, 0.45, 0.63, and 0.71 respectively, and the overall Matthews correlation coefficient of the 20 amino acid types in the 24 interface CDR positions was 0.312. The structure-based computational antibody design methodology was further tested with other antibodies binding to VEGF. The results indicate that the methodology could provide alternatives to the current antibody technologies based on animal immune systems in engineering therapeutic and diagnostic antibodies against predetermined antigen epitopes