In this study, we investigated the stability of the fully activated conformation of the orexin receptor 2 (OX2R) embedded in a pure POPC bilayer using MD simulations. Various thermodynamic ensembles (i.e., NPT, NVT, NVE, NPAT, μVT, and NPγT) were employed to explore the dynamical heterogeneity of the system in a comprehensive way. In addition, informational similarity metrics (e.g., Jensen-Shannon divergence) as well as Markov state modeling approaches were utilized to elucidate the receptor kinetics. Special attention was paid to assessing surface tension within the simulation box, particularly under NPγT conditions, where 21 nominal surface tension constants were evaluated. Our findings suggest that traditional thermodynamic ensembles such as NPT may not adequately control physical properties of the POPC membrane, impacting the plausibility of the OX2R model. In general, the performed study underscores the importance of employing the NPγT ensemble for computational investigations of membrane-embedded receptors, as it effectively maintains zero surface tension in the simulated system. These results offer valuable insights for future research aimed at understanding receptor dynamics and designing targeted therapeutics.

Biomedical Research Centre University Hospital Hradec Kralove Sokolska 581 500 03 Hradec Kralove Czech Republic

Department of Epidemiology 2nd Faculty of Medicine Charles University 5 Uvalu 84 150 06 Prague Czech Republic

References provided by Crossref.org