Molecular Modeling Studies ofYersinia pestisDihydrofolate Reductase

Microwave-assisted Synthesis and Docking Studies of Phenylureas as Candidates for the Drug Design Against the Biological Warfare Agent Yersinia Pestis

Background:

Bubonic plague is amongst the diseases with the highest potential for being used in biological warfare attacks today. This disease, caused by the bacterium Yersina pestis, is highly infectious and can achieve 100% of fatal victims when in its most dangerous form. Besides, there is no effective vaccine, and the chemotherapy available today against plague is ineffective if not administered at the beginning of the infection.

Objective:

Willing to contribute for changing this reality we propose here new phenylureas as candidates for the drug design against plague meant to target the enzyme dihydrofolate reductase from Y. pestis (YpDHFR).

Methods:

Seven phenylureas, four of them new, were synthesized, following synthetic routes adapted from procedures available in the literature, and using microwave irradiation. After, they were submitted to docking studies inside YpDHFR and human DHFR (HssDHFR) in order to check their potential as selective inhibitors.

Results:

Our results revealed four new phenylureas and a new synthetic route for this kind of molecule using microwave irradiation. Also, our docking studies showed that these compounds are capable of binding to both HssDHFR and YpDHFR, with U1 - U4 and U23 showing more selectivity for HssDHFR and U7, U8 being more selective towards YpDHFR.

Conclusion:

We reported the synthesis with good yields of seven phenylureas, following a simple and clean alternative synthetic route using microwave irradiation. Further molecular docking studies of our compounds suggested that two are capable of binding more selectivity to YpDHFR, qualifying as potential candidates for the drug design of new drugs against plague.

A multilayer screening approach toward the discovery of novel Pf-DHFR inhibitors

A small yet diverse xanthone library was build and computationally docked against wild type Pf-DHFR by Molegro Virtual Docker (MolDock). For analysis of results an integrated approach based on re-ranking, scaling (based on heavy atom counts), pose clustering and visual inspection was implemented. Standard methods such as self-docking (for docking), EF analysis, average rank determinations (for size normalization), and cluster quality indices (for pose clustering) were used for validation of results. Three compounds X5, X113A and X164B displayed contact footprints similar to the known inhibitors with good scores. Finally, 16 compounds were extracted from ZINC data base by similarity based screening, docking score and drug/lead likeness. Out of these 16 compounds, 11 displayed very close contact footprints to experimentally known inhibitors, indicating there potential utility in further drug discovery efforts.

Analyses of the Binding between Water Soluble C60 Derivatives and Potential Drug Targets through a Molecular Docking Approach

Fullerene C60, a unique sphere-shaped molecule consisting of carbon, has been proved to have inhibitory effects on many diseases. However, the applications of C60 in medicine have been severely hindered by its complete insolubility in water and low solubility in almost all organic solvents. In this study, the water-soluble C60 derivatives and the C60 binding protein’s structures were collected from the literature. The selected proteins fall into several groups, including acetylcholinesterase, glutamate racemase, inosine monophosphate dehydrogenase, lumazine synthase, human estrogen receptor alpha, dihydrofolate reductase and N-myristoyltransferase. The C60 derivatives were docked into the binding sites in the proteins. The binding affinities of the C60 derivatives were calculated. The bindings between proteins and their known inhibitors or native ligands were also characterized in the same way. The results show that C60 derivatives form good interactions with the binding sites of different protein targets. In many cases, the binding affinities of C60 derivatives are better than those of known inhibitors and native ligands. This study demonstrates the interaction patterns of C60 derivatives and their binding partners, which will have good impact on the fullerene-based drug discovery.

Molecular dynamics of interactions between rigid and flexible antifolates and dihydrofolate reductase from pyrimethamine-sensitive and pyrimethamine-resistant Plasmodium falciparum

Currently, the usefulness of antimalarials such as pyrimethamine (PYR) is drastically reduced due to the emergence of resistant Plasmodium falciparum (Pf) caused by its dihydrofolate reductase (PfDHFR) mutations, especially the quadruple N51I/C59R/S108N/I164L mutations. The resistance was due to the steric conflict of PYR with S108N. WR99210 (WR), a dihydrotriazine antifolate with a flexible side chain that can avoid such conflict, can overcome this resistance through tight binding with the mutant. To understand factors contributing to different binding affinities of PYR/WR to the wild type (WT) and quadruple mutant (QM), we performed simulations on WR-WT, WR-QM, PYR-WT, and PYR-QM complexes and found that Ile14 and Asp54 were crucial for PYR/WR binding to PfDHFR due to strong hydrogen bonds. The quadruple mutations cause PYR to form, on average, fewer hydrogen bonds with Ile14 and Leu164, and to be displaced from its optimal orientation for Asp54 interaction. The predicted binding affinity ranking (WR-QM ≈ WR-WT ≈ PYR-WT >> PYR-QM) reasonably agrees with the inhibition constant (K(i)) ranking. Our results reveal important residues for tight binding of PYR/WR to WT/QM, which may be used to evaluate the inhibition effectiveness of antimalarials and to provide fundamental information for designing new drugs effective against drug-resistant P. falciparum.

Insight into the molecular mechanism about lowered dihydrofolate binding affinity to dihydrofolate reductase-like 1 (DHFRL1)

Human dihydrofolate reductase-like 1 (DHFRL1) has been identified as a second human dihydrofolate reductase (DHFR) enzyme. Although DHFRL1 have high sequence homology with human DHFR, dihydrofolate (DHF) exhibits a lowered binding affinity to DHFRL1 and the corresponding molecular mechanism is still unknown. To address this question, we studied the binding of DHF to DHFRL1 and DHFR by using molecular dynamics simulation. Moreover, to investigate the role the 24th residue of DHFR/DHFRL1 plays in DHF binding, R24W DHFRL1 mutant was also studied. The van der Waals interaction are more crucial for the total DHF binding energies, while the difference between the DHF binding energies of human DHFR and DHFRL1 can be attributed to the electrostatic interaction and the polar desolvation free energy.More specifically, lower DHF affinity to DHFRL1 can be mainly attributed to the reduction of net electrostatic interactions of residues Arg32 and Gln35 of DHFRL1 with DHF as being affected by Arg24. The side chain of Arg24 in DHFRL1 can extend deeply into the binding sites of DHF and NADPH, and disturb the DHF binding by steric effect, which rarely happens in human DHFR and R24W DHFRL1 mutant. Additionally, the conformation of loop I in DHFRL1 was also studied in this work. Interestingly, the loop conformation resemble to normal closed state of Escherichia coli DHFR other than the closed state of human DHFR. We hope this work will be useful to understand the general characteristics of DHFRL1.

Molecular modeling study of dihydrofolate reductase inhibitors. Molecular dynamics simulations, quantum mechanical calculations, and experimental corroboration

A molecular modeling study on dihydrofolate reductase (DHFR) inhibitors was carried out. By combining molecular dynamics simulations with semiempirical (PM6), ab initio, and density functional theory (DFT) calculations, a simple and generally applicable procedure to evaluate the binding energies of DHFR inhibitors interacting with the human enzyme is reported here, providing a clear picture of the binding interactions of these ligands from both structural and energetic viewpoints. A reduced model for the binding pocket was used. This approach allows us to perform more accurate quantum mechanical calculations as well as to obtain a detailed electronic analysis using the quantum theory of atoms in molecules (QTAIM) technique. Thus, molecular aspects of the binding interactions between inhibitors and the DHFR are discussed in detail. A significant correlation between binding energies obtained from DFT calculations and experimental IC₅₀ values was obtained, predicting with an acceptable qualitative accuracy the potential inhibitor effect of nonsynthesized compounds. Such correlation was experimentally corroborated synthesizing and testing two new inhibitors reported in this paper.

ITP Adjuster 1.0: A New Utility Program to Adjust Charges in the Topology Files Generated by the PRODRG Server

The suitable computation of accurate atomic charges for the GROMACS topology *.itp files of small molecules, generated in the PRODRG server, has been a tricky task nowadays because it does not calculate atomic charges using an ab initio method. Usually additional steps of structure optimization and charges calculation, followed by a tedious manual replacement of atomic charges in the *.itp file, are needed. In order to assist this task, we report here the ITP Adjuster 1.0, a utility program developed to perform the replacement of the PRODRG charges in the *.itp files of small molecules by ab initio charges.

Structural analysis of Pla protein from the biological warfare agent Yersinia pestis: docking and molecular dynamics of interactions with the mammalian plasminogen system

Yersinia pestis protein Pla is a plasmid-coded outer membrane protein with aspartic-protease activity. Pla exhibits a plasminogen (Plg) activator activity (PAA) that promotes the cleavage of Plg to the active serine-protease form called plasmin. Exactly how Pla activates Plg into plasmin remains unclear. To investigate this event, we performed the interactions between the predicted Plg and Pla protein structures by rigid-body docking with the HEX program and evaluated the complex stability by molecular dynamics (MD) using the GROMACS package programs. The predicted docked complex of Plg-Pla shows the same interaction site predicted by experimental site-direct mutagenesis in other studies. After a total of 8 ns of MD simulation, we observed the relaxation of the beta-barrel structure of Pla and the progressive approximation and stabilization between the cleavage site of Plg into the extracellular loops of Pla, followed by the increase in the number of H bonds. We also report here the aminoacids that participate in the active site and the sub sites of interaction. The total understanding of these interactions can be an important tool for drug design against bacterial proteases.