g-Lung Chen

Homology modeling, docking, and molecular dynamics reveal HR1039 as a potent inhibitor of 2009 A(H1N1) influenza neuraminidase

The neuraminidase of the influenza virus is the target of antiviral drugs oseltamivir and zanamivir. Clinical practices have shown that zanamivir and oseltamivir are effective in treating the 2009 A(H1N1) influenza virus. However, drug resistance strains are also emerging. Herein, we report the findings from homology modeling and molecular simulations of 2009 A(H1N1) neuraminidase complexed with zanamivir, oseltamivir, and several herb extracts with potential activities. Our docked oseltamivir and zanamivir results are consistent with previous studies. Based on the same procedure, the docked results of herb extracts HR1039 and HR1040 suggest that they are potential potent inhibitors of neuraminidase. Also, the binding modes of HR1039/HR1040 are different from those of oseltmivir and zanamivir, and may be effective in treating oseltamivir-resistant influenza virus strains.

Potential of mean force for Syrian hamster prion epitope protein – Monoclonal fab 3f4 antibody interaction studies

Simulating antigen-antibody interactions are crucial for understanding antigen-antibody associations in immunology. To shed further light on this question, we study a dissociation of the Syrian hamster prion epitope protein-fab 3f4 antibody complex structure. The stretching, that is, the distance between the center of mass of the prion epitope protein and the fab 3f4 antibody, has been studied using potential of mean force (PMF) calculations based on molecular dynamics (MD) and the implicit water model. For the complex structure, there are four important intermediates, U-shaped groove on the antibodies, and two inter-protein molecular hydrogen bonds in the stretching process. Use of our simulations may help in understanding the binding mechanics of the complex structure, and thus of significance in the design of antibodies against prion disease.

Potential of mean force of the hepatitis C virus core protein-monoclonal 19D9D6 antibody interaction.

Antigen-antibody interactions are critical for understanding antigen-antibody associations in immunology. To shed further light on this question, we studied a dissociation of the 19D9D6-HCV core protein antibody complex structure. However, forced separations in single molecule experiments are difficult, and therefore molecular simulation techniques were applied in our study. The stretching, that is, the distance between the center of mass of the HCV core protein and the 19D9D6 antibody, has been studied using the potential of mean force calculations based on molecular dynamics and the explicit water model. Our simulations indicate that the 7 residues Gly70, Gly72, Gly134, Gly158, Glu219, Gln221 and Tyr314, the interaction region (antibody), and the 14 interprotein molecular hydrogen bonds might play important roles in the antigen-antibody interaction, and this finding may be useful for protein engineering of this antigen-antibody structure. In addition, the 3 residues Gly134, Gly158 and Tyr314 might be more important in the development of bioactive antibody analogs.

Theoretical investigation on reaction of sulbactam with wild-type SHV-1 β-lactamase: acylation, tautomerization, and deacylation.

Molecular dynamics (MD) simulation and quantum mechanical (QM) calculations were used to investigate the reaction mechanism of sulbactam with class A wild-type SHV-1 β-lactamase including acylation, tautomerization, and deacylation. Five different sulbactam-enzyme configurations were investigated by MD simulations. In the acylation step, we found that Glu166 cannot activate Ser70 directly for attacking on the carbonyl carbon, and Lys73 would participate in the reaction acting as a relay. Additionally, we found that sulbactam carboxyl can also act as a general base. QM calculations were performed on the formation mechanism of linear intermediates. We suggest that both imine and trans-enamine intermediates can be obtained in the opening of a five-membered thiazolidine ring. By MD simulation, we found that imine intermediate can exist in two conformations, which can generate subsequent trans- and cis-enamine intermediates, respectively. The QM calculations revealed that trans-enamine intermediate is much more stable than other intermediates. The deacylation mechanism of three linear intermediates (imine, trans-enamine, cis-enamine) was investigated separately. It is remarkably noted that, in cis-enamine intermediate, Glu166 cannot activate water for attacking on the carbonyl carbon directly. This leads to a decreasing of the deacylation rate of cis-enamine. These findings will be potentially useful in the development of new inhibitors.

Structural insights into ligand binding of PGRP1 splice variants in Chinese giant salamander (Andrias davidianus) from molecular dynamics and free energy calculations

Peptidoglycan (PGN) recognition proteins (PGRPs) are important pattern recognition receptors of the innate immune system. A number of PGRP splicing variants produced by alternative splicing of PGRP genes have been reported. However, several important aspects of interactions between PGRP splice variants and their ligands are still unclear. In the present study, three dimensional models of salamander PGRP1 (adPGRP1) and its splice variant (adPGRP1a) were constructed, and their key amino acids involved in interacting with PGNs were analyzed. The results revealed that adPGRP1a has a typical PGRPs structure containing five β-sheets and four α-helices, while adPGRP1 contained five β-sheets and only one α-helix due to the lack of 51 amino acids at its C-terminus. Molecular docking revealed that van der Waals and Coulombic interactions contributed to interactions in the protein–ligand complex. Further binding energy of adPGRP-PGNs computed by the MM-PBSA method revealed that adPGRP1a and adPGRP1 might selectively bind to different PGNs; the former might selectively bind Dap-type PGNs and the latter both types of PGNs. In addition, the binding energy of each residue of adPGRP1a and adPGRP1 was also calculated, revealing that residues involved in the interaction of protein–ligand complexes were different in adPGRP1a and adPGRP1. These results provided a first insight into the potential basis for interaction between PGRPs generated by alternative splicing and PGN derivatives.

Why tazobactam and sulbactam have different intermediates population with SHV-1 β-lactamase: a molecular dynamics study

AbstractThe imine intermediates of tazobactam and sulbactam bound to SHV-1 β-lactamase were investigated by molecular dynamics (MD) simulation respectively. Hydrogen bond networks around active site were found different between tazobactam and sulbactam acyl-enzymes. In tazobactam imine intermediate, it was observed that the triazolyl ring formed stable hydrogen bonds with Asn170 and Thr167. The results suggest that conformation of imine determined the population of intermediates. In imine intermediate of tazobactam, the triazolyl ring is trapped in Thr_Asn pocket, and it restricts the rotation of C5-C6 bond so that tazobactam can only generate trans enamine intermediate. Further, conformational cluster analyses are performed to substantiate the results. These findings provide an explanation for the corresponding experimental results, and will be potentially useful in the development of new inhibitors. FigureThe distribution of dihedral angle N4-C5-C6-C7 in two systems (imine_taz and imine_sul) along MD simulations

COMPUTATIONAL MODELING STUDY ON METABOLISM MECHANISM OF OSELTAMIVIR

Oseltamivir (OTV) is widely used in the treatment of both influenza virus A and B infections. Additionally, OTV is an effective antiviral drug in treating the 2009 A (H1N1) influenza virus. Clinical studies concluded that OTV is readily extensively converted to the active carboxylate metabolite after oral administration. In order to investigate the metabolism mechanism of OTV, we carried out density functional theory (DFT) quantum mechanical calculations. The molecule orbital (MO) theory and natural population analysis (NPA) were also employed to help understanding the reaction mechanism. All possible reaction pathways for OTV metabolism are considered, involving hydrolysis of ester and amide. Two mechanisms were considered in this work, viz. concerted mechanism and stepwise mechanism. Our results indicate the stepwise mechanism is more favorable in both hydrolysis reactions and the rate-determining stage is the formation of the tetrahedral intermediate. In addition, the hydrolysis reactions can be assist...

Computer modeling on the tautomerization of sulbactam intermediate in SHV-1 β-lactamases: E166A mutant vs. wild type.

We present a theoretical study for the tautomerization of sulbactam intermediates in different SHV-1 β-lactamases: E166A and wild-type (WT). Molecular dynamics (MD) simulations were employed and hydrogen bonds network around active site was found different between the WT and E166A acyl-enzymes. In E166A, Asn170 restricts the C5--C6 bond rotation, thus stabilizes the dihedral angle N4--C5--C6--C7 of imine to a trans conformation. The DFT calculations (B3LYP/6-31+G* and B3LYP/6-31++G**) were performed on tautomerization reactions. Two mechanisms including direct and stepwise proton transfer reactions were proposed based on the MD results. In E166A, the substrate carboxyl group acts as a relay station which assists the proton transfer with a very low energy barrier. However, in WT, such stepwise mechanism is difficult to proceed because of the large separation between C6 and substrate carboxyl group. Our results explain why E166A SHV-1 β-lactamases forms greater population of trans-enamine than WT.