Wei-Wei Han

The metabolism of CYP2C9 and CYP2C19 for gliclazide by homology modeling and docking study.

With homology modeling techniques, a 3D structure model of CYP2C19 was built and refined with molecular mechanics and molecular dynamics simulations. The refined model was assessed to be reasonable by Profile-3D and PROCHECK programs. With the aid of the automatic molecular docking, one substrate and two inhibitors were docked to CYP2C19 by InsightII/Affinity program. The docking results, which are in well agreement with the reported results, demonstrate that the refined model of CYP2C19 is reliable. Then, with the refined model of CYP2C19 and the crystal structure of CYP2C9, the metabolisms of them for gliclazide in two different metabolic pathways were studied and the results show that both enzymes have more favorable interaction energies and stronger affinity with gliclazide in methylhydroxylation pathway than in 6beta-hydroxylation pathway. It is exciting that substrate inhibition phenomenon can be found in metabolisms of CYP2C9 and CYP2C19 for gliclazide in two metabolic pathways. Gliclazide can change the conformation of the active sites and decrease obviously the affinities between gliclazide in the active site and enzymes when it is docked in the second active sites in CYP2C9 and CYP2C19. These results are in well agreement with the kinetic experimental results.

Experimental and computational studies indicate the mutation of Glu12 to increase the thermostability of oligomeric protease from Pyrococcus horikoshii

AbstractThe intracellular protease from Pyrococcus horikoshii (PhpI) is a member of the DJ-1/ThiJ/PfpI superfamily, which is suggested to be involved in cellular protection against environmental stresses. In this study, flexible docking approach was employed to dock the ligand into the active site of PhpI. By analyzing the results, active site architecture and certain key residues responsible for substrate specificity were identified on the enzyme. Our docking result indicates that Glu12 plays an important role in substrate binding. The kinetic experiment conducted by Zhan shows that the E12T mutant is more stable than that of the wild-type. We also predict that Glu15, Lys43, and Tyr46 may be important in the catalytic efficiency and thermostability of enzyme. The new structural and mechanistic insights obtained from computational study should be valuable for detailed structures and mechanisms of the member of the DJ-1 superfamily. Figurea AAFR-AMC in the a E12T b WT

Computational study on new natural polycyclic compounds of H1N1 influenza virus neuraminidase

A new strain of influenza A (H1N1) virus is a major cause of morbidity and mortality around the world. The neuraminidase of the influenza virus has been the most potential target for the anti-influenza drugs such as oseltamivir and zanamivir. However, the emergence of drug-resistant variants of these drugs makes a pressing need for the development of new neuraminidase inhibitors for controlling illness and transmission. Here a 3D structure model of H1N1 avian influenza virus neuraminidase type 1 (N1) was constructed based on the structure of the template H5N1 avian influenza virus N1. Upon application of virtual screening technique for N1 inhibitors, two novel compounds (ZINC database ID: ZINC02128091, ZINC02098378) were found as the most favorable interaction energy with N1. Docking results showed that the compounds bound not only in the active pocket, but also in a new hydrophobic cave which contains Arg368, Trp399, Ile427, Pro431 and Lys432 of N1. Our result suggested that both of the screened compounds containing the hydrophobic group bring a strong conjugation effect with Arg293, Arg368 Lys432 of N1 by pi-pi interaction. However, the control inhibitors zanamivir and oseltamivir do not have this effect. The details of N1-compound binding structure obtained will be valuable for the development of a new anti-influenza virus agent.

CATALYTIC REACTION MECHANISM OF HUMAN PHOTORECEPTOR RETINOL DEHYDROGENASE: A THEORETICAL STUDY

Human photoreceptor retinol dehydrogenase (hRDH8) catalyzes the reduction of all-trans-retinal to all-trans-retinol with NADPH as a rate-limiting step in the visual cycle. Based on the docking results of the substrate to the 3D structure of hRDH8 which is generated by homology modeling method, three quantum chemical calculation models with different sizes were used to investigate the catalytic reaction mechanism of hRDH8 with the aid of density functional theory. The calculations indicate that hRDH8 employs a general acid/base mechanism that a proton is transferred to the keto oxygen of the substrate after the pro-S hydride of NADPH transfer to keto carbon of the substrate. The H-transfer order is converse to that in the proposed mechanism of 17s-hydroxysteroid dehydrogenase 1, which is highly related to the hRDH8 sequence. Tyr155 always provides the proton to the keto oxygen of the substrate whether unprotonated Lys159 is considered or not in the calculation models. However, protonated Lys159 changes the...

Molecular docking study of the affinity of CYP2C9 and CYP2D6 for imrecoxib

With the aid of the automatic molecular docking, the affinity of CYP2C9 and CYP2D6 for imrecoxib was studied by InsightII/Affinity program. The results indicate that CYP2C9–imrecoxib complex has higher stability and stronger affinity because CYP2C9 has more favorable interaction energy (-62.72 kcal/mol) and higher Ludi score (610) with imrecoxib than CYP2D6 (-50.22 kcal/mol and 551) and this is consistent with the results of the kinetic experiments by Li et al. By analyzing the theoretical results combined with the experimental ones, we suggest that the affinity difference is caused by the difference of the structure between CYP2C9 and CYP2D6, and the most important residues for enzyme–substrate complexes, such as Phe476, Asn204, Phe100, Leu366 and Arg108 of CYP2C9 and Phe120, Glu216, and Phe483 of CYP2D6 were also identified.

Characterization of the PH1704 protease from Pyrococcus horikoshii OT3 and the critical functions of Tyr120.

The PH1704 protease from hyperthermophilic archaean Pyrococcus horikoshii OT3 is a member of DJ-1/ThiJ/PfpI superfamily with diverse functional subclasses. The recombinant PH1704 was efficiently purified and was systematically characterized by a combination of substrate specificity analysis, steady-state kinetics study and molecular docking research. The homogeneous protease was obtained as a presumed dodecamer with molecular weight of ∼240 kDa. Iodoacetamide strongly inhibited the peptidase activity, confirming that Cys100 is a nucleophilic residue. The recombinant protein was identified as both an aminopeptidase and an endopeptidase. Experimental data showed that L-R-amc was the best substrate of PH1704. Structural interaction fingerprint analysis (SIFt) indicated the binding pose of PH1704 and showed that Tyr120 is important in substrate binding. Kinetic parameters K cat and K cat /K m of the Y120P mutant with L-R-amc was about 7 and 7.8 times higher than that of the wild type (WT). For the endopeptidase Y120P with AAFR-amc, K cat and K cat /K m is 10- and 21- fold higher than that of WT. Experimental data indicate the important functions of Tyr120: involvement in enzyme activity to form a hydrogen bond with Cys100 and as an entrance gate of the substrate with Lys43. The results of this study can be used to investigate the DJ-1/ThiJ/PfpI superfamily.

Theoretical study on the allosteric regulation of an oligomeric protease from Pyrococcus horikoshii by Cl- Ion.

The thermophilic intracellular protease (PH1704) from Pyrococcus horikoshii that functions as an oligomer (hexamer or higher forms) has proteolytic activity and remarkable stability. PH1704 is classified as a member of the C56 family of peptidases. This study is the first to observe that the use of Cl− as an allosteric inhibitor causes appreciable changes in the catalytic activity of the protease. Theoretical methods were used for further study. Quantum mechanical calculations indicated the binding mode of Cl− with Arg113. A molecular dynamics simulation explained how Cl− stabilized distinct contact species and how it controls the enzyme activity. The new structural insights obtained from this study are expected to stimulate further biochemical studies on the structures and mechanisms of allosteric proteases. It is clear that the discovery of new allosteric sites of the C56 family of peptidases may generate opportunities for pharmaceutical development and increases our understanding of the basic biological processes of this peptidase family.

The 3D structure of the defense-related rice protein Pir7b predicted by homology modeling and ligand binding studies

AbstractTo better understand the ligand-binding mechanism of protein Pir7b, important part in detoxification of a pathogen-derived compound against Pyricularia oryzae, a 3D structure model of protein Pir7b was constructed based on the structure of the template SABP2. Three substrates were docking to this protein, two of them were proved to be active, and some critical residues are identified, which had not been confirmed by the experiments. His87 and Leu17 considered as ‘oxyanion hole’ contribute to initiating the Ser86 nucleophilic attack. Gln187 and Asp139 can form hydrogen bonds with the anilid group to maintain the active binding orientation with the substrates. The docking model can well interpret the specificity of protein Pir7b towards the anilid moiety of the substrates and provide valuable structure information about the ligand binding to protein Pir7b. FigureLigand binding analysis based on the refined Pir7b model. Magenta dash line, hydrogen bond; Red dash line, distance label. (a) Docking of 2-naphthol AS-acetate to Pir7b model. A 3D figure of 2-naphthol AS-acetate-Pir7b complex is also attached (b) Docking of 2-naphthol AS-2-chlor-propionate to Pir7b model. (c) Docking of 2-naphthol-acetate to Pir7b model.

Investigating the substrate binding mechanism of sulfotransferase 2A1 based on substrate tunnel analysis: a molecular dynamics simulation study

Cytosolic sulfotransferases (SULTs) catalyze the transfer of a sulfonate group from the unique cofactor 3′-phosphoadenosine 5′-phosphosulfate (PAPS) to a large number of diverse substrates. In this work, tunnels that facilitate the transport of substrates in the enzyme were studied, with and without bound cofactor, using extensive molecular dynamics simulations. Residues making up tunnels, as well as residues forming bottlenecks to the tunnels, were identified. Conformation analysis of the active-site cap was also performed. We found that binding of cofactor could significantly narrow the tunnel based on the closing of the active-site cap to the enzyme. The roles of the key residues identified in this work deserve further exploration experimentally.

TOWARD A BLUEPRINT FOR β-PRIMEVEROSIDASE FROM TEA LEAVES STRUCTURE/FUNCTION PROPERTIES: HOMOLOGY MODELING STUDY

By means of the Homology modeling and the known structure of cyannogenic β-glycosidase from white clover (1CBG, EC 3.2.1.21), we construct a 3D model of the β-primeverosidase (EC 3.2.1.149) and search for the binding site of substrate. The 3D model is then refined by using molecular mechanics (optimization and molecular dynamics) simulation. Finally, the refined model is further assessed by Profile-3D and PROCHECK, and the results showed that the final model is reliable. Furthermore, the docking of the substrates into the active site of the protein indicates that β-primeverosidase is able to hydrolyze β-primeverosides, but not act on 2-phenylethyl β-D-glucopyranoside. These results suggest that β-primeverosidase shows broad substrate specificity with respect to the disaccharide glycon moiety (subsite -2). This is consistent with the experimental observation. Thr271 and Thr415 play important roles in subsite -2 of β-primeverosidase. Our results may be helpful for further experimental investigations.