pyMBE: the Python-based Molecule Builder for ESPResSo

pyMBE provides tools to facilitate building up molecules with complex architectures in the Molecular Dynamics software ESPResSo. Some examples of molecules that can be set up with pyMBE are polyelectrolytes, peptides and proteins. pyMBE bookkeeps all the information about the molecule topology, permitting to link each particle to its corresponding residue and molecule. pyMBE uses the Pint library to enable input parameters in any arbitrary unit system, which is later transformed in the reduced unit system used in ESPResSo.

Dependencies

• ESPResSo =4.2.1
• Pint >=0.20.01
• Pandas >=1.5.3
• Pint-Pandas >=0.3
• Numpy >=1.23
• SciPy
• pdoc (for building the docs)
• CMake (for running the testsuite)

Contents

• figs/: folder with various images used in the tutorials of pyMBE.
• lib/: folder with various libraries.
• maintainer/: folder with various scripts used by the maintainers.
• parameters/: folder with various sets of parameters from previous works.
• samples/: folder with various sample scripts showcasing how to use pyMBE to setup different systems.
• testsuite/: folder with several test scripts and data for continous integration of the library.
• tutorials/: folder with the available tutorials on pyMBE.
• visualization/: folder with helper scripts to aid the visualization of vtf trajectories from constant pH and Grand reaction simulations with VMD.
• AUTHORS.md: list of authors and contributors of pyMBE.
• CONTRIBUTING.md: rules on how to contribute to pyMBE.
• LICENSE.txt: license of pyMBE.
• pyMBE.py: source code of pyMBE
• requirements.txt: list of required libraries to use pyMBE.

Usage

Set up the pyMBE virtual environment

To use pyMBE in your simulations, first clone this repository locally:

git clone [email protected]:pyMBE-dev/pyMBE.git

Please, be aware that pyMBE is intended to be a supporting tool to setup simulations with ESPResSo. Thus, for most of its functionalities ESPResSo must also be available. Following the NEP29 guidelines, we recommend the users of pyMBE to use Python3.10+ when using our module.

The pyMBE module uses its own Python virtual enviroment to avoid incompatibility issues when loading its requirements from other libraries. The Python module venv is needed to set up pyMBE. If venv is not in the Python distribution of the user, the user will need to first install 'venv' before setting up pyMBE. For Ubuntu users, this can be done as follows:

sudo apt install python3-venv

To set up pyMBE, the users need to install its virtual environment, install its Python dependencies and configure the path to the ESPResSo build folder as follows:

python3 -m venv pymbe  # creates a local folder named pymbe, which contains the virtual environment
source pymbe/bin/activate  # activates the pymbe venv
python3 maintainer/configure_venv.py --espresso_path=/home/user/espresso/build # please, adapt the espresso path accordingly
python3 -m pip install -r requirements.txt
python3 simulation_script.py # run the espresso simulation script
deactivate  # deactivate the virtual environment

We highlight that the path /home/user/espresso/build is just an example of a possible path to the ESPResSo build folder. The user should change this path to match the local absolute path were ESPResSo was installed. For more details on how to install ESPResSo, please consult the ESPResSo installation guide.

The pyMBE virtual enviroment can be deactivated at any moment:

deactivate

Cluster users who rely on module files to load dependencies should opt for the following alternative:

module load ESPResSo/4.2.1-foss-2022a # adapt module name
python3 -m venv --system-site-packages pymbe
source pymbe/bin/activate
python3 maintainer/configure_venv.py
python3 -m pip install -r requirements.txt
deactivate
module purge

We highlight that the module files need to be loaded before every activation of the virtual environment.

Now you can use pyMBE and ESPResSo by activating the virtual environment:

$ source pymbe/bin/activate
(pymbe) $ python3 -c "import espressomd.version; print(espressomd.version.friendly())"
4.2
(pymbe) $ python3 -c "import pyMBE; print(pyMBE.__file__)"
/home/user/Documents/pyMBE/pyMBE.py
$ deactivate

To use pyMBE in JupyterLab, register the virtual environment in a new kernel:

source pymbe/bin/activate
python3 -m pip install ipykernel "jupyterlab>=4.0.8" "PyOpenGL>=3.1.5"
python3 -m ipykernel install --user --name=pyMBE
deactivate

Please be aware the pyMBE kernel will be registered outside the environment, typically in your home folder. You can later inspect the list of registered kernels and delete unwanted ones with the following commands:

jupyter kernelspec list
jupyter kernelspec uninstall pymbe

The JupyterLab main menu will now show a new Python kernel called "pyMBE" that uses the virtual environment.

Use pyMBE in your simulation scripts

source pymbe/bin/activate
python3 samples/peptide.py
deactivate

Run the tutorial of pyMBE

You can run the interactive tutorial of pyMBE with the command:

source pymbe/bin/activate
jupyter-lab tutorials/pyMBE_tutorial.ipynb
deactivate

Be sure to use the pyMBE kernel instead of the default Python3 kernel. The currently active kernel is usually displayed in the top right corner of the notebook.

Run the testsuite

To make sure your code is valid, please run the testsuite before submitting your contribution:

source pymbe/bin/activate
make tests -j4
deactivate

Here, -j4 instructs CTest to run the test cases in parallel using 4 CPU cores. This number can be adjusted depending on your hardware specifications. You can use make unit_tests -j4 to run the subset of fast tests, but keep in mind those won't be able to detect more serious bugs that only manifest themselves in long simulations. You can also run individual test cases directly, for example with python3 testsuite/parameter_test.py.

When contributing new features, consider adding a unit test in the testsuite/ folder and a corresponding line in the testsuite/CTestTestfile.cmake file.

Every contribution is automatically tested in CI using EESSI (https://www.eessi.io) and the EESSI GitHub Action.

References

Check out the corresponding paper to learn more about pyMBE. If you use pyMBE in your research, please cite our paper:

@article{beyer2024pymbe,
  author = {Beyer, David and Torres, Paola B. and Pineda, Sebastian P. and
            Narambuena, Claudio F. and Grad, Jean-No{\"e}l and Ko{\v{s}}ovan,
            Peter and Blanco, Pablo M.},
  title = {{pyMBE}: The {P}ython-based molecule builder for {ESPResSo}},
  journal = {The Journal of Chemical Physics},
  volume = {161},
  number = {2},
  pages = {022502},
  year = {2024},
  month = jul,
  issn = {0021-9606},
  doi = {10.1063/5.0216389},
}

When using a released version of pyMBE, we recommend citing the corresponding Zenodo record in addition to the pyMBE paper, for example: "We set up our coarse-grained models using pyMBE v0.8.0 [@beyer2024pymbe; @zenodo2024pymbe]".

Please also make sure to properly cite the original authors if you use the resources provided in the parameters/ folder. The relevant references are provided as metadata in the corresponding files.

License

Copyright (C) 2023-2024 pyMBE-dev team

pyMBE is free software: you can redistribute it and/or modify it under the terms of the GNU General Public License as published by the Free Software Foundation, either version 3 of the License, or (at your option) any later version.

You should have received a copy of the GNU General Public License along with this program. If not, see https://www.gnu.org/licenses/.