Hey everyone, I am coding a stochastic model in python and I have hit a wall and don't know what to do or how to code something.
This model is to stochastically simulate the interactions between 2 molecules, molecule1 and molecule 2, taking into account their transcription, splicing, and translation rates (and degradation). It is meant to have an ON/OFF switch that activates/deactivates the transcription of the genes in question. My issue is with the probability of the switch being ON. I want it to have a base rate and increase with the amount of protein present (say, increases 10% every 10 proteins with a cap on 30%). There is a loop of positive autofeedback for each molecule, as well as cross-feedback.
If anyone has any idea on how to code for this I would appreciate a comment. I can share the full code if requested (don't want to share the complete thing here, as it would make the post massive).
Thanks a lot!
I would imagine this is something one of the AI code tools could handle.