Hi, all:<br>I&#39;m trying to embed cholesterol molecules into a DOPC bilayer in a coarse grained model. <br>The force field and coordinates for a pure DOPC bilayer and for cholesterol were obtained on the MARTINI&#39;s site. I randomly replace a DOPC molecule with a cholesterol molecule and  do the energy minimization. After embedded a certain amount of cholesterols, I run the simulation by coupling DOPC and CHOL together to a thermostat. However, the cholesterols turn to be very unstable and flip-flop a lot between the monolayers, which I think is not expected.
<br>  Any suggestions will be greatly appreciated.Thank you.<br>dj<br>