Aromatic-aromatic interactions in the formation of the MDM2-p53 complex

Although the molecular interaction of MDM2 with the transactivation domain of p53 has been thoroughly studied, there is very limited information regarding the steps involved in the recognition mechanism between these proteins. On this basis, we performed four high-temperature molecular dynamics simulations in explicit solvent to gain insight into the interactions involved in the fist contact toward the formation of the complex. We found that the presence of specific intermolecular aromatic pairs at the interface of the complex, around the native-like state of MDM2, is consistent among independent molecular dynamics runs. This observation suggests that aromatic-aromatic interactions are closely related to the first contact between MDM2 and p53. Thus, we propose that aromatic-aromatic interactions are an important, and probably essential, requirement for the formation and stabilization of the MDM2-p53 complex.     

2D-QSAR study,molecular docking,and molecular dynamics simulation studies of interaction mechanism between inhibitors and transforming growth factor-beta receptor I (ALK5)

Transforming growth factor type 1 receptor (ALK5) is kinase associated with a wide variety of pathological processes, and inhibition of ALK5 is a good strategy to treat many kinds of cancer and fibrotic diseases. Recently, a series of compounds have been synthesized as ALK5 inhibitors. However, the study of their selectivity against other potential targets remains elusive. In this research, a data-set of ALK5 inhibitors were collected and studied based on the combination of 2D-QSAR, molecular docking and molecular dynamics simulation. The quality of QSAR models were assessed statistically by F, R², and R²_ADJ, proved to be credible. The cross-validations for the models (q²_LOO = 0.571 and 0.629, respectively) showed their robustness, while the external validations (r²_test = 0.703 and 0.764, respectively) showed their predictive power. Besides, the predicted binding free energy results calculated by MM/GBSA method were in accordance with the experimental data, and the van der Waals energy term was the factor that had the most significant impact on ligand binding. What is more, several important residues were found to significantly affect the binding affinity. Finally, based on our analyses above, a proposed series of molecules were designed.     

Proapoptotic modification of substituted isoindolinones as MDM2-p53 inhibitors

A series of novel amino acid ester derivatives of 2,3-substituted isoindolinones was synthesized and evaluated for p53-mediated apoptotic activity. The rationale for augmentation of the target activity of 2,3-substituted isoindolinones was based on the introduction of new fragments in the structure of the inhibitor that would provide additional binding sites in the hydrophobic cavity of MDM2. To select for the anticipated modifications we employed molecular docking. Synthesized molecules were evaluated for their ability to induce apoptosis in two cancer cell lines and their derivatives with different status of p53 (colorectal HCT116 and osteosarcoma U2OS cells) by Annexin V staining. The target activity was estimated using high-content imaging system Operetta. Valine and phenylglycine ester derivatives were identified as potentially active MDM2-p53 inhibitors.     

Targeting the conformational transitions of MDM2 and MDMX: Insights into key residues affecting p53 recognition

The oncogenic proteins MDM2 and MDMX have distinct and critical roles in the control of the activity of the p53 tumor suppressor protein. Recently, we have used spatial coarse graining simulations to analyze the conformational transitions manifest in the p53 recognition of MDM2 and MDMX. These conformational movements are different between MDM2 and MDMX and unveil the presence of conserved and nonconserved interactions in the p53 binding cleft that may be exploited in the design of selective and dual modulators of the oncogenic proteins. In this study, we investigate the conformational profiles of apo‐ and p53‐bound states of MDM2 and MDMX using molecular dynamic simulations along a time scale of 60 ns. The analysis of the trajectories is instrumental to discuss energetical and conformational aspects of p53 recognition and to point out specific key residues whose conformational shifts have crucial roles in affecting the apo‐ and p53‐bound states of MDM2 and MDMX. Among these, in particular, linear discriminant analyses identify diverse conformations of Y99/Y100 (MDMX/MDM2) as markers of the apo‐ and p53‐bound states of the oncogenic proteins. The results of this study shed further light on different p53 recognition in MDM2 and MDMX and may prove useful for the design and identification of new potent and selective synthetic modulators of p53‐MDM2/MDMX interactions. Proteins 2009. © 2009 Wiley‐Liss, Inc.     

A spatiotemporal characterization of the effect of p53 phosphorylation on its interaction with MDM2

The interaction of p53 and MDM2 is modulated by the phosphorylation of p53. This mechanism is key to activating p53, yet its molecular determinants are not fully understood. To study the spatiotemporal characteristics of this molecular process we carried out Brownian dynamics simulations of the interactions of the MDM2 protein with a p53 peptide in its wild type state and when phosphorylated at Thr18 (pThr18) and Ser20 (pSer20). We found that p53 phosphorylation results in concerted changes in the topology of the interaction landscape in the diffusively bound encounter complex domain. These changes hinder phosphorylated p53 peptides from binding to MDM2 well before reaching the binding site. The underlying mechanism appears to involve shift of the peptide away from the vicinity of the MDM2 protein, peptide reorientation, and reduction in peptide residence time relative to wild-type p53 peptide. pThr18 and pSr20 p53 peptides experience reduction in residence times by factors of 13.6 and 37.5 respectively relative to the wild-type p53 peptide, indicating a greater role for Ser20 phosphorylation in abrogating p53 MDM2 interactions. These detailed insights into the effect of phosphorylation on molecular interactions are not available from conventional experimental and theoretical approaches and open up new avenues that incorporate molecular interaction dynamics, for stabilizing p53 against MDM2, which is a major focus of anticancer drug lead development.     

A spatiotemporal characterization of the effect of p53 phosphorylation on its interaction with MDM2

The interaction of p53 and MDM2 is modulated by the phosphorylation of p53. This mechanism is key to activating p53, yet its molecular determinants are not fully understood. To study the spatiotemporal characteristics of this molecular process we carried out Brownian dynamics simulations of the interactions of the MDM2 protein with a p53 peptide in its wild type state and when phosphorylated at Thr18 (pThr18) and Ser20 (pSer20). We found that p53 phosphorylation results in concerted changes in the topology of the interaction landscape in the diffusively bound encounter complex domain. These changes hinder phosphorylated p53 peptides from binding to MDM2 well before reaching the binding site. The underlying mechanism appears to involve shift of the peptide away from the vicinity of the MDM2 protein, peptide reorientation, and reduction in peptide residence time relative to wild-type p53 peptide. pThr18 and pSr20 p53 peptides experience reduction in residence times by factors of 13.6 and 37.5 respectively relative to the wild-type p53 peptide, indicating a greater role for Ser20 phosphorylation in abrogating p53 MDM2 interactions. These detailed insights into the effect of phosphorylation on molecular interactions are not available from conventional experimental and theoretical approaches and open up new avenues that incorporate molecular interaction dynamics, for stabilizing p53 against MDM2, which is a major focus of anticancer drug lead development.     

Computational investigation of the HIV-1 Rev multimerization using molecular dynamics simulations and binding free energy calculations

The HIV Rev protein mediates the nuclear export of viral mRNA, and is thereby essential for the production of late viral proteins in the replication cycle. Rev forms a large organized multimeric protein-protein complex for proper functioning. Recently, the three-dimensional structures of a Rev dimer and tetramer have been resolved and provide the basis for a thorough structural analysis of the binding interaction. Here, molecular dynamics (MD) and binding free energy calculations were performed to elucidate the forces thriving dimerization and higher order multimerization of the Rev protein. It is found that despite the structural differences between each crystal structure, both display a similar behavior according to our calculations. Our analysis based on a molecular mechanics-generalized Born surface area (MM/GBSA) and a configurational entropy approach demonstrates that the higher order multimerization site is much weaker than the dimerization site. In addition, a quantitative hot spot analysis combined with a mutational analysis reveals the most contributing amino acid residues for protein interactions in agreement with experimental results. Additional residues were found in each interface, which are important for the protein interaction. The investigation of the thermodynamics of the Rev multimerization interactions performed here could be a further step in the development of novel antiretrovirals using structure based drug design. Moreover, the variability of the angle between each Rev monomer as measured during the MD simulations suggests a role of the Rev protein in allowing flexibility of the arginine rich domain (ARM) to accommodate RNA binding.     

Probing the interaction mechanism of small molecule inhibitors with matriptase based on molecular dynamics simulation and free energy calculations

Matriptase is a serine protease associated with a wide variety of human tumors and carcinoma progression. Up to now, many promising anti-cancer drugs have been developed. However, the detailed structure–function relationship between inhibitors and matriptase remains elusive. In this work, molecular dynamics simulation and binding free energy studies were performed to investigate the biochemistry behaviors of two class inhibitors binding to matriptase. The binding free energies predicted by MM/GBSA methods are in good agreement with the experimental bioactivities, and the analysis of the individual energy terms suggests that the van der Waals interaction is the major driving force for ligand binding. The key residues responsible for achieving strong binding have been identified by the MM/GBSA free energy decomposition analysis. Especially, Trp215 and Phe99 had an important impact on active site architecture and ligand binding. The results clearly identify the two class inhibitors exist different binding modes. Through summarizing the two different modes, we have mastered some important and favorable interaction patterns between matriptase and inhibitors. Our findings would be helpful for understanding the interaction mechanism between the inhibitor and matriptase and afford important guidance for the rational design of potent matriptase inhibitors.     

Computational investigation of binding mechanism of substituted pyrazinones targeting corticotropin releasing factor-1 receptor deliberated for anti-depressant drug design

In spite of various research investigations towards anti-depressant drug discovery program, no one drug has not yet launched last 20 years. Corticotropin-releasing factor-1 (CRF-1) is one of the most validated targets for the development of antagonists against depression, anxiety and post-traumatic stress disorders. Various research studies suggest that pyrazinone based CRF-1 receptor antagonists were found to be highly potent and efficacious. In this research investigation, we identified the pharmacophore and binding pattern through 2D and 3D-QSAR and molecular docking respectively. Molecular dynamics studies were also performed to explore the binding pattern recognition. We establish the relationship between activity and pharmacophoric features to design new potent compounds. The best 2D-QSAR model was generated through multiple linear regression method with r² value of 0.97 and q² value of 0.89. Also 3D-QSAR model was obtained through k-nearest neighbor molecular field analysis method with q² value of 0.52 and q²_se value of 0.36. Molecular docking and binding energy were also evaluated to define binding patterns and pharmacophoric groups, including (i) hydrogen bond with residue Asp284, Glu305 and (ii) π–π stacking with residue Trp9. Compound 11i has the highest binding affinity compared to reference compounds, so this compound could be a potent drug for stress related disorders. Most of the compounds, including reference compounds were found within acceptable range of physicochemical parameters. These observations could be provided the leads for the design and optimization of novel CRF-1 receptor antagonists.

Communicated by Ramaswamy H. Sarma     

In silico design of novel hERG-neutral sildenafil-like PDE5 inhibitors

Cyclic nucleotide phosphodiesterase enzymes (PDEs) have functions in regulating the levels of intracellular second messengers, 3′, 5′-cyclic adenosine monophosphate (cAMP) and 3′, 5′-cyclic guanosine monophosphate (cGMP), via hydrolysis and decomposing mechanisms in cells. They take essential roles in modulating various cellular activities such as memory and smooth muscle functions. PDE type 5 (PDE5) inhibitors enhance the vasodilatory effects of cGMP in the corpus cavernosum and they are used to treat erectile dysfunction. Patch clamp experiments showed that the IC₅₀ values of the human ether-à-go-go-related gene (hERG1) potassium (K) ion channel blocking affinity of PDE5 inhibitors sildenafil, vardenafil, and tadalafil as 33, 12, and 100 μM, respectively. hERG1 channel is responsible for the regulation of the action potential of human ventricular myocyte by contributing the rapid component of delayed rectifier K⁺ current (I_Kr) component of the cardiac action potential. In this work, interaction patterns and binding affinity predictions of selected PDE5 inhibitors against the hERG1 channel are studied. It is attempted to develop PDE5 inhibitor analogs with lower binding affinity to hERG1 ion channel while keeping their pharmacological activity against their principal target PDE5 using in silico methods. Based on detailed analyses of docking poses and predicted interaction energies, novel analogs of PDE5 inhibitors with lower predicted binding affinity to hERG1 channels without loosing their principal target activity were proposed. Moreover, molecular dynamics (MD) simulations and post-processing MD analyses (i.e. Molecular Mechanics/Generalized Born Surface Area calculations) were performed. Detailed analysis of molecular simulations helped us to better understand the PDE5 inhibitor–target binding interactions in the atomic level. Results of this study can be useful for designing of novel and safe PDE5 inhibitors with enhanced activity and other tailored properties.     

羽扇豆醇介导MDM2-p53通路对胃癌细胞生物学行为的影响

目的:探讨羽扇豆醇介导鼠双微基因2(Mouse double microgene 2,MDM2)-p53通路对胃癌细胞生物学行为的影响及相关机制。方法:对数生长期的胃癌小鼠MFC细胞株随机分为三组。实验1组与实验2组给予10 mg/L和20 mg/L的羽扇豆醇处理,对照组以等体积的1×磷酸盐缓冲液处理。对比三组MFC细胞细胞增殖、凋亡、迁移与侵袭,及MDM2-p53通路蛋白表达。结果:细胞处理后6 h与12 h,实验1组与实验2组的细胞增殖指数、细胞迁移与侵袭指数、MDM2蛋白相对表达水平显著低对于对照组,实验2组也低于实验1组,对比差异都有统计学意义(P<0.05)。细胞处理后6 h与12 h,实验1组与实验2组的细胞凋亡指数、p53蛋白相对表达水平显著高于对照组,实验2组也高于实验1组,对比差异都有统计学意义(P<0.05)。结论:羽扇豆醇能促进胃癌细胞p53蛋白的表达,抑制MDM2蛋白的表达,从而促进细胞凋亡,抑制胃癌的增殖、侵袭与转移,且具有剂量依赖性。     

Identification of potent inhibitors against snake venom metalloproteinase (SVMP) using molecular docking and molecular dynamics studies

Snake venom metalloproteinase (SVMP) (Echis coloratus (Carpet viper) is a multifunctional enzyme that is involved in producing several symptoms that follow a snakebite, such as severe local hemorrhage, nervous system effects and tissue necrosis. Because the three-dimensional (3D) structure of SVMP is not known, models were constructed, and the best model was selected based on its stereo-chemical quality. The stability of the modeled protein was analyzed through molecular dynamics (MD) simulation studies. Structure-based virtual screening was performed, and 15 potential molecules with the highest binding energies were selected. Further analysis was carried out with induced fit docking, Prime/MM–GBSA (ΔG_Bind calculations), quantum-polarized ligand docking, and density functional theory calculations. Further, the stability of the lead molecules in the SVMP-active site was examined using MD simulation. The results showed that the selected lead molecules were highly stable in the active site of SVMP. Hence, these molecules could potentially be selective inhibitors of SVMP. These lead molecules can be experimentally validated, and their backbone structural scaffold could serve as building blocks in designing drug-like molecules for snake antivenom.     

Fluorescence polarization assay and inhibitor design for MDM2/p53 interaction

MDM2 is an important negative regulator of the tumor suppressor protein p53 which regulates the expression of many genes including MDM2. The delicate balance of this autoregulatory loop is crucial for the maintenance of the genome and control of the cell cycle and apoptosis. MDM2 hyperactivity, due to amplification/overexpression or mutational inactivation of the ARF locus, inhibits the function of wild-type p53 and can lead to the development of a wide variety of cancers. Thus, the development of anti-MDM2 therapies may restore normal p53 function in tumor cells and induce growth suppression and apoptosis. We report here a novel high-throughput fluorescence polarization binding assay and its application in rank ordering small-molecule inhibitors that block the binding of MDM2 to a p53-derived fluorescent peptide.     

Molecular dynamics simulations and MM/GBSA methods to investigate binding mechanisms of aminomethylpyrimidine inhibitors with DPP-IV

The aminomethylpyrimidines were investigated as a novel class of DPP-IV inhibitors. In this Letter, the binding mechanisms of how slight change of substitution or position influences the binding affinity of five representative analogs was investigated by molecular dynamics simulation, free energy calculations and energy decomposition analysis. The conserved hydrogen bonds with Glu205 and Glu206 slightly favor the inhibitor binding; the van der Waals interactions, especially the two key contacts with Tyr547 and Tyr666, dominate in the binding free energy and play a crucial role on distinguishing the high active inhibitors from the low ones.