Dielectric relaxation in proteins: the computational perspective

In photoexcitation and electron transfer, a new dipole or charge is introduced, and the structure is adjusted. This adjustment represents dielectric relaxation, which is the focus of this review. We concentrate on a few selected topics. We discuss linear response theory, as a unifying framework and a tool to describe non-equilibrium states. We review recent, molecular dynamics simulation studies that illustrate the calculation of dynamic and thermodynamic properties, such as Stokes shifts or reorganization free energies. We then turn to the macroscopic, continuum electrostatic view. We recall the physical definition of a dielectric constant and revisit the decomposition of the free energy into a reorganization and a static term. We review some illustrative continuum studies and discuss some difficulties that can arise with the continuum approach. In conclusion, we consider recent developments that will increase the accuracy and broaden the scope of all these methods.     

The fuzzy polynucleotide space: basic properties

MOTIVATION: Any triplet codon may be regarded as a 12-dimensional fuzzy code. Sufficient information about a particular sequence may not be available in certain situations. The investigator will be confronted with imprecise sequences, yet want to make comparisons of sequences. Fuzzy polynucleotides can be compared by using geometrical interpretation of fuzzy sets as points in a hypercube. RESULTS: We introduce the space of fuzzy polynucleotides and a means of measuring dissimilitudes between them. We establish mathematical principles to measure dissimilarities between fuzzy polynucleotides and present several examples in this metric space. We calculate the frequencies of the nucleotides at the three base sites of a codon in the coding sequences of Escherichia coli K-12 and Mycobacterium tuberculosis H37Rv, and consider them as points in that fuzzy space. We compute the distance between the genomes of E.coli and M.tuberculosis.     

Actin: its cumbersome pilgrimage through cellular compartments

In this article, we follow the history of one of the most abundant, most intensely studied proteins of the eukaryotic cells: actin. We report on hallmarks of its discovery, its structural and functional characterization and localization over time, and point to present days' knowledge on its position as a member of a large family. We focus on the rather puzzling number of diverse functions as proposed for actin as a dual compartment protein. Finally, we venture on some speculations as to its origin.     

DNA hairpins: fuel for autonomous DNA devices

We present a study of the hybridization of complementary DNA hairpin loops, with particular reference to their use as fuel for autonomous DNA devices. The rate of spontaneous hybridization between complementary hairpins can be reduced by increasing the neck length or decreasing the loop length. Hairpins with larger loops rapidly form long-lived kissed complexes. Hairpin loops may be opened by strand displacement using an opening strand that contains the same sequence as half of the neck and a "toehold" complementary to a single-stranded domain adjacent to the neck. We find loop opening via an external toehold to be 10-100 times faster than via an internal toehold. We measure rates of loop opening by opening strands that are at least 1000 times faster than the spontaneous interaction between hairpins. We discuss suitable choices for loop, neck, and toehold length for hairpin loops to be used as fuel for autonomous DNA devices.     

Condor-G: A Computation Management Agent for Multi-Institutional Grids

In recent years, there has been a dramatic increase in the number of available computing and storage resources. Yet few tools exist that allow these resources to be exploited effectively in an aggregated form. We present the Condor-G system, which leverages software from Globus and Condor to enable users to harness multi-domain resources as if they all belong to one personal domain. We describe the structure of Condor-G and how it handles job management, resource selection, security, and fault tolerance. We also present results from application experiments with the Condor-G system. We assert that Condor-G can serve as a general-purpose interface to Grid resources, for use by both end users and higher-level program development tools.     

Evolution of mate-choice imprinting: competing strategies

Mate‐choice imprinting, the determination of mating preferences at an early age based on an individual's observation of adults, plays a role in mate choice in a wide variety of animals. Theoretical work has thus far been focused either on the effects of mate‐choice imprinting on the evolution of the male trait used as a mating cue, or on the evolution of imprinting against a nonimprinting background. We ask the question: if multiple types of imprinting are possible in a species, which is likely to evolve? We develop a haploid population genetic model to compare the evolution of three forms of imprinting: paternal, maternal, and oblique (nonparental adult) imprinting. We find that paternal imprinting is the most likely to evolve, whereas maternal and oblique are nearly equivalent. We identify two factors that determine a strategy's success: its “imprinting set,” the set of individuals imprinted upon, and phenogenotypic disequilibrium, the association between imprinted preferences and mating cues. We assess the predictive power of these factors, and find that the imprinting set is the primary determinant of a strategy's success. We suggest that the imprinting set concept may be generalized to predict the success of additional imprinting strategies, such as mate‐choice copying.     

Microsatellites: consensus and controversy

Microsatellite DNA loci have recently been adopted for many biological applications. Comparative studies across a wide range of species has revealed many details of their mutational properties and evolutionary life cycles. Experience shows that a full understanding of these processes is essential to ensure the effective use of microsatellites as analytical tools. In this article, we review the controversies that have arisen as biologists have taken up this new technology and the emerging consensus that has resulted from their debates. We point to the need for comparative DNA sequencing studies to produce input data for a new generation of theoretical models of microsatellite behaviour. We conclude by presenting our own conceptual model, ‘Snakes and Ladders’, as an aid to theory development.     

The canid genome: behavioral geneticists' best friend?

We review a range of studies on the genetic contribution to behavior in canid species. We begin by identifying factors that make canids a promising model in behavioral genetics and proceed to review research over the last decade that has used canids to identify genetic contributions to behavior. We first review studies that have selectively bred dogs to identify genetic contributions to behavior and then review studies that estimate heritability from populations of non‐laboratory bred dogs. We subsequently review studies that used molecular genetics to identify gene–behavior associations and note associations that have been uncovered. We then note challenges in canid behavioral genetics research that require further consideration. We finish by suggesting alternative phenotyping methods and identify areas in which canids may have as yet unexploited advantages, such as in gene–environment interaction studies where genetic factors are found to moderate the effects of environmental variables.     

Saltatory search: a theoretical analysis

Many animal search in a saltatory fashion: they move forward,pause briefly, and move forward again. Although many optimal-foragingmodels have been developed, most do not address how an animalsearches for food. We view search strategies as "time-distance"functions to allow not only for the possibility of oscillationsin body speed, as implied by saltatory search, but other movementpatterns as well, including cruise search. The key feature ofour models is distinguishing between the body position and thescan position (where the forager is looking). We see the varyingmovement of saltatory search as a consequence of the curvaturein the functions that relate body speed to benefits (Jensen'sinequality)     

Phylogeny of Polygonia, Nymphalis and related butterflies (Lepidoptera: Nymphalidae): a total-evidence analysis

We investigated the phylogeny of butterflies in the tribe Nymphalini sensu Harvey 1991, comprising the genera Vanessa, Cynthia, Bassaris, Aglais, Inachis, Nymphalis, Polygonia, Kaniska, Antanartia, Hypanartia, Symbrenthia, Mynes and Araschnia . Evidence from the mitochondrial gene ndl, the nuclear gene 'wingless' and from morphology/ ecology/behaviour were used separately and combined to analyse relationships. Phylogenies based on the different types of data agreed in many aspects of basic topology. We show that an analysis of only wing pattern characters (based on Nijhout's homology system) results in a topology broadly similar to the one resulting from analysis of the complete matrix. We found support for a monophyletic Nymphalini, where Hypanartia may be the sister clade to all other genera. Mynes, Symbrenthia and Araschnia together seem to form another basal clade. Evidence presented gives only moderate support for a monophyletic Vanessa in the wide sense, including also Cynthia and Bassaris , but strong support for the monophyly of the largely holarctic clade Aglais + Inachis + Nymphalis + Polygonia + Kaniska + Roddia . Within the latter group there is strong support for a clade consisting of Aglais + Inachis and for a second clade which includes Nymphalis, Polygonia (and its sister clade, the monotypic Kaniska) as well as Roddia l-album (= Nymphalis vaualbum ). As a consequence of this topology, Aglais is recognized as a taxon separate from Nymphalis . We present a hypothesis of species relationships within the focal group of genera. We also analyse and discuss the implications of excluding or including ecological data in phylogenetic tree construction, when the tree is to be used for studies in phylogenetic ecology.     

Extremotolerance in fungi: evolution on the edge

Our planet offers many opportunities for life on the edge: high and low temperatures, high salt concentrations, acidic and basic conditions and toxic environments, to name but a few extremes. Recent studies have revealed the diversity of fungi that can occur in stressful environments that are hostile to most eukaryotes. We review these studies here, with the additional purpose of proposing some mechanisms that would allow for the evolutionary adaptation of eukaryotic microbial life under extreme conditions. We focus, in particular, on life in ice and life at high salt concentrations, as there is a surprising similarity between the fungal populations in these two kinds of environments, both of which are characterized by low water activity. We propose steps of evolution of generalist species towards the development of specialists in extreme habitats. We argue that traits present in some fungal groups, such as asexuality, synthesis of melanin-like pigments and a flexible morphology, are preadaptations that facilitate persistence and eventual adaptation to conditions on the ecological edge, as well as biotope switches. These processes are important for understanding the evolution of extremophiles; moreover, they have implications for the emergence of novel fungal pathogens.     

W::Neo: a novel dual-selection marker for high efficiency gene targeting in Drosophila

We have recently developed a so-called genomic engineering approach that allows for directed, efficient and versatile modifications of Drosophila genome by combining the homologous recombination (HR)-based gene targeting with site-specific DNA integration. In genomic engineering and several similar approaches, a "founder" knock-out line must be generated first through HR-based gene targeting, which can still be a potentially time and resource intensive process. To significantly improve the efficiency and success rate of HR-based gene targeting in Drosophila, we have generated a new dual-selection marker termed W::Neo, which is a direct fusion between proteins of eye color marker White (W) and neomycin resistance (Neo). In HR-based gene targeting experiments, mutants carrying W::Neo as the selection marker can be enriched as much as fifty times by taking advantage of the antibiotic selection in Drosophila larvae. We have successfully carried out three independent gene targeting experiments using the W::Neo to generate genomic engineering founder knock-out lines in Drosophila.     

Basic networks: Definition and applications

We define basic networks as the undirected subgraphs with minimal number of units in which the distances (geodesics, minimal path lengths) among a set of selected nodes, which we call seeds, in the original graph are conserved. The additional nodes required to draw the basic network are called connectors. We describe a heuristic strategy to find the basic networks of complex graphs. We also show how the characterization of these networks may help to obtain relevant biological information from highly complex protein-protein interaction data.     

Edge detection in vegetation: Jornada revisited

Abstract. We return to the problem of pattern detection with focus on locating vegetation edges. We assume transect sampling with contiguous quadrats and we postulate a vegetation model in which pattern confounds responses to random, serial, and possibly edge-generating effects. Realizing that strong serial effects will enhance the appearance of vegetation continuity, we probe the residuals for evidence of edges. Two independently derived structures furnish the evidence: a deviations profile and an ordination map. The profile displays compositional changes along the transect as deviations from a theoretical continuum state. The ordination traces the transect as an irregular line in vegetation hyperspace. On this irregular line, the transect segments appear as apices and the compositional dissimilarities of the transect segments as inter-apex distances. Sharp apices (in full dimension) isolated by long inter-apex distances (also in full dimension) indicate possible edges. While the subtending angle to an apex reflects changes in population proportions, the criterion for measuring edge sharpness, the height of the deviations profile expresses quantitative effects, the criterion for measuring edge intensity. The deviations profile does not seem to have direct parallel in past edge detection. The inter-apex distance has a parallel in the Ludwig-Cornelius compositional distance. We describe the details of sampling and analysis, and present numerical results. We also include a short glossary in App. 5.     

Adherence and drug resistance: predictions for therapy outcome

We combine standard pharmacokinetics with an established model of viral replication to predict the outcome of therapy as a function of adherence to the drug regimen. We consider two types of treatment failure: failure to eliminate the wild-type virus, and the emergence of drug-resistant virus. Specifically, we determine the conditions under which resistance dominates as a result of imperfect adherence. We derive this result for both single- and triple-drug therapies, with attention to conditions which favour the emergence of viral strains that are resistant to one or more drugs in a cocktail. Our analysis provides quantitative estimates of the degree of adherence necessary to prevent resistance. We derive results specific to the treatment of human immunodeficiency virus infection, but emphasize that our method is applicable to a range of viral or other infections treated by chemotherapy.     

Cultural complexity and demography: The case of folktales

We investigate the relationship between cultural complexity and population size in a non-technological cultural domain for which we have suitable quantitative records: folktales. We define three levels of complexity for folk narratives: the number of tale types, the number of narrative motifs, and, finally, the number of traits in variants of the same type, for two well-known tales for which we have data from previous studies. We found a positive relationship between number of tale types and population size, a negative relationship for the number of narrative motifs, and no relationship for the number of traits. The absence of a consistent relationship between population size and complexity in folktales provides a novel perspective on the current debates in cultural evolution. We propose that the link between cultural complexity and demography could be domain dependent: in some domains (e.g. technology) this link is important, whereas in others, such as folktales, complex traditions can be easily maintained in small populations as well as large ones, as they may appeal to universal cognitive biases.     

Global cardiac function: mechano-energetico-informatics

This review on the global cardiac function covers cardiac mechanics, energetics, and informatics that I have developed with my collaborators over the last 30 years in Japan and USA. We first established E(max) (end-systolic maximum elastance or pressure/volume ratio) as a new index of ventricular contractility using canine hearts. We then expanded the E(max) concept to PVA (systolic pressure-volume area consisting of external mechanical work and mechanical potential energy) as an innovative measure of total mechanical energy of ventricular contraction and discovered it to be a reliable determinant of ventricular energetics or O(2) consumption (V(O(2))). We have discovered that E(max) shifts the V(O(2))-PVA relation and the E(max) dependency (O(2) cost of E(max)) varies among different pathophysiological hearts. We also searched for the basis of E(max) in crossbridge behavior information contained in an X-ray diffraction of papillary muscle. Recently, we established a new integrative analysis to estimate total Ca(2+) recruited for excitation-contraction coupling in a beating heart using the E(max)-PVA-V(O(2)) information. These global, mechano-energetico-informatic approaches seem to facilitate better understanding of cardiac function, as required in the present post-genomic era when more physiomic knowledge is required not only in cardiac function but also in all other physiologic functions.     

Using the general-purpose reactivity indicator: challenging examples

We elucidate the regioselectivity of nucleophilic attack on substituted benzenesulfonates, quinolines, and pyridines using a general-purpose reactivity indicator (GPRI) for electrophiles. We observe that the GPRI is most accurate when the incoming nucleophile resembles a point charge. We further observe that the GPRI often chooses reactive “dead ends” as the most reactive sites as well as sterically hindered reactive sites. This means that care must be taken to remove sites that are inherently unreactive. Generally, among sites where reactions actually occur, the GPRI identifies the sites in the molecule that lead to the kinetically favored product(s). Furthermore, the GPRI can discern which sites react with hard reagents and which sites react with soft reagents. Because it is currently impossible to use the mathematical framework of conceptual DFT to identify sterically inaccessible sites and reactive dead ends, the GPRI is primarily useful as an interpretative, not a predictive, tool.