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

Update instructions for hydrogenation at pH 7 #28

Open
wants to merge 1 commit into
base: master
Choose a base branch
from

Conversation

naveenarun
Copy link

I noticed that the command obabel mymols.sdf -O outputfile.smi -p gives the wrong output, and should be obabel mymols.sdf -O outputfile.smi -p 7.4 instead.

@ghutchis
Copy link
Member

Really? The default is 7.4. Can you give me an example on how the "wrong output" occurs?

@naveenarun
Copy link
Author

Sure, here is the relevant info to reproduce. I used 5kh7.pdb (attached here as 5kh7.txt since Github restricts file types) and used the following two commands with their corresponding outputs:

obabel 5kh7.pdb -O 5kh7_nop.mol2 -p: gives 5kh7_nop.txt

obabel 5kh7.pdb -O 5kh7_withp.mol2 -p 7.4: gives 5kh7_withp.txt

The main difference is that the carboxyl group on 6T7 is protonated when -p is used but not when -p 7.4 is used. This can be verified using grep -c 'UNL16' *.mol2 in which an extra hydrogen is present in 5kh7_nop.mol2.

I tried doing a scan of number of hydrogens vs. pH using grep -c ' H ' *.mol2, it seems that -p corresponds to a pH of 0 in this case.

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

Successfully merging this pull request may close these issues.

2 participants