Hi OpenMM!
I am simulating an antiparallel beta sheet that likes to fall apart during equilibration. I need to prevent this. I am wondering if anyone on this forum has any advice for keeping the beta sheet together. I was thinking to use CustomNonbondedForce between a few particular atoms in the adjacent strands to keep them together. Do you recommend a flat bottom or harmonic potential? What other advice do you have for someone trying to do this? I realize what I'm simulating is somewhat artificial but that's fine for my purposes right now. Please advise anything you can think of!
Thank you!
Dan
retaining antiparallel beta sheet
- Daniel Konstantinovsky
- Posts: 77
- Joined: Tue Jun 11, 2019 12:21 pm
- Peter Eastman
- Posts: 2588
- Joined: Thu Aug 09, 2007 1:25 pm
Re: retaining antiparallel beta sheet
If you want a constraint between a few specific atoms, you probably want a CustomBondForce, not CustomNonbonedForce. As for questions about potential shape, strength, etc. I'll leave that to people with more experience doing that sort of thing.
- Daniel Konstantinovsky
- Posts: 77
- Joined: Tue Jun 11, 2019 12:21 pm
Re: retaining antiparallel beta sheet
Oh ok, thank you, I will use a bond force.