Abstract:CHARMM is rich in methodology and functionality as one of the first programs addressing problems of molecular dynamics and modeling of biological macromolecules and their partners, e.g., small molecule ligands. When combined with the highly developed CHARMM parameters for proteins, nucleic acids, small molecules, lipids, sugars, and other biologically relevant building blocks, and the versatile CHARMM scripting language, CHARMM has been a trendsetting platform for modeling studies of biological macromolecules. To further enhance the utility of accessing and using CHARMM functionality in increasingly complex workflows associated with modeling biological systems, we introduce pyCHARMM, Python bindings, functions, and modules to complement and extend the extensive set of modeling tools and methods already available in CHARMM. These include access to CHARMM function-generated variables associated with the system (psf), coordinates, velocities and forces, atom selection variables, and force field related parameters. The ability to augment CHARMM forces and energies with energy terms or methods derived from machine learning or other sources, written in Python, CUDA, or OpenCL and expressed as Python callable routines is introduced together with analogous functions callable during dynamics calculations. Integration of Python-based graphical engines for visualization of simulation models and results is also accessible. Loosely coupled parallelism is available for workflows such as free energy calculations, using MBAR/TI approaches or high-throughput multisite λ-dynamics (MSλD) free energy methods, string path optimization calculations, replica exchange, and molecular docking with a new Python-based CDOCKER module. CHARMM accelerated platform kernels through the CHARMM/OpenMM API, CHARMM/DOMDEC, and CHARMM/BLaDE API are also readily integrated into this Python framework. We anticipate that pyCHARMM will be a robust platform for the development of comprehensive and complex workflows utilizing Python and its extensive functionality as well as an optimal platform for users to learn molecular modeling methods and practices within a Python-friendly environment such as Jupyter Notebooks.

CHARMM: A program for macromolecular energy, minimization, and dynamics calculations

CHARMM: The biomolecular simulation program

ORAL: All purpose molecular mechanics simulator and energy minimizer

CHARMM at 45: Enhancements in Accessibility, Functionality, and Speed

Amber Free Energy Tools: Interoperable Software for Free Energy Simulations Using Generalized Quantum Mechanical/molecular Mechanical and Machine Learning Potentials.

Simultaneous Optimization of Biomolecular Energy Functions on Features from Small Molecules and Macromolecules.

The Amber biomolecular simulation programs

CHARMM-GUI Input Generator for NAMD, GROMACS, AMBER, OpenMM, and CHARMM/OpenMM Simulations Using the CHARMM36 Additive Force Field

TAMkin: A Versatile Package for Vibrational Analysis and Chemical Kinetics

Accurate Free Energy Calculation via Multiscale Simulations Driven by Hybrid Machine Learning and Molecular Mechanics Potentials

CHARMM-GUI Free Energy Calculator for Absolute and Relative Ligand Solvation and Binding Free Energy Simulations

CASM Monte Carlo: Calculations of the thermodynamic and kinetic properties of complex multicomponent crystals

pyCHARMM: Embedding CHARMM Functionality in a Python Framework

CHARMM-GUI PDB manipulator for advanced modeling and simulations of proteins containing nonstandard residues.

MOLCAS: a program package for computational chemistry

Automation of the CHARMM General Force Field (CGenFF) II: Assignment of Bonded Parameters and Partial Atomic Charges

GROMACS: Fast, flexible, and free

charmm2gmx: An Automated Method to Port the CHARMM Additive Force Field to GROMACS

Molecular Dynamics Simulations Studies Based GROMACS

Integrated Modeling Program, Applied Chemical Theory (IMPACT)

CHARMM-GUI ligand reader and modeler for CHARMM force field generation of small molecules