Structural and Dynamical Insight into PPARγ Antagonism: In Silico Study of the Ligand-Receptor Interactions of Non-Covalent Antagonists

The structural and dynamical properties of the peroxisome proliferator-activated receptor γ (PPARγ) nuclear receptor have been broadly studied in its agonist state but little is known about the key features required for the receptor antagonistic activity. Here we report a series of molecular dynamics (MD) simulations in combination with free energy estimation of the recently discovered class of non-covalent PPARγ antagonists. Their binding modes and dynamical behavior are described in details. Two key interactions have been detected within the cavity between helices H3, H11 and the activation helix H12, as well as with H12. The strength of the ligand-amino acid residues interactions has been analyzed in relation to the specificity of the ligand dynamical and antagonistic features. According to our results, the PPARγ activation helix does not undergo dramatic conformational changes, as seen in other nuclear receptors, but rather perturbations that occur through a significant ligand-induced reshaping of the ligand-receptor and the receptor-coactivator binding pockets. The H12 residue Tyr473 and the charge clamp residue Glu471 play a central role for the receptor transformations. Our results also demonstrate that MD can be a helpful tool for the compound phenotype characterization (full agonists, partial agonists or antagonists) when insufficient experimental data are available.     

Potential Activity of Fevicordin-A from Phaleria macrocarpa (Scheff) Boerl. Seeds as Estrogen Receptor Antagonist Based on Cytotoxicity and Molecular Modelling Studies

Fevicordin-A (FevA) isolated from Phaleria macrocarpa (Scheff) Boerl. seeds was evaluated for its potential anticancer activity by in vitro and in silico approaches. Cytotoxicity studies indicated that FevA was selective against cell lines of human breast adenocarcinoma (MCF-7) with an IC₅₀ value of 6.4 μM. At 11.2 μM, FevA resulted in 76.8% cell death of T-47D human breast cancer cell lines. Critical pharmacophore features amongst human Estrogen Receptor-α (hERα) antagonists were conserved in FevA with regard to a hypothesis that they could make notable contributions to its pharmacological activity. The binding stability as well as the dynamic behavior of FevA towards the hERα receptor in agonist and antagonist binding sites were probed using molecular dynamics (MD) simulation approach. Analysis of MD simulation suggested that the tail of FevA was accountable for the repulsion of the C-terminal of Helix-11 (H11) in both agonist and antagonist receptor forms. The flexibility of loop-534 indicated the ability to disrupt the hydrogen bond zipper network between H3 and H11 in hERα. In addition, MM/GBSA calculation from the molecular dynamic simulations also revealed a stronger binding affinity of FevA in antagonistic action as compared to that of agonistic action. Collectively, both the experimental and computational results indicated that FevA has potential as a candidate for an anticancer agent, which is worth promoting for further preclinical evaluation.     

Exploring the Binding Mechanism of PF-07321332 SARS-CoV-2 Protease Inhibitor through Molecular Dynamics and Binding Free Energy Simulations

The novel coronavirus disease, caused by severe acute respiratory coronavirus 2 (SARS-CoV-2), rapidly spreading around the world, poses a major threat to the global public health. Herein, we demonstrated the binding mechanism of PF-07321332, α-ketoamide, lopinavir, and ritonavir to the coronavirus 3-chymotrypsin-like-protease (3CL^pro) by means of docking and molecular dynamic (MD) simulations. The analysis of MD trajectories of 3CL^pro with PF-07321332, α-ketoamide, lopinavir, and ritonavir revealed that 3CL^pro–PF-07321332 and 3CL^pro–α-ketoamide complexes remained stable compared with 3CL^pro–ritonavir and 3CL^pro–lopinavir. Investigating the dynamic behavior of ligand–protein interaction, ligands PF-07321332 and α-ketoamide showed stronger bonding via making interactions with catalytic dyad residues His41–Cys145 of 3CL^pro. Lopinavir and ritonavir were unable to disrupt the catalytic dyad, as illustrated by increased bond length during the MD simulation. To decipher the ligand binding mode and affinity, ligand interactions with SARS-CoV-2 proteases and binding energy were calculated. The binding energy of the bespoke antiviral PF-07321332 clinical candidate was two times higher than that of α-ketoamide and three times than that of lopinavir and ritonavir. Our study elucidated in detail the binding mechanism of the potent PF-07321332 to 3CL^pro along with the low potency of lopinavir and ritonavir due to weak binding affinity demonstrated by the binding energy data. This study will be helpful for the development and optimization of more specific compounds to combat coronavirus disease.     

Co-Incubation with PPARβ/δ Agonists and Antagonists Modeled Using Computational Chemistry: Effect on LPS Induced Inflammatory Markers in Pulmonary Artery

Peroxisome proliferator activated receptor beta/delta (PPARβ/δ) is a nuclear receptor ubiquitously expressed in cells, whose signaling controls inflammation. There are large discrepancies in understanding the complex role of PPARβ/δ in disease, having both anti- and pro-effects on inflammation. After ligand activation, PPARβ/δ regulates genes by two different mechanisms; induction and transrepression, the effects of which are difficult to differentiate directly. We studied the PPARβ/δ-regulation of lipopolysaccharide (LPS) induced inflammation (indicated by release of nitrite and IL-6) of rat pulmonary artery, using different combinations of agonists (GW0742 or L−165402) and antagonists (GSK3787 or GSK0660). LPS induced release of NO and IL-6 is not significantly reduced by incubation with PPARβ/δ ligands (either agonist or antagonist), however, co-incubation with an agonist and antagonist significantly reduces LPS-induced nitrite production and Nos2 mRNA expression. In contrast, incubation with LPS and PPARβ/δ agonists leads to a significant increase in Pdk−4 and Angptl−4 mRNA expression, which is significantly decreased in the presence of PPARβ/δ antagonists. Docking using computational chemistry methods indicates that PPARβ/δ agonists form polar bonds with His287, His413 and Tyr437, while antagonists are more promiscuous about which amino acids they bind to, although they are very prone to bind Thr252 and Asn307. Dual binding in the PPARβ/δ binding pocket indicates the ligands retain similar binding energies, which suggests that co-incubation with both agonist and antagonist does not prevent the specific binding of each other to the large PPARβ/δ binding pocket. To our knowledge, this is the first time that the possibility of binding two ligands simultaneously into the PPARβ/δ binding pocket has been explored. Agonist binding followed by antagonist simultaneously switches the PPARβ/δ mode of action from induction to transrepression, which is linked with an increase in Nos2 mRNA expression and nitrite production.     

Comparative Molecular Dynamics Simulations of Mitogen-Activated Protein Kinase-Activated Protein Kinase 5

The mitogen-activated protein kinase-activated protein kinase MK5 is a substrate of the mitogen-activated protein kinases p38, ERK3 and ERK4. Cell culture and animal studies have demonstrated that MK5 is involved in tumour suppression and promotion, embryogenesis, anxiety, cell motility and cell cycle regulation. In the present study, homology models of MK5 were used for molecular dynamics (MD) simulations of: (1) MK5 alone; (2) MK5 in complex with an inhibitor; and (3) MK5 in complex with the interaction partner p38α. The calculations showed that the inhibitor occupied the active site and disrupted the intramolecular network of amino acids. However, intramolecular interactions consistent with an inactive protein kinase fold were not formed. MD with p38α showed that not only the p38 docking region, but also amino acids in the activation segment, αH helix, P-loop, regulatory phosphorylation region and the C-terminal of MK5 may be involved in forming a very stable MK5-p38α complex, and that p38α binding decreases the residual fluctuation of the MK5 model. Electrostatic Potential Surface (EPS) calculations of MK5 and p38α showed that electrostatic interactions are important for recognition and binding.     

The inverse type II β-turn on D-Trp-Phe, a pharmacophoric motif for MOR agonists

Herein we propose the D ‐Trp‐Phe sequence within an inverse type II β‐turn as a new kind of pharmacophoric motif for μ‐opioid receptor (MOR) cyclopeptide agonists. Initially, we observed that c[Tyr‐D ‐Pro‐D ‐Trp‐Phe‐Gly] ( 4 ), an analogue of endomorphin‐1 (H‐Tyr‐Pro‐Trp‐Phe‐NH₂) lacking the crucial protonatable amino group of Tyr 1, is a MOR agonist with 10^?8 M affinity. Molecular docking analysis suggested that the relevant interactions with the receptor involve D ‐Trp‐Phe. The bioactive conformation of this region was investigated by selected derivatives of 4 designed to adopt an inverse type II β‐turn. These efforts led to c[Tyr‐Gly‐D ‐Trp‐Phe‐Gly] ( 14 ) and to the cyclotetrapeptide c[D ‐Asp‐1‐amide‐β‐Ala‐D ‐Trp‐Phe] ( 15 ), showing improved nanomolar affinity. Both 14 and 15 selectively bind MOR, as they have negligible affinity for the κ‐ and δ‐opioid receptors. Both 14 and 15 behave as partial MOR agonists in functional assays. Conformational and docking analyses confirm the role of the inverse β‐turn in binding. These results indicate that the D ‐Trp‐Phe inverse β‐turn structure can be used for designing non‐endomorphin‐like peptidomimetic opioid agonists in general, characterized by an atypical mechanism of interaction between ligand and receptor.     

ERα36–GPER1 Collaboration Inhibits TLR4/NFκB-Induced Pro-Inflammatory Activity in Breast Cancer Cells

Inflammation is important for the initiation and progression of breast cancer. We have previously reported that in monocytes, estrogen regulates TLR4/NFκB-mediated inflammation via the interaction of the Erα isoform ERα36 with GPER1. We therefore investigated whether a similar mechanism is present in breast cancer epithelial cells, and the effect of ERα36 expression on the classic 66 kD ERα isoform (ERα66) functions. We report that estrogen inhibits LPS-induced NFκB activity and the expression of downstream molecules TNFα and IL-6. In the absence of ERα66, ERα36 and GPER1 are both indispensable for this effect. In the presence of ERα66, ERα36 or GPER1 knock-down partially inhibits NFκB-mediated inflammation. In both cases, ERα36 overexpression enhances the inhibitory effect of estrogen on inflammation. We also verify that ERα36 and GPER1 physically interact, especially after LPS treatment, and that GPER1 interacts directly with NFκB. When both ERα66 and ERα36 are expressed, the latter acts as an inhibitor of ERα66 via its binding to estrogen response elements. We also report that the activation of ERα36 leads to the inhibition of breast cancer cell proliferation. Our data support that ERα36 is an inhibitory estrogen receptor that, in collaboration with GPER1, inhibits NFκB-mediated inflammation and ERα66 actions in breast cancer cells.     

Computational Analysis of the Interactions between the S100B Extracellular Chaperone and Its Amyloid β Peptide Client

S100B is an astrocytic extracellular Ca²⁺-binding protein implicated in Alzheimer’s disease, whose role as a holdase-type chaperone delaying Aβ₄₂ aggregation and toxicity was recently uncovered. Here, we employ computational biology approaches to dissect the structural details and dynamics of the interaction between S100B and Aβ₄₂. Driven by previous structural data, we used the Aβ_25–35 segment, which recapitulates key aspects of S100B activity, as a starting guide for the analysis. We used Haddock to establish a preferred binding mode, which was studied with the full length Aβ using long (1 μs) molecular dynamics (MD) simulations to investigate the structural dynamics and obtain representative interaction complexes. From the analysis, Aβ-Lys28 emerged as a key candidate for stabilizing interactions with the S100B binding cleft, in particular involving a triad composed of Met79, Thr82 and Glu86. Binding constant calculations concluded that coulombic interactions, presumably implicating the Lys28(Aβ)/Glu86(S100B) pair, are very relevant for the holdase-type chaperone activity. To confirm this experimentally, we examined the inhibitory effect of S100B over Aβ aggregation at high ionic strength. In agreement with the computational predictions, we observed that electrostatic perturbation of the Aβ-S100B interaction decreases anti-aggregation activity. Altogether, these findings unveil features relevant in the definition of selectivity of the S100B chaperone, with implications in Alzheimer’s disease.     

Molecular modelling and site-directed mutagenesis of human GALR1 galanin receptor defines determinants of receptor subtype specificity

Human galanin is a 30 amino acid neuropeptide that elicits arange of biological activities by interaction with G protein-coupledreceptors. We have generated a model of the human GALR1 galaninreceptor subtype (hGALR1) based on the alpha carbon maps offrog rhodopsin and investigated the significance of potentialcontact residues suggested by the model using site-directedmutagenesis. Mutation of Phe186 within the second extracellularloop to Ala resulted in a 6-fold decrease in affinity for galanin,representing a change in free energy consistent with hydrophobicinteraction. Our model suggests interaction between Phe186 ofhGALR1 and Ala7 or Leu11 of galanin. Receptor subtype specificitywas investigated by replacement of residues in hGALR1 with thecorresponding residues in hGALR2 and use of the hGALR2-specificligands hGalanin(2–30) and [D-Trp2]hGalanin(1–30).The His267Ile mutant receptor exhibited a pharmacological profilecorresponding to that of hGALR1, suggesting that His267 is notinvolved in a receptor–ligand interaction. The mutationPhe115Ala resulted in a decreased binding affinity for hGalaninand for hGALR2-specific analogues, indicating Phe115 to be ofstructural importance to the ligand binding pocket of hGALR1but not involved in direct ligand interaction. Analysis of Glu271Trpsuggested that Glu271 of hGALR1 interacts with the N-terminusof galanin and that the Trp residue in the corresponding positionin hGALR2 is involved in receptor subtype specificity of binding.Our model supports previous reports of Phe282 of hGALR1 interactingwith Trp2 of galanin and His264 of hGALR1 interacting with Tyr9of galanin.     

Effect of pH on the Aggregation of α-syn12 Dimer in Explicit Water by Replica-Exchange Molecular Dynamics Simulation

The dimeric structure of the N-terminal 12 residues drives the interaction of α-synuclein protein with membranes. Moreover, experimental studies indicated that the aggregation of α-synuclein is faster at low pH than neutral pH. Nevertheless, the effects of different pH on the structural characteristics of the α-syn12 dimer remain poorly understood. We performed 500 ns temperature replica exchange molecular dynamics (T-REMD) simulations of two α-syn12 peptides in explicit solvent. The free energy surfaces contain ten highly populated regions at physiological pH, while there are only three highly populated regions contained at acidic pH. The anti-parallel β-sheet conformations were found as the lowest free energy state. Additionally, these states are nearly flat with a very small barrier which indicates that these states can easily transit between themselves. The dimer undergoes a disorder to order transition from physiological pH to acidic pH and the α-syn12 dimer at acidic pH involves a faster dimerization process. Further, the Lys6–Asp2 contact may prevent the dimerization.     

In Silico Screening of Novel α1-GABAA Receptor PAMs towards Schizophrenia Based on Combined Modeling Studies of Imidazo [1,2-a]-Pyridines

The ionotropic GABA_A receptor (GABA_AR) has been proven to be an important target of atypical antipsychotics. A novel series of imidazo [1,2-a]-pyridine derivatives, as selective positive allosteric modulators (PAMs) of α1-containing GABA_ARs with potent antipsychotic activities, have been reported recently. To better clarify the pharmacological essentiality of these PAMs and explore novel antipsychotics hits, three-dimensional quantitative structure–activity relationships (3D-QSAR), molecular docking, pharmacophore modeling, and molecular dynamics (MD) were performed on 33 imidazo [1,2-a]-pyridines. The constructed 3D-QSAR models exhibited good predictive abilities. The dockings results and MD simulations demonstrated that hydrogen bonds, π–π stackings, and hydrophobic interactions play essential roles in the binding of these novel PAMs in the GABA_AR binding pocket. Four hit compounds (DS01–04) were then screened out by the combination of the constructed models and computations, including the pharmacophore model, Topomer Search, molecular dockings, ADME/T predictions, and MD simulations. The compounds DS03 and DS04, with higher docking scores and better predicted activities, were also found to be relatively stable in the binding pocket by MD simulations. These results might provide a significant theoretical direction or information for the rational design and development of novel α1-GABA_AR PAMs with antipsychotic activities.     

Structure,Dynamics, and Ligand Recognition of Human-Specific CHRFAM7A (Dupα7) Nicotinic Receptor Linked to Neuropsychiatric Disorders

Cholinergic α7 nicotinic receptors encoded by the CHRNA7 gene are ligand-gated ion channels directly related to memory and immunomodulation. Exons 5–7 in CHRNA7 can be duplicated and fused to exons A-E of FAR7a, resulting in a hybrid gene known as CHRFAM7A, unique to humans. Its product, denoted herein as Dupα7, is a truncated subunit where the N-terminal 146 residues of the ligand binding domain of the α7 receptor have been replaced by 27 residues from FAM7. Dupα7 negatively affects the functioning of α7 receptors associated with neurological disorders, including Alzheimer’s diseases and schizophrenia. However, the stoichiometry for the α7 nicotinic receptor containing dupα7 monomers remains unknown. In this work, we developed computational models of all possible combinations of wild-type α7 and dupα7 pentamers and evaluated their stability via atomistic molecular dynamics and coarse-grain simulations. We assessed the effect of dupα7 subunits on the Ca²⁺ conductance using free energy calculations. We showed that receptors comprising of four or more dupα7 subunits are not stable enough to constitute a functional ion channel. We also showed that models with dupα7/α7 interfaces are more stable and are less detrimental for the ion conductance in comparison to dupα7/dupα7 interfaces. Based on these models, we used protein–protein docking to evaluate how such interfaces would interact with an antagonist, α-bungarotoxin, and amyloid Aβ₄₂. Our findings show that the optimal stoichiometry of dupα7/α7 functional pentamers should be no more than three dupα7 monomers, in favour of a dupα7/α7 interface in comparison to a homodimer dupα7/dupα7 interface. We also showed that receptors bearing dupα7 subunits are less sensitive to Aβ₄₂ effects, which may shed light on the translational gap reported for strategies focused on nicotinic receptors in ‘Alzheimer’s disease research.     

Role and Regulation of Mechanotransductive HIF-1α Stabilisation in Periodontal Ligament Fibroblasts

Orthodontic tooth movement (OTM) creates compressive and tensile strain in the periodontal ligament, causing circulation disorders. Hypoxia-inducible factor 1α (HIF-1α) has been shown to be primarily stabilised by compression, but not hypoxia in periodontal ligament fibroblasts (PDLF) during mechanical strain, which are key regulators of OTM. This study aimed to elucidate the role of heparan sulfate integrin interaction and downstream kinase phosphorylation for HIF-1α stabilisation under compressive and tensile strain and to which extent downstream synthesis of VEGF and prostaglandins is HIF-1α-dependent in a model of simulated OTM in PDLF. PDLF were subjected to compressive or tensile strain for 48 h. In various setups HIF-1α was experimentally stabilised (DMOG) or destabilised (YC-1) and mechanotransduction was inhibited by surfen and genistein. We found that HIF-1α was not stabilised by tensile, but rather by compressive strain. HIF-1α stabilisation had an inductive effect on prostaglandin and VEGF synthesis. As expected, HIF-1α destabilisation reduced VEGF expression, whereas prostaglandin synthesis was increased. Inhibition of integrin mechanotransduction via surfen or genistein prevented stabilisation of HIF-1α. A decrease in VEGF expression was observed, but not in prostaglandin synthesis. Stabilisation of HIF-1α via integrin mechanotransduction and downstream phosphorylation of kinases seems to be essential for the induction of VEGF, but not prostaglandin synthesis by PDLF during compressive (but not tensile) orthodontic strain.