# How to add an organic molecule in an enzym's active pocket?

**URL:** https://gromacs.bioexcel.eu/t/how-to-add-an-organic-molecule-in-an-enzyms-active-pocket/7272
**Category:** User discussions
**Tags:** analysis-tools, simulation-setup
**Created:** [September 26, 2023, 9:05pm UTC](https://gromacs.bioexcel.eu/t/how-to-add-an-organic-molecule-in-an-enzyms-active-pocket/7272 "2023-09-26T21:05:47Z")
**Posts on this page:** 3
**Page:** 1

<div class="post-metadata">

### Author: ![Ali\_Jafari](https://avatars.discourse-cdn.com/v4/letter/a/ee59a6/32.png) [@Ali\_Jafari](https://gromacs.bioexcel.eu/u/Ali_Jafari)
#### Post date: [September 26, 2023, 9:05pm UTC](https://gromacs.bioexcel.eu/t/how-to-add-an-organic-molecule-in-an-enzyms-active-pocket/7272/1 "2023-09-26T21:05:47Z")

</div>

GROMACS version: 2020.7  
GROMACS modification: No  
I’ve been trying to run a MD simulation on a ligand-enzyme (in my case CALB) complex in acetonitrile as solvent but the problem is that I don’t know how to insert my ligand in the active pocket since it’s not a natural ligand for CALB.  
Can I make this complex using molecular docking softwares like Glide in Schrödinger suite and then use its results to run my calculation? (Would it be reliable?) I appretiate any help in advance.

P. S. : My main reason for confusion is that I’m using CHARMM36-jul2022.ff for MD simulations but molecular docking softwares like Glide in Schrödinger suite employ FFs like OPLS3

---

<div class="post-metadata">

### Author: ![MagnusL](https://dub1.discourse-cdn.com/flex017/user_avatar/gromacs.bioexcel.eu/magnusl/32/2380_2.png) [@MagnusL](https://gromacs.bioexcel.eu/u/MagnusL)
#### Post date: [September 27, 2023, 6:58am UTC](https://gromacs.bioexcel.eu/t/how-to-add-an-organic-molecule-in-an-enzyms-active-pocket/7272/2 "2023-09-27T06:58:25Z")

</div>

Using output coordinates from a docking software is almost certainly the most reliable way to go. Of course, the quality depends on how accurate the docking poses are, but if there are no experimental structures available then docking poses will give you good starting positions.

That the force fields, for docking and MD simulations, are different should, in most cases, not make a dramatic difference, as long as both are good enough at representing the molecules you are studying. You will need to run energy minimization, and equilibration on your complex, before starting any thorough analyses. Since you need to add water between docking and MD simulations that would be necessary anyhow.

---

<div class="post-metadata">

### Author: ![Ali\_Jafari](https://avatars.discourse-cdn.com/v4/letter/a/ee59a6/32.png) [@Ali\_Jafari](https://gromacs.bioexcel.eu/u/Ali_Jafari)
#### Post date: [September 27, 2023, 7:16am UTC](https://gromacs.bioexcel.eu/t/how-to-add-an-organic-molecule-in-an-enzyms-active-pocket/7272/3 "2023-09-27T07:16:59Z")

</div>

Thank you so much for your answer 🙏  
It was very informative
