Dataset Information

Relative stability of different DNA guanine quadruplex stem topologies derived using large-scale quantum-chemical computations.

ABSTRACT: We provide theoretical predictions of the intrinsic stability of different arrangements of guanine quadruplex (G-DNA) stems. Most computational studies of nucleic acids have applied Molecular Mechanics (MM) approaches using simple pairwise-additive force fields. The principle limitation of such calculations is the highly approximate nature of the force fields. In this study, we for the first time apply accurate QM computations (DFT-D3 with large atomic orbital basis sets) to essentially complete DNA building blocks, seven different folds of the cation-stabilized two-quartet G-DNA stem, each having more than 250 atoms. The solvent effects are approximated by COSMO continuum solvent. We reveal sizable differences between MM and QM descriptions of relative energies of different G-DNA stems, which apparently reflect approximations of the DNA force field. Using the QM energy data, we propose correction to earlier free energy estimates of relative stabilities of different parallel, hybrid, and antiparallel G-stem folds based on classical simulations. The new energy ranking visibly improves the agreement between theory and experiment. We predict the 5'-anti-anti-3' GpG dinucleotide step to be the most stable one, closely followed by the 5'-syn-anti-3' step. The results are in good agreement with known experimental structures of 2-, 3-, and 4-quartet G-DNA stems. Besides providing specific results for G-DNA, our study highlights basic limitations of force field modeling of nucleic acids. Although QM computations have their own limitations, mainly the lack of conformational sampling and the approximate description of the solvent, they can substantially improve the quality of calculations currently relying exclusively on force fields.

SUBMITTER: Sponer J

PROVIDER: S-EPMC3775466 | biostudies-literature | 2013 Jul

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

Relative stability of different DNA guanine quadruplex stem topologies derived using large-scale quantum-chemical computations.

Šponer Jiří J Mládek Arnošt A Špačková Naďa N Cang Xiaohui X Cheatham Thomas E TE Grimme Stefan S

Journal of the American Chemical Society 20130619 26

We provide theoretical predictions of the intrinsic stability of different arrangements of guanine quadruplex (G-DNA) stems. Most computational studies of nucleic acids have applied Molecular Mechanics (MM) approaches using simple pairwise-additive force fields. The principle limitation of such calculations is the highly approximate nature of the force fields. In this study, we for the first time apply accurate QM computations (DFT-D3 with large atomic orbital basis sets) to essentially complete ...[more]

PMID: 23742743

Similar Datasets

Project description:Two major techniques have been used to determine the three-dimensional structures of proteins: X-ray diffraction and NMR spectroscopy. In particular, the validation of NMR-derived protein structures is one of the most challenging problems in NMR spectroscopy. Therefore, researchers have proposed a plethora of methods to determine the accuracy and reliability of protein structures. Despite these proposals, there is a growing need for more sophisticated, physics-based structure validation methods. This approach will enable us to (a) characterize the "quality" of the NMR-derived ensemble as a whole by a single parameter, (b) unambiguously identify flaws in the sequence at a residue level, and (c) provide precise information, such as sets of backbone and side-chain torsional angles, that we can use to detect local flaws. Rather than reviewing all of the existing validation methods, this Account describes the contributions of our research group toward a solution of the long-standing problem of both global and local structure validation of NMR-derived protein structures. We emphasize a recently introduced physics-based methodology that makes use of observed and computed (13)C(alpha) chemical shifts (at the density functional theory (DFT) level of theory) for an accurate validation of protein structures in solution and in crystals. By assessing the ability of computed (13)C(alpha) chemical shifts to reproduce observed (13)C(alpha) chemical shifts of a single structure or ensemble of structures in solution and in crystals, we accomplish a global validation by using the conformationally averaged root-mean-square deviation, ca-rmsd, as a scoring function. In addition, the method enables us to provide local validation by identifying a set of individual amino acid conformations for which the computed and observed (13)C(alpha) chemical shifts do not agree within a certain error range and may represent a nonreliable fold of the protein model. Although it is computationally intensive, our validation method has several advantages, which we illustrate through a series of applications. This method makes use of the (13)C(alpha) chemical shifts, not shielding, that are ubiquitous to proteins and can be computed precisely from the phi, psi, and chi torsional angles. There is no need for a priori knowledge of the oligomeric state of the protein, and no knowledge-based information or additional NMR data are required. The primary limitation at this point is the computational cost of such calculations. However, we anticipate that enhancements in the speed of calculating these chemical shifts coupled with the ever-increasing computational power should soon make this a standard method accessible to the general NMR community.

Project description:Many real-world situations require navigating decisions for both reward and threat. While there has been significant progress in understanding mechanisms of decision-making and mediating neurocircuitry separately for reward and threat, there is limited understanding of situations where reward and threat contingencies compete to create approach-avoidance conflict (AAC). Here, we leverage computational learning models, independent component analysis (ICA), and multivariate pattern analysis (MVPA) approaches to understand decision-making during a novel task that embeds concurrent reward and threat learning and manipulates congruency between reward and threat probabilities. Computational modeling supported a modified reinforcement learning model where participants integrated reward and threat value into a combined total value according to an individually varying policy parameter, which was highly predictive of decisions to approach reward vs avoid threat during trials where the highest reward option was also the highest threat option (i.e., approach-avoidance conflict). ICA analyses demonstrated unique roles for salience, frontoparietal, medial prefrontal, and inferior frontal networks in differential encoding of reward vs threat prediction error and value signals. The left frontoparietal network uniquely encoded degree of conflict between reward and threat value at the time of choice. MVPA demonstrated that delivery of reward and threat could accurately be decoded within salience and inferior frontal networks, respectively, and that decisions to approach reward vs avoid threat were predicted by the relative degree to which these reward vs threat representations were active at the time of choice. This latter result suggests that navigating AAC decisions involves generating mental representations for possible decision outcomes, and relative activation of these representations may bias subsequent decision-making towards approaching reward or avoiding threat accordingly.

Dataset Information

Relative stability of different DNA guanine quadruplex stem topologies derived using large-scale quantum-chemical computations.

Publications

Relative stability of different DNA guanine quadruplex stem topologies derived using large-scale quantum-chemical computations.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets