# how to  distribute 20 polymers (ligand)  randomly  around the target protein to study interactions

**URL:** <https://gromacs.bioexcel.eu/t/how-to-distribute-20-polymers-ligand-randomly-around-the-target-protein-to-study-interactions/4700>\
**Category:** User discussions\
**Created:** [August 29, 2022, 2:33pm UTC](https://gromacs.bioexcel.eu/t/how-to-distribute-20-polymers-ligand-randomly-around-the-target-protein-to-study-interactions/4700 "2022-08-29T14:33:50Z")\
**Posts on this page:** 7\
**Page:** 1

<div class="post-metadata">

**Author:** ![mehdi60](https://avatars.discourse-cdn.com/v4/letter/m/a698b9/32.png) [@mehdi60](https://gromacs.bioexcel.eu/u/mehdi60)\
**Post date:** [August 29, 2022, 2:33pm UTC](https://gromacs.bioexcel.eu/t/how-to-distribute-20-polymers-ligand-randomly-around-the-target-protein-to-study-interactions/4700/1 "2022-08-29T14:33:50Z")

</div>

GROMACS version:2021.4  
GROMACS modification: No  
Hello everyone, I am new to GROMACS and trying to run protein-ligand simulations to consider protein-polymer interaction for the drug delivery study. So I want to randomly distribute 20 polymers around the target protein. How can I do it? Please kindly guide me in this regard. (I’m using gromos96 54a7 force field via gromacs 2021.4)

Thank you

---

<div class="post-metadata">

**Author:** ![jalemkul](https://dub1.discourse-cdn.com/flex017/user_avatar/gromacs.bioexcel.eu/jalemkul/32/18_2.png) [@jalemkul](https://gromacs.bioexcel.eu/u/jalemkul)\
**Post date:** [August 30, 2022, 4:57pm UTC](https://gromacs.bioexcel.eu/t/how-to-distribute-20-polymers-ligand-randomly-around-the-target-protein-to-study-interactions/4700/2 "2022-08-30T16:57:57Z")

</div>

Use `gmx insert-molecules`.

---

<div class="post-metadata">

**Author:** ![mehdi60](https://avatars.discourse-cdn.com/v4/letter/m/a698b9/32.png) [@mehdi60](https://gromacs.bioexcel.eu/u/mehdi60)\
**Post date:** [April 17, 2023, 7:21pm UTC](https://gromacs.bioexcel.eu/t/how-to-distribute-20-polymers-ligand-randomly-around-the-target-protein-to-study-interactions/4700/3 "2023-04-17T19:21:24Z")

</div>

Thank you very much for your kind reply.  
I am going to add 3 types of molecules (A, B and C) to my simulation box, so that first A and B molecules are added to the box, and after they interact with each other, C is added, and C interacts with the A&B mixture. For this purpose, I plan to do the following process in `gmx insert-molecules` .

My procedure: With each call of `gmx insert-molecules` I insert one type molecule, then use the output from that step as the input for the next to insert another type of molecule. Is it correct?

Thanks, regards

---

<div class="post-metadata">

**Author:** ![jalemkul](https://dub1.discourse-cdn.com/flex017/user_avatar/gromacs.bioexcel.eu/jalemkul/32/18_2.png) [@jalemkul](https://gromacs.bioexcel.eu/u/jalemkul)\
**Post date:** [April 18, 2023, 2:35am UTC](https://gromacs.bioexcel.eu/t/how-to-distribute-20-polymers-ligand-randomly-around-the-target-protein-to-study-interactions/4700/4 "2023-04-18T02:35:46Z")

</div>

> [@mehdi60](#):
>
> My procedure: With each call of `gmx insert-molecules` I insert one type molecule, then use the output from that step as the input for the next to insert another type of molecule. Is it correct?

This is correct.

---

<div class="post-metadata">

**Author:** ![mehdi60](https://avatars.discourse-cdn.com/v4/letter/m/a698b9/32.png) [@mehdi60](https://gromacs.bioexcel.eu/u/mehdi60)\
**Post date:** [April 18, 2023, 10:25am UTC](https://gromacs.bioexcel.eu/t/how-to-distribute-20-polymers-ligand-randomly-around-the-target-protein-to-study-interactions/4700/5 "2023-04-18T10:25:27Z")

</div>

Thanks a lot…

---

<div class="post-metadata">

**Author:** ![mehdi60](https://avatars.discourse-cdn.com/v4/letter/m/a698b9/32.png) [@mehdi60](https://gromacs.bioexcel.eu/u/mehdi60)\
**Post date:** [June 16, 2023, 11:54pm UTC](https://gromacs.bioexcel.eu/t/how-to-distribute-20-polymers-ligand-randomly-around-the-target-protein-to-study-interactions/4700/6 "2023-06-16T23:54:52Z")

</div>

I want to generate a topology file for a copolymer like PLGA with 1100 atoms to use in charmm forcefield. As you know, cgenff calculates topology files for molecules with less than 384 atoms. How can I solve my problem and which software should I use for my purpose?

---

<div class="post-metadata">

**Author:** ![ddesh5208660](https://avatars.discourse-cdn.com/v4/letter/d/8797f3/32.png) [@ddesh5208660](https://gromacs.bioexcel.eu/u/ddesh5208660)\
**Post date:** [September 26, 2023, 11:09am UTC](https://gromacs.bioexcel.eu/t/how-to-distribute-20-polymers-ligand-randomly-around-the-target-protein-to-study-interactions/4700/7 "2023-09-26T11:09:44Z")

</div>

You have to create RTP file in the required ff as well as take help from other software packages like Swiss param to fill the parameters in RTP file.Make sure don’t blinded accept parametrized values!
