Intramolecular vibrational dephasing obeys a power law at intermediate times

Experimental intramolecular vibrational dephasing transients for several large organic molecules are reanalyzed. Fits to the experimental data, as well as full numerical quantum calculations with a factorized potential surface for all active degrees of freedom of fluorene indicate that power law decays, not exponentials, occur at intermediate times. The results support a proposal that power law decays describe vibrational dephasing dynamics in large molecules at intermediate times because of the local nature of energy flow.

A power law for cells

Darwin observed that multiple, lowly organized, rudimentary, or exaggerated structures show increased relative variability. However, the cellular basis for these laws has never been investigated. Some animals, such as the nematode Caenorhabditis elegans, are famous for having organs that possess the same number of cells in all individuals, a property known as eutely. But for most multicellular creatures, the extent of cell number variability is unknown. Here we estimate variability in organ cell number for a variety of animals, plants, slime moulds, and volvocine algae. We find that the mean and variance in cell number obey a power law with an exponent of 2, comparable to Taylor's law in ecological processes. Relative cell number variability, as measured by the coefficient of variation, differs widely across taxa and tissues, but is generally independent of mean cell number among homologous tissues of closely related species. We show that the power law for cell number variability can be explained by stochastic branching process models based on the properties of cell lineages. We also identify taxa in which the precision of developmental control appears to have evolved. We propose that the scale independence of relative cell number variability is maintained by natural selection.

A universal scaling law between gray matter and white matter of cerebral cortex

Neocortex, a new and rapidly evolving brain structure in mammals, has a similar layered architecture in species over a wide range of brain sizes. Larger brains require longer fibers to communicate between distant cortical areas; the volume of the white matter that contains long axons increases disproportionally faster than the volume of the gray matter that contains cell bodies, dendrites, and axons for local information processing, according to a power law. The theoretical analysis presented here shows how this remarkable anatomical regularity might arise naturally as a consequence of the local uniformity of the cortex and the requirement for compact arrangement of long axonal fibers. The predicted power law with an exponent of 4/3 minus a small correction for the thickness of the cortex accurately accounts for empirical data spanning several orders of magnitude in brain sizes for various mammalian species, including human and nonhuman primates.

A simple explanation for taxon abundance patterns

For taxonomic levels higher than species, the abundance distributions of the number of subtaxa per taxon tend to approximate power laws but often show strong deviations from such laws. Previously, these deviations were attributed to finite-time effects in a continuous-time branching process at the generic level. Instead, we describe herein a simple discrete branching process that generates the observed distributions and find that the distribution's deviation from power law form is not caused by disequilibration, but rather that it is time independent and determined by the evolutionary properties of the taxa of interest. Our model predicts—with no free parameters—the rank-frequency distribution of the number of families in fossil marine animal orders obtained from the fossil record. We find that near power law distributions are statistically almost inevitable for taxa higher than species. The branching model also sheds light on species-abundance patterns, as well as on links between evolutionary processes, self-organized criticality, and fractals.

Predicting species diversity in tropical forests

A fundamental question in ecology is how many species occur within a given area. Despite the complexity and diversity of different ecosystems, there exists a surprisingly simple, approximate answer: the number of species is proportional to the size of the area raised to some exponent. The exponent often turns out to be roughly 1/4. This power law can be derived from assumptions about the relative abundances of species or from notions of self-similarity. Here we analyze the largest existing data set of location-mapped species: over one million, individually identified trees from five tropical forests on three continents. Although the power law is a reasonable, zeroth-order approximation of our data, we find consistent deviations from it on all spatial scales. Furthermore, tropical forests are not self-similar at areas ≤50 hectares. We develop an extended model of the species-area relationship, which enables us to predict large-scale species diversity from small-scale data samples more accurately than any other available method.

PartsList: a web-based system for dynamically ranking protein folds based on disparate attributes, including whole-genome expression and interaction information

As the number of protein folds is quite limited, a mode of analysis that will be increasingly common in the future, especially with the advent of structural genomics, is to survey and re-survey the finite parts list of folds from an expanding number of perspectives. We have developed a new resource, called PartsList, that lets one dynamically perform these comparative fold surveys. It is available on the web at http://bioinfo.mbb.yale.edu/partslist and http://www.partslist.org. The system is based on the existing fold classifications and functions as a form of companion annotation for them, providing ‘global views’ of many already completed fold surveys. The central idea in the system is that of comparison through ranking; PartsList will rank the approximately 420 folds based on more than 180 attributes. These include: (i) occurrence in a number of completely sequenced genomes (e.g. it will show the most common folds in the worm versus yeast); (ii) occurrence in the structure databank (e.g. most common folds in the PDB); (iii) both absolute and relative gene expression information (e.g. most changing folds in expression over the cell cycle); (iv) protein–protein interactions, based on experimental data in yeast and comprehensive PDB surveys (e.g. most interacting fold); (v) sensitivity to inserted transposons; (vi) the number of functions associated with the fold (e.g. most multi-functional folds); (vii) amino acid composition (e.g. most Cys-rich folds); (viii) protein motions (e.g. most mobile folds); and (ix) the level of similarity based on a comprehensive set of structural alignments (e.g. most structurally variable folds). The integration of whole-genome expression and protein–protein interaction data with structural information is a particularly novel feature of our system. We provide three ways of visualizing the rankings: a profiler emphasizing the progression of high and low ranks across many pre-selected attributes, a dynamic comparer for custom comparisons and a numerical rankings correlator. These allow one to directly compare very different attributes of a fold (e.g. expression level, genome occurrence and maximum motion) in the uniform numerical format of ranks. This uniform framework, in turn, highlights the way that the frequency of many of the attributes falls off with approximate power-law behavior (i.e. according to V–b, for attribute value V and constant exponent b), with a few folds having large values and most having small values.

Slip complexity in earthquake fault models.

We summarize studies of earthquake fault models that give rise to slip complexities like those in natural earthquakes. For models of smooth faults between elastically deformable continua, it is critical that the friction laws involve a characteristic distance for slip weakening or evolution of surface state. That results in a finite nucleation size, or coherent slip patch size, h*. Models of smooth faults, using numerical cell size properly small compared to h*, show periodic response or complex and apparently chaotic histories of large events but have not been found to show small event complexity like the self-similar (power law) Gutenberg-Richter frequency-size statistics. This conclusion is supported in the present paper by fully inertial elastodynamic modeling of earthquake sequences. In contrast, some models of locally heterogeneous faults with quasi-independent fault segments, represented approximately by simulations with cell size larger than h* so that the model becomes "inherently discrete," do show small event complexity of the Gutenberg-Richter type. Models based on classical friction laws without a weakening length scale or for which the numerical procedure imposes an abrupt strength drop at the onset of slip have h* = 0 and hence always fall into the inherently discrete class. We suggest that the small-event complexity that some such models show will not survive regularization of the constitutive description, by inclusion of an appropriate length scale leading to a finite h*, and a corresponding reduction of numerical grid size.

Kinetics and mechanism of polyamide ("peptide") nucleic acid binding to duplex DNA.

To elucidate the mechanism of recognition of double-stranded DNA (dsDNA) by homopyrimidine polyamide ("peptide") nucleic acid (PNA) leading to the strand-displacement, the kinetics of the sequence-specific PNA/DNA binding have been studied. The binding was monitored with time by the gel retardation and nuclease S1 cleavage assays. The experimental kinetic curves obey pseudo-first-order kinetics and the dependence of the pseudo-first-order rate constant, kps, on PNA concentration, P, obeys a power law kps approximately P gamma with 2 < gamma < 3. The kps values for binding of decamer PNA to dsDNA target sites with one mismatch are hundreds of times slower than for the correct site. A detailed kinetic scheme for PNA/DNA binding is proposed that includes two major steps of the reaction of strand invasion: (i) a transient partial opening of the PNA binding site on dsDNA and incorporation of one PNA molecule with the formation of an intermediate PNA/DNA duplex and (ii) formation of a very stable PNA2/DNA triplex. A simple theoretical treatment of the proposed kinetic scheme is performed. The interpretation of our experimental data in the framework of the proposed kinetic scheme leads to the following conclusions. The sequence specificity of the recognition is essentially provided at the "search" step of the process, which consists in the highly reversible transient formation of duplex between one PNA molecule and the complementary strand of duplex DNA while the other DNA strand is displaced. This search step is followed by virtually irreversible "locking" step via PNA2/DNA triplex formation. The proposed mechanism explains how the binding of homopyrimidine PNA to dsDNA meets two apparently mutually contradictory features: high sequence specificity of binding and remarkable stability of both correct and mismatched PNA/DNA complexes.