PMCC PMCC

Search tips
Search criteria

Advanced
Results 1-25 (1226439)

Clipboard (0)
None

Related Articles

1.  Multiscale image-based modeling and simulation of gas flow and particle transport in the human lungs 
Improved understanding of structure and function relationships in the human lungs in individuals and sub-populations is fundamentally important to the future of pulmonary medicine. Image-based measures of the lungs can provide sensitive indicators of localized features, however to provide a better prediction of lung response to disease, treatment and environment, it is desirable to integrate quantifiable regional features from imaging with associated value-added high-level modeling. With this objective in mind, recent advances in computational fluid dynamics (CFD) of the bronchial airways - from a single bifurcation symmetric model to a multiscale image-based subject-specific lung model - will be reviewed. The interaction of CFD models with local parenchymal tissue expansion - assessed by image registration - allows new understanding of the interplay between environment, hot spots where inhaled aerosols could accumulate, and inflammation. To bridge ventilation function with image-derived central airway structure in CFD, an airway geometrical modeling method that spans from the model ‘entrance’ to the terminal bronchioles will be introduced. Finally, the effects of turbulent flows and CFD turbulence models on aerosol transport and deposition will be discussed.
CFD simulation of airflow and particle transport in the human lung has been pursued by a number of research groups, whose interest has been in studying flow physics and airways resistance, improving drug delivery, or investigating which populations are most susceptible to inhaled pollutants. The three most important factors that need to be considered in airway CFD studies are lung structure, regional lung function, and flow characteristics. Their correct treatment is important because the transport of therapeutic or pollutant particles is dependent on the characteristics of the flow by which they are transported; and the airflow in the lungs is dependent on the geometry of the airways and how ventilation is distributed to the peripheral tissue. The human airway structure spans more than 20 generations, beginning with the extra-thoracic airways (oral or nasal cavity, and through the pharynx and larynx to the trachea), then the conducting airways, the respiratory airways, and to the alveoli. The airways in individuals and sub-populations (by gender, age, ethnicity, and normal vs. diseased states) may exhibit different dimensions, branching patterns and angles, and thickness and rigidity. At the local level, one would like to capture detailed flow characteristics, e.g. local velocity profiles, shear stress, and pressure, for prediction of particle transport in an airway (lung structure) model that is specific to the geometry of an individual, to understand how inter-subject variation in airway geometry (normal or pathological) influences the transport and deposition of particles. In a systems biology – or multiscale modeling – approach, these local flow characteristics can be further integrated with epithelial cell models for the study of mechanotransduction. At the global (organ) level, one would like to match regional ventilation (lung function) that is specific to the individual, thus ensuring that the flow that transports inhaled particles is appropriately distributed throughout the lung model. Computational models that do not account for realistic distribution of ventilation are not capable of predicting realistic particle distribution or targeted drug deposition. Furthermore, the flow in the human lung can be transitional or turbulent in the upper and proximal airways, and becomes laminar in the distal airways. The flows in the laminar, transitional and turbulent regimes have different temporal and spatial scales. Therefore, modeling airway structure and predicting gas flow and particle transport at both local and global levels require image-guided multiscale modeling strategies.
In this article, we will review the aforementioned three key aspects of CFD studies of the human lungs: airway structure (conducting airways), lung function (regional ventilation and boundary conditions), and flow characteristics (modeling of turbulent flow and its effect on particle transport). For modeling airway structure, we will focus on the conducting airways, and review both symmetric vs. asymmetric airway models, idealized vs. CT-based airway models, and multiscale subject-specific airway models. Imposition of physiological subject-specific boundary conditions (BCs) in CFD is essential to match regional ventilation in individuals, which is also critical in studying preferential deposition of inhaled aerosols in sub-populations, e.g. normals vs. asthmatics that may exhibit different ventilation patterns. Subject-specific regional ventilation defines flow distributions and characteristics in airway segments and bifurcations, which subsequently determines the transport and deposition of aerosols in the entire lungs. Turbulence models are needed to capture the transient and turbulent nature of the gas flow in the human lungs. Thus, the advantages and disadvantages of different turbulence models as well as their effects on particle transport will be discussed. The ultimate goal of the development is to identify sensitive structural and functional variables in sub-populations of normal and diseased lungs for potential clinical applications.
doi:10.1002/wsbm.1234
PMCID: PMC3763693  PMID: 23843310
2.  Cell-based multiscale computational modeling of small molecule absorption and retention in the lungs 
Pharmaceutical research  2010;27(3):457-467.
Purpose
For optimizing the local, pulmonary targeting of inhaled medications, it is important to analyze the relationship between the physicochemical properties of small molecules and their absorption, retention and distribution in the various cell types of the airways and alveoli.
Methods
A computational, multiscale, cell-based model was constructed to facilitate analysis of pulmonary drug transport and distribution. The relationship between the physicochemical properties and pharmacokinetic profile of monobasic molecules was explored. Experimental absorption data of compounds with diverse structures were used to validate this model. Simulations were performed to evaluate the effect of active transport and organelle sequestration on the absorption kinetics of compounds.
Results
Relating the physicochemical properties to the pharmacokinetic profiles of small molecules reveals how the absorption half-life and distribution of compounds are expected to vary in different cell types and anatomical regions of the lung. Based on logP, pKa and molecular radius, the absorption rate constants (Ka) calculated with the model were consistent with experimental measurements of pulmonary drug absorption.
Conclusions
The cell-based mechanistic model developed herein is an important step towards the rational design of local, lung-targeted medications, facilitating the design and interpretation of experiments aimed at optimizing drug transport properties in lung.
doi:10.1007/s11095-009-0034-9
PMCID: PMC2907074  PMID: 20099073
3.  A CELL-BASED MOLECULAR TRANSPORT SIMULATOR FOR PHARMACOKINETIC PREDICTION AND CHEMINFORMATIC EXPLORATION 
Molecular pharmaceutics  2006;3(6):704-716.
In the body, cell monolayers serve as permeability barriers, determining transport of drug molecules from one organ or tissue compartment to another. After oral administration, for example, drug transport across the epithelial cell monolayer lining the lumen of the intestine determines the fraction of drug in the gut that is absorbed by the body. By modeling passive transcellular transport properties in the presence of an apical to basolateral concentration gradient, we demonstrate how a computational, cell-based molecular transport simulator can be used to define a physicochemical property space occupied by molecules with desirable permeability and intracellular retention characteristics. Considering extracellular domains of cell surface receptors located on the opposite side of a cell monolayer as a drug’s desired site-of-action, simulation of transcellular transport can be used to define the physicochemical properties of molecules with maximal transcellular permeability but minimal intracellular retention. Arguably, these molecules would possess very desirable features: least likely to exhibit non-specific toxicity, metabolism and side effects associated with high (undesirable) intracellular accumulation; and, most likely to exhibit favorable bioavailability and efficacy associated with maximal rates of transport across cells and minimal intracellular retention, resulting in (desirable) accumulation at the extracellular site-of-action. Calculated permeability predictions showed good correlations with PAMPA, Caco2, and intestinal permeability measurements, without “training” the model and without resorting to statistical regression techniques to “fit” the data. Therefore, cell-based molecular transport simulators could be useful in silico screening tools for chemical genomics and drug discovery.
doi:10.1021/mp060046k
PMCID: PMC2710883  PMID: 17140258
Metoprolol; permeability; chemical space; computer aided drug design; virtual screening; chemical genomics; cellular pharmacokinetics; cheminformatics; drug transport; PAMPA; Biopharmaceutics Classification System
4.  Visualizing chemical structure-subcellular localization relationships using fluorescent small molecules as probes of cellular transport 
Background
To study the chemical determinants of small molecule transport inside cells, it is crucial to visualize relationships between the chemical structure of small molecules and their associated subcellular distribution patterns. For this purpose, we experimented with cells incubated with a synthetic combinatorial library of fluorescent, membrane-permeant small molecule chemical agents. With an automated high content screening instrument, the intracellular distribution patterns of these chemical agents were microscopically captured in image data sets, and analyzed off-line with machine vision and cheminformatics algorithms. Nevertheless, it remained challenging to interpret correlations linking the structure and properties of chemical agents to their subcellular localization patterns in large numbers of cells, captured across large number of images.
Results
To address this challenge, we constructed a Multidimensional Online Virtual Image Display (MOVID) visualization platform using off-the-shelf hardware and software components. For analysis, the image data set acquired from cells incubated with a combinatorial library of fluorescent molecular probes was sorted based on quantitative relationships between the chemical structures, physicochemical properties or predicted subcellular distribution patterns. MOVID enabled visual inspection of the sorted, multidimensional image arrays: Using a multipanel desktop liquid crystal display (LCD) and an avatar as a graphical user interface, the resolution of the images was automatically adjusted to the avatar’s distance, allowing the viewer to rapidly navigate through high resolution image arrays, zooming in and out of the images to inspect and annotate individual cells exhibiting interesting staining patterns. In this manner, MOVID facilitated visualization and interpretation of quantitative structure-localization relationship studies. MOVID also facilitated direct, intuitive exploration of the relationship between the chemical structures of the probes and their microscopic, subcellular staining patterns.
Conclusion
MOVID can provide a practical, graphical user interface and computer-assisted image data visualization platform to facilitate bioimage data mining and cheminformatics analysis of high content, phenotypic screening experiments.
doi:10.1186/1758-2946-5-44
PMCID: PMC3852740  PMID: 24093553
Machine vision; Cheminformatics; Virtual reality; Data mining; Optical probes; Multivariate analysis; Human-computer interaction; Graphical user interface
5.  A Bidirectional Coupling Procedure Applied to Multiscale Respiratory Modeling☆ 
Journal of computational physics  2013;244:10.1016/j.jcp.2012.10.021.
In this study, we present a novel multiscale computational framework for efficiently linking multiple lower-dimensional models describing the distal lung mechanics to imaging-based 3D computational fluid dynamics (CFD) models of the upper pulmonary airways in order to incorporate physiologically appropriate outlet boundary conditions. The framework is an extension of the Modified Newton’s Method with nonlinear Krylov accelerator developed by Carlson and Miller [1, 2, 3]. Our extensions include the retention of subspace information over multiple timesteps, and a special correction at the end of a timestep that allows for corrections to be accepted with verified low residual with as little as a single residual evaluation per timestep on average. In the case of a single residual evaluation per timestep, the method has zero additional computational cost compared to uncoupled or unidirectionally coupled simulations. We expect these enhancements to be generally applicable to other multiscale coupling applications where timestepping occurs. In addition we have developed a “pressure-drop” residual which allows for stable coupling of flows between a 3D incompressible CFD application and another (lower-dimensional) fluid system. We expect this residual to also be useful for coupling non-respiratory incompressible fluid applications, such as multiscale simulations involving blood flow.
The lower-dimensional models that are considered in this study are sets of simple ordinary differential equations (ODEs) representing the compliant mechanics of symmetric human pulmonary airway trees. To validate the method, we compare the predictions of hybrid CFD-ODE models against an ODE-only model of pulmonary airflow in an idealized geometry. Subsequently, we couple multiple sets of ODEs describing the distal lung to an imaging-based human lung geometry. Boundary conditions in these models consist of atmospheric pressure at the mouth and intrapleural pressure applied to the multiple sets of ODEs. In both the simplified geometry and in the imaging-based geometry, the performance of the method was comparable to that of monolithic schemes, in most cases requiring only a single CFD evaluation per time step. Thus, this new accelerator allows us to begin combining pulmonary CFD models with lower-dimensional models of pulmonary mechanics with little computational overhead. Moreover, because the CFD and lower-dimensional models are totally separate, this framework affords great flexibility in terms of the type and breadth of the adopted lower-dimensional model, allowing the biomedical researcher to appropriately focus on model design. Research funded by the National Heart and Blood Institute Award 1RO1HL073598.
doi:10.1016/j.jcp.2012.10.021
PMCID: PMC3856712  PMID: 24347680
computational fluid dynamics; multiscale coupling; pulmonary airflows; Krylov subspace; modified Newton-Raphson
6.  Machine vision assisted analysis of structure-localization relationships in a combinatorial library of prospective bioimaging probes 
With a combinatorial library of bioimaging probes, it is now possible to use machine vision to analyze the contribution of different building blocks of the molecules to their cell-associated visual signals. For athis purpose, cell-permeant, fluorescent styryl molecules were synthesized by condensation of 168 aldehyde with 8 pyridinium/quinolinium building blocks. Images of cells incubated with fluorescent molecules were acquired with a high content screening instrument. Chemical and image feature analysis revealed how variation in one or the other building block of the styryl molecules led to variations in the molecules' visual signals. Across each pair of probes in the library, chemical similarity was significantly associated with spectral and total signal intensity similarity. However, chemical similarity was much less associated with similarity in subcellular probe fluorescence patterns. Quantitative analysis and visual inspection of pairs of images acquired from pairs of styryl isomers confirm that many closely-related probes exhibit different subcellular localization patterns. Therefore, idiosyncratic interactions between styryl molecules and specific cellular components greatly contribute to the subcellular distribution of the styryl probes' fluorescence signal. These results demonstrate how machine vision and cheminformatics can be combined to analyze the targeting properties of bioimaging probes, using large image data sets acquired with automated screening systems.
doi:10.1002/cyto.a.20713
PMCID: PMC2692593  PMID: 19243023
Cheminformatics; machine vision; bioimaging; fluorescence; styryl; high content screening; image cytometry; combinatorial chemistry
7.  Impact of hydrogel nanoparticle size and functionalization on in vivo behavior for lung imaging and therapeutics 
Molecular pharmaceutics  2009;6(6):1891-1902.
Polymer chemistry offers the possibility of synthesizing multifunctional nanoparticles which incorporate moieties that enhance diagnostic and therapeutic targeting of cargo delivery to the lung. However, since rules for predicting particle behavior following modification are not well defined, it is essential that probes for tracking fate in vivo are also included. Accordingly, we designed polyacrylamide-based hydrogel particles of differing sizes, functionalized with a nona-arginine cell-penetrating peptide (Arg9), and labeled with imaging components to assess lung retention and cellular uptake after intratracheal administration. Radiolabeled microparticles (1–5 µm diameter) and nanoparticles (20–40 nm diameter) without and with Arg9 showed diffuse airspace distribution by positron emission tomography imaging. Biodistribution studies revealed that particle clearance and extrapulmonary distribution was, in part, size dependent. Microparticles were rapidly cleared by mucociliary routes but unexpectedly, also through the circulation. In contrast, nanoparticles had prolonged lung retention enhanced by Arg9 and were significantly restricted to the lung. For all particle types, uptake was predominant in alveolar macrophages, and, to a lesser extent, lung epithelial cells. In general, particles did not induce local inflammatory responses, with the exception of microparticles bearing Arg9. Whereas microparticles may be advantageous for short-term applications, nano-sized particles constitute an efficient high-retention and non-inflammatory vehicle for the delivery of diagnostic imaging agents and therapeutics to lung airspaces and alveolar macrophages that can be enhanced by Arg9. Importantly, our results show that minor particle modifications may significantly impact in vivo behavior within the complex environments of the lung, underscoring the need for animal modeling.
doi:10.1021/mp900215p
PMCID: PMC2804872  PMID: 19852512
cell penetrating peptide; lung; macrophage; mice; microparticle; nanoparticle; polyacrylamide; polyarginine; positron emission tomography
8.  A Physiologically-Motivated Compartment-Based Model of the Effect of Inhaled Hypertonic Saline on Mucociliary Clearance and Liquid Transport in Cystic Fibrosis 
PLoS ONE  2014;9(11):e111972.
Background
Cystic Fibrosis (CF) lung disease is characterized by liquid hyperabsorption, airway surface dehydration, and impaired mucociliary clearance (MCC). Herein, we present a compartment-based mathematical model of the airway that extends the resolution of functional imaging data.
Methods
Using functional imaging data to inform our model, we developed a system of mechanism-motivated ordinary differential equations to describe the mucociliary clearance and absorption of aerosolized radiolabeled particle and small molecules probes from human subjects with and without CF. We also utilized a novel imaging metric in vitro to gauge the fraction of airway epithelial cells that have functional ciliary activity.
Results
This model, and its incorporated kinetic rate parameters, captures the MCC and liquid dynamics of the hyperabsorptive state in CF airways and the mitigation of that state by hypertonic saline treatment.
Conclusions
We postulate, based on the model structure and its ability to capture clinical patient data, that patients with CF have regions of airway with diminished MCC function that can be recruited with hypertonic saline treatment. In so doing, this model structure not only makes a case for durable osmotic agents used in lung-region specific treatments, but also may provide a possible clinical endpoint, the fraction of functional ciliated airway.
doi:10.1371/journal.pone.0111972
PMCID: PMC4226497  PMID: 25383714
9.  Phenotype Recognition with Combined Features and Random Subspace Classifier Ensemble 
BMC Bioinformatics  2011;12:128.
Background
Automated, image based high-content screening is a fundamental tool for discovery in biological science. Modern robotic fluorescence microscopes are able to capture thousands of images from massively parallel experiments such as RNA interference (RNAi) or small-molecule screens. As such, efficient computational methods are required for automatic cellular phenotype identification capable of dealing with large image data sets. In this paper we investigated an efficient method for the extraction of quantitative features from images by combining second order statistics, or Haralick features, with curvelet transform. A random subspace based classifier ensemble with multiple layer perceptron (MLP) as the base classifier was then exploited for classification. Haralick features estimate image properties related to second-order statistics based on the grey level co-occurrence matrix (GLCM), which has been extensively used for various image processing applications. The curvelet transform has a more sparse representation of the image than wavelet, thus offering a description with higher time frequency resolution and high degree of directionality and anisotropy, which is particularly appropriate for many images rich with edges and curves. A combined feature description from Haralick feature and curvelet transform can further increase the accuracy of classification by taking their complementary information. We then investigate the applicability of the random subspace (RS) ensemble method for phenotype classification based on microscopy images. A base classifier is trained with a RS sampled subset of the original feature set and the ensemble assigns a class label by majority voting.
Results
Experimental results on the phenotype recognition from three benchmarking image sets including HeLa, CHO and RNAi show the effectiveness of the proposed approach. The combined feature is better than any individual one in the classification accuracy. The ensemble model produces better classification performance compared to the component neural networks trained. For the three images sets HeLa, CHO and RNAi, the Random Subspace Ensembles offers the classification rates 91.20%, 98.86% and 91.03% respectively, which compares sharply with the published result 84%, 93% and 82% from a multi-purpose image classifier WND-CHARM which applied wavelet transforms and other feature extraction methods. We investigated the problem of estimation of ensemble parameters and found that satisfactory performance improvement could be brought by a relative medium dimensionality of feature subsets and small ensemble size.
Conclusions
The characteristics of curvelet transform of being multiscale and multidirectional suit the description of microscopy images very well. It is empirically demonstrated that the curvelet-based feature is clearly preferred to wavelet-based feature for bioimage descriptions. The random subspace ensemble of MLPs is much better than a number of commonly applied multi-class classifiers in the investigated application of phenotype recognition.
doi:10.1186/1471-2105-12-128
PMCID: PMC3098787  PMID: 21529372
10.  Optimized Delivery System Achieves Enhanced Endomyocardial Stem Cell Retention 
Background
Regenerative cell-based therapies are associated with limited myocardial retention of delivered stem cells. The objective of this study is to develop an endocardial delivery system for enhanced cell retention.
Methods and Results
Stem cell retention was simulated in silico using one and three-dimensional models of tissue distortion and compliance associated with delivery. Needle designs, predicted to be optimal, were accordingly engineered using nitinol – a nickel and titanium alloy displaying shape memory and super-elasticity. Biocompatibility was tested with human mesenchymal stem cells. Experimental validation was performed with species-matched cells directly delivered into Langendorff-perfused porcine hearts or administered percutaneously into the endocardium of infarcted pigs. Cell retention was quantified by flow cytometry and real time quantitative polymerase chain reaction methodology. Models, computing optimal distribution of distortion calibrated to favor tissue compliance, predicted that a 75°-curved needle featuring small-to-large graded side holes would ensure the highest cell retention profile. In isolated hearts, the nitinol curved needle catheter (C-Cath) design ensured 3-fold superior stem cell retention compared to a standard needle. In the setting of chronic infarction, percutaneous delivery of stem cells with C-Cath yielded a 37.7±7.1% versus 10.0±2.8% retention achieved with a traditional needle, without impact on biocompatibility or safety.
Conclusions
Modeling guided development of a nitinol-based curved needle delivery system with incremental side holes achieved enhanced myocardial stem cell retention.
doi:10.1161/CIRCINTERVENTIONS.112.000422
PMCID: PMC4273747  PMID: 24326777
stem cell retention; endocardial delivery; stem cell; catheter; regeneration; myocardial infarction
11.  Metabolic modeling of endosymbiont genome reduction on a temporal scale 
This study explores the order in which individual metabolic genes are lost in an in silico evolutionary process leading from the metabolic network of Eschericia coli to that of the genome-reduced endosymbiont Buchnera aphidicola.
Simulating the reductive evolutionary process under several growth conditions, a remarkable correlation between in silico and phylogenetically reconstructed gene loss time is obtained.A gene's k-robustness (its depth of backups) is prime determinant of its loss time.In silico gene loss time is a better predictor of their actual loss times than genomic features and network properties.Simulating the reductive evolutionary process by the loss of large blocks followed by single-gene deletions, as known to occur in evolution, yields a remarkable correspondence with the phylogenetic reconstruction and the block loss reported in the literature.
An open fundamental challenge in Systems Biology is whether a genome-scale model can predict patterns of genome evolution by realistically accounting for the associated biochemical constraints. In this study, we explore the order in which individual genes are lost in an in silico evolutionary process, leading from the metabolic network of Eschericia coli to that of the endosymbiont Buchnera aphidicola.
To evaluate the in silico gene loss time, we repeated the reductive evolutionary process introduced by Pál et al (2006), denoting the in silico deletion time of a gene in a single run of the reductive evolutionary process as the number of genes deleted before its own deletion occurred. By comparing the in silico evaluations of the gene loss time to that obtained by a phylogenetic reconstruction (Figure 1), we could evaluate the ability of an in silico process to predict temporal patterns of genome reduction. Applying this procedure on a literature-based viable media, we obtained a mean Spearman's correlation of 0.46 (53% of the maximal correlation, empirical P-value <9.9e−4) between in silico and phylogenetically reconstructed loss times. In order to provide an upper bound on evolutionary necessity stemming from metabolic constraints, we searched the space of potential growth media and biomass functions via a simulated annealing search algorithm aimed at identifying an environment/biomass function that maximizes the target correlation between in silico and reconstructed loss times. Simulating the reductive evolutionary process under the growth conditions and biomass function obtained in this process, we managed to improve the correlation between in silico and reconstructed loss times to a mean Spearman's correlation of 0.54 (63% of the maximal correlation, empirical P-value <9.9e−4, Figure 3).
Examining the dependency of the predicted loss time of each gene on its intrinsic network-level properties we find a very strong inverse Spearman's correlation of −0.84 (empirical P-value <9.9e−4) between the order of gene loss predicted in silico and the k-robustness levels of the genes, the latter denoting the depth of their functional backups in the network (Deutscher et al, 2006). Moreover, in order to examine whether the relative loss time of a gene is influenced by its functional dependencies with other genes, we performed a flux-coupling analysis and identified pairs of reactions whose activities asymmetrically depend on each other, i.e., are directionally coupled (Burgard et al, 2004). We find that genes encoding reactions whose activity is needed for activating the other reaction (and not vice versa) have a tendency to be lost later, as one would expect (binomial P-value <1e−14).
To assess the scale of these results, we examined as a control how well genomic features and network properties predict the phylogenetically reconstructed gene loss times. We examined the dependency of the latter on several factors that are known be inversely correlated with the propensity of a gene to be lost (Brinza et al, 2009; Delmotte et al, 2006; Tamames et al, 2007), including the genes' mRNA levels, tAI values (Covert et al, 2004; Reis et al, 2004; Sharp and Li, 1987; Tuller et al, 2010a) and the number of partners the gene products have in a protein–protein interaction network. Remarkably, these genomic features yield considerably lower Spearman's correlation than that obtained by the in silico simulations. Moreover, multiply regressing the loss times from the phylogenetic reconstruction on the in silico gene loss time predictions and the genomic and network variables, we found that the (normalized) coefficient of the in silico predictions in the regression is much higher than those of the genomic features, further testifying to the considerable independent predictive power of the metabolic model.
Finally, simulating the evolutionary process as large block deletions at first followed by single-gene deletions as is thought to occur in evolution (Moran and Mira, 2001; van Ham et al, 2003), a remarkable correspondence with the phylogenetic reconstruction was found. Namely, we find that after a certain amount of genes are deleted from the genome, no further block deletions can occur due to the increasing density of essential genes. Notably, the maximum amount of genes that can be deleted in blocks (i.e., until no more blocks can be deleted) corresponds to the number of genes appearing in our phylogenetic reconstruction from the LCA (last common ancestor of Buchnera and E. coli) to the LCSA (last common symbiotic ancestor, nodes 1–3 in Figure 1A), as described in the literature.
A fundamental challenge in Systems Biology is whether a cell-scale metabolic model can predict patterns of genome evolution by realistically accounting for associated biochemical constraints. Here, we study the order in which genes are lost in an in silico evolutionary process, leading from the metabolic network of Eschericia coli to that of the endosymbiont Buchnera aphidicola. We examine how this order correlates with the order by which the genes were actually lost, as estimated from a phylogenetic reconstruction. By optimizing this correlation across the space of potential growth and biomass conditions, we compute an upper bound estimate on the model's prediction accuracy (R=0.54). The model's network-based predictive ability outperforms predictions obtained using genomic features of individual genes, reflecting the effect of selection imposed by metabolic stoichiometric constraints. Thus, while the timing of gene loss might be expected to be a completely stochastic evolutionary process, remarkably, we find that metabolic considerations, on their own, make a marked 40% contribution to determining when such losses occur.
doi:10.1038/msb.2011.11
PMCID: PMC3094061  PMID: 21451589
constraint-based modeling; endosymbiont; evolution; metabolism
12.  Protective role of the lung collectins surfactant protein A and surfactant protein D in airway inflammation 
The acute inflammatory airway response is characterized by a time-dependent onset followed by active resolution. Emerging evidence suggests that epithelial cells of the proximal and distal air spaces release host defense mediators that can facilitate both the initiation and the resolution part of inflammatory airway changes. These molecules, also known as the hydrophilic surfactant proteins (surfactant protein [SP]–A and SP-D) belong to the class of collagenous lectins (collectins). The collectins are a small family of soluble pattern recognition receptors containing collagenous regions and C-type lectin domains. SP-A and SP-D are most abundant in the lung. Because of their structural uniqueness, specific localization, and functional versatility, lung collectins are important players of the pulmonary immune responses. Recent studies in our laboratory and others indicated significant associations of lung collectin levels with acute and chronic airway inflammation in both animal models and patients, suggesting the usefulness of these molecules as disease biomarkers. Research on wild-type and mutant recombinant molecules in vivo and in vitro showed that SP-A and SP-D bind carbohydrates, lipids, and nucleic acids with a broad-spectrum specificity and initiate phagocytosis of inhaled pathogens as well as apoptotic cells. Investigations on gene-deficient and conditional overexpresser mice indicated that lung collectins also directly modulate innate immune cell function and T-cell–dependent inflammatory events. Thus, these molecules have a unique, dual-function capacity to induce pathogen elimination and control proinflammatory mechanisms, suggesting a potential suitability for therapeutic prevention and treatment of chronic airway inflammation. This article reviews evidence supporting that the lung collectins play an immune-protective role and are essential for maintenance of the immunologic homeostasis in the lung.
doi:10.1016/j.jaci.2008.10.014
PMCID: PMC4097097  PMID: 19000577
Macrophage; dendritic cell; surfactant; SP-D; SP-A; innate immune regulation
13.  Bioassays to Monitor Taspase1 Function for the Identification of Pharmacogenetic Inhibitors 
PLoS ONE  2011;6(5):e18253.
Background
Threonine Aspartase 1 (Taspase1) mediates cleavage of the mixed lineage leukemia (MLL) protein and leukemia provoking MLL-fusions. In contrast to other proteases, the understanding of Taspase1's (patho)biological relevance and function is limited, since neither small molecule inhibitors nor cell based functional assays for Taspase1 are currently available.
Methodology/Findings
Efficient cell-based assays to probe Taspase1 function in vivo are presented here. These are composed of glutathione S-transferase, autofluorescent protein variants, Taspase1 cleavage sites and rational combinations of nuclear import and export signals. The biosensors localize predominantly to the cytoplasm, whereas expression of biologically active Taspase1 but not of inactive Taspase1 mutants or of the protease Caspase3 triggers their proteolytic cleavage and nuclear accumulation. Compared to in vitro assays using recombinant components the in vivo assay was highly efficient. Employing an optimized nuclear translocation algorithm, the triple-color assay could be adapted to a high-throughput microscopy platform (Z'factor = 0.63). Automated high-content data analysis was used to screen a focused compound library, selected by an in silico pharmacophor screening approach, as well as a collection of fungal extracts. Screening identified two compounds, N-[2-[(4-amino-6-oxo-3H-pyrimidin-2-yl)sulfanyl]ethyl]benzenesulfonamide and 2-benzyltriazole-4,5-dicarboxylic acid, which partially inhibited Taspase1 cleavage in living cells. Additionally, the assay was exploited to probe endogenous Taspase1 in solid tumor cell models and to identify an improved consensus sequence for efficient Taspase1 cleavage. This allowed the in silico identification of novel putative Taspase1 targets. Those include the FERM Domain-Containing Protein 4B, the Tyrosine-Protein Phosphatase Zeta, and DNA Polymerase Zeta. Cleavage site recognition and proteolytic processing of these substrates were verified in the context of the biosensor.
Conclusions
The assay not only allows to genetically probe Taspase1 structure function in vivo, but is also applicable for high-content screening to identify Taspase1 inhibitors. Such tools will provide novel insights into Taspase1's function and its potential therapeutic relevance.
doi:10.1371/journal.pone.0018253
PMCID: PMC3102056  PMID: 21647428
14.  Chemical Address Tags of Fluorescent Bioimaging Probes 
Chemical address tags can be defined as specific structural features shared by a set of bioimaging probes having a predictable influence on cell-associated visual signals obtained from these probes. Here, using a large image dataset acquired with a high content screening instrument, machine vision and cheminformatics analysis have been applied to reveal chemical address tags. With a combinatorial library of fluorescent molecules, fluorescence signal intensity, spectral, and spatial features characterizing each one of the probes' visual signals were extracted from images acquired with the three different excitation and emission channels of the imaging instrument. With multivariate regression, the additive contribution from each one of the different building blocks of the bioimaging probes towards each measured, cell-associated image-based feature was calculated. In this manner, variations in the chemical features of the molecules were associated with the resulting staining patterns, facilitating quantitative, objective analysis of chemical address tags. Hierarchical clustering and paired image-cheminformatics analysis revealed key structure-property relationships amongst many building blocks of the fluorescent molecules. The results point to different chemical modifications of the bioimaging probes that can exert similar (or different) effects on the probes' visual signals. Inspection of the clustered structures suggests intramolecular charge migration or partial charge distribution as potential mechanistic determinants of chemical address tag behavior.
doi:10.1002/cyto.a.20847
PMCID: PMC2907078  PMID: 20104576
Cheminformatics; machine vision; bioimaging; fluorescence; high content screening; image cytometry; combinatorial chemistry
15.  Emergence of tissue polarization from synergy of intracellular and extracellular auxin signaling 
Here, we provide a novel mechanistic framework for cell polarization during auxin-driven plant development that combines intracellular auxin signaling for regulation of expression of PINFORMED (PIN) auxin efflux transporters and the theoretical assumption of extracellular auxin signaling for regulation of PIN subcellular dynamics.The competitive utilization of auxin signaling component in the apoplast might account for the elusive mechanism for cell-to-cell communication for tissue polarization.Computer model simulations faithfully and robustly recapitulate experimentally observed patterns of tissue polarity and asymmetric auxin distribution during formation and regeneration of vascular systems, and during the competitive regulation of shoot branching by apical dominance.Our model generated new predictions that could be experimentally validated, highlighting a mechanistically conceivable explanation for the PIN polarization and canalization of the auxin flow in plants.
A key question of developmental biology relates to a fundamental issue in cell and tissue polarities, namely, how an individual cell in a polarized tissue senses the polarities of its neighbors and its position within tissue. In plant development, this issue is of pronounced importance, because plants have a remarkable ability to redefine cell and tissue polarities in different developmental programs, such as embryogenesis, postembryonic organogenesis, vascular tissue formation, and tissue regeneration (Kleine-Vehn and Friml, 2008).
A polar, cell-to-cell transport of the small signaling molecule auxin in conjunction with local auxin biosynthesis determines auxin gradients during embryonic and postembryonic development, giving positional cues for primordia formation, organ patterning, and tropistic growth (Friml et al, 2002; Benková et al, 2003; Reinhardt et al, 2003; Heisler et al, 2005; Scarpella et al, 2006; Dubrovsky et al, 2008). Over the past decades, theoretical models proposed that auxin acts as a polarizing cue in the center of a positive feedback mechanisms for auxin transport that has a key role in synchronized polarity rearrangements. However, the mechanistic basis for such a feedback loop between auxin and its own transport remains to a large extent elusive.
The direction of auxin transport largely depends on the polar subcellular localization of PINFORMED (PIN) proteins at the plasma membrane (Petrášek et al, 2006; Wiśniewska et al, 2006). These proteins recycle between the plasma membrane and intracellular endosomal compartments (Geldner et al, 2001; Dhonukshe et al, 2007), and their recycling modulates PIN-dependent auxin efflux rates and enable rapid changes in PIN polarity (Dubrovsky et al, 2008; Kleine-Vehn et al, 2008a). Nevertheless, the molecular basis for PIN polarization in plants remains unknown.
To gain new mechanistic insights in the hypothetical feedback mechanisms governing PIN polarization, several theoretical studies (Mitchison, 1980; Sachs, 1981; Rolland-Lagan and Prusinkiewicz, 2005; Jönsson et al, 2006; Smith et al, 2006; Merks et al, 2007; Bayer et al, 2009; Kramer, 2009) have been carried out. These models suggest that auxin promotes its own transport by modulating the amount of PIN proteins at the plasma membrane by incorporating either not yet identified flux gradient-based component (Mitchison, 1980; Rolland-Lagan and Prusinkiewicz, 2005; Bayer et al, 2009; Kramer, 2009) or an unknown short-range intercellular signal-transmitting auxin concentrations of its direct neighbors (Jönsson et al, 2006; Smith et al, 2006; Merks et al, 2007; Bayer et al, 2009; Sahlin et al, 2009).
Here, we propose a feedback driven, biologically plausible model for PIN polarization and auxin transport that introduces the combination of intracellular and extracellular auxin signaling pathways as a unified approach for tissue polarization in plants. Our computer model is based on chemiosmotic hypothesis (Goldsmith et al, 1981; Figure 1A) and integrates up-to-date experimental data, such as auxin feedback on PIN expression (Peer et al, 2004; Heisler et al, 2005) via a nuclear auxin signaling pathway (Chapman and Estelle, 2009; Figure 1B), auxin carrier recycling auxin (Dubrovsky et al, 2008; Kleine-Vehn et al, 2008a; Figure 1C), and auxin feedback on PIN endocytosis (Paciorek et al, 2005) via novel hypothetical, yet plausible, assumption of extracellular auxin perception (Figure 1D).
The heart of our extracellular receptor-based polarization (ERP) mechanism is the competitive utilization of auxin receptors in the intercellular space that allows a direct and simple cell-to-cell communication scheme. In our model, auxin binds to its extracellular receptor in the concentration-dependent manner and induces signal to modulate PIN protein abundance at the plasma membrane (Figure 1D). The direct mode of the signal transfer involves temporal immobilization of recruited receptors to the plasma membrane, which is reflected by reduced diffusion of receptors involved in auxin signaling (Figure 1D). This competitive utilization mechanism enables cell-to-cell communication in our model, leading to receptor enrichment at the site of higher auxin concentration (Figure 1D). The PIN polarization and polar auxin transport in our model both depend on and contribute to the establishment of differential auxin signaling in the cell wall. This feedback loop leads ultimately to the alignment of PIN polarization within a tissue.
We demonstrated the plausibility of the ERP model for various processes, including de novo vascularization, venation patterning, and tissue regeneration in computer simulations performed with only minimal initial assumptions, a discrete auxin source, and a distal sink. The ERP model reproduces the very detailed PIN polarization events that occur during primary vein initiation (Scarpella et al, 2006), such as basal PIN1 polarity in provascular cells, transient adverse PIN1 polarization in neighboring cells during the alignment of tissue polarization, and inner-lateral polarity displayed by the tissues surrounding a conductive auxin channel (Figure 3). Additionally, the ERP model generates high auxin concentration and high auxin flux simultaneously in emerging veins, revising the classical canalization models (Mitchison, 1980; Rolland-Lagan and Prusinkiewicz, 2005). Importantly, all our model simulations support the claim that the ERP model represents the first single approach that faithfully reproduces PIN polarization, both with the auxin gradient (basal PIN1 polarity in provascular cells) and against the auxin gradient (transient adverse PIN1 polarization in neighboring cells surrounding the provascular bundle), as well as producing the corresponding auxin distribution patterns during auxin canalization.
The proposed model introduces the extracellular auxin signaling pathway, which is crucial to account for coordinated PIN polarization and auxin distribution during venation patterning in plants. The putative candidate for extracellular auxin receptor is auxin-binding protein 1 (ABP1), which resides in the lumen of the endoplasmic reticulum and is secreted to the cell wall (Napier et al, 2002; Tromas et al, 2009) where it is physiologically active (Leblanc et al, 1999; Steffens et al, 2001). Additionally, auxin inhibits clathrin-dependent PIN internalization via binding to ABP1 (Robert et al, 2010). Thus, we speculate that the extracellular fraction of ABP1 (or additionally yet to be identified ABPs) could correspond to the common pool of extracellular auxin receptors in the ERP model. A future challenge will be to test whether the ERP model unifies complex PIN polarization and auxin distribution patterns in embryogenesis, root system maintenance, and de novo organ formation.
Plant development is exceptionally flexible as manifested by its potential for organogenesis and regeneration, which are processes involving rearrangements of tissue polarities. Fundamental questions concern how individual cells can polarize in a coordinated manner to integrate into the multicellular context. In canalization models, the signaling molecule auxin acts as a polarizing cue, and feedback on the intercellular auxin flow is key for synchronized polarity rearrangements. We provide a novel mechanistic framework for canalization, based on up-to-date experimental data and minimal, biologically plausible assumptions. Our model combines the intracellular auxin signaling for expression of PINFORMED (PIN) auxin transporters and the theoretical postulation of extracellular auxin signaling for modulation of PIN subcellular dynamics. Computer simulations faithfully and robustly recapitulated the experimentally observed patterns of tissue polarity and asymmetric auxin distribution during formation and regeneration of vascular systems and during the competitive regulation of shoot branching by apical dominance. Additionally, our model generated new predictions that could be experimentally validated, highlighting a mechanistically conceivable explanation for the PIN polarization and canalization of the auxin flow in plants.
doi:10.1038/msb.2010.103
PMCID: PMC3018162  PMID: 21179019
auxin; canalization; cell polarity; PIN proteins
16.  Intratracheal Bleomycin Causes Airway Remodeling and Airflow Obstruction in Mice 
Experimental lung research  2012;38(3):135-146.
Introduction
In addition to parenchymal fibrosis, fibrotic remodeling of the distal airways has been reported in interstitial lung diseases. Mechanisms of airway wall remodeling, which occurs in a variety of chronic lung diseases, are not well defined and current animal models are limited.
Methods
We quantified airway remodeling in lung sections from subjects with idiopathic pulmonary fibrosis (IPF) and controls. To investigate intratracheal bleomycin as a potential animal model for fibrotic airway remodeling, we evaluated lungs from C57BL/6 mice after bleomycin treatment by histologic scoring for fibrosis and peribronchial inflammation, morphometric evaluation of subepithelial connective tissue volume density, TUNEL assay, and immunohistochemistry for transforming growth factor β1 (TGFβ1), TGFβ2, and the fibroblast marker S100A4. Lung mechanics were determined at 3 weeks post-bleomycin.
Results
IPF lungs had small airway remodeling with increased bronchial wall thickness compared to controls. Similarly, bleomycin treated mice developed dose-dependent airway wall inflammation and fibrosis and greater airflow resistance after high dose bleomycin. Increased TUNEL+ bronchial epithelial cells and peribronchial inflammation were noted by 1 week, and expression of TGFβ1 and TGFβ2 and accumulation of S100A4+ fibroblasts correlated with airway remodeling in a bleomycin dose-dependent fashion.
Conclusions
IPF is characterized by small airway remodeling in addition to parenchymal fibrosis, a pattern also seen with intratracheal bleomycin. Bronchial remodeling from intratracheal bleomycin follows a cascade of events including epithelial cell injury, airway inflammation, pro-fibrotic cytokine expression, fibroblast accumulation, and peribronchial fibrosis. Thus, this model can be utilized to investigate mechanisms of airway remodeling.
doi:10.3109/01902148.2012.658595
PMCID: PMC4046254  PMID: 22394287
airway; inflammation; remodeling; mouse model; fibroblast
17.  A Multiscale Agent-Based in silico Model of Liver Fibrosis Progression 
Chronic hepatic inflammation involves a complex interplay of inflammatory and mechanical influences, ultimately manifesting in a characteristic histopathology of liver fibrosis. We created an agent-based model (ABM) of liver tissue in order to computationally examine the consequence of liver inflammation. Our liver fibrosis ABM (LFABM) is comprised of literature-derived rules describing molecular and histopathological aspects of inflammation and fibrosis in a section of chemically injured liver. Hepatocytes are modeled as agents within hexagonal lobules. Injury triggers an inflammatory reaction, which leads to activation of local Kupffer cells and recruitment of monocytes from circulation. Portal fibroblasts and hepatic stellate cells are activated locally by the products of inflammation. The various agents in the simulation are regulated by above-threshold concentrations of pro- and anti-inflammatory cytokines and damage-associated molecular pattern molecules. The simulation progresses from chronic inflammation to collagen deposition, exhibiting periportal fibrosis followed by bridging fibrosis, and culminating in disruption of the regular lobular structure. The ABM exhibited key histopathological features observed in liver sections from rats treated with carbon tetrachloride (CCl4). An in silico “tension test” for the hepatic lobules predicted an overall increase in tissue stiffness, in line with clinical elastography literature and published studies in CCl4-treated rats. Therapy simulations suggested differential anti-fibrotic effects of neutralizing tumor necrosis factor alpha vs. enhancing M2 Kupffer cells. We conclude that a computational model of liver inflammation on a structural skeleton of physical forces can recapitulate key histopathological and macroscopic properties of CCl4-injured liver. This multiscale approach linking molecular and chemomechanical stimuli enables a model that could be used to gain translationally relevant insights into liver fibrosis.
doi:10.3389/fbioe.2014.00018
PMCID: PMC4126446  PMID: 25152891
cirrhosis; computer simulation; inflammation; elastography; hepatocyte
18.  A model of yeast cell-cycle regulation based on multisite phosphorylation 
Multisite phosphorylation of CDK target proteins provides the requisite nonlinearity for cell cycle modeling using elementary reaction mechanisms.Stochastic simulations, based on Gillespie's algorithm and using realistic numbers of protein and mRNA molecules, compare favorably with single-cell measurements in budding yeast.The role of transcription–translation coupling is critical in the robust operation of protein regulatory networks in yeast cells.
Progression through the eukaryotic cell cycle is governed by the activation and inactivation of a family of cyclin-dependent kinases (CDKs) and auxiliary proteins that regulate CDK activities (Morgan, 2007). The many components of this protein regulatory network are interconnected by positive and negative feedback loops that create bistable switches and transient pulses (Tyson and Novak, 2008). The network must ensure that cell-cycle events proceed in the correct order, that cell division is balanced with respect to cell growth, and that any problems encountered (in replicating the genome or partitioning chromosomes to daughter cells) are corrected before the cell proceeds to the next phase of the cycle. The network must operate robustly in the context of unavoidable molecular fluctuations in a yeast-sized cell. With a volume of only 5×10−14 l, a yeast cell contains one copy of the gene for each component of the network, a handful of mRNA transcripts of each gene, and a few hundreds to thousands of protein molecules carrying out each gene's function. How large are the molecular fluctuations implied by these numbers, and what effects do they have on the functioning of the cell-cycle control system?
To answer these questions, we have built a new model (Figure 1) of the CDK regulatory network in budding yeast, based on the fact that the targets of CDK activity are typically phosphorylated on multiple sites. The activity of each target protein depends on how many sites are phosphorylated. The target proteins feedback on CDK activity by controlling cyclin synthesis (SBF's role) and degradation (Cdh1's role) and by releasing a CDK-counteracting phosphatase (Cdc14). Every reaction in Figure 1 can be described by a mass-action rate law, with an accompanying rate constant that must be estimated from experimental data. As the transcription and translation of mRNA molecules have major effects on fluctuating numbers of protein molecules (Pedraza and Paulsson, 2008), we have included mRNA transcripts for each protein in the model.
To create a deterministic model, the rate laws are combined, according to standard principles of chemical kinetics, into a set of 60 differential equations that govern the temporal dynamics of the control system. In the stochastic version of the model, the rate law for each reaction determines the probability per unit time that a particular reaction occurs, and we use Gillespie's stochastic simulation algorithm (Gillespie, 1976) to compute possible temporal sequences of reaction events. Accurate stochastic simulations require knowledge of the expected numbers of mRNA and protein molecules in a single yeast cell. Fortunately, these numbers are available from several sources (Ghaemmaghami et al, 2003; Zenklusen et al, 2008). Although the experimental estimates are not always in good agreement with each other, they are sufficiently reliable to populate a stochastic model with realistic numbers of molecules.
By simulating thousands of cells (as in Figure 5), we can build up representative samples for computing the mean and s.d. of any measurable cell-cycle property (e.g. interdivision time, size at division, duration of G1 phase). The excellent fit of simulated statistics to observations of cell-cycle variability is documented in the main text and Supplementary Information.
Of particular interest to us are observations of Di Talia et al (2007) of the timing of a crucial G1 event (export of Whi5 protein from the nucleus) in a population of budding yeast cells growing at a specific growth rate α=ln2/(mass-doubling time). Whi5 export is a consequence of Whi5 phosphorylation, and it occurs simultaneously with the release (activation) of SBF (see Figure 1). Using fluorescently labeled Whi5, Di Talia et al could easily measure (in individual yeast cells) the time, T1, from cell birth to the abrupt loss of Whi5 from the nucleus. Correlating T1 to the size of the cell at birth, Vbirth, they found that, for a sample of daughter cells, αT1 versus ln(Vbirth) could be fit with two straight lines of slope −0.7 and −0.3. Our simulation of this experiment (Figure 7 of the main text) compares favorably with Figure 3d and e in Di Talia et al (2007).
The major sources of noise in our model (and in protein regulatory networks in yeast cells, in general) are related to gene transcription and the small number of unique mRNA transcripts. As each mRNA molecule may instruct the synthesis of dozens of protein molecules, the coefficient of variation of molecular fluctuations at the protein level (CVP) may be dominated by fluctuations at the mRNA level, as expressed in the formula (Pedraza and Paulsson, 2008) where NM, NP denote the number of mRNA and protein molecules, respectively, and ρ=τM/τP is the ratio of half-lives of mRNA and protein molecules. For a yeast cell, typical values of NM and NP are 8 and 800, respectively (Ghaemmaghami et al, 2003; Zenklusen et al, 2008). If ρ=1, then CVP≈25%. Such large fluctuations in protein levels are inconsistent with the observed variability of size and age at division in yeast cells, as shown in the simplified cell-cycle model of Kar et al (2009) and as we have confirmed with our more realistic model. The size of these fluctuations can be reduced to a more acceptable level by assuming a shorter half-life for mRNA (say, ρ=0.1).
There must be some mechanisms whereby yeast cells lessen the protein fluctuations implied by transcription–translation coupling. Following Pedraza and Paulsson (2008), we suggest that mRNA gestation and senescence may resolve this problem. Equation (3) is based on a simple, one-stage, birth–death model of mRNA turnover. In Supplementary Appendix 1, we show that a model of mRNA processing, with 10 stages each of mRNA gestation and senescence, gives reasonable fluctuations at the protein level (CVP≈5%), even if the effective half-life of mRNA is 10 min. A one-stage model with τM=1 min gives comparable fluctuations (CVP≈5%). In the main text, we use a simple birth–death model of mRNA turnover with an ‘effective' half-life of 1 min, in order to limit the computational complexity of the full cell-cycle model.
In order for the cell's genome to be passed intact from one generation to the next, the events of the cell cycle (DNA replication, mitosis, cell division) must be executed in the correct order, despite the considerable molecular noise inherent in any protein-based regulatory system residing in the small confines of a eukaryotic cell. To assess the effects of molecular fluctuations on cell-cycle progression in budding yeast cells, we have constructed a new model of the regulation of Cln- and Clb-dependent kinases, based on multisite phosphorylation of their target proteins and on positive and negative feedback loops involving the kinases themselves. To account for the significant role of noise in the transcription and translation steps of gene expression, the model includes mRNAs as well as proteins. The model equations are simulated deterministically and stochastically to reveal the bistable switching behavior on which proper cell-cycle progression depends and to show that this behavior is robust to the level of molecular noise expected in yeast-sized cells (∼50 fL volume). The model gives a quantitatively accurate account of the variability observed in the G1-S transition in budding yeast, which is governed by an underlying sizer+timer control system.
doi:10.1038/msb.2010.55
PMCID: PMC2947364  PMID: 20739927
bistability; cell-cycle variability; size control; stochastic model; transcription–translation coupling
19.  A Novel Small Molecule Antagonist of Choline Kinase-α That Simultaneously Suppresses MAPK and PI3K/AKT Signaling 
Oncogene  2011;30(30):3370-3380.
Choline kinase-α expression and activity are increased in multiple human neoplasms as a result of growth factor stimulation and activation of cancer-related signaling pathways. The product of choline kinase-α, phosphocholine, serves as an essential metabolic reservoir for the production of phosphatidylcholine, the major phospholipid constituent of membranes and substrate for the production of lipid second messengers. Using in silico screening for small molecules that may interact with the choline kinase-α substrate binding domain, we identified a novel competitive inhibitor, N-(3,5-dimethylphenyl)-2-[[5-(4-ethylphenyl)-1H-1,2,4-triazol-3-yl]sulfanyl] acetamide (termed CK37) that inhibited purified recombinant human choline kinase-α activity, reduced the steady-state concentration of phosphocholine in transformed cells, and selectively suppressed the growth of neoplastic cells relative to normal epithelial cells. Choline kinase-α activity is required for the downstream production of phosphatidic acid, a promoter of several Ras signaling pathways. CK37 suppressed MAPK and PI3K/AKT signaling, disrupted actin cytoskeletal organization, and reduced plasma membrane ruffling. Finally, administration of CK37 significantly decreased tumor growth in a lung tumor xenograft mouse model, suppressed tumor phosphocholine, and diminished activating phosphorylations of ERK and AKT in vivo. Together, these results further validate choline kinase-α as a molecular target for the development of agents that interrupt Ras signaling pathways, and indicate that receptor-based computational screening should facilitate the identification of new classes of choline kinase-α inhibitors.
doi:10.1038/onc.2011.51
PMCID: PMC3136659  PMID: 21423211
Chemotherapy; Choline Kinase; Metabolism; In silico; Phosphocholine
20.  Signature-Based Small Molecule Screening Identifies Cytosine Arabinoside as an EWS/FLI Modulator in Ewing Sarcoma 
PLoS Medicine  2007;4(4):e122.
Background
The presence of tumor-specific mutations in the cancer genome represents a potential opportunity for pharmacologic intervention to therapeutic benefit. Unfortunately, many classes of oncoproteins (e.g., transcription factors) are not amenable to conventional small-molecule screening. Despite the identification of tumor-specific somatic mutations, most cancer therapy still utilizes nonspecific, cytotoxic drugs. One illustrative example is the treatment of Ewing sarcoma. Although the EWS/FLI oncoprotein, present in the vast majority of Ewing tumors, was characterized over ten years ago, it has never been exploited as a target of therapy. Previously, this target has been intractable to modulation with traditional small-molecule library screening approaches. Here we describe a gene expression–based approach to identify compounds that induce a signature of EWS/FLI attenuation. We hypothesize that screening small-molecule libraries highly enriched for FDA-approved drugs will provide a more rapid path to clinical application.
Methods and Findings
A gene expression signature for the EWS/FLI off state was determined with microarray expression profiling of Ewing sarcoma cell lines with EWS/FLI-directed RNA interference. A small-molecule library enriched for FDA-approved drugs was screened with a high-throughput, ligation-mediated amplification assay with a fluorescent, bead-based detection. Screening identified cytosine arabinoside (ARA-C) as a modulator of EWS/FLI. ARA-C reduced EWS/FLI protein abundance and accordingly diminished cell viability and transformation and abrogated tumor growth in a xenograft model. Given the poor outcomes of many patients with Ewing sarcoma and the well-established ARA-C safety profile, clinical trials testing ARA-C are warranted.
Conclusions
We demonstrate that a gene expression–based approach to small-molecule library screening can identify, for rapid clinical testing, candidate drugs that modulate previously intractable targets. Furthermore, this is a generic approach that can, in principle, be applied to the identification of modulators of any tumor-associated oncoprotein in the rare pediatric malignancies, but also in the more common adult cancers.
Todd Golub and colleagues show that a gene expression-based screen of small-molecule libraries can identify candidate drugs that modulate cancer-associated oncoproteins.
Editors' Summary
Background.
Cancer occurs when cells accumulate genetic changes (mutations) that allow them to divide uncontrollably and to travel throughout the body (metastasize). Chemotherapy, a mainstay of cancer treatments, works by killing rapidly dividing cells. Because some normal tissues also contain dividing cells and are therefore sensitive to chemotherapy drugs, it is hard to treat cancer without causing serious side effects. In recent years, however, researchers have identified some of the mutations that drive the growth of cancer cells. This raises the possibility of designing drugs that kill only cancer cells by specifically targeting “oncoproteins” (the abnormal proteins generated by mutations that transform normal cells into cancer cells). Some “targeted” drugs have already reached the clinic, but unfortunately medicinal chemists do not know how to inhibit the function of many classes of oncoproteins with the small organic molecules that make the best medicines. One oncoprotein in this category is EWS/FLI. This contains part of a protein called EWS fused to part of a transcription factor (a protein that controls cell behavior by telling the cell which proteins to make) called FLI. About 80% of patients with Ewing sarcoma (the second commonest childhood cancer of bone and soft tissue) have the mutation responsible for EWS/FLI expression. Localized Ewing sarcoma can be treated with nontargeted chemotherapy (often in combination with surgery and radiotherapy), but treatment for recurrent or metastatic disease remains very poor.
Why Was This Study Done?
Researchers have known for years that EWS/FLI expression drives the development of Ewing sarcoma by activating the expression of target genes needed for tumor formation. However, EWS/FLI has never been exploited as a target for therapy of this cancer—mainly because traditional approaches used to screen libraries of small molecules do not identify compounds that modulate the activity of transcription factors. In this study, the researchers have used a new gene expression–based, high-throughput screening (GE-HTS) approach to identify compounds that modulate the activity of EWS/FLI.
What Did the Researchers Do and Find?
The researchers used a molecular biology technique called microarray expression profiling to define a 14-gene expression signature that differentiates between Ewing sarcoma cells in which the EWS/FLI fusion protein is active and those in which it is inactive. They then used this signature to screen a library of about 1,000 chemicals (many already approved for other clinical uses) in a “ligation-mediated amplification assay.” For this, the researchers grew Ewing sarcoma cells with the test chemicals, extracted RNA from the cells, and generated a DNA copy of the RNA. They then added two short pieces of DNA (probes) specific for each signature gene to the samples. In samples that expressed a given signature gene, both probes bound and were then ligated (joined together) and amplified. Because one of each probe pair also contained a unique “capture sequence,” the signature genes expressed in each sample were finally identified by adding colored fluorescent beads, each linked to DNA complementary to a different capture sequence. The most active modulator of EWS/FLI activity identified by this GE-HTS approach was cytosine arabinoside (ARA-C). At levels achievable in people, this compound reduced the abundance of EWS/FLI protein in and the viability and cancer-like behavior of Ewing sarcoma cells growing in test tubes. ARA-C treatment also slowed the growth of Ewing sarcoma cells transplanted into mice.
What Do These Findings Mean?
These findings identify ARA-C, which is already used to treat children with some forms of leukemia, as a potent modulator of EWS/FLI activity. More laboratory experiments are needed to discover how ARA-C changes the behavior of Ewing sarcoma cells. Nevertheless, given the poor outcomes currently seen in many patients with Ewing sarcoma and the historical reluctance to test new drugs in children, these findings strongly support the initiation of clinical trials of ARA-C in children with Ewing sarcoma. These results also show that the GE-HTS approach is a powerful way to identify candidate drugs able to modulate the activity of some of the oncoproteins (including transcription factors and other previously intractable targets) that drive cancer development.
Additional Information.
Please access these Web sites via the online version of this summary at http://dx.doi.org/10.1371/journal.pmed.0040122.
Cancerquest from Emory University, provides information on cancer biology (also includes information in Spanish, Chinese and Russian)
The MedlinePlus encyclopedia has pages on Ewing sarcoma
Information for patients and health professionals on Ewing sarcoma is available from the US National Cancer Institute
Cancerbackup offers information for patients and their parents on Ewing sarcoma
Wikipedia has pages on DNA microarrays and expression profiling (note that Wikipedia is a free online encyclopedia that anyone can edit)
doi:10.1371/journal.pmed.0040122
PMCID: PMC1851624  PMID: 17425403
21.  A Versatile Technique for the In Vivo Imaging of Human Tumor Xenografts Using Near-Infrared Fluorochrome-Conjugated Macromolecule Probes 
PLoS ONE  2013;8(12):e82708.
Here, we present a versatile method for detecting human tumor xenografts in vivo, based on the enhanced permeability and retention (EPR) effect, using near-infrared (NIR) fluorochrome-conjugated macromolecule probes. Bovine serum albumin (BSA) and two immunoglobulins—an anti-human leukocyte antigen (HLA) monoclonal antibody and isotype control IgG2a—were labeled with XenoLight CF770 fluorochrome and used as NIR-conjugated macromolecule probes to study whole-body imaging in a variety of xenotransplantation mouse models. NIR fluorescent signals were observed in subcutaneously transplanted BxPC-3 (human pancreatic cancer) cells and HCT 116 (colorectal cancer) cells within 24 h of NIR-macromolecule probe injection, but the signal from the fluorochrome itself or from the NIR-conjugated small molecule (glycine) injection was not observed. The accuracy of tumor targeting was confirmed by the localization of the NIR-conjugated immunoglobulin within the T-HCT 116 xenograft (in which the orange-red fluorescent protein tdTomato was stably expressed by HCT 116 cells) in the subcutaneous transplantation model. However, there was no significant difference in the NIR signal intensity of the region of interest between the anti-HLA antibody group and the isotype control group in the subcutaneous transplantation model. Therefore, the antibody accumulation within the tumor in vivo is based on the EPR effect. The liver metastasis generated by an intrasplenic injection of T-HCT 116 cells was clearly visualized by the NIR-conjugated anti-HLA probe but not by the orange-red fluorescent signal derived from the tdTomato reporter. This result demonstrated the superiority of the NIR probes over the tdTomato reporter protein at enhancing tissue penetration. In another xenograft model, patient-derived xenografts (PDX) of LC11-JCK (human non-small cell lung cancer) were successfully visualized using the NIR-conjugated macromolecule probe without any genetic modification. These results suggested that NIR-conjugated macromolecule, preferably, anti-HLA antibody probe is a valuable tool for the detection of human tumors in experimental metastasis models using whole-body imaging.
doi:10.1371/journal.pone.0082708
PMCID: PMC3866180  PMID: 24358218
22.  Development of fluorescent probes for bioimaging applications 
Fluorescent probes, which allow visualization of cations such as Ca2+, Zn2+ etc., small biomolecules such as nitric oxide (NO) or enzyme activities in living cells by means of fluorescence microscopy, have become indispensable tools for clarifying functions in biological systems. This review deals with the general principles for the design of bioimaging fluorescent probes by modulating the fluorescence properties of fluorophores, employing mechanisms such as acceptor-excited Photoinduced electron Transfer (a-PeT), donor-excited Photoinduced electron Transfer (d-PeT), and spirocyclization, which have been established by our group. The a-PeT and d-PeT mechanisms are widely applicable for the design of bioimaging probes based on many fluorophores and the spirocyclization process is also expected to be useful as a fluorescence off/on switching mechanism. Fluorescence modulation mechanisms are essential for the rational design of novel fluorescence probes for target molecules. Based on these mechanisms, we have developed more than fifty bioimaging probes, of which fourteen are commercially available. The review also describes some applications of the probes developed by our group to in vitro and in vivo systems.
doi:10.2183/pjab.86.837
PMCID: PMC3037519  PMID: 20948177
probe; bioimaging; photoinduced electron transfer; fluorescence; spirocyclization
23.  Developing a multiscale, multi-resolution agent-based brain tumor model by graphics processing units 
Multiscale agent-based modeling (MABM) has been widely used to simulate Glioblastoma Multiforme (GBM) and its progression. At the intracellular level, the MABM approach employs a system of ordinary differential equations to describe quantitatively specific intracellular molecular pathways that determine phenotypic switches among cells (e.g. from migration to proliferation and vice versa). At the intercellular level, MABM describes cell-cell interactions by a discrete module. At the tissue level, partial differential equations are employed to model the diffusion of chemoattractants, which are the input factors of the intracellular molecular pathway. Moreover, multiscale analysis makes it possible to explore the molecules that play important roles in determining the cellular phenotypic switches that in turn drive the whole GBM expansion. However, owing to limited computational resources, MABM is currently a theoretical biological model that uses relatively coarse grids to simulate a few cancer cells in a small slice of brain cancer tissue. In order to improve this theoretical model to simulate and predict actual GBM cancer progression in real time, a graphics processing unit (GPU)-based parallel computing algorithm was developed and combined with the multi-resolution design to speed up the MABM. The simulated results demonstrated that the GPU-based, multi-resolution and multiscale approach can accelerate the previous MABM around 30-fold with relatively fine grids in a large extracellular matrix. Therefore, the new model has great potential for simulating and predicting real-time GBM progression, if real experimental data are incorporated.
doi:10.1186/1742-4682-8-46
PMCID: PMC3312859  PMID: 22176732
24.  A multiscale MDCT image-based breathing lung model with time-varying regional ventilation 
A novel algorithm is presented that links local structural variables (regional ventilation and deforming central airways) to global function (total lung volume) in the lung over three imaged lung volumes, to derive a breathing lung model for computational fluid dynamics simulation. The algorithm constitutes the core of an integrative, image-based computational framework for subject-specific simulation of the breathing lung. For the first time, the algorithm is applied to three multi-detector row computed tomography (MDCT) volumetric lung images of the same individual. A key technique in linking global and local variables over multiple images is an in-house mass-preserving image registration method. Throughout breathing cycles, cubic interpolation is employed to ensure C1 continuity in constructing time-varying regional ventilation at the whole lung level, flow rate fractions exiting the terminal airways, and airway deformation. The imaged exit airway flow rate fractions are derived from regional ventilation with the aid of a three-dimensional (3D) and one-dimensional (1D) coupled airway tree that connects the airways to the alveolar tissue. An in-house parallel large-eddy simulation (LES) technique is adopted to capture turbulent-transitional-laminar flows in both normal and deep breathing conditions. The results obtained by the proposed algorithm when using three lung volume images are compared with those using only one or two volume images. The three-volume-based lung model produces physiologically-consistent time-varying pressure and ventilation distribution. The one-volume-based lung model under-predicts pressure drop and yields un-physiological lobar ventilation. The two-volume-based model can account for airway deformation and non-uniform regional ventilation to some extent, but does not capture the non-linear features of the lung.
doi:10.1016/j.jcp.2012.12.007
PMCID: PMC3685439  PMID: 23794749
Multiscale; pulmonary air flow; boundary condition; regional ventilation; image registration; MDCT
25.  Theoretical models for the simulation of particle deposition and tracheobronchial clearance in lungs of patients with chronic bronchitis 
Introduction
Based upon theoretical models particle deposition and clearance in human respiratory systems affected by chronic bronchitis can be approximated reliably. As a consequence of those hypothetical results, optimal frame conditions (e.g., inhalation time and volume, particle properties) for inhalation therapies can be determined.
Methods
Simulation of particle deposition was conducted by modelling a partly or fully obstructed tracheobronchial architecture. Bronchitis-induced reductions of the airway calibres were computed by application of specific scaling factors. Three different scenarios of chronic bronchitis were modelled. Brownian motion, inertial impaction, interception, and gravitational settling were assumed as main deposition forces influencing inhaled particular mass. Tracheobronchial clearance was approximated by application of generation-specific mucus velocities as well as the consideration of a slow bronchial clearance phase, whose half-time varied between 5 and 20 days.
Results
Under different breathing conditions (i.e., sitting and light-work breathing) deposition of submicron and µm-sized particles is significantly enhanced within the bronchial lung region, but also alveolar deposition becomes partly enhanced. By changing the inhalation conditions target sites of therapeutic aerosols may be reached with rather high accuracy. Based on the data of this modified models, particle retention in lung airways of patients suffering from chronic bronchitis may be noticeably prolonged, with 24-hour retention values being increased by up to 50%.
Discussion and conclusions
As exhibited by the results, particle deposition behaviour in lungs affected by chronic bronchitis differs remarkably from that in healthy lungs. These theoretical finds are mostly supported by experimental data. Further, experimental and theoretical deposition results may be used for an estimation of the grade of disease. Tracheobronchial clearance reduces its efficiency with each progress of the disease which increases the probability of bacterial infections in the airways.
doi:10.3978/j.issn.2305-5839.2012.11.02
PMCID: PMC4200694  PMID: 25332949
Stochastic lung model; particle deposition; bronchial clearance; chronic bronchitis

Results 1-25 (1226439)