MOLECULAR DOCKING AND MOLECULAR DYNAMICS STUDIES OF THE INTERACTION BETWEEN COUMARIN-NEUROTRANSMITTER DERIVATIVES AND CARBONIC ANHYDRASE IX

2021 ◽  
Author(s):  
Dušan Dimić ◽  
◽  
Dejan Milenković ◽  
Edina Avdović ◽  
Goran Kaluđerović ◽  
...  

Novel biologically active compounds can be obtained by the structural modification of coumarins. In this contribution, five new derivatives of 4-hydroxycoumarin with tyramine, octopamine, norepinephrine, 3-methoxytyramine, and dopamine were obtained. Their structures were optimized based on the previously obtained crystal structure of the 4-hydroxycoumarin-dopamine derivative. The special emphasis was put on the effect of various substituents on the structure of obtained compounds and intramolecular interactions governing the stability. To investigate their possible antitumor activity, molecular docking and molecular dynamics simulations were performed with Carbonic anhydrase, a prognostic factor in several cancers, and compared to the native ligand, 5-acetamido-1,3,4-thiadiazole- 2-sulfonamide. The results have shown that all of the coumarin-neurotransmitter derivatives bind to the active pocket of protein with the binding energies higher than for the native ligand. The main contributions to the binding energies were discussed. The Root Mean Square Deviation (RMSD), Root Mean Square Fluctuation (RMSF), and Radius of gyration (Rg), as results of MD simulations, were used to predict the activity of compounds towards chosen protein. The highest MD binding energies were obtained for the derivatives with dopamine and 3-methoxytyramine, with the van der Waals interaction and hydrogen bonds being the most important contributors.

2021 ◽  
Vol 45 (1) ◽  
Author(s):  
Mohammad Mahfuz Ali Khan Shawan ◽  
Sajal Kumar Halder ◽  
Md. Ashraful Hasan

Abstract Background At present, the entire world is in a war against COVID-19 pandemic which has gradually led us toward a more compromised “new normal” life. SARS-CoV-2, the pathogenic microorganism liable for the recent COVID-19 outbreak, is extremely contagious in nature resulting in an unusual number of infections and death globally. The lack of clinically proven therapeutic intervention for COVID-19 has dragged the world’s healthcare system into the biggest challenge. Therefore, development of an efficient treatment scheme is now in great demand. Screening of different biologically active plant-based natural compounds could be a useful strategy for combating this pandemic. In the present research, a collection of 43 flavonoids of 7 different classes with previously recorded antiviral activity was evaluated via computational and bioinformatics tools for their impeding capacity against SARS-CoV-2. In silico drug likeness, pharmacophore and Absorption, Distribution, Metabolism, Excretion and Toxicity (ADMET) profile analysis of the finest ligands were carried out using DataWarrior, DruLiTo and admetSAR programs, respectively. Molecular docking was executed by AutoDock Vina, while molecular dynamics simulation of the target protein–ligand bound complexes was done using nanoscalable molecular dynamics and visual molecular dynamics software package. Finally, the molecular target analysis of the selected ligands within Homo sapiens was conducted with SwissTargetPredcition web server. Results Out of the forty-three flavonoids, luteolin and abyssinone II were found to develop successful docked complex within the binding sites of target proteins in terms of lowest binding free energy and inhibition constant. The root mean square deviation and root mean square fluctuation values of the docked complex displayed stable interaction and efficient binding between the ligands and target proteins. Both of the flavonoids were found to be safe for human use and possessed good drug likeness properties and target accuracy. Conclusions Conclusively, the current study proposes that luteolin and abyssinone II might act as potential therapeutic candidates for SARS-CoV-2 infection. In vivo and in vitro experiments, however, should be taken under consideration to determine the efficiency and to demonstrate the mechanism of action.


Molecules ◽  
2021 ◽  
Vol 26 (21) ◽  
pp. 6538
Author(s):  
Md. Mominur Rahman ◽  
Md. Junaid ◽  
S. M. Zahid Hosen ◽  
Mohammad Mostafa ◽  
Lei Liu ◽  
...  

Inflammation plays an important role in different chronic diseases. Brominated indoles derived from the Australian marine mollusk Dicathais orbita (D. orbita) are of interest for their anti-inflammatory properties. This study evaluates the binding mechanism and potentiality of several brominated indoles (tyrindoxyl sulfate, tyrindoleninone, 6-bromoisatin, and 6,6′-dibromoindirubin) against inflammatory mediators cyclooxygenases-1/2 (COX-1/2) using molecular docking, followed by molecular dynamics simulation, along with physicochemical, drug-likeness, pharmacokinetic (pk), and toxicokinetic (tk) properties. Molecular docking identified that these indole compounds are anchored, with the main amino acid residues, positioned in the binding pocket of the COX-1/2, required for selective inhibition. Moreover, the molecular dynamics simulation based on root mean square deviation (RMSD), radius of gyration (Rg), solvent accessible surface area (SASA), and root mean square fluctuation (RMSF) analyses showed that these natural brominated molecules transit rapidly to a progressive constant configuration during binding with COX-1/2 and seem to accomplish a consistent dynamic behavior by maintaining conformational stability and compactness. The results were comparable to the Food and Drug Administration (FDA)-approved selective COX inhibitor, aspirin. Furthermore, the free energy of binding for the compounds assessed by molecular mechanics–Poisson–Boltzmann surface area (MM–PBSA) confirmed the binding capacity of indoles towards COX-1/2, with suitable binding energy values except for the polar precursor tyrindoxyl sulfate (with COX-1). The physicochemical and drug-likeness analysis showed zero violations of Lipinski’s rule, and the compounds are predicted to have excellent pharmacokinetic profiles. These indoles are projected to be non-mutagenic and free from hepatotoxicity, with no inhibition of human ether-a-go–go gene (hERG) I inhibitors, and the oral acute toxicity LD50 in rats is predicted to be similar or lower than aspirin. Overall, this work has identified a plausible mechanism for selective COX inhibition by natural marine indoles as potential therapeutic candidates for the mitigation of inflammation.


2016 ◽  
Vol 15 ◽  
pp. CIN.S38044 ◽  
Author(s):  
Merlin Lopus ◽  
D. Meshach Paul ◽  
R. Rajasekaran

Tumor suppressor gene, STK11, encodes for serine-threonine kinase, which has a critical role in regulating cell growth and apoptosis. Mutations of the same lead to the inactivation of STK11, which eventually causes different types of cancer. In this study, we focused on identifying those driver mutations through analyzing structural variations of mutants, viz., D194N, E199K, L160P, and Y49D. Native and the mutants were analyzed to determine their geometrical deviations such as root-mean-square deviation, root-mean-square fluctuation, radius of gyration, potential energy, and solvent-accessible surface area using conformational sampling technique. Additionally, the global minimized structure of native and mutants was further analyzed to compute their intramolecular interactions and distribution of secondary structure. Subsequently, simulated thermal denaturation and docking studies were performed to determine their structural variations, which in turn alter the formation of active complex that comprises STK11, STRAD, and MO25. The deleterious effect of the mutants would result in a comparative loss of enzyme function due to variations in their binding energy pertaining to spatial conformation and flexibility. Hence, the structural variations in binding energy exhibited by the mutants, viz., D194N, E199K, L160P, and Y49D, to that of the native, consequently lead to pathogenesis.


2020 ◽  
Author(s):  
Sk. Md Na ◽  
M. Srinivasa R

Abstract Molecular Dynamics simulation using Gromacs with OPLS-AA force field is performed for 100ns between SARS-CoV-2 main protease and Dexamethasone / Umifenovir drugs at 300 K/1 atm pressure. The trajectory of Root Mean Square Deviation (RMSD) and Radius of Gyration(Rg) emphasized the achievement of equilibrium and compactness. The drug-binding affinities on SARS-CoV-2 main protease are estimated via MM/PBSA method. The sign with magnitude of computed Gibbs free energy indicated the presence of strong interactions between SARS-CoV-2 and drugs of Dexamethasone / Umifenovir. The study revealed that the drug Dexamethasone is more effective over Umifenovir in binding SARS-CoV-2 main protease.


Molecules ◽  
2021 ◽  
Vol 26 (17) ◽  
pp. 5304
Author(s):  
Mohammad G. Al-Thiabat ◽  
Amirah Mohd Gazzali ◽  
Noratiqah Mohtar ◽  
Vikneswaran Murugaiyah ◽  
Ezatul Ezleen Kamarulzaman ◽  
...  

Drug targeting is a progressive area of research with folate receptor alpha (FRα) receiving significant attention as a biological marker in cancer drug delivery. The binding affinity of folic acid (FA) to the FRα active site provides a basis for recognition of FRα. In this study, FA was conjugated to beta-cyclodextrin (βCD) and subjected to in silico analysis (molecular docking and molecular dynamics (MD) simulation (100 ns)) to investigate the affinity and stability for the conjugated system compared to unconjugated and apo systems (ligand free). Docking studies revealed that the conjugated FA bound into the active site of FRα with a docking score (free binding energy < −15 kcal/mol), with a similar binding pose to that of unconjugated FA. Subsequent analyses from molecular dynamics (MD) simulations, root mean square deviation (RMSD), root mean square fluctuation (RMSF), and radius of gyration (Rg) demonstrated that FA and FA–βCDs created more dynamically stable systems with FRα than the apo-FRα system. All systems reached equilibrium with stable RMSD values ranging from 1.9–2.4 Å and the average residual fluctuation values of the FRα backbone atoms for all residues (except for terminal residues ARG8, THR9, THR214, and LEU215) were less than 2.1 Å with a consistent Rg value of around 16.8 Å throughout the MD simulation time (0–100 ns). The conjugation with βCD improved the stability and decreased the mobility of all the residues (except residues 149–151) compared to FA–FRα and apo-FRα systems. Further analysis of H-bonds, binding free energy (MM-PBSA), and per residue decomposition energy revealed that besides APS81, residues HIS20, TRP102, HIS135, TRP138, TRP140, and TRP171 were shown to have more favourable energy contributions in the holo systems than in the apo-FRα system, and these residues might have a direct role in increasing the stability of holo systems.


Molecules ◽  
2020 ◽  
Vol 25 (18) ◽  
pp. 4094
Author(s):  
Yi Zhang ◽  
Ting-jian Zhang ◽  
Shun Tu ◽  
Zhen-hao Zhang ◽  
Fan-hao Meng

Src plays a crucial role in many signaling pathways and contributes to a variety of cancers. Therefore, Src has long been considered an attractive drug target in oncology. However, the development of Src inhibitors with selectivity and novelty has been challenging. In the present study, pharmacophore-based virtual screening and molecular docking were carried out to identify potential Src inhibitors. A total of 891 molecules were obtained after pharmacophore-based virtual screening, and 10 molecules with high docking scores and strong interactions were selected as potential active molecules for further study. Absorption, distribution, metabolism, elimination and toxicity (ADMET) property evaluation was used to ascertain the drug-like properties of the obtained molecules. The proposed inhibitor–protein complexes were further subjected to molecular dynamics (MD) simulations involving root-mean-square deviation and root-mean-square fluctuation to explore the binding mode stability inside active pockets. Finally, two molecules (ZINC3214460 and ZINC1380384) were obtained as potential lead compounds against Src kinase. All these analyses provide a reference for the further development of novel Src inhibitors.


2016 ◽  
Vol 14 (1) ◽  
pp. 188-196 ◽  
Author(s):  
Syed Lal Badshah ◽  
Abdul Naeem Khan ◽  
Yahia Nasser Mabkhot

AbstractA molecular dynamics (MD) simulation study of the enzymatic portion of cholera toxin; cholera toxin A-1 polypeptide (CTA1) was performed at 283, 310 and 323 K. From total energy analysis it was observed that this toxin is stable thermodynamically and these outcomes were likewise confirmed by root mean square deviations (RMSD) investigations. The Cα root mean square fluctuation (RMSF) examinations revealed that there are a number of residues inside CTA1, which can be used as target for designing and synthesizing inhibitory drugs, in order to inactivate cholera toxin inside the human body. The fluctuations in the radius of gyration and hydrogen bonding in CTA1 proved that protein unfolding and refolding were normal routine phenomena in its structure at all temperatures. Solvent accessible surface area study identified the hydrophilic nature of the CTA1, and due to this property it can be a potential biological weapon. The structural identification (STRIDE) algorithm for proteins was successfully used to determine the partially disordered secondary structure of CTA1. On account of this partially disordered secondary structure, it can easily deceive the proteolytic enzymes of the endoplasmic reticulum of host cells.


2021 ◽  
Vol 13 (1) ◽  
Author(s):  
Andrew T. McNutt ◽  
Paul Francoeur ◽  
Rishal Aggarwal ◽  
Tomohide Masuda ◽  
Rocco Meli ◽  
...  

AbstractMolecular docking computationally predicts the conformation of a small molecule when binding to a receptor. Scoring functions are a vital piece of any molecular docking pipeline as they determine the fitness of sampled poses. Here we describe and evaluate the 1.0 release of the Gnina docking software, which utilizes an ensemble of convolutional neural networks (CNNs) as a scoring function. We also explore an array of parameter values for Gnina 1.0 to optimize docking performance and computational cost. Docking performance, as evaluated by the percentage of targets where the top pose is better than 2Å root mean square deviation (Top1), is compared to AutoDock Vina scoring when utilizing explicitly defined binding pockets or whole protein docking. Gnina, utilizing a CNN scoring function to rescore the output poses, outperforms AutoDock Vina scoring on redocking and cross-docking tasks when the binding pocket is defined (Top1 increases from 58% to 73% and from 27% to 37%, respectively) and when the whole protein defines the binding pocket (Top1 increases from 31% to 38% and from 12% to 16%, respectively). The derived ensemble of CNNs generalizes to unseen proteins and ligands and produces scores that correlate well with the root mean square deviation to the known binding pose. We provide the 1.0 version of Gnina under an open source license for use as a molecular docking tool at https://github.com/gnina/gnina.


2021 ◽  
Vol 12 (6) ◽  
pp. 7239-7248

The novel coronavirus, recognized as COVID-19, is the cause of an infection outbreak in December 2019. The effect of temperature and pH changes on the main protease of SARS-CoV-2 were investigated using all-atom molecular dynamics simulation. The obtained results from the root mean square deviation (RMSD) and root mean square fluctuations (RMSF) analyses showed that at a constant temperature of 25℃ and pH=5, the conformational change of the main protease is more significant than that of pH=6 and 7. Also, by increasing temperature from 25℃ to 55℃ at constant pH=7, a remarkable change in protein structure was observed. The radial probability of water molecules around the main protease was decreased by increasing temperature and decreasing pH. The weakening of the binding energy between the main protease and water molecules due to the increasing temperature and decreasing pH has reduced the number of hydrogen bonds between the main protease and water molecules. Finding conditions that alter the conformation of the main protease could be fundamental because this change could affect the virus’s functionality and its ability to impose illness.


2020 ◽  
Vol 2 ◽  
pp. 382
Author(s):  
C. G. Koutroulos

Using as an example potentials of the form U±(r)=-D±(cosh^2(r/R))^-1 the binding energies as well as the root mean square radii of the orbits of the Λ particle in hypernuclei in the ground and excited states were calculated in the relativistic and non-relativistic cases and the results are compared.


Sign in / Sign up

Export Citation Format

Share Document