A molecular dynamics approach to ligand-receptor interaction in the aspirin-human serum albumin complex

Investigation of the binding behavior of bioactive 7-methoxyflavone to human serum albumin by coupling multi-spectroscopic with computational approaches

The natural flavonoids with bioactivity as secondary plant metabolites are mostly found in fruits, vegetables, tea and herbs, the distribution and bioavailability of which in vivo depends on the interaction and successive binding with carrier proteins in the systemic circulation. In this paper, the binding behavior of bioactive 7-methoxyflavone (7-MF) with human serum albumin (HSA) was studied with the aid of the combination of multi-spectroscopic methods, molecular docking and molecular dynamic simulation. The results of multi-spectroscopic experiments revealed that 7-MF interacted with HSA predominantly via fluorescence static quenching and the microenvironment around the fluorophore Trp residues in HSA became more hydrophilicity with the binding of 7-MF. Thermodynamic analysis demonstrated that hydrogen bonds and van der Waals forces played a dominant role in stabilizing the HSA-7-MF complex. Moreover, the docking experiment and molecular dynamic simulation further confirmed that 7-MF could enter the active cavity of HSA and caused more stable conformation and change of secondary structure of HSA through forming hydrogen bond. The exploration of the mechanism of 7-MF binding to HSA lights a new avenue to understand the stability, transport and distribution of 7-MF and 7-MF may hold great potential to be extended as a promising alternative of dietary supplements or pharmaceutical agents.

Synthesis, evaluation, molecular dynamics simulation and targets identification of novel pyrazole-containing imide derivatives

A new series of novel pyrazole-containing imide derivatives were synthesized and evaluated for their anticancer activities against A-549, Bel7402, and HCT-8 cell lines. Among these compounds A2, A4, A11 and A14 possessed high inhibition activity against A-549 cell lines with IC₅₀ values at 4.91, 3.22, 27.43 and 18.14 μM, respectively, better than that of 5-fluorouracil (IC₅₀=59.27 μM). A2, A4, and A11 also exhibited significant inhibitory activity towards HCT-8 and Bel7402 cell lines. Interestingly, the Heat Shock Protein 90α (Hsp90α, PDB ID: 1UYK) was found to be the potential drug target of these synthesized compounds with the aid of PharmMapper server (http://lilab.ecust.edu.cn/pharmmapper/) and docking module of Schrödinger (Maestro 10.2). Additionally, molecular dynamics simulation was performed out to explore the most likely binding mode of compound A2 with Hsp90α.Communicated by Ramaswamy H. Sarma.

Identification of the potential dual inhibitor of protein tyrosine phosphatase sigma and leukocyte common antigen-related phosphatase by virtual screen, molecular dynamic simulations and post-analysis

Owing to their inhibitory role in regulating oligodendrocyte differentiation and apoptosis, protein tyrosine phosphatase sigma (PTPσ) and leukocyte common antigen-related phosphatase (LAR) play a crucial potential role in treating spinal cord injury (SCI) disease. In this research, the computer aided drug design (CADD) methods were applied to discover the potential dual-target drug involving virtual screen, molecular docking and molecular dynamic simulation. Initially, the top 20 compounds with higher docking score than the positive controls (ZINC13749892, ZINC14516161) were virtually screened out from NCI and ZINC databases, and then were submitted in ADMET to predict their drug properties. Among these potential compounds, ZINC72417086 showed a higher docking score and satisfied Lipinski's rule of five. In addition, the post-analysis demonstrated that when ZINC72417086 bound to PTPσ and LAR, it could stable proteins conformations and destroy the residues interactions between P-loop and other loop regions in active pocket. Meanwhile, residue ARG1595 and ARG1528 could play a crucial role in in the inhibition of PTPσ and LAR, respectively. This research offered a novel approach for rapid discovery of dual-target leads compounds to treat SCI.Communicated by Ramaswamy H. Sarma.

Exploring the cause of the inhibitor 4AX attaching to binding site disrupting protein tyrosine phosphatase 4A1 trimerization by molecular dynamic simulation

Ectopic overexpression of protein tyrosine phosphatase of liver regeneration-1 (PTP4A1, also called PRL-1) markedly enhanced hepatocellular carcinoma (HCC) cells migration and invasion. The PTP4A1 trimerization played a vital role in mediating cell proliferation and motility. Biochemical and structural studies have proved that the compound 4AX, a well-known inhibitor for PRL1, directly binds to the PTP4A1 trimer interface and obstructs trimer formation of PTP4A1. However, the molecular basis of the ligand-4AX inhibition on PTP4A1 trimer conformations remains unclear. In this study, the docking analysis and the molecular dynamics simulation (MD simulation) study were performed to investigate how the molecule binding at each interface disrupted the trimer formation. The results suggested that the ligand-4AX attaching to the binding site changed the conformation of A:Q131, A:Q135 in the AC interface, C:R18, C:P96 in the CA interface and B:Q131 in the BA interface, leading to the weak interactions between subunits and thus resulting in the disruption of the PTP4A1 trimerization.

Combined comparative molecular field analysis, comparative molecular similarity indices analysis, molecular docking and molecular dynamics studies of histone deacetylase 6 inhibitors

Human histone deacetylase isoform 6 (HDAC6) has been shown to have an immense role in cell motility and aggresome formation and is being an attractive selective target for the treatment of multiple tumour types and neurodegenerative conditions. The discovery of selective HDAC6 inhibitors with new chemical functionalities is therefore of utmost interest to researchers. In order to examine the structural requirements for HDAC6-specific inhibitors and to derive predictive model which can be used for designing new selective HDAC6 inhibitors, a three-dimensional quantitative structure-activity relationship study was carried out on a diverse set of ligands using common feature-based pharmacophore alignment followed by employing comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA) techniques. The models displayed high correlation of 0.978 and 0.991 for best CoMFA and CoMSIA models, respectively, and a good statistical significance. The model could be used for predicting activities of the test set compounds as well as for deriving useful information regarding steric, electrostatic, hydrophobic properties of the molecules used in this study. Further, the training and test set molecules were docked into the HDAC6 binding site and molecular dynamics was carried out to suggest structural requirements for design of new inhibitors.

Telomerase Inhibition by a New Synthetic Derivative of the Aporphine Alkaloid Boldine

Telomerase, the enzyme responsible for cell immortality, is an important target in anti-cancer drug discovery. Boldine, an abundant aporphine alkaloid of Peumus boldus, is known to inhibit telomerase at non-toxic concentrations. Cytotoxicity of N-benzylsecoboldine hydrochloride (BSB), a synthetic derivative of boldine, was determined using the MTT method in MCF7 and MDA-MB231 cells. Aliquots of cell lysates were incubated with various concentrations of BSB in qTRAP (quantitative telomere repeat amplification protocol)-ligand experiments before substrate elongation by telomerase or amplification by hot-start Taq polymerase. The crystal structure of TERT, the catalytic subunit of telomerase from Tribolium castaneum, was used for docking and molecular dynamics analysis. The qTRAP-ligand data gave an IC₅₀ value of about 0.17 ± 0.1 µM for BSB, roughly 400 times stronger than boldine, while the LD₅₀ in the cytotoxicity assays were 12.5 and 21.88 µM, respectively, in cells treated for 48 h. Although both compounds interacted well with the active site, MD analysis suggests a second binding site with which BSB interacts via two hydrogen bonds, much more strongly than boldine. Theoretical analyses also evaluated the IC₅₀ for BSB as submicromolar. BSB, with greater hydrophobicity and flexibility than boldine, represents a promising structure to inhibit telomerase at non-toxic concentrations.

A tunable coarse-grained model for ligand-receptor interaction

Cell-surface receptors are the most common target for therapeutic drugs. The design and optimization of next generation synthetic drugs require a detailed understanding of the interaction with their corresponding receptors. Mathematical approximations to study ligand-receptor systems based on reaction kinetics strongly simplify the spatial constraints of the interaction, while full atomistic ligand-receptor models do not allow for a statistical many-particle analysis, due to their high computational requirements. Here we present a generic coarse-grained model for ligand-receptor systems that accounts for the essential spatial characteristics of the interaction, while allowing statistical analysis. The model captures the main features of ligand-receptor kinetics, such as diffusion dependence of affinity and dissociation rates. Our model is used to characterize chimeric compounds, designed to take advantage of the receptor over-expression phenotype of certain diseases to selectively target unhealthy cells. Molecular dynamics simulations of chimeric ligands are used to study how selectivity can be optimized based on receptor abundance, ligand-receptor affinity and length of the linker between both ligand subunits. Overall, this coarse-grained model is a useful approximation in the study of systems with complex ligand-receptor interactions or spatial constraints.

The current importance of cell surface receptors as primary targets for drug treatment explains the increasing interest in a mathematical and quantitative description of the process of ligand-receptor interaction. Recently, a new generation of synthetic chimeric ligands has been developed to selectively target unhealthy cells, without harming healthy tissue. To understand these and other types of complex ligand-receptor systems, conventional chemical interaction models often rely on simplifications and assumptions about the spatial characteristics of the system, while full atomistic molecular dynamics simulations are too computationally demanding to perform many particle statistical analysis. In this paper, we describe a novel approach to model the interaction between ligands and receptors based on a coarse grained approximation that includes explicitly both spatial and kinetic details of the interaction, while allowing many-particle simulations and therefore, statistical analysis at biologically relevant time scales. The model is used to study the binding properties of generic chimeric ligands to understand how cell specificity is achieved, its dependence on receptor concentration and the influence of the distance between subunits of the chimera. Overall, this approach proves optimal to study other ligand-receptor systems with complex spatial regulation, such as receptor clustering, multimerization and multivalent asymmetric ligand binding.