01
Identify donor and acceptor
Find a donor atom bonded to hydrogen and a compatible acceptor atom with available electron density.
Computational Analyses · Basic to Advanced | AMRA-LAB
Hydrogen-Bond Analysis
What hydrogen bonds are, how geometrical criteria detect them, why they matter in proteins and protein–ligand complexes, how to interpret their number, occupancy, distance, angle, persistence, and when H-bond analysis can mislead you.
02
Apply geometry
Measure donor–acceptor distance and hydrogen-bond angle for every trajectory frame.
03
Interpret persistence
Study counts, occupancy, distributions, lifetimes, residue identity, and structural context together.
1. What is a hydrogen bond?
A hydrogen bond is a directional non-covalent interaction involving a hydrogen donor and an electronegative acceptor.
Simple definition
A typical hydrogen bond is written as D–H···A. The donor atom D is covalently bonded to hydrogen H, while the acceptor A provides electron density that interacts with the polarized hydrogen.
In proteins and biomolecular simulations, common donors include O–H and N–H groups. Common acceptors include suitable oxygen and nitrogen atoms, but the exact classification depends on protonation state, chemical environment, force-field topology, and software rules.
বাংলায় সহজ করে: Donor atom-এর সঙ্গে যুক্ত hydrogen যখন কাছাকাছি acceptor atom-এর সঙ্গে নির্দিষ্ট distance ও angle-এ interaction করে, তখন তাকে hydrogen bond বলা হয়।
Important: Hydrogen bonds are not identified from distance alone. Their direction and chemical donor/acceptor identity also matter.
2. Why is H-bond analysis important?
Protein architecture
Backbone hydrogen bonds support helices, sheets, turns, and many local structural motifs.
বাংলা: Secondary structure বজায় রাখতে সাহায্য করে।
Ligand recognition
Direct and water-mediated H-bonds can help orient a ligand and contribute to binding specificity.
বাংলা: Ligand pocket-এ সঠিক orientation ধরে রাখতে সাহায্য করতে পারে।
Solvent network
Water can compete with, replace, or bridge protein–ligand and intraprotein hydrogen bonds.
বাংলা: Water direct bond ভাঙতে বা bridge তৈরি করতে পারে।
What this analysis can reveal
- Whether key donor–acceptor pairs remain geometrically compatible during simulation.
- Which residues repeatedly participate in direct or water-mediated interactions.
- Whether a binding pose maintains expected polar contacts.
- Whether mutation, ligand binding, pH-related protonation, or conformational change modifies the H-bond network.
- Whether secondary-structure changes coincide with loss or gain of backbone hydrogen bonds.
বাংলা: কোন residue বা ligand atom বারবার H-bond তৈরি করছে এবং interaction সময়ের সঙ্গে টিকে আছে কি না তা বোঝা যায়।
3. How is a hydrogen bond detected?
Trajectory software usually applies a geometric rule to each possible donor–hydrogen–acceptor triplet.
Identify donor atoms
Select chemically valid donor atoms that carry bonded hydrogens.
বাংলা: প্রথমে hydrogen-যুক্ত donor atom চিহ্নিত করা হয়।
Identify acceptor atoms
Select atoms capable of accepting a hydrogen bond according to chemical and software definitions.
বাংলা: এরপর সম্ভাব্য acceptor atom নেওয়া হয়।
Measure donor–acceptor distance
Calculate the distance between D and A. If it exceeds the chosen cutoff, that frame is not counted as bonded.
বাংলা: Donor ও acceptor বেশি দূরে থাকলে H-bond ধরা হবে না।
Measure direction or angle
Evaluate the D–H···A geometry. More linear arrangements generally represent better directional alignment.
বাংলা: শুধু distance নয়, hydrogen-এর direction-ও সঠিক হতে হবে।
Assign bond state per frame
Store 1 when the triplet satisfies all selected criteria and 0 when it does not.
বাংলা: প্রতিটি frame-এ bond আছে বা নেই হিসেবে record করা হয়।
Summarize the trajectory
Calculate total count, occupancy, donor–acceptor distribution, angle distribution, persistence, and pair identity.
বাংলা: সব frame মিলিয়ে count, occupancy ও distribution বের করা হয়।
Typical GROMACS convention: donor–acceptor distance ≤ 0.35 nm and hydrogen–donor–acceptor angle ≤ 30°. Always report the exact definition used because angle conventions differ across programs.
4. Animated hydrogen-bond geometry
Move the donor–acceptor distance and angle to see when a candidate interaction satisfies the selected geometric rule.
Donor–hydrogen–acceptor visualizer
The bond turns green only when both criteria are satisfied.
D–H···A geometry
PassDistance criterion
PassAngle criterion
BondedFrame assignment
বাংলা: Distance cutoff-এর ভেতরে এবং angle cutoff-এর ভেতরে থাকলেই এই frame-এ H-bond গণনা করা হবে।
How occupancy is calculated
Every frame contributes either bonded = 1 or not bonded = 0.
Occupancy equation
Occupancy (%) = (bonded frames ÷ analyzed frames) × 100
Example: a bond present in 750 of 1000 frames has 75% occupancy.
75.0%Occupancy
750Bonded frames
250Unbonded frames
Frequently present interaction
This interaction appears in most analyzed frames, but its biological importance still depends on residue identity, geometry, alternatives, water competition, and replicate consistency.
বাংলা: Occupancy দেখায় মোট analyzed frame-এর কত শতাংশে নির্দিষ্ট H-bond ছিল। বেশি occupancy মানেই binding energy বেশি—এমন নয়।
5. Which hydrogen bonds should be analyzed?
Intraprotein hydrogen bonds
Measure H-bonds within the protein to study secondary-structure support, domain contacts, loop closure, salt-bridge-associated polar networks, and mutation-induced rearrangement.
বাংলা: Protein-এর নিজের ভেতরের H-bond network বোঝায়।
Protein–ligand hydrogen bonds
Measure direct contacts between protein donors/acceptors and ligand acceptors/donors. Report pair identity and occupancy, not only the total number.
বাংলা: Ligand ও pocket residue-এর direct polar contact দেখা হয়।
Protein–water and ligand–water H-bonds
Useful for solvation, exposed polar groups, desolvation, hydration changes, and competition between water and direct binding contacts.
বাংলা: Water কোন polar group-কে stabilize করছে তা বোঝা যায়।
Water-mediated bridges
A bridging water simultaneously connects two partners through two hydrogen bonds. Direct pair analysis alone does not automatically identify complete bridge networks.
বাংলা: একটি water protein ও ligand-এর মাঝে bridge হিসেবে কাজ করতে পারে।
Backbone–backbone H-bonds
These are central to helices and beta sheets. Combine with secondary-structure analysis because simple counts do not uniquely define structural motifs.
বাংলা: Helix ও sheet বজায় রাখার H-bond বোঝায়।
Side-chain and catalytic H-bond networks
Analyze active-site residues, proton relays, catalytic triads, conserved waters, and residues controlling substrate positioning.
বাংলা: Catalytic mechanism বা active-site network বোঝার জন্য প্রয়োজন।
6. How should H-bond plots be interpreted?
Stable average count
The total number fluctuates around a similar range. This suggests a broadly persistent network, but individual bonds may still exchange rapidly.
Decreasing count
May indicate structural expansion, loss of contacts, solvent competition, ligand movement, or a conformational transition. It is not automatically evidence of instability.
Increasing count
May reflect compaction, formation of new polar contacts, water-network reorganization, or ligand adaptation. More bonds are not always better.
Strong switching
Frequent formation and breakage may indicate flexible contacts, alternate donor–acceptor partners, solvent exchange, or multiple conformational states.
Total count limitation: Two systems can have the same average H-bond count but completely different residue pairs, occupancies, lifetimes, and structural functions.
7. Occupancy, lifetime, and persistence
| Metric | Meaning | Useful for | Major limitation |
|---|---|---|---|
| Instantaneous count | Number of H-bonds in each frame. | Overall time evolution. | Does not identify which bonds are present. |
| Pair occupancy | Percentage of frames satisfying the criterion. | Persistent residue or ligand contacts. | Depends on cutoff, sampling, and intermittent rebinding. |
| Continuous lifetime | Duration before a bond first breaks. | Uninterrupted persistence. | Very sensitive to brief geometric interruptions. |
| Intermittent lifetime | Allows temporary breaking and rebinding under a defined model. | Dynamic contact memory. | Requires explicit methodological definition. |
| Distance distribution | Distribution of donor–acceptor separation. | Geometry quality and alternate states. | Distance alone is insufficient. |
| Angle distribution | Distribution of directional alignment. | Bond geometry and orientation. | Angle convention varies among tools. |
বাংলায় সহজ করে: Count মোট bond সংখ্যা দেখায়, occupancy নির্দিষ্ট bond কত সময় ছিল তা দেখায়, আর lifetime bond কতক্ষণ একটানা টিকে ছিল তা বোঝায়।
8. When is H-bond analysis necessary?
Especially useful when
- Testing whether a docked polar interaction survives MD.
- Comparing apo and ligand-bound protein states.
- Studying mutation effects on local networks.
- Analyzing secondary-structure changes.
- Investigating catalytic residue arrangements.
- Examining solvent exposure and water bridges.
- Comparing multiple ligands using identical definitions and sampling.
Not sufficient by itself when
- Claiming binding strength or binding free energy.
- Proving that a ligand remains bound.
- Determining global protein stability.
- Explaining hydrophobic or ionic interactions.
- Assessing entropic contributions.
- Comparing simulations with different protonation or atom naming without validation.
- Interpreting metal coordination as a normal H-bond.
Combine with: ligand RMSD, contact analysis, minimum distances, salt bridges, SASA, secondary structure, water occupancy, free-energy methods, visual inspection, and replicate simulations.
9. Common mistakes and misleading conclusions
“More hydrogen bonds always mean stronger binding”
Incorrect. Binding depends on desolvation, geometry, competing interactions, entropy, electrostatics, conformational strain, and the surrounding environment.
বাংলা: বেশি H-bond মানেই বেশি binding affinity নয়।
Ignoring protonation and tautomer states
A wrong protonation or tautomer assignment changes donor and acceptor identity and may create or remove apparent hydrogen bonds.
বাংলা: ভুল protonation হলে donor/acceptor definition ভুল হবে।
Reporting only the average count
Average count hides pair exchange. Report important donor–acceptor identities, occupancy, distributions, and representative structures.
বাংলা: শুধু average number দিলে কোন bond টিকে ছিল তা জানা যায় না।
Using incompatible groups
For GROMACS group-based analysis, selected groups must meet the tool requirements; identical or non-overlapping groups are expected.
বাংলা: Group selection overlap করলে analysis ভুল বা invalid হতে পারে।
Confusing angle conventions
Some programs use deviation from linearity; others report the D–H···A angle directly. A cutoff of 30° in one convention may correspond to 150° in another representation.
বাংলা: Software angle কীভাবে define করে তা না জানলে cutoff ভুল বুঝবেন।
Ignoring periodic boundary conditions
Broken molecules or complexes across the box can generate false distances or missed interactions. Process and visually inspect the trajectory.
বাংলা: PBC correction না করলে bond count ভুল হতে পারে।
Treating all hydrogen bonds as equivalent
A catalytic H-bond, a surface water contact, and a transient terminal contact do not have equal structural significance.
বাংলা: সব H-bond-এর biological importance একই নয়।
10. GROMACS hydrogen-bond analysis
Commands vary between the current selection-based tool and the legacy implementation. Confirm the syntax supported by your installed GROMACS version.
Step 1 — Correct the trajectory
gmx trjconv -s md.tpr -f md.xtc -o md_whole.xtc -pbc mol -center
Choose a suitable centring group and output group. Verify membrane proteins and multimers visually because generic processing may be inappropriate.
Step 2 — Current selection-based intraprotein analysis
gmx hbond -s md.tpr -f md_whole.xtc -r 'group "Protein"' -t 'group "Protein"' -num hbonds_protein.xvg -dist hbdist_protein.xvg -ang hbang_protein.xvg
The current tool supports reference and target selections. The selections must be identical or non-overlapping.
Step 3 — Protein–ligand H-bonds
gmx hbond -s md.tpr -f md_whole.xtc -r 'group "Protein"' -t 'resname LIG' -num hbonds_protein_ligand.xvg -dist hbdist_protein_ligand.xvg -ang hbang_protein_ligand.xvg
Replace LIG with the actual ligand residue name. Validate donor/acceptor assignment from the topology and chemistry.
Step 4 — Explicit criteria
gmx hbond -s md.tpr -f md_whole.xtc -r 'group "Protein"' -t 'resname LIG' -hbr 0.35 -hba 30 -num hbonds_PL.xvg
Here -hbr controls the H-bond distance cutoff and -hba controls the angular cutoff in the current GROMACS definition.
Legacy installation example
gmx hbond-legacy -s md.tpr -f md_whole.xtc -n index.ndx -num hbnum.xvg -dist hbdist.xvg -ang hbang.xvg
Legacy output can include autocorrelation and lifetime-related analyses. Do not mix current and legacy syntax without checking your version.
Version warning: The H-bond command interface has changed across GROMACS releases. Run
gmx hbond -h and record your software version before analysis.11. What should be reported?
Software name and exact version.
Trajectory processing and PBC correction.
Reference and target atom selections.
Donor and acceptor definitions.
Distance cutoff and measured atom pair.
Angle definition and cutoff.
Time range and saved-frame interval.
Whether direct, water-mediated, or both were analyzed.
Occupancy definition and continuity model.
Replicates, uncertainty, and pair-selection rules.
Protonation and tautomer treatment.
Companion analyses supporting interpretation.
Example reporting sentence
“Protein–ligand hydrogen bonds were analyzed over the production trajectory after periodic-boundary correction using a donor–acceptor distance cutoff of 0.35 nm and the stated angular criterion. The number of hydrogen bonds, donor–acceptor distance distributions, angle distributions, and pair occupancies were evaluated.”
বাংলা: Method section-এ group, cutoff, angle convention, occupancy definition, time range ও software version স্পষ্টভাবে লিখুন।
12. H-bond analysis FAQ
Is a high H-bond occupancy always good?
No. It means the chosen geometry is frequently satisfied. Importance depends on location, chemistry, alternatives, desolvation, function, and reproducibility.
Can one atom form multiple hydrogen bonds?
Depending on chemistry and geometry, a donor or acceptor can participate in alternative or simultaneous interactions. Interpret unrealistic over-coordination carefully.
Why does the H-bond count fluctuate rapidly?
Hydrogen bonds are dynamic and geometric cutoffs convert continuous motion into a binary state. Small movements near a cutoff can cause frequent switching.
Should water be included?
Include water when hydration, competition, or bridging matters. Exclude it when you specifically want direct protein–ligand contacts, but state that choice.
Can H-bond analysis prove binding?
No. A ligand may maintain H-bonds while drifting, or remain bound through hydrophobic and ionic interactions with few hydrogen bonds.
Why do different programs give different counts?
They may use different donor/acceptor rules, atom naming, distance definitions, angle conventions, cutoffs, protonation states, and PBC handling.
Should I compare raw counts between different-sized proteins?
Not directly. Larger systems have more possible donors and acceptors. Normalize or focus on equivalent regions and matched atom selections.
Final interpretation rule
Hydrogen-bond analysis describes how often chemically and geometrically defined donor–hydrogen–acceptor arrangements occur. It does not independently measure binding free energy, global stability, or biological activity.
বাংলায় মূল কথা: H-bond count-এর সঙ্গে pair identity, occupancy, distance, angle, water competition, protonation, structure visualization এবং অন্যান্য interaction analysis মিলিয়ে conclusion দিতে হবে।