Hideyoshi Fuji

High-performance drug discovery: computational screening by combining docking and molecular dynamics simulations.

Virtual compound screening using molecular docking is widely used in the discovery of new lead compounds for drug design. However, this method is not completely reliable and therefore unsatisfactory. In this study, we used massive molecular dynamics simulations of protein-ligand conformations obtained by molecular docking in order to improve the enrichment performance of molecular docking. Our screening approach employed the molecular mechanics/Poisson-Boltzmann and surface area method to estimate the binding free energies. For the top-ranking 1,000 compounds obtained by docking to a target protein, approximately 6,000 molecular dynamics simulations were performed using multiple docking poses in about a week. As a result, the enrichment performance of the top 100 compounds by our approach was improved by 1.6–4.0 times that of the enrichment performance of molecular dockings. This result indicates that the application of molecular dynamics simulations to virtual screening for lead discovery is both effective and practical. However, further optimization of the computational protocols is required for screening various target proteins.

Derivatives of 5-nitro-furan-2-carboxylic acid carbamoylmethyl ester inhibit RNase H activity associated with HIV-1 reverse transcriptase.

The RNase H activity associated with human immunodeficiency virus type 1 (HIV-1) is an attractive target for an antiretroviral drug development. We screened 20000 small-molecular-weight compounds for RNase H inhibitors and identified a novel RNase H-inhibiting structure characterized by a 5-nitro-furan-2-carboxylic acid carbamoylmethyl ester (NACME) moiety. Two NACME derivatives, 5-nitro-furan-2-carboxylic acid adamantan-1-carbamoylmethyl ester (compound 1) and 5-nitro-furan-2-carboxylic acid [[4-(4-bromo-phenyl)-thiazol-2-yl]-(tetrahydro-furan-2-ylmethyl)-carbamoyl]-methyl ester (compound 2), effectively blocked HIV-1 and MLV RT-associated RNase H activities with IC(50)s of 3-30 microM but had little effect on bacterial RNase H activity in vitro. Additionally, 20-25 microM compound 2 effectively inhibited HIV-1 replication. An in silico docking simulation indicated that the conserved His539 residue, and two metal ions in the RNase H catalytic center are involved in RNase H inhibition by NACME derivatives. Taken together, these data suggest that NACME derivatives may be potent lead compounds for development of a novel class of antiretroviral drugs.

Mechanism of drug resistance of hemagglutinin of influenza virus and potent scaffolds inhibiting its function.

Highly pathogenic influenza viruses have become a global threat to humans. It is important to select an effective therapeutic option suitable for the subtypes in an epidemic or pandemic. To increase the options, the development of novel antiviral agents acting on targets different from those of the currently approved drugs is required. In this study, we performed molecular dynamics simulations on a spike protein on the viral envelop, hemagglutinin, for the wild-type and three kinds of mutants using a model system consisting of a trimeric hemagglutinin complex, viral lipid membrane, solvation waters, and ions. A natural product, stachyflin, which shows a high level of antiviral activity specific to some subtypes of influenza viruses, was examined on binding to the wild-type hemagglutinin by docking simulation. The compound potency of stachyflin is, however, easily lost due to resistant mutations. From a comparison of simulation results between the wild-type and the resistant mutants, the reason for the drug resistance of hemagglutinin was clarified. Next, 8 compounds were selected from a chemical database by in silico screening, considering the findings from the simulations. Inhibitory activities to suppress the proliferation of influenza virus were measured by cell-based antiviral assays, and two chemical scaffolds were found to be potent for an inhibitor. More than 30 derivatives bearing either of these two chemical scaffolds were synthesized, and cell culture assays were carried out to evaluate the compound potency. Several derivatives displayed a high compound potency, and 50% effective concentrations of two synthesized compounds were below 1 μM.

Structural insight of DNA topoisomerases I from camptothecin-producing plants revealed by molecular dynamics simulations.

DNA topoisomerase I (Top1) catalyzes changes in DNA topology by cleaving and rejoining one strand of the double stranded (ds)DNA. Eukaryotic Top1s are the cellular target of the plant-derived anticancer indole alkaloid camptothecin (CPT), which reversibly stabilizes the Top1-dsDNA complex. However, CPT-producing plants, including Camptotheca acuminata, Ophiorrhiza pumila and Ophiorrhiza liukiuensis, are highly resistant to CPT because they possess point-mutated Top1. Here, the adaptive convergent evolution is reported between CPT production ability and mutations in their Top1, as a universal resistance mechanism found in all tested CPT-producing plants. This includes Nothapodytes nimmoniana, one of the major sources of CPT. To obtain a structural insight of the resistance mechanism, molecular dynamics simulations of CPT- resistant and -sensitive plant Top1s complexed with dsDNA and topotecan (a CPT derivative) were performed, these being compared to that for the CPT-sensitive human Top1. As a result, two mutations, Val617Gly and Asp710Gly, were identified in O. pumila Top1 and C. acuminata Top1, respectively. The substitutions at these two positions, surprisingly, are the same as those found in a CPT derivative-resistant human colon adenocarcinoma cell line. The results also demonstrated an increased linker flexibility of the CPT-resistant Top1, providing an additional explanation for the resistance mechanism found in CPT-producing plants. These mutations could reflect the long evolutionary adaptation of CPT-producing plant Top1s to confer a higher degree of resistance.

Ab initio protein structure prediction with force field parameters derived from water-phase quantum chemical calculation

Molecular dynamics (MD) simulations are extensively used in the study of the structures and functions of proteins. Ab initio protein structure prediction is one of the most important subjects in computational biology, and many trials have been performed using MD simulation so far. Since the results of MD simulations largely depend on the force field, reliable force field parameters are indispensable for the success of MD simulation. In this work, we have modified atom charges in a standard force field on the basis of water‐phase quantum chemical calculations. The modified force field turned out appropriate for ab initio protein structure prediction by the MD simulation with the generalized Born method. Detailed analysis was performed in terms of the conformational stability of amino acid residues, the stability of secondary structure of proteins, and the accuracy for prediction of protein tertiary structure, comparing the modified force field with a standard one. The energy balance between α‐helix and β‐sheet structures was significantly improved by the modification of charge parameters.

Discovery of novel scaffolds for γ-secretase modulators without an arylimidazole moiety

Gamma-secretase modulators (GSMs) selectively inhibit the production of amyloid-β 42 (Aβ42) and may therefore be useful in the management of Alzheimers disease. Most heterocyclic GSMs that are not derived from nonsteroidal anti-inflammatory drugs contain an arylimidazole moiety that potentially inhibits cytochrome P450 (CYP) activity. Here, we discovered imidazopyridine derivatives that represent a new class of scaffold for GSMs, which do not have a strongly basic end group such as arylimidazole. High-throughput screening identified 2-methyl-8-[(2-methylbenzyl)oxy]-3-(pyridin-4-yl)imidazo[1,2-a]pyridine (3a), which inhibited the cellular production of Aβ42 (IC50 = 7.1 µM) without changing total production of Aβ. Structural optimization of this series of compounds identified 5-[8-(benzyloxy)-2-methylimidazo[1,2-a]pyridin-3-yl]-2-ethylisoindolin-1-one (3m) as a potent inhibitor of Aβ42 (IC50 = 0.39 µM) but not CYP3A4. Further, 3m demonstrated a sustained pharmacokinetic profile in mice and sufficiently penetrated the brain.

Prediction of Ligand Binding Affinity to Target Proteins by Molecular Mechanics Theoretical Calculation

Accurate estimation of ligand-receptor binding affinity is indispensable for computer-assisted drug discovery and structure-based drug design. Many computational scoring functions for estimating binding affinity have been proposed. Every scoring function reported so far, however, has strengths and weaknesses depending on the chemical properties of ligands and the feature of the binding site of the receptor. Hence, potential functions that can be used for many kinds of target proteins are required. In this work, we developed a software program based on Morse-type potential functions that enables evaluation of binding affinity and geometry optimization. Eight different kinds of proteins were used as test data, and ligand chemicals for which the binding pose to the protein and inhibitory constant are known were selected for evaluation. The calculated binding score and the experimentally measured inhibitory constant showed good compatibilities for six target proteins but poor correlation for one target. These compatibilities were compared with the results obtained by using two other software programs. The comparison suggested that the performance of the software developed in this work is good. Since the software can be handled in a computer facility with a many-core system, the software will be effective for search for an active compound from a chemical database and for assistance in chemical modification of the active compound in the pharmaceutical research field.