Browse
Submit Data
Databases
API
Help

Dataset Information

0 Views

0 Connections

0 Citations

0 Reanalyses

0 Downloads

Omics score: 0

Exploring the counterion atmosphere around DNA: what can be learned from molecular dynamics simulations?

ABSTRACT: The counterion distribution around a DNA dodecamer (5'-CGCGAATTCGCG-3') is analyzed using both standard and novel techniques based on state of the art molecular dynamics simulations. Specifically, we have explored the population of Na(+) in the minor groove of DNA duplex, and whether or not a string of Na(+) can replace the spine of hydration in the narrow AATT minor groove. The results suggest that the insertion of Na(+) in the minor groove is a very rare event, but that when once the ion finds specific sites deep inside the groove it can reside there for very long periods of time. According to our simulation the presence of Na(+) inside the groove does not have a dramatic influence in the structure or dynamics of the duplex DNA. The ability of current MD simulations to obtain equilibrated pictures of the counterion atmosphere around DNA is critically discussed.

SUBMITTER: Rueda M

PROVIDER: S-EPMC1304490 | biostudies-literature | 2004 Aug

REPOSITORIES: biostudies-literature

ACCESS DATA

Json Xml

Publications

Exploring the counterion atmosphere around DNA: what can be learned from molecular dynamics simulations?

Rueda Manuel M Cubero Elena E Laughton Charles A CA Orozco Modesto M

Biophysical journal 20040801 2

The counterion distribution around a DNA dodecamer (5'-CGCGAATTCGCG-3') is analyzed using both standard and novel techniques based on state of the art molecular dynamics simulations. Specifically, we have explored the population of Na(+) in the minor groove of DNA duplex, and whether or not a string of Na(+) can replace the spine of hydration in the narrow AATT minor groove. The results suggest that the insertion of Na(+) in the minor groove is a very rare event, but that when once the ion finds ...[more]

PMID: 15298889

Similar Datasets

Towards exact molecular dynamics simulations with machine-learned force fields.

Project description:Molecular dynamics (MD) simulations employing classical force fields constitute the cornerstone of contemporary atomistic modeling in chemistry, biology, and materials science. However, the predictive power of these simulations is only as good as the underlying interatomic potential. Classical potentials often fail to faithfully capture key quantum effects in molecules and materials. Here we enable the direct construction of flexible molecular force fields from high-level ab initio calculations by incorporating spatial and temporal physical symmetries into a gradient-domain machine learning (sGDML) model in an automatic data-driven way. The developed sGDML approach faithfully reproduces global force fields at quantum-chemical CCSD(T) level of accuracy and allows converged molecular dynamics simulations with fully quantized electrons and nuclei. We present MD simulations, for flexible molecules with up to a few dozen atoms and provide insights into the dynamical behavior of these molecules. Our approach provides the key missing ingredient for achieving spectroscopic accuracy in molecular simulations.

| S-EPMC6155327 | biostudies-literature

Exploring kainate receptor pharmacology using molecular dynamics simulations.

Project description:Ionotropic glutamate receptors (iGluRs) are enticing targets for pharmaceutical research; however, the search for selective ligands is a laborious experimental process. Here we introduce a purely computational procedure as an approach to evaluate ligand-iGluR pharmacology. The ligands are docked into the closed ligand-binding domain and during the molecular dynamics (MD) simulation the bi-lobed interface either opens (partial agonist/antagonist) or stays closed (agonist) according to the properties of the ligand. The procedure is tested with closely related set of analogs of the marine toxin dysiherbaine bound to GluK1 kainate receptor. The modeling is set against the abundant binding data and electrophysiological analyses to test reproducibility and predictive value of the procedure. The MD simulations produce detailed binding modes for analogs, which in turn are used to define structure-activity relationships. The simulations suggest correctly that majority of the analogs induce full domain closure (agonists) but also distinguish exceptions generated by partial agonists and antagonists. Moreover, we report ligand-induced opening of the GluK1 ligand-binding domain in free MD simulations. The strong correlation between in silico analysis and the experimental data imply that MD simulations can be utilized as a predictive tool for iGluR pharmacology and functional classification of ligands.

| S-EPMC2876335 | biostudies-literature

Exploring Structures and Dynamics of Protamine Molecules through Molecular Dynamics Simulations.

Project description:Protamines are arginine-rich proteins that condense DNA in sperm. Despite their importance in reproduction, information on protamine structure is scarce. We, therefore, used molecular dynamics to examine the structures of salmon, bull P1, and human P1 protamines. The sizes and shapes of each protamine varied widely, indicating that they were disordered with structures covering a broad conformational landscape, from hairpin loop structures to extended coils. Despite their general disorder, the protamines did form secondary structures, including helices and hairpin loops. In eutherians, hairpins may promote disulfide bonding that facilitates protamine-DNA condensation, but the specifics of this bonding is not well established. We examined inter-residue distances in the simulations to predict residue pairs likely to form intramolecular bonds, leading to the identification of bonding pairs consistent with previous results in bull and human. These results support a model for eutherian protamine structures where a highly charged center is surrounded by disulfide-bond-stabilized loops.

| S-EPMC9685783 | biostudies-literature

Detecting Counterion Dynamics in DNA-Protein Association.

Project description:Due to a high density of negative charges on its surface, DNA condenses cations as counterions, forming the so-called "ion atmosphere". Although the release of counterions upon DNA-protein association has been postulated to have a major contribution to the binding thermodynamics, this release remains to be confirmed through a direct observation of the ions. Herein, we report the characterization of the ion atmosphere around DNA using NMR spectroscopy and directly detect the release of counterions upon DNA-protein association. NMR-based diffusion data reveal the highly dynamic nature of counterions within the ion atmosphere around DNA. Counterion release is observed as an increase in the apparent ionic diffusion coefficient, which directly provides the number of counterions released upon DNA-protein association.

| S-EPMC6980997 | biostudies-literature

Molecular dynamics simulations of p53 DNA-binding domain.

Project description:We have studied room-temperature structural and dynamic properties of the p53 DNA-binding domain in both DNA-bound and DNA-free states. A cumulative 55 ns of explicit solvent molecular dynamics simulations with the particle mesh Ewald treatment of electrostatics was performed. It was found that the mean structures in the production portions of the trajectories agree well with the crystal structure: backbone root-mean-square deviations are in the range of 1.6 and 2.0 A. In both simulations, noticeable backbone deviations from the crystal structure are observed only in loop L6, due to the lack of crystal packing in the simulations. More deviations are observed in the DNA-free simulation, apparently due to the absence of DNA. Computed backbone B-factor is also in qualitative agreement with the crystal structure. Interestingly, little backbone structural change is observed between the mean simulated DNA-bound and DNA-free structures. A notable difference is observed only at the DNA-binding interface. The correlation between native contacts and inactivation mechanisms of tumor mutations is also discussed. In the H2 region, tumor mutations at sites D281, R282, E285, and E286 may weaken five key interactions that stabilize H2, indicating that their inactivation mechanisms may be related to the loss of local structure around H2, which in turn may reduce the overall stability to a measurable amount. In the L2 region, tumor mutations at sites Y163, K164, E171, V173, L194, R249, I251, and E271 are likely to be responsible for the loss of stability in the protein. In addition to apparent DNA contacts that are related to DNA binding, interactions R175/S183, S183/R196, and E198/N235 are highly occupied only in the DNA-bound form, indicating that they are more likely to be responsible for DNA binding.

| S-EPMC2522240 | biostudies-literature

Molecular dynamics simulations of DNA-polycation complex formation.

Project description:Complexes formed from DNA and polycations are of interest because of their potential use in gene therapy; however, there remains a lack of understanding of the structure and formation of DNA-polycation complexes at atomic scale. In this work, molecular dynamics simulations of the DNA duplex d(CGCGAATTCGCG) in the presence of polycation chains are carried out to shed light on the specific atomic interaction that result in complex formation. The structures of complexes formed from DNA with polyethylenimine, which is considered one of the most promising DNA vector candidates, and a second polycation, poly-L-lysine, are compared. After an initial separation of approximately 50 A, the DNA and polycation come together and form a stable complex within 10 ns. The DNA does not undergo any major structural changes on complexation and remains in the B-form. In the formed complex, the charged amine groups of the polycation mainly interact with DNA phosphate groups, with polycation intrusion into the major and minor grooves dependent on the identity and charge state of the polycation. The ability of the polycation to effectively neutralize the charge of the DNA phosphate groups and the resulting influence on the DNA helix interaction are discussed.

| S-EPMC2756393 | biostudies-literature

Exploring TRPC3 Interaction with Cholesterol through Coarse-Grained Molecular Dynamics Simulations.

Project description:Transient receptor potential canonical 3 (TRPC3) channel belongs to the superfamily of transient receptor potential (TRP) channels which mediate Ca2+ influx into the cell. These channels constitute essential elements of cellular signalling and have been implicated in a wide range of diseases. TRPC3 is primarily gated by lipids and its surface expression has been shown to be dependent on cholesterol, yet a comprehensive exploration of its interaction with this lipid has thus far not emerged. Here, through 80 µs of coarse-grained molecular dynamics simulations, we show that cholesterol interacts with multiple elements of the transmembrane machinery of TRPC3. Through our approach, we identify an annular binding site for cholesterol on the pre-S1 helix and a non-annular site at the interface between the voltage-sensor-like domain and pore domains. Here, cholesterol interacts with exposed polar residues and possibly acts to stabilise the domain interface.

| S-EPMC9313397 | biostudies-literature

Molecular dynamics simulations of DNA-free and DNA-bound TAL effectors.

Project description:TAL (transcriptional activator-like) effectors (TALEs) are DNA-binding proteins, containing a modular central domain that recognizes specific DNA sequences. Recently, the crystallographic studies of TALEs revealed the structure of DNA-recognition domain. In this article, molecular dynamics (MD) simulations are employed to study two crystal structures of an 11.5-repeat TALE, in the presence and absence of DNA, respectively. The simulated results indicate that the specific binding of RVDs (repeat-variable diresidues) with DNA leads to the markedly reduced fluctuations of tandem repeats, especially at the two ends. In the DNA-bound TALE system, the base-specific interaction is formed mainly by the residue at position 13 within a TAL repeat. Tandem repeats with weak RVDs are unfavorable for the TALE-DNA binding. These observations are consistent with experimental studies. By using principal component analysis (PCA), the dominant motions are open-close movements between the two ends of the superhelical structure in both DNA-free and DNA-bound TALE systems. The open-close movements are found to be critical for the recognition and binding of TALE-DNA based on the analysis of free energy landscape (FEL). The conformational analysis of DNA indicates that the 5' end of DNA target sequence has more remarkable structural deformability than the other sites. Meanwhile, the conformational change of DNA is likely associated with the specific interaction of TALE-DNA. We further suggest that the arrangement of N-terminal repeats with strong RVDs may help in the design of efficient TALEs. This study provides some new insights into the understanding of the TALE-DNA recognition mechanism.

| S-EPMC3794935 | biostudies-literature

Exploring thermophysical properties of CoCrFeNiCu high entropy alloy via molecular dynamics simulations.

Project description:High entropy alloys (HEAs) are alloys composed of five or more primary elements in equal or nearly equal proportions of atoms. In the present study, the thermophysical properties of the CoCrFeNiCu high entropy alloy (HEA) were investigated by a molecular dynamics (MD) method at nanoscale. The effects of the content of individual elements on lattice thermal conductivity k p were revealed, and the results suggested that adjusting the atomic content can be a way to control the lattice thermal conductivity of HEAs. The effects of temperature on k p were investigated quantitively, and a power-law relationship of k p with T -0.419 was suggested, which agrees with previous findings. The effects of temperature and the content of individual elements on volumetric specific heat capacity C v were also studied: as the temperature increases, the C v of all HEAs slightly decreases and then increases. The effects of atomic content on C v varied with the comprising elements. To further understand heat transfer mechanisms in the HEAs, the phonon density of states (PDOS) at different temperatures and varying atomic composition was calculated: Co and Ni elements facilitate the high-frequency vibration of phonons and the Cu environment weakens the heat transfer via low-frequency vibration of photons. As the temperature increases, the phonon mean free path (MFP) in the equiatomic CoCrFeNiCu HEA decreases, which may be attributed to the accelerated momentum of atoms and intensified collisions of phonons. The present research provides theoretical foundations for alloy design and have implications for high-performance alloy smelting.

| S-EPMC11369437 | biostudies-literature

DNA basepair step deformability inferred from molecular dynamics simulations.

Project description:The sequence-dependent DNA deformability at the basepair step level was investigated using large-scale atomic resolution molecular dynamics simulation of two 18-bp DNA oligomers: d(GCCTATAAACGCCTATAA) and d(CTAGGTGGATGACTCATT). From an analysis of the structural fluctuations, the harmonic potential energy functions for all 10 unique steps with respect to the six step parameters have been evaluated. In the case of roll, three distinct groups of steps have been identified: the flexible pyrimidine-purine (YR) steps, intermediate purine-purine (RR), and stiff purine-pyrimidine (RY). The YR steps appear to be the most flexible in tilt and partially in twist. Increasing stiffness from YR through RR to RY was observed for rise, whereas shift and slide lack simple trends. A proposed measure of the relative importance of couplings identifies the slide-rise, twist-roll, and twist-slide couplings to play a major role. The force constants obtained are of similar magnitudes to those based on a crystallographic ensemble. However, the current data have a less complicated and less pronounced sequence dependence. A correlation analysis reveals concerted motions of neighboring steps and thus exposes limitations in the dinucleotide model. The comparison of DNA deformability from this and other studies with recent quantum-chemical stacking energy calculations suggests poor correlation between the stacking and flexibility.

| S-EPMC1303568 | biostudies-literature

OmicsDI is part of the ELIXIR infrastructure

OmicsDI is an Elixir interoperability service. Learn more ›

Tweets

OmicsDI Databases

PRIDE
PeptideAtlas
MassIVE
JPOST Repository
Physiome Model Repository

EGA
EVA
ENA
LINCS
PAXDB
Cell Collective

MetaboLights
Metabolomics Workbench
MetabolomeExpress
GNPS
BioModels
FAIRDOMHub

ArrayExpress
dbGaP
ExpressionAtlas
GEO
NODE

Information

Databases
Help
API
Contact us
Code on GitHub
Terms of use
Submit Data