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.     

Guarana (Paullinia cupana var. sorbilis), an anciently consumed stimulant from the Amazon rain forest: the seeded-fruit transcriptome

Guarana (Paullinia cupana var. sorbilis) is a plant native to the central Amazon basin. Roasted seed extracts have been used as medicinal beverages since pre-Colombian times, due to their reputation as stimulants, aphrodisiacs, tonics, as well as protectors of the gastrointestinal tract. Guarana plants are commercially cultivated exclusively in Brazil to supply the national carbonated soft-drink industry and natural product stores around the world. In this report, we describe and discuss the annotation of 15,387 ESTs from guarana seeded-fruits, highlighting sequences from the flavonoid and purine alkaloid pathways, and those related to biotic stress avoidance. This is the largest set of sequences registered for the Sapindaceae family.     

An NMR structural study of nickel-substituted rubredoxin

The Ni(II) and Zn(II) derivatives of Desulfovibrio vulgaris rubredoxin (DvRd) have been studied by NMR spectroscopy to probe the structure at the metal centre. The βCH₂ proton pairs from the cysteines that bind the Ni(II) atom have been identified using 1D nuclear Overhauser enhancement (NOE) difference spectra and sequence specifically assigned via NOE correlations to neighbouring protons and by comparison with the published X-ray crystal structure of a Ni(II) derivative of Clostridium pasteurianum rubredoxin. The solution structures of DvRd(Zn) and DvRd(Ni) have been determined and the paramagnetic form refined using pseudocontact shifts. The determination of the magnetic susceptibility anisotropy tensor allowed the contact and pseudocontact contributions to the observed chemical shifts to be obtained. Analysis of the pseudocontact and contact chemical shifts of the cysteine Hβ protons and backbone protons close to the metal centre allowed conclusions to be drawn as to the geometry and hydrogen-bonding pattern at the metal binding site. The importance of NH–S hydrogen bonds at the metal centre for the delocalization of electron spin density is confirmed for rubredoxins and can be extrapolated to metal centres in Cu proteins: amicyanin, plastocyanin, stellacyanin, azurin and pseudoazurin.     

Predictors of Beta-Blocker Intolerance and Mortality in Patients After Acute Coronary Syndrome

Purpose

To investigate the predictors of intolerance to beta-blockers treatment and the 6-month mortality in hospitalized patients with acute coronary syndrome (ACS).

Methods

This was a single-center, prospective, and longitudinal study including 370 consecutive ACS patients in Killip class I or II. BBs were prescribed according to international guidelines and withdrawn if intolerance occurred. The study was approved by the institutional ethics committee of our university. Statistics: the clinical parameters evaluated at admission, and the related intolerance to BBs and death at 6 months were analyzed using logistic regression (p<0.05)in PATIENTS.

Results

BB intolerance was observed in 84 patients and was associated with no prior use of statins (OR: 2.16, 95%CI: 1.26–3.69, p= 0.005) and Killip class II (OR: 2.5, 95%CI: 1.30-4.75, p=0.004) in the model adjusted for age, sex, blood pressure, and renal function. There was no association with ST-segment alteration or left anterior descending coronary artery plaque. Intolerance to BB was associated with the greatest risk of death (OR: 4.5, 95%CI: 2.15–9.40, p<0.001).

Conclusions

After ACS, intolerance to BBs in the first 48 h of admission was associated to non previous use of statin and Killip class II and had a high risk of death within 6 months.     

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.     

Reproductive activity of Coelomera lanio (Coleoptera: Chrysomelidae)

The Cecropia spp. (Cecropiaceae) trees are attacked by several insects; among them, Coelomera lanio (Dalman) (Coleoptera: Chrysomelidae). The reproductive behavior of C. lanio was studied under laboratory conditions (12 hour photoperiod, 24.1 +/- 0.1 degrees C and mean relative humidity 67.7 +/- 0.6%), in Vi?osa, Minas Gerais, Brazil. The insects were reared in cages and in Petri dishes and fed leaves of Cecropia pachystachya Trec. Reproductive activity began 5.8 +/- 0.2 days after adult emergence and mean copulation time was 2.5 +/- 0.1 min. The female began oviposition only after 25.7 +/- 0.7 days. Each female laid a mean of 4.7 +/- 0.4 times (range 1-9). The mean number of eggs per oviposition and female was 129.2 +/- 2.4 and 587.4 +/- 92.1 respectively, and the time between egg-layings averaged 16.3 +/- 0.8 days.     

Quality determination and the repair of poor quality spots in array experiments

Background  

A common feature of microarray experiments is the occurence of missing gene expression data. These missing values occur for a variety of reasons, in particular, because of the filtering of poor quality spots and the removal of undefined values when a logarithmic transformation is applied to negative background-corrected intensities. The efficiency and power of an analysis performed can be substantially reduced by having an incomplete matrix of gene intensities. Additionally, most statistical methods require a complete intensity matrix. Furthermore, biases may be introduced into analyses through missing information on some genes. Thus methods for appropriately replacing (imputing) missing data and/or weighting poor quality spots are required.     

The role of local versus biogeographical processes in influencing diversity and body‐size variation in mammal assemblages

Our objective was to estimate and analyze the body‐size distribution parameters of terrestrial mammal assemblages at different spatial scales, and to determine whether these parameters are controlled by local ecological processes or by larger‐scale ones. Based on 93 local assemblages, plus the complete mammal assemblage from three continents (Africa, North, and South America), we estimated three key distribution parameters (diversity/size slope, skewness, and modal size) and compared the values to those expected if size distributions are mainly controlled by local interactions. Mammal diversity decreased much faster as body size increased than predicted by fractal niche theory, both at continental and at local scales, with continental distributions showing steeper slopes than the localities within them. South America showed a steeper slope (after controlling for species diversity), compared to Africa and North America, at local and continental scales. We also found that skewness and modal body size can show strikingly different correlations with predictor variables, such as species richness and median size, depending on the use of untransformed versus log‐transformed data, due to changes in the distribution density generated by log‐transformation. The main differences in slope, skewness, and modal size between local and continental scales appear to arise from the same biogeographical process, where small‐sized species increase in diversity much faster (due to higher spatial turnover rates) than large‐sized species. This process, which can operate even in the absence of competitive saturation at local scales, generates continental assemblages with steeper slopes, smaller modal sizes, and higher right skewness (toward small‐sized species) compared to local communities. In addition, historical factors can also affect the size distribution slopes, which are significantly steeper, in South American mammal assemblages (probably due to stronger megafauna extinction events in South America) than those in North America and Africa.