Skip to content

Commit 8ef28ae

Browse files
authored
Update README.md
1 parent ddda87e commit 8ef28ae

1 file changed

Lines changed: 23 additions & 18 deletions

File tree

README.md

Lines changed: 23 additions & 18 deletions
Original file line numberDiff line numberDiff line change
@@ -1,31 +1,30 @@
1-
# Graph Attentional Protein Parametrization (GrAPPa)
1+
# Grappa (Graph Attentional Protein Parametrization)
22

33
_A machine-learned molecular mechanics force field using a deep graph attentional network <br>(code supporting [https://arxiv.org/abs/2404.00050](https://arxiv.org/abs/2404.00050))_
44

55
## Abstract
66

7-
Simulating large molecular systems over long timescales requires force fields that are both accurate and efficient. In recent years, E(3) equivariant neural networks have lifted the tension between computational efficiency and accuracy of force fields, but they are still several orders of magnitude more expensive than classical molecular mechanics (MM) force fields.
8-
9-
Here, we propose a novel machine learning architecture to predict MM parameters from the molecular graph, employing a graph attentional neural network and a transformer with symmetry-preserving positional encoding. The resulting force field, Grappa, outperforms established and other machine-learned MM force fields in terms of accuracy at the same computational efficiency and can be used in existing Molecular Dynamics (MD) engines like GROMACS and OpenMM. It predicts energies and forces of small molecules, peptides, RNA and - showcasing its extensibility to uncharted regions of chemical space - radicals at state-of-the-art MM accuracy. We demonstrate Grappa's transferability to macromolecules in MD simulations, during which large protein are kept stable and small proteins can fold. Our force field sets the stage for biomolecular simulations close to chemical accuracy, but with the same computational cost as established protein force fields.
7+
Simulating large molecular systems over long timescales requires force fields that are both accurate and efficient.
8+
In recent years, E(3) equivariant neural networks have lifted the tension between computational efficiency and accuracy of force fields, but they are still several orders of magnitude more expensive than established molecular mechanics (MM) force fields.
9+
Here, we propose Grappa, a machine learning framework to predict MM parameters from the molecular graph, employing a graph attentional neural network and a transformer with symmetry-preserving positional encoding.
10+
The resulting Grappa force field outperformstabulated and machine-learned MM force fields in terms of accuracy at the same computational efficiency and can be used in existing Molecular Dynamics (MD) engines like GROMACS and OpenMM.
11+
It predicts energies and forces of small molecules, peptides, RNA and - showcasing its extensibility to uncharted regions of chemical space - radicals at state-of-the-art MM accuracy.
12+
We demonstrate Grappa's transferability to macromolecules in MD simulations from a small fast folding protein up to a whole virus particle. Our force field sets the stage for biomolecular simulations closer to chemical accuracy, but with the same computational cost as established protein force fields.
1013

1114
<details open>
1215
<summary>Grappa Overview</summary>
1316
<p align="center">
1417
<img src="docs/figures/grappa_overview.png" width="50%" style="max-width: 200px; display: block; margin: auto;">
1518
</p>
16-
<p><i>Grappa first predicts node embeddings from the molecular graph. In a second step, it predicts MM parameters for each n-body interaction from the embeddings of the contributing nodes, respecting the necessary permutation symmetry.</i></p>
17-
</details>
18-
19-
<details>
20-
<summary><b>Performance on MM Benchmark Datasets</b></summary>
21-
<p align="center">
22-
<img src="docs/figures/table_benchmark.png" width="100%" style="max-width: 200px; display: block; margin: auto;">
23-
</p>
24-
<p><i>Grappa's energy and force-component RMSE in kcal/mol and kcal/mol/Å on the test dataset (trained with the same train-val-test partition) from Espaloma [<a href="https://arxiv.org/abs/2307.07085v4">Takaba et al. 2023</a>], compared with classical forcefields [<a href="https://pubs.aip.org/aip/jcp/article/153/11/114502/199591/A-fast-and-high-quality-charge-model-for-the-next">He et al.</a>], [<a href="https://doi.org/10.1021/acs.jctc.5b00255">Maier et al.</a>, <a href="https://pubs.acs.org/doi/10.1021/ct200162x">Zgarbova et al.</a>]</i></p>
19+
<p><i>
20+
Grappa predicts MM parameters in two steps.
21+
First, atom embeddings are predicted from the molecular graph with a graph neural network.
22+
Then, transformers with symmetric positional encoding followed by permutation invariant pooling maps the embeddings to MM parameters with desired permutation symmetries.
23+
Once the MM parameters are predicted, the potential energy surface can be evaluated with MM-efficiency for different spatial conformations.
24+
</i></p>
2525
</details>
2626

2727

28-
2928
<details open><summary><b>Table of contents</b></summary>
3029

3130
- [Usage](#usage)
@@ -108,11 +107,15 @@ conda create -n grappa python=3.10 -y
108107
conda activate grappa
109108
```
110109

111-
In cpu mode, Grappa is available on pypi:
110+
In cpu mode, Grappa is available on PyPi:
112111
```{bash}
113112
pip install grappa-ff
114113
```
115114

115+
Depending on the platform used, installation of OpenMM or GROMACS and Kimmdy is needed (see below).
116+
117+
### Installation from source
118+
116119
To install Grappa from source, clone the repository and install requirements and the package itself with pip:
117120

118121
```{bash}
@@ -130,15 +133,15 @@ python tests/test_installation.py
130133

131134
### GROMACS
132135

133-
The creation of custom GROMACS topology files is handled by [Kimmdy](https://github.com/hits-mbm-dev/kimmdy), which can be installed in the same environment as grappa via pip,
136+
The creation of custom GROMACS topology files is handled by [Kimmdy](https://github.com/hits-mbm-dev/kimmdy), which can be installed in the same environment as Grappa via pip,
134137

135138
```{bash}
136139
pip install kimmdy==6.8.3
137140
```
138141

139142
### OpenMM
140143

141-
Unfortunately, OpenMM is not available on pip and has to be installed via conda in the same environment as grappa,
144+
OpenMM is not available on pip and has to be installed via conda in the same environment as Grappa,
142145

143146
```{bash}
144147
conda install -c conda-forge openmm
@@ -149,7 +152,9 @@ Since the resolution of package dependencies can be slow in conda, it is recomme
149152

150153
### GPU mode
151154

152-
For training grappa models, neither OpenMM nor Kimmdy ar needed, only an environment with a working installation of [PyTorch](https://pytorch.org/) and [DGL](https://www.dgl.ai/) for the cuda version of choice. Instructions for installing dgl with cuda can be found at `installation/README.md`.
155+
For training Grappa models, neither OpenMM nor Kimmdy ar needed, only an environment with a working installation of [PyTorch](https://pytorch.org/) and [DGL](https://www.dgl.ai/) for the cuda version of choice.
156+
Note that installing Grappa in GPU mode is only recommended if training a model is intended.
157+
Instructions for installing dgl with cuda can be found at `installation/README.md`.
153158
In this environment, Grappa can be installed by
154159

155160
```{bash}

0 commit comments

Comments
 (0)