Quantum and molecular dynamics study for binding of macrocyclic inhibitors to human alpha-thrombin

Molecular dynamics simulations followed by quantum mechanical calculation and Molecular Mechanics Poisson-Boltzmann Surface Area (MM-PBSA) analysis have been carried out to study binding of proline- and pyrazinone-based macrocyclic inhibitors (L86 and T76) to human alpha-thrombin. Detailed binding interaction energies between these inhibitors and individual protein fragments are calculated using DFT method based on a new quantum mechanical approach for computing protein-ligand interaction energy. The analysis of detailed interaction energies provides insight on the protein-ligand binding mechanism. Study shows that T76 and L86 bind to thrombin in a very similar "inhibition mode" except that T76 has relatively weaker binding interaction with Glu(217). The analysis from quantum calculation of binding interaction is consistent with the MM-PBSA calculation of binding free energy, and the calculated free energies for L86/T76-thrombin binding agree well with the experimental data.     

Novel tacrine-pyridinium hybrid reactivators of organophosphorus-inhibited acetylcholinesterase: Synthesis,molecular docking,and in vitro reactivation study

First-line medical treatment against nerve agents consists of co-administration of anticholinergic agents and oxime reactivators, which reactivate inhibited AChE. Pralidoxime, a commonly used oxime reactivator, is effective against some nerve agents but not against others; thus, new oxime reactivators are needed. Novel tacrine-pyridinium hybrid reactivators in which 4-pyridinealdoxime derivatives are connected to tacrine moieties by linear carbon chains of different lengths (C2–C7) were prepared (Scheme 1, 5a–f). Their binding affinities to electric eel AChE were tested because oximes can inhibit free AChE, and the highest AChE activity (95%, 92%, and 90%) was observed at 1?μM concentrations of the oximes (5a, 5b, and 5c, respectively). Based on their inhibitory affinities towards free AChE, 1?μM concentrations of the oxime derivatives (5) were used to examine reactivation of paraoxon-inhibited AChE. Reactivation ability increased as the carbon linker chains lengthened (n?=?2–5), and 5c and 5d showed remarkable reactivation ability (41%) compared to that of 2-PAM (16%) and HI-6 (4%) against paraoxon-inhibited electric eel AChE at 1?μM concentrations. Molecular docking simulation showed that the most stable binding free energy was observed in 5c at 73.79?kcal?mol^?1, and the binding mode of 5c is acceptable for the oxygen atom of oximate to attack the phosphorus atom of paraoxon and reactivate paraoxon-inhibited eel AChE model structure.     

Insight into the binding modes and inhibition mechanisms of adamantyl-based 1,3-disubstituted urea inhibitors in the active site of the human soluble epoxide hydrolase

Soluble epoxide hydrolase (sEH) is a promising new target for treating hypertension and inflammation. Considerable efforts have been devoted to develop novel inhibitors. In this study, the binding modes and interaction mechanisms of a series of adamantyl-based 1,3-disubstituted urea inhibitors were investigated by molecular docking, molecular dynamics simulations, binding free energy calculations, and binding energy decomposition analysis. Based on binding affinity, the most favorable binding mode was determined for each inhibitor. The calculation results indicate that the total binding free energy (ΔG_TOT, the sum of enthalpy ΔG_MM-GB/SA, and entropy ?TΔS) presents a good correlation with the experimental inhibitory activity (IC₅₀, r²?=?.99). The van der Waals energy contributes most to the total binding free energy (ΔG_TOT). A detailed discussion on the interactions between inhibitors and those residues located in the active pocket is made based on hydrogen bond and binding modes analysis. According to binding energy decomposition, the residues Asp333 and Trp334 contribute the most to binding free energy in all systems. Furthermore, Hip523 plays a major role in determining this class of inhibitor-binding orientations. Combined with the results of hydrogen bond analysis and binding free energy, we believe that the conserved hydrogen bonds play a role only in anchoring the inhibitors to the exact site for binding and the number of hydrogen bonds may not directly relate to the binding free energy. The results we obtained will provide valuable information for the design of high potency sEH inhibitors.     

Mapping the HSP90 binding region of the glucocorticoid receptor

In animal cells, unliganded steroid receptors are complexed with a 90-kDa heat shock protein, HSP90; hormone binding by the receptor leads to the release of HSP90. We found that the 795-amino acid rat glucocorticoid receptor protein formed oligomeric complexes in vitro upon synthesis in rabbit reticulocyte lysates; these oligomers also dissociated in the presence of hormone. Similar complexes formed when X795, a receptor derivative containing only the C-terminal half (amino acids 407-795) of the protein, was translated in vitro. Moreover, X795 was co-immunoadsorbed from the reticulocyte lysates together with HSP90 by three different anti-HSP90 monoclonal antibodies, indicating that the in vitro translated receptor binds HSP90 and that the interaction occurs within the C-terminal half of the receptor. To localize the HSP90 binding region in greater detail, various deletion mutants of X795 were translated in vitro and assayed for oligomer formation and for co-immunoadsorption with HSP90. The results indicated that HSP90 interacted with the receptor within a subregion of the hormone binding domain, between amino acids 568 and 616. These findings are consistent with the proposal that HSP90 may participate in the mechanism of signal transduction by steroid receptors.     

Fatty acid binding to serum albumin: Molecular simulation approaches

Background

Binding affinity for human serum albumin (HSA) is one of the most important factors affecting the distribution and free blood concentration of many ligands. The effect of fatty acids (FAs) on HSA-ligand binding has long been studied. Since the elucidation of the 3-dimensional structure of HSA, molecular simulation approaches have been applied to studies of the structure–function relationship of HSA–FA binding.

Scope of review

We review current insights into the effects of FA binding on HSA, focusing on the biophysical insights obtained using molecular simulation approaches such as docking, molecular dynamics (MD), and binding free energy calculations.

Major conclusions

Possible conformational changes on binding of FA molecules to HSA have been observed through MD simulations. High- and low-affinity FA-binding sites on HSA have been identified based on binding free energy calculations. The relationship between the warfarin binding affinity of HSA and FA molecules has been clarified based on the results of simulations of multi-site FA binding that cannot be experimentally observed.

General significance

Molecular simulation approaches have great potentials to provide detailed biophysical insights into HSA as well as the effects of the binding of FAs or other ligands to HSA. This article is part of a Special Issue entitled Serum Albumin.     

Unusual binding mode of scorpion toxin BmKTX onto potassium channels relies on its distribution of acidic residues

Besides classical scorpion toxin–potassium channel binding modes, novel modes remain unknown. Here, we report a novel binding mode of native toxin BmKTX towards Kv1.3 channel. The combined experimental and computational data indicated that BmKTX-D33H analog used the classical anti-parallel β-sheet domain as the channel-interacting interface together with the conserved channel pore-blocking Lys²⁶. However, the wild-type BmKTX was found to use Arg²³ rather than Lys²⁶ as the new pore-blocking residue, and mainly adopt the turn motif between the α-helix and antiparallel β-sheet domains to recognize K_v1.3 channel. Together, these findings not only reveal that scorpion toxin–potassium channel interaction modes are more diverse than thought, but also highlight the functional role of toxin acidic residues in mediating diverse toxin–potassium channel binding modes.     

S100A4: a novel partner for heat shock protein 47 in antler stem cells and insight into the calcium ion-induced conformational changes

Abstract

S100A4 is a multiple-function protein highly expressed in tumor or stem cells. We found S100A4 was a novel protein partner for heat shock protein 47 (HSP47) in deer antlerogenic periosteum cells (AP cells), indicating that S100A4 could bind with HSP47. S100A4 had both calcium-dependent and calcium-independent patterns (labeled as SCd and SCi, respectively) to execute different biological activities. Homology models of HSP47, SCd and SCi were constructed. HSP47:collagen model, HSP47:collagen I-V, HSP47:SCd and HSP47:SCi complexes were built using ZDOCK software. Together with free SCd and SCi, 200?ns molecular dynamic (MD) simulations were performed to analyze binding free energies and SCi/SCd conformational changes. The energetic results showed that SCi had the strongest affinity to HSP47, and followed by collagens. SCd had little interaction with HSP47. Decomposition energy results showed that collagen model interacted with HSP47 mainly though neutral amino acids. When SCi bound with HSP47, the majority of mediated amino acids were charged. These results indicated that SCi could compete with collagen on the binding site of HSP47. Root mean square fluctuation (RMSF) values and cross-correlation matrices of principal component analysis (PCA) were calculated to evaluate the SCi/SCd structural variation during MD simulation. Both HSP47 and Ca²⁺ could stabilize the conformation of SCi/SCd. The loops interacting with Ca²⁺s and linking the two EF-hand motifs were impacted particularly. The relative moving directions of α-helices in EF-hands were distinct by the binding effect of HSP47 and Ca²⁺. We found that SCi may regulate the differentiation of AP cells by disturbing the interaction between HSP47 and collagen.

Communicated by Ramaswamy H. Sarma     

Crystal structure of the single-stranded RNA binding protein HutP from Geobacillus thermodenitrificans

RNA binding proteins control gene expression by the attenuation/antitermination mechanism. HutP is an RNA binding antitermination protein. It regulates the expression of hut operon when it binds with RNA by modulating the secondary structure of single-stranded hut mRNA. HutP necessitates the presence of l-histidine and divalent metal ion to bind with RNA. Herein, we report the crystal structures of ternary complex (HutP–l-histidine–Mg²⁺) and EDTA (0.5 M) treated ternary complex (HutP–l-histidine–Mg²⁺), solved at 1.9 Å and 2.5 Å resolutions, respectively, from Geobacillus thermodenitrificans. The addition of 0.5 M EDTA does not affect the overall metal-ion mediated ternary complex structure and however, the metal ions at the non-specific binding sites are chelated, as evidenced from the results of structural features.     

Blinded predictions of binding modes and energies of HSP90-α ligands for the 2015 D3R grand challenge

In the framework of the 2015 D3R inaugural grand challenge, blind binding pose and affinity predictions were performed for a set of 180 ligands of the Heat Shock Protein HSP90-α protein, a relevant cancer target. Spectral clustering was used to rapidly identify alternative binding site conformations in publicly available crystallographic HSP90-α structures. Subsequently, multiple docking and scoring protocols employing the software Autodock Vina and rDock were applied to predict binding modes and rank order ligands. Alchemical free energy calculations were performed with the software FESetup and Sire/OpenMM to predict binding affinities for three congeneric series subsets. Some of the protocols used here were ranked among the top submissions according to most of the evaluation metrics. Docking performance was excellent, but the scoring results were disappointing. A critical assessment of the results is reported, as well as suggestions for future similar competitions.     

Ligand binding by antibody IgE Lb4: assessment of binding site preferences using microcalorimetry, docking, and free energy simulations.

Antibody IgE Lb4 interacts favorably with a large number of different compounds. To improve the current understanding of the structural basis of this vast cross-reactivity, the binding of three dinitrophenyl (DNP) amino acids (DNP-alanine, DNP-glycine, and DNP-serine) is investigated in detail by means of docking and molecular dynamics free energy simulations. Experimental binding energies obtained by isothermal titration microcalorimetry are used to judge the results of the computational studies. For all three ligands, the docking procedure proposes two plausible subsites within the binding region formed by the antibody CDR loops. By subsequent molecular dynamics simulations and calculations of relative free energies of binding, one of these subsites, a tyrosine-surrounded pocket, is revealed as the preferred point of complexation. For this subsite, results consistent with experimental observations are obtained; DNP-glycine is found to bind better than DNP-serine, and this, in turn, is found to bind better than DNP-alanine. The suggested binding mode makes it possible to explain both the moderate binding affinity and the differences in binding energy among the three ligands.     

Calculation of substrate binding affinities for a bacterial GH78 rhamnosidase through molecular dynamics simulations

Ram2 from Pediococcus acidilactici is a rhamnosidase from the glycoside hydrolase family 78. It shows remarkable selectivity for rutinose rather than para-nitrophenyl-alpha-l-rhamnopyranoside (p-NPR). Molecular dynamics simulations were performed using a homology model of this enzyme, in complex with both substrates. Free energy calculations lead to predicted binding affinities of −34.4 and −30.6 kJ mol⁻¹ respectively, agreeing well with an experimentally estimated relative free energy of 5.4 kJ mol⁻¹. Further, the most relevant binding poses could be determined. While p-NPR preferably orients its rhamnose moiety toward the active site, rutinose interacts most strongly with its glucose moiety. A detailed hydrogen bond analysis confirms previously implicated residues in the active site (Asp217, Asp222, Trp226, Asp229 and Glu488) and quantifies the importance of individual residues for the binding. The most important amino acids are Asp229 and Phe339 which are involved in many interactions during the simulations. While Phe339 was observed in more simulations, Asp229 was involved in more persistent interactions (forming an average of at least 2 hydrogen bonds during the simulation). These analyses directly suggest mutations that could be used in a further experimental characterization of the enzyme. This study shows once more the strength of computer simulations to rationalize and guide experiments at an atomic level.     

Calmodulin-regulated binding of the 90-kDa heat shock protein to actin filaments

We have found that the 90-kDa heat shock protein (HSP90) prepared from a mouse lymphoma exists in homodimeric form under physiological conditions and has the ability to bind to F-actin (Koyasu, S., Nishida, E., Kadowaki, T., Matsuzaki, F., Iida, K., Harada, F., Kasuga, M., Sakai, H., and Yahara, I. (1986) Proc. Natl. Acad. Sci. U.S.A., in press). Here we show that calmodulin regulates the binding of HSP90 to F-actin in a Ca2+-dependent manner. The binding of HSP90 to F-actin occurred optimally under physiological solution conditions, i.e. in 2 mM MgCl2 + 100 mM KCl. The binding was saturable in a molar ratio of about 1 HSP90 (dimer) to 10 actins. HSP90 was dissociated from F-actin by the binding of tropomyosin to F-actin. Calmodulin was found to inhibit the binding of HSP90 to F-actin in a Ca2+-dependent manner. Moreover, the equilibrium gel filtration demonstrated that calmodulin binds to HSP90 in the presence of Ca2+, but not in the absence of Ca2+. These data indicate that HSP90 complexed with Ca2+-calmodulin is unable to bind to F-actin. Ca2+-dependent interaction of HSP90 with calmodulin as well as calmodulin-regulated binding of HSP90 to F-actin revealed here may provide new insight into the function of HSP90 and the regulation of actin structure in cells.     

Macrocyclic inhibitors of the malarial aspartic proteases plasmepsin I, II, and IV

The first macrocyclic inhibitor of the Plasmodium falciparum aspartic proteases plasmepsin I, II, and IV with considerable selectivity over the human aspartic protease cathepsin D has been identified. A series of macrocyclic compounds were designed and synthesized. Cyclizations were accomplished using ring-closing metathesis with the second generation Grubbs catalyst. These compounds contain either a 13-membered or a 16-membered macrocycle and incorporate a 1,2-dihydroxyethylene as transition state mimicking unit. The binding mode of this new class of compounds was predicted with automated docking and molecular dynamics simulations, with an estimation of the binding affinities through the linear interaction energy (LIE) method.     

Radiosynthesis,biological evaluation and preliminary microPET study of 18F-labeled 5-resorcinolic triazolone derivative based on ganetespib targeting HSP90

Heat-shock protein 90 (HSP90) is a molecular chaperone that activates oncogenic transformation in several solid tumors, including lung and breast cancers. Ganetespib, a most promising candidate among several HSP90 inhibitors under clinical trials, has entered Phase III clinical trials for cancer therapy. Despite numerous evidences validating HSP90 as a target of anticancer, there are few studies on PET agents targeting oncogenic HSP90. In this study, we synthesized and biologically evaluated a novel ¹⁸F-labeled 5-resorcinolic triazolone derivative (1, [¹⁸F]PTP-Ganetespib) based on ganetespib. [¹⁸F]PTP-Ganetespib was labeled by click chemistry of Ganetespib-PEG-Alkyne (10) and [¹⁸F]PEG-N₃ (11) with 37.3?±?5.11% of radiochemical yield and 99.7?±?0.09% of radiochemical purity. [¹⁸F]PTP-Ganetespib showed proper LogP (0.96?±?0.06) and good stability in human serum over 97% for 2?h. [¹⁸F]PTP-Ganetespib showed high uptakes in breast cancer cells containing triple negative breast cancer (TNBC) MDA-MB-231 and Her2-negative MCF-7 cells, which are target breast cancer cell lines of HSP90 inhibitor, ganetespib, as an anticancer. Blocking of HSP90 by the pretreatment of ganetespib exhibited significantly decreased accumulation of [¹⁸F]PTP-Ganetespib in MDA-MB-231 and MCF-7 cells, indicating the specific binding of [¹⁸F]PTP-Ganetespib to MDA-MB-231 and MCF-7 cells with high HSP90 expression. In the biodistribution and microPET imaging studies, the initial uptake into tumor was weaker than in other thoracic and abdominal organs, but [¹⁸F]PTP-Ganetespib was retained relatively longer in the tumor than other organs. The uptake of [¹⁸F]PTP-Ganetespib in tumors was not sufficient for further development as a tumor-specific PET imaging agent by itself, but this preliminary PET imaging study of [¹⁸F]PTP-Ganetespib can be basis for developing new PET imaging agents based on HSP90 inhibitor, ganetespib.     

Investigation on the binding mechanism of loratinib with the c-ros oncogene 1 (ROS1) receptor tyrosine kinase via molecular dynamics simulation and binding free energy calculations

The c-ros oncogene 1 (ROS1) has proven to be an important cancer target for the treatment of various human cancers. The anaplastic lymphoma kinase inhibitor crizotinib has been granted approval for the treatment of patients with ROS1 positive metastatic non-small-cell lung cancer by the Food and Drug Administration on 2016. However, serious resistance due to the secondary mutation of glycine 2032 to arginine (G2032R) was developed in clinical studies. Loratinib (PF-06463922), a macrocyclic analog of crizotinib, showed significantly improved inhibitory activity against wild–type (WT) ROS1 and ROS1^G2032R mutant. To provide insights into the inhibition mechanism, molecular dynamics simulations and free energy calculations were carried out for the complexes of loratinib with WT and G2032R mutated ROS1. The apo-ROS1^WT and apo-ROS1^G2032R systems showed similar RMSF distributions, while ROS1^G2032R-loratinib showed significantly higher than that of WT ROS1-loratinib, which revealed that the binding of loratinib to ROS1^G2032R significantly interfered the ?uctuation of protein. Calculations of binding free energies indicate that G2032R mutation significantly reduces the binding affinity of loratinib for ROS1, which arose mostly from the increase of conformation entropy and the decrease of solvation energy. Furthermore, detailed per-residue binding free energies highlighted the increased and decreased contributions of some residues in the G2032R mutated systems. The present study revealed the detailed inhibitory mechanism of loratinib as potent WT and G2032R mutated ROS1 inhibitor, which was expected to provide a basis for rational drug design.     

Effect of the 90 kDa heat shock protein, HSP90, on glucocorticoid receptor binding to DNA-cellulose

Glucocorticoid receptors in the IM-9 human lymphoblastoid cell line were affinity labeled with [3H]dexamethasone 21-mesylate and activated to a DNA-binding form by filtration through a Bio-Gel A-1.5m column. The 90 kDa heat shock protein, HSP90, was identified by labeling IM-9 cells with 35S-methionine at both 37 degrees C and 42 degrees C and purified to near homogeneity by sequential chromatography through DE52 and hydroxyapatite. Addition of purified HSP90 to activated, affinity labeled glucocorticoid receptors in a molecular ratio of 16 to 1 inhibited the binding of the receptors to DNA-cellulose. HSP90 did not affect the binding of other proteins to DNA-cellulose, indicating that the inhibitory effect of HSP90 was specific for the glucocorticoid receptor. These results suggest that HSP90 may associate with the glucocorticoid receptor, masking its DNA-binding site and thereby inhibiting receptor interaction with DNA.     

Focal adhesion protein abnormalities in myelodysplastic mesenchymal stromal cells

Direct cell–cell contact between haematopoietic progenitor cells (HPCs) and their cellular microenvironment is essential to maintain ‘stemness’. In cancer biology, focal adhesion (FA) proteins are involved in survival signal transduction in a wide variety of human tumours. To define the role of FA proteins in the haematopoietic microenvironment of myelodysplastic syndromes (MDS), CD73-positive mesenchymal stromal cells (MSCs) were immunostained for paxillin, pFAK [Y³⁹⁷], and HSP90α/β and p130CAS, and analysed for reactivity, intensity and cellular localisation. Immunofluorescence microscopy allowed us to identify qualitative and quantitative differences, and subcellular localisation analysis revealed that in pathological MSCs, paxillin, pFAK [Y³⁹⁷], and HSP90α/β formed nuclear molecular complexes. Increased expression of paxillin, pFAK [Y³⁹⁷], and HSP90α/β and enhanced nuclear co-localisation of these proteins correlated with a consistent proliferative advantage in MSCs from patients with refractory anaemia with excess blasts (RAEB) and negatively impacted clonogenicity of HPCs. These results suggest that signalling via FA proteins could be implicated in HPC–MSC interactions. Further, because FAK is an HSP90α/β client protein, these results suggest the utility of HSP90α/β inhibition as a target for adjuvant therapy for myelodysplasia.     

A single mutation within a Ca binding loop increases proteolytic activity,thermal stability,and surfactant stability

We improved the enzymatic properties of the oxidatively stable alkaline serine protease KP-43 through protein engineering to make it more suitable for use in laundry detergents. To enhance proteolytic activity, the gene encoding KP-43 was mutagenized by error-prone PCR. Screening identified a Tyr195Cys mutant enzyme that exhibited increased specific activity toward casein between pH 7 and 11. At pH 10, the mutant displayed 1.3-fold higher specific activity for casein compared to the wild-type enzyme, but the activity of the mutant was essentially unchanged toward several synthetic peptides. Furthermore, the Tyr195Cys mutation significantly increased thermal stability and surfactant stability of the enzyme under oxidizing conditions. Examination of the crystal structure of KP-43 revealed that Tyr195 is a solvent exposed residue that forms part of a flexible loop that binds a Ca^2 + ion. This residue lies 15–20 Å away from the residues comprising the catalytic triad of the enzyme. These results suggest that the substitution at position 195 does not alter the structure of the active center, but instead may affect a substrate–enzyme interaction. We propose that the Tyr195Cys mutation enhances the interaction with Ca^2 + and affects the packing of the Ca^2 + binding loop, consequently increasing protein stability. The simultaneously increased proteolytic activity, thermal stability, and surfactant stability of the Tyr195Cys mutant enzyme make the protein an ideal candidate for laundry detergent application.     

RNA targeting by DNA binding drugs: Structural,conformational and energetic aspects of the binding of quinacrine and DAPI to A-form and H-form of poly(rC)·poly(rG)

A key step in the rational design of new RNA binding small molecules necessitates a complete elucidation of the molecular aspects of the binding of existing molecules to RNA structures. This work focuses towards the understanding of the interaction of a DNA intercalator, quinacrine and a minor groove binder 4′,6-diamidino-2-phenylindole (DAPI) with the right handed Watson–Crick base paired A-form and the left-handed Hoogsteen base paired H^L-form of poly(rC)·poly(rG) evaluated by multifaceted spectroscopic and viscometric techniques. The energetics of their interaction has also been elucidated by isothermal titration calorimetry. Results of this study converge to suggest that (i) quinacrine intercalates to both A-form and H^L-form of poly(rC)·poly(rG); (ii) DAPI shows both intercalative and groove-binding modes to the A-form of the RNA but binds by intercalative mode to the H^L-form. Isothermal calorimetric patterns of quinacrine binding to both the forms of RNA and of DAPI binding to the H^L-form are indicative of single binding while the binding of DAPI to the A-form reveals two kinds of binding. The binding of both the drugs to both conformations of RNA is exothermic; while the binding of quinacrine to both conformations and DAPI to the A-form (first site) is entropy driven, the binding of DAPI to the second site of A-form and H^L-conformation is enthalpy driven. Temperature dependence of the binding enthalpy revealed that the RNA–ligand interaction reactions are accompanied by small heat capacity changes that are nonetheless significant. We conclude that the binding affinity characteristics and energetics of interaction of these DNA binding molecules to the RNA conformations are significantly different and may serve as data for the development of effective structure selective RNA-based antiviral drugs.