Exploring key features of selectivity in somatostatin receptors through molecular dynamics simulations

Somatostatin receptors (SSTRs) are widely distributed throughout the human body and play crucial roles in various physiological processes. They are recognized as key targets for both radiotherapy and radiodiagnosis due to their overexpression in several cancer types. However, the discovery and design of selective drugs for each of the five isoforms have been significantly hindered by the lack of complete structural information. In this study, we conducted a systematic computational analysis of all five SSTRs in complex with the endogenous ligand somatostatin to elucidate their structural and dynamic features. We thoroughly characterized each isoform using available experimental structures for SSTR2 and SSTR4, as well as AlphaFold2 models for SSTR1, SSTR3, and SSTR5. By performing multi-copy μs-long molecular dynamics simulations, we examined the differences and similarities in dynamical behaviour and somatostatin binding among all SSTRs. Our analysis focused on understanding the opening and closing movements of the extracellular loop 2, which are crucial for ligand binding and recognition. Interestingly, we observed a unique conformation of somatostatin within the binding pocket of SSTR5 in which the loop can partially close, as compared to the other isoforms. Fingerprint analyses provided distinct interaction patterns of somatostatin with all receptors, thus enabling precise guidelines for the discovery and development of more selective somatostatin-based pharmaceuticals tailored for precision medicine therapies.