Virtual screening,docking, and dynamics of potential new inhibitors of dihydrofolate reductase from Yersinia pestis

In the present work, we propose to design drugs that target the enzyme dihydrofolate redutase (DHFR) as a means of a novel drug therapy against plague. Potential inhibitors of DHFR from Yersinia pestis (YpDHFR) were selected by virtual screening and subjected to docking, molecular dynamics (MD) simulations, and Poisson–Boltzmann surface area method, in order to evaluate their interactions in the active sites of YpDHFR and human DHFR (HssDHFR). The results suggested selectivity for three compounds that were further used to propose the structures of six new potential selective inhibitors for YpDHFR.     

Docking and molecular dynamics studies of peripheral site ligand–oximes as reactivators of sarin-inhibited human acetylcholinesterase

In the present work, we performed docking and molecular dynamics simulations studies on two groups of long-tailored oximes designed as peripheral site binders of acetylcholinesterase (AChE) and potential penetrators on the blood brain barrier. Our studies permitted to determine how the tails anchor in the peripheral site of sarin-inhibited human AChE, and which aminoacids are important to their stabilization. Also the energy values obtained in the docking studies corroborated quite well with the experimental results obtained before for these oximes.     

RECENT EVOLUTION AND DIVERGENCE AMONG POPULATIONS OF A RARE MEXICAN ENDEMIC,CHIHUAHUA SPRUCE,FOLLOWING HOLOCENE CLIMATIC WARMING

Fragmentation and reduction in population size are expected to reduce genetic diversity. However, examples from natural populations of forest trees are scarce. The range of Chihuahua spruce retreated northward and fragmented coincident with the warming climate that marked the early Holocene. The isolated populations vary from 15 to 2441 trees, which provided an opportunity to test whether census number is a good predictor of genetic diversity. Mean expected heterozygosity, H_e, based on 24 loci in 16 enzyme systems, was 0.093 for 10 sampled populations, which is within the range reported for conifers. However, estimates varied more than twofold among populations and H_e was closely related to the logarithm of the number of mature trees in the population (r_He,N = 0.93). Diversity among populations, F_ST, was 24.8% of the total diversity, which is higher than that observed in almost all conifer species studied. Nei's genetic distance, D, was not related to geographic distance between populations, and D? was 0.033, which is higher than estimates for most wide-ranging species. Most populations had excess homozygosity and the fixation index, F_IS, was higher than that reported for all but one species of conifer. Nm, the number of migrants per generation, was 0.43 to 0.76, depending on estimation procedure, and is the smallest observed in conifers. The data suggest that populations of Chihuahua spruce have differentiated by drift and that they are effectively isolated. The results illustrate how a combination of paleontological observation and molecular markers can be used to illuminate recent evolutionary events. Multilocus estimates of outcrossing for two small populations were zero (complete selfing) and 0.153, respectively, which are in striking contrast to the near complete outcrossing observed in most conifers. The high fixation index and a high proportion of empty seeds (45%) suggest that inbreeding may be a serious problem for conservation of Chihuahua spruce.     

Tyrosyl-DNA Phosphodiesterase I Catalytic Mutants Reveal an Alternative Nucleophile That Can Catalyze Substrate Cleavage

Tyrosyl-DNA phosphodiesterase I (Tdp1) catalyzes the repair of 3′-DNA adducts, such as the 3′-phosphotyrosyl linkage of DNA topoisomerase I to DNA. Tdp1 contains two conserved catalytic histidines: a nucleophilic His (His^nuc) that attacks DNA adducts to form a covalent 3′-phosphohistidyl intermediate and a general acid/base His (His^gab), which resolves the Tdp1-DNA linkage. A His^nuc to Ala mutant protein is reportedly inactive, whereas the autosomal recessive neurodegenerative disease SCAN1 has been attributed to the enhanced stability of the Tdp1-DNA intermediate induced by mutation of His^gab to Arg. However, here we report that expression of the yeast His^nucAla (H182A) mutant actually induced topoisomerase I-dependent cytotoxicity and further enhanced the cytotoxicity of Tdp1 His^gab mutants, including H432N and the SCAN1-related H432R. Moreover, the His^nucAla mutant was catalytically active in vitro, albeit at levels 85-fold less than that observed with wild type Tdp1. In contrast, the His^nucPhe mutant was catalytically inactive and suppressed His^gab mutant-induced toxicity. These data suggest that the activity of another nucleophile when His^nuc is replaced with residues containing a small side chain (Ala, Asn, and Gln), but not with a bulky side chain. Indeed, genetic, biochemical, and mass spectrometry analyses show that a highly conserved His, immediately N-terminal to His^nuc, can act as a nucleophile to catalyze the formation of a covalent Tdp1-DNA intermediate. These findings suggest that the flexibility of Tdp1 active site residues may impair the resolution of mutant Tdp1 covalent phosphohistidyl intermediates and provide the rationale for developing chemotherapeutics that stabilize the covalent Tdp1-DNA intermediate.     

Analysis of the structure and dynamics of human serum albumin

A computational approach to designing a peptide-based ligand for the purification of human serum albumin (HSA) was undertaken using molecular docking and molecular dynamics (MD) simulation. A three-step procedure was performed to design a specific ligand for HSA. Based on the candidate pocket structure of HSA (warfarin binding site), a peptide library was built. These peptides were then docked into the pocket of HSA using the GOLD program. The GOLDscore values were used to determine the affinity of peptides for HSA. Consequently, the dipeptide Trp–Trp, which shows a high GOLDscore value, was selected and linked to a spacer arm of Lys[CO(CH₂)₅NH] on the surface of ECH-lysine sepharose 4 gel. For further evaluation, the Autodock Vina program was used to dock the linked compound into the pocket of HSA. The docking simulation was performed to obtain a first guess of the binding structure of the spacer–Trp–Trp–HSA complex and subsequently analyzed by MD simulations to assess the reliability of the docking results. These MD simulations indicated that the ligand–HSA complex remains stable, and water molecules can bridge between the ligand and the protein by hydrogen bonds. Finally, absorption spectroscopic studies were performed to illustrate the appropriateness of the binding affinity of the designed ligand toward HSA. These studies demonstrate that the designed dipeptide can bind preferentially to the warfarin binding site. Graphical Abstract

Three-step computational approach to the design of a dipeptide ligand for human serum albumin purification exploiting structure-based docking and molecular dynamics simulation     

The role of the oximes HI-6 and HS-6 inside human acetylcholinesterase inhibited with nerve agents: a computational study

The oximes 4-carbamoyl-1-[({2-[(E)-(hydroxyimino) methyl] pyridinium-1-yl} methoxy) methyl] pyridinium (known as HI-6) and 3-carbamoyl-1-[({2-[(E)-(hydroxyimino) methyl] pyridinium-1-yl} methoxy) methyl] pyridinium (known as HS-6) are isomers differing from each other only by the position of the carbamoyl group on the pyridine ring. However, this slight difference was verified to be responsible for big differences in the percentual of reactivation of acetylcholinesterase (AChE) inhibited by the nerve agents tabun, sarin, cyclosarin, and VX. In order to try to find out the reason for this, a computational study involving molecular docking, molecular dynamics, and binding energies calculations, was performed on the binding modes of HI-6 and HS-6 on human AChE (HssAChE) inhibited by those nerve agents.     

A molecular dynamics study of components of the ginger (Zingiber officinale) extract inside human acetylcholinesterase: implications for Alzheimer disease

Components of ginger (Zingiber officinale) extracts have been described as potential new drug candidates against Alzheimer disease (AD), able to interact with several molecular targets related to the AD treatment. However, there are very few theoretical studies in the literature on the possible mechanisms of action by which these compounds can work as potential anti-AD drugs. For this reason, we performed here docking, molecular dynamic simulations and mmpbsa calculations on four components of ginger extracts former reported as active inhibitors of human acetylcholinesterase (HssAChE), and compared our results to the known HssAChE inhibitor and commercial drug in use against AD, donepezil (DNP). Our findings points to two among the compounds studied: (E)-1,7-bis(4-hydroxy-3-methoxyphenyl)hept-4-en-3-on and 1-(3,4-dihydroxy-5-methoxyphenyl)-7-(4-hydroxy-3- ethoxyphenyl) heptane-3,5-diyl diacetate, as promising new HssAChE inhibitors that could be as effective as DNP. We also mapped the binding of the studied compounds in the different binding pockets inside HssAChE and established the preferred interactions to be favored in the design of new and more efficient inhibitors.     

EVIDENCE FOR AN EXTREME BOTTLENECK IN A RARE MEXICAN PINYON: GENETIC DIVERSITY,DISEQUILIBRIUM, AND THE MATING SYSTEM IN PINUS MAXIMARTINEZII

Maxipiñon (Pinus maximartinezii Rzedowski), which is confined to a single population of approximately 2000 to 2500 mature trees, covers about 400 ha in southern Zacatecas, Mexico. Genetic diversity measured by expected heterozygosity was 0.122, which is moderate for pines. However, percentage polymorphic loci was low, 30.3%. The fixation index (F) of 0.081 indicated only slight heterozygote deficiency. Mating system analysis indicated a significant but low level of selling; the multilocus outcrossing rate, t_m, was 0.816. The mean of single locus estimates, t_s, was smaller (0.761), perhaps suggesting mating among relatives, although the difference between t_m and t_s was not statistically significant. The most striking features of maxipiñon's genetic structure were that no polymorphic locus had more than two alleles and most alleles at polymorphic loci were at intermediate frequencies. This is in contrast to other pines, which often have three to five or more alleles at some loci and in which the distribution of allele frequencies is U-shaped, most alleles being present at frequencies less than 10% or greater than 90%. A population with only two alleles per locus and at intermediate frequencies could occur if the population had been reduced to an extreme bottleneck and then expanded rapidly before random drift modified allele frequencies. A novel origin from a hybridization event would also explain the results. Significant gametic disequilibrium was detected at several pairs of loci in both maternal and paternal gametes. The presence of disequilibrium is in agreement with an origin from an extreme bottleneck, perhaps even a single seed. Furthermore, it demands that the event be relatively recent. The number of generations, as calculated from the observed mean disequilibrium, suggested that maxipiñon derived from an extreme bottleneck four to five generations ago, which is less than 1000 years in this species.     

Behavior of uncharged oximes compared to HI6 and 2-PAM in the human AChE-tabun conjugate: a molecular modeling approach

Tabun is one of the most dangerous nerve agents because it has deleterious effects like inhibition of the essential enzymes acetylcholinesterase (AChE) and butyrylcholinesterase. Some oximes such HI6 as 2-PAM are nucleophiles that are capable to reactivate inhibited human AChE under some conditions. Zwitterionic and cationic species have the best chance of productive action on inhibited AChE. However uncharged oximes can give important interaction information. In order to investigate the interaction and behavior of cationic and uncharged oximes, we performed molecular docking simulations and molecular dynamics and calculated binding energies of complexes of these compounds with human AChE. The uncharged oximes of larger structure were more susceptible to the influence of the substituents on the phosphorus atom and presented low binding energies. In contrast, HI 6 and 2-PAM showed high binding energy values with great contribution of the amino acid Asp74, demonstrating the importance of the quaternary nitrogen to the affinity and interaction of the oximes/AChE tabun-inhibited complexes.     

Dynamics of heme complexed with human serum albumin: a theoretical approach

Human serum albumin (HSA) is the most abundant protein in the blood serum. It binds several ligands and has an especially strong affinity for heme, hence becoming a natural candidate for oxygen transport. In order to analyze the interaction of HSA-heme, molecular dynamics simulations of HSA with bound heme were performed. Based on the results of X-ray diffraction, the binding site of the heme, localized in subdomain IB, was considered. We analyzed the fluctuations and their correlations along trajectories to detect collective motions. The role of H bonds and salt bridges in the stabilization of heme in its pocket was also investigated. Complementarily, the localization of water molecules in the hydrophobic pocket and the interaction with heme were discussed.