Partitioned likelihood support and the evaluation of data set conflict

In simultaneous analyses of multiple data partitions, the trees relevant when measuring support for a clade are the optimal tree, and the best tree lacking the clade (i.e., the most reasonable alternative). The parsimony-based method of partitioned branch support (PBS) forces each data set to arbitrate between the two relevant trees. This value is the amount each data set contributes to clade support in the combined analysis, and can be very different to support apparent in separate analyses. The approach used in PBS can also be employed in likelihood: a simultaneous analysis of all data retrieves the maximum likelihood tree, and the best tree without the clade of interest is also found. Each data set is fitted to the two trees and the log-likelihood difference calculated, giving partitioned likelihood support (PLS) for each data set. These calculations can be performed regardless of the complexity of the ML model adopted. The significance of PLS can be evaluated using a variety of resampling methods, such as the Kishino-Hasegawa test, the Shimodiara-Hasegawa test, or likelihood weights, although the appropriateness and assumptions of these tests remains debated.

A new principle for tight junction modulation based on occludin peptides

The aim of this study was to investigate whether peptides from the extracellular loops of the tight junction protein occludin could be used as a new principle for tight junction modulation. Peptides of 4 to 47 amino acids in length and covering the two extracellular loops of the tight junction protein occludin were synthesized, and their effect on the tight junction permeability in Caco-2 cells was investigated using [C-14] mannitol as a paracellular marker. Lipopeptide derivatives of one of the active occludin peptides (OPs), synthesized by adding a lipoamino acid containing 14 carbon atoms (C-14-) to the N terminus of the peptide, were also investigated. Peptides corresponding to the N terminus of the first extracellular loop of occludin increased the permeability of the tight junctions without causing short-term toxicity. However, the peptides had an effect only when added to the basolateral side of the cells, which could be partly explained by degradation by apical peptidases and aggregate formation. By contrast, the lipopeptide C-14-OP90-103, which protects the peptide from degradation and aggregation, displayed a rapid apical effect. The L- and D-diastereomers of C-14-OP90-103 had distinctly different effects. The D-isomer, which releases intact OP90-103 from the lipoamino acid, displayed a rapid and transient increase in tight junction permeability. The L- isomer, which releases OP90-103 more rapidly, gave a more sustained increase in tight junction permeability. In conclusion, C-14-OP90-103 represents a prototype of a new class of tight junction modulators that act on the extracellular domains of tight junction proteins.

The transform likelihood ratio method for rare event simulation with heavy tails

We present a novel method, called the transform likelihood ratio (TLR) method, for estimation of rare event probabilities with heavy-tailed distributions. Via a simple transformation ( change of variables) technique the TLR method reduces the original rare event probability estimation with heavy tail distributions to an equivalent one with light tail distributions. Once this transformation has been established we estimate the rare event probability via importance sampling, using the classical exponential change of measure or the standard likelihood ratio change of measure. In the latter case the importance sampling distribution is chosen from the same parametric family as the transformed distribution. We estimate the optimal parameter vector of the importance sampling distribution using the cross-entropy method. We prove the polynomial complexity of the TLR method for certain heavy-tailed models and demonstrate numerically its high efficiency for various heavy-tailed models previously thought to be intractable. We also show that the TLR method can be viewed as a universal tool in the sense that not only it provides a unified view for heavy-tailed simulation but also can be efficiently used in simulation with light-tailed distributions. We present extensive simulation results which support the efficiency of the TLR method.

Estimating disease likelihood: a case of rubbery figures

In diagnosis and prognosis, we should avoid intuitive “guesstimates” and seek a validated numerical aid

First-principle density-functional calculation of the Raman spectra of BEDT-TTF

We present a first-principles density-functional calculation for the Raman spectra of a neutral BEDT-TTF molecule. Our results are in excellent agreement with experimental results. We show that a planar Structure is not a stable state of a neutral BEDT-TTF molecule. We consider three possible conformations and discuss their relation to disorder in these systems.

Modelling cell lifespan and proliferation: is likelihood to die or to divide independent of age?

In cell lifespan studies the exponential nature of cell survival curves is often interpreted as showing the rate of death is independent of the age of the cells within the population. Here we present an alternative model where cells that die are replaced and the age and lifespan of the population pool is monitored until a, steady state is reached. In our model newly generated individual cells are given a determined lifespan drawn from a number of known distributions including the lognormal, which is frequently found in nature. For lognormal lifespans the analytic steady-state survival curve obtained can be well-fit by a single or double exponential, depending on the mean and standard deviation. Thus, experimental evidence for exponential lifespans of one and/or two populations cannot be taken as definitive evidence for time and age independence of cell survival. A related model for a dividing population in steady state is also developed. We propose that the common adoption of age-independent, constant rates of change in biological modelling may be responsible for significant errors, both of interpretation and of mathematical deduction. We suggest that additional mathematical and experimental methods must be used to resolve the relationship between time and behavioural changes by cells that are predominantly unsynchronized.

Coadaptation: A unifying principle in evolutionary thermal biology

Over the last 50 yr, thermal biology has shifted from a largely physiological science to a more integrated science of behavior, physiology, ecology, and evolution. Today, the mechanisms that underlie responses to environmental temperature are being scrutinized at levels ranging from genes to organisms. From these investigations, a theory of thermal adaptation has emerged that describes the evolution of thermoregulation, thermal sensitivity, and thermal acclimation. We review and integrate current models to form a conceptual model of coadaptation. We argue that major advances will require a quantitative theory of coadaptation that predicts which strategies should evolve in specific thermal environments. Simply combining current models, however, is insufficient to understand the responses of organisms to thermal heterogeneity; a theory of coadaptation must also consider the biotic interactions that influence the net benefits of behavioral and physiological strategies. Such a theory will be challenging to develop because each organism's perception of and response to thermal heterogeneity depends on its size, mobility, and life span. Despite the challenges facing thermal biologists, we have never been more pressed to explain the diversity of strategies that organisms use to cope with thermal heterogeneity and to predict the consequences of thermal change for the diversity of communities.

Complex numbers and symmetries in quantum mechanics, and a nonlinear superposition principle for Wigner functions

Complex numbers appear in the Hilbert space formulation of quantum mechanics, but not in the formulation in phase space. Quantum symmetries are described by complex, unitary or antiunitary operators defining ray representations in Hilbert space, whereas in phase space they are described by real, true representations. Equivalence of the formulations requires that the former representations can be obtained from the latter and vice versa. Examples are given. Equivalence of the two formulations also requires that complex superpositions of state vectors can be described in the phase space formulation, and it is shown that this leads to a nonlinear superposition principle for orthogonal, pure-state Wigner functions. It is concluded that the use of complex numbers in quantum mechanics can be regarded as a computational device to simplify calculations, as in all other applications of mathematics to physical phenomena.

First-principle study of adsorption of hydrogen on Ti-doped Mg(0001) surface

Ab initio density functional theory (DFT) calculations are performed to study the adsorption of H-2 molecules on a Ti-doped Mg(0001) surface. We find that two hydrogen molecules are able to dissociate on top of the Ti atom with very small activation barriers (0.103 and 0.145 eV for the first and second H-2 molecules, respectively). Additionally, a molecular adsorption state of H-2 above the Ti atom is observed for the first time and is attributed to the polarization of the H-2 molecule by the Ti cation. Our results parallel recent findings for H-2 adsorption on Ti-doped carbon nanotubes or fullerenes. They provide new insight into the preliminary stages of hydrogen adsorption onto Ti-incorporated Mg surfaces.