I am attempting to find potential Hydrogen bonds between Hydrogen donors and aromatic ring acceptors. I do this by predicting the location of Hydrogens on residues and then calculating how far these Hydrogens are from aromatic rings. If a certain Hydrogen is <7.0 Angstroms from a certain aromatic ring, then I take it under consideration: I form the N-H vector, which is the vector created by the Hydrogen under question and the Nitrogen in the backbone of the residue that the Hydrogen belongs to. I test that this N-H vector is pointing toward the plane of the aromatic ring, and I also test that the point of intersection between the plane of the aromatic and the N-H vector is within 6 Angstroms of the center of the aromatic ring.
If all of these conditions are met, then I consider it a Hydrogen bond between the Hydrogen and the aromatic ring. However, my data must be incorrect, because I am seeing situations where a Hydrogen is < 1.0 Angstrom from the plane of the aromatic. Atoms should not be getting that close to each other.
I thoroughly tested my method by hand using an example situation where my code identified one of the sidechain Hydrogens on an ASN is 0.3 Angstroms from the plane of the aromatic of a TRP. Unfortunately, I could not find any bugs. You can find a PDF of this verification here.
Any suggestions on how my method might be flawed would be greatly appreciated.
No comments:
Post a Comment