Replies: 1 comment · 12 replies
-
|
Hi @diaghil, does your PDB file contain explicit hydrogen atoms for both the protein and the ligand? If it does, can you also check if there are bonds defined in the PDB file, and if so is it all bonds or only a subset? If it's a subset, you can add the |
Beta Was this translation helpful? Give feedback.
-
|
@diaghil I haven't received anything (also not in my spam folder), could you double check? |
Beta Was this translation helpful? Give feedback.
-
|
@cbouy Let me know if you received my new email. |
Beta Was this translation helpful? Give feedback.
-
|
I did, will have a look in the coming days |
Beta Was this translation helpful? Give feedback.
-
|
Please let me know. |
Beta Was this translation helpful? Give feedback.
-
|
Okay the problem was indeed with only having partial bonds, but I managed to solve it by calling import prolif as plf
import MDAnalysis as mda
from MDAnalysis.topology.tables import vdwradii
u = mda.Universe("file.pdb")
lig = u.select_atoms("resname UNK")
# I restrict the protein selection to 20 around the ligand but for the trajectory analysis you could
# only keep the `not group ligand` for the selection in case the ligand moves around a lot
prot = u.select_atoms("(byres around 20 group ligand) and not group ligand", ligand=lig)
# guess bonds for the salts
prot.guess_bonds(
vdwradii={
"Na": vdwradii["NA"],
"Mg": vdwradii["MG"],
"Cl": vdwradii["CL"],
})
# run
fp = plf.Fingerprint()
fp.run(u.trajectory, lig, prot)
df = fp.to_dataframe()and you get those interaction types: df.groupby(level="interaction", axis=1).sum()
Tell me if this works for you, |
Beta Was this translation helpful? Give feedback.
Uh oh!
There was an error while loading. Please reload this page.
-
I am trying to see if my fingerprint displays any interactions other than VdWContact. My code is pretty much referencing the documentation directly, and I know the PDB file contains other kinds of interactions.
After removing HOH residues from the PDB input file and making sure none of the atom IDs in the PDB were faulty, I'm still failing to see a full range of results, even without the trajectory file being passed as an argument.
Beta Was this translation helpful? Give feedback.
All reactions