S-Leaping: An Adaptive, Accelerated Stochastic Simulation Algorithm, Bridging [Formula: see text]-Leaping and R-Leaping.

We propose the S-leaping algorithm for the acceleration of Gillespie's stochastic simulation algorithm that combines the advantages of the two main accelerated methods; the [Formula: see text]-leaping and R-leaping algorithms. These algorithms are known to be efficient under different conditions; the [Formula: see text]-leaping is efficient for non-stiff systems or systems with partial equilibrium, while the R-leaping performs better in stiff system thanks to an efficient sampling procedure. However, even a small change in a system's set up can critically affect the nature of the simulated system and thus reduce the efficiency of an accelerated algorithm. The proposed algorithm combines the efficient time step selection from the [Formula: see text]-leaping with the effective sampling procedure from the R-leaping algorithm. The S-leaping is shown to maintain its efficiency under different conditions and in the case of large and stiff systems or systems with fast dynamics, the S-leaping outperforms both methods. We demonstrate the performance and the accuracy of the S-leaping in comparison with the [Formula: see text]-leaping and R-leaping on a number of benchmark systems involving biological reaction networks.

Price of disorder in the lac repressor hinge helix.

The Lac system of genes has been pivotal in understanding gene regulation. When the lac repressor protein binds to the correct DNA sequence, the hinge region of the protein goes through a disorder to order transition. The structure of this region of the protein is well understood when it is in this bound conformation, but less so when it is not. Structural studies show that this region is flexible. Our simulations show this region is extremely flexible in solution; however, a high concentration of salt can help kinetically trap the hinge helix. Thermodynamically, disorder is more favorable without the DNA present.

Integrated In Vivo Genotoxicity Assessment of Procarbazine Hydrochloride Demonstrates Induction of Pig-a and LacZ Mutations, and Micronuclei, in MutaMouse Hematopoietic Cells.

Procarbazine hydrochloride (PCH) is a DNA-reactive hematopoietic carcinogen with potent and well-characterized clastogenic activity. However, there is a paucity of in vivo mutagenesis data for PCH, and in vitro assays often fail to detect the genotoxic effects of PCH due to the complexity of its metabolic activation. We comprehensively evaluated the in vivo genotoxicity of PCH on hematopoietic cells of male MutaMouse transgenic rodents using a study design that facilitated assessments of micronuclei and Pig-a mutation in circulating erythrocytes, and lacZ mutant frequencies in bone marrow. Mice were orally exposed to PCH (0, 6.25, 12.5, and 25 mg/kg/day) for 28 consecutive days. Blood samples collected 2 days after cessation of treatment exhibited significant dose-related induction of micronuclei in both immature and mature erythrocytes. Bone marrow and blood collected 3 and 70 days after cessation of treatment also showed significantly elevated mutant frequencies in both the lacZ and Pig-a assays even at the lowest dose tested. PCH-induced lacZ and Pig-a (immature and mature erythrocytes) mutant frequencies were highly correlated, with R2 values ≥0.956, with the exception of lacZ vs. Pig-a mutants in mature erythrocytes at the 70-day time point (R2 = 0.902). These results show that PCH is genotoxic in vivo and demonstrate that the complex metabolism and resulting genotoxicity of PCH is best evaluated in intact animal models. Our results further support the concept that multiple biomarkers of genotoxicity, especially hematopoietic cell genotoxicity, can be readily combined into one study provided that adequate attention is given to manifestation times. Environ. Mol. Mutagen. 60:505-512, 2019. © 2018 Her Majesty the Queen in Right of Canada.

Integration of sperm DNA damage assessment into OECD test guidelines for genotoxicity testing using the MutaMouse model.

The Organisation for Economic Co-operation and Development (OECD) endorses test guidelines (TG) for identifying chemicals that are genotoxic, such as the transgenic rodent gene mutation assay (TG 488). Current OECD TG do not include assays for sperm DNA damage resulting in a critical testing gap. We evaluated the performance of the Sperm Chromatin Structure Assay (SCSA) and the Terminal Deoxynucleotidyl Transferase-Mediated Deoxyuridine Triphosphate Nick end Labeling (TUNEL) assay to detect sperm DNA damage within the recommended TG 488 protocol. MutaMouse males received 0, 0.5, 1, or 2 mg/kg/day triethylenemelamine (TEM), a multifunctional alkylating agent, for 28 days orally and tissues were collected two (blood) and three (sperm and bone marrow) days later. TEM significantly increased the frequency of lacZ mutants in bone marrow, and of micronuclei (MN) in both reticulocytes (%MN-RET) and normochromatic erythrocytes (%MN-NCE) in a dose-dependent manner (P < 0.05). The percentage of DNA fragmentation index (%DFI) and %TUNEL positive cells demonstrated dose-related increases in sperm (P < 0.05), and the two assay results were strongly correlated (R = 0.9298). Within the same animal, a good correlation was observed between %MN-NCE and %DFI (R = 0.7189). Finally, benchmark dose modelling (BMD) showed comparable BMD10 values among the somatic and germ cell assays. Our results suggest that sperm DNA damage assays can be easily integrated into standard OECD designs investigating genotoxicity in somatic tissues to provide key information on whether a chemical is genotoxic in germ cells and impact its risk assessment.

Coarse-grained analysis of stochastically simulated cell populations with a positive feedback genetic network architecture.

Among the different computational approaches modelling the dynamics of isogenic cell populations, discrete stochastic models can describe with sufficient accuracy the evolution of small size populations. However, for a systematic and efficient study of their long-time behaviour over a wide range of parameter values, the performance of solely direct temporal simulations requires significantly high computational time. In addition, when the dynamics of the cell populations exhibit non-trivial bistable behaviour, such an analysis becomes a prohibitive task, since a large ensemble of initial states need to be tested for the quest of possibly co-existing steady state solutions. In this work, we study cell populations which carry the lac operon network exhibiting solution multiplicity over a wide range of extracellular conditions (inducer concentration). By adopting ideas from the so-called "equation-free" methodology, we perform systems-level analysis, which includes numerical tasks such as the computation of coarse steady state solutions, coarse bifurcation analysis, as well as coarse stability analysis. Dynamically stable and unstable macroscopic (population level) steady state solutions are computed by means of bifurcation analysis utilising short bursts of fine-scale simulations, and the range of bistability is determined for different sizes of cell populations. The results are compared with the deterministic cell population balance model, which is valid for large populations, and we demonstrate the increased effect of stochasticity in small size populations with asymmetric partitioning mechanisms.

Quantitation of interactions between two DNA loops demonstrates loop domain insulation in E. coli cells.

Eukaryotic gene regulation involves complex patterns of long-range DNA-looping interactions between enhancers and promoters, but how these specific interactions are achieved is poorly understood. Models that posit other DNA loops--that aid or inhibit enhancer-promoter contact--are difficult to test or quantitate rigorously in eukaryotic cells. Here, we use the well-characterized DNA-looping proteins Lac repressor and phage λ CI to measure interactions between pairs of long DNA loops in E. coli cells in the three possible topological arrangements. We find that side-by-side loops do not affect each other. Nested loops assist each other's formation consistent with their distance-shortening effect. In contrast, alternating loops, where one looping element is placed within the other DNA loop, inhibit each other's formation, thus providing clear support for the loop domain model for insulation. Modeling shows that combining loop assistance and loop interference can provide strong specificity in long-range interactions.

A simple and sensitive biosensor strain for detecting toxoflavin using ß-galactosidase activity.

In this study, we developed a simple and sensitive biosensor for the determination of toxoflavin (which is toxic to various plants, fungi, animals, and bacteria) in natural samples based on ß-galactosidase activity. The proposed toxoflavin detection method for toxin-producing bacteria or toxin-contaminated foods is simple and cost effective. Burkholderia glumae, a species known to cause rice grain rot and wilt in various field crops, produces toxoflavin under the control of a LysR-type transcriptional regulator ToxR and its ligand toxoflavin. As the expression of toxoflavin biosynthetic genes requires toxoflavin as a co-activator of ToxR, a novel biosensor stain was constructed based on lacZ reporter gene integration into the first gene of the toxoflavin biosynthesis operon, toxABCDE of B. glumae. The biosensor was composed of a sensor strain (COK71), substrates (X-gal or ONPG), and culture medium, without any complex preparation process. We demonstrated that the biosensor strain is highly specific to toxoflavin, and can quantify relative amounts of toxoflavin compared with known concentrations of toxoflavin. The proposed method was reliable and simple; samples containing 50-500 nM of toxoflavin could be analyzed. More importantly, the proposed biosensor strain could identify toxoflavin-producing bacteria in real samples. The excellent performance of this biosensor is useful for diagnostic purposes, such as detecting toxoflavin-contaminated foods and environmental samples.

Interplay of protein and DNA structure revealed in simulations of the lac operon.

The E. coli Lac repressor is the classic textbook example of a protein that attaches to widely spaced sites along a genome and forces the intervening DNA into a loop. The short loops implicated in the regulation of the lac operon suggest the involvement of factors other than DNA and repressor in gene control. The molecular simulations presented here examine two likely structural contributions to the in-vivo looping of bacterial DNA: the distortions of the double helix introduced upon association of the highly abundant, nonspecific nucleoid protein HU and the large-scale deformations of the repressor detected in low-resolution experiments. The computations take account of the three-dimensional arrangements of nucleotides and amino acids found in crystal structures of DNA with the two proteins, the natural rest state and deformational properties of protein-free DNA, and the constraints on looping imposed by the conformation of the repressor and the orientation of bound DNA. The predicted looping propensities capture the complex, chain-length-dependent variation in repression efficacy extracted from gene expression studies and in vitro experiments and reveal unexpected chain-length-dependent variations in the uptake of HU, the deformation of repressor, and the folding of DNA. Both the opening of repressor and the presence of HU, at levels approximating those found in vivo, enhance the probability of loop formation. HU affects the global organization of the repressor and the opening of repressor influences the levels of HU binding to DNA. The length of the loop determines whether the DNA adopts antiparallel or parallel orientations on the repressor, whether the repressor is opened or closed, and how many HU molecules bind to the loop. The collective behavior of proteins and DNA is greater than the sum of the parts and hints of ways in which multiple proteins may coordinate the packaging and processing of genetic information.

Detection of ß-galactosidase activity: X-gal staining.

X-gal staining is a rapid and convenient histochemical technique used to detect reporter gene expression. A prerequisite is the creation or acquisition of transgenic reporter mouse lines, in which the bacterial LacZ gene has been knocked into the gene of interest or placed under the control of regulatory elements corresponding to the gene of interest. Expression is marked by a dark blue stain and can be detected at the single cell level, providing a robust visual readout of gene expression in the developing kidney. Here, we describe the methodology, applications, and limitations of this technique.

Cost-benefit tradeoffs in engineered lac operons.

Cells must balance the cost and benefit of protein expression to optimize organismal fitness. The lac operon of the bacterium Escherichia coli has been a model for quantifying the physiological impact of costly protein production and for elucidating the resulting regulatory mechanisms. We report quantitative fitness measurements in 27 redesigned operons that suggested that protein production is not the primary origin of fitness costs. Instead, we discovered that the lac permease activity, which relates linearly to cost, is the major physiological burden to the cell. These findings explain control points in the lac operon that minimize the cost of lac permease activity, not protein expression. Characterizing similar relationships in other systems will be important to map the impact of cost/benefit tradeoffs on cell physiology and regulation.

Proliferating primary hepatocytes from the pUR288 lacZ plasmid mouse are valuable tools for genotoxicity assessment in vitro.

Safety assessments of substances with regard to genotoxicity are generally based on a combination of in vitro and in vivo tests. These tests are performed according to a (tiered) test strategy whereby a positive result in vitro usually triggers further testing in vivo. A low specificity and high frequency of irrelevant positive results associated with most in vitro mammalian cell genotoxicity assays necessitates the design and validation of suitable alternatives. As such, we examined the feasibility of culturing primary hepatocytes from the pUR288 lacZ reporter mouse, and moreover, using established cultures to reliably assess genotoxic activity in vitro. Initial studies characterizing the metabolic capacity of proliferating lacZ primary hepatocytes indicated that these cells retained at least some activities important for xenobiotic metabolism: cytochrome P450 1A1 enzyme activities were markedly increased in the hepatocytes after exposure to benzo[a]pyrene, and also UDP-glucuronosyl transferase and glutathione-S-transferase activities, both Phase II enzymes, were detected. Increasing levels of phosphorylated p53 at residue serine 389 after ultraviolet treatment indicated a properly functioning p53, one of the criteria for an effective new test system. Four genotoxic substances with different mechanisms of genotoxicity, i.e., benzo[a]pyrene, bleomycin, etoposide, and cyclophosphamide, were tested in the lacZ rescue assay. For etoposide and cyclophosphamide, the induction of mutant colonies was rather low. Exposure to benzo[a]pyrene and bleomycin, however, yielded a clear concentration-dependent induction of the lacZ mutant frequency. Based on our preliminary observations, proliferating lacZ primary hepatocytes are a promising new tool for the assessment of genotoxic hazard.

Deterministic and stochastic population-level simulations of an artificial lac operon genetic network.

BACKGROUND: The lac operon genetic switch is considered as a paradigm of genetic regulation. This system has a positive feedback loop due to the LacY permease boosting its own production by the facilitated transport of inducer into the cell and the subsequent de-repression of the lac operon genes. Previously, we have investigated the effect of stochasticity in an artificial lac operon network at the single cell level by comparing corresponding deterministic and stochastic kinetic models. RESULTS: This work focuses on the dynamics of cell populations by incorporating the above kinetic scheme into two Monte Carlo (MC) simulation frameworks. The first MC framework assumes stochastic reaction occurrence, accounts for stochastic DNA duplication, division and partitioning and tracks all daughter cells to obtain the statistics of the entire cell population. In order to better understand how stochastic effects shape cell population distributions, we develop a second framework that assumes deterministic reaction dynamics. By comparing the predictions of the two frameworks, we conclude that stochasticity can create or destroy bimodality, and may enhance phenotypic heterogeneity. CONCLUSIONS: Our results show how various sources of stochasticity act in synergy with the positive feedback architecture, thereby shaping the behavior at the cell population level. Further, the insights obtained from the present study allow us to construct simpler and less computationally intensive models that can closely approximate the dynamics of heterogeneous cell populations.

The rate of the molecular clock and the cost of gratuitous protein synthesis.

BACKGROUND: The nature of the protein molecular clock, the protein-specific rate of amino acid substitutions, is among the central questions of molecular evolution. Protein expression level is the dominant determinant of the clock rate in a number of organisms. It has been suggested that highly expressed proteins evolve slowly in all species mainly to maintain robustness to translation errors that generate toxic misfolded proteins. Here we investigate this hypothesis experimentally by comparing the growth rate of Escherichia coli expressing wild type and misfolding-prone variants of the LacZ protein. RESULTS: We show that the cost of toxic protein misfolding is small compared to other costs associated with protein synthesis. Complementary computational analyses demonstrate that there is also a relatively weaker, but statistically significant, selection for increasing solubility and polarity in highly expressed E. coli proteins. CONCLUSIONS: Although we cannot rule out the possibility that selection against misfolding toxicity significantly affects the protein clock in species other than E. coli, our results suggest that it is unlikely to be the dominant and universal factor determining the clock rate in all organisms. We find that in this bacterium other costs associated with protein synthesis are likely to play an important role. Interestingly, our experiments also suggest significant costs associated with volume effects, such as jamming of the cellular environment with unnecessary proteins.

Using deep sequencing to characterize the biophysical mechanism of a transcriptional regulatory sequence.

Cells use protein-DNA and protein-protein interactions to regulate transcription. A biophysical understanding of this process has, however, been limited by the lack of methods for quantitatively characterizing the interactions that occur at specific promoters and enhancers in living cells. Here we show how such biophysical information can be revealed by a simple experiment in which a library of partially mutated regulatory sequences are partitioned according to their in vivo transcriptional activities and then sequenced en masse. Computational analysis of the sequence data produced by this experiment can provide precise quantitative information about how the regulatory proteins at a specific arrangement of binding sites work together to regulate transcription. This ability to reliably extract precise information about regulatory biophysics in the face of experimental noise is made possible by a recently identified relationship between likelihood and mutual information. Applying our experimental and computational techniques to the Escherichia coli lac promoter, we demonstrate the ability to identify regulatory protein binding sites de novo, determine the sequence-dependent binding energy of the proteins that bind these sites, and, importantly, measure the in vivo interaction energy between RNA polymerase and a DNA-bound transcription factor. Our approach provides a generally applicable method for characterizing the biophysical basis of transcriptional regulation by a specified regulatory sequence. The principles of our method can also be applied to a wide range of other problems in molecular biology.

In vivo assessment of a tissue-engineered vascular graft combining a biodegradable elastomeric scaffold and muscle-derived stem cells in a rat model.

Limited autologous vascular graft availability and poor patency rates of synthetic grafts for bypass or replacement of small-diameter arteries remain a concern in the surgical community. These limitations could potentially be improved by a tissue engineering approach. We report here our progress in the development and in vivo testing of a stem-cell-based tissue-engineered vascular graft for arterial applications. Poly(ester urethane)urea scaffolds (length = 10 mm; inner diameter = 1.2 mm) were created by thermally induced phase separation (TIPS). Compound scaffolds were generated by reinforcing TIPS scaffolds with an outer electrospun layer of the same biomaterial (ES-TIPS). Both TIPS and ES-TIPS scaffolds were bulk-seeded with 10 x 10(6) allogeneic, LacZ-transfected, muscle-derived stem cells (MDSCs), and then placed in spinner flask culture for 48 h. Constructs were implanted as interposition grafts in the abdominal aorta of rats for 8 weeks. Angiograms and histological assessment were performed at the time of explant. Cell-seeded constructs showed a higher patency rate than the unseeded controls: 65% (ES-TIPS) and 53% (TIPS) versus 10% (acellular TIPS). TIPS scaffolds had a 50% mechanical failure rate with aneurysmal formation, whereas no dilation was observed in the hybrid scaffolds. A smooth-muscle-like layer of cells was observed near the luminal surface of the constructs that stained positive for smooth muscle alpha-actin and calponin. LacZ+ cells were shown to be engrafted in the remodeled construct. A confluent layer of von Willebrand Factor-positive cells was observed in the lumen of MDSC-seeded constructs, whereas acellular controls showed platelet and fibrin deposition. This is the first evidence that MDSCs improve patency and contribute to the remodeling of a tissue-engineered vascular graft for arterial applications.

In vivo genotoxicity of EMS: statistical assessment of the dose response curves.

EMS induced micronuclei and lacZ mutations in in vivo studies in mice with a clearly sublinear dose dependency. As reported elsewhere in this issue, NOEL dose values of between 25 mg/kg/day and 80 mg/kg/day were observed for the different endpoints and tissues analysed. Here we show that statistical assessment of the data provides solid support that the induction of mutagenic and clastogenic effects after in vivo treatment with the directly DNA damaging mutagen EMS adheres to a thresholded dose response relation. These data corroborate similar evidence obtained in in vitro studies. We conclude that cells are fully capable of repairing large amounts of DNA ethylations induced by EMS without experiencing elevated mutation frequencies. The stochastic, linear risk assessment model generally employed for DNA damaging genotoxins can therefore be refuted for EMS. While presently this conclusion cannot be generalized to other genotoxins a change of paradigm appears to be indicated at least for alkylating agents inducing a comparable type and spectrum of DNA lesions as EMS.

Safety assessment of intradiscal gene therapy II: effect of dosing and vector choice.

STUDY DESIGN: Clinical, biochemical, and histologic analysis was performed after in vivo delivery of cDNA encoding various anabolic cytokines and marker genes to the lumbar epidural space of New Zealand white rabbits, using both adenoviral and adeno-associated viral vectors. OBJECTIVE: To mimic errant or misplaced doses of gene therapy to better ascertain the potential risks associated with alternative vectors and transgene products with regard to their application to problems of the intervertebral disc. SUMMARY OF BACKGROUND DATA: Work done with several anabolic cytokines including bone morphogenic proteins and transforming growth factors, has demonstrated the potential of gene therapy. Recently, data has been published demonstrating that improperly dosed or delivered adenoviral-mediated gene therapy within the subarachnoid space can result in significant morbidity in rabbits. There are currently no studies examining the effect of these errors within the epidural space or using an adeno-associated viral (AAV) vector. METHODS: Using either adenoviral or AAV vectors, complementary DNA (cDNA) encoding anabolic cytokines bone morphogenic protein-2 (BMP-2) and transforming growth factor-beta 1 and marker proteins LacZ and green fluorescent protein were injected into the epidural space of 37 New Zealand white rabbits at the L5/6 level. Rabbits were then observed clinically for up to 6 weeks, after which the rabbits were sacrificed in order to perform a comprehensive biochemical and histologic analysis. RESULTS: Following adenoviral-mediated delivery of anabolic cytokine cDNA, up to eighty percent of rabbits demonstrated significant clinical, biochemical, and histologic morbidity. Conversely, AAV-mediated delivery of any cDNA and adenoviral-mediated delivery of marker protein cDNA resulted in no clinical, histologic, or biochemical morbidity. CONCLUSION: Properly dosed and directed gene therapy seems to be both safe and potentially efficacious. This study suggests that side effects of gene therapy may be due to a combination of dosing, transgene product, and vector choice, and that newer AAV vectors may reduce these side-effects and decrease the risk of this technology.

The cost of expression of Escherichia coli lac operon proteins is in the process, not in the products.

Transcriptional regulatory networks allow bacteria to express proteins only when they are needed. Adaptive hypotheses explaining the evolution of regulatory networks assume that unneeded expression is costly and therefore decreases fitness, but the proximate cause of this cost is not clear. We show that the cost in fitness to Escherichia coli strains constitutively expressing the lactose operon when lactose is absent is associated with the process of making the lac gene products, i.e., associated with the acts of transcription and/or translation. These results reject the hypotheses that regulation exists to prevent the waste of amino acids in useless protein or the detrimental activity of unnecessary proteins. While the cost of the process of protein expression occurs in all of the environments that we tested, the expression of the lactose permease could be costly or beneficial, depending on the environment. Our results identify the basis of a single selective pressure likely acting across the entire E. coli transcriptome.

Cost-benefit theory and optimal design of gene regulation functions.

Cells respond to the environment by regulating the expression of genes according to environmental signals. The relation between the input signal level and the expression of the gene is called the gene regulation function. It is of interest to understand the shape of a gene regulation function in terms of the environment in which it has evolved and the basic constraints of biological systems. Here we address this by presenting a cost-benefit theory for gene regulation functions that takes into account temporally varying inputs in the environment and stochastic noise in the biological components. We apply this theory to the well-studied lac operon of E. coli. The present theory explains the shape of this regulation function in terms of temporal variation of the input signals, and of minimizing the deleterious effect of cell-cell variability in regulatory protein levels. We also apply the theory to understand the evolutionary tradeoffs in setting the number of regulatory proteins and for selection of feed-forward loops in genetic circuits. The present cost-benefit theory can be used to understand the shape of other gene regulatory functions in terms of environment and noise constraints.

Stochastic and deterministic simulations of heterogeneous cell population dynamics.

A Monte Carlo algorithm, which can accurately simulate the dynamics of entire heterogeneous cell populations, was developed. The algorithm takes into account the random nature of cell division as well as unequal partitioning of cellular material at cell division. Moreover, it is general in the sense that it can accommodate a variety of single-cell, deterministic reaction kinetics as well as various stochastic division and partitioning mechanisms. The validity of the algorithm was assessed through comparison of its results with those of the corresponding deterministic cell population balance model in cases where stochastic behavior is expected to be quantitatively negligible. Both algorithms were applied to study: (a) linear intracellular kinetics and (b) the expression dynamics of a genetic network with positive feedback architecture, such as the lac operon. The effects of stochastic division as well as those of different division and partitioning mechanisms were assessed in these systems, while the comparison of the stochastic model with a continuum model elucidated the significance of cell population heterogeneity even in cases where only the prediction of average properties is of primary interest.