Multi-scale computational modelling in biology and physiology

Recent advances in biotechnology and the availability of ever more powerful computers have led to the formulation of increasingly complex models at all levels of biology. One of the main aims of systems biology is to couple these together to produce integrated models across multiple spatial scales and physical processes. In this review, we formulate a definition of multi-scale in terms of levels of biological organisation and describe the types of model that are found at each level. Key issues that arise in trying to formulate and solve multi-scale and multi-physics models are considered and examples of how these issues have been addressed are given for two of the more mature fields in computational biology: the molecular dynamics of ion channels and cardiac modelling. As even more complex models are developed over the coming few years, it will be necessary to develop new methods to model them (in particular in coupling across the interface between stochastic and deterministic processes) and new techniques will be required to compute their solutions efficiently on massively parallel computers. We outline how we envisage these developments occurring.     

Northern richness and southern poverty: contrasting genetic footprints of glacial refugia in the relictual tree Sciadopitys verticillata (Coniferales: Sciadopityaceae)

Sciadopitys verticillata is amongst the most relictual of all plants, being the last living member of an ancient conifer lineage, the Sciadopityaceae, and is distributed in small and disjunct populations in high rainfall regions of Japan. Although mega‐fossils indicate the persistence of the species within Japan through the Pleistocene glacial–interglacial cycles, how the species withstood the colder and drier climates of the glacials is not well known. The present study utilized phylogeography and palaeodistribution modelling to test whether the species survived within pollen‐based coastal temperate forest glacial refugia or within previously unidentified refugia close to its current range. Sixteen chloroplast haplotypes were found that displayed significant geographical structuring. Unexpectedly, northern populations in central Honshu most distant from coastal refugia had the highest chloroplast diversity and were differentiated from the south, a legacy of glacial populations possibly in inland river valleys close to its current northern range. By contrast, populations near putative coastal refugia in southern Japan, harboured the lower chloroplast diversity and were dominated by a single haplotype. Fragment size polymorphism at a highly variable and homoplasious mononucleotide repeat region in the trnT‐trnL intergenic spacer reinforced the contrasting patterns of diversity observed between northern and southern populations. The divergent histories of northern and southern populations revealed in the present study will inform the management of this globally significant conifer. © 2012 The Linnean Society of London, Biological Journal of the Linnean Society, 2012, 108 , 263–277.     

ACCA phosphopeptide recognition by the BRCT repeats of BRCA1

The tumour suppressor gene BRCA1 encodes a 220 kDa protein that participates in multiple cellular processes. The BRCA1 protein contains a tandem of two BRCT repeats at its carboxy-terminal region. The majority of disease-associated BRCA1 mutations affect this region and provide to the BRCT repeats a central role in the BRCA1 tumour suppressor function. The BRCT repeats have been shown to mediate phospho-dependant protein-protein interactions. They recognize phosphorylated peptides using a recognition groove that spans both BRCT repeats. We previously identified an interaction between the tandem of BRCA1 BRCT repeats and ACCA, which was disrupted by germ line BRCA1 mutations that affect the BRCT repeats. We recently showed that BRCA1 modulates ACCA activity through its phospho-dependent binding to ACCA. To delineate the region of ACCA that is crucial for the regulation of its activity by BRCA1, we searched for potential phosphorylation sites in the ACCA sequence that might be recognized by the BRCA1 BRCT repeats. Using sequence analysis and structure modelling, we proposed the Ser1263 residue as the most favourable candidate among six residues, for recognition by the BRCA1 BRCT repeats. Using experimental approaches, such as GST pull-down assay with Bosc cells, we clearly showed that phosphorylation of only Ser1263 was essential for the interaction of ACCA with the BRCT repeats. We finally demonstrated by immunoprecipitation of ACCA in cells, that the whole BRCA1 protein interacts with ACCA when phosphorylated on Ser1263.     

The response of stream periphyton to Pacific salmon: using a model to understand the role of environmental context

相似文献   

Modelling intravascular delivery from drug-eluting stents with biodurable coating: investigation of anisotropic vascular drug diffusivity and arterial drug distribution

In-stent restenosis occurs in coronary arteries after implantation of drug-eluting stents with non-uniform restenosis thickness distribution in the artery cross section. Knowledge of the spatio-temporal drug uptake in the arterial wall is useful for investigating restenosis growth but may often be very expensive/difficult to acquire experimentally. In this study, local delivery of a hydrophobic drug from a drug-eluting stent implanted in a coronary artery is mathematically modelled to investigate the drug release and spatio-temporal drug distribution in the arterial wall. The model integrates drug diffusion in the coating and drug diffusion with reversible binding in the arterial wall. The model is solved by the finite volume method for both high and low drug loadings relative to its solubility in the stent coating with varied isotropic–anisotropic vascular drug diffusivities. Drug release profiles in the coating are observed to depend not only on the coating drug diffusivity but also on the properties of the surrounding arterial wall. Time dependencies of the spatially averaged free- and bound-drug levels in the arterial wall on the coating and vascular drug diffusivities are discussed. Anisotropic vascular drug diffusivities result in slightly different average drug levels in the arterial wall but with very different spatial distributions. Higher circumferential vascular diffusivity results in more uniform drug loading in the upper layers and is potentially beneficial in reducing in-stent restenosis. An analytical expression is derived which can be used to determine regions in the arterial with higher free-drug concentration than bound-drug concentration.     

作物生长模拟模型研究和应用

作物生长模拟模型研究和应用宇振荣（北京农业大学农业生态环境系,１０００９４）ＳｔｕｄｉｅｓｏｎＣｒｏｐＧｒｏｗｔｈＭｏｄｅｌｌｉｎｇａｎｄＩｔｓＡｐｐｌｉｃａｔｉｏｎ．￥ＹｕＺｈｅｎｒｏｎｇ（ＤｅｐａｒｔｍｅｎｔｏｆＡｇｒｏ－Ｅ－ｃｏｌｏｇｉｃａｌＥ...     

Palaeodistribution modelling does not support disjunct Pleistocene refugia in several Central American plant taxa

Aim We combine evidence from palaeoniche modelling studies of several tree species to estimate the extent of Central American forest during the Last Glacial Maximum (LGM). In particular, we ask whether the distributions of these species are likely to have changed since the LGM, and whether LGM distributions coincide with previously proposed Pleistocene refugia in this area. Location Central American wet and seasonally dry forests. Methods We developed ecological niche models using two simulations of Pleistocene climate and occurrence data for 15 Neotropical plant species. We focused on palaeodistribution models of three ‘focal’ tree species that occur in wet and seasonally dry Central American forests, where recent phylogeographic data suggest Pleistocene differentiation coincident with previously proposed refugia. We added predictions from six wet‐forest and six seasonally dry‐forest obligate plant species to gauge whether Pleistocene range shifts were specific to habitat type. Correlation analyses were performed between projected LGM and present distributions, LGM distributions and previously proposed refugia. We also asked whether modelled palaeodistributions were smaller than their current extents. Results According to our models, the ranges of the study species were not reduced during the LGM, and did not correlate with refugial models, regardless of habitat type. Relative range sizes between present and LGM distributions did not indicate significant range changes since the LGM. However, relative range sizes differed overall between the two palaeoclimate models. Main conclusions Many of the modelled palaeodistributions of study species were not restricted to refugia during the LGM, regardless of forest type. While constrained from higher elevations, most species found suitable habitat at coastal margins and on newly exposed land due to lowered sea levels during the LGM. These results offer no corroboration for Pleistocene climate change as a driver of genetic differentiation in the ‘focal’ species. We offer alternative explanations for genetic differentiation found in plant species in this area.     

The eutrophication of shallow coastal lakes in Southwest England — understanding and recommendations for restoration,based on palaeolimnology,historical records,and the modelling of changing phosphorus loads

Palaeolimnological studies of sediments from Slapton Ley and Loe Pool, two coastal freshwater lakes in Southwest England, show that in the period since 1945, they have been eutrophicated by nutrient inputs from intensification of agriculture, but also from sewage effluent. Two simple models have been used to identify the main sources of catchment outputs, and in the case of Slapton Ley, to evaluate historical changes in land use, and their likely effect on lake trophic status.Restoration strategies may also be evaluated using the same models. They suggest that in order to reduce loads on either lake to within OECD permissible limits, not only will all sewage inputs need to be prevented, and non-phosphate detergents used, but also losses from agricultural land must be reduced. This could take the form of the keeping of fewer cattle (the main source of organic nitrogen and phosphorus in both catchments), or the zoning of the respective catchments so that steep slopes close to riparian zones are not used, as at present, for the grazing of livestock.A better option, however, would appear to be the establishment along most of the rivers draining into these lakes, of buffer strips of woodland at least 15 m wide. According to the models, this measure, along with treatment or diversion of sewage effluent, would reduce phosphorus loads upon the lakes to within acceptable limits.