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

Ligand with Charge #53

Open
pealwl opened this issue Aug 13, 2024 · 1 comment
Open

Ligand with Charge #53

pealwl opened this issue Aug 13, 2024 · 1 comment

Comments

@pealwl
Copy link

pealwl commented Aug 13, 2024

I would like to know what I should do if my ligand is charged in my protein-ligand complex. For example, the ligand has a -COOH group in it, and it is deprotonated in my ligand and becomes -COO⁻. If I run Uni-GBSA normally, it will give the following error:

Running: /home/username/anaconda3/envs/gbsa/bin/atomtype -i ANTECHAMBER_AC.AC0 -o ANTECHAMBER_AC.AC -p gaff
Info: Total number of electrons: 71; net charge: 0
Info: The number of electrons is odd (71).
Please check the total charge (-nc flag) and spin multiplicity (-m flag).

which means that it assumes that the molecule is neutral and the molecule is a radical (-COO·). In this case, how can I specify the negative charge in Uni-GBSA?

Any help would be appreciated.

@Aunity
Copy link
Collaborator

Aunity commented Oct 16, 2024

You can try the command unigbsa-validate to fix this problem. But please double check the output file.

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