scholarly journals Sequence-dependent Three Interaction Site (TIS) Model for Single and Double-stranded DNA

2018 ◽  
Author(s):  
Debayan Chakraborty ◽  
Naoto Hori ◽  
D. Thirumalai

AbstractWe develop a robust coarse-grained model for single and double stranded DNA by representing each nucleotide by three interaction sites (TIS) located at the centers of mass of sugar, phosphate, and base. The resulting TIS model includes base-stacking, hydrogen bond, and electrostatic interactions as well as bond-stretching and bond angle potentials that account for the polymeric nature of DNA. The choices of force constants for stretching and the bending potentials were guided by a Boltzmann inversion procedure using a large representative set of DNA structures extracted from the Protein Data Bank. Some of the parameters in the stacking interactions were calculated using a learning procedure, which ensured that the experimentally measured melting temperatures of dimers are faithfully reproduced. Without any further adjustments, the calculations based on the TIS model reproduces the experimentally measured salt and sequence dependence of the size of single stranded DNA (ssDNA), as well as the persistence lengths of poly(dA) and poly(dT) chains. Interestingly, upon application of mechanical force the extension of poly(dA) exhibits a plateau, which we trace to the formation of stacked helical domains. In contrast, the force-extension curve (FEC) of poly(dT) is entropic in origin, and could be described by a standard polymer model. We also show that the persistence length of double stranded DNA, formed from two complementary ssDNAs with one hundred and thirty base pairs, is consistent with the prediction based on the worm-like chain. The persistence length, which decreases with increasing salt concentration, is in accord with the Odijk-Skolnick-Fixman theory intended for stiff polyelectrolyte chains near the rod limit. The range of applications, which did not require adjusting any parameter after the initial construction based solely on PDB structures and melting profiles of dimers, attests to the transferability and robustness of the TIS model for ssDNA and dsDNA.

2021 ◽  
Author(s):  
Salvatore Assenza ◽  
Rubén Pérez

AbstractWe introduce MADna, a sequence-dependent coarse-grained model of double-stranded DNA (dsDNA), where each nucleotide is described by three beads localized at the sugar and base moieties, and at the phosphate group. The sequence dependence is included by considering a step-dependent parameterization of the bonded interactions, which are tuned in order to reproduce the values of key observables obtained from exhaustive atomistic simulations from literature. The predictions of the model are benchmarked against an independent set of all-atom simulations, showing that it captures with high fidelity the sequence dependence of conformational and elastic features beyond the single step considered in its formulation. A remarkably good agreement with experiments is found for both sequence-averaged and sequence-dependent conformational and elastic features, including the stretching and torsion moduli, the twist-stretch and twist-bend couplings, the persistence length and the helical pitch. Overall, for the inspected quantities, the model has a precision comparable to atomistic simulations, hence providing a reliable coarse-grained description for the rationalization of singlemolecule experiments and the study of cellular processes involving dsDNA. Owing to the simplicity of its formulation, MADna can be straightforwardly included in common simulation engines.


2021 ◽  
Author(s):  
Sidney Carvalho ◽  
Ralf Metzler ◽  
Andrey Cherstvy ◽  
Daniel Caetano

Several applications arise from the confinement of proteins on surfaces since their stability and biological activity are enhanced. It is also known that the way a protein adsorbs on the surface is important for its biological function since its active sites should not be obstructed. In this study, the adsorption properties of hen egg-white Lysozyme, HEWL, into a negatively charged silica pore is examined employing a coarse-grained model and constant-pH Monte Carlo simulations. The role of electrostatic interactions is taken into account when including the Debye-Hueckel potentials into the Ca structure-based model. We evaluate the effects of pH, salt concentration, and pore radius on the protein preferential orientation and spatial distribution of its residues regarding the pore surface. By mapping the residues that stay closer to the pore surface, we find the increase of pH leads to orientational changes of the adsorbed protein when the solution pH gets closer to the HEWL isoelectric point. At these conditions, the pKa shift of these important residues caused by the adsorption into the charged confining surface results in a HEWL charge distribution that stabilizes the adsorption in the observed protein orientation. We compare our observations to the results of pKa shift for HEWL available in the literature and to some experimental data.


2019 ◽  
Vol 9 (3) ◽  
pp. 20180066 ◽  
Author(s):  
S. Pasquali ◽  
E. Frezza ◽  
F. L. Barroso da Silva

Electrostatic interactions play a pivotal role in many biomolecular processes. The molecular organization and function in biological systems are largely determined by these interactions. Owing to the highly negative charge of RNA, the effect is expected to be more pronounced in this system. Moreover, RNA base pairing is dependent on the charge of the base, giving rise to alternative secondary and tertiary structures. The equilibrium between uncharged and charged bases is regulated by the solution pH, which is therefore a key environmental condition influencing the molecule’s structure and behaviour. By means of constant-pH Monte Carlo simulations based on a fast proton titration scheme, coupled with the coarse-grained model HiRE-RNA, molecular dynamic simulations of RNA molecules at constant pH enable us to explore the RNA conformational plasticity at different pH values as well as to compute electrostatic properties as local p K a values for each nucleotide.


2016 ◽  
Vol 72 (3) ◽  
pp. 324-337 ◽  
Author(s):  
A. Janner

Considered is the coarse-grained modeling of icosahedral viruses in terms of a three-dimensional lattice (the digital modeling lattice) selected among the projected points in space of a six-dimensional icosahedral lattice. Backbone atomic positions (Cα's for the residues of the capsid and phosphorus atoms P for the genome nucleotides) are then indexed by their nearest lattice point. This leads to a fine-grained lattice point characterization of the full viral chains in the backbone approximation (denoted as digital modeling). Coarse-grained models then follow by a proper selection of the indexed backbone positions, where for each chain one can choose the desired coarseness. This approach is applied to three viruses, the Satellite tobacco mosaic virus, the bacteriophage MS2 and the Pariacoto virus, on the basis of structural data from the Brookhaven Protein Data Bank. In each case the various stages of the procedure are illustrated for a given coarse-grained model and the corresponding indexed positions are listed. Alternative coarse-grained models have been derived and compared. Comments on related results and approaches, found among the very large set of publications in this field, conclude this article.


2020 ◽  
Author(s):  
Raju Lunkad ◽  
Anastasiia Murmiliuk ◽  
Pascal Hebbeker ◽  
Milan Boublík ◽  
Zdeněk Tošner ◽  
...  

Weak ampholytes are ubiquitous in nature and commonly found in artificial pH-responsive systems. However, our limited understanding of their ionisation response and the lack of predictive capabilities hinder the bottom-up design of such systems. Here, we used a coarse-grained model of a flexible polymer with weakly ionisable monomer units to quantitatively analyse the ionisation behaviour of two oligopeptides. Differences in ionisation response between oligopeptides and monomeric amino acids showed that electrostatic interactions between weak acid and base side chains play a key role in oligopeptide ionisation, as predicted by our model. Moreover, by comparing our simulations with experimental results from potentiometric titration, capillary zone electrophoresis and NMR, we demonstrated that our model reliably predicts the ionisation response and electrophoretic mobilities of various peptide sequences. Ultimately, our model is the first step towards using predictive bottom-up design of responsive ampholytes to tailor their properties as a function of charge and pH.<br>


Soft Matter ◽  
2018 ◽  
Vol 14 (22) ◽  
pp. 4530-4540 ◽  
Author(s):  
Tyler J. Weyer ◽  
Alan R. Denton

Swelling and structural properties of ionic microgel suspensions are described by a coarse-grained model that includes elastic and electrostatic interactions.


2020 ◽  
Author(s):  
Raju Lunkad ◽  
Anastasiia Murmiliuk ◽  
Pascal Hebbeker ◽  
Milan Boublík ◽  
Zdeněk Tošner ◽  
...  

Weak ampholytes are ubiquitous in nature and commonly found in artificial pH-responsive systems. However, our limited understanding of their ionisation response and the lack of predictive capabilities hinder the bottom-up design of such systems. Here, we used a coarse-grained model of a flexible polymer with weakly ionisable monomer units to quantitatively analyse the ionisation behaviour of two oligopeptides. Differences in ionisation response between oligopeptides and monomeric amino acids showed that electrostatic interactions between weak acid and base side chains play a key role in oligopeptide ionisation, as predicted by our model. Moreover, by comparing our simulations with experimental results from potentiometric titration, capillary zone electrophoresis and NMR, we demonstrated that our model reliably predicts the ionisation response and electrophoretic mobilities of various peptide sequences. Ultimately, our model is the first step towards using predictive bottom-up design of responsive ampholytes to tailor their properties as a function of charge and pH.<br>


2017 ◽  
Author(s):  
Olivier Mailhot ◽  
Vincent Frappier ◽  
François Major ◽  
Rafael Najmanovich

ABSTRACTMotivationThe use of Normal Mode Analysis (NMA) methods to study both protein and nucleic acid dynamics is well established. However, the most widely used coarse-grained methods are based on backbone geometry alone and do not take into account the chemical nature of the residues. Elastic Network Contact Model (ENCoM) is a coarse-grained NMA method that includes a pairwise atom-type non-bonded interaction term, which makes it sensitive to the sequence of the studied molecule. We adapted ENCoM to simulate the dynamics of ribonucleic acid (RNA) molecules.ResultsENCoM outperforms the most commonly used coarse-grained model on RNA, Anisotropic Network Model (ANM), in the prediction of b-factors, in the prediction of conformational change as measured by overlap (a measure of effective prediction of structural transitions) and in the prediction of structural variance from NMR ensembles. These benchmarks were derived from the set of all RNA structures available from the Protein Data Bank (PDB) and contain more total cases than previous studies applying NMA to RNA. We thus established ENCoM as an attractive tool for fast and accurate exploration of the conformational space of RNA molecules.AvailabilityENCoM is open source software available at https://github.com/NRGlab/ENCoM


2019 ◽  
Vol 5 (6) ◽  
pp. eaav1697 ◽  
Author(s):  
Min Ju Shon ◽  
Sang-Hyun Rah ◽  
Tae-Young Yoon

Submicrometer elasticity of double-stranded DNA (dsDNA) governs nanoscale bending of DNA segments and their interactions with proteins. Single-molecule force spectroscopy, including magnetic tweezers (MTs), is an important tool for studying DNA mechanics. However, its application to short DNAs under 1 μm is limited. We developed an MT-based method for precise force-extension measurements in the 100-nm regime that enables in situ correction of the error in DNA extension measurement, and normalizes the force variability across beads by exploiting DNA hairpins. The method reduces the lower limit of tractable dsDNA length down to 198 base pairs (bp) (67 nm), an order-of-magnitude improvement compared to conventional tweezing experiments. Applying this method and the finite worm-like chain model we observed an essentially constant persistence length across the chain lengths studied (198 bp to 10 kbp), which steeply depended on GC content and methylation. This finding suggests a potential sequence-dependent mechanism for short-DNA elasticity.


Sign in / Sign up

Export Citation Format

Share Document