Evaluation of Urtica dioica Phytochemicals against Therapeutic Targets of Allergic Rhinitis Using Computational Studies.

Allergic rhinitis (AR) is a prevalent inflammatory condition affecting millions globally, with current treatments often associated with significant side effects. To seek safer and more effective alternatives, natural sources like Urtica dioica (UD) are being explored. However, UD's mechanism of action remains unknown. Therefore, to elucidate it, we conducted an in silico evaluation of UD phytochemicals' effects on known therapeutic targets of allergic rhinitis: histamine receptor 1 (HR1), neurokinin 1 receptor (NK1R), cysteinyl leukotriene receptor 1 (CLR1), chemoattractant receptor-homologous molecule expressed on type 2 helper T cells (CRTH2), and bradykinin receptor type 2 (BK2R). The docking analysis identified amentoflavone, alpha-tocotrienol, neoxanthin, and isorhamnetin 3-O-rutinoside as possessing a high affinity for all the receptors. Subsequently, molecular dynamics (MD) simulations were used to analyze the key interactions; the free energy of binding was calculated through Generalized Born and Surface Area Solvation (MMGBSA), and the conformational changes were evaluated. Alpha-tocotrienol exhibited a high affinity while also inducing positive conformational changes across all targets. Amentoflavone primarily affected CRTH2, neoxanthin targeted NK1R, CRTH2, and BK2R, and isorhamnetin-3-O-rutinoside acted on NK1R. These findings suggest UD's potential to treat AR symptoms by inhibiting these targets. Notably, alpha-tocotrienol emerges as a promising multi-target inhibitor. Further in vivo and in vitro studies are needed for validation.

Interactions of a boron-containing levodopa derivative on D2 dopamine receptor and its effects in a Parkinson disease model.

Levodopa is a cornerstone in Parkinson's disease treatment. Beneficial effects are mainly by binding on D2 receptors. Docking simulations of a set of compounds including well-known D2-ligands and a pool of Boron-Containing Compounds (BCC), particularly boroxazolidones with a tri/tetra-coordinated boron atom, were performed on the D2 Dopamine receptor (D2DR). Theoretical results yielded higher affinity of the compound DPBX, a Dopaboroxazolidone, than levodopa on D2DR. Essential interactions with residues in the third and sixth transmembrane domains of the D2DR appear to be crucial to induce and stabilize interactions in the active receptor state. Results from a motor performance evaluation of a murine model of Parkinson's disease agree with theoretical results, as DPBX showed similar efficacy to that of levodopa for diminishing MPTP-induced parkinsonism. This beneficial effect was disrupted with prior Risperidone (D2DR antagonist) administration, supporting the role of D2DR in the biological effect of DPBX. In addition, DPBX limited neuronal loss in substantia nigra in a similar manner to that of levodopa administration.

Molecular Basis of Inhibitory Mechanism of Naltrexone and Its Metabolites through Structural and Energetic Analyses.

Naltrexone is a potent opioid antagonist with good blood-brain barrier permeability, targeting different endogenous opioid receptors, particularly the mu-opioid receptor (MOR). Therefore, it represents a promising candidate for drug development against drug addiction. However, the details of the molecular interactions of naltrexone and its derivatives with MOR are not fully understood, hindering ligand-based drug discovery. In the present study, taking advantage of the high-resolution X-ray crystal structure of the murine MOR (mMOR), we constructed a homology model of the human MOR (hMOR). A solvated phospholipid bilayer was built around the hMOR and submitted to microsecond (µs) molecular dynamics (MD) simulations to obtain an optimized hMOR model. Naltrexone and its derivatives were docked into the optimized hMOR model and submitted to µs MD simulations in an aqueous membrane system. The MD simulation results were submitted to the molecular mechanics-generalized Born surface area (MMGBSA) binding free energy calculations and principal component analysis. Our results revealed that naltrexone and its derivatives showed differences in protein-ligand interactions; however, they shared contacts with residues at TM2, TM3, H6, and TM7. The binding free energy and principal component analysis revealed the structural and energetic effects responsible for the higher potency of naltrexone compared to its derivatives.

Microsecond MD Simulations to Explore the Structural and Energetic Differences between the Human RXRα-PPARγ vs. RXRα-PPARγ-DNA.

The heterodimeric complex between retinoic X receptor alpha (RXRα) and peroxisome proliferator-activated receptor gamma (PPARγ) is one of the most important and predominant regulatory systems, controlling lipid metabolism by binding to specific DNA promoter regions. X-ray and molecular dynamics (MD) simulations have revealed the average conformation adopted by the RXRα-PPARγ heterodimer bound to DNA, providing information about how multiple domains communicate to regulate receptor properties. However, knowledge of the energetic basis of the protein-ligand and protein-protein interactions is still lacking. Here we explore the structural and energetic mechanism of RXRα-PPARγ heterodimer bound or unbound to DNA and forming complex with co-crystallized ligands (rosiglitazone and 9-cis-retinoic acid) through microsecond MD simulations, molecular mechanics generalized Born surface area binding free energy calculations, principal component analysis, the free energy landscape, and correlated motion analysis. Our results suggest that DNA binding alters correlated motions and conformational mobility within RXRα-PPARγ system that impact the dimerization and the binding affinity on both receptors. Intradomain correlated motions denotes a stronger correlation map for RXRα-PPARγ-DNA than RXRα-PPARγ, involving residues at the ligand binding site. In addition, our results also corroborated the greater role of PPARγ in regulation of the free and bound DNA state.

131I-C19 Iodide Radioisotope and Synthetic I-C19 Compounds as K-Ras4B-PDE6δ Inhibitors: A Novel Approach against Colorectal Cancer-Biological Characterization, Biokinetics and Dosimetry.

In 40-50% of colorectal cancer (CRC) cases, K-Ras gene mutations occur, which induce the expression of the K-Ras4B oncogenic isoform. K-Ras4B is transported by phosphodiesterase-6δ (PDE6δ) to the plasma membrane, where the K-Ras4B-PDE6δ complex dissociates and K-Ras4B, coupled to the plasma membrane, activates signaling pathways that favor cancer aggressiveness. Thus, the inhibition of the K-Ras4B-PDE6δ dissociation using specific small molecules could be a new strategy for the treatment of patients with CRC. This research aimed to perform a preclinical proof-of-concept and a therapeutic potential evaluation of the synthetic I-C19 and 131I-C19 compounds as inhibitors of the K-Ras4B-PDE6δ dissociation. Molecular docking and molecular dynamics simulations were performed to estimate the binding affinity and the anchorage sites of I-C19 in K-Ras4B-PDE6δ. K-Ras4B signaling pathways were assessed in HCT116, LoVo and SW620 colorectal cancer cells after I-C19 treatment. Two murine colorectal cancer models were used to evaluate the I-C19 therapeutic effect. The in vivo biokinetic profiles of I-C19 and 131I-C19 and the tumor radiation dose were also estimated. The K-Ras4B-PDE6δ stabilizer, 131I-C19, was highly selective and demonstrated a cytotoxic effect ten times greater than unlabeled I-C19. I-C19 prevented K-Ras4B activation and decreased its dependent signaling pathways. The in vivo administration of I-C19 (30 mg/kg) greatly reduced tumor growth in colorectal cancer. The biokinetic profile showed renal and hepatobiliary elimination, and the highest radiation absorbed dose was delivered to the tumor (52 Gy/74 MBq). The data support the idea that 131I-C19 is a novel K-Ras4B/PDE6δ stabilizer with two functionalities: as a K-Ras4B signaling inhibitor and as a compound with radiotherapeutic activity against colorectal tumors.

Exploring the inhibitory activity of valproic acid against the HDAC family using an MMGBSA approach.

Valproic acid (VPA) is a compound currently used in clinical practice for the treatment of epilepsy as well as bipolar and mood disorders. VPA targets histone deacetylases (HDACs), which participate in the removal of acetyl groups from lysine in several proteins, regulating a wide variety of functions within the organism. An imbalance or malfunction of these enzymes is associated with the development and progression of several diseases, such as cancer and neurodegenerative diseases. HDACs are divided into four classes, but VPA only targets Class I (HDAC1-3 and 8) and Class IIa (HDAC4-5, 7 and 9) HDACs; however, structural and energetic information regarding the manner by which VPA inhibits these HDACs is lacking. Here, the structural and energetic features that determine this recognition were studied using molecular docking and molecular dynamics (MD) simulation. It was found that VPA reaches the catalytic site in HDAC1-3 and 7, whereas in HDAC6, VPA only reaches the catalytic tunnel. In HDAC4, VPA was bound adjacent to L1 and L2, a zone that participates in corepressor binding, and in HDAC8, VPA was bound to the hydrophobic active site channel (HASC), in line with previous reports.

Dissecting the molecular recognition of dual lapatinib derivatives for EGFR/HER2.

Abnormalities in the expression levels of EGFR/HER2 are found in many different types of human cancer; therefore, the design of dual inhibitors of EGFR/HER2 is a recognized anti-cancer strategy. Some lapatinib derivatives have been previously synthesized by modification at the methylsulfonylethylaminomethylfuryl group and biologically evaluated, demonstrating that the 2i compound shows potent inhibitory activity against EGFR/HER2-overexpressing cancer cells. In the present study, we explored the structural and energetic features that guide the molecular recognition of 2i using various EGFR/HER2 states. Molecular dynamics (MD) simulation with an MMPB(GB)SA approach was used to generate the inactive EGFR/HER2-ligand complexes. Our results corroborate that slight modification of lapatinib contributes to an increase in the affinity of the 2i compound for inactive EGFR/HER2 as compared with lapatinib compound, which is in accordance with experimental results. Comparison with previous results reveals that lapatinib and its derivative bind more strongly to the inactive than the intermediate active-inactive HER2 state. Principal component analysis allowed the observation that coupling of 2i to EGFR/HER2 is linked to a reduction in the conformational mobility, which may also contribute to the improvement in affinity observed for this compound as compared with lapatinib.

Impact of tetramerization on the ligand recognition of N1 influenza neuraminidase via MMGBSA approach.

Influenza virus neuraminidase (NA) is a homotetrameric surface protein that, in contrast to other non-influenza NAs, requires a quaternary assembly to exhibit enzymatic activity, suggesting that the oligomeric state significantly impacts the active site of influenza NA. Nevertheless, most structure-based drug design studies have been reported by employing the monomeric state in the closed or open-loop due to the computational cost of employing the tetrameric NA. In this work, we present MD simulations coupled to the MMGBSA approach of avian N1 type NA in its monomeric and tetrameric closed and open-loop state both with and without the inhibitor oseltamivir and its natural substrate, sialic acid. Structural and energetic analyses revealed that the tetrameric state impacts flexibility as well as the map of interactions participating in stabilizing the protein-ligand complexes with respect to the monomeric state. It was observed that the tetrameric state exerts dissimilar effects in binding affinity, characteristic of positive and negative cooperativity for oseltamivir and sialic acid, respectively. Based on our results, to perform a confident structure-based drug design, as well as to evaluate the impact of key mutations through MD simulations, it is important to consider the tetrameric state closed-loop state.

Molecular mechanism of the association and dissociation of Deltarasin from the heterodimeric KRas4B-PDEδ complex.

The formation of the KRas4B-PDEδ complex activates different signaling pathways required for the development and maintenance of cancer. Previous experimental and theoretical studies have allowed researchers to design an inhibitor of the KRas4B-PDEδ complex, "Deltarasin." This inhibitor binds to the prenyl-binding pocket of PDEδ and subsequently inhibits the proliferation of human pancreatic ductal adenocarcinoma cells that depend on oncogenic KRas4B. Nevertheless, structural and energetic information about the inhibitory effects of Deltarasin on the KRas4B-PDEδ complex are not available. In this study, we explore the properties of Deltarasin in inhibiting the formation of wild-type and mutant KRas4B-PDEδ complexes present in different cell lines expressing mutant RAS genes (G12D, G12C, G12V, G13D, Q61L, and Q61R) using 1.7 µs molecular dynamics (MD) simulations in combination with the MMGBSA approach. Our results revealed the energetic and structural mechanisms that suggest a higher affinity of Deltarasin for PDEδ than the farnesylated HVR. Moreover, Deltarasin exerts another dissociative effect by binding to the protein-protein dimeric interface of wild-type KRas4B-PDEδ, whereas associative and dissociative effects were observed for mutant KRas4B-PDEδ, providing a mechanistic explanation for the inhibitory effects of Deltarasin on different cancer cell lines.

Cell-based assays and molecular dynamics analysis of a boron-containing agonist with different profiles of binding to human and guinea pig beta2 adrenoceptors.

The design of beta2 adrenoceptor (ß2AR) agonists is attractive because of their wide-ranging applications in medicine, and the details of agonist interactions with ß2AR are interesting because it is considered a prototype for G-protein coupled receptors. Preclinical studies for agonist development have involved biological assays with guinea pigs due to a similar physiology to humans. Boron-containing Albuterol derivatives (BCADs) designed as bronchodilators have improved potency and efficacy compared with their boron-free precursor on guinea pig ß2ARs (gpß2ARs), and two of the BCADs (BR-AEA and boronterol) conserve these features on cells expressing human ß2ARs (hß2ARs). The aim of this study was to test the BCAD Politerol on gpß2ARs and hß2ARs in vitro and in silico. Politerol displayed higher potency and efficacy on gpß2AR than on hß2AR in experimental assays, possible explanations are provided based on molecular modeling, and molecular dynamics simulations of about 0.25 µs were performed for the free and bound states adding up to 2 µs in total. There were slight differences, particularly in the role of the boron atom, in the interactions of Politerol with gpß2ARs and hß2ARs, affecting movements of transmembrane domains 5-7, known to be pivotal in receptor activation. These findings could be instrumental in the design of compounds selective for hß2ARs.

KRas4B-PDE6δ complex stabilization by small molecules obtained by virtual screening affects Ras signaling in pancreatic cancer.

BACKGROUND: The GTPase KRas4B has been utilized as a principal target in the development of anticancer drugs. PDE6δ transports KRas4B to the plasma membrane, where it is released to activate various signaling pathways required for the initiation and maintenance of cancer. Therefore, identifying new small molecules that prevent activation of this GTPase by stabilizing the KRas4B-PDE6δ molecular complex is a practical strategy to fight against cancer. METHODS: The crystal structure of the KRas4B-PDE6δ heterodimer was employed to locate possible specific binding sites at the protein-protein interface region. Virtual screening of Enamine-database compounds was performed on the located potential binding sites to identify ligands able to simultaneously bind to the KRas4B-PDE6δ heterodimer. A molecular dynamics approach was used to estimate the binding free-energy of the complex. Cell viability and apoptosis were measured by flow cytometry. G-LISA was used to measure Ras inactivation. Western blot was used to measure AKT and ERK activation. MIA PaCa-2 cells implanted subcutaneously into nude mice were treated with D14 or C22 and tumor volumes were recorded. RESULTS: According to the binding affinity estimation, D14 and C22 stabilized the protein-protein interaction in the KRas4B-PDE6δ complex based on in vitro evaluation of the 38 compounds showing antineoplastic activity against pancreatic MIA PaCa-2 cancer cells. In this work, we further investigated the antineoplastic cellular properties of two of them, termed D14 and C22, which reduced the viability in the human pancreatic cancer cells lines MIA PaCa-2, PanC-1 and BxPC-3, but not in the normal pancreatic cell line hTERT-HPNE. Compounds D14 and C22 induced cellular death via apoptosis. D14 and C22 significantly decreased Ras-GTP activity by 33% in MIA PaCa-2 cells. Moreover, D14 decreased AKT phosphorylation by 70% and ERK phosphorylation by 51%, while compound C22 reduced AKT phosphorylation by 60% and ERK phosphorylation by 36%. In addition, compounds C22 and D14 significantly reduced tumor growth by 88.6 and 65.9%, respectively, in a mouse xenograft model. CONCLUSIONS: We identified two promising compounds, D14 and C22, that might be useful as therapeutic drugs for pancreatic ductal adenocarcinoma treatment.

The small organic molecule C19 binds and strengthens the KRAS4b-PDEδ complex and inhibits growth of colorectal cancer cells in vitro and in vivo.

BACKGROUND: Colorectal cancer is the third most common cancer worldwide; and in 40% of all cases, KRAS4b-activating mutations occur. KRAS4b is transported by phosphodiesterase-6δ (PDEδ) to the plasma membrane, where it gets activated. PDEδ downregulation prevents redistribution and activation of KRAS4b. Thus, targeting the KRAS4b-PDEδ complex is a treatment strategy for colorectal cancer. METHODS: Using docking and molecular dynamics simulations coupled to molecular mechanics, the generalized born model and solvent accessibility (MMGBSA) approach to explore protein-ligand stability, we found that the compound ((2S)-N-(2,5-diclorofenil)-2-[(3,4-dimetoxifenil)metilamino]-propanamida), termed C19, bound and stabilized the KRAS4b-PDEδ complex. We investigated whether C19 decreases the viability and proliferation of colorectal cancer cells, in addition to knowing the type of cell death that it causes and if C19 decreases the activation of KRAS4b and their effectors. RESULTS: C19 showed high cytotoxicity in the colorectal cancer cell lines HCT116 and LoVo, with a stronger effect in KRAS-dependent LoVo cells. Importantly, C19 significantly decreased tumor size in a xenograft mouse model and showed lower side effects than 5-fluorouracil that is currently used as colorectal cancer treatment. CONCLUSIONS: Mechanistically, the cytotoxic effect was due to increased apoptosis of tumor cells and decreased phosphorylation of Erk and Akt. Therefore, our results suggest that C19 may serve as a promising new treatment for colorectal cancer.

Exploring the conformational and binding properties of unphosphorylated/phosphorylated monomeric and trimeric Bcl-2 through docking and molecular dynamics simulations.

B-cell lymphoma (Bcl-2) is commonly associated with the progression and preservation of cancer and certain lymphomas; therefore, it is considered as a biological target against cancer. Nevertheless, evidence of all its structural binding sites has been hidden because of the lack of a complete Bcl-2 model, given the presence of a flexible loop domain (FLD), which is responsible for its complex behavior. FLD region has been implicated in phosphorylation, homotrimerization, and heterodimerization associated with Bcl-2 antiapoptotic function. In this contribution, homology modeling, molecular dynamics (MD) simulations in the microsecond (µs) time-scale and docking calculations were combined to explore the conformational complexity of unphosphorylated/phosphorylated monomeric and trimeric Bcl-2 systems. Conformational ensembles generated through MD simulations allowed for identifying the most populated unphosphorylated/phosphorylated monomeric conformations, which were used as starting models to obtain trimeric complexes through protein-protein docking calculations, also submitted to µs MD simulations. Principal component analysis showed that FLD represents the main contributor to total Bcl-2 mobility, and is affected by phosphorylation and oligomerization. Subsequently, based on the most representative unphosphorylated/phosphorylated monomeric and trimeric Bcl-2 conformations, docking studies were initiated to identify the ligand binding site of several known Bcl-2 inhibitors to explain their influence in homo-complex formation and phosphorylation. Docking studies showed that the different conformational states experienced by FLD, such as phosphorylation and oligomerization, play an essential role in the ability to make homo and hetero-complexes. © 2016 Wiley Periodicals, Inc. Biopolymers 105: 393-413, 2016.

Targeting quorum sensing by designing azoline derivatives to inhibit the N-hexanoyl homoserine lactone-receptor CviR: Synthesis as well as biological and theoretical evaluations.

To counteract bacterial resistance, we investigated the interruption of quorum sensing mediated by non-classical bioisosteres of the N-hexanoyl homoserine lactone with an azoline core. For this purpose, a set of selected 2-substituted azolines was synthesized, establishing the basis for a new protocol to synthesize 2-amino imidazolines. The synthesized compounds were evaluated as inhibitors of violacein production in Chromobacterium violaceum. Theoretical studies on bioisostere-protein interactions were performed using CviR. The results show that some azolines decreased violacein production, suggesting an antiquorum sensing profile against Gram-negative bacteria. Docking and molecular dynamic simulations together with binding free energy calculations revealed the exact binding and inhibitory profiles. These theoretical results show relationship with the in vitro activity of the azoline series.

Binding free energy calculations between bovine ß-lactoglobulin and four fatty acids using the MMGBSA method.

The bovine dairy protein ß-lactoglobulin (ßlg) is a promiscuous protein that has the ability to bind several hydrophobic ligands. In this study, based on known experimental data, the dynamic interaction mechanism between bovine ßlg and four fatty acids was investigated by a protocol combining molecular dynamics (MD) simulations and molecular mechanics generalized Born surface area (MMGBSA) binding free energy calculations. Energetic analyses revealed binding free energy trends that corroborated known experimental findings; larger ligand size corresponded to greater binding affinity. Finally, binding free energy decomposition provided detailed information about the key residues stabilizing the complex.

Ligand entry into the calyx of ß-lactoglobulin.

Although the thermodynamic principles that control the binding of drug molecules to their protein targets are well understood, the detailed process of how a ligand reaches a protein binding site has been an intriguing question over decades. The short time interval between the encounter between a ligand and its receptor to the formation of the stable complex has prevented experimental observations. Bovine ß-lactoglobulin (ßlg) is a lipocalin member that carries fatty acids (FAs) and other lipids in the cellular environment. Βlg accommodates a FA molecule in its highly hydrophobic cavity and exhibits the capability of recognizing a wide variety of hydrophobic ligands. To elucidate the ligand entry process on ßlg, we report molecular dynamics simulations of the encounter between palmitate (PA) or laurate (LA) and ßlg. Our results show that residues localized in loops at the cavity entrance play an important role in the ligand penetration process. Analysis of the short-term interaction energies show that the forces operating on the systems lead to average conformations very close to the crystallographic holo-forms. Whereas the binding free energy analysis using the molecular mechanics Generalized Born surface area method shows that these conformations were thermodynamically favorable.

Seeking potential anticonvulsant agents that target GABAA receptors using experimental and theoretical procedures.

The aim of this study was to identify compounds that possess anticonvulsant activity by using a pentylenetetrazol (PTZ)-induced seizure model. Theoretical studies of a set of ligands, explored the binding affinities of the ligands for the GABA(A) receptor (GABA(A)R), including some benzodiazepines. The ligands satisfy the Lipinski rules and contain a pharmacophore core that has been previously reported to be a GABA(A)R activator. To select the ligands with the best physicochemical properties, all of the compounds were analyzed by quantum mechanics and the energies of the highest occupied molecular orbital and lowest unoccupied molecular orbital were determined. Docking calculations between the ligands and the GABA(A)R were used to identify the complexes with the highest Gibbs binding energies. The identified compound D1 (dibenzo(b,f)(1,4)diazocine-6,11(5H,12H)-dione) was synthesized, experimentally tested, and the GABA(A)R-D1 complex was submitted to 12-ns-long molecular dynamics (MD) simulations to corroborate the binding conformation obtained by docking techniques. MD simulations were also used to analyze the decomposition of the Gibbs binding energy of the residues involved in the stabilization of the complex. To validate our theoretical results, molecular docking and MD simulations were also performed for three reference compounds that are currently in commercial use: clonazepam (CLZ), zolpidem and eszopiclone. The theoretical results show that the GABA(A)R-D1, and GABA(A)R-CLZ complexes bind to the benzodiazepine binding site, share a similar map of binding residues, and have similar Gibbs binding energies and entropic components. Experimental studies using a PTZ-induced seizure model showed that D1 possesses similar activity to CLZ, which corroborates the predicted binding free energy identified by theoretical calculations.

Evaluation of structural and thermodynamic insight of ERß with DPN and derivatives through MMGBSA/MMPBSA methods.

Estrogen receptors (ERs) are nuclear factors that exist as two subtypes: ERα and ERß. Among the different selective ERß agonist ligands, the widely used ERß-selective agonist DPN (diarylpropionitrile) is highlighted. Recent experimental and thermodynamic information between R-DPN and S-DPN enantiomers with ERß is important for evaluating further the ability of MD simulations combined with end-point methods to reproduce experimental binding affinity and generate structural insight not provided through crystallographic data. In this research, starting from crystallographic data and experimental binding affinities, we explored the structural and thermodynamic basis of the molecular recognition of ERß with DPN and derivatives through triplicate MD simulations combined with end-point methods. Conformational analysis showed some regions with the highest mobility linked to ligand association that, at the time, impacted the total protein fluctuation. Binding free energy (ΔG) analysis revealed that the Molecular Mechanics Generalized-Born Surface Area (MMGBSA) approach was able to reproduce the experimental tendency with a strong correlation (R = 0.778), whereas per-residue decomposition analysis revealed that all the systems interacted strongly with eight residues (L298, E305, L339, M340, L343, F356, H475, and L476). The comparison between theoretical studies using the MMGBSA approach with experimental results provides new insights for drug designing of new DPN derivatives.

Unveiling the Mechanisms of EGCG-p53 Interactions through Molecular Dynamics Simulations.

Green tea consumption is associated with protective and preventive effects against various types of cancer. These effects are attributed to polyphenols, particularly epigallocatechin-3-gallate (EGCG). EGCG acts by directly inhibiting tumor suppressor protein p53. The binding mechanism by which EGCG inhibits p53 activity is associated with residues Trp23-Lys24 and Pro47-Thr55 within the p53 N-terminal domain (NTD). However, the structural and thermodynamic aspects of the interaction between EGCG and p53 are poorly understood. Therefore, based on crystallographic data, we combine docking, molecular dynamics (MD) simulations, and molecular mechanics generalized Born surface area approaches to explore the intricacies of the EGCG-p53 binding mechanism. A triplicate microsecond MD simulation for each system is initially performed to capture diverse p53 NTD conformations. From the start, the most populated cluster of the second run (R2-1) stands out due to a unique opening between Trp23 and Trp53. During MD simulations, this conformation allows EGCG to sustain a high level of stability and affinity while interacting with both regions of interest and deepening the binding pocket. Structural analysis emphasizes the significance of pyrogallol motifs in EGCG binding. Therefore, the conformational shift in this gap is pivotal, enabling EGCG to impede p53 interactions and manifest its anticancer properties. These findings enhance the present comprehension of the anticancer properties of green tea polyphenols and pave the way for future therapeutic developments.

Determining Structural Changes for Ligand Recognition between Human and Rat Phosphorylated BACE1 in Silico and Its Phosphorylation by GSK3ß at Thr252 by in Vitro Studies.

Alzheimer's disease (AD) is a neurodegenerative disease affecting older adults. AD pathogenesis involves the production of the highly neurotoxic amyloid-ß peptide 1-42 (Aß1-42) from ß-site amyloid precursor protein cleaving enzyme 1 (BACE1). The phosphorylation of BACE1 at Thr252 increases its enzymatic activity. This study examined the phosphorylation of BACE1 from human and rat BACE1 in silico through phosphorylation predictors. Besides, we explored how phosphorylation at various sites affected the BACE1 structure and its affinity with amyloid precursor protein (APP) and six BACE1 inhibitors. Additionally, we evaluated the phosphorylation of Thr252-BACE1 by glycogen synthase kinase 3 ß (GSK3ß) in vitro. The phosphorylation predictors showed that Thr252, Ser59, Tyr76, Ser71, and Ser83 could be phosphorylated. Also, Ser127 in rat BACE1 can be phosphorylated, but human BACE1 has a Gly at this position. Molecular dynamics simulations showed that Ser127 plays an important role in the open and closed BACE1 conformational structures. Docking studies and the molecular mechanics generalized Born surface area (MMGBSA) approach showed that human BACE1 phosphorylated at Thr252 and rat BACE1 phosphorylated at Ser71 have the best binding and free energy with APP, forming hydrogen bonds with Asp672. Importantly, inhibitors have a higher affinity for the phosphorylated rat BACE1 than for its human counterpart, which could explain their failure during clinical trials. Finally, in vitro experiments showed that GSK3ß could phosphorylate BACE1. In conclusion, BACE1 phosphorylation influences the BACE1 conformation and its recognition of ligands and substrates. Thus, these features should be carefully considered in the design of BACE1 inhibitors.