/alchemy

Alchemical tools for OpenMM

Primary LanguagePython

Build Status

Alchemical tools for OpenMM

This package contains factories for creating alchemically-modified versions of OpenMM systems.

Installation

The easiest way to install is through the conda package manager through the omnia binstar repository:

conda install -c omnia alchemy

Examples

Create alchemical intermediates for p-xylene binding to T4 lysozyme L99A in implicit solvent.

# Create a reference system for T4 lysozyme L99A
from openmmtools import testsystems
complex = testsystems.LysozymeImplicit()
[reference_system, positions] = [complex.system, complex.positions]
receptor_atoms = range(0,2603) # T4 lysozyme L99A
ligand_atoms = range(2603,2621) # p-xylene

# Create a factory to produce alchemical intermediates.
factory = AbsoluteAlchemicalFactory(reference_system, ligand_atoms=ligand_atoms)

# Get the default protocol for 'denihilating' in complex in implicit solvent.
protocol = factory.defaultComplexProtocolImplicit()

# Create the perturbed systems using this protocol.
systems = factory.createPerturbedSystems(protocol)