Molecular dynamics simulation studies for DNA sequence recognition by reactive metabolites of anticancer compounds

The discovery of novel anticancer molecules 5F‐203 (NSC703786) and 5‐aminoflavone (5‐AMF, NSC686288) has addressed the issues of toxicity and reduced efficacy by targeting over expressed Cytochrome P450 1A1 (CYP1A1) in cancer cells. CYP1A1 metabolizes these compounds into their reactive metabolites, which are proven to mediate their anticancer effect through DNA adduct formation. However, the drug metabolite–DNA binding has not been explored so far. Hence, understanding the binding characteristics and molecular recognition for drug metabolites with DNA is of practical and fundamental interest. The present study is aimed to model binding preference shown by reactive metabolites of 5F‐203 and 5‐AMF with DNA in forming DNA adducts. To perform this, three different DNA crystal structures covering sequence diversity were selected, and 12 DNA‐reactive metabolite complexes were generated. Molecular dynamics simulations for all complexes were performed using AMBER 11 software after development of protocol for DNA‐reactive metabolite system. Furthermore, the MM‐PBSA/GBSA energy calculation, per‐nucleotide energy decomposition, and Molecular Electrostatic Surface Potential analysis were performed. The results obtained from present study clearly indicate that minor groove in DNA is preferable for binding of reactive metabolites of anticancer compounds. The binding preferences shown by reactive metabolites were also governed by specific nucleotide sequence and distribution of electrostatic charges in major and minor groove of DNA structure. Overall, our study provides useful insights into the initial step of mechanism of reactive metabolite binding to the DNA and the guidelines for designing of sequence specific DNA interacting anticancer agents. Copyright © 2014 John Wiley & Sons, Ltd.     

毛竹的异质性空间点格局分析

毛竹是一种重要的经济用材,研究其空间分布格局有益于对其空间发生进行预测。选取了一块空间肥量分布不均的毛竹样地,使用异质性空间点格局分析方法对毛竹的空间分布进行了研究,结果发现毛竹在给定的距离尺度下呈现出显著的泊松分布,而并不是聚集分布。同时选用二元二次方程拟合观察的数据,发现对毛竹强度的拟合优度较高。研究揭示出异质性空间环境中毛竹的分布特点,潜在地表明土壤肥量的空间分布对毛竹空间分布的影响,对毛竹的集约化经营具有重要的指导意义;同时,鉴于国内生态研究中较少采用异质性空间点格局分析方法,研究对于推动这种方法在国内生态学研究中的普及还具有一定的参考价值。     

The contribution of the mitochondrial genome to sex‐specific fitness variance

Maternal inheritance of mitochondrial DNA (mtDNA) facilitates the evolutionary accumulation of mutations with sex‐biased fitness effects. Whereas maternal inheritance closely aligns mtDNA evolution with natural selection in females, it makes it indifferent to evolutionary changes that exclusively benefit males. The constrained response of mtDNA to selection in males can lead to asymmetries in the relative contributions of mitochondrial genes to female versus male fitness variation. Here, we examine the impact of genetic drift and the distribution of fitness effects (DFE) among mutations—including the correlation of mutant fitness effects between the sexes—on mitochondrial genetic variation for fitness. We show how drift, genetic correlations, and skewness of the DFE determine the relative contributions of mitochondrial genes to male versus female fitness variance. When mutant fitness effects are weakly correlated between the sexes, and the effective population size is large, mitochondrial genes should contribute much more to male than to female fitness variance. In contrast, high fitness correlations and small population sizes tend to equalize the contributions of mitochondrial genes to female versus male variance. We discuss implications of these results for the evolution of mitochondrial genome diversity and the genetic architecture of female and male fitness.     

Investigation of PDE5/PDE6 and PDE5/PDE11 selective potent tadalafil-like PDE5 inhibitors using combination of molecular modeling approaches,molecular fingerprint-based virtual screening protocols and structure-based pharmacophore development

The essential biological function of phosphodiesterase (PDE) type enzymes is to regulate the cytoplasmic levels of intracellular second messengers, 3′,5′-cyclic guanosine monophosphate (cGMP) and/or 3′,5′-cyclic adenosine monophosphate (cAMP). PDE targets have 11 isoenzymes. Of these enzymes, PDE5 has attracted a special attention over the years after its recognition as being the target enzyme in treating erectile dysfunction. Due to the amino acid sequence and the secondary structural similarity of PDE6 and PDE11 with the catalytic domain of PDE5, first-generation PDE5 inhibitors (i.e. sildenafil and vardenafil) are also competitive inhibitors of PDE6 and PDE11. Since the major challenge of designing novel PDE5 inhibitors is to decrease their cross-reactivity with PDE6 and PDE11, in this study, we attempt to identify potent tadalafil-like PDE5 inhibitors that have PDE5/PDE6 and PDE5/PDE11 selectivity. For this aim, the similarity-based virtual screening protocol is applied for the “clean drug-like subset of ZINC database” that contains more than 20 million small compounds. Moreover, molecular dynamics (MD) simulations of selected hits complexed with PDE5 and off-targets were performed in order to get insights for structural and dynamical behaviors of the selected molecules as selective PDE5 inhibitors. Since tadalafil blocks hERG1 K channels in concentration dependent manner, the cardiotoxicity prediction of the hit molecules was also tested. Results of this study can be useful for designing of novel, safe and selective PDE5 inhibitors.     

Estimation of uncertainties in X-ray refinement results by use of perturbed structures

The uncertainties in the refined parameters for a 1.5-A X-ray structure of carbon-monoxy (FeII) myoglobin are estimated by combining energy minimization with least-squares refinement against the X-ray data. The energy minimizations, done without reference to the X-ray data, provide perturbed structures which are used to restart conventional X-ray refinement. The resulting refined structures have the same, or better, R-factor and stereochemical parameters as the original X-ray structure, but deviate from it by 0.13 A rms for the backbone atoms and 0.31 A rms for the sidechain atoms. Atoms interacting with a disordered sidechain, Arg 45 CD3, are observed to have larger positional uncertainties. The uncertainty in the B-factors, within the isotropic harmonic motion approximation, is estimated to be 15%. The resulting X-ray structures are more consistent with the energy parameters used in simulations.     

AIPAR: ab initio parametrization of intermolecular potentials for computer simulations

An unambiguous, fully ab initio and automated technique denoted AIPAR (ab initio parametrization) implemented in the SJBR program has been proposed to yield intermolecular interaction potentials between polar molecules and water. The AIPAR procedure has been applied to several organic molecules covering a wide range of structure and functional groups, namely methanol, acetone (propanone), methanethiol (methyl mercaptan), imidazole (1,3-diazole), oxazole and furan. The AIPAR-derived sets of parameters compare well with the empirical OPLS ones, mainly when the all-atoms model is employed in the OPLS procedure. Monte Carlo simulations were performed for an aqueous solution of methanol and for an equimolar binary mixture methanol–water using the AIPAR and OPLS parameters. The thermodynamic and geometric results obtained with the parameters obtained with the AIPAR procedure compare favorably with the OPLS simulations, even for the binary mixture, demonstrating the precision, robustness and transferability of the parameters obtained with the AIPAR procedure.Figure Superimposed configurations of water (without the hydrogen atoms) around the methanol molecule obtained with the AGOA procedure.     

Physical scoring function based on AMBER force field and Poisson-Boltzmann implicit solvent for protein structure prediction

A well-behaved physics-based all-atom scoring function for protein structure prediction is analyzed with several widely used all-atom decoy sets. The scoring function, termed AMBER/Poisson-Boltzmann (PB), is based on a refined AMBER force field for intramolecular interactions and an efficient PB model for solvation interactions. Testing on the chosen decoy sets shows that the scoring function, which is designed to consider detailed chemical environments, is able to consistently discriminate all 62 native crystal structures after considering the heteroatom groups, disulfide bonds, and crystal packing effects that are not included in the decoy structures. When NMR structures are considered in the testing, the scoring function is able to discriminate 8 out of 10 targets. In the more challenging test of selecting near-native structures, the scoring function also performs very well: for the majority of the targets studied, the scoring function is able to select decoys that are close to the corresponding native structures as evaluated by ranking numbers and backbone Calpha root mean square deviations. Various important components of the scoring function are also studied to understand their discriminative contributions toward the rankings of native and near-native structures. It is found that neither the nonpolar solvation energy as modeled by the surface area model nor a higher protein dielectric constant improves its discriminative power. The terms remaining to be improved are related to 1-4 interactions. The most troublesome term is found to be the large and highly fluctuating 1-4 electrostatics term, not the dihedral-angle term. These data support ongoing efforts in the community to develop protein structure prediction methods with physics-based potentials that are competitive with knowledge-based potentials.     

Proton microbeam radiotherapy with scanned pencil-beams – Monte Carlo simulations

Irradiation, delivered by a synchrotron facility, using a set of highly collimated, narrow and parallel photon beams spaced by 1 mm or less, has been termed Microbeam Radiation Therapy (MRT). The tolerance of healthy tissue after MRT was found to be better than after standard broad X-ray beams, together with a more pronounced response of malignant tissue. The microbeam spacing and transverse peak-to-valley dose ratio (PVDR) are considered to be relevant biological MRT parameters. We investigated the MRT concept for proton microbeams, where we expected different depth-dose profiles and PVDR dependences, resulting in skin sparing and homogeneous dose distributions at larger beam depths, due to differences between interactions of proton and photon beams in tissue. Using the FLUKA Monte Carlo code we simulated PVDR distributions for differently spaced 0.1 mm (sigma) pencil-beams of entrance energies 60, 80, 100 and 120 MeV irradiating a cylindrical water phantom with and without a bone layer, representing human head. We calculated PVDR distributions and evaluated uniformity of target irradiation at distal beam ranges of 60–120 MeV microbeams. We also calculated PVDR distributions for a 60 MeV spread-out Bragg peak microbeam configuration. Application of optimised proton MRT in terms of spot size, pencil-beam distribution, entrance beam energy, multiport irradiation, combined with relevant radiobiological investigations, could pave the way for hypofractionation scenarios where tissue sparing at the entrance, better malignant tissue response and better dose conformity of target volume irradiation could be achieved, compared with present proton beam radiotherapy configurations.     

In silico investigation of PARP-1 catalytic domains in holo and apo states for the design of high-affinity PARP-1 inhibitors

The rational design of high-affinity inhibitors of poly-ADP-ribose polymerase-1 (PARP-1) is at the heart of modern anti-cancer drug design. While relevance of enzyme to DNA repair processes in cellular environment is firmly established, the structural and functional understanding of the main determinants for high-affinity ligands controlling PARP-1 activity is still lacking. The conserved active site of PARP-1 represents an ideal target for inhibitors and may offer a novel target at the treatment of breast cancer. To fill the gap in the structural knowledge, we report on the combination of molecular dynamics (MD) simulations, principal component analysis (PCA), and conformational analysis that analyzes in great details novel binding mode for a number of inhibitors at the PARP-1. While optimization of the binding affinity for original target is an important goal in the drug design, many of the promising molecules for treatment of the breast cancer are plagued by significant cardiotoxicity. One of the most common side-effects reported for a number of polymerase inhibitors is its off-target interactions with cardiac ion channels and hERG1 channel, in particular. Thus, selected candidate PARP-1 inhibitors were also screened in silico at the central cavities of hERG1 potassium ion channel.     

Prerequisite for His4 in deltorphin A for highδ opioid receptor selectivity

Summary Analysis of deltorphin A position 4 analogues included: backbone constrained N MeHis, spinacine (Spi), N MePhe and the tetrahydroisoquinoline-3-carboxylic acid (Tic); spatially confined side-chain (Phg); and imidazole alkylation ofl- andd-His⁴ enantiomers. High selectivity was lost with the following replacements: N MeHis⁴, N MePhe⁴ and Phg⁴ reduced binding and the constrained residues also increasedµ binding; ring closure between the side-chain and amino group to yield Spi⁴ or Tic⁴ increasedµ affinity. Imidazole methylation of His⁴ marginally affected opioid binding and doubled selectivity; alkylatedd-His⁴-derivatives generally maintained selectivity in spite of decreased affinities. Thus, His⁴ imidazole preserves selectivity by facilitating high binding and by repulsion at theµ receptor. Several low energy conformers of deltorphin A indicated that the His⁴ imidazole preferred a spatial orientation parallel to the phenolic side-chain of Tyr¹ suggestive that this conformation might contribute to high affinity and selectivity.