PyRod: Tracing Water Molecules in Molecular Dynamics Simulations

Area of Science:

• Computational chemistry and structural biology
• Drug discovery and molecular modeling

Background:

• Water molecules play a crucial role in ligand binding within protein pockets, competing with ligands for interaction sites.
• Understanding water molecule location and thermodynamics is vital for effective drug design and high-throughput virtual screening.

Purpose of the Study:

• To introduce PyRod, a novel Python software for visualizing protein binding pockets and generating pharmacophore features.
• To enable structure-based virtual screening by analyzing water molecule behavior in molecular dynamics simulations.

Main Methods:

• PyRod analyzes water molecule environments in molecular dynamics simulations, differentiating between hydrogen-bonded and hydrophobic/charged/aromatic interactions.
• It generates dynamic molecular interaction fields (dMIFs) for visualization and pharmacophore features for virtual screening.
• The software was tested on five drug targets, with pharmacophores evaluated using DUD-E benchmarking sets.

Main Results:

• PyRod successfully identified hot spots for ligand binding and generated pharmacophore features.
• The software demonstrated strong performance, particularly for HIV1 protease, achieving an early enrichment factor of 54.6.
• Generated pharmacophores are purely protein-based, derived solely from water molecule analysis.

Conclusions:

• PyRod offers a new approach to structure-based screening by utilizing water molecule dynamics.
• The software facilitates virtual screening campaigns when limited or no cocrystallized ligand information is available.
• PyRod provides an accessible tool for drug design, enhancing the identification of potential therapeutic agents.