pablo-arantes / making-it-rain

Cloud-based molecular simulations for everyone
MIT License
385 stars 105 forks source link

Gromacs_inputs #38

Closed cwilliamszf closed 2 years ago

cwilliamszf commented 2 years ago

I am new to molecular dynamics and I am trying to do some simulations of mutant GPCRs in membranes. To create a system for your Colab notebook for simulation I used GPCRmodsim to get a short simulation and necessary GROMACS files. Unfortunately, when I run the simulation, I get the following error:

Too few fields in [ dihedraltypes ] line: SI OS 1 0.000 3.766 3

Any help would be appreciated.

I have included the gro and top files and also for good measure the FF files. WTactive.zip FFfiles.zip

cwilliamszf commented 2 years ago

Addendum here is the screenshot of the error. image

pablo-arantes commented 2 years ago

HI @cwilliamszf,

Sorry for my delay. Regarding your question, I've checked your files and seems to me that topologies from GPCRmodsim are a little bit different from GROMACS topologies. I believe this is your major issue here. If you could create the topology directly in GROMACS and try again, it would be better. I can understand that GPCRmodsim topologies work well in GROMACS, however we are using OpenMM engine to run MD simulation. With that in mind, OpenMM is very restricted in the file format.

Sorry for that, I hope you understand.

Thank you.

Best,

Pablo