Ecosyste.ms: Awesome

An open API service indexing awesome lists of open source software.

Awesome Lists | Featured Topics | Projects

https://github.com/rasbt/hbind

Calculates hydrogen-bond interaction tables for protein-small molecule complexes, based on protein PDB and protonated ligand MOL2 structure input. Raschka et al. (2018) J. Computer-Aided Molec. Design
https://github.com/rasbt/hbind

bioinformatics computational-biology data-science hydrogen-bonds protein-ligand-interfaces

Last synced: 25 days ago
JSON representation

Host: GitHub
URL: https://github.com/rasbt/hbind
Owner: rasbt
License: apache-2.0
Created: 2017-10-20T22:32:48.000Z (about 7 years ago)
Default Branch: master
Last Pushed: 2019-11-06T14:57:07.000Z (about 5 years ago)
Last Synced: 2024-10-04T12:57:26.109Z (about 1 month ago)
Topics: bioinformatics, computational-biology, data-science, hydrogen-bonds, protein-ligand-interfaces
Language: C
Homepage: http://psa-lab.github.io/Hbind
Size: 1.1 MB
Stars: 26
Watchers: 5
Forks: 14
Open Issues: 3
Metadata Files:
- Readme: README.md
- License: LICENSE.txt

Awesome Lists containing this project

README

# Hbind -- Identifying hydrogen bonds by donor/acceptor chemistry and geometric constraints

![](docs/sources/images/hbind-logo.png)

Software to rigorously define intermolecular H-bonds by donor/acceptor chemistry and geometric constraints, which was developed, used, and described in detail in

- Raschka, Sebastian, Alex Wolf, Joseph Bemister‐Buffington, and Leslie A. Kuhn (2018) “Protein-ligand interfaces are polarized: discovery of a strong trend for intermolecular hydrogen bonds to favor donors on the protein side with implications for predicting and designing ligand complexes” Journal of Computer-Aided Molecular Design. *Journal of Computer-Aided Molecular Design* (accepted manuscript) [[biorxiv preprint](https://www.biorxiv.org/content/early/2018/02/05/260612)]

Code Repository: https://github.com/psa-lab/Hbind
Documentation: [https://psa-lab.github.io/Hbind](https://psa-lab.github.io/Hbind)

---

![](docs/sources/images/1kpf_interact.png)

[Interactions between an PKCI-substrate analog ([1KPF](https://www.rcsb.org/pdb/explore.do?structureId=1kpf)) with its ligand (adenosine monophsophate) via [HbindViz](https://github.com/rasbt/HbindViz) and [PyMOL](https://pymol.org); hydrogen atoms not shown]

---

## Installation

Installing Hbind requires [Perl](https://www.perl.org) and the [GCC](https://gcc.gnu.org) compiler, both of which come pre-installed with most Unix- and Linux-based operating systems.

To install Hbind, simply download this repository, unpack it, and navigate into the main Hbind folder. Then, execute the following command in your terminal:

perl install_hbind.pl

Upon successful installation, the Hbind software will be ready to use from the `bin/` subdirectory. To show the help and usage menu execute execute the following command in your terminal:

./bin/hbind -h

```
HBIND Version: 1.0.0

Documentation: http://psa-lab.github.io/Hbind
Raschka, Wolf, Bemister-Buffington, Kuhn (2018)
Protein Structure and Analysis Lab, MSU (http://kuhnlab.bmb.msu.edu)

USAGE:
-p STRING Path to protein PDB file
-l STRING Path to ligand mol2 file (in docked conformation)
-s Include saltbridges in the output
-t Include a summary table in the output
```

Please see the "Usage" section below for more information on how to use Hbind.

## Usage

To execute Hbind, please navigate to its main directory from where you installed it and run the following command to identify the hydrogen bonds in a protein ligand interface:

./bin/hbind -p ./example_files/1KPF.pdb \
-l ./example_files/1KPF_AMP.mol2

Here, 1KPF.pdb is an excerpt from the PDB structure [1KPF](https://www.rcsb.org/pdb/explore.do?structureId=1kpf) (PKCI-substrate analog) with its ligand `1KPF_AMP.mol2` (adenosine monophsophate). The structures are provided in Hbind's subdirectory "[example_files/](./example_files)".

***Before you run Hbind, please make sure that you provide the protein structure in PDB format and the ligand structure in MOL2 format, and please ensure that the ligand's protonation state (hydrogen assignment) is correct.***

To protonate a ligand, tools such as [OpenEye MolCharge](https://docs.eyesopen.com/quacpac/molchargeusage.html) can be used. However, it is recommended to protonate a ligand in presence of its protein environment, for example, by using tools such as [What If](http://swift.cmbi.ru.nl/whatif/) or [Yasara OptHyd](http://www.yasara.org/features.htm).

After running the command from the beginning of this section, the following output should be displayed:

```
HBIND Version: 1.0.0

Documentation: http://psa-lab.github.io/Hbind
Raschka, Wolf, Bemister-Buffington, Kuhn (2018)
Protein Structure and Analysis Lab, MSU (http://kuhnlab.bmb.msu.edu)

Ligand file: /Users/sebastian/code/hbind/example_files/1KPF_AMP.mol2
Protein file: /Users/sebastian/code/hbind/example_files/1KPF.pdb

++++++++++++++++++++++++++++++++ HBind Interaction Table ++++++++++++++++++++++++++++++++
# | Ligand Atom -- Protein Atom | Bond D-H-A Lig.-Prot.
# | # TYPE -- RES CH_ID RES_NO A_TYPE | DIST. ANGLE INTERACTION
| hbond 1 2 O.2 -- ASN A 99 ND2 2.911 174.0 Acceptor - Donor
| hbond 2 3 O.3 -- SER A 107 N 2.802 166.8 Acceptor - Donor
| hbond 3 3 O.3 -- VAL A 108 N 3.031 178.8 Acceptor - Donor
| hbond 4 4 O.3 -- HIS A 112 NE2 2.593 163.1 Acceptor - Donor
| hbond 5 4 O.3 -- HIS A 114 NE2 2.581 140.4 Acceptor - Donor
| hbond 6 10 O.3 -- ASP A 43 OD2 2.714 167.1 Donor - Acceptor
| hbond 7 12 O.3 -- ASP A 43 OD1 2.607 154.1 Donor - Acceptor
| hbond 8 22 N.ar -- ILE A 44 N 3.156 154.9 Acceptor - Donor
```

To include salt bridges in the listing, you can use the `-s` option as shown below:

./bin/hbind -p ./example_files/1KPF.pdb \
-l ./example_files/1KPF_AMP.mol2 -s

```
HBIND Version: 1.0.0

Documentation: http://psa-lab.github.io/Hbind
Raschka, Wolf, Bemister-Buffington, Kuhn (2018)
Protein Structure and Analysis Lab, MSU (http://kuhnlab.bmb.msu.edu)

Ligand file: /Users/sebastian/code/hbind/example_files/1KPF_AMP.mol2
Protein file: /Users/sebastian/code/hbind/example_files/1KPF.pdb

For additional information, i.e. a summary talbe that includes the number of hydrophobic contacts or [SLIDE](http://www.kuhnlab.bmb.msu.edu/software/slide/) docking scores, provide the `-t` flag:

./bin/hbind -p ./example_files/1KPF.pdb \
-l ./example_files/1KPF_AMP.mol2 -s -t

```
HBIND Version: 1.0.0

Documentation: http://psa-lab.github.io/Hbind
Raschka, Wolf, Bemister-Buffington, Kuhn (2018)
Protein Structure and Analysis Lab, MSU (http://kuhnlab.bmb.msu.edu)

Ligand file: /Users/sebastian/code/hbind/example_files/1KPF_AMP.mol2
Protein file: /Users/sebastian/code/hbind/example_files/1KPF.pdb

## Hydrogen bond rules

Hbind's rules for identifying hydrogen bonds protein-ligand interfaces are based on the criteria by Ippolito et al. [1] and McDonald and Thornton [2].

#### Hydrogen bond criteria:

- Hydrogen to acceptor distance: 1.5-2.5 Å
- Donor to acceptor distance: 2.4-3.5 Å
- Donor-H-acceptor angle (θ): 120-180°
- Pre-acceptor–acceptor–H angle (φ): 90-180°

(Hydrogen bonds must meet all of these criteria.)

#### Criteria for metal interactions:

- Lone pair atom distance to K or Na: 2.0-2.9 Å
- Lone pair atom distance to Ca, Co, Cu, Fe, Mg, Mn, Ni, or Zn: 1.7-2.6 Å

Additional command line options are available to list longer-range salt bridge interactions (up to 4.5 Å) and the number of hydrophobic contancts between protein and ligand atoms.

[1] Ippolito, Joseph A, Richard S Alexander & David W Christianson. 1990. Hydrogen bond stereo-chemistry in protein structure and function. *Journal of Molecular Biology* 215(3). 457–471.
[2] McDonald, Ian & Janet M Thornton. 1994. Atlas of side-chain and main-chain hydrogen bonding. http://www.biochem.ucl.ac.uk/bsm/atlas: Biochemistry and Molecular Biology Department, University College London.

## Hbind tools

### Protein Recognition Index (PRI)

Software for assessing the native-likeness of designed or predicted protein-ligand interfaces, which can be used to guide protein mutagenesis and ligand design.

- [https://psa-lab.github.io/protein-recognition-index/](https://psa-lab.github.io/protein-recognition-index/)

### Hbind Interaction Viz

A program for creating PyMOL H-bond interaction views from Hbind tables.

- [http://psa-lab.github.io/HbindViz](http://psa-lab.github.io/HbindViz)