key: cord-0005479-bsycsuoa authors: Cele, Favourite N.; Kumalo, Hezekiel; Soliman, Mahmoud E. S. title: Mechanism of Inhibition of Hsp90 Dimerization by Gyrase B Inhibitor Coumermycin A1 (C–A1) Revealed by Molecular Dynamics Simulations and Thermodynamic Calculations date: 2016-07-04 journal: Cell Biochem Biophys DOI: 10.1007/s12013-016-0743-8 sha: 33e73a8f5255ef56f603d50190e854539930220c doc_id: 5479 cord_uid: bsycsuoa Heat shock protein (Hsp) 90 an emerging and attracting target in the anti-HIV drug discovery process due to the key role it plays in the pathogenicity of HIV-1 virus. In this research study, long-range all-atom molecular dynamics simulations were engaged for the bound and the unbound proteins to enhance the understanding of the molecular mechanisms of the Hsp90 dimerization and inhibition. Results evidently showed that coumermycin A1 (C–A1), a recently discovered Hsp90 inhibitor, binds at the dimer’s active site of the Hsp90 protein and leads to a substantial parting between dimeric opposed residues, which include Arg591.B, Lys594.A, Ser663.A, Thr653.B, Ala665.A, Thr649.B, Leu646.B and Asn669.A. Significant differences in magnitudes were observed in radius of gyration, root-mean-square deviation and root-mean-square fluctuation, which confirms a reasonably more flexible state in the apo conformation associated with it dimerization. In contrast, the bound conformer of Hsp90 showed less flexibility. This visibly highpoints the inhibition process resulting from the binding of the ligand. These findings were further validated by principal component analysis. We believe that the detailed dynamic analyses of Hsp90 presented in this study, would give an imperative insight and better understanding to the function and mechanisms of inhibition. Furthermore, information obtained from the binding mode of the inhibitor would be of great assistance in the design of more potent inhibitors against the HIV target Hsp90. ELECTRONIC SUPPLEMENTARY MATERIAL: The online version of this article (doi:10.1007/s12013-016-0743-8) contains supplementary material, which is available to authorized users. Acquired immunodeficiency syndrome (AIDS) is the most challenging disease outbreak thus far and is caused by the human immunodeficiency virus (HIV) [21] . HIV/AIDS affected over 35 million people globally, and Sub-Saharan Africa is the most affected region with approximately 68 % people affected (https://en.wikipedia.org/wiki/Epidemio Electronic supplementary material The online version of this article (doi:10.1007/s12013-016-0743-8) contains supplementary material, which is available to authorized users. logy_of_HIV/AIDS#Sub-Saharan_Africa). Highly active antiretroviral therapy (HAART), approved by the Food and Drug Administration (FDA), is currently the most effective therapeutic regimen decreasing the viral load, which prolongs the survival of AID-free patients [4, 40, 58] . These drugs include inhibitors of three essential HIV enzymes protease, integrase (IN), reverse transcriptase (RT) as well as entry and progression enzyme inhibitors of in the HIV life cycle [11, 41] . In spite of all currently existing anti-HIV therapies, resistant strains still poses a problem. This creates a demanding necessity for the identification of new antivirals that do not induce drug resistance. Hsp90 is an ATP dependent chaperone protein crucial for maintaining active forms of other proteins which are known as client proteins [28] . Hsp90 has also been reported to be responsible for the replication of nearly all kinds of viruses such as DNA, RNA viruses, and double-stranded RNA viruses [23] . Hsp90 is a large protein, forming a homodimeric complex. Each monomer shares a common domain organization which consists of a C-terminal dimerization domain (CTD), a middle domain (MD) and an N-terminal ATPase domain (NTD) (Fig. 1) , which is found mostly in almost every part of eukaryotic cells [15, 35, 46] . Amongst other proteins, Hsp90 is one of the most profuse proteins in the cytoplasm, where it makes up approximately 1-2 % of the total protein levels (Anon n.d.). It is responsible for the controls of activities, such as the maturation, localization of a selected large number of substrates called ''client proteins'' in the cytoplasm. These clients proteins are part of several processes including transcription, translation, mitochondrial function, kinetochore assembly, inflammation, immunology, cellular antiviral defense pathways, apoptosis, centrosome function and cell cycle [18, 30, 39] . Studies have demonstrated that the replication of viruses to be sensitive to Hsp90 inhibitors at very low concentrations without affecting cellular viability, due to their dependence in Hsp90 for their survival. Hsp90 is regarded as an essential protein for malignant transformation and progression, thus target protein for anti-cancer drugs [45, 54] . A previous study in our laboratory gave a computational perspective of Hsp90 as an anti-cancer target [36] . Recently it has been reported as a potential target for anti-HIV therapy [60] . Studies have recently revealed that the CTD of Hsp90 protein has a crucial role in the dimerization of the chaperones with an additional binding site for inhibitors [52] . It has also been reported that the nucleotide binding site situated at the CTD is more favorable to guanosine-5 0triphosphate (GTP) over ATP, which would possibly enhance the selectivity of the new inhibitors against this site [9, 17, 59] . Herein we focus on the C-terminal dimerization inhibitors as this is the vital domain for the biological activity of Hsp90 [11] . To date, a broad-spectrum of antivirals have been identified, including Hsp90 inhibitors [5, 16, 23, 53] . Numerous antiviral activity of Hsp90 inhibitors have been proved against tissue culture of ill health conditions such as picornaviruses poliovirus, coxsackievirus, rhinovirus, influenza virus, paramyxoviruses, hepatitis C virus (HCV), Ebola virus, vesicular stomatitis virus, severe acquired respiratory syndrome (SARS), feline viral rhinotracheitis (FHV), HIV, vaccinia virus, and herpes viruses [14, 16, 32, 33, 44, 57, 60, 61, 64] . Gryase B inhibitors (novobiocin, coumermycin A1 (C-A1) and daunorubicin) have been reported to impair HIV-1 replication the most, by targeting Hsp90 [60] . It has been shown that they inhibit the formation Hsp90 dimer by binding in the C-terminal domain, which impairs viral gene expression [60] . The present study looked at C-A1 (Fig. 2) , which has previously been screened for its antiretroviral activity by means of a chemical genetics (CG) approach [60] . This study confirmed the impairment of viral gene expression by C-A1 binding to Hsp90. C-A1 weakens both the HIV-1 integration and gene expression possesses. Molecular dynamic simulations and post-dynamics (PD) analysis appeared as the closest methods to the actual experimental work and is of great help in the understanding of the complex biological phenomenas [12] . An application of continuous MD simulations to reveal folding and the unfolding activities of biological enzymes has made it convenient to fathom the dimerization process complexity of macromolecules [8, 20, 25, 37, 55] . The heightened utilities of the PD methods demonstrated to be effectual in comprehending the conformational background of biological macromolecules. The most extensively improved PD methodologies used to seek structural variations within distinctive biological phenomenas include principal component analysis PCA or essential dynamics (ED) [11] . Although PCA has attested its competency to be implement in the likening of motions of different macromolecules there are additional techniques, which were applied to understand such biological conformational behaviors [9, 49] . Such may include dynamic cross-correlation (DCC), which was valuable to gain insights on macromolecular motions in different biological systems [22] . To understand both the dimerization and the inhibition processes of Hsp90, virtual MD simulations were employed for the bound and the unbound conformations. To our understanding, this is the first report of such adequate computational insight on Hsp90 as a vital HIV target. Thus, we trust that this account aids as a significant computational study towards an improved understanding of the Hsp90 protein structure, dynamics and its mechanism of inhibition. Herein, we presented the structural and dynamic insight that can be executed in the discovery of drugs and the development of more effective HIV inhibitors against Hsp90. There is no available crystal structure of the human C-terminal portion of Hsp90, however structures of yeast (2CG9) [2] and E. coli (1SF8) [26] Hsp90 orthologues are available. 2CG9 was chosen for this study due to the fact there is very high sequence and structural conservation among Hsp90 proteins from prokaryotes to Homo sapien s [2] . The ligand and protein systems were prepared in the similar manner as mentioned in the previous reports [42, 43] . Auto-Dock Vina [56] was used to generate docked conformation of C-A1 inside Hsp90's active site using the same procedure previously reported [60] . The unbound and C-A1 bound models of Hsp90 proteins were exposed to all-atom, continuous MD simulations in precise solvent using the GPU version PMEMD engine integrated with Amber 14 [3, 6, 7] . The ff99sb force field [19] applied with Amber 14 was employed in the description of the protein systems. To neutralize the system, missing hydrogen atoms and counter ions were added using the Leap module, integrated with Amber 14. All the systems were submerged into an orthorhombic box with TIP3P [29] water molecules in a manner that none of the protein atoms were contained in 12 Å of any box edge prior to the set up. The particle mesh Ewald (PME) [27] technique was employed to calculate the continuous electrostatic interactions and van de Waals cut-off of 12 Å . Initial system minimization for 500 steps using steepest descent was followed by 100,000 steps using conjugate gradients algorithm was performed to relax the system and remove all unfavorable clashes. CPU version of Amber14 was used to perform all the minimization but before the minimization step, a gradual heating from 0 to 300 K by means of Langevin thermostat [62] with a collision frequency of 1/ps using an official ensemble (NVT) was used. All systems were successively equilibrated at 300 K in a NPT ensemble for 100,000 ps with no restrained and to uphold the pressure at 1 bar Berendsen barostat [10] was used. To restrict the bonds of all covalently bonded atoms to hydrogen atoms for all MD simulations runs the SHAKE algorithm [34] was used. The MD was run continuously in a NPT ensemble for 100 ns with a target pressure of 1 bar and a pressure-coupling constant of 2 ps. In every 1 ps the trajectories were examined and more analyses namely RMSD, RMSF, Rg, inter atomic distances, DCC and PCA were performed using PTRAJ and CPPTRAJ modules intergrated in Amber 14. All plots and visualization were performed in Chimera molecular modeling tool [47] and Origin data analysis software (www.originlab.com) respectively. The docked complex of Hsp90-C-A1 was calculated to predict the binding efficiencies and confirm the HCAAR C-A1 binding to Hsp90. The binding affinity predictions were executed using the MMPBSA with Eq. (1) [38] and MMGBSA method that with Eq. (2) [24] . Where, DE MM is the differentiation between the minimized energies of the Hsp90-C-A1 complex and the total energies of the Hsp90 and Hsp90 inhibitor including the electrostatic and the van der Waals energies. The change in entropy of the ligand binding conformations is presented by TDS, the variance in the P/GBSA solvation energies of the Hsp90-C-A1 complex and the summation of the solvation energies of the Hsp90 and Hsp90 inhibitor is represented by DG solv and DG SA is the difference in the surface area energies for the Hsp90 and Hsp90 inhibitor. Both MMPBSA and MMGBSA methods have been recognized to ensure the inhibitors are ranked accurately in accordance with their binding energy and thus can aid as a potent tool in the drug design research [1, 50, 63] . Principle component analysis (PCA) is extensively used to effortlessly examine and identify data generated from MD simulations to highpoint principal modes accountable for changes in conformation [11] . The PTRAJ and CPPTRAJ modules were used in performing PCA C-alpha atoms [48] integrated with Amber 14 [13] The plots analysis showing the central conformational motions descriptive of each structure were created using Origin data analyses programme (www.originlab.com). The DCCM between fluctuations based on the residues throughout the simulation was analyzed using the module called CPPTRAJ assimilated with Amber 14. The DCCM is best defined by the equation below [31] : Where r i (t) defines the vector of the atom's coordinates as a function of time t, \ Á [ t measures the time ensemble average and Dr i (t) = r i(t) -\ r i (t) [ t . The backbone Ca atomic fluctuations were considered during the DCCM analysis. The lig-plot (Anon n.d.) analysis is best suited for displaying the 2D interaction between the ligand and the amino acid residue contributing to the ligand binding to the receptor. The structure of inhibitor C-A1 docked within the Hsp90 dimer's binding site gave essential evidence on the binding of the ligand inside the active site (S1). It was observed that C-A1 formed hydrogen bonds interactions with Arg591.B, Ala595.A, Glu660.A, Lys480.B, Ser663.A and Ser666.A (S1) residues of both Hsp90 dimer subunits. Interestingly, the binding of C-A1 is at the dimer of Hsp90 possessing multiple common residues from both sub-unit A and B in its active site. The averaged values of the calculated free energy, using MM/PBSA approach; over the 100 ns MD simulations for each CAAR to the binding of C-A1 to Hsp90 are illustrated in Fig. 8 and S3 . The CAAR formed hydrogen bond interactions, which is one of the key interactions favorable and vital for binding of a ligand to its receptor, possess higher binding affinities. The plot of the interaction (Fig. 3) high spots the C-A1 location within the Hsp90 dimer and reciprocal residues' active site, together with the hydrogen bond interactions shown in red dotted lines, from each subunit. The existence of hydrophobic/aromatic interactions, which promotes the appropriate binding of inhibitor into the dimer, consisting of key active site residues involved in dimer packing, including Ser595.A and Glu660.A, are visibly illustrated in the binding mode. It is assumed that binding of C-A1 in the Hsp90 dimer's active site results to the dimer inhibition of Hsp90. The dynamic insight of the dimerization mechanism process of Hsp90 through inhibition by C-A1 was obtained from the continuous MD simulations using the structure of Hsp90-C-A1 complex (docked) (S2). Protein dynamics, structure, and function are highly interrelated. The conformational changes of Hsp90 are particularly critical for enzyme function [51, 65] . Studying internal atomistic motions can unravel the dynamic nature of Hsp90, which is responsible for numerous biological functions as enzymes. Molecular dynamic simulations were performed for each enzyme (Apo and bound) to establish the drug-binding pattern and to better understand the interaction of C-A1 with the CTD, to ensure the stability of the MD trajectories dynamic and the stabilities dissimilarities of MD simulations. The values of the RMSD for the protein backbone atoms and side chain atoms in relation to the initial minimized structure throughout the simulations were calculated for the Ca atoms for the dimer (S4) and whole protein (Fig. 4) . The average structure for the complete simulation was measured by calculating values for all residues of both in Apo and bound conformations (Fig. 5) by means of the RMSF in order to obtain the system's subset movement. Residues between 400 and 600 have the highest RMSF values, which corresponds to CTD. The most stable areas are the areas with the lowest RMSF values. Overall the bound enzyme has slightly lower RMSF values of 1.88 Å compared to the Apo with the value of 2.64 Å , which indicates that the backbone of the bound conformation is slightly more stable than that of the free enzyme. The mass-weighted root mean square distance of a collection of atoms from their common center of mass best describes the R g . To measure for the compactness of a structure, the R g was calculated. The analysis shown in Fig. 6 and S5 provides us with an understanding into the dimensions of the dimeric cleavage and of the overall Hsp90 protein respectively. We observed the major fluctuations in both systems between 0 and 100,000 ps. The R g indicated higher structural deviations for the apo conformation in comparison to the bound conformation. The mean value of the bound is lower than that of the apo, which indicates that the bound protein is stable. This could be due to C-A1 inhibitory effect on the Hsp90 dimerization. The trajectories' projection obtained at 300 K onto the first two principal components (PC1, PC2) presented the motion of two proteins in phase space. Clusters were better defined in the bound conformation than the apo conformation. Furthermore, the bound conformation covered a larger region of phase space particularly along PC1 and PC2 plane than the apo conformation as depicted in Fig. 7 . In an effort to offer more insight into the involvement of each amino acid residue toward the binding of C-A1, PRED was computed for the system. This helps us to observe the residues with the most contribution to the ligand binding, as seen from Fig. 8 the residues from both monomers showed interaction favorable to the ligand via different interactions. This further advocates that C-A1 inhibit dimerization. This information will ideally assist medicinal chemists in the design of inhibitors that precisely interrelate with these amino acid residues. The DCCM is shown in Fig. 9 and S7, which has been calculated for both bound and unbound conformations of Hsp90 dimer subunits and Hsp90 as a whole. The free state of Hsp90 displayed a more correlated motion, which further justifies the occurrence of the dimerization process resulting in a more correlated residue-residue interaction. The bound state, on the contrary, showed a much greater reduction in correlated motions during simulations. With these findings we can thus conclude that C-A1 inhibited the dimerization mechanism of Hsp90. The suitability of C-A1 to inhibit dimerization mechanism is hence confirmed as Hsp90 ground state is stabilized. The distance between the dimer residues in each dimeric subunit were calculated by measuring the interatomic distances between the Ca atoms along the simulations to gain an in-depth insight of dynamics of the dimerization process. Figure 10 and S3 highlights these residues. The distances between the orthogonal residues for the apo conformation were comparatively higher than that of the bound conformation as the average Ca distances specified on Table 1 and S6. The standard deviation of distances for the apo conformation was higher than that of the bound conformation, indicating that the latter enzyme is more stable (agreeing with the RMSF, R g and PCA studies). These substantial dissimilarities in distances between residues in the arranging of the dimer guaranteed the inhibitions of Hsp90 by C-A1 through the formation of a stable interaction with residues involved in the dimerization, preventing the dimer from opening and allowing the substrate from entering. This fact further correlates with previous biochemical analyses aimed at mapping the binding site of C-A1 in Hsp90 suggesting that the viral gene expression inhibited by C-A1 through the interference with Hsp90 dimerization [60] . Snapshots for the MD simulations for both the apo and the C-A1-bound conformations of Hsp90 were taken at specific time interval corresponding to the distances of orthogonal residues (S6) where it is evident for both conformations opening and closing. This further substantiates the lack of dimerization in the Hsp90-C-A1 complex (Fig. 11 ). Molecular dynamic simulations revealed difference events of dimer mechanism of Hsp90 in its unbound and C-A1bound conformations. Micro-molecule inhibitors such as C-A1 prevent the dimerization process of by targeting the dimer of Hsp90. We believe a more conformational rigid system was due a dimerization process being mired. This was indicated after the analysis of the RMSF and R g proving a more conformational supple nature of unbound Hsp90, whereas the C-A1-bound conformation of Hsp90 was found to be more stable with a declined dimerization throughout the MD simulation. The binding of C-A1 to the Hsp90 dimer subunits inhibits the process of dimerization. Observing the distance between the residues that are involved in the dimerization process further validates this. This is the first report emphasizing the important computational insight of the features of a crucial HIV target, which would also serve as an appropriate foundation in the process of designing new compounds against Hsp90 as future anti-HIV drugs. Fig. 10 The dimeric region residues from both subunits assumed to be part of the Hsp90 dimerization process Burger's medicinal chemistry and drug discovery Crystal structure of an Hsp90-nucleotide-p23/Sba1 closed chaperone complex Amber14 Identification of broad-spectrum antiviral compounds and assessment of the druggability of their target for efficacy against respiratory syncytial virus (RSV) The Global HIV/AIDS epidemic The henry J. kaiser family foundation Methods for molecular dynamics simulations of protein folding/unfolding in solution A perspective on enzyme catalysis Molecular dynamics with coupling to an external bath An integrated molecular dynamics, principal component analysis and residue interaction network approach reveals the impact of M184V mutation on HIV reverse transcriptase resistance to lamivudine Challenges in computational studies of enzyme structure, function and dynamics Amber 14 reference manual A functional heat shock protein 90 chaperone is essential for efficient flock house virus RNA polymerase synthesis in Drosophila cells Dimerization and N-terminal domain proximity underlie the function of the molecular chaperone heat shock protein 90 Antiviral activity and RNA polyinerase degradation following Hsp90 inhibition in a range of negative strand viruses Probing molecular mechanisms of the Hsp90 chaperone: biophysical modeling identifies key regulators of functional dynamics Targeting Hsp90 for the treatment of cancer A point-charge force field for molecular mechanics simulations of proteins based on condensed-phase quantum mechanical calculations Molecular dynamics simulations and drug discovery HIV-1 reverse transcriptase still remains a new drug target: structure, function, classical inhibitors, and new inhibitors with innovative mechanisms of actions Molecular dynamics simulation of human LOX-1 provides an explanation for the lack of OxLDL binding to the Trp150Ala mutant Broad action of Hsp90 as a host chaperone required for viral replication Improving docking results via reranking of ensembles of ligand poses in multiple X-ray protein conformations with mm-gbsa Molecular dynamics simulations of protein folding from the transition state The crystal structure of the carboxy-terminal dimerization domain of htpG, the Escherichia coli Hsp90, reveals a potential substrate binding site An implementation of the smooth particle-mesh Ewald (PME) method on GPU hardware Advances in the clinical development of heat shock protein 90 (Hsp90) inhibitors in cancers Comparison of simple potential functions for simulating liquid water Regulation of tumor cell mitochondrial Homeostasis by an organelle-specific Hsp90 chaperone network A novel approach of dynamic cross correlation analysis on molecular dynamics simulations and its application to Ets1 dimer-DNA complex Destabilization of PDK1 by Hsp90 inactivation suppresses hepatitis C virus replication through inhibition of PRK2-mediated viral RNA polymerase phosphorylation Heat shock protein and heat shock factor 1 expression and localization in vaccinia virus infected human monocyte derived macrophages A fast SHAKE algorithm to solve distance constraint equations for small molecules in molecular dynamics simulations Conformational dynamics of the molecular chaperone Hsp90 Heat-shock protein 90 (Hsp90) as anticancer target for drug discovery: an ample computational perspective How fast-folding proteins fold Accurate prediction of the relative potencies of members of a series of kinase inhibitors using molecular docking and MM-GBSA scoring Diverse cellular functions of the Hsp90 molecular chaperone uncovered using systems approaches Insights into the reasons for discontinuation of the first highly active antiretroviral therapy (HAART) regimen in a cohort of antiretroviral naive patients Dynamic features of apo and bound HIV-Nef protein reveal the anti-HIV dimerization inhibition mechanism Dual acting HIV inhibitors: integrated rational in silico design strategy Computer-aided perspective for the design of flexible HIV non-nucleoside reverse transcriptase inhibitors (NNRTIs): de-novo design, virtual screening and molecular dynamics simulations Hsp90 inhibitors suppress HCV replication in replicon cells and humanized liver mice Hsp90 molecular chaperone inhibitors: are we there yet? Structure and mechanism of the Hsp90 molecular chaperone machinery UCSF Chimera-a visualization system for exploratory research and analysis PTRAJ and CPPTRAJ: software for processing and analysis of molecular dynamics trajectory data Enzymatic effect on reactant and transition states. Chalcone isomerase Accurate ensemble molecular dynamics binding free energy ranking of multidrug-resistant HIV-1 proteases Activation of Hsp90 enzymatic activity and conformational dynamics through rationally designed allosteric ligands. Chemistry (Weinheim an der Bergstrasse, Germany) Structural models and binding site prediction of the C-terminal domain of human Hsp90: a new target for anticancer drugs HSP90 inhibitors: current development and potential in cancer therapy. Recent patents on anti-cancer drug discovery Hsp90 inhibitors as anti-cancer agents, from basic discoveries to clinical development Describing protein folding kinetics by molecular dynamics simulations. 1. Theory . The AutoDock Vina: improving the speed and accuracy of docking with a new scoring function, efficient optimization, and multithreading 17-AAG, an Hsp90 inhibitor, suppress hepatitis C virus (HCV) replication Highly active antiretroviral therapy Structural and computational biology of the molecular chaperone Hsp90: from understanding molecular mechanisms to computer-based inhibitor design. Current topics in medicinal chemistry Gyrase B inhibitor impairs HIV-1 replication by targeting Hsp90 and the capsid protein A targeted analysis of cellular chaperones reveals contrasting roles for heat shock protein 70 in flock house virus RNA replication Self-guided Langevin dynamics simulation method Using molecular modelling to study interactions between molecules with biological activity Heat-shock protein 90 promotes nuclear transport of herpes simplex virus 1 capsid protein by interacting with acetylated tubulin Hsp90 and co-chaperones twist the functions of diverse client proteins Acknowledgments The authors would like to acknowledge the School of Health Science, University of KwaZulu-Natal for financial support. Computational resources granted by the Center for High Performance Computing (http://www.chpc.ac.za) and San Diego is highly appreciated. Conflict of interests The authors avow no conflict of interest to the above article.