Modulation of DNA loop lifetimes by the free energy of loop formation.
Ontology highlight
ABSTRACT:
SUBMITTER: Chen YJ
PROVIDER: S-EPMC4267329 | biostudies-literature | 2014 Dec
REPOSITORIES: biostudies-literature
Publications
Modulation of DNA loop lifetimes by the free energy of loop formation.
Chen Yi-Ju YJ Johnson Stephanie S Mulligan Peter P Spakowitz Andrew J AJ Phillips Rob R
Proceedings of the National Academy of Sciences of the United States of America 20141119 49
Storage and retrieval of the genetic information in cells is a dynamic process that requires the DNA to undergo dramatic structural rearrangements. DNA looping is a prominent example of such a structural rearrangement that is essential for transcriptional regulation in both prokaryotes and eukaryotes, and the speed of such regulations affects the fitness of individuals. Here, we examine the in vitro looping dynamics of the classic Lac repressor gene-regulatory motif. We show that both loop assoc ...[more]
Similar Datasets
Project description:DNA/DNA duplex formation is the basic mechanism that is used in genome tiling arrays and SNP arrays manufactured by Affymetrix. However, detailed knowledge of the physical process is still lacking. In this study, we show a free energy analysis of DNA/DNA duplex formation these arrays based on the positional-dependent nearest-neighbor (PDNN) model, which was developed previously for describing DNA/RNA duplex formation on expression microarrays. Our results showed that the two ends of a probe contribute less to the stability of the duplexes and that there is a microarray surface effect on binding affinities. We also showed that free energy cost of a single mismatch depends on the bases adjacent to the mismatch site and obtained a comprehensive table of the cost of a single mismatch under all possible combination of adjacent bases. The mismatch costs were found to be correlated with those determined in aqueous solution. We further demonstrate that the DNA copy number estimated from the SNP array correlates negatively with the target length; this is presumably caused by inefficient PCR amplification for long fragments. These results provide important insights into the molecular mechanisms of microarray technology and have implications for microarray design and the interpretation of observed data.
| S-EPMC1807971 | biostudies-literature
Project description:The lengths of G-tracts and their connecting loop sequences determine G-quadruplex folding and stability. Complete understanding of the sequence-structure relationships remains elusive. Here, single-loop G-quadruplexes were investigated using explicit solvent molecular dynamics (MD) simulations to characterize the effect of loop length, loop sequence, and G-tract length on the folding topologies and stability of G-quadruplexes. Eight loop types, including different variants of lateral, diagonal, and propeller loops, and six different loop sequences [d0 (i.e., no intervening residues in the loop), dT, dT(2), dT(3), dTTA, and dT(4)] were considered through MD simulation and free energy analysis. In most cases the free energetic estimates agree well with the experimental observations. The work also provides new insight into G-quadruplex folding and stability. This includes reporting the observed instability of the left propeller loop, which extends the rules for G-quadruplex folding. We also suggest a plausible explanation why human telomere sequences predominantly form hybrid-I and hybrid-II type structures in K(+) solution. Overall, our calculation results indicate that short loops generally are less stable than longer loops, and we hypothesize that the extreme stability of sequences with very short loops could possibly derive from the formation of parallel multimers. The results suggest that free energy differences, estimated from MD and free energy analysis with current force fields and simulation protocols, are able to complement experiment and to help dissect and explain loop sequence, loop length, and G-tract length and orientation influences on G-quadruplex structure.
| S-EPMC3168932 | biostudies-literature
Project description:Limitations in protein homology modeling often arise from the inability to adequately model loops. In this paper we focus on the selection of loop conformations. We present a complete computational treatment that allows the screening of loop conformations to identify those that best fit a molecular model. The stability of a loop in a protein is evaluated via computations of conformational free energies in solution, i.e., the free energy difference between the reference structure and the modeled one. A thermodynamic cycle is used for calculation of the conformational free energy, in which the total free energy of the reference state (i.e., gas phase) is the CHARMm potential energy. The electrostatic contribution of the solvation free energy is obtained from solving the finite-difference Poisson-Boltzmann equation. The nonpolar contribution is based on a surface area-based expression. We applied this computational scheme to a simple but well-characterized system, the antibody hypervariable loop (complementarity-determining region, CDR). Instead of creating loop conformations, we generated a database of loops extracted from high-resolution crystal structures of proteins, which display geometrical similarities with antibody CDRs. We inserted loops from our database into a framework of an antibody; then we calculated the conformational free energies of each loop. Results show that we successfully identified loops with a "reference-like" CDR geometry, with the lowest conformational free energy in gas phase only. Surprisingly, the solvation energy term plays a confusing role, sometimes discriminating "reference-like" CDR geometry and many times allowing "non-reference-like" conformations to have the lowest conformational free energies (for short loops). Most "reference-like" loop conformations are separated from others by a gap in the gas phase conformational free energy scale. Naturally, loops from antibody molecules are found to be the best models for long CDRs (> or = 6 residues), mainly because of a better packing of backbone atoms into the framework of the antibody model.
| S-EPMC1181139 | biostudies-other
Project description:Circular permutations usually retain the native structure and function of a protein while inevitably perturbing its folding dynamics. By using simulations with a structure-based model and a rigorous methodology to determine free-energy surfaces from trajectories, we evaluate the effect of a circular permutation on the free-energy landscape of the protein T4 lysozyme. We observe changes which, although subtle, largely affect the cooperativity between the two subdomains. Such a change in cooperativity has been previously experimentally observed and recently also characterized using single molecule optical tweezers and the Crooks relation. The free-energy landscapes show that both the wild type and circular permutant have an on-pathway intermediate, previously experimentally characterized, in which one of the subdomains is completely formed. The landscapes, however, differ in the position of the rate-limiting step for folding, which occurs before the intermediate in the wild type and after in the circular permutant. This shift of transition state explains the observed change in the cooperativity. The underlying free-energy landscape thus provides a microscopic description of the folding dynamics and the connection between circular permutation and the loss of cooperativity experimentally observed.
| S-EPMC3821731 | biostudies-literature
Project description:MOTIVATION:Predicting the secondary structure of RNA is a fundamental task in bioinformatics. Algorithms that predict secondary structure given only the primary sequence, and a model to evaluate the quality of a structure, are an integral part of this. These algorithms have been updated as our model of RNA thermodynamics changed and expanded. An exception to this has been the treatment of multi-loops. Although more advanced models of multi-loop free energy change have been suggested, a simple, linear model has been used since the 1980s. However, recently, new dynamic programing algorithms for secondary structure prediction that could incorporate these models were presented. Unfortunately, these models appear to have lower accuracy for secondary structure prediction. RESULTS:We apply linear regression and a new parameter optimization algorithm to find better parameters for the existing linear model and advanced non-linear multi-loop models. These include the Jacobson-Stockmayer and Aalberts & Nandagopal models. We find that the current linear model parameters may be near optimal for the linear model, and that no advanced model performs better than the existing linear model parameters even after parameter optimization. AVAILABILITY AND IMPLEMENTATION:Source code and data is available at https://github.com/maxhwardg/advanced_multiloops. SUPPLEMENTARY INFORMATION:Supplementary data are available at Bioinformatics online.
| S-EPMC6821347 | biostudies-literature
Project description:Estimation of the energy from a given Boltzmann sample is straightforward since one just has to average the contribution of the individual configurations. On the other hand, calculation of the absolute entropy, S (hence the absolute free energy F) is difficult because it depends on the entire (unknown) ensemble. We have developed a new method called "the hypothetical scanning molecular dynamics" (HSMD) for calculating the absolute S from a given sample (generated by any simulation technique). In other words, S (like the energy) is "written" on the sample configurations, where HSMD provides a prescription of how to "read" it. In practice, each sample conformation, i, is reconstructed with transition probabilities, and their product leads to the probability of i, hence to the entropy. HSMD is an exact method where all interactions are considered, and the only approximation is due to insufficient sampling. In previous studies HSMD (and HS Monte CarloHSMC) has been extended systematically to systems of increasing complexity, where the most recent is the seven-residue mobile loop, 304-310 (Gly-His-Gly-Ala-Gly-Gly-Ser) of the enzyme porcine pancreatic alpha-amylase modeled by the AMBER force field and AMBER with the implicit solvation GB/SA (paper I, Cheluvaraja, S.; Meirovitch, H. J. Chem. Theory Comput. 2008, 4, 192). In the present paper we make a step further and extend HSMD to the same loop capped with TIP3P explicit water at 300 K. As in paper I, we are mainly interested in entropy and free energy differences between the free and bound microstates of the loop, which are obtained from two separate MD samples of these microstates. The contribution of the loop to S and F is calculated by HSMD and that of water by a particular thermodynamic integration procedure. As expected, the free microstate is more stable than the bound microstate by a total free energy difference, Ffree-Fbound=-4.8+/-1, as compared to -25.5 kcal/mol obtained with GB/SA. We find that relatively large systematic errors in the loop entropies, Sfree(loop) and Sbound(loop) are cancelled in their difference which is thus obtained efficiently and with high accuracy, i.e., with a statistical error of 0.1 kcal/mol. This cancellation, which has been observed in previous HSMD studies, is in accord with theoretical arguments given in paper I.
| S-EPMC2671085 | biostudies-literature
Project description:For the first time, an expansive study into the concentration and extended decay behavior of environmentally persistent free radicals in PM2.5 was performed. Results from this study revealed three types of radical decay-a fast decay, slow decay, and no decay-following one of four decay patterns: a relatively fast decay exhibiting a 1/e lifetime of 1-21 days accompanied by a slow decay with a 1/e lifetime of 21-5028 days (47% of samples); a single slow decay including a 1/e lifetime of 4-2083 days (24% of samples); no decay (18% of samples); and a relatively fast decay displaying an average 1/e lifetime of 0.25-21 days followed by no decay (11% of samples). Phenol correlated well with the initial radical concentration and fast decay rate. Other correlations for common atmospheric pollutants (ozone, NOx, SO2, etc.) as well as meteorological conditions suggested photochemical processes impact the initial radical concentration and fast decay rate. The radical signal in PM2.5 was remarkably similar to semiquinones in cigarette smoke. Accordingly, radicals inhaled from PM2.5 were related to the radicals inhaled from smoking cigarettes, expressed as the number of equivalent cigarettes smoked. This calculated to 0.4-0.9 cigarettes per day for nonextreme air quality in the United States.
| S-EPMC4309365 | biostudies-literature
Project description:BackgroundEstimation of DNA duplex hybridization free energy is widely used for predicting cross-hybridizations in DNA computing and microarray experiments. A number of software programs based on different methods and parametrizations are available for the theoretical estimation of duplex free energies. However, significant differences in free energy values are sometimes observed among estimations obtained with various methods, thus being difficult to decide what value is the accurate one.ResultsWe present in this study a quantitative comparison of the similarities and differences among four published DNA/DNA duplex free energy calculation methods and an extended Nearest-Neighbour Model for perfect matches based on triplet interactions. The comparison was performed on a benchmark data set with 695 pairs of short oligos that we collected and manually curated from 29 publications. Sequence lengths range from 4 to 30 nucleotides and span a large GC-content percentage range. For perfect matches, we propose an extension of the Nearest-Neighbour Model that matches or exceeds the performance of the existing ones, both in terms of correlations and root mean squared errors. The proposed model was trained on experimental data with temperature, sodium and sequence concentration characteristics that span a wide range of values, thus conferring the model a higher power of generalization when used for free energy estimations of DNA duplexes under non-standard experimental conditions.ConclusionsBased on our preliminary results, we conclude that no statistically significant differences exist among free energy approximations obtained with 4 publicly available and widely used programs, when benchmarked against a collection of 695 pairs of short oligos collected and curated by the authors of this work based on 29 publications. The extended Nearest-Neighbour Model based on triplet interactions presented in this work is capable of performing accurate estimations of free energies for perfect match duplexes under both standard and non-standard experimental conditions and may serve as a baseline for further developments in this area of research.
| S-EPMC2837027 | biostudies-literature
Project description:The Cre-recombination system has become an important tool for genetic manipulation of higher organisms and a model for site-specific DNA-recombination mechanisms employed by the λ-Int superfamily of recombinases. We report a novel quantitative approach for characterizing the probability of DNA-loop formation in solution using time-dependent ensemble Förster resonance energy transfer measurements of intra- and inter-molecular Cre-recombination kinetics. Our method uses an innovative technique for incorporating multiple covalent modifications at specific sites in covalently closed DNA. Because the mechanism of Cre recombinase does not conform to a simple kinetic scheme, we employ numerical methods to extract rate constants for fundamental steps that pertain to Cre-mediated loop closure. Cre recombination does not require accessory proteins, DNA supercoiling or particular metal-ion cofactors and is thus a highly flexible system for quantitatively analyzing DNA-loop formation in vitro and in vivo.
| S-EPMC3424569 | biostudies-literature
Project description:The loop 287-290 (Ile, Phe, Arg, and Phe) of the protein acetylcholinesterase (AChE) changes its structure upon interaction of AChE with diisopropylphosphorofluoridate (DFP). Reversible dissociation measurements suggest that the free-energy (F) penalty for the loop displacement is DeltaF=Ffree-Fbound approximately -4 kcal/mol. Therefore, this loop has been the target of two studies by Olson's group for testing the efficiency of procedures for calculating F. In this paper, we test for the first time the performance of our "hypothetical scanning molecular dynamics" (HSMD) method and the validity of the related modeling for a loop with bulky side chains in explicit water. Thus, we consider only atoms of the protein that are the closest to the loop (they constitute the "template"), where the rest of the atoms are ignored. The template's atoms are fixed in the X-ray coordinates of the free protein, and the loop is capped with a sphere of TIP3P water molecules; also, the X-ray structure of the bound loop is attached to the free template. We carry out two separate MD simulations starting from the free and bound X-ray structures, where only the atoms of the loop and water are allowed to move while the template-water and template-loop (AMBER) interactions are considered. The absolute Ffree and Fbound (of the loop and water) are calculated from the corresponding trajectories. A main objective of this paper is to assess the reliability of this model, and for this several template sizes are studied capped with 80-220 water molecules. We find that consistent results for the free energy (which also agree with the experimental data above) require a template larger than a minimal size and a number of water molecules approximately equal to the experimental density of bulk water. For example, we obtain DeltaFtotal=DeltaFwater+DeltaFloop=-3.1+/-2.5 and -3.6+/-4 kcal/mol for a template consisting of 944 atoms and a sphere containing 160 and 180 waters, respectively. Our calculations demonstrate the important contribution of water to the total free energy. Namely, for water densities close to the experimental value, DeltaFwater is always negative leading thereby to a negative DeltaFtotal (while DeltaFloop is always positive). Also, the contribution of the water entropy TDeltaSwater to DeltaFtotal is significant. Various aspects related to the efficiency of HSMD are tested and improved, and plans for future studies are discussed.
| S-EPMC2747743 | biostudies-literature