Physical quantity of residue electrostatic energy in flavin mononucleotide binding protein dimer.

The electrostatic (ES) energy of each residue was for the first time quantitatively evaluated in a flavin mononucleotide binding protein (FBP). A residue electrostatic energy (RES) was obtained as the sum of the ES energies between atoms in each residue and all other atoms in the FBP dimer using atomic coordinates obtained by a molecular dynamics (MD) simulation. ES is one of the most important energies among the interaction energies in a protein. It is determined from the RES, the residues which mainly contribute to stabilize the structure of each subunit, and the binding energy between two subunits can be estimated. The RES of all residues in subunit A (Sub A) and subunit B (Sub B) were attractive forces, even though the residues contain net negative or positive charges. This reveals that the ES energies of any of the residues can contribute to stabilize the protein structure. The total binding ES energy over all residues among the subunits was distributed between -0.2 to -1.2 eV (mean = -0.67 eV) from the MD simulation time.

Dynamics of the protein structure of T169S pyranose 2-oxidase in solution: Molecular dynamics simulation.

Pyranose 2-oxidase (P2O) from Trametes multicolor contains FAD as cofactor, and forms a tetramer. The protein structure of a mutated P2O, T169S (Thr169 is replaced by Ser), in solution was studied by means of molecular dynamics simulation and analyses of photoinduced electron transfer (ET) from Trp168 to excited isoalloxazine (Iso*), and was compared with wild type (WT) P2O. Hydrogen bonding between Iso and nearby amino acids was very similar as between T169S and WT protein. Distances between Iso and Tyr456 were extremely heterogeneous among the subunits, 1.7 (1.5 in WT) in subunit A (Sub A), 0.97 (2.2 in WT) in Sub B, 1.3 (2.1 in WT) in Sub C, 1.3 nm (2.0 in WT) in Sub D. Mean values of root of mean square fluctuation over all residues were greater by four times than those in WT. This suggests that the protein structure of T169S is much more flexible than that of WT. Electrostatic (ES) energies between Iso anion in one subunit and ionic groups in the entire protein were evaluated. It was found that more than 50% of the total ES energy in each subunit is contributed from other subunits. Reported fluorescence decays were analyzed by a method as WT, previously reported. Electron affinities of Iso* in T169S were appreciably higher than those in WT. Static dielectric constants near Iso and Trp168 were also quite higher in T169S than those in WT.

Conformational difference between two subunits in flavin mononucleotide binding protein dimers from Desulfovibrio vulgaris (MF): molecular dynamics simulation.

The structural and dynamical properties of five FMN binding protein (FBP) dimers, WT (wild type), E13K (Glu13 replaced by Lys), E13R (Glu13 replaced by Arg), E13T (Glu13 replaced by Thr) and E13Q (Glu13 replaced by Gln), were investigated using a method of molecular dynamics simulation (MDS). In crystal structures, subunit A (Sub A) and subunit B (Sub B) were almost completely equivalent in all of the five FBP dimers. However, the predicted MDS structures of the two subunits were not equivalent in solution, revealed by the distances and inter-planar angles between isoalloxazine (Iso) and aromatic amino acids (Trp32, Tyr35 and Trp106) as well as the hydrogen bonding pairs between Iso and nearby amino acids. Residue root of mean square fluctuations (RMSF) also displayed considerable differences between Sub A and Sub B and in the five FBP dimers. The dynamics of the whole protein structures were examined with the distance (RNN) between the peptide N atom of the N terminal (Met1) and the peptide N atom of the C terminal (Leu122). Water molecules were rarely accessible to Iso in all FBP dimers which are in contrast with other flavoenzymes.

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

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.

Anti-HIV-1 integrase effect of compounds from Aglaia andamanica leaves and molecular docking study with acute toxicity test in mice.

CONTEXT: Acquired immunodeficiency syndrome (AIDS) is a serious health problem worldwide. It has been reported that Aglaia andamanica Hiern (Meliaceae) leaves possessed an antiviral effect. Therefore, a search of anti-HIV-1 integrase (HIV-1 IN) agents from A. andamanica is a promising target. OBJECTIVE: The objective of this study is to evaluate anti-HIV-1 IN activity of isolated compounds from A. andamanica using an in vitro assay and molecular docking study as well as testing acute toxicity in mice using the up and down method. MATERIALS AND METHODS: The leaves and compounds (3-100 µg/mL) from A. andamanica were determined for the anti-HIV-1 IN effect using the multiplate integration assay (MIA) by detection the absorbance of the final product, p-nitrophenol, at 405 nm. The molecular docking with the HIV-1 IN of the active compound N-methyl-trans-4-hydroxy-l-proline (10) was also studied. The Swiss albino mice were used for an acute toxicity test. RESULTS AND DISCUSSION: Among the isolated compounds, 10 showed marked anti-HIV-1 IN effect with an IC50 value of 11.8 µg/mL, whereas other compounds were inactive (IC50 value > 100 µg/mL). The molecular docking of compound 10 with an HIV-1 IN enzyme was also studied. The result revealed that this compound formed the hydrogen bonding with the Thr66, Asn155, and Lys159 of the HIV-1 IN binding site. The acute toxicity of the A. andamanica extract was not observed at the dose 2000 mg/kg mice. This is the first report of A. andamanica for anti-HIV-1 IN activity.

Classification of the mechanisms of photoinduced electron transfer from aromatic amino acids to the excited flavins in flavoproteins.

In many flavoproteins photoinduced electron transfer (ET) efficiently takes place from aromatic amino acids such as tryptophan or tyrosine to the excited isoalloxazine, so that the fluorescence lifetimes of isoalloxazine in some flavoproteins become ultrashort. The mechanism of ET in the flavoproteins was classified into four classes from the relationship between logarithmic ET rates (ln Rate) and the donor-acceptor distances (Rc), using reported data. The physical quantity, GT, is defined as the sum of solvent reorganization energy, electrostatic energy between a donor cation and an Iso anion, the standard free energy gap between the photoproducts and reactants, and net electrostatic energy between the photoproducts and other ionic groups in the flavoproteins (NetES). When GT fluctuates around zero with Rc, the ET rate becomes fastest (faster than 1 ps(-1)) in Kakitani and Mataga rates. In the ultrafast ET processes, the ln Rate becomes a parabolic function (category 1) of Rc as in FMN binding proteins and pyranose 2-oxidase at the shorter emission wavelengths, when NetES is negligible compared to the other quantities in the GT function. In the ultrafast ET processes, the ln Rate does not display any clear function of Rc (category 2) when NetES is dominant in the GT function, because of no direct relation between NetES and Rc. ET in flavodoxin from Helicobacter pylori may be classified into category 2. When GT linearly varies with Rc around a certain positive value, the ET rates become much slower (<1 ps(-1)). In this case the ln Rate linearly decreases with Rc (category 3), as Tyr224 in d-amino acid oxidase dimers. It is also conceivable that the ln Rate decreases with much scattered function of Rc (category 4), when NetES is dominant in the GT function, as Tyr314 in d-amino acid oxidase dimers. In ET processes of category 1, ET rates decrease as Rc becomes shorter than the distance at the maximum values of ln Rates, where GT is negative. Conditions and physical meanings were discussed for the GT-negative region.

Anti-HIV-1 integrase activity and molecular docking of compounds from Albizia procera bark.

CONTEXT: Albizia procera (Roxb.) Benth. (Mimosaceae) has been traditionally used in Thai longevity preparations. Thus, searching for HIV-1 integrase (HIV-1 IN) agents from natural sources is of interest. OBJECTIVE: The objective of this study is to examine the inhibitory activity against HIV-1 IN of compounds isolated from the stem bark of Albizia procera. MATERIALS AND METHODS: The EtOH extract and isolated compounds of Albizia procera bark were examined for anti-HIV-1 IN activity at various concentrations (10-100 µg/mL and 10-100 µM) using the multiplate integration assay and molecular docking. RESULTS AND DISCUSSIONS: The results showed that the ethanol extract had good anti-HIV-1 IN activity with an IC50 value of 19.5 µg/mL, whereas ethyl acetate fraction exhibited the most potent with an IC50 value of 19.1 µg/mL, followed by water fraction (IC50 value = 21.3 µg/mL), hexane and chloroform fractions (IC50 value > 100 µg/mL), respectively. From bioassay-guided isolation, the ethyl acetate fraction was further separated to give two compounds which are (+)-catechin (1) and protocatechuic acid (2), respectively. Of the tested samples, (+)-catechin (1) exhibited appreciable activity against HIV-1 IN with an IC50 value of 46.3 µM, whereas protocatechuic acid (2) showed mild activity with 46.0% inhibition at concentration of 100 µM. (+)-Catechin (1) could interact with Thr66, Gly148, and Glu152 in the core domain of IN enzyme, whereas protocatechuic acid (2) could bind with Thr66, His67, Glu152, Asn155, and Lys159. This is the first report on anti-HIV-1 IN activity of Albizia procera bark. These results may suggest that Albizia procera bark has potential as anti-HIV-1 IN agent.

Anti-HIV-1 Integrase Activity and Molecular Docking Study of Compounds from Caesalpinia sappan L.

Caesalpinia sappan L. (Caesalpiniaceae) has been traditionally used as blood tonic, expectorant, and astringent by boiling with water. Searching for HIV-1 integrase (IN) inhibitors from this plant is a promising approach. The EtOH extract of C. sappan and its isolated compounds were tested for their anti-HIV-1 IN effect using the multiplate integration assay, and the active compounds were determined for their mechanisms by molecular docking technique. Extraction from the heartwoods and roots of C. sappan led to the isolation of nine compounds. Among the compounds tested, sappanchalcone (2) displayed the strongest effect against HIV-1 IN with an IC50 value of 2.3 µM followed by protosappanin A (9, IC50 = 12.6 µM). Structure-activity relationships of compounds from C. sappan were found, in which the vicinal hydroxyl moiety were essential for anti-HIV-1 IN effect of compounds 2 and 9 by binding with the amino acid residues Gln148 and Thr66 in the core domain of the HIV- 1 IN enzyme, respectively.

Use of a hexasubstituted benzene scaffold in the development of multivalent HIV-1 integrase inhibitors.

The highly directional hexasubstituted benzene moiety was used as the central scaffold to create new human immunodeficiency virus (HIV)-1 integrase inhibitors through the attachment of multiple active groups. A series of potential inhibitors having substituted polyhydroxylated mono, bis and tris-cinnamoyl derivatives connected on the scaffold were prepared through Claisen-Schmidt condensations with substituted benzaldehydes, followed by partial demethylation to uncover the active phenolic groups required for the interactions with the integrase enzyme active sites. Using a multiplate integration assay method, four compounds carrying at least two sets of interacting moieties were found to be relatively potent integrase inhibitors with IC50 values in the low micromolar range. The results confirmed that multiple polyhydroxylated groups were required on the platform in order to effectively interact with the enzyme. The results from molecular docking studies consistently complemented the experimental results and revealed the nature of the potential key binding interactions responsible for the apparent activity of the active compounds.

The C-terminal region of xylanase domain in Xyn11A from Paenibacillus curdlanolyticus B-6 plays an important role in structural stability.

Paenibacillus curdlanolyticus B-6 produces an extracellular multienzyme complex containing a major xylanase subunit, designated Xyn11A, which includes two functional domains belonging to glycosyl hydrolase family-11 (GH11) and carbohydrate binding module family-36 (CBM36) and possesses a glycine and asparagine-rich linker (linker). To clarify the roles of each functional domain, recombinant proteins XynXL and XynX (CBM36 deleted and CBM36 and linker deleted, respectively) were constructed. Their xylanase activities were similar toward soluble xylan, whereas XynXL showed decreased hydrolysis activity toward insoluble xylan while XynX had no xylanase activity. To determine the significance of the linker and its neighbor region, XynX was subjected to secondary structural alignments using circular dichroism (CD) spectroscopy and three-dimensional (3D) structural analysis. A seven amino acid (NTITIGG) neighbor linker sequence was highly conserved among GH11 xylanases of Paenibacillus species. Although XynX exhibited a typical GH11 xylanase structure, conformational gaps were observed in the ß6- and ß12-sheets and in CD spectra. Flipping of the Arg163 side chains in the subsite was also observed upon analysis of superimposed models. Docking analysis using xylohexaose indicated that flipping of the Arg163 side chains markedly affected substrate binding in the subsite. To identify the amino acids related to stabilizing the substrate binding site, XynX with an extended C-terminal region was designed. At least seven amino acids were necessary to recover substrate binding and xylanase activity. These results indicated that the seven amino acid neighbor Xyn11A linker plays an important role in the activity and conformational stability of the xylanase domain.

Non-equivalent conformations of D-amino acid oxidase dimer from porcine kidney between the two subunits. Molecular dynamics simulation and photoinduced electron transfer.

The structural difference between two subunits of D-amino acid oxidase dimer from porcine kidney was studied by molecular dynamics simulation (MDS) and rate of photoinduced electron transfer (ET) from aromatic amino acids as tyrosines (Tyr) and tryptophanes (Trp) to the excited isoalloxazine (Iso*). The donor-acceptor distances (Rc) between isoalloxazine (Iso) and the donors were shortest in Tyr224 (0.74 nm) in Sub A at 10 °C (Sub A10), in Tyr224 (0.79 nm) in Sub B at 10 °C (Sub B10), in Tyr228 (0.85 nm) in Sub A at 30 °C (Sub A30), and in Tyr224 (0.72 nm) in Sub B at 30 °C (Sub B30). The Rcs were mostly shorter in the dimer than those in the monomer. Hydrogen bonding (H-bond) pairs between Iso and surrounding amino acids varied with the subunit and temperature. O2 of the Iso ring formed an H-bond exclusively with Thr317OG1 (side chain) in both Sub A10 and Sub A30, while it formed with Gly315N (peptide), Leu316N and Thr317N in Sub B10 and Sub B30. N3H of Iso formed an H-bond with Leu51O (peptide) in Sub A10 and Sub A30, but not in Sub B10 and Sub B30. Electron affinity of Iso* was appreciably lower in Sub A10 compared to Sub B10, while it was opposite at 30 °C. ET rate to Iso* was fastest from Tyr224 in Sub A10, while it was fastest from Tyr314 in Sub B10. The ET rate was fastest from Tyr314 in Sub A30, while it was fastest from Tyr224 in Sub B30. The greater ET rates in the dimer as compared to those in the monomer were elucidated with shorter Rc in the dimer as compared to the monomer. The static dielectric constants inside the subunits and the static dielectric constant between Iso and Tyr224 or Tyr228 were not different appreciably. A few water molecules and sometimes an amino acid were located between Iso and Tyr224, which may be the reason why the dielectric constant of the entire subunits did not differ from that between Iso and Tyr224.

Photoinduced electron transfer modeling to simulate flavoprotein fluorescence decay.

A method of analysis is described on the photoinduced electron transfer (PET) from aromatic amino acids as tryptophans (Trp) and tyrosines (Tyr) to the excited isoalloxazine (Iso*) in FMN-binding proteins (FBP) from Desulfovibrio vulgaris (strain, Miyazaki F). Time-dependent geometrical factors as the donor-acceptor distances are determined by means of a molecular dynamics simulation (MDS) of the proteins. Fluorescence decays of the single mutated isoforms of FBP are used as experimental data. The electrostatic (ES) energy between the photoproducts and ionic groups in the proteins is introduced into the Kakitani and Mataga (KM) model, which is modeled for an electron transfer process in solution. The PET parameters contained in the KM rate are determined by means of a nonlinear least square method, according to the Marquardt algorithm. The agreement between the observed and calculated decays is quite good, but not optimal. Characteristics on PET in flavoproteins, obtained by the present method, are described. Possible improvements of the method are discussed.

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.

Simultaneous analyses of photoinduced electron transfer in the wild type and four single substitution isomers of the FMN binding protein from Desulfovibrio vulgaris, Miyazaki F.

The mechanism of photoinduced electron transfer (PET) from the aromatic amino acids (Trp32, Tyr35 and Trp106) to the excited flavin mononucleotide (FMN) in the wild type (WT) and four single amino acid substitution isomers (E13T, E13Q, W32A and W32Y) of FMN binding protein (FBP) from the Desulfovibrio vulgaris (Miyazaki F) were simultaneously analyzed (Method A) with the Marcus-Hush (MH) theory and Kakitani-Mataga (KM) theory using ultrafast fluorescence dynamics of these proteins. In addition, the PET mechanism of the WT, E13T and E13Q FBP systems (Method B) were also analyzed with both MH and KM theories. The KM theory could describe all of the experimental fluorescence decays better than the MH theory by both Methods A and B. The PET rates were found to largely depend on the electrostatic energies between photo-products, isoalloxazine (Iso) anion and the PET donor cations, and the other ionic groups, and hence on static dielectric constants. The dielectric constant (ε(0)(DA)) around the PET donors and acceptor was separately determined from those (ε(0)(j), j = WT, E13T, E13Q, W32Y and W32A) in the domain between the Iso anion or the donor cations and the other ionic groups in the proteins. The values of ε(0)(DA) were always lower than those of ε(0)(j), which is reasonable because no amino acid exists between the PET donors and acceptor in all systems. The values of the dielectric constants ε(0)(j) (j = WT, E13T and E13Q) were similar to those obtained previously from the analysis of the crystal structures and the average lifetimes of these FBP proteins. Energy gap law in the FBP systems was examined. An excellent parabolic function of the logarithms of the PET rates was obtained against the total free energy gap. The PET in these FBP isomers mostly took place in the so-called normal region, and partly in the inverted region.

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.

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.

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.

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.

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.

Hybrid quantum mechanical/molecular mechanical molecular dynamics simulations of HIV-1 integrase/inhibitor complexes.

Human immunodeficiency virus (HIV)-1 integrase (IN) is an attractive target for development of acquired immunodeficiency syndrome chemotherapy. In this study, conventional and coupled quantum mechanical and molecular mechanical (QM/MM) molecular dynamics (MD) simulations of HIV-1 IN complexed with 5CITEP (IN-5CITEP) were carried out. In addition to differences in the bound position of 5CITEP, significant differences at the two levels of theory were observed in the metal coordination geometry and the areas involving residues 116-119 and 140-166. In the conventional MD simulation, the coordination of Mg(2+) was found to be a near-perfect octahedral geometry whereas a distorted octahedral complex was observed in QM/MM. All of the above reasons lead to a different pattern of protein-ligand salt link formation that was not observed in the classical MD simulation. Furthermore to provide a theoretical understanding of inhibition mechanisms of 5CITEP and its derivative (DKA), hybrid QM/MM MD simulations of the two complexes (IN-5CITEP and IN-DKA) have been performed. The results reveal that areas involving residues 60-68, 116-119, and 140-149 were substantially different among the two systems. The two systems show similar pattern of metal coordination geometry, i.e., a distorted octahedron. In IN-DKA, both OD1 and OD2 of Asp-64 coordinate the Mg(2+) in a monodentate fashion whereas only OD1 is chelated to the metal as observed in IN-5CITEP. The high potency of DKA as compared to 5CITEP is supported by a strong salt link formed between its carboxylate moiety and the ammonium group of Lys-159. Detailed comparisons between HIV-1 IN complexed with DKA and with 5CITEP provide information about ligand structure effects on protein-ligand interactions in particular with the Lys-159. This is useful for the design of new selective HIV-1 IN inhibitors.