Skip to content

ATOMScience-org/AMPL

Folders and files

NameName
Last commit message
Last commit date

Latest commit

 

History

2,405 Commits
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 

Develop test Linter Documentation Status codecov

GitHub Release License LinkedIn YouTube

| Install | Docker | Tutorials | Features | Pipeline parameters | Docs |

ATOM Modeling PipeLine (AMPL) for Drug Discovery

An open-source, end-to-end software pipeline for data curation, model building, and molecular property prediction to advance in silico drug discovery.

Created by the Accelerating Therapeutics for Opportunities in Medicine (ATOM) Consortium

AMPL logo

The ATOM Modeling PipeLine (AMPL) extends the functionality of DeepChem and supports an array of machine learning and molecular featurization tools to predict key potency, safety and pharmacokinetic-relevant parameters. AMPL has been benchmarked on a large collection of pharmaceutical datasets covering a wide range of parameters. This is a living software project with active development. Check back for continued updates. Feedback is welcomed and appreciated, and the project is open to contributions! An article describing the AMPL project was published in JCIM. The AMPL pipeline documentation is available here.

Check out our new tutorial series that walks through AMPL's end-to-end modeling pipeline to build a machine learning model! View them in our docs or as Jupyter notebooks in our repo.

Static Badge

In addition to our written tutorials, we now provide a series of video tutorials on our YouTube channel, ATOMScience-org. These videos are created to assist users in exploring and leveraging AMPL's robust capabilities. We provided a playlist for easy streamlined Learning:

AMPL Tutorial Playlist


Table of contents

Useful links


Installation

AMPL supports Python 3.10 and provides platform-specific environments for cpu, cuda, rocm, and mchip.

For dependency-specific installation details, see:

AMPL uses uv for Python environment and dependency management.

Install the published package from PyPI

Use this option if you want to install the released atomsci-ampl package without cloning this repository.

If you are working from a local checkout of this repo, use the local uv environment workflow above instead.

Note: AMPL requires Python 3.10. For LLNL users on LC, run module load python/3.10.8.

Get the latest release version from PyPI

Create or activate a Python 3.10 environment, then install:

pip install atomsci-ampl

Or use uv install:

Example:

uv venv --python 3.10 <.venv_name>
source <.venv_name>/bin/activate
uv pip install atomsci-ampl

Set up an uv environment for local development

Use this workflow if you want to run or develop AMPL from a local clone of this repository.

What is uv?

uv is the Python environment and dependency tool used by this project. It creates virtual environments and installs the packages defined in pyproject.toml.

Requirements

Item Value
Python 3.10
Supported range >=3.10,<3.11
Platforms cpu, cuda, rocm, mchip

Install uv

Install uv with one of the following:

curl -LsSf https://astral.sh/uv/install.sh | sh

or

pip install uv

Verify the installation:

uv --version

Choose your platform

Platform Use this if... Environment
cpu you want a CPU-only environment .venv-cpu
cuda you are using NVIDIA GPUs on Linux .venv-cuda
rocm you are using AMD GPUs .venv-rocm
mchip you are on Apple Silicon .venv-mchip

Create or refresh your environment

For most users, use the make sync-<platform> commands:

make sync-cpu

Other options:

make sync-cuda
make sync-rocm
make sync-mchip

Then activate the environment:

source .venv-cpu/bin/activate

Replace cpu with your platform as needed.

Example:

make sync-cpu
source .venv-cpu/bin/activate

Daily use

If the environment is already working, you can usually just activate it:

source .venv-cpu/bin/activate

If commands like pytest are missing, or imports fail, refresh the environment:

make sync-cpu
source .venv-cpu/bin/activate

Platform lockfiles

This project keeps one lockfile per platform:

Platform Lockfile Environment
cpu uv.lock.cpu .venv-cpu
cuda uv.lock.cuda .venv-cuda
rocm uv.lock.rocm .venv-rocm
mchip uv.lock.mchip .venv-mchip

The make sync-<platform> commands use the matching lockfile to create or refresh .venv-<platform>.

For maintainers only

If dependencies change in pyproject.toml, you must regenerate the platform lockfile:

make update-lock-cpu

Other supported targets:

make update-lock-cuda
make update-lock-rocm
make update-lock-mchip

These commands create the lockfile if it does not already exist, or update it if it does.

Any changes to pyproject.toml or uv.lock.<platform> must be committed and pushed to the AMPL GitHub repository. Do not leave regenerated lockfiles uncommitted.

Most users do not need to run make update-lock-<platform>.

Troubleshooting

uv not found
uv --version

If this fails, install uv and start a new shell.

Wrong Python or missing pytest

Check:

which python
python --version
which pytest

They should point into the active environment, for example:

.venv-cpu/bin/python
.venv-cpu/bin/pytest

If not, refresh the environment:

make sync-cpu
source .venv-cpu/bin/activate
Import or package problems

Refresh the environment:

make sync-cpu
source .venv-cpu/bin/activate

Building AMPL for Local Development

If you want to develop AMPL locally from a repository checkout:

git clone https://github.com/ATOMScience-org/AMPL.git
cd AMPL
./install-dev.sh

This installs the package from the local source tree, so code changes are available without reinstalling.

Build and install from a local clone

If you want to build the package locally and then install the built artifact:

git clone https://github.com/ATOMScience-org/AMPL.git
cd AMPL
./build.sh
./install.sh

Optional, for LLNL LC only

If you use model_tracker, install atomsci.clients:

pip install -r clients_requirements.txt

Create jupyter notebook kernel (optional)

To run AMPL from Jupyter Notebook, first activate your environment and then run:

python -m ipykernel install --user --name atomsci-env

(Optional) LLNL LC only: if you use model_tracker, install atomsci.clients

# LLNL only: required for ATOM model_tracker
pip install -r clients_requirements.txt

Create jupyter notebook kernel (optional)

To run AMPL from Jupyter Notebook. To setup a new kernel, first activate your environment and then run the following command:

python -m ipykernel install --user --name atomsci-env

Install with Docker

  • Download and install Docker Desktop.
  • Create a workspace folder to mount with Docker environment and transfer files.
  • Get the Docker image and run it. Since 1.6.3, there are some changes with the AMPL Docker.

To retrieve, run version 1.6.2 or earlier, please specify the desired version tag:

docker pull atomsci/atomsci-ampl:v1.6.2
docker run -it -p 8888:8888 -v </local_workspace_folder>:</directory_in_docker> atomsci/atomsci-ampl:v1.6.2

For AMPL versions 1.6.3 and later, we offer downloadable images for various platforms (CPU, GPU or Linux/ARM64). To run a Docker container, be sure to append bash at the end of the command to open a bash session.

docker pull atomsci/atomsci-ampl:latest-<platform> # can be cpu, gpu, or arm (for arm64 chip)
docker run -it -p 8888:8888 -v </local_workspace_folder>:</directory_in_docker> atomsci/atomsci-ampl:latest-<platform> bash
#inside docker environment
jupyter-notebook --ip=0.0.0.0 --allow-root --port=8888 &
# -OR-
jupyter-lab --ip=0.0.0.0 --allow-root --port=8888 &
  • Visit the provided URL in your browser, ie
  • From the notebook, you may need to set the kernel that atomsci is installed ("atomsci-venv") in order to acccess the atomsci package.

For additional options related to building, running, and other Docker development tasks, please refer to Makefile.md.


To remove an entire virtual environment named "atomsci-env":

rm -rf $ENVROOT/atomsci-env

AMPL Features

AMPL enables tasks for modeling and prediction from data ingestion to data analysis and can be broken down into the following stages:

1. Data curation

  • Generation of RDKit molecular SMILES structures
  • Processing of qualified or censored data processing
  • Curation of activity and property values

2. Featurization

  • Extended connectivity fingerprints (ECFP)
  • Graph convolution latent vectors from DeepChem
  • Chemical descriptors from Mordred package
  • Descriptors generated by MOE (requires MOE license)

3. Model training and tuning

  • Test set selection
  • Cross-validation
  • Uncertainty quantification

4. Supported models

  • scikit-learn random forest models
  • XGBoost models
  • Fully connected neural networks
  • Graph convolution models

5. Visualization and analysis

  • Visualization and analysis tools
Details of running specific features are within the [parameter (options) documentation](atomsci/ddm/docs/PARAMETERS.md). More detailed documentation is in the [library documentation](https://ampl.readthedocs.io/en/latest/).

Running AMPL

AMPL can be run from the command line or by importing into Python scripts and Jupyter notebooks.

Python scripts and Jupyter notebooks

AMPL can be used to fit and predict molecular activities and properties by importing the appropriate modules. See the examples for more descriptions on how to fit and make predictions using AMPL.

Pipeline parameters

AMPL includes many parameters to run various model fitting and prediction tasks.

  • Pipeline options (parameters) can be set within JSON files containing a parameter list.
  • The parameter list with detailed explanations of each option can be found at atomsci/ddm/docs/PARAMETERS.md.
  • Example pipeline JSON files can be found in the tests directory and the example directory.

Library documentation

AMPL includes detailed docstrings and comments to explain the modules. Full HTML documentation of the Python library is available with the package at https://ampl.readthedocs.io/en/latest/.

More information on AMPL usage


Tests

AMPL includes a suite of software tests. This section explains how to run a very simple test that is fast to run. The Python test fits a random forest model using Mordred descriptors on a set of compounds from Delaney, et al with solubility data. A molecular scaffold-based split is used to create the training and test sets. In addition, an external holdout set is used to demonstrate how to make predictions on new compounds.

To run the Delaney Python script that curates a dataset, fits a model, and makes predictions, run the following commands:

source $ENVROOT/atomsci-env/bin/activate # activate your pip environment.
cd atomsci/ddm/test/integrative/delaney_RF
pytest

Note: This test generally takes a few minutes on a modern system

The important files for this test are listed below:

  • test_delany_RF.py: This script loads and curates the dataset, generates a model pipeline object, and fits a model. The model is reloaded from the filesystem and then used to predict solubilities for a new dataset.
  • config_delaney_fit_RF.json: Basic parameter file for fitting
  • config_delaney_predict_RF.json: Basic parameter file for predicting

More example and test information


Advanced AMPL usage

Command line

AMPL can fit models from the command line with:

python model_pipeline.py --config_file filename.json # [filename].json is the name of the config file

To get more info on an AMPL config file, please refer to:

Hyperparameter optimization

Hyperparameter optimization for AMPL model fitting is available to run on SLURM clusters or with [Optuna](https://optuna.readthedocs.io/) (Bayesian Optimization). To run Bayesian Optimization, the following steps can be followed.
  1. (Optional) Install Optuna with "pip install optuna"

  2. Pre-split your dataset with computed_descriptors if you want to use Mordred/MOE/RDKit descriptors.

  3. In the config JSON file, set the following parameters.

    • "hyperparam": "True"
    • "search_type": "optuna"
    • "descriptor_type": "mordred_filtered,rdkit_raw" (use comma to separate multiple values)
    • "model_type": "RF|20" (the number after | is the number of evaluations of Bayesian Optimization)
    • "featurizer": "ecfp,computed_descriptors" (use comma if you want to try multiple featurizers, note the RF and graphconv are not compatible)
    • "result_dir": "/path/to/save/the/final/results,/temp/path/to/save/models/during/optimization" (Two paths separated by a comma)

    RF model specific parameters:

    • "rfe": "uniformint|8,512", (RF number of estimators)
    • "rfd": "uniformint|8,512", (RF max depth of the decision tree)
    • "rff": "uniformint|8,200", (RF max number of features)

    Use the following schemes to define the searching domains

    method|parameter1,parameter2...

    method: supported searching schemes in Optuna include: choice, uniform, loguniform, uniformint. For details, see the Optuna documentation.

    parameters:

    • choice: all values to search from, separated by comma, e.g. choice|0.0001,0.0005,0.0002,0.001
    • uniform: low and high bound of the interval to search, e.g. uniform|0.00001,0.001
    • loguniform: low and high bound (in natural log) of the interval to search, e.g. loguniform|-13.8,-6.9
      • Note: For backwards compatibility, loguniform values are submitted to Optuna in log scale. Although Optuna supports non-log-scaled value ranges for log-uniform distributions, we maintain the original log-scaled specification format to ensure consistency with existing configurations.
    • uniformint: low and high bound of the interval as integers, e.g. uniformint|8,256

    NN model specific parameters:

    • "lr": "loguniform|-13.8,-6.9", (learning rate)
    • "ls": "uniformint|3|8,512", (layer_sizes)
      • The number between two bars (|) is the number of layers, namely 3 layers, each one with 8~512 nodes
      • Note that the number of layers (number between two |) can not be changed during optimization, if you want to try different number of layers, just run several optimizations.
    • "dp": "uniform|3|0,0.4", (dropouts)
      • 3 layers, each one has a dropout range from 0 to 0.4
      • Note that the number of layers (number between two |) can not be changed during optimization, if you want to try different number of layers, just run several optimizations.

    XGBoost model specific parameters:

    • "xgbg": "uniform|0,0.4", (xgb_gamma, Minimum loss reduction required to make a further partition on a leaf node of the tree)
    • "xgbl": "loguniform|-6.9,-2.3", (xgb_learning_rate, Boosting learning rate (xgboost's "eta"))
  4. Run hyperparameter search in batch mode or submit a slurm job.

    python hyperparam_search_wrapper.py --config_file filename.json
    
  5. Save a checkpoint to continue it later.

    To save a checkpoint file of the hyperparameter search job, you want to set the following two parameters.

    • "hp_checkpoint_save": "/path/to/the/checkpoint/file.pkl"
    • "hp_checkpoint_load": "/path/to/the/checkpoint/file.pkl"

    If the "hp_checkpoint_load" is provided, the hyperparameter search will continue from the checkpoint.


Advanced testing

Running all tests

To run the full set of tests, use Pytest from the test directory:

source $ENVROOT/atomsci-env/bin/activate # activate your pip environment. "atomsci" is an example here.
cd atomsci/ddm/test
pytest

Running SLURM tests

Several of the tests take some time to fit. These tests can be submitted to a SLURM cluster as a batch job. Example general SLURM submit scripts are included as `pytest_slurm.sh`.
source $ENVROOT/atomsci-env/bin/activate # activate your pip environment. "atomsci-env" is an example here.
cd atomsci/ddm/test/integrative/delaney_NN
sbatch pytest_slurm.sh
cd ../../../..
cd atomsci/ddm/test/integrative/wenzel_NN
sbatch pytest_slurm.sh

Running tests without internet access

AMPL works without internet access. Curation, fitting, and prediction do not require internet access.

However, the public datasets used in tests and examples are not included in the repo due to licensing concerns. These are automatically downloaded when the tests are run.

If a system does not have internet access, the datasets will need to be downloaded before running the tests and examples. From a system with internet access, run the following shell script to download the public datasets. Then, copy the AMPL directory to the offline system.

cd atomsci/ddm/test
bash download_datset.sh
cd ../../..
# Copy AMPL directory to offline system

AMPL tutorials

Please follow link, "atomsci/ddm/examples/tutorials", to access a collection of AMPL tutorial notebooks. The tutorial notebooks give an exhaustive coverage of AMPL features. The AMPL team has prepared the tutorials to help beginners understand the basics to advanced AMPL features, and a reference for advanced AMPL users.


Development

Installing the AMPL for development

Using "pip install -e ." will create a namespace package in your environment directory that points back to your git working directory, so every time you reimport a module you'll be in sync with your working code. Since site-packages is already in your sys.path, you won't have to fuss with PYTHONPATH or setting sys.path in your notebooks.

Code Push Policy

It's recommended to use a development branch to do the work. After each release, there will be a branch opened for development.

The policy is

  1. Create a branch based off a development ("1.6.0 "for example) or "master" branch
  2. Create a pull request. Assign a reviewer to approve the code changes

Note: Step 2 is required for pushing directly to "master". For a development branch, this step is recommended but not required.

Docstring format

The "Google docstring" format is used in the AMPL code. When writing new code, please use the same Docstring style. Refer here and here for examples.

Versioning

Versions are managed through GitHub tags on this repository.

Built with

  • DeepChem: A rich repository of chemistry-specific model types and utilities
  • RDKit: Molecular informatics library
  • Mordred: Chemical descriptors
  • Other Python package dependencies

Project information

Authors

The Accelerating Therapeutics for Opportunities in Medicine (ATOM) Consortium

  • Amanda J. Minnich (1)
  • Kevin McLoughlin (1)
  • Margaret Tse (2)
  • Jason Deng (2)
  • Andrew Weber (2)
  • Neha Murad (2)
  • Benjamin D. Madej (3)
  • Bharath Ramsundar (4)
  • Tom Rush (2)
  • Stacie Calad-Thomson (2)
  • Jim Brase (1)
  • Jonathan E. Allen (1)  

Contributors

  • Amanda Paulson (5)
  • Stewart He (1)
  • Da Shi (6)
  • Ravichandran Sarangan (7)
  • Jessica Mauvais (1)

1. Lawrence Livermore National Laboratory
2. GlaxoSmithKline Inc.
3. Frederick National Laboratory for Cancer Research
4. Computable
5. University of California, San Francisco
6. Schrodinger
7. Leidos  

Support, Suggestions or Report Issues

  • If you have suggestions or like to report issues, please click here.  

Contributing

Thank you for contributing to AMPL!

  • Contributions must be submitted through pull requests.
  • All new contributions must adhere to the MIT license.  

Release

AMPL is distributed under the terms of the MIT license. All new contributions must be made under this license.

See MIT license and NOTICE for more details.

  • LLNL-CODE-795635
  • CRADA TC02264

About

The ATOM Modeling PipeLine (AMPL) is an open-source, modular, extensible software pipeline for building and sharing models to advance in silico drug discovery.

Topics

Resources

License

Code of conduct

Contributing

Stars

Watchers

Forks

Packages

 
 
 

Contributors