Somsak Pianwanit

How amantadine and rimantadine inhibit proton transport in the M2 protein channel

To understand how antiviral drugs inhibit the replication of influenza A virus via the M2 ion channel, molecular dynamics simulations have been applied to the six possible protonation states of the M2 ion channel in free form and its complexes with two commercial drugs in a fully hydrated lipid bilayer. Among the six different states of free M2 tetramer, water density was present in the pore of the systems with mono-protonated, di-protonated at adjacent position, tri-protonated and tetra-protonated systems. In the presence of inhibitor, water density in the channel was considerably better reduced by rimantadine than amantadine, agreed well with the experimental IC(50) values. With the preferential position and orientation of the two drugs in all states, two mechanisms of action, where the drug binds to the opening pore and the histidine gate, were clearly explained, i.e., (i) inhibitor was detected to localize slightly closer to the histidine gate and can facilitate the orientation of His37 imidazole rings to lie in the close conformation and (ii) inhibitor acts as a blocker, binding at almost above the opening pore and interacts slightly with the three pore-lining residues, Leu26, Ala30 and Ser31. Here, the inhibitors were found to bind very weakly to the channel due to their allosteric hindrance while theirs side chains were strongly solvated.

Homology modeling and molecular dynamics simulations of Dengue virus NS2B/NS3 protease: insight into molecular interaction.

The pathogenic West Nile virus (WNV) and Dengue virus (DV) are growing global threats for which there are no specific treatments. Both viruses possess a two component NS2B/NS3 protease which cleaves viral precursor proteins. Whereas for the WNV protease two crystal structures in complex with an inhibitor have been solved recently, no such information is available for the DV protease. Here, we report the generation of a homology model of DV NS2B/NS3 protease. Since it is known from the related WNV protease that it adopts a distinct conformation in free and in inhibitor‐complexed form, a special emphasis was given to the analysis of the protease flexibility. Therefore, several models of DV NS2B/NS3 protease complexed with the peptidic inhibitor (Bz‐Nle(P4)‐Lys(P3)‐Arg(P2)‐Arg(P1)‐H) were generated. The first DV protease model (DV‐1) was constructed using the available crystal structure of the apo DV NS2B/NS3 protease. The second model (DV‐2) was built taking the WNV NS3/NS2B protease in the inhibitor‐complexed form as the template structure. Molecular dynamics simulations which were carried out for the WNV crystal structures as well as for the DV models provided an understanding of the role of NS2B for maintaining the protease in the active conformation. It was also demonstrated that NS2B is not only important for maintaining NS3 in the active form, but is also essential for establishing the interaction between residues from the S2 pocket and the peptidic inhibitor. The DV NS2B/NS3 model in the productive conformation can now be used for structure‐based design purposes. Copyright

Source of High Pathogenicity of an Avian Influenza Virus H5N1: Why H5 Is Better Cleaved by Furin

The origin of the high pathogenicity of an emerging avian influenza H5N1 due to the -RRRKK- insertion at the cleavage loop of the hemagglutinin H5, was studied using the molecular dynamics technique, in comparison with those of the noninserted H5 and H3 bound to the furin (FR) active site. The cleavage loop of the highly pathogenic H5 was found to bind strongly to the FR cavity, serving as a conformation suitable for the proteolytic reaction. With this configuration, the appropriate interatomic distances were found for all three reaction centers of the enzyme-substrate complex: the arrangement of the catalytic triad, attachment of the catalytic Ser(368) to the reactive S1-Arg, and formation of the oxyanion hole. Experimentally, the--RRRKK--insertion was also found to increase in cleavage of hemagglutinin by FR. The simulated data provide a clear answer to the question of why inserted H5 is better cleaved by FR than the other subtypes, explaining the high pathogenicity of avian influenza H5N1.

Postprocessing of protein-ligand docking poses using linear response MM-PB/SA: application to Wee1 kinase inhibitors.

Prediction of the binding strength of untested ligands is a central issue in structure-based drug design. In order to rapidly screen large compound databases, simple scoring schemes are often used in target-based virtual screening. The resulting scores often correlate poorly with biological affinities. More rigorous scoring methods, such as MM-PB/SA, correlate better with biological data by considering solvation effects and protein flexibility in the calculation of the binding free energy of a ligand. Here we describe the performance of a modified MM-PB/SA method on 222 Wee1 kinase inhibitors (48 pyridopyrimidine and 174 pyrrolocarbazole derivatives). Docking of these inhibitors into the available Wee1 kinase crystal structure yielded a consistent binding mode, and the derived MM-PB/SA models showed a significant correlation between calculated and experimental data (r(2) values between 0.64 and 0.67). Further study of these models on external test sets of Wee1 kinase inhibitors and structurally related decoys showed that a model based on a single kinase-inhibitor conformation can discriminate the active inhibitors from decoys. We also tested whether the linear interaction energy method with continuum electrostatics (LIECE) yields comparable results to MM-PB/SA and whether the LIECE and MM-PB/SA models can be applied for virtual screening of compound libraries.

Structural basis for the temperature-induced transition of D-amino acid oxidase from pig kidney revealed by molecular dynamic simulation and photo-induced electron transfer.

The structural basis for the temperature-induced transition in the D-amino acid oxidase (DAAO) monomer from pig kidney was studied by means of molecular dynamic simulations (MDS). The center to center (Rc) distances between the isoalloxazine ring (Iso) and all aromatic amino acids (Trp and Tyr) were calculated at 10 °C and 30 °C. Rc was shortest in Tyr224 (0.82 and 0.88 nm at 10 and 30 °C, respectively), and then in Tyr228. Hydrogen bonding (H-bond) formed between the Iso N1 and Gly315 N (peptide), between the Iso N3H and Leu51 O (peptide) and between the Iso N5 and Ala49 N (peptide) at 10 °C, whilst no H-bond was formed at the Iso N1 and Iso N3H at 30 °C. The H-bond of Iso O4 with Leu51 N (peptide) at 10 °C switched to that with Ala49 N (peptide) at 30 °C. The reported fluorescence lifetimes (228 and 182 ps at 10 and 30 °C, respectively) of DAAO were analyzed with Kakitani and Mataga (KM) ET theory. The calculated fluorescence lifetimes displayed an excellent agreement with the observed lifetimes. The ET rate was fastest from Tyr224 to the excited Iso (Iso*) at 10 °C and from Tyr314 at 30 °C, despite the fact that the Rc was shortest between Iso and Tyr224 at both temperatures. This was explained by the electrostatic energy in the protein. The differences in the observed fluorescence lifetimes at 10 and 30 °C were ascribed to the differences in electron affinity of the Iso* at both temperatures, in which the free energies of the electron affinity of Iso* at 10 and 30 °C were -8.69 eV and -8.51 eV respectively. The other physical quantities related to ET did not differ appreciably at both temperatures. The electron affinities at both temperatures were calculated with a semi-empirical molecular orbital method (MO) of PM6. Mean calculated electron affinities over 100 snapshots with 0.1 ps intervals were -7.69 eV at 10 °C and -7.59 eV at 30 °C. The difference in the calculated electron affinities, -0.11 eV, was close to the observed difference in the free energies, -0.18 eV. The present quantitative analysis predicts that the highest ET rate can occur from a donor with longer donor-acceptor distance, which was explained by differences in electrostatic energy.

Receptor-based 3D-QSAR studies of checkpoint Wee1 kinase inhibitors

One hundred and seventy-four pyrrolo[3,4-c]carbazole-1,3(2H,6H)-dione derivatives reported as inhibitors of the kinase Wee1 were used for a molecular docking and three-dimensional quantitative structure-activity relationship (3D-QSAR) study. Due to the availability of the three-dimensional structure of the Wee1 kinase a receptor-based alignment strategy was applied. Six available Wee1-inhibitor crystal structures were analyzed using the docking program GOLD resulting in a good reproduction of the experimentally derived position and interaction of the cocrystallized inhibitors. Since only a low correlation between docking scores and inhibitory activities was obtained for the series of 174 inhibitors a receptor-based 3D-QSAR study was performed, dividing the data set into 144 training set molecules and an external test set of 30 compounds. Besides the ligand alignment derived from the docking study we tested several other alignment procedures as basis for the 3D-QSAR analysis. The most predictive model was obtained using the alignment from the GOLD docking study. The CoMFA model was found to be robust (q(LOO)(2)=0.764 and r(2)=0.870). The predictive ability of the model was further examined by carrying out leave-20%-out and leave-50%-out cross-validation (q(2)=0.747 for leave-20%-out and 0.737 for leave-50%-out) and predicting the activities of 30 inhibitors used as external test set (r(pred)(2)=0.790). The graphical analysis of the CoMFA contour plot together with the key residues of the binding pocket provided important insight into the relevant interactions of the inhibitors. The results not only provide information about the essential features of potent Wee1 inhibitors but also show the advantage of using receptor-based alignment for 3D-QSAR analysis.

Homology modelling and molecular dynamics simulations of wild type and mutated flavodoxins from Desulfovibrio vulgaris (Miyazaki F): insight into FMN–apoprotein interactions

The flavodoxin from Desulfovibrio vulgaris, strain Miyazaki F (FD-DvMF), binds one molecule of flavin mononucleotide (FMN) as a cofactor and is considered to associate with electron transport reactions. However, although the 3D structure of the related FD from D. vulgaris strain Hildenborough has been determined, that for FD-DvMF has not. In this study, we have predicted the protein structures of the wild type and the W59F, Y97F and W59F-Y97F substitutional mutants of FD-DvMF by a homology modelling approach. Subsequently, the dynamic properties of these four FD-DvMF variants were investigated by molecular dynamics simulations. The results revealed that peptide O of Trp59 formed H-bond with Tyr99OH only in Y97F, leading to FMN being buried deeper inside the protein than in the other three variants and reducing the accessibility of water to FMN. The phosphate oxygen atoms formed extensive H-bonds with amino acid residues in the 10-loop region in all variants resulting in the highest degree of stabilisation. The OH groups of the ribityl chain and the isoalloxazine ring formed H-bonds with amino acid residues of the 60- and 90-loop regions, respectively. The decomposition free energy calculations suggest that the greatest contribution come from the 10-loop region, which is compatible with the published data. That the calculated binding and decomposition free energies were both greatest in Y97F is proposed to be due to the H-bond between peptide O of Trp59 and Tyr99OH.

Theoretical analyses of the fluorescence lifetimes of the D-amino acid oxidase–benzoate complex dimer from porcine kidney: molecular dynamics simulation and photoinduced electron transfer

The mechanism of photoinduced electron transfer (ET) from benzoate (Bz) and aromatic amino acids to the excited isoalloxazine (Iso*) in the D-amino acid oxidase–benzoate complex (DAOB) dimer from porcine kidney was studied using molecular dynamics simulation (MDS) and an electron transfer theory, and compared with that in the DAOB monomer. The DAOB dimer displayed two fluorescent lifetime components of 0.85 ps and 4.8 ps, as reported. The ET parameters contained in the Kakitani and Mataga (KM) model were determined so as to reproduce these lifetimes with MDS atomic coordinates. The Bz–isoalloxazine (Iso) distances were 0.66 nm in subunit A (Sub A), 0.68 nm in subunit B (Sub B) and 0.61 nm in the monomer. The fluorescent lifetimes of 4.8 ps and 0.85 ps were found to originate from Sub A and Sub B, respectively. In Sub A, Tyr228 was the fastest ET donor followed by Bz and Tyr55, while Bz was followed by Tyr228 and Tyr314 in Sub B. The ET rate from Bz was fastest in Sub B, followed by that in Sub A and the DAOB monomer. The static dielectric constants obtained near Iso were 2.4–2.6 in the DAOB dimer and monomer and 5.8–5.9 in holo D-amino oxidase (DAAO). The different dielectric constants could account for the experimental fluorescence peak observed for DAOB (524 nm) and DAAO (530 nm). Logarithmic ET rates decreased linearly with the donor–acceptor distance expressed by both center to center distance (Rc) and edge to edge distance (Re) in Sub A and Sub B of DAOB dimer and monomer, which reveals that the conventional Dutton rule holds in the ET processes in DAOB. The logarithmic ET rates were decomposed into the electronic coupling (EC), square root (SQ) and exponential (GTRAM) terms. It was found that both the EC term and the GTRAM term also decreased linearly with Rc. The sum of the slopes in the EC and GTRAM vs. Rc plots coincided with the slopes in the logarithmic ET rate vs. Rc functions, suggesting that the GTRAM term makes a significant contribution to the linear relations between logarithmic ET rate and Rc.

Computational Studies of HIV-1 Integrase and its Inhibitors

Integration of the genome of the human immunodeficiency virus (HIV) into that of the host genome is catalyzed by HIV integrase (IN) and is an essential step in HIV-1 life cycle. Therefore, drug discovery efforts have been undertaken to identify selective IN inhibitors with the goal of improving the outcome of AIDS therapy using Highly Active Anti Retroviral Therapy (HAART). As computational technology has grown rapidly and is increasingly being used worldwide to accelerate the drug discovery processes, the aim of this review is to summarize the applications of the computer-aided drug design (CADD) techniques to HIV-1 IN and its inhibitors. The following applications are emphasized, including two- and three-dimensional quantitative structure activity relationships (2D/3D-QSAR), pharmacophore modeling, database searching, molecular docking, molecular dynamics simulations, and de novo methodologies.

Anti-HIV-1 integrase compounds from Dioscorea bulbifera and molecular docking study

Abstract Context Dioscorea bulbifera L. (Dioscoreaceae) has been used in a traditional Thai longevity medicine preparation. Isolation of inhibitors from natural products is a potential source for continuous development of new HIV-1 integrase (IN) inhibitors. Objective The objective of this study is to isolate the compounds and evaluate their anti-HIV-1 IN activity, as well as to predict the potential interactions of the compounds with an IN. Materials and methods The ethyl acetate and water fractions (1–100 μg/mL) of Dioscorea bulbifera bulbils were isolated and tested for their anti-HIV-1 IN activity using the multiplate integration assay (MIA). The interactions of the active compounds with IN were investigated using a molecular docking method. Results and discussions The ethyl acetate and water fractions of Dioscorea bulbifera bulbils afforded seven compounds. Among these, allantoin (1), 2,4,3′,5′-tetrahydroxybibenzyl (2), and 5,7,4′-trihydroxy-2-styrylchromone (5) were isolated for the first time from this plant. Myricetin (4) exhibited the most potent activity with an IC50 value of 3.15 μM, followed by 2,4,6,7-tetrahydroxy-9,10-dihydrophenanthrene (3, IC50 value= 14.20 μM), quercetin-3-O-β-d-glucopyranoside (6, IC50 value = 19.39 μM) and quercetin-3-O-β-d-galactopyranoside (7, IC50 value = 21.80 μM). Potential interactions of the active compounds (3, 4, 6, and 7) with the IN active site were additionally investigated. Compound 4 showed the best binding affinity to IN and formed strong interactions with various amino acid residues. These compounds interacted with Asp64, Thr66, His67, Glu92, Asp116, Gln148, Glu152, Asn155, and Lys159, which are involved in both the 3′-processing and strand transfer reactions of IN. In particular, galloyl, catechol, and sugar moieties were successful inhibitors for HIV-1 IN.