Hi, all:<br>I'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'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>