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

Correct the judgment of hydrogen bonds #370

Open
wants to merge 10 commits into
base: master
Choose a base branch
from

Conversation

Linsastar
Copy link

Errors of adding the essential hydrogen bond interaction have been corrected. In particular, OE was divided into OE1 and OE2.
Besides, OE1, OE2, OD1, OD2 and O were defined as the hydrogen acceptors and the rest hydrogen donors.

Linsastar and others added 2 commits March 2, 2024 11:02
Corrected the error of not distinguishing between hydrogen receptors and donors, distinguishing OE as OE1 and OE2
@a-r-j
Copy link
Owner

a-r-j commented Mar 3, 2024

Thanks for the PR @Linsastar!

Could you clarify why the atoms are different to those specified in the existing Hbond acceptors and donors

@Linsastar
Copy link
Author

Thanks for the PR @Linsastar!

Could you clarify why the atoms are different to those specified in the existing Hbond acceptors and donors
Hi, Arian, My name is Chenyang Gong and I'm an undergraduate student of Dr. Yandong Huang. Comparing the protocols by you and us, there are mainly three differences. First, backbone O (acceptor) and N (donor) atoms are ignored by you. In addition, S atoms on Met, Ser or Thr are not considered by you. But we think these S are likely donors. Finally, we take HSE as the default state of a neutral His where imidazole N-epsilon is protonated and N-delta is deprotonated. Thus, N-epsilon can be a donor and N-delta an acceptor.

@a-r-j
Copy link
Owner

a-r-j commented Mar 7, 2024

Hi Chenyang, thanks for clarifying, that all makes sense to me. Appreciate it!

To cover your points:

  1. Backbone atoms are a clear omission by us here. I think it would be useful however from a code POV to allow users to choose between sc-sc, bb-sc, bb-bb H bonds. What do you think?
  2. Re S atoms, my understanding is that the MET thiother group does not exhibit strong electronegativity. This also seems to be reflected in the current reference for the selection (and some cursory reading). We'll need some concrete references to back up this change,
  3. This change makes sense to me.

@Linsastar
Copy link
Author

Hi Arian, Thanks for your feedback! (1) It makes sense to me that users have three hbond options, sc-sc, bb-sc, bb-bb, according to the locations of donors and acceptors. (2) Indeed, Sulfur involved hydrogen bonds are typically weak when compared with others and thus can be discarded.

Linsastar and others added 2 commits March 16, 2024 15:54
… chain can be identified

Corrected hydrogen bond determination, now hydrogen bonds on the main chain can be identified
Copy link

sonarcloud bot commented Mar 16, 2024

Quality Gate Passed Quality Gate passed

Issues
0 New issues
0 Accepted issues

Measures
0 Security Hotspots
No data about Coverage
No data about Duplication

See analysis details on SonarCloud

if rgroup_df is None:
rgroup_df = G.graph["rgroup_df"]
# if rgroup_df is None:
# rgroup_df = G.graph["rgroup_df"]
Copy link
Owner

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

What is the purpose of not using the supplied rgroup_df?

Copy link
Author

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Yes, you're right. I commented out the code like that just to facilitate debugging and determine the hydrogen bonds on the main chain.

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