Computational characterization of how the VX nerve agent binds human serum paraoxonase 1

Human serum paraoxonase 1 (HuPON1) is an enzyme that can hydrolyze various chemical warfare nerve agents including VX. A previous study has suggested that increasing HuPON1’s VX hydrolysis activity one to two orders of magnitude would make the enzyme an effective countermeasure for in vivo use against VX. This study helps facilitate further engineering of HuPON1 for enhanced VX-hydrolase activity by computationally characterizing HuPON1’s tertiary structure and how HuPON1 binds VX. HuPON1’s structure is first predicted through two homology modeling procedures. Docking is then performed using four separate methods, and the stability of each bound conformation is analyzed through molecular dynamics and solvated interaction energy calculations. The results show that VX’s lone oxygen atom has a strong preference for forming a direct electrostatic interaction with HuPON1’s active site calcium ion. Various HuPON1 residues are also detected that are in close proximity to VX and are therefore potential targets for future mutagenesis studies. These include E53, H115, N168, F222, N224, L240, D269, I291, F292, and V346. Additionally, D183 was found to have a predicted pKa near physiological pH. Given D183’s location in HuPON1’s active site, this residue could potentially act as a proton donor or accepter during hydrolysis. The results from the binding simulations also indicate that steered molecular dynamics can potentially be used to obtain accurate binding predictions even when starting with a closed conformation of a protein’s binding or active site.     

Functional and structural insights revealed by molecular dynamics simulations of an essential RNA editing ligase in Trypanosoma brucei

RNA editing ligase 1 (TbREL1) is required for the survival of both the insect and bloodstream forms of Trypanosoma brucei, the parasite responsible for the devastating tropical disease African sleeping sickness. The type of RNA editing that TbREL1 is involved in is unique to the trypanosomes, and no close human homolog is known to exist. In addition, the high-resolution crystal structure revealed several unique features of the active site, making this enzyme a promising target for structure-based drug design. In this work, two 20 ns atomistic molecular dynamics (MD) simulations are employed to investigate the dynamics of TbREL1, both with and without the ATP substrate present. The flexibility of the active site, dynamics of conserved residues and crystallized water molecules, and the interactions between TbREL1 and the ATP substrate are investigated and discussed in the context of TbREL1's function. Differences in local and global motion upon ATP binding suggest that two peripheral loops, unique to the trypanosomes, may be involved in interdomain signaling events. Notably, a significant structural rearrangement of the enzyme's active site occurs during the apo simulations, opening an additional cavity adjacent to the ATP binding site that could be exploited in the development of effective inhibitors directed against this protozoan parasite. Finally, ensemble averaged electrostatics calculations over the MD simulations reveal a novel putative RNA binding site, a discovery that has previously eluded scientists. Ultimately, we use the insights gained through the MD simulations to make several predictions and recommendations, which we anticipate will help direct future experimental studies and structure-based drug discovery efforts against this vital enzyme.     

Analysis of active-site amino-acid residues of human serum paraoxonase using competitive substrates

Serum paraoxonase (PON1) is a calcium-dependent six-fold beta-propeller protein structurally similar to the di-isopropylfluorophosphatase (DFPase) found in the squid Loligo vulgaris. Human serum paraoxonase (HuPON1) has been shown to hydrolyze an array of substrates even though relatively little is known about its physiological role(s) or its catalytic mechanism. Through site-directed mutagenesis studies, designed from a DFPase-like homology model, and from a crystal structure of a hybrid PON1 molecule, amino-acid residues essential for enzyme function, including H115 and F222, have been identified. It was shown previously that, when H115 is replaced with tryptophan, the resulting enzyme hydrolyzes paraoxon but not phenyl acetate. This study shows that, when present simultaneously, phenyl acetate competitively inhibits paraoxon hydrolysis by H115W. Conversely, when F222 is replaced with tyrosine, mutant F222Y can hydrolyze phenyl acetate but not paraoxon. The presence of DFP, an inhibitor of both arylesterase and paraoxonase activities of wild-type HuPON1 (mean Ki=0.48+/-0.15 mM), has no effect on the ability of F222Y to catalyze the hydrolysis of phenyl acetate, suggesting that the F222Y mutant is unable to bind DFP. Together, the results suggest that, in wild-type HuPON1, H115 and F222 are important in determining substrate binding and specificity, but are not likely to be directly involved in substrate hydrolysis.     

Identification of Crucial Amino Acids in Mouse Aldehyde Oxidase 3 That Determine Substrate Specificity

In order to elucidate factors that determine substrate specificity and activity of mammalian molybdo-flavoproteins we performed site directed mutagenesis of mouse aldehyde oxidase 3 (mAOX3). The sequence alignment of different aldehyde oxidase (AOX) isoforms identified variations in the active site of mAOX3 in comparison to other AOX proteins and xanthine oxidoreductases (XOR). Based on the structural alignment of mAOX3 and bovine XOR, differences in amino acid residues involved in substrate binding in XORs in comparison to AOXs were identified. We exchanged several residues in the active site to the ones found in other AOX homologues in mouse or to residues present in bovine XOR in order to examine their influence on substrate selectivity and catalytic activity. Additionally we analyzed the influence of the [2Fe-2S] domains of mAOX3 on its kinetic properties and cofactor saturation. We applied UV-VIS and EPR monitored redox-titrations to determine the redox potentials of wild type mAOX3 and mAOX3 variants containing the iron-sulfur centers of mAOX1. In addition, a combination of molecular docking and molecular dynamic simulations (MD) was used to investigate factors that modulate the substrate specificity and activity of wild type and AOX variants. The successful conversion of an AOX enzyme to an XOR enzyme was achieved exchanging eight residues in the active site of mAOX3. It was observed that the absence of the K889H exchange substantially decreased the activity of the enzyme towards all substrates analyzed, revealing that this residue has an important role in catalysis.     

Structure/function analyses of human serum paraoxonase (HuPON1) mutants designed from a DFPase-like homology model

Human serum paraoxonase (HuPON1) is a calcium-dependent enzyme that hydrolyzes esters, including organophosphates and lactones, and exhibits anti-atherogenic properties. A few amino acids have been shown to be essential for the enzyme's arylesterase and organophosphatase activities. Until very recently, a three-dimensional model was not available for HuPON1, so functional roles have not been assigned to those residues. Based on sequence-structure alignment studies, we have folded the amino acid sequence of HuPON1 onto the sixfold beta-propeller structure of squid diisopropylfluorophosphatase (DFPase). We tested the validity of this homology model by circular dichroism (CD) spectroscopy and site-directed mutagenesis. Consistent with predictions from the homology model, CD data indicated that the structural composition of purified HuPON1 consists mainly of beta-sheets. Mutants of HuPON1 were assayed for enzymatic activity against phenyl acetate and paraoxon. Substitution of residues predicted to be important for substrate binding (L69, H134, F222, and C284), calcium ion coordination (D54, N168, N224, and D269), and catalytic mechanism of HuPON1 (H285) led to enzyme inactivation. Mutants F222Y and H115W exhibited substrate-binding selectivity towards phenyl acetate and paraoxon, respectively. The homology model presented here is very similar to the recently obtained PON1 crystal structure, and has allowed identification of several residues within the enzyme active site.     

Kinetic and crystallographic analysis of mutant Escherichia coli aminopeptidase P: insights into substrate recognition and the mechanism of catalysis

Aminopeptidase P (APPro) is a manganese-dependent enzyme that cleaves the N-terminal amino acid from polypeptides where the second residue is proline. APPro shares a similar fold, substrate specificity, and catalytic mechanism with methionine aminopeptidase and prolidase. To investigate the roles of conserved residues at the active site, seven mutant forms of APPro were characterized kinetically and structurally. Mutation of individual metal ligands selectively abolished binding of either or both Mn(II) atoms at the active site, and none of these metal-ligand mutants had detectable catalytic activity. Mutation of the conserved active site residues His243 and His361 revealed that both are required for catalysis. We propose that His243 stabilizes substrate binding through an interaction with the carbonyl oxygen of the requisite proline residue of a substrate and that His361 stabilizes substrate binding and the gem-diol catalytic intermediate. Sequence, structural, and kinetic analyses reveal that His350, conserved in APPro and prolidase but not in methionine aminopeptidase, forms part of a hydrophobic binding pocket that gives APPro its proline specificity. Further, peptides in which the required proline residue is replaced by N-methylalanine or alanine are cleaved by APPro, but they are extremely poor substrates due to a loss of interactions between the prolidyl ring of the substrate and the hydrophobic proline-binding pocket.     

Structural basis for broad substrate specificity in higher plant beta-D-glucan glucohydrolases

Family 3 beta-D-glucan glucohydrolases are distributed widely in higher plants. The enzymes catalyze the hydrolytic removal of beta-D-glucosyl residues from nonreducing termini of a range of beta-D-glucans and beta-D-oligoglucosides. Their broad specificity can be explained by x-ray crystallographic data obtained from a barley beta-D-glucan glucohydrolase in complex with nonhydrolyzable S-glycoside substrate analogs and by molecular modeling of enzyme/substrate complexes. The glucosyl residue that occupies binding subsite -1 is locked tightly into a fixed position through extensive hydrogen bonding with six amino acid residues near the bottom of an active site pocket. In contrast, the glucosyl residue at subsite +1 is located between two Trp residues at the entrance of the pocket, where it is constrained less tightly. The relative flexibility of binding at subsite +1, coupled with the projection of the remainder of bound substrate away from the enzyme's surface, means that the overall active site can accommodate a range of substrates with variable spatial dispositions of adjacent beta-D-glucosyl residues. The broad specificity for glycosidic linkage type enables the enzyme to perform diverse functions during plant development.     

The mechanism of inhibition of the cyclin-dependent kinase-2 as revealed by the molecular dynamics study on the complex CDK2 with the peptide substrate HHASPRK

Molecular dynamics (MD) simulations were used to explain structural details of cyclin-dependent kinase-2 (CDK2) inhibition by phosphorylation at T14 and/or Y15 located in the glycine-rich loop (G-loop). Ten-nanosecond-long simulations of fully active CDK2 in a complex with a short peptide (HHASPRK) substrate and of CDK2 inhibited by phosphorylation of T14 and/or Y15 were produced. The inhibitory phosphorylations at T14 and/or Y15 show namely an ATP misalignment and a G-loop shift (~5 A) causing the opening of the substrate binding box. The biological functions of the G-loop and GxGxxG motif evolutionary conservation in protein kinases are discussed. The position of the ATP gamma-phosphate relative to the phosphorylation site (S/T) of the peptide substrate in the active CDK2 is described and compared with inhibited forms of CDK2. The MD results clearly provide an explanation previously not known as to why a basic residue (R/K) is preferred at the P(2) position in phosphorylated S/T peptide substrates.     

Homology modeling of human serum paraoxonase1 and its molecular interaction studies with aspirin and cefazolin

Human serum paraoxonase1 (HuPON1) belongs to the family of A-esterases (EC.3.1.8.1). It is associated with HDL particle and prevents atherosclerosis by cleaving lipid hydroperoxides and other proatherogenic molecules of oxidized low density lipoproteins (LDL). Since the precise structure of HuPON1 is not yet available, the structure-function relationship between HuPON1 and activators/inhibitors is still unknown. Therefore, a theoretical model of HuPON1 was generated using homology modelling and precise molecular interactions of an activator aspirin and an inhibitor cefazolin with PON1 were studied using Autodock software. The ligand binding residues were found to be similar to the predicted active site residues. Both cefazolin and aspirin were found to dock in the vicinity of the predicted active sites of PON1; cefazolin bound at residues N166, S193 and Y71, while aspirin at residues N309, I310 and L311. Binding region in the PON1 by prediction (3D2GO server) and docking studies provide useful insight into mechanism of substrate and inhibitor binding to the enzyme active site.     

Stereoselectivity of phosphotriesterase with paraoxon derivatives: a computational study

The bacterial enzyme phosphotriesterase (PTE) exhibits stereoselectivity toward hydrolysis of chiral substrates with a preference for the S_p enantiomer. In this work, docking analysis and two explicit-solvent molecular dynamics (MD) simulations were performed to characterize and differentiate the structural dynamics of PTE bound to the S_p and R_p paraoxon derivative enantiomers (R_p-1 and S_p-1) hydrolyzed with distinct catalytic efficiencies. Comparative analysis of the molecular trajectories for PTE bound to R_p-1 and S_p-1 suggested that substrate binding induced conformational changes in the loops near the active site. After 100 ns of MD simulation, the Zn β²⁺ metal ion formed hexacoordinated- and tetracoordinated geometries in the S_p-1-PTE and R_p-1-PTE ensembles, respectively. Simulation results further showed that the hydrogen bond between Asp301 and His254 occurred with a higher probability after S_p-1 binding to PTE (47.5%) than that after R_p-1 binding (22.2%). These results provide a qualitative and molecular-level explanation for the 10 orders of magnitude increase in the catalytic efficiency of PTE toward the S_p enantiomer of paraoxon.     

Dynamic profiling of double-stranded RNA binding proteins

Double-stranded (ds) RNA is a key player in numerous biological activities in cells, including RNA interference, anti-viral immunity and mRNA transport. The class of proteins responsible for recognizing dsRNA is termed double-stranded RNA binding proteins (dsRBP). However, little is known about the molecular mechanisms underlying the interaction between dsRBPs and dsRNA. Here we examined four human dsRBPs, ADAD2, TRBP, Staufen 1 and ADAR1 on six dsRNA substrates that vary in length and secondary structure. We combined single molecule pull-down (SiMPull), single molecule protein-induced fluorescence enhancement (smPIFE) and molecular dynamics (MD) simulations to investigate the dsRNA-dsRBP interactions. Our results demonstrate that despite the highly conserved dsRNA binding domains, the dsRBPs exhibit diverse substrate specificities and dynamic properties when in contact with different RNA substrates. While TRBP and ADAR1 have a preference for binding simple duplex RNA, ADAD2 and Staufen1 display higher affinity to highly structured RNA substrates. Upon interaction with RNA substrates, TRBP and Staufen1 exhibit dynamic sliding whereas two deaminases ADAR1 and ADAD2 mostly remain immobile when bound. MD simulations provide a detailed atomic interaction map that is largely consistent with the affinity differences observed experimentally. Collectively, our study highlights the diverse nature of substrate specificity and mobility exhibited by dsRBPs that may be critical for their cellular function.     

Oligomeric states of the detergent-solubilized human serum paraoxonase (PON1)

Human plasma paraoxonase (HuPON1) is a high density lipoprotein (HDL)-bound enzyme exhibiting antiatherogenic properties. The molecular basis for the binding specificity of HuPON1 to HDL has not been established. Isolation of HuPON1 from HDL requires the use of detergents. We have determined the activity, dispersity, and oligomeric states of HuPON1 in solutions containing mild detergents using nondenaturing electrophoresis, size exclusion chromatography, and cross-linking. HuPON1 was active whatever its oligomeric state. In nonmicellar solutions, HuPON1 was polydisperse. In contrast, HuPON1 exhibited apparent homogeneity in micellar solutions, except with CHAPS. The enzyme apparent hydrodynamic radius varied with the type of detergent and protein concentration. In C(12)E(8) micellar solutions, from sedimentation velocity, equilibrium analytical ultracentrifugation, and radioactive detergent binding, HuPON1 was described as monomers and dimers in equilibrium. A decrease of the detergent concentration shifted this equilibrium toward the formation of dimers. About 100 detergent molecules were associated per monomer and dimer. The assembly of amphiphilic molecules, phospholipids in vivo, in sufficiently large aggregates could be a prerequisite for anchoring of HuPON1 and then allowing stabilization of the enzyme activity. Changes of HDL size and shape could strongly affect the binding affinity and stability of HuPON1 and result in reduced antioxidative capacity of the lipoprotein.     

Regulatory phosphorylation of cyclin-dependent kinase 2: insights from molecular dynamics simulations

The structures of fully active cyclin-dependent kinase-2 (CDK2) complexed with ATP and peptide substrate, CDK2 after the catalytic reaction, and CDK2 inhibited by phosphorylation at Thr14/Tyr15 were studied using molecular dynamics (MD) simulations. The structural details of the CDK2 catalytic site and CDK2 substrate binding box were described. Comparison of MD simulations of inhibited complexes of CDK2 was used to help understand the role of inhibitory phosphorylation at Thr14/Tyr15. Phosphorylation at Thr14/Tyr15 causes ATP misalignment for the phosphate-group transfer, changes in the Mg²⁺ coordination sphere, and changes in the H-bond network formed by CDK2 catalytic residues (Asp127, Lys129, Asn132). The inhibitory phosphorylation causes the G-loop to shift from the ATP binding site, which leads to opening of the CDK2 substrate binding box, thus probably weakening substrate binding. All these effects explain the decrease in kinase activity observed after inhibitory phosphorylation at Thr14/Tyr15 in the G-loop. Interaction of the peptide substrate, and the phosphorylated peptide product, with CDK2 was also studied and compared. These results broaden hypotheses drawn from our previous MD studies as to why a basic residue (Arg/Lys) is preferred at the P₊₂ substrate position. Figure View of the substrate binding site of the fully active cyclin-dependent kinase-2 (CDK2) (pT160-CDK2/cyclin A/ATP). The pThr160 activation site is located in the T-loop (yellow secondary structure). The G-loop, which partly forms the ATP binding site, is shown in blue. The Thr14 and Tyr15 inhibitory phosphorylation sites located in the G-loop are shown in licorice representation     

Docking and molecular dynamics simulation studies on glycation-induced conformational changes of human paraoxonase 1

Human paraoxonase 1 (huPON1) is a calcium-dependent esterase responsible for hydrolysis of a wide variety of substrates including organophosphates, esters, lactones, and paraoxon. Although its natural substrate is unknown, the action of PON as an antioxidant is well documented. Because recent reports have suggested glycation may induce reduced PON activity in diabetes, we investigated the structural features of huPON1 and its glycated mutant by template-based modeling, docking, and molecular dynamics (MD) simulations. Our results corroborated the importance of the His115–His134 dyad in both the lactonase and paraoxonase activity of huPON1. Structural alterations in the glycated model reflected weak interactions between the docked substrate and the active site cleft. We also used MD simulation to gain insight into glycation-induced conformational changes of huPON1 and the implication of this on depleted enzymatic activity. The catalytic calcium found on the surface interacts with the side chain oxygen of residues, including Asn224, Asn270, Asn168, Asp269, and Glu53, and this interaction with the respective residues undergoes minor displacement on glycation. The root-mean-square fluctuation had high motional flexibility in the non-glycated model whereas the conformation of the glycated structure was comparatively stable. Our findings emphasize the consequence of glycation-induced alterations and their effect on overall enzymatic activity.     

MD simulations reveal alternate conformations of the oxyanion hole in the Zika virus NS2B/NS3 protease

Recent crystallography studies have shown that the binding site oxyanion hole plays an important role in inhibitor binding, but can exist in two conformations (active/inactive). We have undertaken molecular dynamics (MD) calculations to better understand oxyanion hole dynamics and thermodynamics. We find that the Zika virus (ZIKV) NS2B/NS3 protease maintains a stable closed conformation over multiple 100-ns conventional MD simulations in both the presence and absence of inhibitors. The S1, S2, and S3 pockets are stable as well. However, in two of eight simulations, the A132-G133 peptide bond in the binding pocket of S1' spontaneously flips to form a 3₁₀-helix that corresponds to the inactive conformation of the oxyanion hole, and then maintains this conformation until the end of the 100-ns conventional MD simulations without inversion of the flip. This conformational change affects the S1' pocket in ZIKV NS2B/NS3 protease active site, which is important for small molecule binding. The simulation results provide evidence at the atomic level that the inactive conformation of the oxyanion hole is more favored energetically when no specific interactions are formed between substrate/inhibitor and oxyanion hole residues. Interestingly, however, transition between the active and inactive conformation of the oxyanion hole can be observed by boosting the valley potential in accelerated MD simulations. This supports a proposed induced-fit mechanism of ZIKV NS2B/NS3 protease from computational methods and provides useful direction to enhance inhibitor binding predictions in structure-based drug design.     

Molecular dynamics simulations of lignin peroxidase in solution

The dynamical and structural properties of lignin peroxidase and its Trp171Ala mutant have been investigated in aqueous solution using molecular dynamics (MD) simulations. In both cases, the enzyme retained its overall backbone structure and all its noncovalent interactions in the course of the MD simulations. Very interestingly, the analysis of the MD trajectories showed the presence of large fluctuations in correspondence of the residues forming the heme access channel; these movements enlarge the opening and facilitate the access of substrates to the enzyme active site. Moreover, steered molecular dynamics docking simulations have shown that lignin peroxidase natural substrate (veratryl alcohol) can easily approach the heme edge through the access channel.     

Structural basis for a change in substrate specificity: crystal structure of S113E isocitrate dehydrogenase in a complex with isopropylmalate, Mg2+, and NADP

Isocitrate dehydrogenase (IDH) catalyzes the oxidative decarboxylation of isocitrate and has negligible activity toward other (R)-malate-type substrates. The S113E mutant of IDH significantly improves its ability to utilize isopropylmalate as a substrate and switches the substrate specificity (k(cat)/K(M)) from isocitrate to isopropylmalate. To understand the structural basis for this switch in substrate specificity, we have determined the crystal structure of IDH S113E in a complex with isopropylmalate, NADP, and Mg(2+) to 2.0 A resolution. On the basis of a comparison with previously determined structures, we identify distinct changes caused by the amino acid substitution and by the binding of substrates. The S113E complex exhibits alterations in global and active site conformations compared with other IDH structures that include loop and helix conformational changes near the active site. In addition, the angle of the hinge that relates the two domains was altered in this structure, which suggests that the S113E substitution and the binding of substrates act together to promote catalysis of isopropylmalate. Ligand binding results in reorientation of the active site helix that contains residues 113 through 116. E113 exhibits new interactions, including van der Waals contacts with the isopropyl group of isopropylmalate and a hydrogen bond with N115, which in turn forms a hydrogen bond with NADP. In addition, the loop and helix regions that bind NADP are altered, as is the loop that connects the NADP binding region to the active site helix, changing the relationship between substrates and enzyme. In combination, these interactions appear to provide the basis for the switch in substrate specificity.     

Insights into the binding specificity of wild type and mutated wheat germ agglutinin towards Neu5Acα(2‐3)Gal: a study by in silico mutations and molecular dynamics simulations

Wheat germ agglutinin (WGA) is a plant lectin, which specifically recognizes the sugars NeuNAc and GlcNAc. Mutated WGA with enhanced binding specificity can be used as biomarkers for cancer. In silico mutations are performed at the active site of WGA to enhance the binding specificity towards sialylglycans, and molecular dynamics simulations of 20 ns are carried out for wild type and mutated WGAs (WGA1, WGA2, and WGA3) in complex with sialylgalactose to examine the change in binding specificity. MD simulations reveal the change in binding specificity of wild type and mutated WGAs towards sialylgalactose and bound conformational flexibility of sialylgalactose. The mutated polar amino acid residues Asn114 (S114N), Lys118 (G118K), and Arg118 (G118R) make direct and water mediated hydrogen bonds and hydrophobic interactions with sialylgalactose. An analysis of possible hydrogen bonds, hydrophobic interactions, total pair wise interaction energy between active site residues and sialylgalactose and MM‐PBSA free energy calculation reveals the plausible binding modes and the role of water in stabilizing different binding modes. An interesting observation is that the binding specificity of mutated WGAs (cyborg lectin) towards sialylgalactose is found to be higher in double point mutation (WGA3). One of the substituted residues Arg118 plays a crucial role in sugar binding. Based on the interactions and energy calculations, it is concluded that the order of binding specificity of WGAs towards sialylgalactose is WGA3 > WGA1 > WGA2 > WGA. On comparing with the wild type, double point mutated WGA (WGA3) exhibits increased specificity towards sialylgalactose, and thus, it can be effectively used in targeted drug delivery and as biological cell marker in cancer therapeutics. Copyright © 2014 John Wiley & Sons, Ltd.     

Comparative molecular dynamics simulations of HIV-1 integrase and the T66I/M154I mutant: binding modes and drug resistance to a diketo acid inhibitor

HIV-1 IN is an essential enzyme for viral replication and an interesting target for the design of new pharmaceuticals for use in multidrug therapy of AIDS. L-731,988 is one of the most active molecules of the class of beta-diketo acids. Individual and combined mutations of HIV-1 IN at residues T66, S153, and M154 confer important degrees of resistance to one or more inhibitors belonging to this class. In an effort to understand the molecular mechanism of the resistance of T66I/M154I IN to the inhibitor L-731,988 and its specific binding modes, we have carried out docking studies, explicit solvent MD simulations, and binding free energy calculations. The inhibitor was docked against different protein conformations chosen from prior MD trajectories, resulting in 2 major orientations within the active site. MD simulations have been carried out for the T66I/M154I DM IN, DM IN in complex with L-731,988 in 2 different orientations, and 1QS4 IN in complex with L-731,988. The results of these simulations show a similar dynamical behavior between T66I/M154I IN alone and in complex with L-731,988, while significant differences are observed in the mobility of the IN catalytic loop (residues 138-149). Water molecules bridging the inhibitor to residues from the active site have been identified, and residue Gln62 has been found to play an important role in the interactions between the inhibitor and the protein. This work provides information about the binding modes of L-731,988, as well as insight into the mechanism of inhibitor-resistance in HIV-1 integrase.     

Structural insights into the binding mode of flavonols with the active site of matrix metalloproteinase-9 through molecular docking and molecular dynamic simulations studies

Cartilage degradation in rheumatoid arthritis is mediated principally by the collagenases and gelatinases. Gelatinase B (also called matrix metalloproteinase 9 – MMP-9), is a valid target molecule which is known to participate in cartilage degradation as well as angiogenesis associated with the disease and inhibition of its activity shall prevent cartilage damage and angiogenesis. The focus of this study is to investigate the possibilities of MMP-9 inhibition by flavonol class of bioflavonoids by studying their crucial binding interactions at the active site of MMP 9 using molecular docking (Glide XP and QPLD) and further improvisation by post-docking MM-GBSA and molecular dynamic (MD) simulations. The results show that flavonols can convincingly bind to active site of MMP-9 as demonstrated by their stable interactions at the S1′ specificity pocket and favourable binding energies. Gossypin has emerged as a promising candidate with a docking score of ?14.618 kcal/mol, binding energy of ?79.97 kcal/mol and a stable MD pattern over 15 ns. In addition, interaction mechanisms with respect to catalytic site zinc are also discussed. Further, the drug-like characters of the ligands were also analysed using ADME analysis.