Molecular Dynamics Simulation of Chemokine Receptors in Lipid Bilayer: A Case Study on C-C Chemokine Receptor Type 2

Cholesterol Dependence of the Conformational Changes in Metabotropic Glutamate Receptor 1

Metabotropic glutamate receptors (mGluRs) are class C G protein-coupled receptors that function as obligate dimers in regulating neurotransmission and synaptic plasticity in the central nervous system. The mGluR1 subtype has been shown to be modulated by the membrane lipid environment, particularly cholesterol, though the molecular mechanisms remain elusive. In this study, we employed all-atom molecular dynamics simulations to investigate the effects of cholesterol on the conformational dynamics of the mGluR1 seven-transmembrane (7TM) domain in an inactive state model. Simulations were performed with three different cholesterol concentrations (0%, 10%, and 25%) in a palmitoyl-oleoyl phosphatidylcholine (POPC) lipid bilayer system. Our results demonstrate that cholesterol induces conformational changes in the mGluR1 dimer more significantly than in the individual protomers. Notably, cholesterol modulates the dynamics and conformations of the TM1 and TM2 helices at the dimer interface. Interestingly, an intermediate cholesterol concentration of 10% elicits more pronounced conformational changes compared to both cholesterol-depleted (0%) and cholesterol-enriched (25%) systems. Specific electrostatic interaction unique to the 10% cholesterol system further corroborate these conformational differences. Given the high sequence conservation of the 7TM domains across mGluR subtypes, the cholesterol-dependent effects observed in mGluR1 are likely applicable to other members of this receptor family. Our findings provide atomistic insights into how cholesterol modulates the conformational landscape of mGluRs, which could impact their function and signaling mechanisms.

A unique ligand-steered strategy for CC chemokine receptor 2 homology modeling to facilitate structure-based virtual screening

CC chemokine receptor 2 (CCR2) antagonists that disrupt CCR2/MCP-1 interaction are expected to treat a variety of inflammatory and autoimmune diseases. The lack of CCR2 crystal structure limits the application of structure-based drug design (SBDD) to this target. Although a few three-dimensional theoretical models have been reported, their accuracy remains to be improved in terms of templates and modeling approaches. In this study, we developed a unique ligand-steered strategy for CCR2 homology modeling. It starts with an initial model based on the X-ray structure of the closest homolog so far, that is, CXCR4. Then, it uses Elastic Network Normal Mode Analysis (EN-NMA) and flexible docking (FD) by AutoDock Vina software to generate ligand-induced fit models. It selects optimal model(s) as well as scoring function(s) via extensive evaluation of model performance based on a unique benchmarking set constructed by our in-house tool, that is, MUBD-DecoyMaker. The model of 81_04 presents the optimal enrichment when combined with the scoring function of PMF04, and the proposed binding mode between CCR2 and Teijin lead by this model complies with the reported mutagenesis data. To highlight the advantage of our strategy, we compared it with the only reported ligand-steered strategy for CCR2 homology modeling, that is, Discovery Studio/Ligand Minimization. Lastly, we performed prospective virtual screening based on 81_04 and CCR2 antagonist bioassay. The identification of two hit compounds, that is, E859-1281 and MolPort-007-767-945, validated the efficacy of our model and the ligand-steered strategy.

A novel amino cellulose derivative using ATRP method: Preparation, characterization, and investigation of its antibacterial activity

In this study, we prepared a novel amino cellulose derivative (benzyl cellulose-g-poly [2-(N,N-Dimethylamino)ethyl methacrylate]) via a homogeneous ATRP method. The successful synthesis of the novel amino cellulose was confirmed by FT-IR and ¹H NMR. This study addressed the different characteristics of the prepared polymer including the thermal stability, solubility, and X-ray diffraction pattern. The antibacterial activity of the synthesized cellulose derivative was investigated using the diffusion disk method against both gram-negative (Escherichia coli, Salmonella enterica) and gram-positive (Staphylococcus aureus, Bacillus subtilis) bacteria. Based on the inhibition zone, it was confirmed that the prepared benzyl cellulose-g-PDMAEMA possesses acceptable antibacterial activity against Escherichia coli, Salmonella enterica, and Staphylococcus aureus while Bacillus subtilis is resistant to the prepared polymer. Also according to the inhibition zone, it was shown that benzyl cellulose-g-PDMAEMA has more impact on E. coli and Salmonella enterica than Staphylococcus aureus. Molecular dynamics simulation was also used to study the interaction of the synthesized cellulose derivative with a model membrane which presented atomistic details of the polymer-lipid interactions. According to the results obtained from the molecular dynamics simulation, the polymer was able to destabilize the structure of the membrane and clearly express its signs of degradation.

CCL2/CCR2 Chemokine System in Embryonic Hypothalamus: Involvement in Sexually Dimorphic Stimulatory Effects of Prenatal Ethanol Exposure on Peptide-Expressing Neurons

Maternal consumption of ethanol during pregnancy is known to increase the offspring's risk for developing alcohol use disorders and associated behavioral disturbances. Studies in adolescent and adult animals suggest the involvement of neuroimmune and neurochemical systems in the brain that control these behaviors. To understand the origin of these effects during early developmental stages, we examined in the embryo and neonate the effects of maternal intraoral administration of ethanol (2 g/kg/day) from embryonic day 10 (E10) to E15 on the inflammatory chemokine C-C motif ligand 2 (CCL2) and its receptor CCR2 in a specific, dense population of neurons in the lateral hypothalamus (LH), where they are closely related to an orexigenic neuropeptide, melanin-concentrating hormone (MCH), known to promote ethanol consumption and related behaviors. We found that prenatal ethanol exposure increases the expression and density of CCL2 and CCR2 cells along with MCH neurons in the LH and the colocalization of CCL2 with MCH. We also discovered that these effects are sexually dimorphic, consistently stronger in female embryos, and are blocked by maternal administration of a CCL2 antibody (1 and 5 µg/day, i.p., E10-E15) that neutralizes endogenous CCL2 and of a CCR2 antagonist INCB3344 (1 mg/day, i.p., E10-E15) that blocks CCL2's main receptor. These results, which in the embryo anatomically and functionally link the CCL2/CCR2 system to MCH neurons in the LH, suggest an important role for this neuroimmune system in mediating ethanol's sexually dimorphic, stimulatory effect on MCH neurons that may promote higher level of alcohol consumption described in females.

Structural basis for ligand modulation of the CCR2 conformational landscape

CC chemokine receptor 2 (CCR2) is a part of the chemokine receptor family, an important class of therapeutic targets. These class A G-protein coupled receptors (GPCRs) are involved in mammalian signaling pathways and control cell migration toward endogenous CC chemokine ligands, named for the adjacent cysteine motif on their N terminus. Chemokine receptors and their associated ligands are involved in a wide range of diseases and thus have become important drug targets. CCR2, in particular, promotes the metastasis of cancer cells and is also implicated in autoimmunity-driven type-1 diabetes, diabetic nephropathy, multiple sclerosis, asthma, atherosclerosis, neuropathic pain, and rheumatoid arthritis. Although promising, CCR2 antagonists have been largely unsuccessful to date. Here, we investigate the effect of an orthosteric and an allosteric antagonist on CCR2 dynamics by coupling long-timescale molecular dynamics simulations with Markov-state model theory. We find that the antagonists shift CCR2 into several stable inactive conformations that are distinct from the crystal structure conformation and disrupt a continuous internal water and sodium ion pathway, preventing transitions to an active-like state. Several metastable conformations present a cryptic drug-binding pocket near the allosteric site that may be amenable to targeting with small molecules. Without antagonists, the apo dynamics reveal intermediate conformations along the activation pathway that provide insight into the basal dynamics of CCR2 and may also be useful for future drug design.

Shedding light on the structural properties of lipid bilayers using molecular dynamics simulation: a review study

In recent years, a massive increase has been observed in the number of published articles describing accurate and reliable molecular dynamics simulations of lipid bilayers. This is due to several reasons, including the development of fast and efficient methods for treating long-range electrostatic interactions, significant progress in computer hardware, progress in atomistic simulation algorithms and the development of well-validated empirical molecular mechanical force fields. Although molecular dynamics is an effective approach for investigating different aspects of lipid bilayers, to the best of our knowledge, there is no review in the literature that explains the different analyses that can be carried out with membrane simulation. This review gives an overview about the some of the most important possible analyses, technical challenges, and existing protocols that can be performed on the biological membrane by molecular dynamics simulation. The reviewed analyses include the degree of membrane disruption, average area per lipid, probability distributions for the area per lipid molecule, membrane thickness, membrane area compressibility, lateral diffusion, rotational diffusion, order parameters, head group tilt, electron density profile, mass density profile, electrostatic potential profile, ordering of vicinity waters, number of hydrogen bonds, and radial distribution function.

Conformational analysis of the human chemokine receptor CXCR3

In the last years, some studies showed the patho-genetic role of CXCR3 bound to its ligands in many human inflammatory diseases and cancers. Thus, the blockage of the CXCR3 interaction site to its ligands is seen as a possible therapeutic target for the treatment of cancer. The presence of flexible regions in the chemokine receptors determines their capability to develop specific mechanisms of action. We have recently focused on the features of the N-terminal region of human CXCR3 free in solution, where we demonstrate the presence of numerous conformational ensembles, dynamically stabilized by H-bonds. Since up to now no structure was experimentally determined for CXCR3, we decided to approach the study of its conformational behavior by molecular dynamics simulations, in a lipid bilayer, surrounded of water, at neutral pH and 300K. Furthermore, we modeled the CXCR3/CXCL11 complex, where CXCL11 is one of its natural ligands. The aim of this work is to have a vision as realistic as possible in dynamic terms of the biological mechanism that drives the search for the ligand, its interaction and the formation of a stable complex between CXCR3 and CXCL11. Overall, our approach has been able to describe the structural events which dynamically characterize the molecular mechanisms involved in the binding of CXCR3 to CXCL11 and the critical role exerted by its N-terminal region in "hunting" and capturing the ligand.

Molecular insight into the Grandivitin- matrix metalloproteinase 9 interactions

Grandivitin (GRA), a natural coumarin, can inhibit Matrix metalloproteinase 9 (MMP9). Binding characteristics are therefore of interest for pharmacodynamics of GRA and coumarin derivatives. A combination of spectroscopic methods and molecular modeling techniques was used to characterize interaction of GRA with MMP9. Fluorescence spectroscopy showed that GRA could quench the MMP9 fluorescence spectra. Changes in the UV-Vis and FT-IR spectra were observed upon ligand binding along with a significant degree of tryptophan fluorescence quenching on complex formation. Fluorescence studies showed that GRA has an ability to quench the intrinsic fluorescence of MMP9. Molecular modeling analysis showed that GRA to be bound in the large hydrophobic cavity of MMP9. Further investigation of the binding site of GRA within the MMP9 molecule suggested that hydrophobic contacts, hydrogen bond formation and electrostatic interactions account for the binding of GRA. According molecular dynamics (MD) simulation results the ligand can interact with the protein, with affecting the secondary structure of MMP9 and with a modification of its tertiary structure. The biological significance of this work is evident because MMP9 serves as a potential target protein for anticancer agents. The binding study of GRA with MMP9 is of great importance in pharmacy, pharmacology and biochemistry. This work can provide some key data to clinical research and supply the theoretical basis for the new drug candidate designing.