hbondauto
HBondAuto is a program developed to calculate and automate the hydrogen bond analysis in Molecular Dynamics Simulations.
Science Score: 44.0%
This score indicates how likely this project is to be science-related based on various indicators:
-
✓CITATION.cff file
Found CITATION.cff file -
✓codemeta.json file
Found codemeta.json file -
✓.zenodo.json file
Found .zenodo.json file -
○DOI references
-
○Academic publication links
-
○Academic email domains
-
○Institutional organization owner
-
○JOSS paper metadata
-
○Scientific vocabulary similarity
Low similarity (11.2%) to scientific vocabulary
Keywords
Repository
HBondAuto is a program developed to calculate and automate the hydrogen bond analysis in Molecular Dynamics Simulations.
Basic Info
Statistics
- Stars: 2
- Watchers: 1
- Forks: 1
- Open Issues: 1
- Releases: 1
Topics
Metadata Files
README.md
🔴❗UNDER DEVELOPMENT❗🔴
Hydrogen bond analysis automation (HBondAuto)
HBondAuto is a program developed to calculate and automate the hydrogen bond analysis in Molecular Dynamics Simulations.

Up until now the program is able to: * Calculate the hydrogen bonds in the system through simulation time * Calculate the hydrogen bonds in specific residues of interest: Up until three. * Plot the hydrogen bonds in specific residues of interest
In future updates: * Plot all the hydrogen bonds
Instalation
Download the code and unzip it on the desirable directory. To prepare the environment use the following command:
conda env create -f environment.yml
Be aware to uncomment the sections on the environment.yml file depending on which OS you are using.
How to use
To calculate the hydrogen bonds this program needs the protein file (pdb), with only the protein, and the trajectory file (xtc) with only the protein atoms (please use gmx makendx -f mdprotein-only.pdb -o index.ndx to create an index and gmx trjconv -f md.xtc -o md_protein-only.xtc -n index.ndx to get only the protein trajectory). Do not use the noPBC file in this analysis.
Modify line 47 to add the residues of interest: ```
Modify here the residues of interest
residues = {'donorresidues':[["PHE"],["GLY"],["GLN"]], 'donorindexes':[["28"],["32"],["33"]], 'acceptorresidues':[["ALA"],["ASN"],["ASN"]], 'acceptorindexes':[["23"],["30"],["30"]]} ```
Activate the environment using:
conda activate hbondanalysis
The program uses a shell script to automate the calculation of data in multiple folders. To use this feature go to the root directory and on the terminal use:
bash \hbondauto.sh
If you are interested only on running one folder, you may just add your files to the HBondAuto directory and use:
python analysis.py
Observations
- This script was developed following some steps on the MDAnalysis tutorial for calculating hydrogen bonds. The tutorial can be found here.
Authorship
- Author: Brenda Ferrari (brendaferrari)
Social preview original photo by Brenda Ferrari (brendaferrari)
Owner
- Name: Brenda Ferrari
- Login: brendaferrari
- Kind: user
- Repositories: 6
- Profile: https://github.com/brendaferrari
Chemistry PhD student - Python and Machine Learning student
Citation (CITATION.cff)
cff-version: 0.0.1 message: "If you use this software, please cite it as below." authors: - family-names: "Ferrari" given-names: "Brenda" orcid: "https://orcid.org/0000-0001-5609-0947" title: "HBondAuto" version: 0.0.1 doi: 10.5281/zenodo.7517217 date-released: 2023-01-09 url: "https://github.com/brendaferrari/HBondAuto"
GitHub Events
Total
Last Year
Dependencies
- matplotlib 3.5.2.*
- mdanalysis 2.2.0.*
- numpy 1.23.1.*
- pandas 1.4.2.*
- python 3.8.*