Electrostatic interaction-mediated conformational changes of adipocyte fatty acid binding protein probed by molecular dynamics simulation

Abstract

Adipocyte fatty acid binding protein (A-FABP) is a potential drug target for treatment of diabetes, obesity and atherosclerosis. Molecular dynamics (MD) simulations, principal component (PC) analysis and binding free energy calculations were combined to probe effect of electrostatic interactions of residues R78, R106 and R126 with inhibitors ZGB, ZGC and IBP on structural stability of three inhibitor/A-FABP complexes. The results indicate that mutation R126A produces significant influence on polar interactions of three inhibitors with A-FABP and these interactions are main force for driving the conformational change of A-FABP. Analyses on hydrogen bond interactions show that the decrease in hydrogen bonding interactions of residues R126 and Y128 with three inhibitors and the increase in that of K58 with inhibitors ZGC and IBP in the R126A mutated systems mostly regulate the conformational changes of A-FABP. This work shows that R126A can generate a significant perturbation on structural stability of A-FABP, which implies that R126 is of significance in inhibitor bindings. We expect that this study can provide a theoretical guidance for design of potent inhibitors targeting A-FABP.

Communicated by Ramaswamy H. Sarma     

In silico 3D structure modeling and inhibitor binding studies of human male germ cell-associated kinase

Human male germ cell-associated kinase (hMAK) is an androgen-inducible gene in prostate epithelial cells, and it acts as a coactivator of androgen receptor signaling in prostate cancer. The 3D structure of the hMAK kinase was modeled based on the crystal structure of CDK2 kinase using comparative modeling methods, and the ATP-binding site was characterized. We have collected five inhibitors of hMAK from the literature and docked into the ATP-binding site of the kinase domain. Solvated interaction energies (SIE) of inhibitor binding are calculated from the molecular dynamics simulations trajectories of protein–inhibitor complexes. The contribution from each active site residue in hMAK toward inhibitor binding revealed the nature and extent of interactions between inhibitors and individual residues. The main chain atoms of Met79 invariably form hydrogen bonds with all five inhibitors. The amino acids Leu7, Val15, and Leu129 stabilize the inhibitors via CH–pi interactions. The Asp140 in the active site and Glu77 in hinge region show characteristic hydrogen bonding interactions with inhibitors. From SIE, the residue-wise interactions revealed the nature of non-bonding contacts and modifications required to increase the inhibitor activity. Our work provides 3D model structure of hMAK and molecular basis for the mechanisms of hMAK inhibition at atomic level that aid in designing new potent inhibitors.     

Discovery of promising FtsZ inhibitors by E-pharmacophore, 3D-QSAR,molecular docking study,and molecular dynamics simulation

Filamentous temperature-sensitive protein Z (FtsZ), playing a key role in bacterial cell division, is regarded as a promising target for the design of antimicrobial agent. This study is looking for potential high-efficiency FtsZ inhibitors. Ligand-based pharmacophore and E-pharmacophore, virtual screening and molecular docking were used to detect promising FtsZ inhibitors, and molecular dynamics simulation was used to study the stability of protein-ligand complexes in this paper. Sixty-three inhibitors from published literatures with pIC₅₀ ranging from 2.483 to 5.678 were collected to develop ligand-based pharmacophore model. 4DXD bound with 9PC was selected to develop the E-pharmacophore model. The pharmacophore models validated by test set method and decoy set were employed for virtual screening to exclude inactive compounds against ZINC database. After molecular docking, ADME analysis, IFD docking and MM-GBSA, 8 hits were identified as potent FtsZ inhibitors. A 50?ns molecular dynamics simulation was implemented on the compounds to assess the stability between potent inhibitors and FtsZ. The results indicated that the candidate compounds had a high docking score and were strongly combined with FtsZ by forming hydrogen bonding interactions with key amino acid residues, and van der Waals forces and hydrophobic interactions had significant contribution to the stability of the binding. Molecular dynamics simulation results showed that the protein-ligand compounds performed well in both the stability and flexibility of the simulation process.     

MKP1 reduces neuroinflammation via inhibiting endoplasmic reticulum stress and mitochondrial dysfunction

MAP kinase phosphatase 1 (MKP1) has been identified as an antiapoptotic protein via sustaining mitochondrial function. However, the role of MKP1 in neuroinflammation has not been fully understood. The aim of this study is to figure out the influence of MKP1 in lipopolysaccharide (LPS)-treated microglia BV-2 cells and investigate whether MKP1 reduces BV-2 cell death via modulating endoplasmic reticulum (ER) stress and mitochondrial dysfunction. The results of this study demonstrated that MKP1 was rapidly downregulated after exposure to LPS. However, the transfection of MKP1 adenovirus could reverse cell viability and attenuate LPS-mediated BV-2 cell apoptosis. Mechanistically, MKP1 overexpression alleviated ER stress and corrected LPS-induced calcium overloading. Besides, MKP1 adenovirus transfection also reversed mitochondrial bioenergetics, maintained mitochondrial membrane potential, and blocked mitochondria-initiated apoptosis signals. Furthermore, we found that MKP1 overexpression is associated with inactivation of mitogen-activated protein kinase–c-Jun N-terminal kinase (MAPK–JNK) pathway. Interestingly, the activation of MAPK–JNK pathway could abolish the protective effects of MKP1 on BV-2 cells survival and mitochondrial function in the presence of LPS. Altogether, our results identified MKP1 as a primary defender of neuroinflammation via modulating ER stress and mitochondrial function in a manner dependent on MAPK–JNK pathway. These findings may open a new window for the treatment of neuroinflammation in the clinical setting.     

Intramolecular dephosphorylation of ERK by MKP3

The dual specificity mitogen-activated protein kinase phosphatase MKP3 downregulates mitogenic signaling through dephosphorylation of extracellular signal-regulated kinase (ERK). Like other MKPs, MKP3 consists of a noncatalytic N-terminal domain and a catalytic C-terminal domain. ERK binding to the N-terminal noncatalytic domain of MKP3 has been shown to increase (up to 100-fold) the catalytic activity of MKP3 toward small artificial substrates. Here, we address the function of the N-terminal domain of MKP3 in either inter- or intramolecular dephosphorylation of pERK (phosphorylated ERK) and the stoichiometry of the MKP3/pERK Michaelis complex. These are important mechanistic distinctions given the observation that ERK exists in a monomer/dimer equilibrium that is shifted toward the dimer when phosphorylated and given that MKP3 undergoes catalytic activation toward other substrates when bound to ERK. Wild-type and engineered mutants of ERK and MKP3, binding analyses, reaction kinetics, and chemical cross-linking studies were used to demonstrate that the monomer of MKP3 binds to the monomeric form of pERK and that MKP3 within the resulting heterodimer performs intramolecular dephosphorylation of pERK. This study provides the first direct evidence that MKP3 utilizes intramolecular dephosphorylation between a complex consisting of one molecule each of MKP3 and ERK. Catalytic activation and substrate tethering by MKP3 lead to a >or=4000-fold rate enhancement (k(cat)/K(m)) for dephosphorylation of pERK.     

Computational revelation of binding mechanisms of inhibitors to endocellular protein tyrosine phosphatase 1B using molecular dynamics simulations

Endocellular protein tyrosine phosphatase 1B (PTP1B) is one of the most promising target for designing and developing drugs to cure type-II diabetes and obesity. Molecular dynamics (MD) simulations combined with molecular mechanics generalized Born surface area (MM-GBSA) and solvated interaction energy methods were applied to study binding differences of three inhibitors (ID: 901, 941, and 968) to PTP1B, the calculated results show that the inhibitor 901 has the strongest binding ability to PTP1B among the current inhibitors. Principal component (PC) analysis was also carried out to investigate the conformational change of PTP1B, and the results indicate that the associations of inhibitors with PTP1B generate a significant effect on the motion of the WPD-loop. Free energy decomposition method was applied to study the contributions of individual residues to inhibitor bindings, it is found that three inhibitors can generate hydrogen bonding interactions and hydrophobic interactions with different residues of PTP1B, which provide important forces for associations of inhibitors with PTP1B. This research is expected to give a meaningfully theoretical guidance to design and develop of effective drugs curing type-II diabetes and obesity.     

HIV protease (HIV PR) inhibitor structure-activity-selectivity, and active site molecular modeling of high affinity Leu [CH(OH)CH2]Val modified viral and nonviral substrate analogs.

This report details the structure-activity relationships of the HIV gag substrate analog Val-Ser-Gln-Asn-Leu psi[CH(OH)CH2]Val-Ile-Val (U-85548E), an inhibitor exhibiting subnanomolar affinity towards HIV type-1 aspartic proteinase (HIV-1 PR). Our data show that the P1-P2' tripeptidyl sequence provides the minimal chemical determinant for HIV-1 PR binding. We describe the structure-activity properties of Leu psi[CH(OH)CH2]Val substitution in other peptidyl ligands of nonviral substrate origin (e.g., angiotensinogen, insulin and pepstatin). Furthermore, the aspartic proteinase selectivities of a few key compounds are summarized relative to evaluation against human renin, human pepsin, and the fungal enzyme, rhizopuspepsin. These studies have led to the rational design of nanomolar potent inhibitors of both HIV-1 and HIV-2 PR. Finally, a 2.5 A resolution X-ray crystallographic structure of U-85548E complexed to synthetic HIV-1 PR dimer (Jaskolski et al., Biochemistry 30, 1600 [1991]) provided a 3-D picture of the inhibitor bound to the enzyme active site, and we performed computer-assisted molecular modeling studies to explore the possible binding modes of the above series of Leu psi[CH(OH)CH2]Val substituted HIV-1 PR inhibitors.     

A computational insight into binding modes of inhibitors XD29, XD35, and XD28 to bromodomain-containing protein 4 based on molecular dynamics simulations

In the current work, conformational changes of bromodomain-containing protein 4 (1) (BRD4-1) induced by bindings of inhibitors XD29 (57G), XD35 (57F), and XD28 (L28) were investigated using molecular dynamics (MD) simulations and principal component analysis. The results demonstrate that inhibitor bindings produce significant effect on the motion of ZA loop in BRD4-1. Moreover, to further study binding modes of three inhibitors to BRD4-1, binding free energies of inhibitors to BRD4-1 were also calculated using molecular mechanics Poisson-Boltzmann surface area (MM-PBSA) method. The results indicate that van der Waals interactions are main factors to modulate inhibitor bindings. Energy decomposition and hydrogen bond analysis demonstrate that residues Pro82, Leu92, Asn140, and Ile146 play important roles in binding processes of inhibitors to BRD4-1. This study is not only helpful for better understanding function and internal dynamics of BRD4-1, but also can provide a theoretical basis for rational designs of effective anticancer drugs targeting BRD4-1.     

Insights into the structural function of the complex of HIV-1 protease with TMC-126: molecular dynamics simulations and free-energy calculations

The binding properties of the protein–inhibitor complex of human immunodeficiency virus type 1 (HIV-1) protease with the inhibitor TMC-126 are investigated by combining computational alanine scanning (CAS) mutagenesis with binding free-energy decomposition (BFED). The calculated results demonstrate that the flap region (residues 38–58) and the active site region (residues 23–32) in HIV-1 protease contribute 63.72% of the protease to the binding of the inhibitor. In particular, the mechanisms for the interactions of key residues of these species are fully explored and analyzed. Interestingly, the regression analyses show that both CAS and BFED based on the generalized Born model yield similar results, with a correlation coefficient of 0.94. However, compared to CAS, BFED is faster and can decompose the per-residue binding free-energy contributions into backbone and side-chain contributions. The results obtained in this study are useful for studying the binding mechanism between receptor and ligand and for designing potent inhibitors that can combat diseases.     

The dimerization domain of HIV-1 viral infectivity factor Vif is required to block virion incorporation of APOBEC3G

We have obtained the 1.7 Å crystal structure of FIV protease (PR) in which 12 critical residues around the active site have been substituted with the structurally equivalent residues of HIV PR (12X FIV PR). The chimeric PR was crystallized in complex with the broad-based inhibitor TL-3, which inhibits wild type FIV and HIV PRs, as well as 12X FIV PR and several drug-resistant HIV mutants [1–4]. Biochemical analyses have demonstrated that TL-3 inhibits these PRs in the order HIV PR > 12X FIV PR > FIV PR, with K_i values of 1.5 nM, 10 nM, and 41 nM, respectively [2–4]. Comparison of the crystal structures of the TL-3 complexes of 12X FIV and wild-typeFIV PR revealed theformation of additinal van der Waals interactions between the enzyme inhibitor in the mutant PR. The 12X FIV PR retained the hydrogen bonding interactions between residues in the flap regions and active site involving the enzyme and the TL-3 inhibitor in comparison to both FIV PR and HIV PR. However, the flap regions of the 12X FIV PR more closely resemble those of HIV PR, having gained several stabilizing intra-flap interactions not present in wild type FIV PR. These findings offer a structural explanation for the observed inhibitor/substrate binding properties of the chimeric PR.     

Genetic complementation screen identifies a mitogen-activated protein kinase phosphatase, MKP3, as a regulator of dopamine transporter trafficking

The antidepressant and cocaine sensitive plasma membrane monoamine transporters are the primary mechanism for clearance of their respective neurotransmitters and serve a pivotal role in limiting monoamine neurotransmission. To identify molecules in pathways that regulate dopamine transporter (DAT) internalization, we used a genetic complementation screen in Xenopus oocytes to identify a mitogen-activated protein (MAP) kinase phosphatase, MKP3/Pyst1/DUSP6, as a molecule that inhibits protein kinase C-induced (PKC) internalization of transporters, resulting in enhanced DAT activity. The involvement of MKP3 in DAT internalization was verified using both overexpression and shRNA knockdown strategies in mammalian cell models including a dopaminergic cell line. Although the isolation of MKP3 implies a role for MAP kinases in DAT internalization, MAP kinase inhibitors have no effect on internalization. Moreover, PKC-dependent down-regulation of DAT does not correlate with the phosphorylation state of several well-studied MAP kinases (ERK1/2, p38, and SAPK/JNK). We also show that MKP3 does not regulate PKC-induced ubiquitylation of DAT but acts at a more downstream step to stabilize DAT at the cell surface by blocking dynamin-dependent internalization and delaying the targeting of DAT for degradation. These results indicate that MKP3 can act to enhance DAT function and identifies MKP3 as a phosphatase involved in regulating dynamin-dependent endocytosis.     

Insights into structure and function of SHIP2-SH2: homology modeling, docking, and molecular dynamics study

SRC homology 2 (SH2)-containing inositol 5′-phosphatase protein (SHIP2) is a potential target for type 2 diabetes. Its ability to dephosphorylate the lipid messenger phosphatidylinositol 3,4,5-trisphosphate [PtdIns(3,4,5)P3], important for insulin signaling, makes it an important target against type 2 diabetes. The insulin-induced SHIP2 interaction with Shc is very important for the membrane localization and functioning of SHIP2. There is a bidentate relationship between the two proteins where two domains each from SHIP2 and Shc are involved in mutual binding. However in the present study, the SHIP2-SH2 domain binding with the phosphorylated tyrosine 317 on the collagen-homology (CH) domain of Shc, has been studied due to the indispensability of this interaction in SHIP2 localization. In the absence of the crystal structure of SHIP2-SH2, its structural model was developed followed by tracking its molecular interactions with Shc through molecular docking and dynamics studies. This study revealed much about the structural interactions between the SHIP2-SH2 and Shc-CH. Finally, docking study of a nonpeptide inhibitor into the SHIP2-SH2 domain further confirmed the structural interactions involved in ligand binding and also proposed the inhibitor as a major starting point against SHIP2-SH2 inhibition. The insights gained from the current study should prove useful in the design of more potent inhibitors against type 2 diabetes.     

Structural basis of caspase-7 inhibition by XIAP

The inhibitor of apoptosis (IAP) proteins suppress cell death by inhibiting the catalytic activity of caspases. Here we present the crystal structure of caspase-7 in complex with a potent inhibitory fragment from XIAP at 2.45 A resolution. An 18-residue XIAP peptide binds the catalytic groove of caspase-7, making extensive contacts to the residues that are essential for its catalytic activity. Strikingly, despite a reversal of relative orientation, a subset of interactions between caspase-7 and XIAP closely resemble those between caspase-7 and its tetrapeptide inhibitor DEVD-CHO. Our biochemical and structural analyses reveal that the BIR domains are dispensable for the inhibition of caspase-3 and -7. This study provides a structural basis for the design of the next-generation caspase inhibitors.     

TRIM59, amplified in ovarian cancer,promotes tumorigenesis through the MKP3/ERK pathway

Tripartite motif containing 59 (TRIM59) functions as an oncoprotein in various human cancers including ovarian cancer. In this study, we found that TRIM59 gene amplification was prevalent in ovarian cancer tissues, and its amplification was significantly correlated with poorer overall survival. Moreover, knockdown of TRIM59 in SKOV3 and OVCAR3 cells, which had relatively high level of TRIM59, suppressed glucose uptake and lactate production. TRIM59 knockdown also decreased the expression of c-Myc and lactate dehydrogenase A, and the phosphorylation of extracellular signal-regulated kinase (ERK). TRIM59 overexpression in A2780 cells, which expressed low level of TRIM59, showed reverse effects. Notably, treatment with an ERK inhibitor (PD98059) completely abolished the oncogenic effects of TRIM59 overexpression. Interestingly, TRIM59 increased the ubiquitination of MAP kinase phosphatase 3 (MKP3), which may dephosphorylate and inactivate ERK. Ectopic expression of MKP3 inhibited the promoting effects of TRIM59 on glycolysis and the phosphorylation of ERK. TRIM59 protein expression was negatively correlated with MKP3 protein expression in ovarian cancer tissues. Finally, TRIM59 amplification potently affected the anticancer effect of 3-bromopyruvate, an inhibitor of glycolysis, in ovarian cancer cells and patient-derived xenograft. In conclusion, these results suggest that TRIM59 may regulate glycolysis in ovarian cancer via the MKP3/ERK pathway.     

Effect of the active site D25N mutation on the structure, stability, and ligand binding of the mature HIV-1 protease

All aspartic proteases, including retroviral proteases, share the triplet DTG critical for the active site geometry and catalytic function. These residues interact closely in the active, dimeric structure of HIV-1 protease (PR). We have systematically assessed the effect of the D25N mutation on the structure and stability of the mature PR monomer and dimer. The D25N mutation (PR(D25N)) increases the equilibrium dimer dissociation constant by a factor >100-fold (1.3 +/- 0.09 microm) relative to PR. In the absence of inhibitor, NMR studies reveal clear structural differences between PR and PR(D25N) in the relatively mobile P1 loop (residues 79-83) and flap regions, and differential scanning calorimetric analyses show that the mutation lowers the stabilities of both the monomer and dimer folds by 5 and 7.3 degrees C, respectively. Only minimal differences are observed in high resolution crystal structures of PR(D25N) complexed to darunavir (DRV), a potent clinical inhibitor, or a non-hydrolyzable substrate analogue, Ac-Thr-Ile-Nle-r-Nle-Gln-Arg-NH(2) (RPB), as compared with PR.DRV and PR.RPB complexes. Although complexation with RPB stabilizes both dimers, the effect on their T(m) is smaller for PR(D25N) (6.2 degrees C) than for PR (8.7 degrees C). The T(m) of PR(D25N).DRV increases by only 3 degrees C relative to free PR(D25N), as compared with a 22 degrees C increase for PR.DRV, and the mutation increases the ligand dissociation constant of PR(D25N).DRV by a factor of approximately 10(6) relative to PR.DRV. These results suggest that interactions mediated by the catalytic Asp residues make a major contribution to the tight binding of DRV to PR.     

Molecular dynamics (MD) simulations and binding free energy calculation studies between inhibitors and type II dehydroquinase (DHQ2)

Type II dehydroquinase (DHQ2) is the third enzyme of the shikimic acid pathway, and it has been the effective target for tuberculosis (TB). So far, developing multiple potent inhibitors of the DHQ2 of Mycobacterium tuberculosis (DHQ2-Mt) has been considered to be the new therapy to TB. Molecular dynamics simulations followed by molecular mechanics-generalised Born surface area were carried out to calculate the free binding energy and to determine the affinity ability of the four chosen inhibitor molecules, L1, L2, L3 and L4. Energy decomposition analyses show the electrostatic interaction and van der Waals interaction of the ligands to every residue of the DHQ2-Mt. The results suggest that some important residues have different interactions with the four ligands, such as Arg19 and Tyr24. These interactions may have an effect on the ligand binding affinity. The binding affinity of monosubstituted inhibitor is higher than that of disubstituted inhibitor, due to some important interactions with the DHQ2-Mt residues. These computational works will be significant to the theoretical research in the future.     

Crystal Structures of Inhibitor Complexes of Human T-Cell Leukemia Virus (HTLV-1) Protease

Human T-cell leukemia virus type 1 (HTLV-1) is a retrovirus associated with several serious diseases, such as adult T-cell leukemia and tropical spastic paraparesis/myelopathy. For a number of years, the protease (PR) encoded by HTLV-1 has been a target for designing antiviral drugs, but that effort was hampered by limited available structural information. We report a high-resolution crystal structure of HTLV-1 PR complexed with a statine-containing inhibitor, a significant improvement over the previously available moderate-resolution structure. We also report crystal structures of the complexes of HTLV-1 PR with five different inhibitors that are more compact and more potent. A detailed study of structure-activity relationships was performed to interpret in detail the influence of the polar and hydrophobic interactions between the inhibitors and the protease.     

Investigating specificity of the anti-hypertensive inhibitor WNK463 against With-No-Lysine kinase family isoforms via multiscale simulations

Abstract

The With-No-Lysine (WNK) kinase family plays a significant role in regulating cation-chloride cotransporters, blood pressure and body fluid homeostasis. Mutations in the gene of WNK family, especially in WNK1 and WNK4 are responsible for pseudohypoaldosteronism type II (PHAII), characterized by hypertension. The selective inhibition of WNK1 over other isoforms has created an immense challenge in the design of an ATP competitive inhibitor due to their high conservatism. In this work, we have compared the selectivity of the inhibitor WNK463, which was designed for WNK1 with other WNK family isoforms by comprehensive molecular modeling, docking and molecular dynamics simulations in conjunction with the Molecular Mechanics Poisson-Boltzmann Surface Area method. Our calculations show that the affinity of the inhibitor decreases in the order WNK2?>?WNK1?>?WNK3?>?WNK4, in agreement with the experiment. Our study reveals that the inhibitor is most selective to WNK2 due to decreased polar solvation and configurational entropy compared to other isoforms. Furthermore, our analyses indicated that the nonpolar contribution from the hydrophobic residues and hydrogen bonds in the hinge region gatekeeper residue Met³⁰⁴ of WNK1 and its equivalent residue from other kinases played a critical role in stabilizing the inhibitor against WNK kinases. Residues Lys²³³, Met³⁰⁴, Phe³⁵⁶ and Leu³⁶⁹ of WNK1 were the essential residue differences compared to other isoforms that led to specific interactions thereby forming the basis of molecular binding pattern of binding interactions. Overall, we have identified conserved WNK-inhibitor interactions and elucidated isoform-specific interactions that could be exploited in the design of more potent and selective WNK inhibitors.

Communicated by Ramaswamy H. Sarma