A global resource for computational chemistry

A modular distributable system has been built for high-throughput computation of molecular structures and properties. It has been used to process 250,000 compounds from the NCI database and to make the results searchable by structures and properties. The IUPAC/NIST InChI specification and algorithm has been used to index the structures and enforce integrity during computation. A number of novel features of the PM5 Hamiltonian were identified as a result of the high-throughput approach. The system and the data can be redistributed and reused and promote the value of computed data as a primary chemical resource.     

A computational model for cell differentiation

In this paper, we present a word set generating mechanism, called cell-differentiation system, inspired by the tissue process formation in multicellular organisms, which might model some properties of evolving communities of living cells at the syntactical level. The tools utilized to model these biological phenomena belong to the formal language theory. In this context chromosomal mutations are defined as operations on strings and the differentiation according to the control of gene expression is represented by some random-context conditions in formal languages.In the presented formal framework we prove that in a simplified form of this formalism, with only one cell-type which is regular, one single cell and no mitosis involved, the problem of establishing whether or not the set of vectors of integers indicating the number of cells in each population, is finite, linear or semilinear, is recursively undecidable. However, one can algorithmically decide whether or not a cell-differentiation system of finite cell-type can produce a specific generation of cells.     

A stochastic model for the economic management of a renewable animal resource

An optimal economic harvesting policy, which maximizes the present value of an animal population, capable of renewing itself, is discussed. It is assumed that, unhindered, the successive population levels, X_n, form a Markov chain, with transitions

$\text{X}{}_{\mathrm{n+1}}\text{=?(X}{}_{\mathrm{n}}\text{) + ?}{}_{\mathrm{n}}\text{?(X}{}_{\mathrm{n}}\text{)}$

, where f is the recruitment function, and {?_n} is an iid sequence of random shocks. When a positive set-up cost is present an optimal policy is of the (S,s) type. The optimal population level is compared with that of an equivalent deterministic model. Bioeconomic conditions, which imply the optimality of conservation, or extinction are investigated.     

A computational model for telomere-dependent cell-replicative aging

Telomere shortening provides a molecular basis for the Hayflick limit. Recent data suggest that telomere shortening also influence mitotic rate. We propose a stochastic growth model of this phenomena, assuming that cell division in each time interval is a random process which probability decreases linearly with telomere shortening. Computer simulations of the proposed stochastic telomere-regulated model provides good approximation of the qualitative growth of cultured human mesenchymal stem cells.     

A hybrid computational model for collective cell durotaxis

Collective cell migration is regulated by a complex set of mechanical interactions and cellular mechanisms. Collective migration emerges from mechanisms occurring at single cell level, involving processes like contraction, polymerization and depolymerization, of cell–cell interactions and of cell–substrate adhesion. Here, we present a computational framework which simulates the dynamics of this emergent behavior conditioned by substrates with stiffness gradients. The computational model reproduces the cell’s ability to move toward the stiffer part of the substrate, process known as durotaxis. It combines the continuous formulation of truss elements and a particle-based approach to simulate the dynamics of cell–matrix adhesions and cell–cell interactions. Using this hybrid approach, researchers can quickly create a quantitative model to understand the regulatory role of different mechanical conditions on the dynamics of collective cell migration. Our model shows that durotaxis occurs due to the ability of cells to deform the substrate more in the part of lower stiffness than in the stiffer part. This effect explains why cell collective movement is more effective than single cell movement in stiffness gradient conditions. In addition, we numerically evaluate how gradient stiffness properties, cell monolayer size and force transmission between cells and extracellular matrix are crucial in regulating durotaxis.     

A diffuse biosemiotic model for cell-to-tissue computational closure

The adoption of diffuse rationality creates a practical bridge between biosemiotics and computation in formulating local-to-global self-consistent criteria for cellular-to-tissue interfacing and for the emergence of life and consciousness. Nature is always complex, the more so at biological membranic inter-scalar interfaces. We present an evolutionary model of the relationship between autonomy and dependence across scales, and describe the implications of its alternating complex-rational-complex nature.     

Truthful online double auction based dynamic resource provisioning for multi-objective trade-offs in IaaS clouds

Auction designs have recently been adopted for static and dynamic resource provisioning in IaaS clouds, such as Microsoft Azure and Amazon EC2. However, the existing mechanisms are mostly restricted to simple auctions, single-objective, offline setting, one-sided interactions either among cloud users or cloud service providers (CSPs), and possible misreports of cloud user’s private information. This paper proposes a more realistic scenario of online auctioning for IaaS clouds, with the unique characteristics of elasticity for time-varying arrival of cloud user requests under the time-based server maintenance in cloud data centers. We propose an online truthful double auction technique for balancing the multi-objective trade-offs between energy, revenue, and performance in IaaS clouds, consisting of a weighted bipartite matching based winning-bid determination algorithm for resource allocation and a Vickrey–Clarke–Groves (VCG) driven algorithm for payment calculation of winning bids. Through rigorous theoretical analysis and extensive trace-driven simulation studies exploiting Google cluster workload traces, we demonstrate that our mechanism significantly improves the performance while promising truthfulness, heterogeneity, economic efficiency, individual rationality, and has a polynomial-time computational complexity.

     

A computational model for populations of dividing cells.

A mathematical model of the cell movements due to cell division is presented. In the model we assume that every cell is a computational object with a given volume, and that the cell pushes the neighbouring cells in order to acquire the space for this volume. The Force that each cell exerts over the other cells is derived from a harmonic arbitrary Potential. The main parameter of the model is the average distance among the cells, that checks if the system is in spatial equilibrium or not. We show that just changing the physical constraints we can model two different systems, a two-dimensional culture on a plate and a three-dimensional early embryo. In both cases the patterns of the cell populations we obtain are similar to the real ones.     

A new computational growth model for sea urchin skeletons

A new computational model has been developed to simulate growth of regular sea urchin skeletons. The model incorporates the processes of plate addition and individual plate growth into a composite model of whole-body (somatic) growth. A simple developmental model based on hypothetical morphogens underlies the assumptions used to define the simulated growth processes. The data model is based on a Delaunay triangulation of plate growth center points, using the dual Voronoi polygons to define plate topologies. A spherical frame of reference is used for growth calculations, with affine deformation of the sphere (based on a Young-Laplace membrane model) to result in an urchin-like three-dimensional form. The model verifies that the patterns of coronal plates in general meet the criteria of Voronoi polygonalization, that a morphogen/threshold inhibition model for plate addition results in the alternating plate addition pattern characteristic of sea urchins, and that application of the Bertalanffy growth model to individual plates results in simulated somatic growth that approximates that seen in living urchins. The model suggests avenues of research that could explain some of the distinctions between modern sea urchins and the much more disparate groups of forms that characterized the Paleozoic Era.     

A general structured model for a sequential hermaphrodite population

This paper introduces and analyzes a model of sequential hermaphroditism in the framework of continuously structured population models with sexual reproduction. The model is general in the sense that the birth, transition (from one sex to the other) and death processes of the population are given by arbitrary functions according to a biological meaningful hypotheses. The system is reduced to a single equation introducing the intrinsic sex-ratio subspace. The steady states are analyzed and illustrated for several cases. In particular, neglecting the competition for resources we have explicitly found a unique non-trivial equilibrium which is unstable.     

A computational model for collagen fibre remodelling in the arterial wall

As the interaction between tissue adaptation and the mechanical condition within tissues is complex, mathematical models are desired to study this interrelation. In this study, a mathematical model is presented to investigate the interplay between collagen architecture and mechanical loading conditions in the arterial wall. It is assumed that the collagen fibres align along preferred directions, situated in between the principal stretch directions. The predicted fibre directions represent symmetrically arranged helices and agree qualitatively with morphometric data from literature. At the luminal side of the arterial wall, the fibres are oriented more circumferentially than at the outer side. The discrete transition of the fibre orientation at the media-adventitia interface can be explained by accounting for the different reference configurations of both layers. The predicted pressure-radius relations resemble experimentally measured sigma-shaped curves. As there is a strong coupling between the collagen architecture and the mechanical loading condition within the tissue, we expect that the presented model for collagen remodelling is useful to gain further insight into the processes involved in vascular adaptation, such as growth and smooth muscle tone adaptation.     

A computational model of ameboid deformation and locomotion

Traditional continuum models of ameboid deformation and locomotion are limited by the computational difficulties intrinsic in free boundary conditions. A new model using the immersed boundary method overcomes these difficulties by representing the cell as a force field immersed in fluid domain. The forces can be derived from a direct mechanical interpretation of such cell components as the cell membrane, the actin cortex, and the transmembrane adhesions between the cytoskeleton and the substratum. The numerical cytoskeleton, modeled as a dynamic network of immersed springs, is able to qualitatively model the passive mechanical behavior of a shear-thinning viscoelastic fluid (Bottino 1997). The same network is used to generate active protrusive and contractile forces. When coordinated with the attachment-detachment cycle of the cell's adhesions to the substratum, these forces produce directed locomotion of the model ameba. With this model it is possible to study the effects of altering the numerical parameters upon the motility of the model cell in a manner suggestive of genetic deletion experiments. In the context of this ameboid cell model and its numerical implementation, simulations involving multicellular interaction, detailed internal signaling, and complex substrate geometries are tractable. Received: 5 January 1998 / Revised version: 23 March 1998 / Accepted: 26 March 1998     

A real-time computational model for estimating kinematics of ankle ligaments

Background: An accurate assessment of ankle ligament kinematics is crucial in understanding the injury mechanisms and can help to improve the treatment of an injured ankle, especially when used in conjunction with robot-assisted therapy. A number of computational models have been developed and validated for assessing the kinematics of ankle ligaments. However, few of them can do real-time assessment to allow for an input into robotic rehabilitation programs. Method: An ankle computational model was proposed and validated to quantify the kinematics of ankle ligaments as the foot moves in real-time. This model consists of three bone segments with three rotational degrees of freedom (DOFs) and 12 ankle ligaments. This model uses inputs for three position variables that can be measured from sensors in many ankle robotic devices that detect postures within the foot–ankle environment and outputs the kinematics of ankle ligaments. Validation of this model in terms of ligament length and strain was conducted by comparing it with published data on cadaver anatomy and magnetic resonance imaging. Results: The model based on ligament lengths and strains is in concurrence with those from the published studies but is sensitive to ligament attachment positions. Conclusions: This ankle computational model has the potential to be used in robot-assisted therapy for real-time assessment of ligament kinematics. The results provide information regarding the quantification of kinematics associated with ankle ligaments related to the disability level and can be used for optimizing the robotic training trajectory.     

A computational model for previtamin D(3) production in skin

Low levels of vitamin D have been implicated in a wide variety of health issues from calcemic diseases to cancer, diabetes and cardiovascular disease. For most humans, the majority of vitamin D(3) is derived from sunlight. How much vitamin D is produced under given exposure conditions is still widely discussed. We present a computational model for the production of (pre-)vitamin D within the skin. It accounts for spectral irradiance, optical properties of the skin and concentration profile of provitamin D. Results are computed for various sets of these parameters yielding the distribution of produced previtamin D in the skin.     

Structural model of the mAb 806-EGFR complex using computational docking followed by computational and experimental mutagenesis

In this work, we combined computational protein-protein docking with computational and experimental mutagenesis to predict the structure of the complex formed by monoclonal antibody 806 (mAb 806) and the epidermal growth factor receptor (EGFR). We docked mAb 806, an antitumor antibody, to its epitope of EGFR residues 287-302. Potential mAb 806-EGFR orientations were generated, and computational mutagenesis was used to filter them according to their agreement with experimental mutagenesis data. Further computational mutagenesis suggested additional mutations, which were tested to arrive at a final structure that was most consistent with experimental mutagenesis data. We propose that this is the EGFR-mAb 806 structure, in which mAb 806 binds to an untethered form of the receptor, consistent with published experimental results. The steric hindrance created by the antibody near the EGFR dimer interface interferes with receptor dimerization, and we postulate this as the structural origin for the antitumor effect of mAb 806.     

Autonomic resource management in virtualized data centers using fuzzy logic-based approaches

Data centers, as resource providers, are expected to deliver on performance guarantees while optimizing resource utilization to reduce cost. Virtualization techniques provide the opportunity of consolidating multiple separately managed containers of virtual resources on underutilized physical servers. A key challenge that comes with virtualization is the simultaneous on-demand provisioning of shared physical resources to virtual containers and the management of their capacities to meet service-quality targets at the least cost. This paper proposes a two-level resource management system to dynamically allocate resources to individual virtual containers. It uses local controllers at the virtual-container level and a global controller at the resource-pool level. An important advantage of this two-level control architecture is that it allows independent controller designs for separately optimizing the performance of applications and the use of resources. Autonomic resource allocation is realized through the interaction of the local and global controllers. A novelty of the local controller designs is their use of fuzzy logic-based approaches to efficiently and robustly deal with the complexity and uncertainties of dynamically changing workloads and resource usage. The global controller determines the resource allocation based on a proposed profit model, with the goal of maximizing the total profit of the data center. Experimental results obtained through a prototype implementation demonstrate that, for the scenarios under consideration, the proposed resource management system can significantly reduce resource consumption while still achieving application performance targets.