Off-campus UMass Amherst users: To download campus access dissertations, please use the following link to log into our proxy server with your UMass Amherst user name and password.

Non-UMass Amherst users: Please talk to your librarian about requesting this dissertation through interlibrary loan.

Dissertations that have an embargo placed on them will not be available to anyone until the embargo expires.

Author ORCID Identifier


Open Access Dissertation

Document Type


Degree Name

Doctor of Philosophy (PhD)

Degree Program

Chemical Engineering

Year Degree Awarded


Month Degree Awarded


First Advisor

Michael A. Henson

Subject Categories

Partial Differential Equations | Systems Biology | Transport Phenomena


Biofilms are ubiquitous in medical, environmental, and engineered microbial systems. The majority of naturally occurring microbes grow as mixed species biofilms. These complicated biofilm consortia are comprised of many cell phenotypes with complex interactions and self-organized into three-dimensional structures. Approximately 2% of the US population suffers from non-healing chronic wounds infected by a combination of commensal and pathogenic bacteria whereas about 500,000 cases of Clostridium difficile infections (CDI) are reported annually. These polymicrobial infections are often resilient to antibiotic treatment due to the nutrient-rich environments and species interactions that promote community stability and robustness. This thesis focusses on developing metabolic modeling framework to study the interactions and the spatial/temporal organizations in the biofilms. The modeling framework is based on integrating genome scale metabolic reconstructions of considered species in this work, with the nutrient uptake kinetics to predict the species abundances, growth rates and byproduct secretions. The spatiotemporal modeling framework accounts for the nutrient concentration gradients in the biofilm system. Spatiotemporal biofilm metabolic models were formulated by combining genome scale metabolic reconstructions of considered species with uptake kinetics for available nutrients and reaction-diffusion type equations for species biomass, supplied substrates and synthesized metabolic byproducts. The resulting partial differential equations embedded with linear programs were discretized in the space and integrated using a dynamic flux balance method. This framework was used to calculate the spatial and temporal variations in the species, nutrient and byproduct concentrations in biofilms. This framework was used to study the species organization and dynamics in chronic wound infections, CDI and environmental biofilms. The chronic wound biofilm model was comprising of two most dominant species, Pseudomonas aeruginosa and Staphylococcus aureus. The CDI biofilm model was comprising of representative species from three most common phyla in gut Bacteroidetes thetaiotaomicron, Faecalibacterium prausnitzii, Escherichia coli and pathogen C. difficile. The simulation results were used to study the interspecies interactions, the spatial partitioning in the biofilms and important crossfeeding relationships within the community. These predictions would be useful in devising effective antibiotic treatment strategies to cure the biofilm infections associated with chronic wounds and C. difficile. The environmental biofilm model for cyanobacteria and heterotrophs was developed and validated with the experimental results, this model was used to evaluate the community dynamics under extreme environmental conditions The second modeling framework considered biofilm as a well-mixed homogenous system at steady state. Steady state in silico community models were formulated by combining genome scale metabolic reconstructions of the considered species. The community models were solved using SteadyCom method. This method uses community flux balance analysis to calculate the relative abundance of each species with an objective of maximizing the community growth rate. A 12 species chronic wound community metabolic model covering 74% of 16S rDNA pyrosequencing reads of dominant genera from 2,963 chronic wound patients was developed. The community model was used to predict species abundances averaged across this large patient population. The simulation results from this study were used to identify putative mutualistic interactions between bacteria that could be targeted to enhance treatment efficacy. The frameworks developed in this thesis would be useful in developing patient/disease specific therapeutic treatments.