Allosteric modulators (AMs) are considered as a perpetual
hotspot
in research for their higher selectivity and various effects on orthosteric
ligands (OL). They are classified in terms of their functionalities
as positive, negative, or silent allosteric modulators (PAM, NAM,
or SAM, respectively). In the present work, 11 pairs of three-dimensional
(3D) structures of receptor–orthosteric ligand and receptor–orthosteric
ligand–allosteric modulator complexes have been collected for
the studies, including three different systems: GPCR, enzyme, and
ion channel. Molecular dynamics (MD) simulations are applied to quantify
the dynamic interactions in both the orthosteric and allosteric binding
pockets and the structural fluctuation of the involved proteins. Our
results showed that MD simulations of moderately large molecules or
peptides undergo insignificant changes compared to crystal structure
results. Furthermore, we also studied the conformational changes of
receptors that bound with PAM and NAM, as well as the different allosteric
binding sites in a receptor. There should be no preference for the
position of the allosteric binding pocket after comparing the allosteric
binding pockets of these three systems. Finally, we aligned four distinct
β2 adrenoceptor structures and three N-methyl-d-aspartate receptor (NMDAR) structures to investigate conformational
changes. In the β2 adrenoceptor systems, the aligned results
revealed that transmembrane (TM) helices 1, 5, and 6 gradually increased
outward movement from an enhanced inactive state to an improved active
state. TM6 endured the most significant conformational changes (around
11 Å). For NMDAR, the bottom section of NMDAR’s ligand-binding
domain (LBD) experienced an upward and outward shift during the gradually
activating process. In conclusion, our research provides insight into
receptor–orthosteric ligand–allosteric modulator studies
and the design and development of allosteric modulator drugs using
MD simulation.