MOLECULAR MODELING OF THE INTERACTIONS BETWEEN HISTONE DEACETYLASE 8 AND INHIBITORS

2012 ◽  
Vol 11 (04) ◽  
pp. 907-924 ◽  
Author(s):  
DAWEI HUANG ◽  
XIAOHUI LI ◽  
ZHILONG XIU

Inhibitors of histone deacetylases (HDACs) have become an attractive class of anticancer agent. To understand the interaction between HDAC8 and inhibitors, including "pan-" inhibitors that inhibit many HDACs isoforms and selective inhibitors with no linker domain, docking and molecular dynamics simulation were conducted. Docking results showed the presence of π-π interactions between "linkerless" inhibitors and the aromatic amino acid residues of HDAC8 in the active site. Binding between HDAC8 and inhibitors was also stabilized by hydrogen bond and hydrophobic interaction. In molecular dynamics simulations, the zinc ion was shown to coordinate one more atom of HDAC8-"linkerless" inhibitor complexes than HDAC8-"pan-" inhibitor complexes. Persistent hydrogen bonds also existed between Tyr306 of HDAC8 and some inhibitors. When inhibitors with large cap groups bound to the active pocket of HDAC8, Phe152 and Met274 shifted from their initial positions and the entrance of the active pocket became more open, resulting in the formation of sub-pocket. Hydrophobic interactions contributed most favorably to the binding free energy between HDAC8 and inhibitors. Lys33, Asp178, Asp267, Tyr306 and Leu308 of HDAC8 were favorable for binding with all inhibitors.

RSC Advances ◽  
2016 ◽  
Vol 6 (63) ◽  
pp. 58573-58585 ◽  
Author(s):  
Jianzhong Chen

Molecular dynamics simulation and binding free energy calculations were used to probe drug resistance of HIV-1 protease mutations toward inhibitors.


2021 ◽  
Vol 8 ◽  
Author(s):  
Xiaoqian Zhang ◽  
Hua Yu ◽  
Xiangdong Liu ◽  
Chen Song

The calcium release-activated calcium channel, composed of the Orai channel and the STIM protein, plays a crucial role in maintaining the Ca2+ concentration in cells. Previous studies showed that the L138F mutation in the human Orai1 creates a constitutively open channel independent of STIM, causing severe myopathy, but how the L138F mutation activates Orai1 is still unclear. Here, based on the crystal structure of Drosophila melanogaster Orai (dOrai), molecular dynamics simulations for the wild-type (WT) and the L210F (corresponding to L138F in the human Orai1) mutant were conducted to investigate their structural and dynamical properties. The results showed that the L210F dOrai mutant tends to have a more hydrated hydrophobic region (V174 to F171), as well as more dilated basic region (K163 to R155) and selectivity filter (E178). Sodium ions were located deeper in the mutant than in the wild-type. Further analysis revealed two local but essential conformational changes that may be the key to the activation. A rotation of F210, a previously unobserved feature, was found to result in the opening of the K163 gate through hydrophobic interactions. At the same time, a counter-clockwise rotation of F171 occurred more frequently in the mutant, resulting in a wider hydrophobic gate with more hydration. Ultimately, the opening of the two gates may facilitate the opening of the Orai channel independent of STIM.


2021 ◽  
Author(s):  
Xiaoqian Zhang ◽  
Hua Yu ◽  
Xiangdong Liu ◽  
Chen Song

The calcium release-activated calcium (CRAC) channel, composed of the Orai channel and the STIM protein, plays a crucial role in maintaining the Ca2+ concentration in cells. Previous studies showed that the L138F mutation in the human Orai1 creates a constitutively open channel independent of STIM, causing severe myopathy, but how the L138F mutation activates Orai1 is still unclear. Here, based on the crystal structure of Drosophila melanogaster Orai (dOrai), molecular dynamics simulations for the wild-type (WT) and the L210F (corresponding to L138F in the human Orai1) mutant were conducted to investigate their structural and dynamical properties. The results showed that the L210F dOrai mutant tends to have a more hydrated hydrophobic region (V174 to F171), as well as more dilated basic region (K163 to R155) and selectivity filter (E178). Sodium ions were located deeper in the mutant than in the WT. Further analysis revealed two local but essential conformational changes that may be the key to the activation. A rotation of F210, a previously undescribed feature, was found to result in the opening of the K163 gate through hydrophobic interactions. At the same time, a counter-clockwise rotation of F171 occurred more frequently in the mutant, resulting in a wider hydrophobic gate with more hydration. Ultimately, the opening of the two gates may facilitate the opening of the Orai channel independent of STIM.


2020 ◽  
Author(s):  
Dr. Chirag N. Patel ◽  
Dr. Prasanth Kumar S. ◽  
Dr. Himanshu A. Pandya ◽  
Dr. Rakesh M. Rawal

<p>The pandemic outbreak of COVID-19 virus (SARS-CoV-2) has become critical global health issue. The biophysical and structural evidence shows that SARS-CoV-2 spike protein possesses higher binding affinity towards angiotensin-converting enzyme 2 (ACE2) and hemagglutinin-acetylesterase (HE) glycoprotein receptor. Hence, it was selected as a target to generate the potential candidates for the inhibition of HE glycoprotein. The present study focuses on extensive computational approaches which contains molecular docking, ADMET prediction followed by molecular dynamics simulations and free energy calculations. Furthermore, virtual screening of NPACT compounds identified 3,4,5-Trihydroxy-1,8-bis[(2R,3R)-3,5,7-trihydroxy-3,4-dihydro-2H-chromen-2-yl]benzo[7]annulen-6-one, Silymarin, Withanolide D, Spirosolane and Oridonin were interact with high affinity. The ADMET prediction revealed pharmacokinetics and drug-likeness properties of top-ranked compounds. Molecular dynamics simulations and binding free energy calculations affirmed that these five NPACT compounds were robust HE inhibitor.</p>


2020 ◽  
Vol 76 (10) ◽  
pp. 993-1000
Author(s):  
Akihito Ikeda ◽  
Seia Urata ◽  
Tadashi Ando ◽  
Yasuhiro Suzuki ◽  
Yasufumi Sato ◽  
...  

Vasohibins regulate angiogenesis, tumor growth, metastasis and neuronal differentiation. They form a complex with small vasohibin-binding protein (SVBP) and show tubulin tyrosine carboxypeptidase activity. Recent crystal structure determinations of vasohibin–SVBP complexes have provided a molecular basis for complex formation, substrate binding and catalytic activity. However, the regulatory mechanism and dynamics of the complex remain elusive. Here, the crystal structure of the VASH1–SVBP complex and a molecular-dynamics simulation study are reported. The overall structure of the complex was similar to previously reported structures. Importantly, however, the structure revealed a domain-swapped heterotetramer that was formed between twofold symmetry-related molecules. This heterotetramerization was stabilized by the mutual exchange of ten conserved N-terminal residues from the VASH1 structural core, which was intramolecular in other structures. Interestingly, a comparison of this region with previously reported structures revealed that the patterns of hydrogen bonding and hydrophobic interactions vary. In the molecular-dynamics simulations, differences were found between the heterotetramer and heterodimer, where the fluctuation of the N-terminal region in the heterotetramer was suppressed. Thus, heterotetramer formation and flexibility of the N-terminal region may be important for enzyme activity and regulation.


2016 ◽  
Vol 15 (06) ◽  
pp. 1650054 ◽  
Author(s):  
Seifollah Jalili ◽  
Mina Maddah ◽  
Jeremy Schofield

Cisplatin and oxaliplatin are two widely-used anti-cancer drugs which covalently bind to a same location in DNA strands. Platinum agents make intrastrand and interstrand cross-links with the N7 atoms of guanine nucleotides which prevent DNA from polymerization by causing a distortion in the double helix. Molecular dynamics simulations and free energy calculations were carried out to investigate the binding of two platinum-based anti-cancer drugs with DNA. We compared the binding of these drugs which differ in their carrier ligands, and hence their potential interactions with DNA. When a platinum agent binds to nucleotides, it causes a high amount of deformation in DNA structure. To find the extent of deformation, torsion angles and base pair and groove parameters of DNA were considered. These parameters were compared with normal B-DNA which was considered as the undamaged DNA. The formation of hydrogen bonds between drugs and DNA nucleotides was examined in solution. It was shown that oxaliplatin forms more hydrogen bonds than cisplatin. Our results confirm that the structure of the platinated DNA rearranges significantly and cisplatin tries to deform DNA more than oxaliplatin. The binding free energies were also investigated to understand the affinities, types and the contributions of interactions between drugs and DNA. It was concluded that oxaliplatin tendency for binding to DNA is more than cisplatin in solvent environment. The binding free energy was calculated based on the MM/PBSA and MM/GBSA methods and the results of QM/MM calculations verified them.


Biology ◽  
2020 ◽  
Vol 10 (1) ◽  
pp. 2
Author(s):  
Kamrul Hasan Chowdhury ◽  
Md. Riad Chowdhury ◽  
Shafi Mahmud ◽  
Abu Montakim Tareq ◽  
Nujhat Binte Hanif ◽  
...  

Novel coronavirus disease (COVID-19) was identified from China in December 2019 and spread rapidly through human-to-human transmission, affecting so many people worldwide. Until now, there has been no specific treatment against the disease and repurposing of the drug. Our investigation aimed to screen potential inhibitors against coronavirus for the repurposing of drugs. Our study analyzed sequence comparison among SARS-CoV, SARS-CoV-2, and MERS-CoV to determine the identity matrix using discovery studio. SARS-CoV-2 Mpro was targeted to generate an E-pharmacophore hypothesis to screen drugs from the DrugBank database having similar features. Promising drugs were used for docking-based virtual screening at several precisions. Best hits from virtual screening were subjected to MM/GBSA analysis to evaluate binding free energy, followed by the analysis of binding interactions. Furthermore, the molecular dynamics simulation approaches were carried out to assess the docked complex’s conformational stability. A total of 33 drug classes were found from virtual screening based on their docking scores. Among them, seven potential drugs with several anticancer, antibiotic, and immunometabolic categories were screened and showed promising MM/GBSA scores. During interaction analysis, these drugs exhibited different types of hydrogen and hydrophobic interactions with amino acid residue. Besides, 17 experimental drugs selected from virtual screening might be crucial for drug discovery against COVID-19. The RMSD, RMSF, SASA, Rg, and MM/PBSA descriptors from molecular dynamics simulation confirmed the complex’s firm nature. Seven promising drugs for repurposing against SARS-CoV-2 main protease (Mpro), namely sapanisertib, ornidazole, napabucasin, lenalidomide, daniquidone, indoximod, and salicylamide, could be vital for the treatment of COVID-19. However, extensive in vivo and in vitro studies are required to evaluate the mentioned drug’s activity.


2021 ◽  
Author(s):  
Ahmed Ali Al-Karmalawy ◽  
Radwan Alnajjar ◽  
Mohammed Dahab ◽  
Ahmed Metwaly ◽  
Ibrahim Eissa

Background: Drug repurposing is the fastest effective method to provide treatment for coronavirus disease (COVID-19). Drugs that targeting a closely related virus with similar genetic material such as hepatitis C virus (HCV) and more specifically targeting a similar viral protease would be an excellent choice. Methods: In this study, we carried out a virtual screening for fifteen anti HCV drugs against COVID-19 main protease using computational molecular docking techniques. Moreover, Velpatasvir (4) and Sofosbuvir (13) drugs were further evaluated through molecular dynamics simulations followed by calculating the binding free energy using the molecular mechanics generalised born/solvent accessibility (MM-GBSA) approach. Results: The binding affinity descending order was N3 natural inhibitor (1), Velpatasvir (4), Sofosbuvir (13), Ombitasvir (3), Glecaprevir (2), Asunaprevir (8), Paritaprevir (10), Grazoprevir (11), Elbasvir (6), Ledipasvir (5), Daclatasvir (7), Pibrentasvir (9), Simeprevir (12), Dasabuvir (14), Taribavirin (16) and finally Ribavirin (15). Molecular dynamics simulation reveals that Sofosbuvir (13) has exciting properties and it was stable within the active site; it also showed good MM-GBSA compared to the natural inhibitor N3. Conclusion: Our results could be auspicious for fast repurposing of the examined drugs either alone or in combinations with each other for the treatment of the COVID-19. Furthermore, this work provides a clear spot on the structure-activity relationship (SAR) for targeting the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) main protease and helps the design and synthesis of new drugs in the future targeting it as well.


2021 ◽  
Vol 261 ◽  
pp. 02069
Author(s):  
Qiuxian Zhang ◽  
Hecheng Wang

Alzheimer’s disease (AD) is an age-related neurodegenerative disorder. β-amyloid protein (Aβ) is the key protein which involved in AD. But the physiological function of Aβ is needed to be investigated. Many experimental studies have shown that Aβ could bind to glycoproteins D (gD) on the surface of the herpes virus. However the mechanism is still unclear. In the present study, we elucidate the molecular mechanism of the interaction between Aβ and gD of herpes simplex virus type 2 (HSV-2) by molecular docking and molecular dynamics simulation. Molecular dynamics simulations displayed that Aβ could stably bind to the HSV-2 gD owing to the presence of several interactions. Analysis binding free energy by molecular mechanics Poisson–Boltzmann surface area (MM–PBSA) method revealed that hot residues including Glu3, Glu11, Glu22 and Ala42 of Aβ1-42 were involved in binding with HSV-2 gD. Thus, the HSV-2 gD can be entrapped by Aβ which will be utilized for prevent and therapy of AD in future.


2019 ◽  
Vol 4 (3) ◽  
pp. 47-50
Author(s):  
Reihaneh Sabbaghzadeh

Cyclin-dependent kinases (CDKs, play important roles in cell cycle regulation. Since deregulation of cyclins and/or alteration or absence of inhibitors has been associated with many cancers, there is strong interest in CDKs inhibitors that could play an necessary role in the discovery of a new family of antitumor agents. molecular modeling is used to design new materials, which the correct prediction of physical properties of realistic systems is required. Gromacs is an tool to perform molecular dynamics simulations and energy minimization of bimolecular systems which commonly consists of several tens to thousands of amino acid residues. From the simulations the first 100 ps were regarded as equilibration, leaving 900 ps for analysis purposes. Thermodynamic properties density, potential energy, temperature, and pressure are given. In this results tetramethyl((1,4-phenylenebis(azanediyl))bis((4-chlorophenyl)methylene))bis(phosphonate) was found to be the better selective known inhibitor for cyclin-dependent kinase2 because it shown lowest energy. The phosphoric acid moiety is considered to bind to the affected protein more actively than the corresponding carboxylic acid because of its di anionic character.


Sign in / Sign up

Export Citation Format

Share Document