We used two different mathematical modelling approaches to describe and simulate different aspects of our project. We built an agent-based model to understand the behaviour of nitrogen fixing bacteria in response to the chemoattractant naringenin. We used SimBiotics to visualise stochastic simulations via real-time animations. These models guided experimental work and were then informed by data from our chemotaxis experiments and bacterial growth characterisation. Our model provided insight into the biofilm formation process, including biofilm thickness and number of cells of each nitrogen-fixing species present. We also built a kinetic model describing metabolic flux through the naringenin biosynthetic pathway. Our model employed mass action kinetics to describe the behaviour of reactants and products for each step in the pathway. By coupling this information with models describing the rates of production and turnover of the four naringenin biosynthetic enzymes we developed an improved genetic design for our biosynthetic devices.