How to extract a list of hydrogen bonds from a PDB?

Question

I am analyzing the 7VF2 structure, from the RCSB-database. I noticed that some residues belonging to one of its chains interact with the other chains that are present within the assembly. My question is if there is a method to extract a list of interacting residues that are present in a PDB.

EDIT:

Just as an example, we have Lys 1616 belonging to ZC3H13 chain that interacts with Glu 621 of the protein virilizer homolog (the first chain)

Mol* 3D visualizer can be used for major information in the 3D View in the link page I provided

Answer 1

You can do that with GROMACS. You probably need to check the manual for this, but roughly these are the steps:

gmx pdb2gmx -f input.pdb -o output.gro -water spc
gmx solvate -cp output.gro -cs spc216.gro -o output_solv.gro -p output.top
gmx pdb2gmx -f output_solv.gro -o output_processed.gro -p output.top -ignh
gmx make_ndx -f output_processed.gro -o index.ndx
gmx hbond -f output_processed.gro -n index.ndx -s output.tpr -num hbond.xvg -hbn hbond.ndx

First, you convert the PDB file into a GROMACS topology file using the pdb2gmx command. Then you can solvate the system with water molecules and add hydrogen atoms. Then you cCreate an index file containing the groups of atoms you want to analyze. Then you can use the hbond command to calculate the hydrogen bonds.

I have not used GROMACS in a while. Hope that provides a starting point. With GROMACS you can even run a simulation and monitor the h-bonds over time.

You actually might need to run an energy minimziation first, as the PDB might not in an optimal state for the h-bonds detection.