# How to Identify Hydrogen Bonds in Protein-Ligand Complexes: A Step-by-Step Guide Using PyMOL and PLIP

Hydrogen bonds are central to molecular recognition in protein-ligand complexes, and identifying them accurately from a docked structure requires a defined protocol instead of visual inspection alone. This article provides a practical workflow for detecting hydrogen bonds using PyMOL for visualization and PLIP for systematic interaction analysis, including the geometric criteria that define a hydrogen bond, how to interpret the output, and how to avoid common errors in assignment.

The intended reader is a researcher who has completed a molecular docking run and now needs to determine which residues form hydrogen bonds with the ligand, what those interactions mean for binding affinity, and how to report them reproducibly. The protocol assumes a basic working knowledge of protein structure files and command-line usage, but each step is explained with enough detail for a graduate student or laboratory technician to follow without prior experience in interaction analysis.

## The Role of Hydrogen Bonds in Protein-Ligand Binding

Hydrogen bonds arise when a hydrogen atom attached to an electronegative donor atom interacts with another electronegative acceptor atom. In protein-ligand complexes, the common donors are backbone amide nitrogens, side-chain hydroxyls, and amine groups, while the common acceptors are backbone carbonyl oxygens, side-chain carboxylates, and imidazole nitrogens. The strength of an individual hydrogen bond is modest compared to covalent bonds, but the collective network of hydrogen bonds, hydrophobic contacts, and van der Waals interactions determines the overall binding free energy.

For a researcher interpreting a docking result, hydrogen bonds serve two distinct purposes. First, they provide a mechanistic explanation for why a particular ligand conformation is preferred. A ligand that satisfies the hydrogen-bonding potential of the binding pocket is more likely to retain its pose during molecular dynamics simulation. Second, hydrogen bonds guide medicinal chemistry decisions. If a docked ligand forms a hydrogen bond with a catalytic residue, that interaction becomes a design constraint for analog synthesis.

The practical challenge is that docking programs do not always report hydrogen bonds consistently. Some use geometric criteria that are too permissive, while others omit interactions that would be visible in a crystallographic structure. This is why a two-tool approach is recommended. PyMOL provides rapid visual confirmation with user-controlled distance and angle cutoffs, while PLIP applies a standardized geometric model that is suitable for publication and comparison across multiple complexes.

## Geometric Criteria for Hydrogen Bond Assignment

A hydrogen bond in a protein-ligand complex is defined by three geometric parameters: the donor-acceptor distance, the hydrogen-acceptor distance, and the donor-hydrogen-acceptor angle. The distance between the donor and acceptor heavy atoms typically falls between 2.5 and 3.5 angstroms. The distance between the hydrogen atom and the acceptor is shorter, usually between 1.5 and 2.5 angstroms. The angle formed by the donor, hydrogen, and acceptor atoms is ideally close to 180 degrees, with values above 120 degrees generally considered acceptable.

These thresholds are not arbitrary. They reflect the balance between electrostatic attraction and steric repulsion. When the donor-acceptor distance falls below 2.5 angstroms, the atoms are too close and the interaction becomes repulsive. When the distance exceeds 3.5 angstroms, the electrostatic interaction is too weak to be considered a meaningful hydrogen bond. The angle criterion ensures that the hydrogen atom is oriented toward the acceptor lone pair instead of merely being in proximity.

PyMOL allows the user to set these criteria manually through the distance command. The default settings in PyMOL do not automatically identify hydrogen bonds, so the researcher must specify the donor and acceptor selections. PLIP, in contrast, applies its own internal geometric model that has been optimized against known protein-ligand structures. The PLIP model uses a maximum donor-acceptor distance of 3.5 angstroms and a minimum angle of 120 degrees, which aligns with the values commonly reported in structural biology literature.

A common mistake is to rely solely on the distance between heavy atoms without checking the angle. This produces false positives, particularly in flexible side chains where a donor and acceptor may be within 3.5 angstroms but oriented incorrectly. The angle criterion is especially important for hydroxyl groups, which can rotate around the carbon-oxygen bond and present the hydrogen in multiple orientations.

## Preparing the Protein-Ligand Complex for Analysis

Before any hydrogen bond analysis can begin, the protein-ligand complex must be properly prepared. The starting point is usually a docking output file in PDB or SDF format. The quality of the hydrogen bond assignment depends entirely on the quality of the input structure, so the preparation step deserves careful attention.

The first task is to ensure that hydrogen atoms are present in the structure. Many docking programs, including AutoDock Vina, do not include polar hydrogens in the output by default. Without hydrogen atoms, the angle criterion cannot be evaluated, and the analysis is limited to donor-acceptor distances. PyMOL can add hydrogens using the `h_add` command, but the researcher must verify that the protonation states of titratable residues match the physiological pH of the binding assay.

The second task is to check for alternate conformations. If a side chain has multiple occupancies, the hydrogen bond assignment may differ depending on which conformation is selected. The researcher should choose the conformation that is most consistent with the docking pose and note this decision in the methods section.

The third task is to verify the ligand atom naming. PLIP relies on the atom names in the ligand file to assign element types and hybridization states. If the ligand was prepared with a tool that uses nonstandard atom names, PLIP may fail to recognize certain atoms as hydrogen bond donors or acceptors. The ligand should be checked against the SMILES string or the original mol2 file to confirm that all atom names are consistent.

The NCBI provides access to the Protein Data Bank and related structure resources that can be used to obtain reference structures for comparison. When a crystal structure of a related complex is available, the researcher can compare the docked hydrogen bonds against the experimentally observed interactions to validate the docking pose. This comparison is particularly valuable when the docking score is ambiguous or when multiple poses have similar binding energies.

## Loading the Complex in PyMOL

PyMOL is a molecular visualization system that is widely used in structural biology for its scripting capabilities and its ability to render publication-quality images. The first step in the hydrogen bond analysis is to load the protein and ligand into PyMOL as separate objects.

The protein is typically loaded from a PDB file, while the ligand is loaded from a separate file that contains the docked coordinates. The `fetch` command can retrieve a structure directly from the Protein Data Bank if the researcher is working with a known protein, but for a docked complex, the local files should be used to ensure that the coordinates match the docking output.

Once both objects are loaded, the researcher should align the ligand to the protein binding site. This is usually automatic if the docking output already contains the protein coordinates, but if the ligand was docked into a rigid receptor, the ligand coordinates may need to be transformed to match the protein frame. The `align` command in PyMOL can be used to superimpose the docked complex onto the reference structure.

The display should be configured to show the protein as a cartoon or surface representation and the ligand as sticks. The binding site residues should be shown as sticks as well, with the hydrogen bond donors and acceptors colored by element. This visual setup makes it possible to inspect the hydrogen bond geometry before running the automated analysis.

## Using PyMOL Distance Commands for Hydrogen Bond Detection

PyMOL does not have a dedicated hydrogen bond detection command, but the `distance` command can be used to measure the distance between any two atoms or selections. To identify hydrogen bonds, the researcher must create selections for the potential donors and acceptors in both the protein and the ligand.

The `select` command can be used to define the protein donor atoms, which include backbone amide nitrogens and side-chain hydroxyl and amine groups. The protein acceptors include backbone carbonyl oxygens and side-chain carboxylate and imidazole nitrogens. The ligand donors and acceptors must be defined manually based on the ligand structure.

A practical approach is to use the `distance` command with the `mode=2` option, which measures the distance between all pairs of atoms in two selections. The command syntax is:

```
distance hbonds, protein_donors, ligand_acceptors, 3.5, mode=2
```

This command creates a new object called `hbonds` that contains all donor-acceptor pairs within 3.5 angstroms. The researcher can then inspect each pair visually to verify that the angle criterion is satisfied. The `angle` command in PyMOL can measure the donor-hydrogen-acceptor angle for a specific triplet of atoms.

The limitation of this approach is that it requires manual inspection of each candidate pair. For a complex with many potential hydrogen bonds, this can be time-consuming and prone to oversight. The advantage is that the researcher has full control over the criteria and can adjust the distance cutoff to explore weaker interactions.

For a more systematic analysis, the researcher should use PLIP, which automates the geometric evaluation and produces a standardized output that can be compared across multiple complexes.

## Installing and Running PLIP

PLIP, the Protein-Ligand Interaction Profiler, is a command-line tool that detects non-covalent interactions between a protein and a ligand in a structure file. It is available as a standalone tool that can be installed via pip or conda, and it can also be accessed through a web server for single-structure analysis.

The installation process depends on the operating system. On a Linux or macOS system, the recommended approach is to create a virtual environment and install PLIP using pip:

```
python -m venv plip_env
source plip_env/bin/activate
pip install plip
```

On Windows, the same commands work in the Anaconda Prompt or PowerShell with the appropriate activation command. The PLIP documentation provides detailed installation instructions for each platform.

Once PLIP is installed, the analysis is run from the command line. The basic command for a PDB file containing both protein and ligand is:

```
plip -f complex.pdb -o plip_output
```

This command generates a report in the `plip_output` directory that includes a text file with the detected interactions and a PDB file with the interaction features mapped onto the structure. The text report lists each hydrogen bond with the donor residue, the acceptor residue, the distance, and the angle.

PLIP can also analyze a complex where the protein and ligand are in separate files. The `-p` and `-l` options specify the protein and ligand files, respectively. This is useful when the docking output separates the receptor and ligand coordinates.

The output of PLIP is deterministic for a given input structure, which means that the same complex will produce the same interaction list every time. This reproducibility is essential for comparing results across different docking poses or different ligand analogs.

## Interpreting PLIP Hydrogen Bond Output

The PLIP text report uses a standardized format for each hydrogen bond. A typical entry includes the residue number and name of the protein donor or acceptor, the atom name, the distance between the donor and acceptor heavy atoms, and the donor-hydrogen-acceptor angle. The report also indicates whether the interaction is a backbone or side-chain contact.

The first step in interpretation is to check whether the detected hydrogen bonds match the visual inspection in PyMOL. Discrepancies can arise from differences in the geometric criteria or from the treatment of hydrogen atoms. PLIP adds hydrogen atoms to the structure internally, so the protonation state may differ from what the researcher assumed in PyMOL.

The second step is to assess the biological relevance of each hydrogen bond. A hydrogen bond to a backbone carbonyl is often less specific than a hydrogen bond to a side-chain group, because the backbone is present in every residue. A hydrogen bond to a catalytic residue, such as a serine or histidine in an enzyme active site, is more likely to be functionally important.

The third step is to consider the hydrogen bond network as a whole. A single hydrogen bond may be weak, but a network of three or four hydrogen bonds that anchors the ligand in the pocket can be the dominant contributor to binding. The PLIP output should be examined for clusters of interactions that involve the same ligand functional group.

The study of amidine-based BACE-1 inhibitors provides an example of how hydrogen bond analysis informs inhibitor design. In that study, molecular docking identified a compound that formed stable interactions with the catalytic residue ASP32 and key hydrophobic residues TRP115 and PHE108. The hydrogen bond to ASP32 was considered important because it mimics the interaction of the natural substrate with the catalytic dyad. This type of analysis, where a specific hydrogen bond is linked to a functional role, is the goal of the interpretation step.

## Comparing PyMOL and PLIP Results

The two tools serve complementary roles in the hydrogen bond analysis workflow. PyMOL provides interactive visualization that allows the researcher to explore the binding site and understand the three-dimensional context of each interaction. PLIP provides a systematic, reproducible enumeration of interactions that can be used for quantitative comparisons.

The results from the two tools will not always agree. PyMOL with a distance cutoff of 3.5 angstroms will detect more candidate interactions than PLIP, because PyMOL does not apply an angle filter by default. PLIP will exclude interactions where the angle is below 120 degrees, even if the distance is within range. The researcher should treat the PyMOL results as a superset of candidates and the PLIP results as the filtered set of geometrically valid hydrogen bonds.

Disagreements can also arise from the treatment of hydrogen atoms. PyMOL uses the hydrogen atoms present in the structure file, which may be incomplete if the docking program did not add them. PLIP adds hydrogen atoms internally using its own algorithm, which may produce different positions for flexible groups. The researcher should verify that the hydrogen atom positions are chemically reasonable before accepting either result.

When the two tools disagree, the recommended approach is to examine the specific interaction in detail. The researcher should measure the donor-acceptor distance and the donor-hydrogen-acceptor angle manually using the PyMOL `distance` and `angle` commands. If the geometric parameters are borderline, the interaction should be reported as a weak or potential hydrogen bond instead of a definitive one.

## Recording Hydrogen Bond Data for Publication

The hydrogen bond analysis should be recorded in a format that allows other researchers to reproduce the results. The minimum information includes the software versions, the geometric criteria, and the complete list of detected interactions.

The software versions are important because different versions of PyMOL and PLIP may use different algorithms for hydrogen atom placement and geometric evaluation. The researcher should record the version numbers in the methods section of any publication or report.

The geometric criteria should be stated explicitly. If the researcher used a distance cutoff of 3.5 angstroms and an angle cutoff of 120 degrees, this should be reported. If the criteria were adjusted for a specific interaction, the adjustment should be noted.

The interaction list should include the residue number and name for each protein atom involved in a hydrogen bond, the atom name in the ligand, and the measured distance and angle. This information can be presented in a table or as supplementary material. The PLIP output file provides this information in a structured format that can be included directly.

The study of Solanum trilobatum phytochemicals in oral submucous fibrosis provides an example of how hydrogen bond data is used in a broader drug discovery context. In that study, network pharmacology identified six hub genes as potential therapeutic targets, and molecular docking confirmed that a specific phytochemical exhibited significant binding affinities with all six targets. The hydrogen bond analysis would have been used to explain why the phytochemical bound to each target and to guide the selection of the compound for molecular dynamics simulation.

## Common Failure Patterns in Hydrogen Bond Identification

Several recurring errors can compromise the accuracy of hydrogen bond identification. The most common is the failure to add hydrogen atoms before analysis. Without hydrogen atoms, the angle criterion cannot be evaluated, and the researcher may report interactions that are geometrically impossible.

Another common error is the use of an overly permissive distance cutoff. A cutoff of 4.0 angstroms will detect many interactions that are not true hydrogen bonds. The researcher should use the standard cutoff of 3.5 angstroms for donor-acceptor heavy atom distance and verify that the angle is above 120 degrees.

A third error is the misidentification of donor and acceptor atoms. In some functional groups, such as histidine, the donor and acceptor roles depend on the protonation state and tautomeric form. If the researcher assumes the wrong protonation state, the hydrogen bond assignment will be incorrect.

A fourth error is the failure to consider symmetry-related molecules in a crystal structure. If the analysis is performed on a crystal structure with multiple copies of the protein in the asymmetric unit, the hydrogen bonds may differ between copies. The researcher should specify which copy was used for the analysis.

A fifth error is the overinterpretation of hydrogen bonds in docking poses. A docked pose is a prediction, not an experimental observation. The hydrogen bonds in a docked complex should be validated by molecular dynamics simulation or by comparison with a known crystal structure before they are used to guide experimental work.

## Limitations of Computational Hydrogen Bond Detection

The hydrogen bonds identified by PyMOL and PLIP are based on geometric criteria applied to a static structure. This approach has several limitations that the researcher should acknowledge.

First, the structure is a snapshot. In solution, the protein and ligand are in constant motion, and hydrogen bonds are continuously formed and broken. A hydrogen bond that is present in the docked structure may be transient in the actual complex. Molecular dynamics simulation can provide information about the persistence of hydrogen bonds over time, but this is beyond the scope of the basic analysis described here.

Second, the geometric criteria are approximations. The ideal hydrogen bond geometry depends on the electronic structure of the donor and acceptor groups, which is not captured by simple distance and angle cutoffs. Quantum mechanical calculations can provide more accurate hydrogen bond energies, but they are computationally expensive and not practical for routine analysis.

Third, the protonation states are assumed. The researcher must decide which residues are protonated and which are not, based on the pH of the binding assay. If the protonation state is wrong, the hydrogen bond assignment will be wrong. The researcher should use a pKa prediction tool to verify the protonation states of titratable residues in the binding site.

Fourth, the analysis does not account for water-mediated hydrogen bonds. In many protein-ligand complexes, water molecules bridge the protein and ligand through hydrogen bonds. These water-mediated interactions are not detected by the standard analysis, but they can be important for binding. The researcher should inspect the binding site for ordered water molecules in the crystal structure and consider their role in the hydrogen bond network.

The hydrogen/deuterium exchange mass spectrometry method provides an experimental complement to computational hydrogen bond detection. This technique measures the rate at which backbone amide hydrogens exchange with deuterium in solution, and the binding of a ligand can protect certain regions from exchange. The protected regions correspond to the binding interface, which can be compared with the computationally predicted hydrogen bonds.

## Quality Control Checks for the Analysis

Before accepting the hydrogen bond results, the researcher should perform several quality control checks. These checks are designed to catch errors in the input structure and in the analysis parameters.

The first check is to verify that the protein and ligand are in the same coordinate frame. If the ligand was docked into a rigid receptor, the coordinates should match the receptor coordinates exactly. A misalignment will produce false hydrogen bonds or miss true ones.

The second check is to verify the ligand chemistry. The ligand should have the correct protonation state, tautomeric form, and stereochemistry. An incorrect ligand structure will produce incorrect hydrogen bond assignments.

The third check is to verify the residue numbering. If the protein structure has gaps or nonstandard numbering, the hydrogen bond report may be confusing. The researcher should map the residue numbers to the canonical sequence before reporting the results.

The fourth check is to compare the hydrogen bonds with the docking score. A ligand with a favorable docking score should form a reasonable number of hydrogen bonds. If the ligand has no hydrogen bonds but a favorable score, the binding may be driven by hydrophobic interactions, which is possible but should be noted.

The fifth check is to compare the results with known structures. If a crystal structure of a related complex is available, the researcher should compare the hydrogen bonds in the docked complex with those in the crystal structure. This comparison can validate the docking pose and the hydrogen bond assignment.

## Professional Escalation Criteria

The hydrogen bond analysis described here is suitable for routine docking studies and for generating hypotheses about binding determinants. However, there are situations where the researcher should escalate the analysis to a more sophisticated approach.

If the hydrogen bond analysis produces conflicting results between PyMOL and PLIP, and the conflict cannot be resolved by manual inspection, the researcher should consider using a third tool or a more detailed computational method. The conflict may indicate that the docking pose is not reliable, and the researcher should consider re-docking with different parameters.

If the hydrogen bonds are critical for a medicinal chemistry decision, such as the design of a new analog, the researcher should validate the interactions with molecular dynamics simulation. A 100 nanosecond simulation can reveal whether the hydrogen bonds are stable over time or whether they are transient artifacts of the docking pose. The study of Solanum trilobatum phytochemicals used a 100 nanosecond simulation to confirm the stability of the ligand-protein complex, and the BACE-1 inhibitor study used a 200 nanosecond simulation for the same purpose.

If the researcher is working on a drug discovery project that will proceed to experimental validation, the hydrogen bond predictions should be tested experimentally. Techniques such as X-ray crystallography, nuclear magnetic resonance spectroscopy, or hydrogen/deuterium exchange mass spectrometry can provide experimental evidence for the predicted interactions. The hydrogen/deuterium exchange method is particularly useful because it can map the binding interface without requiring a crystal structure.

If the researcher is not confident in the protonation states of the binding site residues, the analysis should be escalated to a tool that can predict pKa values and assign protonation states automatically. The hydrogen bond assignment is only as good as the protonation state assumption.

## At a Glance

| Analysis Step | Tool | Key Parameter | Common Error | Verification Method |
| --- | --- | --- | --- | --- |
| Structure preparation | PyMOL or docking software | Hydrogen atoms present, correct protonation states | Missing hydrogens, wrong tautomers | Check atom count, compare with reference structure |
| Visual inspection | PyMOL | Distance cutoff 3.5 angstroms, angle cutoff 120 degrees | Using distance only, no angle check | Measure angle with `angle` command |
| Systematic detection | PLIP | Internal geometric model, donor-acceptor distance 3.5 angstroms | Ignoring angle filter, misassigned donor/acceptor | Compare with PyMOL candidates |
| Result interpretation | Manual | Biological context, catalytic residues, binding pocket | Overinterpreting backbone contacts | Compare with known crystal structures |
| Validation | Molecular dynamics | Simulation length 100 to 200 nanoseconds | Accepting static pose without dynamics | Check hydrogen bond persistence over trajectory |

## Practical Implementation Steps

The following steps provide a complete workflow for hydrogen bond identification in a docked protein-ligand complex.

Step 1: Prepare the complex. Load the protein and ligand into PyMOL, add hydrogen atoms if they are missing, and verify the protonation states of titratable residues. Save the prepared complex as a new PDB file.

Step 2: Run PLIP. Install PLIP in a virtual environment and run the analysis on the prepared complex. Save the output report and the interaction PDB file.

Step 3: Visualize the PLIP results. Load the PLIP output PDB file into PyMOL and inspect the detected hydrogen bonds. The interaction features are mapped onto the structure, which makes it easy to see which residues are involved.

Step 4: Verify the hydrogen bonds manually. For each hydrogen bond in the PLIP report, measure the donor-acceptor distance and the donor-hydrogen-acceptor angle in PyMOL. Confirm that the geometric parameters satisfy the standard criteria.

Step 5: Identify the biologically relevant interactions. Focus on hydrogen bonds that involve catalytic residues, conserved binding site residues, or functional groups that are known to be important for binding.

Step 6: Record the results. Create a table with the residue number, residue name, atom name, distance, and angle for each hydrogen bond. Note the software versions and the geometric criteria used.

Step 7: Validate with molecular dynamics if needed. If the hydrogen bonds will guide experimental work, run a molecular dynamics simulation to confirm that the interactions are stable over time.

## Records and Measurements

The hydrogen bond analysis produces several types of records that should be preserved for reproducibility. The input structure file, the PLIP output directory, and the PyMOL session file should all be saved with descriptive names that include the date and the complex identifier.

The PLIP output directory contains a text report and a PDB file with interaction features. The text report is the primary record of the detected hydrogen bonds. The PDB file can be loaded into PyMOL for visualization and can be used to generate publication figures.

The PyMOL session file preserves the visualization state, including the selections, the distance measurements, and the display settings. This file is useful for revisiting the analysis and for sharing the results with collaborators.

The measurements that should be recorded for each hydrogen bond are the donor-acceptor distance, the donor-hydrogen-acceptor angle, and the residue and atom identifiers. These measurements should be reported in any publication or presentation that describes the binding interactions.

The software versions should be recorded in the methods section. The PyMOL version can be found in the help menu or by running the `version` command. The PLIP version can be found by running `plip --version` in the command line.

## Common Failure Patterns and How to Avoid Them

The most frequent failure in hydrogen bond identification is the use of a structure that has not been properly prepared. This includes structures with missing hydrogen atoms, incorrect protonation states, or alternate conformations that have not been resolved. The preparation step should never be skipped, even for structures that come from a docking program that claims to add hydrogens.

The second most frequent failure is the reliance on a single tool. PyMOL alone will produce too many candidates because it does not apply an angle filter. PLIP alone will produce a filtered list, but the researcher may not understand why certain interactions were excluded. The two-tool approach provides both the complete candidate list and the filtered result.

The third failure is the overinterpretation of hydrogen bonds in a docking pose. A docked pose is a prediction, and the hydrogen bonds in that pose may not persist in solution. The researcher should always consider the possibility that the docking pose is wrong and that the hydrogen bonds are artifacts.

The fourth failure is the neglect of water-mediated interactions. Many protein-ligand complexes use water molecules as hydrogen bond bridges, and these interactions are not detected by the standard analysis. The researcher should inspect the binding site for ordered water molecules and consider their role.

The fifth failure is the use of inconsistent criteria across different complexes. If the researcher is comparing hydrogen bonds across a series of ligand analogs, the same geometric criteria must be used for all complexes. Changing the distance cutoff between analyses will produce results that cannot be compared.

## Safety and Ethical Context

Hydrogen bond identification is a computational analysis that does not involve hazardous materials or biological samples. The main safety consideration is the proper use of computer resources, particularly when running molecular dynamics simulations that can consume significant CPU time and memory.

The ethical context relates to the responsible use of computational predictions. Hydrogen bonds identified in a docking study are predictions that should be validated experimentally before they are used to guide drug development. Reporting predicted hydrogen bonds as experimentally confirmed interactions is a form of scientific misconduct.

The researcher should also be aware of the limitations of the computational methods and should communicate these limitations in any publication. The geometric criteria used to identify hydrogen bonds are approximations, and the results should be interpreted with appropriate caution.

The training resources provided by the Galaxy Training Network and The Carpentries offer foundational instruction in the computational skills needed for this type of analysis. The Galaxy Training Network provides accessible workflow training that can help researchers understand how to structure reproducible analyses, while The Carpentries lessons cover the shell, Git, and programming fundamentals that are useful for managing analysis pipelines.

## Frequently Asked Questions

### What is the standard distance cutoff for a hydrogen bond in a protein-ligand complex?

The standard distance cutoff for the donor-acceptor heavy atom distance is 3.5 angstroms. This value represents the balance between electrostatic attraction and steric repulsion. Distances below 2.5 angstroms indicate steric clash, while distances above 3.5 angstroms indicate that the electrostatic interaction is too weak to be considered a hydrogen bond. The hydrogen-acceptor distance is shorter, typically between 1.5 and 2.5 angstroms.

### Why does PyMOL not automatically detect hydrogen bonds?

PyMOL is a molecular visualization system, not a dedicated interaction analysis tool. It provides the `distance` command that can measure distances between atoms or selections, but the researcher must define the donor and acceptor selections manually. PyMOL does not apply an angle filter by default, so the researcher must verify the donor-hydrogen-acceptor angle for each candidate interaction. PLIP is recommended for automated detection because it applies a standardized geometric model.

### How does PLIP decide whether an interaction is a hydrogen bond?

PLIP uses a geometric model that evaluates both the donor-acceptor distance and the donor-hydrogen-acceptor angle. The maximum donor-acceptor distance is 3.5 angstroms, and the minimum angle is 120 degrees. PLIP adds hydrogen atoms internally to the structure, so the protonation state may differ from what the researcher assumed. The output report lists each detected hydrogen bond with the residue, atom, distance, and angle.

### What should I do if PyMOL and PLIP give different hydrogen bond results?

Disagreements between PyMOL and PLIP are expected because the two tools use different criteria. PyMOL with a distance cutoff will detect more candidates because it does not apply an angle filter. PLIP will exclude interactions where the angle is below 120 degrees. When the tools disagree, measure the distance and angle manually in PyMOL and decide based on the geometric parameters. If the interaction is borderline, report it as a weak or potential hydrogen bond.

### Can I identify hydrogen bonds without adding hydrogen atoms to the structure?

You can identify potential hydrogen bonds based on donor-acceptor distances alone, but this approach produces false positives because the angle criterion cannot be evaluated. Hydrogen atoms are required for an accurate analysis. Most docking programs do not include polar hydrogens in the output, so you should add them using the `h_add` command in PyMOL or a similar tool before running the analysis.

### How do I know if a hydrogen bond is biologically important?

A hydrogen bond is biologically important if it involves a residue that is known to be critical for binding or catalysis. Hydrogen bonds to catalytic residues, such as the catalytic aspartate in BACE-1, are typically more important than hydrogen bonds to surface residues. The biological importance should be assessed in the context of the binding site and the known structure-activity relationships for the target.

### Should I validate the hydrogen bonds with molecular dynamics simulation?

Molecular dynamics simulation is recommended when the hydrogen bonds will guide experimental work or medicinal chemistry decisions. A simulation of 100 to 200 nanoseconds can reveal whether the hydrogen bonds are stable over time or whether they are transient artifacts of the docking pose. The studies of Solanum trilobatum phytochemicals and BACE-1 inhibitors both used molecular dynamics simulation to confirm the stability of the predicted interactions.

### What are the limitations of computational hydrogen bond detection?

Computational hydrogen bond detection is based on geometric criteria applied to a static structure. The structure is a snapshot, and the hydrogen bonds may be transient in solution. The protonation states are assumed, and water-mediated hydrogen bonds are not detected. The results should be validated experimentally when possible, using techniques such as X-ray crystallography, nuclear magnetic resonance spectroscopy, or hydrogen/deuterium exchange mass spectrometry.

## Related Bioinformatics Guides

- [Gene Set Enrichment Analysis in R: A Practical Tutorial for Interpreting Omics Data](/knowledge/bioinformatics/gene-set-enrichment-analysis-in-r-a-practical-tutorial-for-interpreting-omics-data)
- [Alphafold Protein Ligand Docking: Structural Analysis and Computational Methodologies in Bioinformatics](/knowledge/bioinformatics/alphafold-protein-ligand-docking)
- [Deep Learning for Protein-Ligand Binding Affinity Prediction in Antiviral Drug Design](/knowledge/bioinformatics/deep-learning-protein-ligand-binding-affinity-antiviral-drug-design)
- [Deep Learning in Protein-Ligand Binding Affinity Prediction for Antiviral Drug Design](/knowledge/bioinformatics/deep-learning-protein-ligand-binding-affinity-prediction-antiviral-drug-design)
- [Metabolomics Data Analysis in R: A Practical Workflow](/knowledge/bioinformatics/metabolomics-data-analysis-in-r-a-practical-workflow)

## References and Further Reading

- [NCBI Data Resources](https://www.ncbi.nlm.nih.gov/). National Center for Biotechnology Information.
- [EMBL-EBI Training](https://www.ebi.ac.uk/training). European Bioinformatics Institute.
- [Bioconductor](https://bioconductor.org/). Bioconductor Project.
- [Galaxy Training Network](https://training.galaxyproject.org/). Galaxy Project.
- [nf-core Documentation](https://nf-co.re/docs). nf-core.
- [The Carpentries Lessons](https://carpentries.org/lessons). The Carpentries.
- [Decoding the therapeutic promise of Solanum trilobatum in oral submucous fibrosis: a network pharmacology, docking, and bioactivity study.](https://doi.org/10.1038/s41598-026-46122-8). 2026.
- [<i>In Silico</i> Investigation of Amidine-Based BACE-1 Inhibitors Against Alzheimer's Disease: SAR, Pharmacokinetics, Molecular Docking and Dynamic Simulations.](https://doi.org/10.3390/ph19010005). 2025.
- [Mapping Protein-Ligand Interactions with Proteolytic Fragmentation, Hydrogen/Deuterium Exchange-Mass Spectrometry](https://doi.org/10.1016/bs.mie.2015.08.010). Methods in Enzymology, 2016.

> This article is educational and does not replace validated analysis plans, institutional policy, clinical interpretation, or specialist review.