When developing molecular dynamics simulations, it is common to utilise existing files or assets found from databases. However, when modelling more specific or complex systems such as membrane proteins, it is likely that such models are not built nor characterised well. For our project, we faced many issues with existing bioinformatics tools, which made it imperative to develop our simulation. Fusion proteins, especially surface-expressed proteins containing linkers, often yield unreliable and inconclusive results due to the limitations of models like AlphaFold. So while AlphaFold helps to democratise structure prediction by generating 3D structures from merely a sequence, these outputs lack contextualisation and thus are mostly reliable as snapshots for static proteins and making them less dependable for more complex structures such as membrane proteins.
In the case of bioremediation, this issue is even more exacerbated as conditions are seldom ideal or able to be safely tested in the scope of an iGEM project. Thus, Molecular Dynamics outputs such as GROMACS simulations aim to inch closer towards providing context to engineered proteins. However, these tools face two problems:
- Computational Stress serves as a bottleneck for many teams or young undergraduate researchers
- Complex Terminal Commands that intimidate and complicate the workflow, which can create a steep learning curve
In this part, we hope to highlight our contribution to the second problem by providing our troubleshooting tips and troubleshooting errors we ran into along the way. In our modelling, we developed literature-accurate bacteria membranes for Pseudomonas putida and Escherichia coli and embedded fusion proteins into them using CharmmGUI. We also subsequently show how to import these models into GROMACS and simulate them in metal ion environments, and measure the stability.
We hope our documentation can show other teams how to do the same and overcome the many troubles we encountered during our learning process.
For more information, please visit our Model page