OpenCMISS: A multi-physics & multi-scale computational infrastructure for the VPH/Physiome project

The VPH/Physiome Project is developing the model encoding standards CellML (cellml.org) and FieldML (fieldml.org) as well as web-accessible model repositories based on these standards (models.physiome.org). Freely available open source computational modelling software is also being developed to solve the partial differential equations described by the models and to visualise results. The OpenCMISS code (opencmiss.org), described here, has been developed by the authors over the last six years to replace the CMISS code that has supported a number of organ system Physiome projects.OpenCMISS is designed to encompass multiple sets of physical equations and to link subcellular and tissue-level biophysical processes into organ-level processes. In the Heart Physiome project, for example, the large deformation mechanics of the myocardial wall need to be coupled to both ventricular flow and embedded coronary flow, and the reaction-diffusion equations that govern the propagation of electrical waves through myocardial tissue need to be coupled with equations that describe the ion channel currents that flow through the cardiac cell membranes.In this paper we discuss the design principles and distributed memory architecture behind the OpenCMISS code. We also discuss the design of the interfaces that link the sets of physical equations across common boundaries (such as fluid-structure coupling), or between spatial fields over the same domain (such as coupled electromechanics), and the concepts behind CellML and FieldML that are embodied in the OpenCMISS data structures. We show how all of these provide a flexible infrastructure for combining models developed across the VPH/Physiome community.     

High resolution positioning of intron ends on the nucleosomes

High resolution sequence-directed nucleosome mapping is applied to 36,000 sequences containing splice junctions, from five different species. As it has been also shown in previous studies, the junctions are found to be preferentially located within nucleosomes. Moreover, the orientation of guanine residues at the GT- and AG-ends of introns within the nucleosomes is such that the guanines are positioned nearest to the surface of histone octamers, 3 and 4 bases upstream from the local DNA pseudo-dyads passing through minor grooves oriented outwards. Since the guanine residues are the most vulnerable to spontaneous damage within the cell (primarily, depurination and oxidation) such positioning of the splice junctions minimizes the damage that is caused by free radicals and highly reactive metabolites.     

The effect of Ginkgo biloba extract on 3-nitropropionic acid-induced neurotoxicity in rats

3-Nitropropionic acid (3-NP), an irreversible inhibitor of succinate dehydrogenase enzyme (SDH), induces neurodegeneration similar to that observed in Huntington’s disease (HD). Reduction of prepulse inhibition (PPI) of acoustic startle response, locomotor hypoactivity, bilateral striatal lesions as well as brain oxidative stress are major features of HD. The present study was designed to investigate neuroprotective effect of Ginkgo biloba extract (EGb 761) on 3-NP induced neurobehavioral changes and striatal lesions.Rats administered 3-NP (20 mg/kg, s.c.) for five consecutive days exhibited PPI deficits and locomotor hypoactivity whereas, pretreatment of animals with EGb 761 (100 mg/kg, i.p. for 15 days) ahead of and during the induction of HD by 3-NP (20 mg/kg for 5 days starting at day 8) ameliorated 3-NP-induced neurobehavioral deficits. Administration of 3-NP increased the level of striatal malondialdehyde (MDA). This effect was prevented in animals pre-treated with EGb 761. Changes in the level of apoptotic regulatory gene expressions, following 3-NP treatment, were demonstrated as both an up-regulation and a down-regulation of the expression levels of striatal Bax and Bcl-xl genes, respectively. In addition, an up-regulation of the expression level of striatal glyceraldehyde-3-phosphate dehydrogenase (GAPDH) was also observed. Pre-treatment with EGb 761 caused a down-regulation in striatal GAPDH and Bax together with an up-regulation of striatal Bcl-xl expression level as compared to the 3-NP treated group. Histochemical examination of striatal tissue showed that EGb 761 significantly prevented 3-NP induced inhibition of SDH activity. Histopathological examination further affirmed the neuroprotective effect of EGb 761 against 3-NP toxicity.Taken together, these results suggest that EGb 761 has a neuroprotective role in the current HD paradigm, which may be related to improvement of energy metabolism, antioxidant properties and antiapoptotic effects.     

基于COMSTAT软件对微生物被膜进行定量分析的方法与意义

目的:微生物被膜是一种具有协调性、功能性和高度结构性的膜状复合物,它可以为微生物提供良好的生存环境,免受外界因素的干扰。研究发现,具有产生微生物被膜能力的细菌治病性明显增强,而现今缺乏一种分析微生物被膜的有效手段。本文探讨利用COMSTAT软件对微生物被膜进行定量分析的方法和意义,为对微生物被膜定性定量分析提供支持手段,从而为研究微生物被膜致病性提供方法理论基础。方法:以葡萄球菌为研究模型,利用激光扫描共聚焦显微镜成像技术,结合COMSTAT微生物被膜分析软件对微生物被膜的单位面积生物量、基质覆盖率、平均厚度、粗糙系数等方面进行定量分析,研究了该葡萄球菌的生物被膜生长变化过程,并考察了抗生素对其生物被膜的抑制作用。结果:在葡萄球菌生物被膜生长过程中,生物量、平均厚度以及平均扩散距离等结构指标数值都有明显增加,而粗糙度和表面积与生物量比值呈现降低趋势,表明了微生物被膜由发生向成熟的转化过程。与此同时,经10μg/mL和100μg/mL的卡那霉素处理得到的葡萄球菌微生物被膜生长受到明显抑制,且随着卡那霉素的浓度增加,抑制效果随之增加。结论:本文运用COMSTAT软件的分析方法首次从生物量、平均厚度等结构指标数值的角度描述了葡萄球菌生物被膜,从而有效评价微生物被膜发生、发展、成熟以及崩解的生长过程。该技术在研究微生物被膜形成的理论机制方面存在潜在价值,可以为研究微生物被膜治病性提供理论基础,具有理论指导意义。     

Scalable Nanohelices for Predictive Studies and Enhanced 3D Visualization

Spring-like materials are ubiquitous in nature and of interest in nanotechnology for energy harvesting, hydrogen storage, and biological sensing applications.  For predictive simulations, it has become increasingly important to be able to model the structure of nanohelices accurately.  To study the effect of local structure on the properties of these complex geometries one must develop realistic models.  To date, software packages are rather limited in creating atomistic helical models.  This work focuses on producing atomistic models of silica glass (SiO₂) nanoribbons and nanosprings for molecular dynamics (MD) simulations. Using an MD model of “bulk” silica glass, two computational procedures to precisely create the shape of nanoribbons and nanosprings are presented.  The first method employs the AWK programming language and open-source software to effectively carve various shapes of silica nanoribbons from the initial bulk model, using desired dimensions and parametric equations to define a helix.  With this method, accurate atomistic silica nanoribbons can be generated for a range of pitch values and dimensions.  The second method involves a more robust code which allows flexibility in modeling nanohelical structures.  This approach utilizes a C++ code particularly written to implement pre-screening methods as well as the mathematical equations for a helix, resulting in greater precision and efficiency when creating nanospring models.  Using these codes, well-defined and scalable nanoribbons and nanosprings suited for atomistic simulations can be effectively created.  An added value in both open-source codes is that they can be adapted to reproduce different helical structures, independent of material.  In addition, a MATLAB graphical user interface (GUI) is used to enhance learning through visualization and interaction for a general user with the atomistic helical structures.  One application of these methods is the recent study of nanohelices via MD simulations for mechanical energy harvesting purposes.     

引物间同源性和裂口对PCR扩增的影响

采用ＰＣ／Ｇｅｎｅ软件对能扩增出特异睡不能扩增的ＰＣＲ引物进行同源性比较,同源间“裂口（ｇａｐｓ）”差数的统计,得出源率小于或等于３５％,“裂口”差数大于或等于３,是获得理想ＰＣＲ扩增效果的关键参数之一的结论,为引物设计提供了一条直观的依据。     

Landscape Habitat Suitability Index Software

ABSTRACT Habitat suitability index (HSI) models are traditionally used to evaluate habitat quality for wildlife at a local scale. Rarely have such models incorporated spatial relationships of habitat components. We introduce Landscape HSImodels, a new Microsoft Windows® (Microsoft, Redmond, WA)—based program that incorporates local habitat as well as landscape-scale attributes to evaluate habitats for 21 species of wildlife. Models for additional species can be constructed using the generic model option. At a landscape scale, attributes include edge effects, patch area, distance to resources, and habitat composition. A moving window approach is used to evaluate habitat composition and interspersion within areas typical of home ranges and territories or larger. The software and sample data are available free of charge from the United States Forest Service, Northern Research Station at http:www.nrs.fs.fed.ushsi .     

SAS‐6 coiled‐coil structure and interaction with SAS‐5 suggest a regulatory mechanism in C. elegans centriole assembly

The centriole is a conserved microtubule‐based organelle essential for both centrosome formation and cilium biogenesis. Five conserved proteins for centriole duplication have been identified. Two of them, SAS‐5 and SAS‐6, physically interact with each other and are codependent for their targeting to procentrioles. However, it remains unclear how these two proteins interact at the molecular level. Here, we demonstrate that the short SAS‐5 C‐terminal domain (residues 390–404) specifically binds to a narrow central region (residues 275–288) of the SAS‐6 coiled coil. This was supported by the crystal structure of the SAS‐6 coiled‐coil domain (CCD), which, together with mutagenesis studies, indicated that the association is mediated by synergistic hydrophobic and electrostatic interactions. The crystal structure also shows a periodic charge pattern along the SAS‐6 CCD, which gives rise to an anti‐parallel tetramer. Overall, our findings establish the molecular basis of the specific interaction between SAS‐5 and SAS‐6, and suggest that both proteins individually adopt an oligomeric conformation that is disrupted upon the formation of the hetero‐complex to facilitate the correct assembly of the nine‐fold symmetric centriole.     

Perl and Python codes to estimate Hardy–Weinberg‐ and linkage‐disequilibrium‐related parameters

We have developed two pieces of Perl and Python codes to calculate the natural log of lambda, a quantity that expresses deviations from Hardy–Weinberg equilibrium, and R, a weighted correlation coefficient between alleles of locus pairs as a quantity that expresses linkage disequilibrium. These quantities can be used in an estimation approach as tools to facilitate the characterization and the comparison of populations.