Making-it-rain

Cloud-based molecular simulations for everyone

+ 28Apr2023: We have updated all the notebooks to Python 3.10, everything is working fine
+ 01May2023: We have fixed the Amber, AlphaFol2+MD and Protei-Ligand notebooks.
+ 15Dec2023: We have corrected all the notebooks, and everything should be functioning properly now.

Hello there!

This is a repository where you can find a Jupyter notebook scripts for running Molecular Dynamics (MD) simulations using OpenMM engine and AMBER and CHARMM force fields files on Google Colab. This repository is a supplementary material of the paper "Making it rain: Cloud-based molecular simulations for everyone" and we encourage you to read it before using this pipeline.

The main goal of this work is to demonstrate how to harness the power of cloud-computing to run microsecond-long MD simulations in a cheap and yet feasible fashion.

Important: We've updated the notebooks to CondaColab. Now, all the dependencies will be installed faster than before (less than half of the previous time). You will see a CondaColab cell, just run and wait a few seconds, the session will restart and this is normal and expected. After that, you can continue running the cells like normal. Do not use the Run all option. Run the condacolab cell individually and wait for the kernel to restart. Only then, you can run all cells if you want. Thank you for your support.

1. AMBER - Using AMBER to generate topology and to build the simulation box
2. CHARMM - Using inputs from CHARMM-GUI solution builder
3. AlphaFold2+MD - Using AlphaFold2_mmseqs2 to generate protein model + MD simulation using AMBER to generate topology and to build simulation box

UPDATE (October 2021)

4. Protein-Ligand simulations - Using AMBER to generate topology and to build the simulation box and for the ligand using GAFF2 force field
5. Using AMBER Inputs - Using inputs from AMBER suite of biomolecular simulation program
6. Using GROMACS Inputs - Using inputs from GROMACS biomolecular simulation package (AMBER, CHARMM and OPLS force fields are compatible)

UPDATE (March 2022)

7. RESP Partial Charges - Using a SMILES as input and outputs a mol2 file with RESP derived partial charges. Options for setting method (HF, B3LYP, ...), basis set (3-21G, 6-31G*) and singlepoint or geometry optimization are available
8. Small Molecules MD - Using a SMILES as a input, geometry optimization with TorchANI and topology with AMBER (GAFF2 force field)
9. GLYCAM - Using inputs from GLYCAM server

UPDATE (August 2022)

10. DRUDE - Using inputs from CHARMM-GUI Drude Prepper

UPDATE (March 2024)

11. Protein-Membrane simulations - Using OpenFF to generate the topology and build the simulation box for protein-membrane systems with AMBER force fields.
12. Martini+cg2all - Utilizing Vermouth, the Python library that powers Martinize2, to generate the topology and build the simulation box for protein systems using Martini force fields. Additionally, employing cg2at to predict all-atom trajectories from coarse-grained (CG) representations.
13. AMBER Mutations - Performing mutations on protein/nucleic acid systems and utilizing AMBER to generate the topology and build the simulation box.

Tired to "just" run molecular simulations and want to try something new? Gabriel Monteiro da Silva (@GMondaSilva) and I are thrilled to share a colab notebook for running the subsampled AlphaFold2 approach for predicting protein conformational ensembles. Check it out:

Bugs

• If you encounter any bugs, please report the issue to https://github.com/pablo-arantes/Making-it-rain/issues

Acknowledgments

• We would like to thank the Psi4 team for developing an excellent and open source suite of ab initio quantum chemistry.
• We would like to thank the Roitberg team for developing the fantastic TorchANI.
• We would like to thank the OpenMM team for developing an excellent and open source engine.
• We would like to thank the AlphaFold team for developing an excellent model and open sourcing the software.
• We would like to thank the ChemosimLab (@ChemosimLab) team for their incredible ProLIF (Protein-Ligand Interaction Fingerprints) tool.
• Credits to Sergey Ovchinnikov (@sokrypton), Milot Mirdita (@milot_mirdita) and Martin Steinegger (@thesteinegger) for their fantastic ColabFold
• Making it rain by Pablo R. Arantes (@pablitoarantes), Marcelo D. Polêto (@mdpoleto), Conrado Pedebos (@ConradoPedebos) and Rodrigo Ligabue-Braun (@ligabue_braun).
• Also, credit to David Koes for his awesome py3Dmol plugin.
• Finally, we would like to thank Professor Giulia Palermo for her support and thoughtful comments in the development of the present work.

Do you want to cite this work?

Arantes P.R., Depólo Polêto M., Pedebos C., Ligabue-Braun R. Making it rain: cloud-based molecular simulations for everyone. Journal of Chemical Information and Modeling 2021. DOI: 10.1021/acs.jcim.1c00998.