scholarly journals Behavior of Chemokine Receptor 6 (CXCR6) in Complex with CXCL16 Soluble form Chemokine by Molecular Dynamic Simulations: General Protein‒Ligand Interaction Model and 3D-QSAR Studies of Synthetic Antagonists

Life ◽  
2021 ◽  
Vol 11 (4) ◽  
pp. 346
Author(s):  
Giovanny Aguilera-Durán ◽  
Antonio Romo-Mancillas
Keyword(s):  
Molecular Dynamics ◽  
3D Qsar ◽  
Md Simulations ◽  
Active State ◽  
Activation Mechanism ◽  
Soluble Form ◽  
Ligand Interaction ◽  
Protein Ligand Interaction ◽  
Polar Interactions ◽  
General Protein

The CXCR6‒CXCL16 axis is involved in several pathological processes, and its overexpression has been detected in different types of cancer, such as prostate, breast, ovary, and lung cancer, along with schwannomas, in which it promotes invasion and metastasis. Moreover, this axis is involved in atherosclerosis, type 1 diabetes, primary immune thrombocytopenia, vitiligo, and other autoimmune diseases, in which it is responsible for the infiltration of different immune system cells. The 3D structure of CXCR6 and CXCL16 has not been experimentally resolved; therefore, homology modeling and molecular dynamics simulations could be useful for the study of this signaling axis. In this work, a homology model of CXCR6 and a soluble form of CXCL16 (CXCR6‒CXCL16s) are reported to study the interactions between CXCR6 and CXCL16s through coarse-grained molecular dynamics (CG-MD) simulations. CG-MD simulations showed the two activation steps of CXCR6 through a decrease in the distance between the chemokine and the transmembrane region (TM) of CXCR6 and transmembrane rotational changes and polar interactions between transmembrane segments. The polar interactions between TM3, TM5, and TM6 are fundamental to functional conformation and the meta-active state of CXCR6. The interactions between D77-R280 and T243-TM7 could be related to the functional conformation of CXCR6; alternatively, the interaction between Q195-Q244 and N248 could be related to an inactive state due to the loss of this interaction, and an arginine cage broken in the presence of CXCL16s allows the meta-active state of CXCR6. A general protein‒ligand interaction supports the relevance of TM3‒TM5‒TM6 interactions, presenting three relevant pharmacophoric features: HAc (H-bond acceptor), HDn (H-bond donator), and Hph (hydrophobic), distributed around the space between extracellular loops (ECLs) and TMs. The HDn feature is close to TM3 and TM6; likewise, the HAc and Hph features are close to ECL1 and ECL2 and could block the rotation and interactions between TM3‒TM6 and the interactions of CXCL16s with the ECLs. Tridimensional quantitative structure-activity relationships (3D-QSAR) models show that the positive steric (VdW) and electrostatic fields coincide with the steric and positive electrostatic region of the exo-azabicyclo[3.3.1]nonane scaffold in the best pIC50 ligands. This substructure is close to the E274 residue and therefore relevant to the activity of CXCR6. These data could help with the design of new molecules that inhibit chemokine binding or antagonize the receptor based on the activation mechanism of CXCR6 and provoke a decrease in chemotaxis caused by the CXCR6‒CXCL16 axis.

scholarly journals Computational Study of C-X-C Chemokine Receptor (CXCR)3 Binding with Its Natural Agonists Chemokine (C-X-C Motif) Ligand (CXCL)9, 10 and 11 and with Synthetic Antagonists: Insights of Receptor Activation towards Drug Design for Vitiligo

Molecules ◽  
2020 ◽  
Vol 25 (19) ◽  
pp. 4413
Author(s):  
Giovanny Aguilera-Durán ◽  
Antonio Romo-Mancillas
Keyword(s):  
Molecular Dynamics ◽  
Computational Study ◽  
Md Simulations ◽  
Receptor Activation ◽  
Activation Mechanism ◽  
Coarse Grained ◽  
Dimensional Structure ◽  
Cytotoxic Lymphocytes ◽  
Ligand Interaction ◽  
Α Subunit

Vitiligo is a hypopigmentary skin pathology resulting from the death of melanocytes due to the activity of CD8+ cytotoxic lymphocytes and overexpression of chemokines. These include CXCL9, CXCL10, and CXCL11 and its receptor CXCR3, both in peripheral cells of the immune system and in the skin of patients diagnosed with vitiligo. The three-dimensional structure of CXCR3 and CXCL9 has not been reported experimentally; thus, homology modeling and molecular dynamics could be useful for the study of this chemotaxis-promoter axis. In this work, a homology model of CXCR3 and CXCL9 and the structure of the CXCR3/Gαi/0βγ complex with post-translational modifications of CXCR3 are reported for the study of the interaction of chemokines with CXCR3 through all-atom (AA-MD) and coarse-grained molecular dynamics (CG-MD) simulations. AA-MD and CG-MD simulations showed the first activation step of the CXCR3 receptor with all chemokines and the second activation step in the CXCR3-CXCL10 complex through a decrease in the distance between the chemokine and the transmembrane region of CXCR3 and the separation of the βγ complex from the α subunit in the G-protein. Additionally, a general protein–ligand interaction model was calculated, based on known antagonists binding to CXCR3. These results contribute to understanding the activation mechanism of CXCR3 and the design of new molecules that inhibit chemokine binding or antagonize the receptor, provoking a decrease of chemotaxis caused by the CXCR3/chemokines axis.

scholarly journals Molecular Dynamics Scoring of Protein–Peptide Models Derived from Coarse-Grained Docking

Molecules ◽  
2021 ◽  
Vol 26 (11) ◽  
pp. 3293
Author(s):  
Mateusz Zalewski ◽  
Sebastian Kmiecik ◽  
Michał Koliński
Keyword(s):  
Molecular Dynamics ◽  
Geometry Optimization ◽  
Computational Prediction ◽  
Md Simulations ◽  
Coarse Grained ◽  
Ligand Interaction ◽  
Vast Number ◽  
Docking Simulations ◽  
Peptide Models ◽  
Peptide Complexes

One of the major challenges in the computational prediction of protein–peptide complexes is the scoring of predicted models. Usually, it is very difficult to find the most accurate solutions out of the vast number of sometimes very different and potentially plausible predictions. In this work, we tested the protocol for Molecular Dynamics (MD)-based scoring of protein–peptide complex models obtained from coarse-grained (CG) docking simulations. In the first step of the scoring procedure, all models generated by CABS-dock were reconstructed starting from their original C-alpha trace representations to all-atom (AA) structures. The second step included geometry optimization of the reconstructed complexes followed by model scoring based on receptor–ligand interaction energy estimated from short MD simulations in explicit water. We used two well-known AA MD force fields, CHARMM and AMBER, and a CG MARTINI force field. Scoring results for 66 different protein–peptide complexes show that the proposed MD-based scoring approach can be used to identify protein–peptide models of high accuracy. The results also indicate that the scoring accuracy may be significantly affected by the quality of the reconstructed protein receptor structures.

scholarly journals An Investigation of Small GTPases in relation to Liver Tumorigenesis Using Traditional Chinese Medicine

2014 ◽  
Vol 2014 ◽  
pp. 1-11 ◽  
Author(s):  
Tzu-Chieh Hung ◽  
Wen-Yuan Lee ◽  
Kuen-Bao Chen ◽  
Yueh-Chiu Chan ◽  
Calvin Yu-Chian Chen
Keyword(s):  
Molecular Dynamics ◽  
Chinese Medicine ◽  
Traditional Chinese Medicine ◽  
Hydrophobic Interactions ◽  
Small Gtpases ◽  
Ligand Interaction ◽  
Protein Ligand Interactions ◽  
Ligand Interactions ◽  
Protein Ligand Interaction ◽  
Better Than

Recently, an important topic of liver tumorigenesis had been published in 2013. In this report, Ras and Rho had defined the relation of liver tumorigenesis. The traditional Chinese medicine (TCM) database has been screened for molecular compounds by simulating molecular docking and molecular dynamics to regulate Ras and liver tumorigenesis. Saussureamine C, S-allylmercaptocysteine, and Tryptophan are selected based on the highest docking score than other TCM compounds. The molecular dynamics are helpful in the analysis and detection of protein-ligand interactions. Based on the docking poses, hydrophobic interactions, and hydrogen bond variations, this research surmises are the main regions of important amino acids in Ras. In addition to the detection of TCM compound efficacy, we suggest Saussureamine C is better than the others for protein-ligand interaction.

scholarly journals 1J1500 An analysis of protein-ligand interaction by multicanonical molecular dynamics

2002 ◽  
Vol 42 (supplement2) ◽  
pp. S57
Author(s):  
T. Terada ◽  
Y. Matsuo ◽  
A. Kidera
Keyword(s):  
Molecular Dynamics ◽  
Ligand Interaction ◽  
Protein Ligand Interaction

scholarly journals Saponins as a Potential PFN2 Inhibitor: In silico Approach

2022 ◽  
Author(s):  
Berly Cárdenas-Pillco ◽  
Lizbeth Campos-Olazaval ◽  
Patricia López ◽  
Jorge Alberto Aguilar-Pineda ◽  
Pamela Lily Gamero-Begazo ◽  
...  
Keyword(s):  
Colorectal Cancer ◽  
Binding Affinity ◽  
Oleanolic Acid ◽  
Md Simulations ◽  
New Drugs ◽  
Actin Binding ◽  
Ligand Interaction ◽  
Site Analysis ◽  
Colorectal Cancer Cells ◽  
Protein Ligand Interaction

Abstract Colorectal cancer (CRC) disease has a high mortality rate and has recently involved human profilin II (Pfn2), an actin-binding protein, as a promoter of its invasiveness and progression. This work evaluated the binding affinity of oleanolic acid saponin over Pfn2 and its structural stability. QM and MM techniques were applied to perform geometrical optimization and calculation of the reactive sites from oleanolic acid, whereas molecular docking and MD simulations for protein-ligand interaction under physiological conditions. Oleanolic acid saponin showed a high binding affinity to the Pfn2 PLPbinding site. Analysis of the protein-ligand structure suggests saponin as a molecule with high potential for developing new drugs against Pfn2 in colorectal cancer cells.

scholarly journals 3D-QSAR, Molecular Docking, and MD Simulations of Anthraquinone Derivatives as PGAM1 Inhibitors

2021 ◽  
Vol 12 ◽  
Author(s):  
Yuwei Wang ◽  
Yifan Guo ◽  
Shaojia Qiang ◽  
Ruyi Jin ◽  
Zhi Li ◽  
...  
Keyword(s):  
Molecular Dynamics ◽  
Molecular Docking ◽  
Rational Design ◽  
Binding Free Energy ◽  
3D Qsar ◽  
Md Simulations ◽  
Structure Activity Relationships ◽  
Structure Activity ◽  
Wide Range ◽  
Anthraquinone Derivatives

PGAM1 is overexpressed in a wide range of cancers, thereby promoting cancer cell proliferation and tumor growth, so it is gradually becoming an attractive target. Recently, a series of inhibitors with various structures targeting PGAM1 have been reported, particularly anthraquinone derivatives. In present study, the structure–activity relationships and binding mode of a series of anthraquinone derivatives were probed using three-dimensional quantitative structure–activity relationships (3D-QSAR), molecular docking, and molecular dynamics (MD) simulations. Comparative molecular field analysis (CoMFA, r2 = 0.97, q2 = 0.81) and comparative molecular similarity indices analysis (CoMSIA, r2 = 0.96, q2 = 0.82) techniques were performed to produce 3D-QSAR models, which demonstrated satisfactory results, especially for the good predictive abilities. In addition, molecular dynamics (MD) simulations technology was employed to understand the key residues and the dominated interaction between PGAM1 and inhibitors. The decomposition of binding free energy indicated that the residues of F22, K100, V112, W115, and R116 play a vital role during the ligand binding process. The hydrogen bond analysis showed that R90, W115, and R116 form stable hydrogen bonds with PGAM1 inhibitors. Based on the above results, 7 anthraquinone compounds were designed and exhibited the expected predictive activity. The study explored the structure–activity relationships of anthraquinone compounds through 3D-QSAR and molecular dynamics simulations and provided theoretical guidance for the rational design of new anthraquinone derivatives as PGAM1 inhibitors.

Sign in / Sign up

Export Citation Format

Share Document