Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Using customized force field for small molecules #851

Open
wenchangzhou-qtx opened this issue May 10, 2024 · 8 comments
Open

Using customized force field for small molecules #851

wenchangzhou-qtx opened this issue May 10, 2024 · 8 comments

Comments

@wenchangzhou-qtx
Copy link

Hey @IAlibay @mikemhenry,

For small molecules, can we now using customized force fields provided by users like myself, other than gaff or espaloma?

Thanks

@mikemhenry
Copy link
Contributor

We support custom force fields if they comply with the smirnoff spec so anything that openff makes or you create with bespoke fit should work.

@aminsagar
Copy link

@mikemhenry
I have also been trying to use espaloma as the small molecule force field. #682
It seems it should be possible to use espaloma #229
Can you please give a minimal example of using espaloma?

Thanks.

@mikemhenry
Copy link
Contributor

Ah! Okay yes this is supported, I don't have time right now to test this, but it should be as easy as specifying espaloma-0.3.2 as the forcefeield you want to use, is that what you mean?

OH! I see in #682 it looks like openmmforcefeild needs an update to support espaloma-0.3.2 since it only supports espaloma-0.2.2 currently. Tracking that issue here openmm/openmmforcefields#336

@wenchangzhou-qtx
Copy link
Author

Thanks! @mikemhenry, on this topic, is it possible to change it inside yaml file like this?

protocol:
  settings:
    forcefield_settings:
      small_molecule_forcefield: 'espaloma-0.2.2'
    protocol_repeats: 2

@mikemhenry
Copy link
Contributor

Yes that should work, if it doesn't please let me know!

@wenchangzhou-qtx
Copy link
Author

wenchangzhou-qtx commented May 14, 2024

Hey @mikemhenry, thanks for checking, I just tried but in the output json file I still have the default, openff-2.1.1 for small_molecule_forcefield and 3 for protocol_repeats, my settings.yaml file looks like this:

mapper:
  method: lomap
  # method: kartograf

network:
  method: generate_radial_network
  settings:
    central_ligand: 'lig_ejm_31'
  # method: generate_radial_network
  # method: generate_maximal_network
  # method: generate_minimal_redundant_network

protocol:
  settings:
    forcefield_settings:
      small_molecule_forcefield: 'espaloma-0.2.2'
    protocol_repeats: 2

The command I run is openfe plan-rbfe-network -M tyk2_ligands.sdf -p tyk2_protein.pdb -o network_setup -s settings.yaml

@IAlibay
Copy link
Member

IAlibay commented May 14, 2024

@wenchangzhou-qtx, CLI YAML inputs for Protocols is not currently available in the 1.0 release of openfe. It should make it soon but for now if you want non-default settings you'll have to use the Python API.

@wenchangzhou-qtx
Copy link
Author

wenchangzhou-qtx commented May 14, 2024 via email

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

No branches or pull requests

4 participants