Skip to content

[ WIP ] Add tutorial for membrane RBFE calculation#286

Open
hannahbaumann wants to merge 2 commits intomainfrom
membrane_tutorial
Open

[ WIP ] Add tutorial for membrane RBFE calculation#286
hannahbaumann wants to merge 2 commits intomainfrom
membrane_tutorial

Conversation

@hannahbaumann
Copy link
Contributor

No description provided.

@review-notebook-app
Copy link

Check out this pull request on  ReviewNB

See visual diffs & provide feedback on Jupyter Notebooks.


Powered by ReviewNB

@github-actions
Copy link

github-actions bot commented Feb 6, 2026

Colab 👈 Launch a Colab session on branch membrane_tutorial

@hannahbaumann hannahbaumann changed the title Add tutorial for membrane RBFE calculation [ WIP ] Add tutorial for membrane RBFE calculation Feb 6, 2026
@@ -0,0 +1,1153 @@
{
Copy link
Contributor

@jthorton jthorton Feb 6, 2026

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Line #25.                settings=RelativeHybridTopologyProtocol._adaptive_settings(

Do we not need to add the lipid force field to the settings or is adaptive doing this? Also this is another place we are exposing _adaptive to the world do we need to add a public entry point?


Reply via ReviewNB

@@ -0,0 +1,1153 @@
{
Copy link
Contributor

@jthorton jthorton Feb 6, 2026

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Line #43.    network = openfe.AlchemicalNetwork(transformations)

Woow so many warnings we need to remove a lot of these!


Reply via ReviewNB

Copy link
Member

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Already done - Hannah's branch isn't up to date with main.

@@ -0,0 +1,1153 @@
{
Copy link
Member

@IAlibay IAlibay Feb 7, 2026

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Two out of three of these imports warning seems to be some kind of problem with our typing vendor in gufe. I have opened an issue.


Reply via ReviewNB

@@ -0,0 +1,1153 @@
{
Copy link
Member

@IAlibay IAlibay Feb 7, 2026

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

For this example, rather than separating this in a separate step, I think it will be more useful to the readers if you just do these calls in-place, i.e. modeller = openmm.app.Modeller(protein.to_openmm_topology(), protein.to_openmm_positions().

With a little bit of formatting, that will look a lot less overwhelming.


Reply via ReviewNB

@@ -0,0 +1,1153 @@
{
Copy link
Member

@IAlibay IAlibay Feb 7, 2026

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

This is too unweidly for users in my opinion, I think OpenFreeEnergy/gufe#568 will need to be completed.


Reply via ReviewNB

@@ -0,0 +1,1153 @@
{
Copy link
Member

@IAlibay IAlibay Feb 7, 2026

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Generally would not recommend overwriting an existing file in your example, that's going to lead users to do bad things on their end - I would instead suggest you set a new file name for it.


Reply via ReviewNB

@@ -0,0 +1,1153 @@
{
Copy link
Member

@IAlibay IAlibay Feb 7, 2026

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Could you somewhere make a statement about which of these would be applicable to the "merged PDB" example shown above? (give that most Maestro users would be going down that path).

numpy array using OpenFF units (see below)

This is the wrong level of details at this point, here I would just say "box vectors in OpenMM format". Then in the section below, I would precise that it has to be a united array.


Reply via ReviewNB

@@ -0,0 +1,1153 @@
{
Copy link
Member

@IAlibay IAlibay Feb 7, 2026

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Here mention that you need these to be in an openff-units united array.


Reply via ReviewNB

@@ -0,0 +1,1153 @@
{
Copy link
Member

@IAlibay IAlibay Feb 7, 2026

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Add some details on how this is done? Any caveats that one might need to consider (i.e. making sure that the system is wrapped / not broken across periodic images)?


Reply via ReviewNB

@@ -0,0 +1,1153 @@
{
Copy link
Member

@IAlibay IAlibay Feb 7, 2026

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Add some details on the types of lipid forcefields that are currently supported (i.e. I know there are some that won't work, which currently do work and which should one avoid?)


Reply via ReviewNB

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment

Labels

None yet

Projects

None yet

Development

Successfully merging this pull request may close these issues.

3 participants