Solving the heart mechanics equations with Newton and quasi Newton methods--a comparison

The non-linear elasticity equations of heart mechanics are solved while emulating the effects of a propagating activation wave. The dynamics of a 1 cm(3) slab of active cardiac tissue was simulated as the electrical wave traversed the muscular heart wall transmurally. The regular Newton (Newton-Raphson) method was compared to two modified Newton approaches, and also to a third approach that delayed update only of some selected Jacobian elements. In addition, the impact of changing the time step (0.01, 0.1 and 1 ms) and the relative non-linear convergence tolerance (10(-4), 10(-3) and 10(-2)) was investigated. Updating the Jacobian only when slow convergence occurred was by far the most efficient approach, giving time savings of 83-96%. For each of the four methods, CPU times were reduced by 48-90% when the time step was increased by a factor 10. Increasing the convergence tolerance by the same factor gave time savings of 3-71%. Different combinations of activation wave speed, stress rate and bulk modulus revealed that the fastest method became relatively even faster as stress rate and bulk modulus was decreased, while the activation speed had negligible influence in this respect.     

Unfolding a folding disease: folding, misfolding and aggregation of the marble brain syndrome-associated mutant H107Y of human carbonic anhydrase II

Most loss-of-function diseases are caused by aberrant folding of important proteins. These proteins often misfold due to mutations. The disease marble brain syndrome (MBS), known also as carbonic anhydrase II deficiency syndrome (CADS), can manifest in carriers of point mutations in the human carbonic anhydrase II (HCA II) gene. One mutation associated with MBS entails the His107Tyr substitution. Here, we demonstrate that this mutation is a remarkably destabilizing folding mutation. The loss-of-function is clearly a folding defect, since the mutant shows 64% of CO(2) hydration activity compared to that of the wild-type at low temperature where the mutant is folded. On the contrary, its stability towards thermal and guanidine hydrochloride (GuHCl) denaturation is highly compromised. Using activity assays, CD, fluorescence, NMR, cross-linking, aggregation measurements and molecular modeling, we have mapped the properties of this remarkable mutant. Loss of enzymatic activity had a midpoint temperature of denaturation (T(m)) of 16 degrees C for the mutant compared to 55 degrees C for the wild-type protein. GuHCl-denaturation (at 4 degrees C) showed that the native state of the mutant was destabilized by 9.2kcal/mol. The mutant unfolds through at least two equilibrium intermediates; one novel intermediate that we have termed the molten globule light state and, after further denaturation, the classical molten globule state is populated. Under physiological conditions (neutral pH; 37 degrees C), the His107Tyr mutant will populate the molten globule light state, likely due to novel interactions between Tyr107 and the surroundings of the critical residue Ser29 that destabilize the native conformation. This intermediate binds the hydrophobic dye 8-anilino-1-naphthalene sulfonic acid (ANS) but not as strong as the molten globule state, and near-UV CD reveals the presence of significant tertiary structure. Notably, this intermediate is not as prone to aggregation as the classical molten globule. As a proof of concept for an intervention strategy with small molecules, we showed that binding of the CA inhibitor acetazolamide increases the stability of the native state of the mutant by 2.9kcal/mol in accordance with its strong affinity. Acetazolamide shifts the T(m) to 34 degrees C that protects from misfolding and will enable a substantial fraction of the enzyme pool to survive physiological conditions.     

BESST - Efficient scaffolding of large fragmented assemblies

Background

The use of short reads from High Throughput Sequencing (HTS) techniques is now commonplace in de novo assembly. Yet, obtaining contiguous assemblies from short reads is challenging, thus making scaffolding an important step in the assembly pipeline. Different algorithms have been proposed but many of them use the number of read pairs supporting a linking of two contigs as an indicator of reliability. This reasoning is intuitive, but fails to account for variation in link count due to contig features.We have also noted that published scaffolders are only evaluated on small datasets using output from only one assembler. Two issues arise from this. Firstly, some of the available tools are not well suited for complex genomes. Secondly, these evaluations provide little support for inferring a software’s general performance.

Results

We propose a new algorithm, implemented in a tool called BESST, which can scaffold genomes of all sizes and complexities and was used to scaffold the genome of P. abies (20 Gbp). We performed a comprehensive comparison of BESST against the most popular stand-alone scaffolders on a large variety of datasets. Our results confirm that some of the popular scaffolders are not practical to run on complex datasets. Furthermore, no single stand-alone scaffolder outperforms the others on all datasets. However, BESST fares favorably to the other tested scaffolders on GAGE datasets and, moreover, outperforms the other methods when library insert size distribution is wide.

Conclusion

We conclude from our results that information sources other than the quantity of links, as is commonly used, can provide useful information about genome structure when scaffolding.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2105-15-281) contains supplementary material, which is available to authorized users.     

Real-time viability assay based on 51Cr retention in adherent cells

The chromium (51Cr) release assay has been widely used for viability measurements, even though it has major disadvantages such as high manual workload and poor time resolution. By the use of LigandTracer 51Cr release viability measurements on adherent cells can be significantly simplified and improved. LigandTracer enables a time-resolved detection of 5SCr in target cells, with the result that the effect of toxic material is updated continuously throughout the experiment. Here we explain the principle behind this novel real-time viability assay and show viability curves for known toxic compounds on A431 and U343MGaCl2:6 cell lines.     

Histidine-rich glycoprotein protects from systemic Candida infection

Fungi, such as Candida spp., are commonly found on the skin and at mucosal surfaces. Yet, they rarely cause invasive infections in immunocompetent individuals, an observation reflecting the ability of our innate immune system to control potentially invasive microbes found at biological boundaries. Antimicrobial proteins and peptides are becoming increasingly recognized as important effectors of innate immunity. This is illustrated further by the present investigation, demonstrating a novel antifungal role of histidine-rich glycoprotein (HRG), an abundant and multimodular plasma protein. HRG bound to Candida cells, and induced breaks in the cell walls of the organisms. Correspondingly, HRG preferentially lysed ergosterol-containing liposomes but not cholesterol-containing ones, indicating a specificity for fungal versus other types of eukaryotic membranes. Both antifungal and membrane-rupturing activities of HRG were enhanced at low pH, and mapped to the histidine-rich region of the protein. Ex vivo, HRG-containing plasma as well as fibrin clots exerted antifungal effects. In vivo, Hrg(-/-) mice were susceptible to infection by C. albicans, in contrast to wild-type mice, which were highly resistant to infection. The results demonstrate a key and previously unknown antifungal role of HRG in innate immunity.     

Preliminary results in the immunodiagnosis of tuberculosis in children based on T cell responses to ESAT-6 and PPD antigens

The aim of this work was to study the difference in interferon gamma (IFN-gamma) production by T lymphocytes after early secretory antigen target 6 (ESAT-6) or purified protein derivate (PPD) stimulation in whole blood culture supernatants from children with suspected tuberculosis (TB) disease (n = 21), latent TB infection (n = 16) and negative controls (NC) (n = 22) from an endemic area in Brazil. The concentration of IFN-gamma (pg/ml) was measured by enzyme linked immunosorbent assay and the differences in the IFN-gamma levels for each group were compared and evaluated using an unpaired Student's t-test; p values < 0.05 were considered significant. Measurement of IFN-gamma levels after ESAT-6 stimulation raised the possibility of early diagnosis in the latent TB group (p = 0.0030). Nevertheless, the same group showed similar responses to the NC group (p > 0.05) after PPD stimulation. The IFN-gamma assay using ESAT-6 as an antigenic stimulus has the potential to be used as a tool for the immunodiagnosis of early TB in children.