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.     

Molecular Dynamics Simulations Based on 1-Phenyl-4-Benzoyl-1-Hydro-Triazole ERRα Inverse Agonists

Estrogen-related receptor α (ERRα), which is overexpressed in a variety of cancers has been considered as an effective target for anticancer therapy. ERRα inverse agonists have been proven to effectively inhibit the migration and invasion of cancer cells. As few crystalline complexes have been reported, molecular dynamics (MD) simulations were carried out in this study to deepen the understanding of the interaction mechanism between inverse agonists and ERRα. The binding free energy was analyzed by the MM-GBSA method. The results show that the total binding free energy was positively correlated with the biological activity of an inverse agonist. The interaction of the inverse agonist with the hydrophobic interlayer composed of Phe328 and Phe495 had an important impact on the biological activity of inverse agonists, which was confirmed by the decomposition of energy on residues. As Glu331 flipped and formed a hydrogen bond with Arg372 in the MD simulation process, the formation of hydrogen bond interaction with Glu331 was not a necessary condition for the compound to act as an inverse agonist. These rules provide guidance for the design of new inverse agonists.     

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.     

Dynophore-Based Approach in Virtual Screening: A Case of Human DNA Topoisomerase IIα

In this study, we utilized human DNA topoisomerase IIα as a model target to outline a dynophore-based approach to catalytic inhibitor design. Based on MD simulations of a known catalytic inhibitor and the native ATP ligand analog, AMP-PNP, we derived a joint dynophore model that supplements the static structure-based-pharmacophore information with a dynamic component. Subsequently, derived pharmacophore models were employed in a virtual screening campaign of a library of natural compounds. Experimental evaluation identified flavonoid compounds with promising topoisomerase IIα catalytic inhibition and binding studies confirmed interaction with the ATPase domain. We constructed a binding model through docking and extensively investigated it with molecular dynamics MD simulations, essential dynamics, and MM-GBSA free energy calculations, thus reconnecting the new results to the initial dynophore-based screening model. We not only demonstrate a new design strategy that incorporates a dynamic component of molecular recognition, but also highlight new derivates in the established flavonoid class of topoisomerase II inhibitors.     

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.     

Molecular Docking and Molecular Dynamics Simulation Studies of Triterpenes from Vernonia patula with the Cannabinoid Type 1 Receptor

A molecular docking approach was employed to evaluate the binding affinity of six triterpenes, namely epifriedelanol, friedelin, α-amyrin, α-amyrin acetate, β-amyrin acetate, and bauerenyl acetate, towards the cannabinoid type 1 receptor (CB1). Molecular docking studies showed that friedelin, α-amyrin, and epifriedelanol had the strongest binding affinity towards CB1. Molecular dynamics simulation studies revealed that friedelin and α-amyrin engaged in stable non-bonding interactions by binding to a pocket close to the active site on the surface of the CB1 target protein. The studied triterpenes showed a good capacity to penetrate the blood–brain barrier. These results help to provide some evidence to justify, at least in part, the previously reported antinociceptive and sedative properties of Vernonia patula.     

Docking and Molecular Dynamic of Microalgae Compounds as Potential Inhibitors of Beta-Lactamase

Bacterial resistance is responsible for a wide variety of health problems, both in children and adults. The persistence of symptoms and infections are mainly treated with β-lactam antibiotics. The increasing resistance to those antibiotics by bacterial pathogens generated the emergence of extended-spectrum β-lactamases (ESBLs), an actual public health problem. This is due to rapid mutations of bacteria when exposed to antibiotics. In this case, β-lactamases are enzymes used by bacteria to hydrolyze the beta-lactam rings present in the antibiotics. Therefore, it was necessary to explore novel molecules as potential β-lactamases inhibitors to find antibacterial compounds against infection caused by ESBLs. A computational methodology based on molecular docking and molecular dynamic simulations was used to find new microalgae metabolites inhibitors of β-lactamase. Six 3D β-lactamase proteins were selected, and the molecular docking revealed that the metabolites belonging to the same structural families, such as phenylacridine (4-Ph), quercetin (Qn), and cryptophycin (Cryp), exhibit a better binding score and binding energy than commercial clinical medicine β-lactamase inhibitors, such as clavulanic acid, sulbactam, and tazobactam. These results indicate that 4-Ph, Qn, and Cryp molecules, homologous from microalgae metabolites, could be used, likely as novel β-lactamase inhibitors or as structural templates for new in-silico pharmaceutical designs, with the possibility of combatting β-lactam resistance     

Molecular Dynamics Simulations of Acylpeptide Hydrolase Bound to Chlorpyrifosmethyl Oxon and Dichlorvos

Acylpeptide hydrolases (APHs) catalyze the removal of N-acylated amino acids from blocked peptides. Like other prolyloligopeptidase (POP) family members, APHs are believed to be important targets for drug design. To date, the binding pose of organophosphorus (OP) compounds of APH, as well as the different OP compounds binding and inducing conformational changes in two domains, namely, α/β hydrolase and β-propeller, remain poorly understood. We report a computational study of APH bound to chlorpyrifosmethyl oxon and dichlorvos. In our docking study, Val471 and Gly368 are important residues for chlorpyrifosmethyl oxon and dichlorvos binding. Molecular dynamics simulations were also performed to explore the conformational changes between the chlorpyrifosmethyl oxon and dichlorvos bound to APH, which indicated that the structural feature of chlorpyrifosmethyl oxon binding in APH permitted partial opening of the β-propeller fold and allowed the chlorpyrifosmethyl oxon to easily enter the catalytic site. These results may facilitate the design of APH-targeting drugs with improved efficacy.     

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.     

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.     

Thermodynamic Studies of Interactions between Sertraline Hydrochloride and Randomly Methylated β-Cyclodextrin Molecules Supported by Circular Dichroism Spectroscopy and Molecular Docking Results

The interaction between sertraline hydrochloride (SRT) and randomly methylated β-cyclodextrin (RMβCD) molecules have been investigated at 298.15 K under atmospheric pressure. The method used—Isothermal Titration Calorimetry (ITC) enabled to determine values of the thermodynamic functions like the enthalpy (ΔH), the entropy (ΔS) and the Gibbs free energy (ΔG) of binding for the examined system. Moreover, the stoichiometry coefficient of binding (n) and binding/association constant (K) value have been calculated from the experimental results. The obtained outcome was compared with the data from the literature for other non-ionic βCD derivatives interacting with SRT and the enthalpy-entropy compensation were observed and interpreted. Furthermore, the connection of RMβCD with SRT was characterized by circular dichroism spectroscopy (CD) and complexes of βCD derivatives with SRT were characterized through the computational studies with the use of molecular docking (MD).     

Combinatorial Virtual Library Screening Study of Transforming Growth Factor-β2–Chondroitin Sulfate System

Transforming growth factor-beta (TGF-β), a member of the TGF-β cytokine superfamily, is known to bind to sulfated glycosaminoglycans (GAGs), but the nature of this interaction remains unclear. In a recent study, we found that preterm human milk TGF-β2 is sequestered by chondroitin sulfate (CS) in its proteoglycan form. To understand the molecular basis of the TGF-β2–CS interaction, we utilized the computational combinatorial virtual library screening (CVLS) approach in tandem with molecular dynamics (MD) simulations. All possible CS oligosaccharides were generated in a combinatorial manner to give 24 di- (CS02), 192 tetra- (CS04), and 1536 hexa- (CS06) saccharides. This library of 1752 CS oligosaccharides was first screened against TGF-β2 using the dual filter CVLS algorithm in which the GOLDScore and root-mean-square-difference (RMSD) between the best bound poses were used as surrogate markers for in silico affinity and in silico specificity. CVLS predicted that both the chain length and level of sulfation are critical for the high affinity and high specificity recognition of TGF-β2. Interestingly, CVLS led to identification of two distinct sites of GAG binding on TGF-β2. CVLS also deduced the preferred composition of the high specificity hexasaccharides, which were further assessed in all-atom explicit solvent MD simulations. The MD results confirmed that both sites of binding form stable GAG–protein complexes. More specifically, the highly selective CS chains were found to engage the TGF-β2 monomer with high affinity. Overall, this work present key principles of recognition with regard to the TGF-β2–CS system. In the process, it led to the generation of the in silico library of all possible CS oligosaccharides, which can be used for advanced studies on other protein–CS systems. Finally, the study led to the identification of unique CS sequences that are predicted to selectively recognize TGF-β2 and may out-compete common natural CS biopolymers.     

The Folding of de Novo Designed Protein DS119 via Molecular Dynamics Simulations

As they are not subjected to natural selection process, de novo designed proteins usually fold in a manner different from natural proteins. Recently, a de novo designed mini-protein DS119, with a βαβ motif and 36 amino acids, has folded unusually slowly in experiments, and transient dimers have been detected in the folding process. Here, by means of all-atom replica exchange molecular dynamics (REMD) simulations, several comparably stable intermediate states were observed on the folding free-energy landscape of DS119. Conventional molecular dynamics (CMD) simulations showed that when two unfolded DS119 proteins bound together, most binding sites of dimeric aggregates were located at the N-terminal segment, especially residues 5–10, which were supposed to form β-sheet with its own C-terminal segment. Furthermore, a large percentage of individual proteins in the dimeric aggregates adopted conformations similar to those in the intermediate states observed in REMD simulations. These results indicate that, during the folding process, DS119 can easily become trapped in intermediate states. Then, with diffusion, a transient dimer would be formed and stabilized with the binding interface located at N-terminals. This means that it could not quickly fold to the native structure. The complicated folding manner of DS119 implies the important influence of natural selection on protein-folding kinetics, and more improvement should be achieved in rational protein design.     

Insights into the Interaction of LVV-Hemorphin-7 with Angiotensin II Type 1 Receptor

Hemorphins are known for their role in the control of blood pressure. Recently, we revealed the positive modulation of the angiotensin II (AngII) type 1 receptor (AT1R) by LVV-hemorphin-7 (LVV-H7) in human embryonic kidney (HEK293) cells. Here, we examined the molecular binding behavior of LVV-H7 on AT1R and its effect on AngII binding using a nanoluciferase-based bioluminescence resonance energy transfer (NanoBRET) assay in HEK293FT cells, as well as molecular docking and molecular dynamics (MD) studies. Saturation and real-time kinetics supported the positive effect of LVV-H7 on the binding of AngII. While the competitive antagonist olmesartan competed with AngII binding, LVV-H7 slightly, but significantly, decreased AngII’s k_D by 2.6 fold with no effect on its B_max. Molecular docking and MD simulations indicated that the binding of LVV-H7 in the intracellular region of AT1R allosterically potentiates AngII binding. LVV-H7 targets residues on intracellular loops 2 and 3 of AT1R, which are known binding sites of allosteric modulators in other GPCRs. Our data demonstrate the allosteric effect of LVV-H7 on AngII binding, which is consistent with the positive modulation of AT1R activity and signaling previously reported. This further supports the pharmacological targeting of AT1R by hemorphins, with implications in vascular and renal physiology.     

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.