Deprotonation states of the two active site water molecules regulate the binding of protein phosphatase 5 with its substrate: A molecular dynamics study

Targeting shikimate pathway: In silico analysis of phosphoenolpyruvate derivatives as inhibitors of EPSP synthase and DAHP synthase

The shikimate pathway consists of seven enzymatic steps involved in the conversion of erythrose-4-phosphate and phosphoenolpyruvate to chorismate and also responsible to the production of aromatic amino acids, such as phenylalanine, tyrosine, and tryptophan which are essential to the bacterial metabolism. The 3-deoxy-D-arabino-heptulosonate-7-phosphate synthase (DAHPS) and 5-enolpyruvylshikimate 3-phosphate synthase (EPSPS) catalyze important steps in the shikimate pathway using as substrate the phosphoenolpyruvate (PEP). Due to the importance of PEP in shikimate pathway, its structure has been investigated to develop new bioinspired competitive inhibitors against DAHPS and EPSPS. In the present study, we perform a literature survey of 28 PEP derivatives, then we analyzed the selectivity and affinity of these compounds against the EPSPS and DAHPS structures using consensual molecular docking, pharmacophore prediction, molecular dynamics (MD) simulations, and binding free energy calculations. Here, we propose consistent binding modes of the selected ligands and indicate that their structures show interesting pharmacophoric properties related to multi-targets inhibitors for both enzymes. Our computational results are supported by previous experimental findings related to the interactions of PEP derivatives with DAHPS and EPSPS structures.

Quantum-Based Modeling of Dephosphorylation in the Catalytic Site of Serine/Threonine Protein Phosphatase-5 (PPP5C)

Serine/threonine protein phosphatase-5 (PP5; PPP5C) is a member of the phosphoprotein phosphatase (PPP) gene family. The PPP catalytic domains feature a bimetal system (M1/M2), an associated bridge hydroxide (W1(OH−)), an M1-bound water/hydroxide (W2), and a highly conserved core sequence. The PPPs are presumed to share a common mechanism: The seryl/threonyl phosphoryl group of the phosphoprotein coordinates the metal ions, W1(OH−) attacks the central phosphorous atom, rupturing the antipodal phosphoester bond and releasing the phosphate-free protein. Also, a histidine/aspartate tandem is responsible for protonating the exiting seryl/threonyl alkoxide. Here, we employed quantum-based computations on a large section of the PP5 catalytic site. A 33-residue, ONIOM(UB3LYP/6-31G(d):UPM7) model was built to perform computations using methylphosphate dianion as a stand-in substrate for phosphoserine/phosphothreonine. We present a concerted transition state (TS) in which W1(OH−) attacks the phosphate center at the same time that the exiting seryl/threonyl alkoxide is protonated directly by the His304/Asp274 tandem, with W2 assigned as a water molecule: W2(H2O). Arg275, proximal to M1, stabilizes the substrate and TS by binding both the ester oxygen (Oγ) and a phosphoryl oxygen (O1) in a bidentate fashion; in the product state, Tyr451 aids in decoupling Arg275 from O1 of the product phosphate ion. The reaction is exothermic (ΔH = −2.0 kcal/mol), occurs in a single step, and has a low activation barrier (ΔH‡ = +10.0 kcal/mol). Our work is an improvement over an earlier computational study that also found bond rupture and alkoxide protonation to be concerted, but concluded that Arg275 is deprotonated during the reactant and TS stages of the pathway. In that earlier study, the critical electron-withdrawal role that Arg275 plays during the hydroxide attack was not correctly accounted for.

Discovery of novel serine/threonine protein phosphatase 1 inhibitors from traditional Chinese medicine through virtual screening and biological assays

Protein phosphatase 1 (PP1) is a critical regulator of several processes, such as muscle contraction, neuronal signaling, glycogen synthesis, and cell proliferation. Dysregulation of PP1 has recently been found to be implicated in cardiac dysfunctions, which indicates that PP1 could be an attractive therapeutic target. However, discovery of PP1 inhibitors with satisfied safety and efficiency is still a challenge. Here, in order to discover potential PP1 inhibitors, compounds extracted from traditional Chinese medicine (TCM) were screened by a novel integrated virtual screening protocol including pharmacophore modeling and docking approaches. Combined with protein phosphatase inhibition assay, ZINC43060554 showed strongly inhibitory activity with IC₅₀ values of 26.78 μM. Furthermore, molecular dynamics simulation and Molecular Mechanics/Generalized Born Surface Area binding free-energy analysis were performed to examine the stability of ligand binding modes. These novel scaffolds discovered in the present study can be used for rational design of PP1 inhibitors with high affinity.Communicated by Ramaswamy H. Sarma.

Multiscale simulation of the interaction of calreticulin-thrombospondin-1 complex with a model membrane microdomain

Cell surface calreticulin (CRT) binding to thrombospondin-1 (TSP1), regulates cell adhesion, migration, anoikis resistance, and collagen production. Due to the essential role of membrane microdomains in CRT-mediated focal adhesion disassembly, we previously studied the effect of raft-like bilayers on TSP1-CRT interactions with all-atom molecular dynamics (AAMD) simulations. However, the simulated systems of protein on the surface of the bilayer(s) in the explicit solvent are too large for long timescale AAMD simulations due to computational expense. In this study, we adopted a multiscale modeling approach of combining AAMD, coarse-grained molecule dynamics (CGMD), and reversed AAMD (REV AAMD) simulations to investigate the interactions of single CRT or of the TSP1-CRT complex with a membrane microdomain at microsecond timescale. Results showed that CRT conformational stabilization by binding of TSP1 in AAMD simulation was undetectable in CGMD simulation, but it was recovered in REV AAMD simulation. Similarly, interactions of the CRT N-domain and TSP1 with the membrane microdomain were lost in CGMD simulations but they were re-gained in the REV AAMD simulations. There was the higher coordination of the CRT P-domain in the TSP1-CRT complex with the lipid components of membrane microdomain compared to that of single CRT, which could directly affect the conformation of CRT and further mediate CRT recruitment of LDL receptor-related protein for signaling events. This study provides structural and molecular insights into TSP1-CRT interactions in a membrane microdomain environment and demonstrates the feasibility of using multiscale simulations to investigate the interactions between protein and membrane microdomains at a long timescale.

Molecular insights into the specific recognition between the RNA binding domain qRRM2 of hnRNP F and G-tract RNA: A molecular dynamics study

Heterogeneous nuclear ribonucleoprotein F (hnRNP F) controls the expression of various genes through regulating the alternative splicing of pre-mRNAs in the nucleus. It uses three quasi-RNA recognition motifs (qRRMs) to recognize G-tract RNA which contains at least three consecutive guanines. The structures containing qRRMs of hnRNP F in complex with G-tract RNA have been determined by nuclear magnetic resonance (NMR) spectroscopy, shedding light on the recognition mechanism of qRRMs with G-tract RNA. However, knowledge of the recognition details is still lacking. To investigate how qRRMs specifically bind with G-tract RNA and how the mutations of any guanine to an adenine in the G-tract affect the binding, molecular dynamics simulations with binding free energy analysis were performed based on the NMR structure of qRRM2 in complex with G-tract RNA. Simulation results demonstrate that qRRM2 binds strongly with G-tract RNA, but any mutation of the G-tract leads to a drastic reduction of the binding free energy. Further comparisons of the energetic components reveal that van der Waals and non-polar interactions play essential roles in the binding between qRRM2 and G-tract RNA, but the interactions are weakened by the effect of RNA mutations. Structural and dynamical analyses indicate that when qRRM2 binds with G-tract RNA, both qRRM2 and G-tract maintain stabilized structures and dynamics; however, the stability is disrupted by the mutations of the G-tract. These results provide novel insights into the recognition mechanism of qRRM2 with G-tract RNA that are not elucidated by the NMR technique.