Optimization and Control of Agent-Based Models in Biology: A Perspective

Agent-based models (ABMs) have become an increasingly important mode of inquiry for the life sciences. They are particularly valuable for systems that are not understood well enough to build an equation-based model. These advantages, however, are counterbalanced by the difficulty of analyzing and using ABMs, due to the lack of the type of mathematical tools available for more traditional models, which leaves simulation as the primary approach. As models become large, simulation becomes challenging. This paper proposes a novel approach to two mathematical aspects of ABMs, optimization and control, and it presents a few first steps outlining how one might carry out this approach. Rather than viewing the ABM as a model, it is to be viewed as a surrogate for the actual system. For a given optimization or control problem (which may change over time), the surrogate system is modeled instead, using data from the ABM and a modeling framework for which ready-made mathematical tools exist, such as differential equations, or for which control strategies can explored more easily. Once the optimization problem is solved for the model of the surrogate, it is then lifted to the surrogate and tested. The final step is to lift the optimization solution from the surrogate system to the actual system. This program is illustrated with published work, using two relatively simple ABMs as a demonstration, Sugarscape and a consumer-resource ABM. Specific techniques discussed include dimension reduction and approximation of an ABM by difference equations as well systems of PDEs, related to certain specific control objectives. This demonstration illustrates the very challenging mathematical problems that need to be solved before this approach can be realistically applied to complex and large ABMs, current and future. The paper outlines a research program to address them

Import of cytochrome c into mitochondria

The import of cytochrome c into mitochondria can be resolved into a number of discrete steps. Here we report on the covalent attachment of heme to apocytochrome c by the enzyme cytochrome c heme lyase in mitochondria from Neurospora crassa. A new method was developed to measure directly the linkage of heme to apocytochrome c. This method is independent of conformational changes in the protein accompanying heme attachment. Tryptic peptides of [35S]cysteine-labelled apocytochrome c, and of enzymatically formed holocytochrome c, were resolved by reverse-phase HPLC. The cysteine-containing peptide to which heme was attached eluted later than the corresponding peptide from apocytochrome c and could be quantified by counting 35S radioactivity as a measure of holocytochrome c formation. Using this procedure, the covalent attachment of heme to apocytochrome c, which is dependent on the enzyme cytochrome c heme lyase, could be measured. Activity required heme (as hemin) and could be reversibly inhibited by the analogue deuterohemin. Holocytochrome c formation was stimulated 5–10-fold by NADH > NADPH > glutathione and was independent of a potential across the inner mitochondrial membrane. NADH was not required for the binding of apocytochrome c to mitochondria and was not involved in the reduction of the cysteine thiols prior to heme attachment. Holocytochrome c formation was also dependent on a cytosolic factor that was necessary for the heme attaching step of cytochrome c import. The factor was a heat-stable, protease-insensitive, low-molecular-mass component of unknown function. Cytochrome c heme lyase appeared to be a soluble protein located in the mitochondrial intermembrane space and was distinct from the previously identified apocytochrome c binding protein having a similar location. A model is presented in which the covalent attachment of heme by cytochrome c heme lyase also plays an essential role in the import pathway of cytochrome c

Tumors and tumor-like lesions of the heart valves

Valvular tumors and tumor-like lesions may have similar morphological and clinical characteristics, and may place the patients at a high risk of stroke in different ways. From January 2004 to June 2008, 11 patients underwent surgery for a suspected valvular tumor. Valvular tumor and tumor-like lesions accounted for 0.32% of adult cardiac operations. Five (45.5%) valvular lesions were papillary fibroelastomas, one (9.1%) was myxoma, 2 (18.2%) were organized thrombi, and 3 (27.3%) were calcification lesions. There was a total of 5 (45.5%) atrioventricular valve lesions, 4 arising from the atrial side of the leaflets, and one from the ventricular side. All 5 (45.5%) semilunar valvular lesions were from the aortic valve. One (9.1%) lesion originated from the chorda tendinea of the mitral valve. All leaflet lesions were resected by a simple shave technique, and all the patients recovered favorably. Valvular tumor and tumor-like lesions are rare. Pre-operative differential diagnoses among these valvular lesions pose important clinical implications for appropriate treatment for the underlying diseases. Prompt therapeutic measures in view of the underlying diseases of the valvular lesions are essential to prevent potential embolic events

Computational modelling of NF-κB activation by IL-1RI and its co-receptor TILRR, predicts a role for Cytoskeletal Sequestration of IκBα in inflammatory signalling.

The transcription factor NF-κB (nuclear factor kappa B) is activated by Toll-like receptors and controlled by mechanotransduction and changes in the cytoskeleton. In this study we combine 3-D predictive protein modelling and in vitro experiments with in silico simulations to determine the role of the cytoskeleton in regulation of NF-κB. Simulations used a comprehensive agent-based model of the NF-κB pathway, which includes the type 1 IL-1 receptor (IL-1R1) complex and signalling intermediates, as well as cytoskeletal components. Agent based modelling relies on in silico reproductions of systems through the interactions of its components, and provides a reliable tool in investigations of biological processes, which require spatial considerations and involve complex formation and translocation of regulatory components. We show that our model faithfully reproduces the multiple steps comprising the NF-κB pathway, and provides a framework from which we can explore novel aspects of the system. The analysis, using 3-D predictive protein modelling and in vitro assays, demonstrated that the NF-κB inhibitor, IκBα is sequestered to the actin/spectrin complex within the cytoskeleton of the resting cell, and released during IL-1 stimulation, through a process controlled by the IL-1RI co-receptor TILRR (Toll-like and IL-1 receptor regulator). In silico simulations using the agent-based model predict that the cytoskeletal pool of IκBα is released to adjust signal amplification in relation to input levels. The results suggest that the process provides a mechanism for signal calibration and enables efficient, activation-sensitive regulation of NF-κB and inflammatory responses

The Astropy Project: Building an inclusive, open-science project and status of the v2.0 core package

The Astropy project supports and fosters the development of open-source and openly-developed Python packages that provide commonly-needed functionality to the astronomical community. A key element of the Astropy project is the core package Astropy, which serves as the foundation for more specialized projects and packages. In this article, we provide an overview of the organization of the Astropy project and summarize key features in the core package as of the recent major release, version 2.0. We then describe the project infrastructure designed to facilitate and support development for a broader ecosystem of inter-operable packages. We conclude with a future outlook of planned new features and directions for the broader Astropy project