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

-nt flag not producing expected result #487

Closed
csbrasnett opened this issue Nov 23, 2022 · 2 comments
Closed

-nt flag not producing expected result #487

csbrasnett opened this issue Nov 23, 2022 · 2 comments

Comments

@csbrasnett
Copy link
Collaborator

I can see this might be related to #472, but I'm not 100% sure so am opening a separate issue.

From the aminoacids.ff file for martini, I would expect the following section regarding protein termini to make the N terminus (BB -BB) a P5 bead, and C terminus (BB +BB) a P6 bead:

[ link ]
[ molmeta ]
neutral_termini true
[ atoms ]
BB {"replace": {"atype": "P5"}}
[ non-edges ]
BB -BB

[ link ]
[ molmeta ]
neutral_termini true
[ features ]
neutral_termini
[ atoms ]
BB {"replace": {"atype": "P6"}}
[ non-edges ]
BB +BB

However, if I run martinize2 with -nt, the result is that both BBs are of type P6.

@pckroon
Copy link
Member

pckroon commented Nov 23, 2022

Termini parameters are defined as modifications (the link hacks you found are from a previous era ;) ). You can find them (for martini3001) in https://github.com/marrink-lab/vermouth-martinize/blob/master/vermouth/data/force_fields/martini3001/modifications.ff. The -nt flag is hardcoded to apply the modifications 'COOH-ter' and 'NH2-ter' (https://github.com/marrink-lab/vermouth-martinize/blob/master/bin/martinize2#L606).
Make sure you're running a version of martinize in which #476 and #477 have been merged.

@csbrasnett
Copy link
Collaborator Author

Thanks for the clarification, can see that now!

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

No branches or pull requests

2 participants