Emerging infectious diseases and animal social systems

Emerging infectious diseases threaten a wide diversity of animals, and important questions remain concerning disease emergence in socially structured populations. We developed a spatially explicit simulation model to investigate whether—and under what conditions—disease-related mortality can impact rates of pathogen spread in populations of polygynous groups. Specifically, we investigated whether pathogen-mediated dispersal (PMD) can occur when females disperse after the resident male dies from disease, thus carrying infections to new groups. We also examined the effects of incubation period and virulence, host mortality and rates of background dispersal, and we used the model to investigate the spread of the virus responsible for Ebola hemorrhagic fever, which currently is devastating African ape populations. Output was analyzed using regression trees, which enable exploration of hierarchical and non-linear relationships. Analyses revealed that the incidence of disease in single-male (polygynous) groups was significantly greater for those groups containing an average of more than six females, while the total number of infected hosts in the population was most sensitive to the number of females per group. Thus, as expected, PMD occurs in polygynous groups and its effects increase as harem size (the number of females) increases. Simulation output further indicated that population-level effects of Ebola are likely to differ among multi-male–multi-female chimpanzees and polygynous gorillas, with larger overall numbers of chimpanzees infected, but more gorilla groups becoming infected due to increased dispersal when the resident male dies. Collectively, our results highlight the importance of social system on the spread of disease in wild mammals.     

Parameter estimation for the distribution of single cell lag times

In Quantitative Microbial Risk Assessment, it is vital to understand how lag times of individual cells are distributed over a bacterial population. Such identified distributions can be used to predict the time by which, in a growth-supporting environment, a few pathogenic cells can multiply to a poisoning concentration level.We model the lag time of a single cell, inoculated into a new environment, by the delay of the growth function characterizing the generated subpopulation. We introduce an easy-to-implement procedure, based on the method of moments, to estimate the parameters of the distribution of single cell lag times. The advantage of the method is especially apparent for cases where the initial number of cells is small and random, and the culture is detectable only in the exponential growth phase.     

Simulation of Y-chromosomal haplotype data

The non-recombining nature of the Y-chromosome determines the non-independence of alleles between loci. The evolution of short tandem repeat (STR) loci in the Y-chromosome is the result of different factors such as differential mutation rates, mutation modes, gene conversion, selection and demographic processes. The degree of correlation between loci is dependent on the magnitude of these processes. The simulation of data is a routine tool used for testing hypotheses in population and evolutionary studies. The most basic parameters hitherto used in lineage haplotype simulations are the allele frequency distributions and mutation rates, assuming either full independence or linkage between loci. In this study we introduce use of the Spearman correlation coefficient to estimate the degree of dependence between non-recombining loci. Then, both the interdependence between loci and the allele frequency distributions at multi-allelic loci are incorporated in an algorithm for simulating haplotypes. We illustrate the method using published and unpublished Y-chromosome STR data.     

Conserved water-mediated H-bonding dynamics of catalytic Asn 175 in plant thiol protease

The role of invariant water molecules in the activity of plant cysteine protease is ubiquitous in nature. On analysing the 11 different Protein DataBank (PDB) structures of plant thiol proteases, the two invariant water molecules W1 and W2 (W220 and W222 in the template 1PPN structure) were observed to form H-bonds with the O_b atom of Asn 175. Extensive energy minimization and molecular dynamics simulation studies up to 2 ns on all the PDB and solvated structures clearly revealed the involvement of the H-bonding association of the two water molecules in fixing the orientation of the asparagine residue of the catalytic triad. From this study, it is suggested that H-bonding of the water molecule at the W1 invariant site better stabilizes the Asn residue at the active site of the catalytic triad.     

A Molecular Dynamics Investigation of Chain Conformational Changes In Compressed Bilayers

Abstract

The conformations of the chains constituting the hydrophilic component of alkyl monolayers and bilayers are investigated by performing molecular dynamics atomistic simulations on these systems at different temperatures. Several monitoring techniques are used to reveal the chain conformations, including atom pair radial distribution functions, evolutions of the torsional angles over thousands of timesteps, frequency distributions of the torsionl angles and ‘snapshot’ plots of the atomic configurations. These methods consistently testify to the stability of the trans (fully extended) character of the strain-free alkyl chains up to room temperature. The chains retain much of this conformation even when the layers are compressed by the application of pressure, to which the chains respond by ‘folding’ at the ends attaching them to the substrate planes while maintaining directions which are mainly normal to these planes. A non-zero gap between the layers is also maintained. A pressure of about 50 kbar abruptly causes all motion in the chains to cease, resulting in a highly ordered lattice structure.     

Computational Oral Absorption Simulation for Low‐Solubility Compounds

Bile micelles play an important role in oral absorption of low‐solubility compounds. Bile micelles can affect solubility, dissolution rate, and permeability. For the pH–solubility profile in bile micelles, the Henderson–Hasselbalch equation should be modified to take bile‐micelle partition into account. For the dissolution rate, in the Nernst–Brunner equation, the effective diffusion coefficient in bile‐micelle media should be used instead of the monomer diffusion coefficient. The diffusion coefficient of bile micelles is 8‐ to 18‐fold smaller than that of monomer molecules. For permeability, the effective diffusion coefficient in the unstirred water layer adjacent to the epithelial membrane, and the free fraction at the epithelial membrane surface should be taken into account. The importance of these aspects is demonstrated here using several in vivo and clinical oral‐absorption data of low‐solubility model compounds. Using the theoretical equations, the food effect on oral absorption is further discussed.     

Comparing Stochastically Restricted Linear Estimators in a Regression Model

Conditions for superiority of the minimum dispersion estimator over another with respect to the covariance matrix are derived when the vector parameter of a regression model is subject to competing stochastic restrictions. The restrictions may also consist both of a deterministic part and a stochastic part.     

Application of a photosystem II model for analysis of fluorescence induction curves in the 100 ns to 10 s time domain after excitation with a saturating light pulse

A mathematical model of photosystem II (PSII) events was used to analyze chlorophyll fluorescence transients in the time domain from 100 ns to 10 s after excitation with a saturating 10-ns flash, applied as a part of specialized illumination protocol, using preparations of a thermophilic strain of the unicellular green alga, Chlorella pyrenoidosa Chick (using both intact and diuron-treated cells). Analysis of simulation results has proven that particular attention should be given to flash-induced recombination processes, including nonradiative recombination in PSII, while subsequent charge transfer along the electron transport chain of thylakoid membrane can be adequately described by a single reaction of quinone reoxidation. The PSII model was extended by taking inhibition by diuron of the electron transport in the acceptor side of PSII into account, which allowed simulation of fluorescence induction curves observed in the presence of this inhibitor. The model parameters were determined (stromal pH, rate constants of nonradiative recombination, and the initial reduction state of the quinone pool) which provided adequate simulation of experimentally observed ratios of the maximal and initial fluorescence levels (F _m/F ₀).     

Adsorption and immobilisation of human insulin on graphene monoxide,silicon carbide and boron nitride nanosheets investigated by molecular dynamics simulation

The adsorption and immobilisation of human insulin onto the bio-compatible nanosheets including graphene monoxide, silicon carbide and boron nitride nanosheets were studied by molecular dynamics simulation at the temperature of 310 K. After equilibration, heating and 100 ns production molecular dynamic runs, it was found that the insulin was adsorbed and immobilised onto the considered surfaces in a native folded state. The structural parameters, including root-mean-square deviation and fluctuation, surface accessible solvent area, radius of gyration (R_g) and the distance between the centre of the mass of immobilised protein and the surface of the considered nanosheets, were measured, analysed and discussed. The energetics of the studied systems such as the interaction energy between protein and nanosheet was also measured and addressed. The discussions were centred on the structural and energetic parameters of the protein and nanosheets, including charge density, hydrophobicity, hydrophilicity and residue polarity. The results also showed that the active site of C-termini of chain B played an important role in the adsorption process and this could be helpful in the protection of insulin in its smart delivery and release applications.