Different environments for a realistic simulation of GPCRs-application to the M2 muscarinic receptor

A model of the human M(2) muscarinic receptor was taken as an example for a class A G-protein coupled receptor to explore the influence of different environments in a molecular dynamics simulation (MDS) on the protein structure. The most commonly used environment is the vacuum, although it is very unnatural for a transmembrane protein. As an alternative a membrane-like system, consisting of a lipophilic central layer and two aqueous flanking layers, was tested. The most realistic system that can be applied is a phospholipid bilayer with a surrounding physiological sodium chloride solution. From all systems good protein structures were received, nevertheless clear differences between the systems were detected in the structural comparison of the models. Subsequently it was analyzed whether the observed structural differences influence ligand binding. For this purpose the antagonist (S)-scopolamine was docked into the binding cavity, which is well known by many reported single and multiple point mutations. As expected from the observed structural variations triggered by the type of environment employed in MDS, also differences in the binding mode of (S)-scopolamine were detected, all contacts, however, which are known to be important were found.