Refined homology model of cytochrome bcc complex B subunit for virtual screening of potential anti-tuberculosis agents

Abstract

Cytochrome bcc complex is important for ATP synthesis and cellular activity, as a crucial step in the terminal reduction of oxygen in aerobic electron transport chains. The b subunit of cytochrome bcc complex (QcrB) has been reported as a promising anti-tuberculosis target, with many novel anti-tuberculosis scaffolds reported. However, the 3D structure of mycobacterium tuberculosis (M. tuberculosis) QcrB has not been released, making it hard to understand the interactions between QcrB and its inhibitors as well as to develop novel anti-tuberculosis scaffolds. Herein we built the optimal homology model of M. tuberculosis QcrB using the M. smegmatis QcrB structure as template, which was refined through all-atom molecular dynamics simulation. Then, the binding modes of known inhibitors were predicted through molecular docking method, along with molecular dynamics simulation and binding free energy calculation to verify the accuracy of docking results and stability of the protein-inhibitor complexes. The informative key residues within QcrB site enabled us to perform structure-based virtual library screening to obtain potential M. tuberculosis QcrB inhibitors, which were validated through molecular dynamics simulation and MM-GBSA calculation and analyzed through pharmacokinetic properties prediction. Our research would provide a deeper insight into the interactions between M. tuberculosis QcrB and its inhibitors, which boosts to develop novel therapy against tuberculosis.

Communicated by Ramaswamy H. Sarma     

Combined pharmacophore modeling, 3D-QSAR and docking studies to identify novel HDAC inhibitors using drug repurposing

Abstract

Histone deacetylases (HDACs), a critical family of epigenetic enzymes, has emerged as a promising target for antitumor drugs. Here, we describe our protocol of virtual screening in identification of novel potential HDAC inhibitors through pharmacophore modeling, 3D-QSAR, molecular docking and molecular dynamics (MD) simulation. Considering the limitation of current virtual screening works, drug repurposing strategy was applied to discover druggable HDAC inhibitor. The ligand-based pharmacophore and 3D-QSAR models were established, and their reliability was validated by different methods. Then, the DrugBank database was screened, followed by molecular docking. MD simulation (100?ns) was performed to further study the stability of ligand binding modes. Finally, results indicated the hit DB03889 with high in silico inhibitory potency was suitable for further experimental analysis.

Communicated by Ramaswamy H. Sarma     

Induced fit docking,free energy calculation and molecular dynamics studies on Mycobacterium tuberculosis alanine racemase inhibitor

Alar, a Pyridoxal 5′-phosphate (PLP)-dependent bacterial enzyme is responsible for the racemisation of L-alanine into D-alanine which is essential for the peptidoglycan biosynthesis in both Gram-positive and Gram-negative bacteria. In the present study, we performed induced fit docking, binding free energy calculation and molecular dynamics simulation to elucidate the Alar inhibition potential of 1,2,4-thiadiazolidine-3,5-dione-based inhibitor 1. The inhibitor binds to the hydrophobic groove of Alar and the binding was found to be stable throughout 20-ns MD simulation. Induced fit docking result showed that Lys42, Tyr46, Tyr175 and Tyr364 residues are primarily responsible for the stabilisation of inhibitor–protein complex. Further, high negative van der Waals binding free energy value of –38.88 kcal/mol, indicated it as the main driving force for the inhibitor binding. Based on the information obtained from this study, we designed few molecules as potent Alar inhibitor. In order to gain structural insight and to validate the stability of complex, we performed 20-ns MD simulation of the designed molecule D1. Results obtained from this study can be used for the design of M. tuberculosis Alar potent inhibitors lacking affinity for the co-factor PLP.     

Discovery of new potential triplet acting inhibitor for Alzheimer’s disease via X-ray crystallography,molecular docking and molecular dynamics

Abstract

The most common brain disorder of late life is Alzheimer’s disease (AD), which is highly complicating dementia. There are several drug targets which are reported to control the severe level of AD; notably, acetylcholinesterase, β-Secretase and glycogen synthase kinase enzymes are approached as a good drug targets for AD. Hence, the present study mainly focused to discover newly synthesized molecule (7-propyl-6H-pyrano[3,2-c:5,6-c']dichromene-6,8(7H)-dione) as a potential triplet acting drug for above said enzymes through the analysis of X-ray crystallography, molecular docking, molecular dynamics and quantum chemical calculation. The target drug molecule was crystallized in the monoclinic crystal structure with P21/n space group. The structure was solved by SHELXS and refined by SHELXL. The crystal packing is stabilized by C???H···O type of interactions. Further, the induced fit docking shows that the molecule has high docking score, glide energy, favorable hydrogen bonding and hydrophobic interactions on the protein targets. The molecular dynamics simulation was performed to understand the stability of the molecule in the presence of active site environment. Finally, quantum chemical calculation has been carried out for the molecule in gas phase and for the corresponding molecule lifted from the active site region. The structural comparison between gas phase and active site helps to understand the conformational modification of the molecule in the active site.

Communicated by Ramaswamy H. Sarma     

An explorative study on Staphylococcus aureus MurE inhibitor: induced fit docking,binding free energy calculation,and molecular dynamics

Staphylococcus aureus MurE enzyme catalyzes the addition of l-lysine as third residue of the peptidoglycan peptide moiety. Due to the high substrate specificity and its ubiquitous nature among bacteria, MurE enzyme is considered as one of the potential target for the development of new therapeutic agents. In the present work, induced fit docking (IFD), binding free energy calculation, and molecular dynamics (MD) simulation were carried out to elucidate the inhibition potential of 2-thioxothiazolidin-4-one based inhibitor 1 against S. aureus MurE enzyme. The inhibitor 1 formed majority of hydrogen bonds with the central domain residues Asn151, Thr152, Ser180, Arg187, and Lys219. Binding free-energy calculation by MM-GBSA approach showed that van der Waals (ΔG_vdW, ?57.30?kcal/mol) and electrostatic solvation (ΔG_solv, ?36.86?kcal/mol) energy terms are major contributors for the inhibitor binding. Further, 30-ns MD simulation was performed to validate the stability of ligand–protein complex and also to get structural insight into mode of binding. Based on the IFD and MD simulation results, we designed four new compounds D1–D4 with promising binding affinity for the S. aureus MurE enzyme. The designed compounds were subjected to the extra-precision docking and binding free energy was calculated for complexes. Further, a 30-ns MD simulation was performed for D1/4C13 complex.     

Unravelling the binding affinity between model transport protein and a prospective tuberculosis therapeutic agent: a spectroscopic and theoretical simulation exploration

Abstract

Haloxyfop was reported to exhibit inhibition effect targeting Mycobacterium tuberculosis and pathogenic parasites. To pave its way for drug development, more research is required to determine the affinities interacting with biological receptors in vivo. In this work, the interactions of Haloxyfop with two model transport proteins were investigated by spectroscopic techniques and theoretical simulation. The interaction mechanism, thermodynamic properties and the impact of Haloxyfop-induced conformational change in serum albumins were revealed by series of fluorescence, UV-Vis absorption and circular dichroic spectroscopy. The specific binding sites were determined by site-competitive replacement experiment. Molecular docking and dynamic simulation provided a visual screening in the microscopic binding mode. The structure of Haloxyfop was roughly divided into three parts that exhibit different covalent interaction affinities. The two isomers of Haloxyfop showed a certain degree of affinity difference. Hydrophobic, polar interaction and π-effect were analyzed in detail, and the surface electrostatic potential energy maps were simulated to provide references. The free energy, calculated by the molecular mechanics-generalized born surface area (MM-GBSA) and molecular mechanics-Poisson Boltzmann surface area (MM-PBSA) methods, was decomposed to per-residues, which intuitively revealed relevant contributions in binding process. The role of water existence was explored through molecular dynamic refinement, and the frontier molecular orbital analysis explored the ionic interaction mechanism in electronic level. In general, multiple chemistry method was adopted to fully unravel the properties of Haloxyfop-binding for the sake of rationalizing the applicability as a therapeutic agent.

Communicated by Ramaswamy H. Sarma     

Study on the binding mode of a pyrrolotriazin derivative with JAK2 by docking and MD simulation

The pyrrolotriazin derivative 2-(4-(4-((7-(3-(N-methylmethylsulfonamido)phenyl)pyrrolo [2,1-f][1,2,4]triazin-2-yl)amino)phenyl)piperidin-1-yl)acetamide (PPA) is a potential Janus kinase 2 (JAK2) inhibitor. The binding mode between PPA and JAK2 was investigated by using a combined method of docking, molecular dynamics (MD) simulation and binding free-energy calculation. The docking calculations preliminarily indicated that there were two possible binding modes 1 and 2; MD simulations and binding free-energy calculations identified that binding mode 1 was more stable and favourable, with the lower MM-PBSA binding free energy of ?34.00?±?0.17?kcal/mol. Moreover, some valuable binding information is revealed as follows: the inhibitor PPA is suitably located at the ATP-binding site of JAK2 and the hydrophobic interaction plays an essential role. PPA not only interacts with residues Leu855, Val863, Ala880, Tyr931, Leu932 and Leu983 via hydrophobic interaction but also interacts with Ser936 and Asp994 by hydrogen bonds. These two factors are advantageous for PPA to strongly bind to JAK2. These results help to understand the action mechanisms and designing new compounds with a higher affinity to JAK2.     

Assessment on the binding characteristics of dasatinib,a tyrosine kinase inhibitor to calf thymus DNA: insights from multi-spectroscopic methodologies and molecular docking as well as DFT calculation

Abstract

The binding characteristics of calf thymus DNA (ct-DNA) with dasatinib (DSTN), a tyrosine kinase inhibitor was assessed through multi-spectroscopic methodologies and viscosity measurement combined with molecular docking as well as DFT calculation to understand the binding mechanism, affinity of DSTN onto ct-DNA, effect of DSTN on ct-DNA conformation, and among others. The results confirmed DSTN bound onto ct-DNA, leading to forming the DSTN–ct-DNA complex with the binding constant of 4.82?×?10³ M^?1 at 310?K. DSTN preferentially inserted to the minor groove of ct-DNA with rich A-T region, that was the binding mode of DSTN onto ct-DNA was groove binding. The enthalpic change (ΔH⁰) and entropic change (ΔS⁰) during the binding process of DSTN with ct-DNA were 128.9?kJ mol^?1 and 489.2?J mol^?1 K^?1, respectively, confirming clearly that the association of DSTN with ct-DNA was an endothermic process and the dominative driven-force was hydrophobic interaction. Meanwhile, the results also indicated that there was a certain extent of electrostatic force and hydrogen bonding, but they maybe play an auxiliary role. The CD measurement results confirmed the alteration in the helical configuration of ct-DNA but almost no change in the base stacking after binding DSTN. The results revealed that there was the obvious change in the conformation, the dipole moment, and the atomic charge distribution of DSTN in the B-DNA complexes, compared with free DSTN, to satisfy the conformational adaptation. From the obtained fronitier molecular orbitals of DSTN, it can be inferred that the nature of DSTN alters with the change of the environment around DSTN.

Communicated by Ramaswamy H. Sarma     

Exploring the binding mechanism of saccharin and sodium saccharin to promoter of human p53 gene by theoretical and experimental methods

Abstract

In the past few decades, extensive discussions have been on the impact of artificial sweeteners on the risk of cancer. The present study aimed to evaluate the interaction of saccharin (SA) and sodium saccharin (SSA) with the promoter of the human p53 gene. The binding ability was assessed using the spectroscopic technique, molecular docking and molecular dynamics (MD) simulation methods. Free energy of binding has been calculated using Molecular Mechanics/Poisson–Boltzmann Surface Area (MM/PBSA) method. Fluorescence spectra of mentioned gene with concentration profiles of SA and SSA were obtained in a physiological condition. A gradual increase without any significant spectral shift in the fluorescence intensity of around 350?nm was evident, indicating the presence of an interaction between both compounds and gene. The docking results showed that both compounds were susceptible to bind to 5′-DG56DG57-3′ nucleotide sequence of gene. Furthermore, the MD simulation demonstrated that the binding positions for SA and SSA were 5′-A1T3T4-3′ and 5′-G44T45-3′ sequences of gene, respectively. The binding of these sweeteners to gene made significant conformational changes to the DNA structure. Hydrogen and hydrophobic interactions are the major forces in complexes stability. Through the groove binding mode, the non-interactive DNA-binding nature of SSA and SA has been demonstrated by the results of spectrofluorometric and molecular modeling. This study could provide valuable insight into the binding mechanism of SA and its salt with p53 gene promoter as macromolecule at the molecular level in atomistic details. This work can contribute to the possibility of the potential hazard of carcinogenicity of this sweetener and to design and apply new and safer artificial sweeteners. Abbreviations SA Saccharin

SSA Sodium Saccharin

Pp53g promoter of human p53 gene

MD Molecular dynamics

RMSD Root-mean-square deviation

RMSF Root-mean-square fluctuation

Rg Radius of Gyration

SASA Solvent-Accessible Surface Area

ADI Acceptable daily intake

MM/PBSA Molecular Mechanics/Poisson–Boltzmann Surface Area

Communicated by Ramaswamy H. Sarma     

Evaluating the suitability of RNA intervention mechanism exerted by some flavonoid molecules against dengue virus MTase RNA capping site: a molecular docking,molecular dynamics simulation,and binding free energy study

Abstract

Dengue virus (DENV) is one of the most dangerous mosquito-borne human pathogens known to the mankind. Currently, no vaccines or standard therapy is avaliable to treate DENV infection. This makes the drug development against DENV more significant and challenging. The MTase domain of DENV RNA RdRp NS5 is a promising drug target, because this domain hosts the RNA capping process of DENV RNA to escape from human immune system. In the present study, we have analysed the RNA intervention mechanism exerted by flavoniod molecules against NS5 MTase RNA capping site by using molecular docking, molecular dynamics simulation and the binding free energy calculations. The results from the docking analysis confirmed that the RNA intervention mecanism is exerted by the quercetagetin (QGN) molecule with all necessary intermolecular interactions and high binding affinity. Notably, QGN forms strong hydrogen bonding interactions with Asn18, Leu20 and Ser150 residues and π???π stacking interaction with Phe25 residue. The apo and QGN bound NS5 MTase and QGN-NS5 MTase complex were used for MD simulation. The results of MD simulation reveal that the RMSD and RMSF values of QGN-MTase complex have increased on comparing the apo protein due to the effect of ligand binding. The binding free energy calulation includes prediction of total binding free energy of ligand-protein complex and per-residue free energy decomposition. The QGN binding to NS5 MTase affects it’s native motion, this result is found from Principal component analysis.

Communicated by Ramaswamy H. Sarma     

Designing of benzothiazole derivatives as promising EGFR tyrosine kinase inhibitors: a pharmacoinformatics study

Abstract

Benzothiazole derivatives represent an important class of therapeutic chemical agents and are widely used for interesting biological activities and therapeutic functions including anticancer, antitumor and antimicrobial. In this study, we have performed similarity/substructure-based search of eMolecule database to find out promising benzothiazole derivatives as EGFR tyrosine kinase inhibitors. Several screening criteria that included molecular docking, pharmacokinetics and synthetic accessibility were used on initially derived about 7000 molecules consisting of benzothiazole as major component. Finally, four molecules were found to be promising EGFR tyrosine kinase inhibitors. The best docked pose of each molecule was considered for binding interactions followed by molecular dynamics (MD) and binding energy calculation. Molecular docking clearly showed the final proposed derivatives potential to form a number of binding interactions. MD simulation trajectories undoubtedly indicated that the EGFR protein becomes stable when proposed derivatives bind to the receptor cavity. Strong binding affinity was found for all molecules toward the EGFR which was substantiated by the binding energy calculation using the MM-PBSA approach. Therefore, proposed benzothiazole derivatives may be promising EGFR tyrosine kinase inhibitors for potential application as cancer therapy.

Communicated by Ramaswamy H. Sarma     

A computational study of binding between 3-(4-fluorophenyl)-N-((4-fluorophenyl)sulphonyl)acrylamide and tubulin

3-(4-Fluorophenyl)-N-((4-fluorophenyl)sulphonyl)acrylamide (FFSA) is a potential tubulin polymerisation inhibitor. In this article, a theoretical study of the binding between FFSA and tubulin in colchicine site was carried out by molecular docking, molecular dynamics (MD) simulation and binding free energy calculations. The docking calculations preliminarily indicate that there are three possible binding modes 1, 2 and 3; MD simulations and binding free energy calculations identify that binding mode 2 is the most favourable, with the lowest binding free energy of ? 29.54 kcal/mol. Moreover, our valuable results for the binding are as follows: the inhibitor FFSA is suitably located at the colchicine site of tubulin, where it not only interacts with residues Leu248β, Lys254β, Leu255β, Lys352β, Met259β and Val181a by hydrophilic interaction, but also interacts with Val181α and Thr179α by hydrogen bond interaction. These two factors are both essential for FFSA strongly binding to tubulin. These theoretical results help understanding the action mechanism and designing new compounds with higher affinity to tubulin.     

Structural investigations into the binding mode of novel neolignans Cmp10 and Cmp19 microtubule stabilizers by in silico molecular docking,molecular dynamics,and binding free energy calculations

Microtubule stabilizers provide an important mode of treatment via mitotic cell arrest of cancer cells. Recently, we reported two novel neolignans derivatives Cmp10 and Cmp19 showing anticancer activity and working as microtubule stabilizers at micromolar concentrations. In this study, we have explored the binding site, mode of binding, and stabilization by two novel microtubule stabilizers Cmp10 and Cmp19 using in silico molecular docking, molecular dynamics (MD) simulation, and binding free energy calculations. Molecular docking studies were performed to explore the β-tubulin binding site of Cmp10 and Cmp19. Further, MD simulations were used to probe the β-tubulin stabilization mechanism by Cmp10 and Cmp19. Binding affinity was also compared for Cmp10 and Cmp19 using binding free energy calculations. Our docking results revealed that both the compounds bind at Ptxl binding site in β-tubulin. MD simulation studies showed that Cmp10 and Cmp19 binding stabilizes M-loop (Phe272-Val288) residues of β-tubulin and prevent its dynamics, leading to a better packing between α and β subunits from adjacent tubulin dimers. In addition, His229, Ser280 and Gln281, and Arg278, Thr276, and Ser232 were found to be the key amino acid residues forming H-bonds with Cmp10 and Cmp19, respectively. Consequently, binding free energy calculations indicated that Cmp10 (?113.655 kJ/mol) had better binding compared to Cmp19 (?95.216 kJ/mol). This study provides useful insight for better understanding of the binding mechanism of Cmp10 and Cmp19 and will be helpful in designing novel microtubule stabilizers.     

Molecular Dynamics,Flexible Docking,Virtual Screening,ADMET Predictions,and Molecular Interaction Field Studies to Design Novel Potential MAO-B Inhibitors

Abstract

Monoamine oxidase is a flavoenzyme bound to the mitochondrial outer membranes of the cells, which is responsible for the oxidative deamination of neurotransmitter and dietary amines. It has two distinct isozymic forms, designated MAO-A and MAO-B, each displaying different substrate and inhibitor specificities. They are the well-known targets for antidepressant, Parkinson's disease, and neuroprotective drugs. Elucidation of the x-ray crystallographic structure of MAO-B has opened the way for the molecular modeling studies. In this work we have used molecular modeling, density functional theory with correlation, virtual screening, flexible docking, molecular dynamics, ADMET predictions, and molecular interaction field studies in order to design new molecules with potential higher selectivity and enzymatic inhibitory activity over MAO-B.     

Discovery of potential sclerostin inhibitors from plants with loop2 region of sclerostin inhibition by interacting with residues outside Pro-Asn-Ala-Ile-Gly motif

Abstract

Sclerostin, an antagonist of the Wnt/β-catenin signaling pathway, was discovered as a potential therapeutic target for stimulating bone formation in osteoporosis. In this study, molecular docking was employed to predict the binding of 29 herbal compounds, which were reported as bone formation stimulators, to the loop2 region of sclerostin. Then, the 50 ns molecular dynamics (MD) simulation of the complexes between sclerostin and the top 10 hits obtained from molecular docking were carried out. Root mean square deviations (RMSDs) analysis of MD trajectories pointed out that all ligands-complexes remain stable throughout the duration of MD simulations. In addition, the molecular mechanics/generalized born surface area (MM/GBSA) binding free energy and energy decomposition analyses were determined. The results here suggested that baicalin is the most promising inhibitor of sclerostin. Interestingly, baicalin binds to sclerostin via the hydrophobic interaction with the amino acid residues on loop2 region but outside the Pro-Asn-Ala-Ile-Gly (PNAIG) motif, particularly the Arg-Gly-Lys-Trp-Trp-Arg (RGKWWR) motif. This finding could be a novel strategy for developing new sclerostin inhibitors in the future.

Communicated by Ramaswamy H. Sarma     

Probing the acting mode and advantages of RMC-4550 as an Src-homology 2 domain-containing protein tyrosine phosphatase (SHP2) inhibitor at molecular level through molecular docking and molecular dynamics

Abstract

The over-activation of Ras/mitogen-activated protein kinase (MAPK) signaling pathway associated with a variety of cancers is usually related with abnormal activation of Src-homology 2 domain-containing protein tyrosine phosphatase (SHP2). For this purpose, SHP2 has attracted extensive interest as a potential target for cancer treatment. RMC-4550, as a newly developed selective inhibitor of SHP2, possesses an overwhelming advantage over the previous generation inhibitor SHP099 in terms of in vitro activity. However, the binding mode of SHP2 with RMC-4550 and the reason for the high efficiency of RMC-4550 as SHP2 inhibitor at molecular level are still unclear. Therefore, in this study, the binding mode of RMC-4550 with SHP2 and the superiorities of RMC-4550 as inhibitor at binding affinity and dynamic interactive behavior with SHP2 were probed by molecular docking and molecular dynamics (MD) simulations. By comparing the results of molecular docking, it was found that SHP2 formed more tight interaction with RMC-4550 than that with SHP099. Subsequently, a series of post-dynamic analyses on three simulation trajectories (SHP2^WT, SHP2^SHP099 and SHP2^RMC-4550) were performed and found that the SHP2 protein bound with RMC-4550 maintained a firmer interaction between N-Src-homology 2 (N-SH2) and PTP domain throughout the MD simulation, leading to a more stable protein conformation. The finding here provides new clues for the design of SHP2 inhibitor against the over-activation of Ras/MAPK pathway.

Communicated by Ramaswamy H. Sarma     

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.     

Communication of γ Phage Lysin plyG Enzymes Binding toward SrtA for Inhibition of Bacillus Anthracis: Protein–Protein Interaction and Molecular Dynamics Study

Abstract

Bacillus anthracis is a pathogenic, Gram-positive bacterium which chiefly affects the livestock of animals and humans through acute disease anthrax. All around the globe this bio-threat organism damages millions of lives in every year and also most of the drugs were not responding properly in inhibition against this diseased pathogen. In recent development, phage therapy is considered as alternative solution to treat this serious infectious disease. In this study, we elucidated the binding of γ phage lysin plyG enzymes toward the SrtA along with its activator peptide LPXTG. Through protein–protein docking and molecular dynamics simulation studies, we showed the distinguished structure complementarity of SrtA and plyG complex. Especially, MD simulation relates strong and stable interaction occurs between the protein complex structures. These results suggest that additional experimental studies on our approach will lead to availability of better inhibitor against the SrtA.     

The identification of novel,high affinity AQP9 inhibitors in an intracellular binding site

Abstract

Background: The involvement of aquaporin (AQP) water and small solute channels in the etiology of several diseases, including cancer, neuromyelitis optica and body fluid imbalance disorders, has been suggested previously. Furthermore, results obtained in a mouse model suggested that AQP9 function contributes to hyperglycemia in type-2 diabetes. In addition, the physiological role of several AQP family members remains poorly understood. Small molecule inhibitors of AQPs are therefore desirable to further study AQP physiological and pathophysiological functions. Methods: The binding of recently established AQP9 inhibitors to a homology model of AQP9 was investigated by molecular dynamics simulations and molecular docking. Putative inhibitor binding sites identified with this procedure were modified by site-directed mutagenesis. Active compounds were measured in a mammalian cell water permeability assay of mutated AQP9 isoforms and tested for changes in inhibitory effects. Controls: Three independent cell lines were established for each mutated AQP9 isoform and functionality of mutant isoforms was established. Principal findings: We have identified putative binding sites of recently established AQP9 inhibitors. This information facilitated successful identification of novel AQP9 inhibitors with low micromolar IC50 values in a cell based assay by in silico screening of a compound library targeting specifically this binding site. Significance: We have established a successful strategy for AQP small molecule inhibitor identification. AQP inhibitors may be relevant as experimental tools, to enhance our understanding of AQP function, and in the treatment of various diseases.     

An integrated in silico approach to understand protein–protein interactions: human meprin-β with fetuin-A

Abstract

Human meprin-β, a zinc metalloprotease belonging to the astacin family, have been found to be associated with many pathological conditions like inflammatory bowel disease, fibrosis and neurodegenerative disease. The inhibition of meprin-β by various inhibitors, both macromolecular and small molecules, is crucial in the control of several diseases. Human fetuin-A, a negative acute phase protein involved in inflammatory disease, has recently been identified as an endogenous inhibitor for meprin-β. In this computational study, an integrated in silico approach was performed using existing structural information of meprin-β coupled with ab initio modelling of human fetuin-A to predict a rational model of the complex through protein–protein docking. Further, the models were optimized and validated to generate an ensemble of conformations through extensive molecular dynamics simulation. Virtual alanine scanning mutagenesis was explored to identify hotspot residues on both proteins significant for protein–protein interaction (PPI). The results of the study provide structural insight into PPI between meprin-β and fetuin-A which can be useful in designing molecules to modulate meprin-β activity.

Communicated by Ramaswamy H. Sarma