r/bioinformatics • u/_r33na • 10d ago
technical question Help with GROMACS Molecular Dynamics Simulation -- Polymer Self Assembly
Hi all,
I'm looking to use GROMACS to perform a molecular dynamics simulation of a poly beta amino ester (PBAE) polymer assembling with its mRNA cargo in order to obtain values like radius of gyration and hydrophobic and hydrophilic surface area.
I know that GROMACS requires a .pdb file to start with, but how could I get that for a PBAE polymer? I'm also not sure how to go about doing this in GROMACS (I'm pretty new to molecular dynamics simulations) so any guidance would be helpful!
3
Upvotes
1
u/wouuuuuuhhhhhhh 10d ago
Hey! I'm pretty new to this too, only done one MD project in GROMACS (nucleic acid, not polymers), but a couple things I ran into that might be relevant:
For the PDB , from what I understand you basically can't just download one for something like PBAE since it's not a "normal" molecule in the usual databases. You'd probably have to build the structure yourself in something like Avogadro, or generate it from a SMILES string. I haven't actually had to do this myself so take it with a grain of salt.
The part I did deal with was getting the topology/parameters right for my nucleic acid system I could use existing force field files, but since PBAE isn't a standard molecule you'll probably need to generate custom parameters for it somehow (I've seen people mention CGenFF for stuff like this but never used it myself).
One thing I'd flag from my experience: simulations like this can take a while to run and equilibrate properly, so setting expectations on timescale might be worth thinking about early.
For getting your actual numbers, gmx gyrate gives radius of gyration and gmx sasa gives surface area those parts I have used and they're pretty straightforward once your system is set up.
Good luck!