Insights into the structure and inhibition of Giardia intestinalis arginine deiminase: homology modeling,docking, and molecular dynamics studies

Giardia intestinalis arginine deiminase (GiADI) is an important metabolic enzyme involved in the energy production and defense of this protozoan parasite. The lack of this enzyme in the human host makes GiADI an attractive target for drug design against G. intestinalis. One approach in the design of inhibitors of GiADI could be computer-assisted studies of its crystal structure, such as docking; however, the required crystallographic structure of the enzyme still remains unresolved. Because of its relevance, in this work, we present a three-dimensional structure of GiADI obtained from its amino acid sequence using the homology modeling approximation. Furthermore, we present an approximation of the most stable dimeric structure of GiADI identified through molecular dynamics simulation studies. An in silico analysis of druggability using the structure of GiADI was carried out in order to know if it is a good target for design and optimization of selective inhibitors. Potential GiADI inhibitors were identified by docking of a set of 3196 commercial and 19 in-house benzimidazole derivatives, and molecular dynamics simulation studies were used to evaluate the stability of the ligand–enzyme complexes.     

Development of a potent 2-oxoamide inhibitor of secreted phospholipase A2 guided by molecular docking calculations and molecular dynamics simulations

Inhibition of group IIA secreted phospholipase A₂ (GIIA sPLA₂) has been an important objective for medicinal chemists. We have previously shown that inhibitors incorporating the 2-oxoamide functionality may inhibit human and mouse GIIA sPLA₂s. Herein, the development of new potent inhibitors by molecular docking calculations using the structure of the known inhibitor 7 as scaffold, are described. Synthesis and biological evaluation of the new compounds revealed that the long chain 2-oxoamide based on (S)-valine GK241 led to improved activity (IC₅₀ = 143 nM and 68 nM against human and mouse GIIA sPLA₂, respectively). In addition, molecular dynamics simulations were employed to shed light on GK241 potent and selective inhibitory activity.     

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.     

Exploring blocker binding to a homology model of the open hERG K+ channel using docking and molecular dynamics methods

Binding of blockers to the human voltage-gated hERG potassium channel is studied using a combination of homology modelling, automated docking calculations and molecular dynamics simulations, where binding affinities are evaluated using the linear interaction energy method. A homology model was constructed based on the available crystal structure of the bacterial KvAP channel and the affinities of a series of sertindole analogues predicted using this model. The calculations reproduce the relative binding affinities of these compounds very well and indicate that both polar interactions near the intracellular opening of the selectivity filter as well as hydrophobic complementarity in the region around F656 are important for blocker binding. These results are consistent with recent alanine scanning mutation experiments on the blocking of the hERG channel by other compounds.     

Reduced catalytic activity of human CYP2C9 natural alleles for gliclazide: molecular dynamics simulation and docking studies

Amongst sulfonylureas, gliclazide is one of the mostly prescribed drugs to diabetic patients and is metabolized extensively by P450 CYP2C9. Among 24-CYP2C9 alleles, the *2/*2 and *3/*3 genotypes showed significantly lower gliclazide clearances with reductions of 25 and 57%, respectively. However, the reason for the change in drug-metabolizing activity induced by these natural alleles is unknown. In the present study, we used molecular dynamics simulation and autodocking studies to provide models for gliclazide-bound complexes of CYP2C9*2, *3 and *2/*3 mutants, which give insight into CYP2C9-gliclazide interactions and explain the reduced enzymatic activity seen in these variants. Our data shows that the size of the substrate-access entry site is significantly reduced in mutants, which limits the access of gliclazide to heme and the active site. The distance from the substrate oxidation site and heme is >5 Å in *3 and *2/*3. Therefore, the addition of an active oxygen molecule by heme-Fe is hindered. The absence of F100, F114 and F476 in the interacting amino acid pocket in *3 reduces catalytic efficiency toward gliclazide. In *1, gliclazide is stabilized by the formation of two hydrogen bonds with R108 while it is absent in mutants. Further in *3 and *2/*3, the key heme-stabilizing residue, R97 stabilization is greatly reduced. Therefore, the decreased catalytic activity of these variants can be explained from the reduced access of the gliclazide to heme, and the interaction between heme and substrate is affected due to their instability in the active site.     

Metronidazole aryloxy,carboxy and azole derivatives: Synthesis,anti-tumor activity,QSAR, molecular docking and dynamics studies

A series of novel metronidazole aryloxy, carboxy and azole derivatives has been synthesized and their cytotoxic activities on three cancer cell lines were evaluated by MTT assay. Compounds 4m, 4l and 4d showed the most potent cytotoxic activity (IC₅₀s?less than?100?µg/mL). Apoptosis was also detected for these compounds by flow cytometry. Docking studies were performed in order to propose the probable target protein. In the next step, molecular dynamics simulation was carried out on the proposed target protein, focal adhesion kinase (FAK, PDB code: 2ETM), bound to compound 4m. As, 4m showed a potent cytotoxic activity and an acceptable apoptotic effect, it can be a potential anticancer candidate that may work through inhibition of FAK.     

Regioselectivity of CYP2B6: homology modeling,molecular dynamics simulation,docking

Human cytochrome P450 (CYP) 2B6 activates the anticancer prodrug cyclophosphamide (CPA) by 4-hydroxylation. In contrast, the same enzyme catalyzes N-deethylation of a structural isomer, the prodrug ifosfamide (IFA), thus causing severe adverse drug effects. To model the molecular interactions leading to a switch in regioselectivity, the structure of CYP2B6 was modeled based on the structure of rabbit CYP2C5. We modeled the missing 22-residue loop in CYP2C5 between helices F and G (the F-G loop), which is not resolved in the X-ray structure, by molecular dynamics (MD) simulations using a simulated annealing protocol. The modeled conformation of the loop was validated by unconstrained MD simulations of the complete enzymes (CYP2C5 and CYP2B6) in water for 70 and 120 ps, respectively. The simulations were stable and led to a backbone r.m.s. deviation of 1.7 A between the two CYPs.The shape of the substrate binding site of CYP2B6 was further analyzed. It consists of three well-defined hydrophobic binding pockets adjacent to the catalytic heme. Size, shape and hydrophobicity of these pockets were compared to the shapes of the two structurally isomeric substrates. In their preferred orientation in the binding site, both substrates fill all three binding pockets without repulsive interactions. The distance to the heme iron is short enough for 4-hydroxylation and N-deethylation to occur for CPA and IFA, respectively. However, if the substrates are docked in the non-preferred orientation (such that 4-hydroxylation and N-deethylation would occur for IFA and CPA, respectively), one pocket is left empty, and clashes were observed between the substrates and the enzyme.     

Molecular modeling studies of atorvastatin analogues as HMGR inhibitors using 3D-QSAR,molecular docking and molecular dynamics simulations

3-Hydroxy-3-methylglutaryl coenzyme-A reductase (HMGR) is generally regarded as targets for the treatment of hypercholesterolemia. HMGR inhibitors (more commonly known as statins) are discovered as plasma cholesterol lowering molecules. In this work, 120 atorvastatin analogues were studied using a combination of molecular modeling techniques including three-dimensional quantitative structure–activity relationship (3D-QSAR), molecular docking and molecular dynamics (MD) simulation. The results show that the best CoMFA (comparative molecular field analysis) model has q² = 0.558 and r² = 0.977, and the best CoMSIA (comparative molecular similarity indices analysis) model has q² = 0.582 and r² = 0.919. Molecular docking and MD simulation explored the binding relationship of the ligand and the receptor protein. The calculation results indicated that the hydrophobic and electrostatic fields play key roles in QSAR model. After MD simulation, we found four vital residues (Lys735, Arg590, Asp690 and Asn686) and three hydrophobic regions in HMGR binding site. The calculation results show that atorvastatin analogues obtained by introduction of F atoms or gem-difluoro groups could obviously improve the inhibitory activity. The new HMGR inhibitor analogues design in this Letter had been submitted which is being currently synthesized by our laboratories.     

Insights into the structure of the PmrD protein with molecular dynamics simulations

Resistance to cationic antimicrobial peptide polymyxin B from Gram-negative bacteria is accomplished by two-component systems (TCSs), protein complexes PmrA/PmrB and PhoP/PhoQ. PmrD is the first protein identified to mediate the connectivity between two TCSs. The 3D structure of PmrD has been recently solved by NMR and its unique fold was revealed. Here, a molecular dynamics study is presented started from the NMR structure. Numerous hydrophobic and electrostatic interactions were identified to contribute to PmrD's 3D stability. Moreover, the mobility of the five loops that connect the protein's six β-strands has been explored. Solvent-accessible surface area calculation revealed that a Leucine-rich hydrophobic cluster of the protein stabilized the protein's structure.     

Combined docking and molecular dynamics simulations to enlighten the capacity of Pseudomonas cepacia and Candida antarctica lipases to catalyze quercetin acetylation

A combined docking and molecular dynamics protocol was applied to investigate quercetin binding modes within the catalytic cavity of Candida antarctica lipase B (CALB) and Pseudomonas cepacia lipase (PCL), aiming to explain the difference of specificity of these enzymes in acetylation reaction. For both lipases, docking of quercetin yielded two families of conformers with either the quercetin A or B-ring pointing towards the catalytic residues. Molecular dynamics (MD) calculations were subsequently performed on several complexes of each family. MD trajectories were analyzed focusing on the orientation of the acyl donor bound to the catalytic serine towards the oxyanion hole residues and the proximity of quercetin hydroxyl groups to the catalytic residues. Results showed that with CALB, the acetate was not correctly positioned within the oxyanion hole whatever the orientation of quercetin, suggesting that no product could be obtained. With PCL, the acetate remained within the oxyanion hole during all MD trajectories. Depending on quercetin orientation, either the 7-OH group or the 3, 5, 3′, 4′-OH groups came alternatively near the catalytic residues, suggesting that all of them could be acylated. The capacity of models to explain the regioselectivity of the reaction was discussed. Key residues and interactions involved in quercetin binding modes were identified and related to the reaction feasibility.     

Structural study of biologically significant ligands with major birch pollen allergen Betv1 by docking and molecular dynamics simulation

The major birch pollen allergen, Betv1 of Betula verrucosa is the main causative agent of birch pollen allergy in humans. Betv1 is capable of binding several physiological ligands including fatty acids, flavones, cytokinins and sterols. Until now, no structural information from crystallography or NMR is available regarding binding mode of any of these ligands into the binding pocket of Betv1. In the present study thirteen ligands have been successfully docked into the hydrophobic cavity of Betv1 and binding free energies of the complexes have been calculated using AutoDock 3.0.5. A linear relationship with correlation coefficient (R2) of 0.6 is obtained between ΔG(b)s values plotted against their corresponding IC50 values. The complex formed between Betv1 and the best docking pose for each ligand has been optimized by molecular dynamics simulation. Here, we describe the ligand binding of Betv1, which provides insight into the biological function of this protein. This knowledge is required for structural alteration or inhibition of some of these ligands in order to modify the allergenic properties of this protein.     

Structural insights into human microsomal epoxide hydrolase by combined homology modeling,molecular dynamics simulations,and molecular docking calculations

A new homology model of human microsomal epoxide hydrolase was derived based on multiple templates. The model obtained was fully evaluated, including MD simulations and ensemble‐based docking, showing that the quality of the structure is better than that of only previously known model. Particularly, a catalytic triad was clearly identified, in agreement with the experimental information available. Analysis of intermediates in the enzymatic mechanism led to the identification of key residues for substrate binding, stereoselectivity, and intermediate stabilization during the reaction. In particular, we have confirmed the role of the oxyanion hole and the conserved motif (HGXP) in epoxide hydrolases, in excellent agreement with known experimental and computational data on similar systems. The model obtained is the first one that fully agrees with all the experimental observations on the system. Proteins 2017; 85:720–730. © 2016 Wiley Periodicals, Inc.     

Design of novel PhMTNA inhibitors,targeting neurological disorder through homology modeling,molecular docking,and dynamics approaches

Vanishing white matter (VWM) is a hereditary human disease, mostly prevalent in childhood caused by the defects in the eukaryotic initiation factor beta subunits. It is the first disease involved in the translation initiation factor, eIF2B. There is no specific treatment for VWM which mainly affect the brain and ovaries. The gray matter remains normal in all characteristics while the white matter changes texture, coming to the pathophysiology, many initiation factors are involved in the initiation of translation of mRNAs into polypeptides. In this study, the three-dimensional structure of PhMTNA protein was modeled and the stability ascertained through Molecular dynamic simulation (MDS) for 100?ns. The active site residues are conserved with the reported BsMTNA structure which is also confirmed through sitemap prediction. Through virtual screening and induced fit docking, top five leads against PhMTNA protein was identified based on their binding mode and affinity. ADME properties and DFT (Density Functional Theory) studies of these compounds were studied. In addition to that, computational mutagenesis studies were performed to identify the hotspot residues involved in the protein-ligand interactions. Overall analysis showed that the compound NCI_941 has a highest binding energy of ?46.256?kcal mol^?1 in the Arg57Ala mutant. Thus, the results suggest that NCI_941 would act as a potent inhibitor against PhMTNA protein.     

Coumarinyl pyranopyrimidines as new neuropeptide S receptor antagonists; design,synthesis, homology and molecular docking

In this work, we described the design, synthesis and characterization of a new class of NPSR antagonists bearing the tetracyclic coumarinyl pyranopyrimidine scaffold incorporated with different acyclic and/or heterocyclic moieties. These compounds are highlighted in this study as never being used as NPSR antagonists before which provides a model for the discovery of new bioactive inhibitors that may hold potential for drug development towards anxiety, food, and addiction disorders. Synthetic and medicinal chemistry studies led to the identification of four potent antagonists, compounds 7d, 10, 12 and 13, which were able to significantly inhibit the stimulatory effect of NPS through counteracting the increased intracellular Ca²⁺ accumulation. The target compound 7d was the most active derivative behaving as a pure NPSR antagonist and displaying IC₅₀ value of 2 μM. Homology model of NPSR was built based on bovine rhodopsin structure. Modeling studies were carried out to further rationalize the NPSR binding mode of the target compounds. Moreover, molecular dynamics simulation study was performed for compounds 7d, 10 and 12 which revealed the stability of the ligand-protein complex and the reliability of the docking studies.     

Drug repurposing for SARS-CoV-2 main protease: Molecular docking and molecular dynamics investigations

The current novel corona virus illness (COVID-19) is a developing viral disease that was discovered in 2019. There is currently no viable therapeutic strategy for this illness management. Because traditional medication development and discovery has lagged behind the threat of emerging and re-emerging illnesses like Ebola, MERS-CoV, and, more recently, SARS-CoV-2. Drug developers began to consider drug repurposing (or repositioning) as a viable option to the more traditional drug development method. The goal of drug repurposing is to uncover new uses for an approved or investigational medicine that aren't related to its original use. The main benefits of this strategy are that there is less developmental risk and that it takes less time because the safety and pharmacologic requirements are met. The main protease (Mpro) of corona viruses is one of the well-studied and appealing therapeutic targets. As a result, the current research examines the molecular docking of Mpro (PDB ID: 5R81) conjugated repurposed drugs. 12,432 approved drugs were collected from ChEMBL and drugbank libraries, and docked separately into the receptor grid created on 5R81, using the three phases of molecular docking including high throughput virtual screening (HTVS), standard precision (SP), and extra precision (XP). Based on docking scores and MM-GBSA binding free energy calculation, top three drugs (kanamycin, sulfinalol and carvedilol) were chosen for further analyses for molecular dynamic simulations.     

On the human CYP2C9*13 variant activity reduction: a molecular dynamics simulation and docking study

Cytochrome P450 2C9 (CYP2C9) plays a key role in the metabolism of clinical drugs. CYP2C9 is a genetically polymorphic enzyme and some of its allelic variants have less activity compared to the wild-type form. Drugs with a narrow therapeutic index may cause serious toxicity to the individuals who carry such allele. CYP2C9*13, firstly identified by some of the present authors in a Chinese poor metabolizer of lornoxicam, is characterized by mutation encoding Leu90Pro substitution. Kinetic experiments show that CYP2C9*13 has less catalytic activity in elimination of diclofenac and lornoxicam in vitro. In order to explore the structure-activity relationship of CYP2C9*13, the three-dimensional structure models of the substrate-free CYP2C9*1 and its variant CYP2C9*13 are constructed on the basis of the X-ray crystal structure of human CYP2C9*1 (PDB code 1R9O) by molecular dynamics simulations. The structure change caused by Leu90Pro replacement is revealed and used to explain the dramatic decrease of the enzymatic activity in clearance of the two CYP2C9 substrates: diclofenac and lornoxicam. The trans configuration of the bond between Pro90 and Asp89 in CYP2C9*13 is firstly identified. The backbone of residues 106-108 in CYP2C9*13 turns over and their side chains block the entrance for substrates accessing so that the entrance of *13 shrinks greatly than that in the wild-type, which is believed to be the dominant mechanism of the catalytic activity reduction. Consequent docking study which is consistent with the results of the kinetic experiments by Guo et al. identifies the most important residues for enzyme-substrate complexes.     

Synthesis,molecular docking and cholinesterase inhibitory activity of hydroxylated 2-phenylbenzofuran derivatives

We have designed, synthesized and evaluated a series of hydroxylated 2-phenylbenzofuran derivatives as potential cholinesterase inhibitors. Starting from a series of 2-phenylbenzofurans previously published, in this paper we present a complete synthesis and the influence on the activity of one or two hydroxyl groups located in meta or in meta and para positions respectively of the 2-phenyl ring and highlight the importance of position of hydroxyl groups. Moreover, simultaneous introduction of halogen at position 7 of the benzofuran scaffold resulted in an improved inhibitory activity against the enzyme. To further provide molecular insight and to identify the most probable ligand-binding site of the protein, docking studies were performed for the top-ranked compounds. Docking results revealed conserved ligand-binding residues and supported the role of catalytic site residues in enzyme inhibition.     

3D-QSAR (CoMFA,CoMSIA), molecular docking and molecular dynamics simulations study of 6-aryl-5-cyano-pyrimidine derivatives to explore the structure requirements of LSD1 inhibitors

Recently, Histone Lysine Specific Demethylase 1 (LSD1) was regarded as a promising anticancer target for the novel drug discovery. And several small molecules as LSD1 inhibitors in different structures have been reported. In this work, we carried out a molecular modeling study on the 6-aryl-5-cyano-pyrimidine fragment LSD1 inhibitors using three-dimensional quantitative structure–activity relationship (3D-QSAR), molecular docking and molecular dynamics simulations. Comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA) were used to generate 3D-QSAR models. The results show that the best CoMFA model has q² = 0.802, r²_ncv = 0.979, and the best CoMSIA model has q² = 0.799, r²_ncv = 0.982. The electrostatic, hydrophobic and H-bond donor fields play important roles in the models. Molecular docking studies predict the binding mode and the interactions between the ligand and the receptor protein. Molecular dynamics simulations results reveal that the complex of the ligand and the receptor protein are stable at 300 K. All the results can provide us more useful information for our further drug design.     

Structural studies,homology modeling and molecular docking of novel non-competitive antagonists of GluK1/GluK2 receptors

Non-competitive ligands of kainate receptors have focused significant attention as medicinal compounds because they seem to be better tolerated than competitive antagonists and uncompetitive blocker of these receptors. Here we present structural studies (X-ray structure determination, NMR and MS spectra) of novel indole-derived non-competitive antagonists of GluK1/GluK2 receptors, homology models of GluK1 and GluK2 receptors based on novel AMPA receptor template as well as molecular docking of ligands to their molecular targets. We find that the allosteric site is in the receptor transduction domain, in one receptor subunit, not between the two subunits as it was indicated by our earlier studies.     

Multiple templates-based homology modeling enhances structure quality of AT1 receptor: validation by molecular dynamics and antagonist docking

We present a comparative account on 3D-structures of human type-1 receptor (AT1) for angiotensin II (AngII), modeled using three different methodologies. AngII activates a wide spectrum of signaling responses via the AT1 receptor that mediates physiological control of blood pressure and diverse pathological actions in cardiovascular, renal, and other cell types. Availability of 3D-model of AT1 receptor would significantly enhance the development of new drugs for cardiovascular diseases. However, templates of AT1 receptor with low sequence similarity increase the complexity in straightforward homology modeling, and hence there is a need to evaluate different modeling methodologies in order to use the models for sensitive applications such as rational drug design. Three models were generated for AT1 receptor by, (1) homology modeling with bovine rhodopsin as template, (2) homology modeling with multiple templates and (3) threading using I-TASSER web server. Molecular dynamics (MD) simulation (15 ns) of models in explicit membrane-water system, Ramachandran plot analysis and molecular docking with antagonists led to the conclusion that multiple template-based homology modeling outweighs other methodologies for AT1 modeling.