hbondauto

HBondAuto is a program developed to calculate and automate the hydrogen bond analysis in Molecular Dynamics Simulations.

https://github.com/brendaferrari/hbondauto

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

hydrogen-bonds molecular-dynamics-simulation
Last synced: 6 months ago · JSON representation ·

Repository

HBondAuto is a program developed to calculate and automate the hydrogen bond analysis in Molecular Dynamics Simulations.

Basic Info
  • Host: GitHub
  • Owner: brendaferrari
  • License: mit
  • Language: Python
  • Default Branch: master
  • Homepage:
  • Size: 37.1 KB
Statistics
  • Stars: 2
  • Watchers: 1
  • Forks: 1
  • Open Issues: 1
  • Releases: 1
Topics
hydrogen-bonds molecular-dynamics-simulation
Created over 3 years ago · Last pushed about 3 years ago
Metadata Files
Readme License Citation

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

Social preview original photo by Brenda Ferrari (brendaferrari)

Owner

  • Name: Brenda Ferrari
  • Login: brendaferrari
  • Kind: user

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

environment.yml conda
  • matplotlib 3.5.2.*
  • mdanalysis 2.2.0.*
  • numpy 1.23.1.*
  • pandas 1.4.2.*
  • python 3.8.*