A Python tutorial for 3DRISM solvation calculations of chemical and biological molecules

Area of Science:

• Computational chemistry
• Biophysics
• Chemical physics

Background:

• The 3-Dimensional Reference Interaction Site Model (3DRISM) is a grid-based solvation model balancing accuracy and efficiency for chemical and biological solutes.
• The EPISOL (Expanded Package of Integral Equation Theory-Based Solvation) package was previously developed to facilitate efficient 3DRISM calculations, offering various closures and compatibility with AMBER and GROMACS.

Purpose of the Study:

• Introduce EPIPY, a new Python-based package designed to streamline 3DRISM calculations.
• Provide a user-friendly interface for running and analyzing 3DRISM simulations.
• Demonstrate the utility of EPIPY through practical examples.

Main Methods:

• Leveraging the EPISOL kernel library within a Python environment.
• Implementing step-by-step tutorials for 3DRISM calculations.
• Utilizing EPIPY for generating water distributions, identifying key solvent coordinates, and calculating solvation free energies.

Main Results:

• Successful development and introduction of EPIPY, a Python package for 3DRISM.
• Demonstration of EPIPY's capability to perform diverse solvation analyses, including water distribution and free energy calculations.
• Availability of source code and tutorials for broader scientific adoption.

Conclusions:

• EPIPY significantly simplifies and enhances the accessibility of 3DRISM solvation modeling.
• The package provides a powerful yet user-friendly tool for researchers in chemistry and biology.
• EPIPY is expected to accelerate research involving solvation analysis of chemical and biological systems.