Mathematical and computational modeling in biology at multiple scales

A variety of topics are reviewed in the area of mathematical and computational modeling in biology, covering the range of scales from populations of organisms to electrons in atoms. The use of maximum entropy as an inference tool in the fields of biology and drug discovery is discussed. Mathematical and computational methods and models in the areas of epidemiology, cell physiology and cancer are surveyed. The technique of molecular dynamics is covered, with special attention to force fields for protein simulations and methods for the calculation of solvation free energies. The utility of quantum mechanical methods in biophysical and biochemical modeling is explored. The field of computational enzymology is examined.     

Sympatric Breeding Auks Shift between Dietary and Spatial Resource Partitioning across the Annual Cycle

When species competing for the same resources coexist, some segregation in the way they utilize those resources is expected. However, little is known about how closely related sympatric breeding species segregate outside the breeding season. We investigated the annual segregation of three closely related seabirds (razorbill Alca torda , common guillemot Uria aalge and Brünnich’s guillemot U . lomvia ) breeding at the same colony in Southwest Greenland. By combining GPS and geolocation (GLS) tracking with dive depth and stable isotope analyses, we compared spatial and dietary resource partitioning. During the breeding season, we found the three species to segregate in diet and/or dive depth, but less in foraging area. During both the post-breeding and pre-breeding periods, the three species had an increased overlap in diet, but were dispersed over a larger spatial scale. Dive depths were similar across the annual cycle, suggesting morphological adaptations fixed by evolution. Prey choice, on the other hand, seemed much more flexible and therefore more likely to be affected by the immediate presence of potential competitors.     

Brainstem structures responsible for paradoxical sleep onset and maintenance

This paper is dedicated to our mentor, Michel Jouvet who inspired our career and transmitted to us his passion for the study of the mechanisms responsible for paradoxical sleep genesis and also that of its still mysterious functions. We expose in the following the progresses in the knowledge in this field brought during 40 years by Michel Jouvet and his team and more recently by the members of a new CNRS laboratory in which we aim to pursue in the path opened by Michel Jouvet.     

The paradox of the clumps mathematically explained

The lumpy distribution of species along a continuous one-dimensional niche axis recently found by Scheffer and van Nes (Scheffer and van Ness 2006) is explained mathematically. We show that it emerges simply from the eigenvalue and eigenvectors of the community matrix. Both the transient patterns—lumps and gaps between them—as well as the asymptotic equilibrium are explained. If the species are evenly distributed along the niche axis, the emergence of these patterns can be demonstrated analytically. The more general case, of randomly distributed species, shows only slight deviations and is illustrated by numerical simulation. This is a robust result whenever the finiteness of the niche is taken into account: it can be extended to different analytic dependence of the interaction coefficients with the distance on the niche axis (i.e., different kernel interactions), different boundary conditions, etc. We also found that there is a critical value both for the width of the species distribution σ and the number of species n below which the clusterization disappears.

Sequence of a human immunoglobulin gamma 3 heavy chain constant region gene: comparison with the other human C gamma genes.

We report the first and complete nucleotide sequence of a human gamma 3 heavy chain constant region gene (C gamma 3). This gene displays the same organization than the others C gamma genes and exhibits normal RNA splice and polyadenylation sites. A comparison of its primary sequence with those of C gamma 1, C gamma 2 and C gamma 4 genes confirms the high degree of homology (95%) of the human family in both coding and non-coding regions, and the divergence of the hinge region. The C gamma 3 gene we sequenced codes for a Gm(b) gamma 3 chain (EZZ). Comparison with other known protein sequences reveals that only two specific aminoacids are involved in the Gm(b) and Gm(g) allotypes, which suggests an important part of the spatial configuration in the allotypic specificities.     

Comparative tests of ectoparasite species richness in seabirds

Background  

The diversity of parasites attacking a host varies substantially among different host species. Understanding the factors that explain these patterns of parasite diversity is critical to identifying the ecological principles underlying biodiversity. Seabirds (Charadriiformes, Pelecaniformes and Procellariiformes) and their ectoparasitic lice (Insecta: Phthiraptera) are ideal model groups in which to study correlates of parasite species richness. We evaluated the relative importance of morphological (body size, body weight, wingspan, bill length), life-history (longevity, clutch size), ecological (population size, geographical range) and behavioural (diving versus non-diving) variables as predictors of louse diversity on 413 seabird hosts species. Diversity was measured at the level of louse suborder, genus, and species, and uneven sampling of hosts was controlled for using literature citations as a proxy for sampling effort.     

Various rat adult tissues express only one major mRNA species from the glyceraldehyde-3-phosphate-dehydrogenase multigenic family.

We have isolated and sequenced a full-length cDNA clone encoding rat glyceraldehyde-3-phosphate-dehydrogenase (GAPDH, E.C.1.2.1.12). The entire mRNA is 1269 nucleotides long exclusive of poly(A) and contains respectively 71 and 196 bases of 5' and 3' non-coding regions. Primer extension as well as S1 nuclease protection experiments clearly established that a single (or at least a highly prominent) GAPDH mRNA species is expressed in all rat tissues examined. This sequence allowed the determination of the hitherto unknown primary structure of rat GAPDH which is 333 aminoacids long. Comparison between GAPDH sequences from rat, man and chicken revealed a high degree of sequence conservation at both nucleotide and protein levels.     

Improved Pancreatic Adenocarcinoma Diagnosis in Jaundiced and Non-Jaundiced Pancreatic Adenocarcinoma Patients through the Combination of Routine Clinical Markers Associated to Pancreatic Adenocarcinoma Pathophysiology

Background

There is still no reliable biomarker for the diagnosis of pancreatic adenocarcinoma. Carbohydrate antigen 19–9 (CA 19–9) is a tumor marker only recommended for pancreatic adenocarcinoma follow-up. One of the clinical problems lies in distinguishing between this cancer and other benign pancreatic diseases such as chronic pancreatitis. In this study we will assess the value of panels of serum molecules related to pancreatic cancer physiopathology to determine whether alone or in combination could help to discriminate between these two pathologies.

Methods

CA 19–9, carcinoembryonic antigen (CEA), C-reactive protein, albumin, insulin growth factor-1 (IGF-1) and IGF binding protein-3 were measured using routine clinical analyzers in a cohort of 47 pancreatic adenocarcinoma, 20 chronic pancreatitis and 15 healthy controls.

Results

The combination of CA 19–9, IGF-1 and albumin resulted in a combined area under the curve (AUC) of 0.959 with 93.6% sensitivity and 95% specificity, much higher than CA 19–9 alone. An algorithm was defined to classify the patients as chronic pancreatitis or pancreatic cancer with the above specificity and sensitivity. In an independent validation group of 20 pancreatic adenocarcinoma and 13 chronic pancreatitis patients, the combination of the four molecules classified correctly all pancreatic adenocarcinoma and 12 out of 13 chronic pancreatitis patients.

Conclusions

Although this panel of markers should be validated in larger cohorts, the high sensitivity and specificity values and the convenience to measure these parameters in clinical laboratories shows great promise for improving pancreatic adenocarcinoma diagnosis.