Numerical methods for non-LTE line radiative transfer: Performance and convergence characteristics

Comparison is made between a number of independent computer programs for radiative transfer in molecular rotational lines. The test models are spherically symmetric circumstellar envelopes with a given density and temperature profile. The first two test models have a simple power law density distribution, constant temperature and a fictive 2-level molecule, while the other two test models consist of an inside-out collapsing envelope observed in rotational transitions of HCO+. For the 2-level molecule test problems all codes agree well to within 0.2%, comparable to the accuracy of the individual codes, for low optical depth and up to 2% for high optical depths (tau=4800). The problem of the collapsing cloud in HCO+ has a larger spread in results, ranging up to 12% for the J=4 population. The spread is largest at the radius where the transition from collisional to radiative excitation occurs. The resulting line profiles for the HCO+ J=4-3 transition agree to within 10%, i.e., within the calibration accuracy of most current telescopes. The comparison project and the results described in this paper provide a benchmark for future code development, and give an indication of the typical accuracy of present day calculations of molecular line transfer.Comment: Accepted for publication in A&

Classifying organisms and artefacts by their outline shapes

We often wish to classify objects by their shapes. Indeed, the study of shapes is an important part of many scientific fields, such as evolutionary biology, structural biology, image processing and archaeology. However, mathematical shape spaces are rather complicated and nonlinear. The most widely used methods of shape analysis, geometric morphometrics, treat the shapes as sets of points. Diffeomorphic methods consider the underlying curve rather than points, but have rarely been applied to real-world problems. Using a machine classifier, we tested the ability of several of these methods to describe and classify the shapes of a variety of organic and man-made objects. We find that one method, based on square-root velocity functions (SRVFs), outperforms all others, including a standard geometric morphometric method (eigenshapes), and that it is also superior to human experts using shape alone. When the SRVF approach is constrained to take account of homologous landmarks it can accurately classify objects of very different shapes. The SRVF method identifies a shortest path between shapes, and we show that this can be used to estimate the shapes of intermediate steps in evolutionary series. Diffeomorphic shape analysis methods, we conclude, now provide practical and effective solutions to many shape description and classification problems in the natural and human sciences.</p

Giant Merkel cell carcinoma of the eyelid: a case report and review of the literature

Merkel cell carcinoma (MCC) is a rare cutaneous tumor and cases located in the eyelid have been described, but still its rarity may lead to difficulty in diagnosis and delay in treatment. A 51-year-old female patient that presented with large lesions in the eyelid underwent surgery after the diagnosis of acute chalazion. Following respiratory distress secondary to pulmonary metastasis, the patient's condition deteriorated and was not fit for complete excision treatment. Histopathological investigation of the biopsies, taken from the tumor, revealed that it was undifferentiated small cell carcinoma. Our aim with this paper is to point out that more cases should be reported for more effective diagnosis, histopathological study, clinical investigation, treatment and prognosis of this specific neoplasm

siRNA-Like Double-Stranded RNAs Are Specifically Protected Against Degradation in Human Cell Extract

RNA interference (RNAi) is a set of intracellular pathways in eukaryotes that controls both exogenous and endogenous gene expression. The power of RNAi to knock down (silence) any gene of interest by the introduction of synthetic small-interfering (si)RNAs has afforded powerful insight into biological function through reverse genetic approaches and has borne a new field of gene therapeutics. A number of questions are outstanding concerning the potency of siRNAs, necessitating an understanding of how short double-stranded RNAs are processed by the cell. Recent work suggests unmodified siRNAs are protected in the intracellular environment, although the mechanism of protection still remains unclear. We have developed a set of doubly-fluorophore labeled RNAs (more precisely, RNA/DNA chimeras) to probe in real-time the stability of siRNAs and related molecules by fluorescence resonance energy transfer (FRET). We find that these RNA probes are substrates for relevant cellular degradative processes, including the RNase H1 mediated degradation of an DNA/RNA hybrid and Dicer-mediated cleavage of a 24-nucleotide (per strand) double-stranded RNA. In addition, we find that 21- and 24-nucleotide double-stranded RNAs are relatively protected in human cytosolic cell extract, but less so in blood serum, whereas an 18-nucleotide double-stranded RNA is less protected in both fluids. These results suggest that RNAi effector RNAs are specifically protected in the cellular environment and may provide an explanation for recent results showing that unmodified siRNAs in cells persist intact for extended periods of time

Enumerating Pathways of Proton Abstraction Based on a Spatial and Electrostatic Analysis of Residues in the Catalytic Site

The pathways of proton abstraction (PA), a key aspect of most catalytic reactions, is often controversial and highly debated. Ultrahigh-resolution diffraction studies, molecular dynamics, quantum mechanics and molecular mechanic simulations are often adopted to gain insights in the PA mechanisms in enzymes. These methods require expertise and effort to setup and can be computationally intensive. We present a push button methodology – Proton abstraction Simulation (PRISM) – to enumerate the possible pathways of PA in a protein with known 3D structure based on the spatial and electrostatic properties of residues in the proximity of a given nucleophilic residue. Proton movements are evaluated in the vicinity of this nucleophilic residue based on distances, potential differences, spatial channels and characteristics of the individual residues (polarity, acidic, basic, etc). Modulating these parameters eliminates their empirical nature and also might reveal pathways that originate from conformational changes. We have validated our method using serine proteases and concurred with the dichotomy in PA in Class A β-lactamases, both of which are hydrolases. The PA mechanism in a transferase has also been corroborated. The source code is made available at www.sanchak.com/prism

Acidithiobacillus ferrooxidans metabolism: from genome sequence to industrial applications

<p>Abstract</p> <p>Background</p> <p><it>Acidithiobacillus ferrooxidans </it>is a major participant in consortia of microorganisms used for the industrial recovery of copper (bioleaching or biomining). It is a chemolithoautrophic, γ-proteobacterium using energy from the oxidation of iron- and sulfur-containing minerals for growth. It thrives at extremely low pH (pH 1–2) and fixes both carbon and nitrogen from the atmosphere. It solubilizes copper and other metals from rocks and plays an important role in nutrient and metal biogeochemical cycling in acid environments. The lack of a well-developed system for genetic manipulation has prevented thorough exploration of its physiology. Also, confusion has been caused by prior metabolic models constructed based upon the examination of multiple, and sometimes distantly related, strains of the microorganism.</p> <p>Results</p> <p>The genome of the type strain <it>A. ferrooxidans </it>ATCC 23270 was sequenced and annotated to identify general features and provide a framework for <it>in silico </it>metabolic reconstruction. Earlier models of iron and sulfur oxidation, biofilm formation, quorum sensing, inorganic ion uptake, and amino acid metabolism are confirmed and extended. Initial models are presented for central carbon metabolism, anaerobic metabolism (including sulfur reduction, hydrogen metabolism and nitrogen fixation), stress responses, DNA repair, and metal and toxic compound fluxes.</p> <p>Conclusion</p> <p>Bioinformatics analysis provides a valuable platform for gene discovery and functional prediction that helps explain the activity of <it>A. ferrooxidans </it>in industrial bioleaching and its role as a primary producer in acidic environments. An analysis of the genome of the type strain provides a coherent view of its gene content and metabolic potential.</p

Assessing the functional coherence of modules found in multiple-evidence networks from Arabidopsis

<p>Abstract</p> <p>Background</p> <p>Combining multiple evidence-types from different information sources has the potential to reveal new relationships in biological systems. The integrated information can be represented as a relationship network, and clustering the network can suggest possible functional modules. The value of such modules for gaining insight into the underlying biological processes depends on their functional coherence. The challenges that we wish to address are to define and quantify the functional coherence of modules in relationship networks, so that they can be used to infer function of as yet unannotated proteins, to discover previously unknown roles of proteins in diseases as well as for better understanding of the regulation and interrelationship between different elements of complex biological systems.</p> <p>Results</p> <p>We have defined the functional coherence of modules with respect to the Gene Ontology (GO) by considering two complementary aspects: (i) the fragmentation of the GO functional categories into the different modules and (ii) the most representative functions of the modules. We have proposed a set of metrics to evaluate these two aspects and demonstrated their utility in <it>Arabidopsis thaliana</it>. We selected 2355 proteins for which experimentally established protein-protein interaction (PPI) data were available. From these we have constructed five relationship networks, four based on single types of data: PPI, co-expression, co-occurrence of protein names in scientific literature abstracts and sequence similarity and a fifth one combining these four evidence types. The ability of these networks to suggest biologically meaningful grouping of proteins was explored by applying Markov clustering and then by measuring the functional coherence of the clusters.</p> <p>Conclusions</p> <p>Relationship networks integrating multiple evidence-types are biologically informative and allow more proteins to be assigned to a putative functional module. Using additional evidence types concentrates the functional annotations in a smaller number of modules without unduly compromising their consistency. These results indicate that integration of more data sources improves the ability to uncover functional association between proteins, both by allowing more proteins to be linked and producing a network where modular structure more closely reflects the hierarchy in the gene ontology.</p

Stratification of biological therapies by pathobiology in biologic-naive patients with rheumatoid arthritis (STRAP and STRAP-EU): two parallel, open-label, biopsy-driven, randomised trials

BACKGROUND: Despite highly effective targeted therapies for rheumatoid arthritis, about 40% of patients respond poorly, and predictive biomarkers for treatment choices are lacking. We did a biopsy-driven trial to compare the response to rituximab, etanercept, and tocilizumab in biologic-naive patients with rheumatoid arthritis stratified for synovial B cell status. METHODS: STRAP and STRAP-EU were two parallel, open-label, biopsy-driven, stratified, randomised, phase 3 trials done across 26 university centres in the UK and Europe. Biologic-naive patients aged 18 years or older with rheumatoid arthritis based on American College of Rheumatology (ACR)–European League Against Rheumatism classification criteria and an inadequate response to conventional synthetic disease-modifying antirheumatic drugs (DMARDs) were included. Following ultrasound-guided synovial biopsy, patients were classified as B cell poor or B cell rich according to synovial B cell signatures and randomly assigned (1:1:1) to intravenous rituximab (1000 mg at week 0 and week 2), subcutaneous tocilizumab (162 mg per week), or subcutaneous etanercept (50 mg per week). The primary outcome was the 16-week ACR20 response in the B cell-poor, intention-to-treat population (defined as all randomly assigned patients), with data pooled from the two trials, comparing etanercept and tocilizumab (grouped) versus rituximab. Safety was assessed in all patients who received at least one dose of study drug. These trials are registered with the EU Clinical Trials Register, 2014-003529-16 (STRAP) and 2017-004079-30 (STRAP-EU). FINDINGS: Between June 8, 2015, and July 4, 2019, 226 patients were randomly assigned to etanercept (n=73), tocilizumab (n=74), and rituximab (n=79). Three patients (one in each group) were excluded after randomisation because they received parenteral steroids in the 4 weeks before recruitment. 168 (75%) of 223 patients in the intention-to-treat population were women and 170 (76%) were White. In the B cell-poor population, ACR20 response at 16 weeks (primary endpoint) showed no significant differences between etanercept and tocilizumab grouped together and rituximab (46 [60%] of 77 patients vs 26 [59%] of 44; odds ratio 1·02 [95% CI 0·47–2·17], p=0·97). No differences were observed for adverse events, including serious adverse events, which occurred in six (6%) of 102 patients in the rituximab group, nine (6%) of 108 patients in the etanercept group, and three (4%) of 73 patients in the tocilizumab group (p=0·53). INTERPRETATION: In this biologic-naive population of patients with rheumatoid arthrtitis, the dichotomic classification into synovial B cell poor versus rich did not predict treatment response to B cell depletion with rituximab compared with alternative treatment strategies. However, the lack of response to rituximab in patients with a pauci-immune pathotype and the higher risk of structural damage progression in B cell-rich patients treated with rituximab warrant further investigations into the ability of synovial tissue analyses to inform disease pathogenesis and treatment response. FUNDING: UK Medical Research Council and Versus Arthritis