Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

DNA is not implemented yet? #521

Open
songyongshun opened this issue May 26, 2023 · 2 comments
Open

DNA is not implemented yet? #521

songyongshun opened this issue May 26, 2023 · 2 comments
Labels

Comments

@songyongshun
Copy link

It seems the martinze2 cannot recognize the DNA.
The error is
WARNING - unknown-residue - Can't add bonds based on atom names for residue A-DT1 because 'Residue DT is not known to force field charmm'. Falling back to distance criteria.

@pckroon
Copy link
Member

pckroon commented Jun 6, 2023

@fgrunewald I seem to remember we used to have DNA but it got removed a while ago for a good reason. I just can't remember the reason. Do you remember?

Either way, @songyongshun you can try to add the required files on your local installation. Do be sure to check the resulting ITP very carefully. One of the things I know are an issue is with the elastic network. It's (currently) not possible to use different parameters for the elastic network for different parts of your system (e.g. proteins and dna), while the respective models do need different parameters....

@fgrunewald
Copy link
Member

@pckroon the issue is / was that Martini2 DNA itp files are somewhat different in terms of EN depending on the martinize-dna version. To make things more complicated not all EN options are currently stable - reason unknown. Finally, repair graph takes ages to repair nucleobases unless they already contain all hydrogen-atoms and proper labels, but that's more an implementation detail. I even have all files lying around somewhere. If I find the time I'll add them to the repro again using the conventions we settled on for polyply.

Polyply has a working Martini2 DNA implementation for single stranded and we're currently expanding it to dsDNA. @songyongshun if you need M2 DNA especially for large structures feel free to open a discussion here.

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
Projects
None yet
Development

No branches or pull requests

3 participants