The GP problem: quantifying gene-to-phenotype relationships

In this paper we refer to the gene-to-phenotype modeling challenge as the GP problem. Integrating information across levels of organization within a genotype-environment system is a major challenge in computational biology. However, resolving the GP problem is a fundamental requirement if we are to understand and predict phenotypes given knowledge of the genome and model dynamic properties of biological systems. Organisms are consequences of this integration, and it is a major property of biological systems that underlies the responses we observe. We discuss the E(NK) model as a framework for investigation of the GP problem and the prediction of system properties at different levels of organization. We apply this quantitative framework to an investigation of the processes involved in genetic improvement of plants for agriculture. In our analysis, N genes determine the genetic variation for a set of traits that are responsible for plant adaptation to E environment-types within a target population of environments. The N genes can interact in epistatic NK gene-networks through the way that they influence plant growth and development processes within a dynamic crop growth model. We use a sorghum crop growth model, available within the APSIM agricultural production systems simulation model, to integrate the gene-environment interactions that occur during growth and development and to predict genotype-to-phenotype relationships for a given E(NK) model. Directional selection is then applied to the population of genotypes, based on their predicted phenotypes, to simulate the dynamic aspects of genetic improvement by a plant-breeding program. The outcomes of the simulated breeding are evaluated across cycles of selection in terms of the changes in allele frequencies for the N genes and the genotypic and phenotypic values of the populations of genotypes.

Couple communication and female problem drinking: A behavioral observation study

Couples with alcohol and relationship problems often report poor communication, yet little is known about the communication of maritally distressed couples in which the woman abuses alcohol (MDWA couples). Compared with maritally distressed couples without alcohol problems (MDNA) and couples with neither problem (NDNA), MDWA couples showed a distinctive pattern of negative communication. Similar to MDNA men, MDWA men spoke negatively to their partners but listened positively to their partners much like NDNA men. MDWA women listened negatively, much as MDNA women did, but spoke positively, like NDNA women did. The interactions of MDWA couples can be characterized as a male-demand-female-withdraw pattern, which is a gender reversal of the female-demand-male-withdraw pattern often observed in MDNA couples.

Simulation of the influence of rate- and state-dependent friction on the macroscopic behavior of complex fault zones with the lattice solid model

In order to understand the earthquake nucleation process, we need to understand the effective frictional behavior of faults with complex geometry and fault gouge zones. One important aspect of this is the interaction between the friction law governing the behavior of the fault on the microscopic level and the resulting macroscopic behavior of the fault zone. Numerical simulations offer a possibility to investigate the behavior of faults on many different scales and thus provide a means to gain insight into fault zone dynamics on scales which are not accessible to laboratory experiments. Numerical experiments have been performed to investigate the influence of the geometric configuration of faults with a rate- and state-dependent friction at the particle contacts on the effective frictional behavior of these faults. The numerical experiments are designed to be similar to laboratory experiments by DIETERICH and KILGORE (1994) in which a slide-hold-slide cycle was performed between two blocks of material and the resulting peak friction was plotted vs. holding time. Simulations with a flat fault without a fault gouge have been performed to verify the implementation. These have shown close agreement with comparable laboratory experiments. The simulations performed with a fault containing fault gouge have demonstrated a strong dependence of the critical slip distance D-c on the roughness of the fault surfaces and are in qualitative agreement with laboratory experiments.

Cadherin-directed actin assembly: E-cadherin physically associates with the Arp2/3 complex to direct actin assembly in nascent adhesive contacts

Cadherin cell adhesion molecules are major determinants of tissue patterning which function in cooperation with the actin cytoskeleton [1-4]. In the context of stable adhesion [1], cadherin/catenin complexes are often envisaged to passively scaffold onto cortical actin filaments. However, cadherins also form dynamic adhesive contacts during wound healing and morphogenesis [2]. Here actin polymerization has been proposed to drive cell surfaces together [5], although F-actin reorganization also occurs as cell contacts mature [6]. The interaction between cadherins and actin is therefore likely to depend on the functional state of adhesion. We sought to analyze the relationship between cadherin homophilic binding and cytoskeletal activity during early cadherin adhesive contacts. Dissecting the specific effect of cadherin ligation alone on actin regulation is difficult in native cell-cell contacts, due to the range of juxtacrine signals that can arise when two cell surfaces adhere [7]. We therefore activated homophilic ligation using a specific functional recombinant protein. We report the first evidence that E-cadherin associates with the Arp2/3 complex actin nucleator and demonstrate that cadherin binding can exert an active, instructive influence on cells to mark sites for actin assembly at the cell surface.

Exploring complex fitness surfaces: Multiple ornamentation and polymorphism in male guppies

The sexual ornamentation used by male guppies to attract females comprises many components, each of which varies considerably among males. Although natural and sexual selection have been shown to contribute to divergence among populations in male sexual ornaments, the role of sexual selection in maintaining polymorphism within populations is less clear. We used both parametric quadratic regression and nonparametric projection pursuit regression techniques to reveal the major axes of non-linear sexual selection on male ornaments. We visualized the fitness surfaces defined by these axes using thin-plate splines to allow a direct comparison of the two methodologies. Identification of the major axes of selection and their visualization was critical in determining the form and strength of nonlinear selection. Both types of analysis revealed fitness surfaces comprising three peaks, suggesting that there is more than one way to make an attractive guppy. Disruptive selection may be an important process underlying the presence of multiple sexual ornaments and may contribute to the maintenance of the high levels of polymorphism in male sexual ornaments found in guppy populations.

Bridging the gap between different statistical approaches: An integrated framework for modelling

This paper proposes a template for modelling complex datasets that integrates traditional statistical modelling approaches with more recent advances in statistics and modelling through an exploratory framework. Our approach builds on the well-known and long standing traditional idea of 'good practice in statistics' by establishing a comprehensive framework for modelling that focuses on exploration, prediction, interpretation and reliability assessment, a relatively new idea that allows individual assessment of predictions. The integrated framework we present comprises two stages. The first involves the use of exploratory methods to help visually understand the data and identify a parsimonious set of explanatory variables. The second encompasses a two step modelling process, where the use of non-parametric methods such as decision trees and generalized additive models are promoted to identify important variables and their modelling relationship with the response before a final predictive model is considered. We focus on fitting the predictive model using parametric, non-parametric and Bayesian approaches. This paper is motivated by a medical problem where interest focuses on developing a risk stratification system for morbidity of 1,710 cardiac patients given a suite of demographic, clinical and preoperative variables. Although the methods we use are applied specifically to this case study, these methods can be applied across any field, irrespective of the type of response.

Further evidence for the reliance of catalysis by rabbit muscle pyruvate kinase upon isomerization of the ternary complex between enzyme and products

Isothermal calorimetry has been used to examine the effect of thermodynamic non-ideality on the kinetics of catalysis by rabbit muscle pyruvate kinase as the result of molecular crowding by inert cosolutes. The investigation, designed to detect substrate-mediated isomerization of pyruvate kinase, has revealed a 15% enhancement of maximal velocity by supplementation of reaction mixtures with 0.1 M proline, glycine or sorbitol. This effect of thermodynamic non-ideality implicates the existence of a substrate-induced conformational change that is governed by a minor volume decrease and a very small isomerization constant; and hence, substantiates earlier inferences that the rate-determining step in pyruvate kinase kinetics is isomerization of the ternary enzyme product complex rather than the release of products. (C) 2003 Elsevier Science B.V. All rights reserved.

Database technologies for l-system simulations in virtual plant applications on bioinformatics

One of the most important advantages of database systems is that the underlying mathematics is rich enough to specify very complex operations with a small number of statements in the database language. This research covers an aspect of biological informatics that is the marriage of information technology and biology, involving the study of real-world phenomena using virtual plants derived from L-systems simulation. L-systems were introduced by Aristid Lindenmayer as a mathematical model of multicellular organisms. Not much consideration has been given to the problem of persistent storage for these simulations. Current procedures for querying data generated by L-systems for scientific experiments, simulations and measurements are also inadequate. To address these problems the research in this paper presents a generic process for data-modeling tools (L-DBM) between L-systems and database systems. This paper shows how L-system productions can be generically and automatically represented in database schemas and how a database can be populated from the L-system strings. This paper further describes the idea of pre-computing recursive structures in the data into derived attributes using compiler generation. A method to allow a correspondence between biologists' terms and compiler-generated terms in a biologist computing environment is supplied. Once the L-DBM gets any specific L-systems productions and its declarations, it can generate the specific schema for both simple correspondence terminology and also complex recursive structure data attributes and relationships.

Quasiequilibrium problem in H-spaces

In the present paper, we study the quasiequilibrium problem and generalized quasiequilibrium problem of generalized quasi-variational inequality in H-spaces by a new method. Some new equilibrium existence theorems are given. Our results are different from corresponding given results or contain some recent results as their special cases. (C) 2003 Elsevier Science Ltd. All rights reserved.