Intrinsic Conformational Energetics Associated with the Glycosyl Torsion in DNA: A Quantum Mechanical Study

AbstractThe glycosyl torsion (χ) in nucleic acids has long been recognized to be a major determinant of their conformational properties. χ torsional energetics were systematically mapped in deoxyribonucleosides using high-level quantum mechanical methods, for north and south sugar puckers and with γ in the g+ and trans conformations. In all cases, the syn conformation is found higher in energy than the anti. When γ is changed from g+ to trans, the anti orientation of the base is strongly destabilized, and the energy difference and barrier between anti and syn are significantly decreased. The barrier between anti and syn in deoxyribonucleosides is found to be less than 10kcal/mol and tends to be lower with purines than with pyrimidines. With γ=g+/χ=anti, a south sugar yields a significantly broader energy well than a north sugar with no energy barrier between χ values typical of A or B DNA. Contrary to the prevailing view, the syn orientation is not more stable with south puckers than with north puckers. The syn conformation is significantly more energetically accessible with guanine than with adenine in 5-nucleotides but not in nucleosides. Analysis of nucleic acid crystal structures shows that γ=trans/χ=anti is a minor but not negligible conformation. Overall, χ appears to be a very malleable structural parameter with the experimental χ distributions reflecting, to a large extent, the associated intrinsic torsional energetics

Context-Dependent Cell Cycle Checkpoint Abrogation by a Novel Kinase Inhibitor

Checkpoint kinase 1 and 2 (Chk1/Chk2), and the Aurora kinases play a critical role in the activation of the DNA damage response and mitotic spindle checkpoints. We have identified a novel inhibitor of these kinases and utilized this molecule to probe the functional interplay between these two checkpoints.Fragment screening, structure guided design, and kinase cross screening resulted in the identification of a novel, potent small molecule kinase inhibitor (VER-150548) of Chk1 and Chk2 kinases with IC(50)s of 35 and 34 nM as well as the Aurora A and Aurora B kinases with IC(50)s of 101 and 38 nM. The structural rationale for this kinase specificity could be clearly elucidated through the X-ray crystal structure. In human carcinoma cells, VER-150548 induced reduplication and the accumulation of cells with >4N DNA content, inhibited histone H3 phosphorylation and ultimately gave way to cell death after 120 hour exposure; a phenotype consistent with cellular Aurora inhibition. In the presence of DNA damage induced by cytotoxic chemotherapeutic drugs, VER-150548 abrogated DNA damage induced cell cycle checkpoints. Abrogation of these checkpoints correlated with increased DNA damage and rapid cell death in p53 defective HT29 cells. In the presence of DNA damage, reduplication could not be observed. These observations are consistent with the Chk1 and Chk2 inhibitory activity of this molecule.In the presence of DNA damage, we suggest that VER-150548 abrogates the DNA damage induced checkpoints forcing cells to undergo a lethal mitosis. The timing of this premature cell death induced by Chk1 inhibition negates Aurora inhibition thereby preventing re-entry into the cell cycle and subsequent DNA reduplication. This novel kinase inhibitor therefore serves as a useful chemical probe to further understand the temporal relationship between cell cycle checkpoint pathways, chemotherapeutic agent induced DNA damage and cell death

rDock: A Fast, Versatile and Open Source Program for Docking Ligands to Proteins and Nucleic Acids

Identification of chemical compounds with specific biological activities is an important step in both chemical biology and drug discovery. When the structure of the intended target is available, one approach is to use molecular docking programs to assess the chemical complementarity of small molecules with the target; such calculations provide a qualitative measure of affinity that can be used in virtual screening (VS) to rank order a list of compounds according to their potential to be active. rDock is a molecular docking program developed at Vernalis for high-throughput VS (HTVS) applications. Evolved from RiboDock, the program can be used against proteins and nucleic acids, is designed to be computationally very efficient and allows the user to incorporate additional constraints and information as a bias to guide docking. This article provides an overview of the program structure and features and compares rDock to two reference programs, AutoDock Vina (open source) and Schrodinger's Glide (commercial). In terms of computational speed for VS, rDock is faster than Vina and comparable to Glide. For binding mode prediction, rDock and Vina are superior to Glide. The VS performance of rDock is significantly better than Vina, but inferior to Glide for most systems unless pharmacophore constraints are used; in that case rDock and Glide are of equal performance. The program is released under the Lesser General Public License and is freely available for download, together with the manuals, example files and the complete test sets, at http://rdock.sourceforge.net

Intrinsic flexibility of B-DNA: the experimental TRX scale

B-DNA flexibility, crucial for DNA–protein recognition, is sequence dependent. Free DNA in solution would in principle be the best reference state to uncover the relation between base sequences and their intrinsic flexibility; however, this has long been hampered by a lack of suitable experimental data. We investigated this relationship by compiling and analyzing a large dataset of NMR 31P chemical shifts in solution. These measurements reflect the BI ↔ BII equilibrium in DNA, intimately correlated to helicoidal descriptors of the curvature, winding and groove dimensions. Comparing the ten complementary DNA dinucleotide steps indicates that some steps are much more flexible than others. This malleability is primarily controlled at the dinucleotide level, modulated by the tetranucleotide environment. Our analyses provide an experimental scale called TRX that quantifies the intrinsic flexibility of the ten dinucleotide steps in terms of Twist, Roll, and X-disp (base pair displacement). Applying the TRX scale to DNA sequences optimized for nucleosome formation reveals a 10 base-pair periodic alternation of stiff and flexible regions. Thus, DNA flexibility captured by the TRX scale is relevant to nucleosome formation, suggesting that this scale may be of general interest to better understand protein-DNA recognition

DNA structures from phosphate chemical shifts

For B-DNA, the strong linear correlation observed by nuclear magnetic resonance (NMR) between the 31P chemical shifts (δP) and three recurrent internucleotide distances demonstrates the tight coupling between phosphate motions and helicoidal parameters. It allows to translate δP into distance restraints directly exploitable in structural refinement. It even provides a new method for refining DNA oligomers with restraints exclusively inferred from δP. Combined with molecular dynamics in explicit solvent, these restraints lead to a structural and dynamical view of the DNA as detailed as that obtained with conventional and more extensive restraints. Tests with the Jun-Fos oligomer show that this δP-based strategy can provide a simple and straightforward method to capture DNA properties in solution, from routine NMR experiments on unlabeled samples

How Thioredoxin Dissociates Its Mixed Disulfide

The dissociation mechanism of the thioredoxin (Trx) mixed disulfide complexes is unknown and has been debated for more than twenty years. Specifically, opposing arguments for the activation of the nucleophilic cysteine as a thiolate during the dissociation of the complex have been put forward. As a key model, the complex between Trx and its endogenous substrate, arsenate reductase (ArsC), was used. In this structure, a Cys29Trx-Cys89ArsC intermediate disulfide is formed by the nucleophilic attack of Cys29Trx on the exposed Cys82ArsC-Cys89ArsC in oxidized ArsC. With theoretical reactivity analysis, molecular dynamics simulations, and biochemical complex formation experiments with Cys-mutants, Trx mixed disulfide dissociation was studied. We observed that the conformational changes around the intermediate disulfide bring Cys32Trx in contact with Cys29Trx. Cys32Trx is activated for its nucleophilic attack by hydrogen bonds, and Cys32Trx is found to be more reactive than Cys82ArsC. Additionally, Cys32Trx directs its nucleophilic attack on the more susceptible Cys29Trx and not on Cys89ArsC. This multidisciplinary approach provides fresh insights into a universal thiol/disulfide exchange reaction mechanism that results in reduced substrate and oxidized Trx

Understanding the Sequence-Dependence of DNA Groove Dimensions: Implications for DNA Interactions

BACKGROUND: The B-DNA major and minor groove dimensions are crucial for DNA-protein interactions. It has long been thought that the groove dimensions depend on the DNA sequence, however this relationship has remained elusive. Here, our aim is to elucidate how the DNA sequence intrinsically shapes the grooves. METHODOLOGY/PRINCIPAL FINDINGS: The present study is based on the analysis of datasets of free and protein-bound DNA crystal structures, and from a compilation of NMR (31)P chemical shifts measured on free DNA in solution on a broad range of representative sequences. The (31)P chemical shifts can be interpreted in terms of the BI↔BII backbone conformations and dynamics. The grooves width and depth of free and protein-bound DNA are found to be clearly related to the BI/BII backbone conformational states. The DNA propensity to undergo BI↔BII backbone transitions is highly sequence-dependent and can be quantified at the dinucleotide level. This dual relationship, between DNA sequence and backbone behavior on one hand, and backbone behavior and groove dimensions on the other hand, allows to decipher the link between DNA sequence and groove dimensions. It also firmly establishes that proteins take advantage of the intrinsic DNA groove properties. CONCLUSIONS/SIGNIFICANCE: The study provides a general framework explaining how the DNA sequence shapes the groove dimensions in free and protein-bound DNA, with far-reaching implications for DNA-protein indirect readout in both specific and non specific interactions

Structure-Guided Evolution of Potent and Selective CHK1 Inhibitors through Scaffold Morphing

Pyrazolopyridine inhibitors with low micromolar potency for CHK1 and good selectivity against CHK2 were previously identified by fragment-based screening. The optimization of the pyrazolopyridines to a series of potent and CHK1-selective isoquinolines demonstrates how fragment-growing and scaffold morphing strategies arising from a structure-based understanding of CHK1 inhibitor binding can be combined to successfully progress fragment-derived hit matter to compounds with activity in vivo. The challenges of improving CHK1 potency and selectivity, addressing synthetic tractability, and achieving novelty in the crowded kinase inhibitor chemical space were tackled by multiple scaffold morphing steps, which progressed through tricyclic pyrimido[2,3-b]azaindoles to N-(pyrazin-2-yl)pyrimidin-4-amines and ultimately to imidazo[4,5-c]pyridines and isoquinolines. A potent and highly selective isoquinoline CHK1 inhibitor (SAR-020106) was identified, which potentiated the efficacies of irinotecan and gemcitabine in SW620 human colon carcinoma xenografts in nude mice

Computational Study of the Human Dystrophin Repeats: Interaction Properties and Molecular Dynamics

Dystrophin is a large protein involved in the rare genetic disease Duchenne muscular dystrophy (DMD). It functions as a mechanical linker between the cytoskeleton and the sarcolemma, and is able to resist shear stresses during muscle activity. In all, 75% of the dystrophin molecule consists of a large central rod domain made up of 24 repeat units that share high structural homology with spectrin-like repeats. However, in the absence of any high-resolution structure of these repeats, the molecular basis of dystrophin central domain's functions has not yet been deciphered. In this context, we have performed a computational study of the whole dystrophin central rod domain based on the rational homology modeling of successive and overlapping tandem repeats and the analysis of their surface properties. Each tandem repeat has very specific surface properties that make it unique. However, the repeats share enough electrostatic-surface similarities to be grouped into four separate clusters. Molecular dynamics simulations of four representative tandem repeats reveal specific flexibility or bending properties depending on the repeat sequence. We thus suggest that the dystrophin central rod domain is constituted of seven biologically relevant sub-domains. Our results provide evidence for the role of the dystrophin central rod domain as a scaffold platform with a wide range of surface features and biophysical properties allowing it to interact with its various known partners such as proteins and membrane lipids. This new integrative view is strongly supported by the previous experimental works that investigated the isolated domains and the observed heterogeneity of the severity of dystrophin related pathologies, especially Becker muscular dystrophy