Sensors used in precision agriculture for the detection of heavy metals in irrigation water are generally expensive and sometimes their deployment and maintenance represent a permanent investment to keep them in operation, leaving a lasting polluting footprint in the environment at the end of their lifespan. This represents an area of opportunity to design new biological devices that can replace part, or all of the sensors currently used. In this article, a novel workflow is proposed to fully carry out the complete process of design, modeling, and simulation of reprogrammable microorganisms . As a proof-of-concept, the workflow has been used to design three whole-cell biosensors for the detection of heavy metals in irrigation water, namely arsenic, mercury and lead. These biosensors are in compliance with the concentration limits established by the World Health Organization (WHO). The proposed workflow allows the design of a wide variety of completely biodevices, which aids in solving problems that cannot be easily addressed with classical computing. The workflow is based on two technologies typical of synthetic biology: the design of synthetic genetic circuits, and synthetic engineering, which allows us to address the design of reprogrammable microorganisms using software and hardware to develop theoretical models. These models enable the behavior prediction of complex biological systems. The output of the workflow is then exported in the form of complete genomes in SBOL, GenBank and FASTA formats, enabling their subsequent implementation in a laboratory. The present proposal enables professionals in the area of computer science to collaborate in biotechnological processes from a theoretical perspective previously or complementary to a design process carried out directly in the laboratory by molecular biologists. Therefore, key results pertaining to this work include the fully workflow that leads to designs that can be tested in the lab or , and a proof-of-concept of how the workflow generates synthetic circuits in the form of three whole-cell heavy metal biosensors that were designed, modeled and simulated using the workflow. The simulations carried out show realistic spatial distributions of biosensors reacting to different concentrations (zero, low and threshold level) of heavy metal presence and at different growth phases (stationary and exponential) that are backed up by the whole design and modeling phases of the workflow.