Heterogeneous behavior of metalloproteins toward metal ion binding and selectivity: insights from molecular dynamics studies

About one-third of the existing proteins require metal ions as cofactors for their catalytic activities and structural complexities. While many of them bind only to a specific metal, others bind to multiple (different) metal ions. However, the exact mechanism of their metal preference has not been deduced to clarity. In this study, we used molecular dynamics (MD) simulations to investigate whether a cognate metal (bound to the structure) can be replaced with other similar metal ions. We have chosen seven different proteins (phospholipase A₂, sucrose phosphatase, pyrazinamidase, cysteine dioxygenase (CDO), plastocyanin, monoclonal anti-CD4 antibody Q425, and synaptotagmin 1 C₂B domain) bound to seven different divalent metal ions (Ca²⁺, Mg²⁺, Zn²⁺, Fe²⁺, Cu²⁺, Ba²⁺, and Sr²⁺, respectively). In total, 49 MD simulations each of 50 ns were performed and each trajectory was analyzed independently. Results demonstrate that in some cases, cognate metal ions can be exchanged with similar metal ions. On the contrary, some proteins show binding affinity specifically to their cognate metal ions. Surprisingly, two proteins CDO and plastocyanin which are known to bind Fe²⁺ and Cu²⁺, respectively, do not exhibit binding affinity to any metal ion. Furthermore, the study reveals that in some cases, the active site topology remains rigid even without cognate metals, whereas, some require them for their active site stability. Thus, it will be interesting to experimentally verify the accuracy of these observations obtained computationally. Moreover, the study can help in designing novel active sites for proteins to sequester metal ions particularly of toxic nature.     

Targeted molecular dynamics simulation studies of binding and conformational changes in E. coli MurD

Enzymes involved in the biosynthesis of bacterial peptidoglycan, an essential cell wall polymer unique to prokaryotic cells, represent a highly interesting target for antibacterial drug design. Structural studies of E. coli MurD, a three-domain ATP hydrolysis driven muramyl ligase revealed two inactive open conformations of the enzyme with a distinct C-terminal domain position. It was hypothesized that the rigid body rotation of this domain brings the enzyme to its closed active conformation, a structure, which was also determined experimentally. Targeted molecular dynamics 1 ns-length simulations were performed in order to examine the substrate binding process and gain insight into structural changes in the enzyme that occur during the conformational transitions into the active conformation. The key interactions essential for the conformational transitions and substrate binding were identified. The results of such studies provide an important step toward more powerful exploitation of experimental protein structures in structure-based inhibitor design.     

Molecular insights of benzodipyrazole as CDK2 inhibitors: combined molecular docking,molecular dynamics,and 3D QSAR studies

Abstract

Benzodipyrazoles have been previously evaluated for their in vitro CDK2 inhibitory activity. In the current investigation, we identified a six-feature common pharmacophore model (AADDRR.33) which is predicted to be responsible for CDK2 inhibition. An efficient 3D QSAR (r^2?=?0.98 and q^2?=?0.82) model was also constructed by employing PLS regression analysis. From the molecular docking studies, we examined the binding patterns of compound 7aa with the target protein and also calculated the binding energy using MM-GBSA calculations. Three hydrogen bonds with Lys 33, Glu 81, and Leu 83 are conserved even after 1000?ps run in a molecular dynamics simulation. We identified the slight displacement in bond lengths and the conformational changes occurred during the dynamics. The results also elucidated the protein residue–ligand interaction fractions which clearly explained the involvement of non-H-bond interactions.     

Translocation mechanism of P-glycoprotein and conformational changes occurring at drug-binding site: Insights from multi-targeted molecular dynamics

P-glycoprotein (P-gp) is well known for multidrug resistance in drug therapy. Its over-expression results into the increased efflux of therapeutic agents rendering them inefficacious. A clear understanding of P-gp efflux mechanism and substrate/inhibitor interactions during the course of efflux cycle will be crucial for designing effective P-gp inhibitors, and therapeutic agents that are non-substrate to P-gp. In the present work, we have modeled P-gp in three different catalytic states. These models were utilized for elucidation of P-gp translocation mechanism using multi-targeted molecular dynamics (MTMD). The gradual changes occurring in P-gp structure from inward open to outward open conformation were sampled out. A detailed investigation of conformational changes occurring in trans-membrane domains (TMDs) during the course of catalytic cycle was carried out. Movements of each TM helices in response to pronounced twisting and translatory motion of NBDs were measured quantitatively. The role of intracellular coupling helices (ICHs) during the structural transition of P-gp was studied, and observed as vital links for structural transition. A close observation of displacements and conformational changes in the residues lining drug-binding pocket was also carried out. Further, we have analyzed the molecular interactions of P-gp substrates/inhibitors during the P-gp translocation to find out how stable binding interactions of a compound at drug-binding site(s) in open conformation, becomes highly destabilized in closed conformation. The study revealed striking differences between the molecular interactions of substrate and inhibitor; inhibitors showed a tendency to maintain stable binding interactions during the catalytic transition cycle.     

Structural study of the zinc and cadmium complexes of a type 2 plant (Quercus suber) metallothionein: insights by vibrational spectroscopy

Zn- and Cd-complexes of Quercus suber metallothionein (QsMT) were obtained by in vivo-synthesis, in order to obtain physiologically representative aggregates, and characterized by spectrometric and spectroscopic methods. The secondary structure elements and the coordination environments of the metal binding sites of the two aggregates were determined, as well as the main metal-containing species formed. The results obtained from the analysis of the Raman and IR spectra reveal that these metal-MT complexes predominantly contain beta-sheet elements (about 60%), whereas they lack alpha-helices. These structural features slightly depend on the divalent metal bound. In particular, Cd(II) binding to QsMT induces a slight increase of the beta-sheet percentage, as well as a decrease in beta-turn elements with respect to Zn(II) binding. Conversely, the in vivo capability of QsMT to inglobe metal and sulfide ions is metal-depending. Spectroscopic vibrational data also confirm the presence of sulfide ligands in the metal clusters of both Zn- and Cd-QsMT, while the participation of the spacer His residue in metal coordination was only found in Cd-QsMT, in agreement with the CD results. Overall data suggest different coordination environments for Zn(II) and Cd(II) ions in QsMT.     

Folding and stability of the three-stranded beta-sheet peptide Betanova: insights from molecular dynamics simulations

The dynamics of the three-stranded beta-sheet peptide Betanova has been studied at four different temperatures (280, 300, 350, and 450 K by molecular dynamics simulation techniques, in explicit water. Two 20-ns simulations at 280 K indicate that the peptide remains very flexible under "folding" conditions sampling a range of conformations that together satisfy the nuclear magnetic resonance (NMR)-derived experimental constraints. Two simulations at 300 K (above the experimental folding temperature) of 20 ns each show partial formation of "native"-like structure, which also satisfies most of the NOE constraints at 280 K. At higher temperature, the presence of compact states, in which a series of hydrophobic contacts remain present, are observed. This is consistent with experimental observations regarding the role of hydrophobic contacts in determining the peptide's stability and in initiating the formation of turns and loops. A set of different structures is shown to satisfy NMR-derived distance restraints and a possible mechanism for the folding of the peptide into the NMR-determined structure is proposed.     

Docking and molecular dynamics studies of new potential inhibitors of the human epidermal receptor 2

Compounds similar to lapatinib and gefitinib have been investigated as potential inhibitors of the intracellular receptor tyrosine kinase (RTK) domain of the human epidermal receptor 2 (HER2), which is a promising molecular target to the drug design of new chemotherapies for breast, lung, ovarian and colorectal cancers. In this study, we have searched potential HER2 inhibitors used for treatment of other illnesses such as hepatitis, bacterial infections and sexual impotence screened in the DrugBank. The compounds selected were subjected to virtual screening docking in order to evaluate the main interactions between them and the RTK domain of HER2. The selected compounds were investigated by flexible docking, molecular dynamics studies and ΔG _bind calculations. The results suggest that antrafenine, saprisartan, reserpine, irinotecan and udenafil are potential candidates to inhibit the RTK domain of HER2.     

Probing ligand binding modes of human cytochrome P450 2J2 by homology modeling, molecular dynamics simulation, and flexible molecular docking

Cytochrome P450 (P450) 2J2 catalyzes epoxidation of arachidonic acid to eicosatrienoic acids, which are related to a variety of diseases such as coronary artery disease, hypertension, and carcinogenesis. Recent experimental data also suggest that P450 2J2 could be a novel biomarker and a potential target for cancer therapy. However, the active site topology and substrate specificity of this enzyme remain unclear. In this study, a three-dimensional model of human P450 2J2 was first constructed on the basis of the crystal structure of human P450 2C9 in complex with a substrate using homology modeling method, and refined by molecular dynamics simulation. Flexible docking approaches were then employed to dock four ligands into the active site of P450 2J2 in order to probe the ligand-binding modes. By analyzing the results, active site architecture and certain key residues responsible for substrate specificity were identified on the enzyme, which might be very helpful for understanding the enzyme's biological role and providing insights for designing novel inhibitors of P450 2J2.     

Capecitabine as a minor groove binder of DNA: molecular docking,molecular dynamics,and multi-spectroscopic studies

The interaction mechanism and binding mode of capecitabine with ctDNA was extensively investigated using docking and molecular dynamics simulations, fluorescence and circular dichroism (CD) spectroscopy, DNA thermal denaturation studies, and viscosity measurements. The possible binding mode and acting forces on the combination between capecitabine and DNA had been predicted through molecular simulation. Results indicated that capecitabine could relatively locate stably in the G-C base-pairs-rich DNA minor groove by hydrogen bond and several weaker nonbonding forces. Fluorescence spectroscopy and fluorescence lifetime measurements confirmed that the quenching was static caused by ground state complex formation. This phenomenon indicated the formation of a complex between capecitabine and ctDNA. Fluorescence data showed that the binding constants of the complex were approximately 2 × 10⁴ M^?1. Calculated thermodynamic parameters suggested that hydrogen bond was the main force during binding, which were consistent with theoretical results. Moreover, CD spectroscopy, DNA melting studies, and viscosity measurements corroborated a groove binding mode of capecitabine with ctDNA. This binding had no effect on B-DNA conformation.     

Effect of molecular confinement on internal enzyme dynamics: frequency domain fluorometry and molecular dynamics simulation studies

The tryptophanyl emission decay of the mesophilic beta-galactosidase from Aspergillus oryzae free in buffer and entrapped in agarose gel is investigated as a function of temperature and compared to that of the hyperthermophilic enzyme from Sulfolobus solfataricus. Both enzymes are tetrameric proteins with a large number of tryptophanyl residues, so the fluorescence emission can provide information on the conformational dynamics of the overall protein structure rather than that of the local environment. The tryptophanyl emission decays are best fitted by bimodal Lorentzian distributions. The long-lived component is ascribed to close, deeply buried tryptophanyl residues with reduced mobility; the short-lived one arises from tryptophanyl residues located in more flexible external regions of each subunit, some of which are involved in forming the catalytic site. The center of both lifetime distribution components at each temperature increases when going from the free in solution mesophilic enzyme to the gel-entrapped and hyperthermophilic enzyme, thus indicating that confinement of the mesophilic enzyme in the agarose gel limits the freedom of the polypeptide chain. A more complex dependence is observed for the distribution widths. Computer modeling techniques are used to recognize that the catalytic sites are similar for the mesophilic and hyperthermophilic beta-galactosidases. The effect due to gel entrapment is considered in dynamic simulations by imposing harmonic restraints to solvent-exposed atoms of the protein with the exclusion of those around the active site. The temperature dependence of the tryptophanyl fluorescence emission decay and the dynamic simulation confirm that more rigid structures, as in the case of the immobilized and/or hyperthermophilic enzyme, require higher temperatures to achieve the requisite conformational dynamics for an effective catalytic action and strongly suggest a link between conformational rigidity and enhanced thermal stability.     

Targeted molecular dynamics simulation studies of calcium binding and conformational change in the C-terminal half of gelsolin

Gelsolin consists of six related domains (G1-G6) and the C-terminal half (G4-G6) acts as a calcium sensor during the activation of the whole molecule, a process that involves large domain movements. In this study, we used targeted molecular dynamics simulations to elucidate the conformational transitions of G4-G6 at an atomic level. Domains G4 and G6 are initially ruptured, followed by a rotation of G6 by approximately 90 degrees , which is the dominant conformational change. During this period, local conformational changes occur at the G4 and G5 calcium-binding sites, facilitating large changes in interdomain distances. Alterations in the binding affinities of the calcium ions in these three domains appear to be related to local conformational changes at their binding sites. Analysis of the relative stabilities of the G4-G6-bound calcium ions suggests that they bind first to G6, then to G4, and finally to G5.     

Potential charge transfer probe induced conformational changes of model plasma protein human serum albumin: Spectroscopic, molecular docking, and molecular dynamics simulation study

The nature of binding of specially designed charge transfer (CT) fluorophore at the hydrophobic protein interior of human serum albumin (HSA) has been explored by massive blue-shift (82 nm) of the polarity sensitive probe emission accompanying increase in emission intensity, fluorescence anisotropy, red edge excitation shift, and average fluorescence lifetimes. Thermal unfolding of the intramolecular CT probe bound HSA produces almost opposite spectral changes. The spectral responses of the molecule reveal that it can be used as an extrinsic fluorescent reporter for similar biological systems. Circular dichrosim spectra, molecular docking, and molecular dynamics simulation studies scrutinize this binding process and stability of the protein probe complex more closely.     

Synthesis,molecular docking,binding free energy calculation and molecular dynamics simulation studies of benzothiazol-2-ylcarbamodithioates as Staphylococcus aureus MurD inhibitors

Abstract

A new series of benzothiazol-2-ylcarbamodithioate functional compounds 5a-f has been designed, synthesized and characterized by spectral data. These compounds were screened for their in vitro antibacterial activity against strains of Staphylococcus aureus (NCIM 5021, NCIM 5022 and methicillin-resistant isolate 43300), Bacillus subtilis (NCIM 2545), Escherichia coli (NCIM 2567), Klebsiella pneumoniae (NCIM 2706) and Psudomonas aeruginosa (NCIM 2036). Compounds 5a and 5d exhibited significant activity against all the tested bacterial strains. Specifically, compounds 5a and 5d showed potent activity against K. pneumoniae (NCIM 2706), while compound 5a also displayed potent activity against S. aureus (NCIM 5021). Compound 5d showed minimum IC₅₀ value of 13.37?μM against S. aureus MurD enzyme. Further, the binding interactions of compounds 5a-f in the catalytic pocket have been investigated using the extra-precision molecular docking and binding free energy calculation by MM-GBSA approach. A 30?ns molecular dynamics simulation of 5d/modeled S. aureus MurD enzyme was performed to determine the stability of the predicted binding conformation.     

Impact of deglycosylation and thermal stress on conformational stability of a full length murine igG2a monoclonal antibody: Observations from molecular dynamics simulations

With the rise of antibody based therapeutics as successful medicines, there is an emerging need to understand the fundamental antibody conformational dynamics and its implications towards stability of these medicines. Both deglycosylation and thermal stress have been shown to cause conformational destabilization and aggregation in monoclonal antibodies. Here, we study instabilities caused by deglycosylation and by elevated temperature (400 K) by performing molecular dynamic simulations on a full length murine IgG2a mAb whose crystal structure is available in the Protein Data bank. C^α‐atom root mean square deviation and backbone root mean square fluctuation calculations show that deglycosylation perturbs quaternary and tertiary structures in the C_H2 domains. In contrast, thermal stress pervades throughout the antibody structure and both Fabs and Fc regions are destabilized. The thermal stress applied in this study was not sufficient to cause large scale unfolding within the simulation time and most amino acid residues showed similar average solvent accessible surface area and secondary structural conformations in all trajectories. C_H3 domains were the most successful at resisting the conformational destabilization. The simulations helped identify aggregation prone regions, which may initiate cross‐β motif formation upon deglycosylation and upon applying thermal stress. Deglycosylation leads to increased backbone fluctuations and solvent exposure of a highly conserved APR located in the edge β‐strand A of the C_H2 domains. Aggregation upon thermal stress is most likely initiated by two APRs that overlap with the complementarity determining regions. This study has important implications for rational design of antibody based therapeutics that are resistant towards aggregation. Proteins 2013. © 2012 Wiley Periodicals, Inc.     

Kinetics, in silico docking, molecular dynamics, and MM-GBSA binding studies on prototype indirubins, KT5720, and staurosporine as phosphorylase kinase ATP-binding site inhibitors: the role of water molecules examined

With an aim toward glycogenolysis control in Type 2 diabetes, we have investigated via kinetic experiments and computation the potential of indirubin (IC?? > 50 μM), indirubin-3'-oxime (IC?? = 144 nM), KT5720 (K(i) = 18.4 nM) and staurosporine (K(i) = 0.37 nM) as phosphorylase kinase (PhKγtrnc) ATP-binding site inhibitors, with the latter two revealed as potent inhibitors in the low nM range. Because of lack of structural information, we have exploited information from homologous kinase complexes to direct in silico calculations (docking, molecular dynamics, and MMGBSA) to predict the binding characteristics of the four ligands. All inhibitors are predicted to bind in the same active site area as the ATP adenine ring, with binding dominated by hinge region hydrogen bonds to Asp104:O and Met106:O (all four ligands) and also Met106:NH (for the indirubins). The PhKγtrnc-staurosporine complex has the greatest number of receptor-ligand hydrogen bonds, while for the indirubin-3'-oxime and KT5720 complexes there is an important network of interchanging water molecules bridging inhibitor-enzyme contacts. The MM-GBSA results revealed the source of staurosporine's low nM potency to be favorable electrostatic interactions, while KT5720 has strong van der Waals contributions. KT5720 interacts with the greatest number of protein residues either by direct or 1-water bridged hydrogen bond interactions, and the potential for more selective PhK inhibition based on a KT5720 analogue has been established. Including receptor flexibility in Schr?dinger induced-fit docking calculations in most cases correctly predicted the binding modes as compared with the molecular dynamics structures; the algorithm was less effective when there were key structural waters bridging receptor-ligand contacts.     

Molecular dynamics investigations of structural and functional changes in Bcl-2 induced by the novel antagonist BDA-366

Apoptosis is a fundamental biological phenomenon, in which anti- or proapoptotic proteins of the Bcl-2 family regulate a committed step. Overexpression of Bcl-2, the prototypical antiapoptotic protein in this family, is associated with therapy resistance in various human cancers. Accordingly, Bcl-2 inhibitors intended for cancer therapy have been developed, typically against the BH3 domain. Recent experimental evidences have shown that the antiapoptotic function of Bcl-2 is not immutable, and that BDA-366, a novel antagonist of the BH4 domain, converts Bcl-2 from a survival molecule to an inducer of cell death. In this study, the underlying mechanisms of this functional conversion were investigated by accelerated molecular dynamics simulation. Results revealed that Pro127 and Trp30 in the BH4 domain rotate to stabilize BDA-366 via π-π interactions, and trigger a series of significant conformational changes of the α3 helix. This rearrangement blocks the hydrophobic binding site (HBS) in the BH3 domain and further prevents binding of BH3-only proteins, which consequently allows the BH3-only proteins to activate the proapoptotic proteins. Analysis of binding free energy confirmed that BDA-366 cross-inhibits BH3-only proteins, implying negative cooperative effects across separate binding sites. The newly identified blocked conformation of the HBS along with the open to closed transition pathway revealed by this study advances the understanding of the Bcl-2 transition from antiapoptotic to proapoptotic function, and yielded new structural insights for novel drug design against the BH4 domain.

Communicated by Ramaswamy H. Sarma     

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.