Investigating cyclic peptides inhibiting CD2-CD58 interactions through molecular dynamics and molecular docking methods

Virtual Screening and Binding Analysis of Potential CD58 Inhibitors in Colorectal Cancer (CRC)

Human cell surface receptor CD58, also known as lymphocyte function-associated antigen 3 (LFA-3), plays a critical role in the early stages of immune response through interacting with CD2. Recent research identified CD58 as a surface marker of colorectal cancer (CRC), which can upregulate the Wnt pathway and promote self-renewal of colorectal tumor-initiating cells (CT-ICs) by degradation of Dickkopf 3. In addition, it was also shown that knockdown of CD58 significantly impaired tumor growth. In this study, we developed a structure-based virtual screening pipeline using Autodock Vina and binding analysis and identified a group of small molecular compounds having the potential to bind with CD58. Five of them significantly inhibited the growth of the SW620 cell line in the following in vitro studies. Their proposed binding models were further verified by molecular dynamics (MD) simulations, and some pharmaceutically relevant chemical and physical properties were predicted. The hits described in this work may be considered interesting leads or structures for the development of new and more efficient CD58 inhibitors.

CD58 Immunobiology at a Glance

The glycoprotein CD58, also known as lymphocyte-function antigen 3 (LFA-3), is a costimulatory receptor distributed on a broad range of human tissue cells. Its natural ligand CD2 is primarily expressed on the surface of T/NK cells. The CD2-CD58 interaction is an important component of the immunological synapse (IS) that induces activation and proliferation of T/NK cells and triggers a series of intracellular signaling in T/NK cells and target cells, respectively, in addition to promoting cell adhesion and recognition. Furthermore, a soluble form of CD58 (sCD58) is also present in cellular supernatant in vitro and in local tissues in vivo. The sCD58 is involved in T/NK cell-mediated immune responses as an immunosuppressive factor by affecting CD2-CD58 interaction. Altered accumulation of sCD58 may lead to immunosuppression of T/NK cells in the tumor microenvironment, allowing sCD58 as a novel immunotherapeutic target. Recently, the crucial roles of costimulatory molecule CD58 in immunomodulation seem to be reattracting the interests of investigators. In particular, the CD2-CD58 interaction is involved in the regulation of antiviral responses, inflammatory responses in autoimmune diseases, immune rejection of transplantation, and immune evasion of tumor cells. In this review, we provide a comprehensive summary of CD58 immunobiology.

Structure-based identification of inhibitors disrupting the CD2-CD58 interactions

The immune system has very intricate mechanisms of fighting against the invading infections which are accomplished by a sequential event of molecular interactions in the body. One of the crucial phenomena in this process is the recognition of T-cells by the antigen-presenting cells (APCs), which is initiated by the rapid interaction between both cell surface receptors, i.e., CD2 located on T-cells and CD58 located on APCs. Under various pathological conditions, which involve undesired immune response, inhibiting the CD2-CD58 interactions becomes a therapeutically relevant opportunity. Herein we present an extensive work to identify novel inhibiting agents of the CD2-CD58 interactions. Classical molecular dynamics (MD) simulations of the CD2-CD58 complex highlighted a series of crucial CD58 residues responsible for the interactions with CD2. Based on such results, a pharmacophore map, complementary to the CD2-binding site of CD58, was created and employed for virtual screening of ~ 300,000 available compounds. On the ~ 6000 compounds filtered from pharmacophore mapping, ADME screening leads to ~ 350 molecules. Molecular docking was then performed on these molecules, and fifteen compounds emerged with significant binding energy (< - 50 kcal/mol) for CD58. Finally, short MD simulations were performed in triplicate on each complex (i) to provide a microscopic view of the ligand binding and (ii) to rule out possibly weak binders of CD58 from the identified hits. At last, we suggest eight compounds for in vitro testing that were identified as promising hits to bind CD58 with a high binding affinity.

Consensus Analyses in Molecular Docking Studies Applied to Medicinal Chemistry

The increasing number of computational studies in medicinal chemistry involving molecular docking has put the technique forward as promising in Computer-Aided Drug Design. Considering the main method in the virtual screening based on the structure, consensus analysis of docking has been applied in several studies to overcome limitations of algorithms of different programs and mainly to increase the reliability of the results and reduce the number of false positives. However, some consensus scoring strategies are difficult to apply and, in some cases, are not reliable due to the small number of datasets tested. Thus, for such a methodology to be successful, it is necessary to understand why, when and how to use consensus docking. Therefore, the present study aims to present different approaches to docking consensus, applications, and several scoring strategies that have been successful and can be applied in future studies.

From quantum to subcellular scales: multi-scale simulation approaches and the SIRAH force field

Modern molecular and cellular biology profits from astonishing resolution structural methods, currently even reaching the whole cell level. This is encompassed by the development of computational methods providing a deep view into the structure and dynamics of molecular processes happening at very different scales in time and space. Linking such scales is of paramount importance when aiming at far-reaching biological questions. Computational methods at the interface between classical and coarse-grained resolutions are gaining momentum with several research groups dedicating important efforts to their development and tuning. An overview of such methods is addressed herein, with special emphasis on the SIRAH force field for coarse-grained and multi-scale simulations. Moreover, we provide proof of concept calculations on the implementation of a multi-scale simulation scheme including quantum calculations on a classical fine-grained/coarse-grained representation of double-stranded DNA. This opens the possibility to include the effect of large conformational fluctuations in chromatin segments on, for instance, the reactivity of particular base pairs within the same simulation framework.