Biomedical subjects

Jeremy S Edwards

Publications and source records attributed to Jeremy S Edwards.

At least 19 recordsLinked to original sources

PubMed2006

Spatial modeling of dimerization reaction dynamics in the plasma membrane: Monte Carlo vs. continuum differential equations.

Bimolecular reactions in the plasma membrane, such as receptor dimerization, are a key signaling step for many signaling systems. For receptors to dimerize, they must first diffuse until a collision happens, upon which a dimerization reaction may occur. Therefore, study of the dynamics of cell signaling on the membrane may require the use of a spatial modeling framework. Despite the availability of spatial simulation methods, e.g., stochastic spatial Monte Carlo (MC) simulation and partial differential equation (PDE) based approaches, many biological models invoke well-mixed assumptions without completely evaluating the importance of spatial organization. Whether one is to utilize a spatial or non-spatial simulation framework is therefore an important decision. In order to evaluate the importance of spatial effects a priori, i.e., without performing simulations, we have assessed the applicability of a dimensionless number, known as second Damköhler number (Da), defined here as the ratio of time scales of collision and reaction, for 2-dimensional bimolecular reactions. Our study shows that dimerization reactions in the plasma membrane with Da approximately >0.1 (tested in the receptor density range of 10(2)-10(5)/microm(2)) require spatial modeling. We also evaluated the effective reaction rate constants of MC and simple deterministic PDEs. Our simulations show that the effective reaction rate constant decreases with time due to time dependent changes in the spatial distribution of receptors. As a result, the effective reaction rate constant of simple PDEs can differ from that of MC by up to two orders of magnitude. Furthermore, we show that the fluctuations in the number of copies of signaling proteins (noise) may also depend on the diffusion properties of the system. Finally, we used the spatial MC model to explore the effect of plasma membrane heterogeneities, such as receptor localization and reduced diffusivity, on the dimerization rate. Interestingly, our simulations show that localization of epidermal growth factor receptor (EGFR) can cause the diffusion limited dimerization rate to be up to two orders of magnitude higher at higher average receptor densities reported for cancer cells, as compared to a normal cell.

Algorithms↗

PubMed2006

Designing experiments that aid in the identification of regulatory networks.

Predictive mathematical models of the interactions of a genetic network can provide insight into the mechanisms of gene regulation, the role of various genes within a network and how multiple genes interact leading to complex traits. However, identification of the parameters and interactions is currently a limiting step in the development of such models. This work reviews the state of the art for design of experiments in biological systems and demonstrates the need for improved design of experiments through the use of a model system. Appropriate design of experiments has a profound impact on the ability to identify a model and on the quality of resulting identified model. Key issues include the selection of appropriate input sequences (e.g. random, independent multivariate inputs) and the selection of the sampling frequencies. This work demonstrates that these issues are especially important in the identification of biochemical networks and that the traditional biochemical approach is incapable of truly identifying the behavior present in such networks.

Computational Biology↗

PubMed2005

Parallel analysis of mutant human glucose 6-phosphate dehydrogenase in yeast using PCR colonies.

We demonstrate a highly parallel strategy to analyze the impact of single nucleotide mutations on protein function. Using our method, it is possible to screen a population and quickly identify a subset of functionally interesting mutants. Our method utilizes a combination of yeast functional complementation, growth competition of mutant pools, and polymerase colonies. A defined mutant human glucose-6-phosphate-dehydrogenase library was constructed which contains all possible single nucleotide missense mutations in the eight-residue glucose-6-phosphate binding peptide of the enzyme. Mutant human enzymes were expressed in a zwf1 (gene encoding yeast homologue) deletion strain of Saccharomyces cerevisiae. Growth rates of the 54 mutant strains arising from this library were measured in parallel in conditions selective for active hG6PD. Several residues were identified which tolerated no mutations (Asp200, His201 and Lys205) and two (Ile199 and Leu203) tolerated several substitutions. Arg198, Tyr202, and Gly204 tolerated only 1-2 specific substitutions. Generalizing from the positions of tolerated and non-tolerated amino acid substitutions, hypotheses were generated about the functional role of specific residues, which could, potentially, be tested using higher resolution/lower throughput methods.

Amino Acid Sequence↗

PubMed2005

Computational modeling reveals molecular details of epidermal growth factor binding.

BACKGROUND: The ErbB family of receptors are dysregulated in a number of cancers, and the signaling pathway of this receptor family is a critical target for several anti-cancer drugs. Therefore a detailed understanding of the mechanisms of receptor activation is critical. However, despite a plethora of biochemical studies and recent single particle tracking experiments, the early molecular mechanisms involving epidermal growth factor (EGF) binding and EGF receptor (EGFR) dimerization are not as well understood. Herein, we describe a spatially distributed Monte Carlo based simulation framework to enable the simulation of in vivo receptor diffusion and dimerization. RESULTS: Our simulation results are in agreement with the data from single particle tracking and biochemical experiments on EGFR. Furthermore, the simulations reveal that the sequence of receptor-receptor and ligand-receptor reaction events depends on the ligand concentration, receptor density and receptor mobility. CONCLUSION: Our computer simulations reveal the mechanism of EGF binding on EGFR. Overall, we show that spatial simulation of receptor dynamics can be used to gain a mechanistic understanding of receptor activation which may in turn enable improved cancer treatments in the future.

Computer Simulation↗

PubMed2005

Heterogeneities in EGF receptor density at the cell surface can lead to concave up scatchard plot of EGF binding.

The mechanism responsible for the concave up nature of the Scatchard plot of epidermal growth factor (EGF) binding on EGF receptor (EGFR) has been a controversial issue for more than a decade. Past efforts to mechanistically simulate the concave up nature of the Scatchard plot of EGF binding have shown that negative cooperativity in EGF binding on an EGFR dimer or inclusion of some external site or binding event can describe this behavior. However, herein we show that heterogeneity in the density of EGFR due to localization in certain regions of the plasma membrane, which has been experimentally reported, can also lead to concave up shape of the Scatchard plot of the EGF binding on EGFR.

Cell Line, Tumor↗

PubMed2005

Time accelerated Monte Carlo simulations of biological networks using the binomial tau-leap method.

UNLABELLED: Developing a quantitative understanding of intracellular networks requires simulations and computational analyses. However, traditional differential equation modeling tools are often inadequate due to the stochasticity of intracellular reaction networks that can potentially influence the phenotypic characteristics. Unfortunately, stochastic simulations are computationally too intense for most biological systems. Herein, we have utilized the recently developed binomial tau-leap method to carry out stochastic simulations of the epidermal growth factor receptor induced mitogen activated protein kinase cascade. Results indicate that the binomial tau-leap method is computationally 100-1000 times more efficient than the exact stochastic simulation algorithm of Gillespie. Furthermore, the binomial tau-leap method avoids negative populations and accurately captures the species populations along with their fluctuations despite the large difference in their size. AVAILABILITY: http://www.dion.che.udel.edu/multiscale/Introduction.html. Fortran 90 code available for academic use by email. SUPPLEMENTARY INFORMATION: Details about the binomial tau-leap algorithm, software and a manual are available at the above website.

Algorithms↗

PubMed2005

Estimating optimal profiles of genetic alterations using constraint-based models.

Metabolic engineering involves application of recombinant DNA methods to manipulate metabolic networks to improve cellular properties. It is critical that the genetic alterations be performed in an optimal manner to maximize profit. In addition to the product yield, productivity consideration is also critical, especially for the production of bulk chemicals such as 1,3-propanediol. In this work, we demonstrate that it is suboptimal from the standpoint of productivity to induce genetic alteration at the start of the production process. A bi-level optimization scheme is formulated to determine the optimal temporal flux profile for the manipulated reaction. In the first case study, an optimal flux in the reaction catalyzed by glycerol kinase is determined to maximize the glycerol production at the end of a 6-h batch cultivation of Escherichia coli under aerobic conditions. The final glycerol concentration is 30% higher for the optimal flux profile compared with having an active flux during the entire batch. The effect of the mass transfer coefficient on the optimal profile and the glycerol concentration is also determined. In the second case study, the anaerobic batch fermentation of the ldh(-) strain of Escherichia coli is considered. The optimal flux in the acetate pathway is determined to maximize the final ethanol concentration. The optimal flux results in higher ethanol concentration (11.92 mmol L(-1)) compared to strains with no acetate flux (8.36 mmol L(-1)) and fully active acetate flux (6.22 mmol L(-1)). We also examine the effects of growth inhibition due to high ethanol concentrations and variations in final batch time on ethanol production.

Algorithms↗

PubMed2004

Autoregulation of regulatory proteins is key for dynamic operation of GAL switch in Saccharomyces cerevisiae.

Autoregulation and nucleocytoplasmic shuttling play important roles in the operation of the GAL regulatory system. However, the significance of these mechanisms in the overall operation of the switch is unclear. In this work, we develop a dynamic model for the GAL system and further validate the same using steady-state and dynamic experimental expression data. Next, the model is used to delineate the relevance of shuttling and autoregulation in response to inducing, repressing, and non-inducing-non-repressing media. The analysis indicates that autoregulation of the repressor, Gal80p, is key in obtaining three distinct steady states in response to the three media. In particular, the analysis rationalizes the intuitively paradoxical observation that the concentration of repressor, Gal80p, actually increases in response to an increase in the inducer concentration. On the other hand, although nucleocytoplasmic shuttling does not affect the dynamics of the system, it plays a dominant role in obtaining a sensitive response to galactose. The dynamic model was also used to obtain insights on the preculturing effect on the system behavior.

Binding Sites↗

PubMed2004

MAPK cascade possesses decoupled controllability of signal amplification and duration.

The three important characteristics of the output signal of mitogen activated protein kinase (MAPK) cascade are time delay between stimulus and response, amplitude gain, and duration of the output signal. In this study, we performed a sensitivity analysis on the computational model of epidermal growth factor receptor (EGFR) activated MAPK cascade developed by Schoeberl and co-workers (1) to identify the sensitive steps of the pathway affecting these characteristics. We show that the signaling network is sensitive in a decoupled manner, which provides the ability to control its output amplitude and duration one at a time. Signal duration is found sensitive only to the phosphatase reactions at the MEK level. In contrast, signal amplitude is found most sensitive to the phosphatase reactions at the ERK level. Time delay is found to be a robust characteristic of the system.

Computational Biology↗

PubMed2004

Digital quantitative measurements of gene expression.

One of the primary goals of functional genomics is to provide a quantitative understanding of gene function. However, the success of this enterprise is dependent on the accuracy and precision of the functional genomic data. A novel approach, digital analysis of gene expression (DAGE) described herein, is an accurate and precise technology for measuring digital gene expression on a relative or absolute scale by simply counting the number of transcripts of a gene being expressed at a given time. The result is a greatly improved technology sensitive enough for identifying and quantifying small (but biologically important and statistically relevant) changes in gene expression. Fourteen genes involved in galactose metabolism in Saccharomyces cerevisiae were analyzed for their expression levels in glucose and galactose minimal media. The quantitative expression results were characterized in terms of distributional and accuracy attributes; they were also in general agreement (in terms of direction of change) with corresponding results obtained using microarray technology. DAGE is likely to have profound implications in the field of functional genomics because the gene expression measurements are digital in nature and therefore more accurate than any other technologies.

Algorithms↗

PubMed2004

An integrated mechanistic model for transcription-coupled nucleotide excision repair.

Preferential repair of the transcribed strand of active genes is usually attributed to a coupling protein that dislodges RNA polymerase stalled at a damage site and recruits repair enzymes. Experimental observations of the effect of transcription on preferential repair in Escherichia coli are contradictory and inexplicable by this model. In this study, it is argued that the multiple conformations displayed by a stalled RNA polymerase result in two sub-pathways for repair: Mfd coupled and direct. Together with the fact that RNA polymerase recruits the repair enzymes in a promoter dependent manner, an integrated mechanistic model is proposed that is capable of explaining the effect of transcription on preferential repair reported in literature. The quantitative behavior of the model is illustrated by describing the various reactions using a biochemical network. The implications of the model on the mechanism for transcription-coupled repair in higher organisms are briefly discussed.

DNA Repair↗

PubMed2004

Detection of allelic variations of human gene expression by polymerase colonies.

BACKGROUND: Quantification of variations of human gene expression is complicated by the small differences between different alleles. Recent work has shown that variations do exist in the relative allelic expression levels in certain genes of heterozygous individuals. Herein, we describe the application of an immobilized polymerase chain reaction technique as an alternative approach to measure relative allelic differential expression. RESULTS: Herein, we report a novel assay, based on immobilized polymerase colonies, that accurately quantifies the relative expression levels of two alleles in a given sample. Mechanistically, this was accomplished by PCR amplifying a gene in a cDNA library in a thin polyacrylamide gel. By immobilizing the PCR, it is ensured that each transcript gives rise to only a single immobilized PCR colony, or "polony". Once polony amplified, the two alleles of the gene were differentially labeled by performing in situ sequencing with fluorescently labeled nucleotides. For these sets of experiments, silent single nucleotide polymorphisms (SNPs) were used to discriminate the two alleles. Finally, a simple count was then performed on the differentially labeled polonies in order to determine the relative expression levels of the two alleles. To validate this technique, the relative expression levels of PKD2 in a family of heterozygous patients bearing the 4208G/A SNP were examined and compared to the literature. CONCLUSIONS: We were able to reproduce the results of allelic variation in gene expression using an accurate technology known as polymerase colonies. Therefore, we have demonstrated the utility of this method in human gene expression analysis.

Alleles↗

PubMed2004

Robust global sensitivity in multiple enzyme cascade system explains how the downstream cascade structure may remain unaffected by cross-talk.

A steady-state framework was applied to the ubiquitous tricyclic enzyme cascade structure, as seen in the mitogen-activated protein (MAP) kinase system, to analyze the effect of upstream kinase concentrations on final output response. The results suggest that signal amplification achieved by the cascade structure ensured that the modifying enzymes at various steps of the cascade were nearly saturated. Thus, there was no change in the response sensitivity with increasing upstream kinase concentration. Analysis was also extended to branching of a signaling pathway as an example of cross-talk. It was observed that the cascade structure confers a larger share of the signal transduction properties to its last kinase. This phenomenon in enzyme cascades may explain how the response of the terminal MAP kinase is unaffected by cross-talk of upstream kinases.

Adenosine Triphosphate↗

PubMed2004

RecA mediated initial alignment of homologous DNA molecules displays apparent first order kinetics with little effect of heterology.

The mechanism and determinants of RecA mediated initial alignment of homologous DNA molecules were studied by performing Monte Carlo simulations of the dynamics of DNA molecules. The simulation procedure was used to assess the effect of heterologous DNA and dilution on the rate of formation and yield of homologous alignments. The results show that the apparent first order kinetic behavior and the impact of heterologous DNA, reported in literature [J. Biol. Chem. 261 (1986) 1025], can be observed even if the conversion of the initially aligned molecules into a stable joint is not rate-determining. The present study is the first step towards developing rigorous computational models to describe the process of homologous recombination, and theoretical frameworks to retrieve biophysical parameters of strand pairing and exchange proteins from in vitro assays of joint molecule formation.

Computer Simulation↗

PubMed2004

Assaying gene function by growth competition experiment.

High-throughput screening and analysis is one of the emerging paradigms in biotechnology. In particular, high-throughput methods are essential in the field of functional genomics because of the vast amount of data generated in recent and ongoing genome sequencing efforts. In this report we discuss integrated functional analysis methodologies which incorporate both a growth competition component and a highly parallel assay used to quantify results of the growth competition. Several applications of the two most widely used technologies in the field, i.e., transposon mutagenesis and deletion strain library growth competition, and individual applications of several developing or less widely reported technologies are presented.

Animals↗

PubMed2004

Equilibrium analysis of allosteric interactions shows zero-order effects.

The binding of effector to an allosteric protein exhibits a non-Michaelis-Menten behavior, resulting in either ultrasensitive or subsensitive response. In the present work, a modular approach has been developed to determine the response curve for allosteric systems at higher concentration of allosteric enzyme than that of effector (zero-order sensitivity, as observed in enzyme cascades) by equilibrium analysis. The analysis shows that, in an allosteric system, the zero-order effect can make the response ultrasensitive or subsensitive with respect to the enzyme concentration. The response is dependent on the number of binding sites, cooperativity, and the total effector concentration. The framework was further applied to a well studied allosteric protein, the Escherichia coli aspartate transcarbamoylase. The predictions are found to be consistent with the reported experimental data.

Allosteric Regulation↗

PubMed2003

Parallel competition analysis of Saccharomyces cerevisiae strains differing by a single base using polymerase colonies.

We describe a strategy to analyze the impact of single nucleotide mutations on protein function. Our method utilizes a combination of yeast functional complementation, growth competition of mutant pools and polyacrylamide gel immobilized PCR. A system was constructed in which the yeast PGK1 gene was expressed from a plasmid-borne copy of the gene in a PGK1 deletion strain of Saccharomyces cerevisiae. Using this system, we demonstrated that the enrichment or depletion of PGK1 point mutants from a mixed culture was consistent with the expected results based on the isolated growth rates of the mutants. Enrichment or depletion of individual point mutants was shown to result from increases or decreases, respectively, in the specific activities of the encoded proteins. Further, we demonstrate the ability to analyze the functional effect of many individual point mutations in parallel. By functional complementation of yeast deletions with human homologs, our technique could be readily applied to the functional analysis of single nucleotide polymorphisms in human genes of medical interest.

Biopolymers↗

PubMed2003

Advances in flux balance analysis.

Biology is going through a paradigm shift from reductionist to holistic, systems-based approaches. The complete genome sequence for a number of organisms is available and the analysis of genome sequence data is proving very useful. Thus, genome sequencing projects and bioinformatic analyses are leading to a complete 'parts catalog' of the molecular components in many organisms. The next challenge will be to reconstruct and simulate overall cellular functions based on the extensive reductionist information. Recent advances have been made in the area of flux balance analysis, a mathematical modeling approach often utilized by metabolic engineers to quantitatively simulate microbial metabolism.

Computational Biology↗