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.

No abstract

Multidrug resistance protein-4 (MRP4) belongs to the ABC transporter superfamily and promotes the transport of xenobiotics including drugs. A non-synonymous single nucleotide polymorphisms (nsSNPs) in the ABCC4 gene can promote changes in the structure and function of MRP4. In this work, the interaction of certain endogen substrates, drug substrates, and inhibitors with wild type-MRP4 (WT-MRP4) and its variants G187W and Y556C were studied to determine differences in the intermolecular interactions and affinity related to SNPs using protein threading modeling, molecular docking, all-atom, coarse grained, and umbrella sampling molecular dynamics simulations (AA-MDS and CG-MDS, respectively). The results showed that the three MRP4 structures had significantly different conformations at given sites, leading to differences in the docking scores (DS) and binding sites of three different groups of molecules. Folic acid (FA) had the highest variation in DS on G187W concerning WT-MRP4. WT-MRP4, G187W, Y556C, and FA had different conformations through 25 ns AA-MD. Umbrella sampling simulations indicated that the Y556C-FA complex was the most stable one with or without ATP. In Y556C, the cyclic adenosine monophosphate (cAMP) and ceefourin-1 binding sites are located out of the entrance of the inner cavity, which suggests that both cAMP and ceefourin-1 may not be transported. The binding site for cAMP and ceefourin-1 is quite similar and the affinity (binding energy) of ceefourin-1 to WT-MRP4, G187W, and Y556C is greater than the affinity of cAMP, which may suggest that ceefourin-1 works as a competitive inhibitor. In conclusion, the nsSNPs G187W and Y556C lead to changes in protein conformation, which modifies the ligand binding site, DS, and binding energy.

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.