PubMed Health⌕ Search

Biomedical subjects

Jirí Sponer

Publications and source records attributed to Jirí Sponer.

35 records · Page 2Linked to original sources
PubMed2004

Theoretical study of the guanine --> 6-thioguanine substitution in duplexes, triplexes, and tetraplexes.

Molecular dynamics and thermodynamic integration calculations have been carried out on a set of G-rich single-strand, duplex, triplex, and quadruplex DNAs to study the structural and stability changes connected with the guanine --> 6-thioguanine (G --> S) mutation. The presence of 6-thioguanine leads to a shift of the geometry from the B/A intermediate to the pure B-form in duplex DNA. The G --> S mutation does not largely affect the structure of the antiparallel triplex when it is located at the reverse-Hoogsteen position, but leads to a non-negligible local distortion in the structure when it is located at the Watson-Crick position. The G --> S mutation leads to destabilization of all studied structures: the lowest effect has been observed for the G --> S mutation in the reverse-Hoogsteen strand of the triplex, a medium effect has been observed in the Watson-Crick strand of the triplex and duplex, and the highest influence of the G -->S mutation has been found for the quadruplex structures.

DNA↗
PubMed2004

Long-residency hydration, cation binding, and dynamics of loop E/helix IV rRNA-L25 protein complex.

Molecular dynamics simulations of RNA-protein complex between Escherichia coli loop E/helix IV (LE/HeIV) rRNA and L25 protein reveal a qualitative agreement between the experimental and simulated structures. The major groove of LE is a prominent rRNA cation-binding site. Divalent cations rigidify the LE major groove geometry whereas in the absence of divalent cations LE extensively interacts with monovalent cations via inner-shell binding. The HeIV region shows bistability of its major groove explaining the observed differences between x-ray and NMR structures. In agreement with the experiments, the simulations suggest that helix-alpha1 of L25 is the least stable part of the protein. Inclusion of Mg2+ cations into the simulations causes perturbation of basepairing at the LE/HeIV junction, which does not, however, affect the protein binding. The rRNA-protein complex is mediated by a number of highly specific hydration sites with long-residing water molecules and two of them are bound throughout the entire 24-ns simulation. Long-residing water molecules are seen also outside the RNA-protein contact areas with water-binding times substantially enhanced compared to simulations of free RNA. Long-residency hydration sites thus represent important elements of the three-dimensional structure of rRNA.

Binding Sites↗
PubMed2004

Accurate interaction energies of hydrogen-bonded nucleic acid base pairs.

Hydrogen-bonded nucleic acids base pairs substantially contribute to the structure and stability of nucleic acids. The study presents reference ab initio structures and interaction energies of selected base pairs with binding energies ranging from -5 to -47 kcal/mol. The molecular structures are obtained using the RI-MP2 (resolution of identity MP2) method with extended cc-pVTZ basis set of atomic orbitals. The RI-MP2 method provides results essentially identical with the standard MP2 method. The interaction energies are calculated using the Complete Basis Set (CBS) extrapolation at the RI-MP2 level. For some base pairs, Coupled-Cluster corrections with inclusion of noniterative triple contributions (CCSD(T)) are given. The calculations are compared with selected medium quality methods. The PW91 DFT functional with the 6-31G basis set matches well the RI-MP2/CBS absolute interaction energies and reproduces the relative values of base pairing energies with a maximum relative error of 2.6 kcal/mol when applied with Becke3LYP-optimized geometries. The Becke3LYP DFT functional underestimates the interaction energies by few kcal/mol with relative error of 2.2 kcal/mol. Very good performance of nonpolarizable Cornell et al. force field is confirmed and this indirectly supports the view that H-bonded base pairs are primarily stabilized by electrostatic interactions.

Base Pairing↗
PubMed2004

DNA deformability at the base pair level.

A complete set of harmonic force constants describing the DNA deformation energetics at the base pair level was obtained using unrestrained atomic-resolution molecular dynamics simulations of selected duplex oligonucleotides and subsequent analysis of structural fluctuations from the simulated trajectories. The deformation was described by the six base pair conformational parameters (buckle, propeller, opening, shear, stretch, stagger). The results for 13 AT pairs and 11 GC pairs in different sequence contexts suggest that buckle and propeller are very flexible (more than roll in TA dinucleotide steps), while stretch is exceptionally stiff. Only stretch and opening stiffness were found to depend unambiguously on the base pair identity (AT vs GC). The relationship of the results to a simple plates-and-springs model of base-base interactions is discussed.

Base Pairing↗
PubMed2004

Theoretical calculation of the NMR spin-spin coupling constants and the NMR shifts allow distinguishability between the specific direct and the water-mediated binding of a divalent metal cation to guanine.

The calculated intermolecular and intramolecular indirect NMR spin-spin coupling constants and NMR shifts were used for the discrimination between the inner-shell and the outer-shell binding motif of hydrated divalent cations Mg(2+) or Zn(2+) with a guanine base. The intermolecular coupling constants (1)J(X,O6) and (1)J(X,N7) (X = Mg(2+), Zn(2+)) can be unambiguously assigned to the specific inner-shell binding motif of the hydrated cation either with oxygen O6 or with nitrogen N7 of guanine. The calculated coupling constants (1)J(Mg,O6) and (1)J(Zn,O6) were 6.2 and -17.5 Hz, respectively, for the inner-shell complex of cation directly interacting with oxygen O6 of guanine. For the inner-shell coordination of the cation at nitrogen N7, the calculated coupling constants (1)J(Mg,N7) and (1)J(Zn,N7) were 5.6 and -36.5 Hz, respectively. When the binding of the cation is water-mediated, the coupling constant is zero. To obtain reliable shifts in NMR parameters, hydrated guanine was utilized as the reference state. The calculated change of NMR spin-spin coupling constants due to the hydration and coordination of the cation with guanine is caused mainly by the variation of Fermi-contact coupling contribution while the variation of diamagnetic spin-orbit, paramagnetic spin-orbit, and spin-dipolar coupling contributions is small. The change of s-character of guanine sigma bonding, sigma antibonding, and lone pair orbitals upon the hydration and cation coordination (calculated using the Natural Bond Orbital analysis) correlates with the variation of the Fermi-contact term. The calculated NMR shifts delta(N7) of -15.3 and -12.2 ppm upon the coordination of Mg(2+) and Zn(2+) ion are similar to the NMR shift of 19.6 ppm toward the high field measured by Tanaka for N7 of guanine upon the coordination of the Cd(2+) cation (Tanaka, Y.; Kojima, C.; Morita, E. H.; Kasai. Y.; Yamasaki, K.; Ono, A.; Kainosho, M.; Taira, K. J. Am. Chem. Soc. 2002, 124, 4595-4601). The present data indicate that measurements of NMR intermolecular coupling constants may be used to discriminate between the specific inner- and outer-shell binding of divalent cations to nucleobases in DNA and RNA.

Cations, Divalent↗
PubMed2004

Ribosomal RNA kink-turn motif--a flexible molecular hinge.

Ribosomal RNA K-turn motifs are asymmetric internal loops characterized by a sharp bend in the phosphodiester backbone resulting in "V" shaped structures, recurrently observed in ribosomes and showing a high degree of sequence conservation. We have carried out extended explicit solvent molecular dynamics simulations of selected K-turns, in order to investigate their intrinsic structural and dynamical properties. The simulations reveal an unprecedented dynamical flexibility of the K-turns around their X-ray geometries. The K-turns sample, on the nanosecond timescale, different conformational substates. The overall behavior of the simulations suggests that the sampled geometries are essentially isoenergetic and separated by minimal energy barriers. The nanosecond dynamics of isolated K-turns can be qualitatively considered as motion of two rigid helix stems controlled by a very flexible internal loop which then leads to substantial hinge-like motions between the two stems. This internal dynamics of K-turns is strikingly different for example from the bacterial 5S rRNA Loop E motif or BWYV frameshifting pseudoknot which appear to be rigid in the same type of simulations. Bistability and flexibility of K-turns was also suggested by several recent biochemical studies. Although the results of MD simulations should be considered as a qualitative picture of the K-turn dynamics due to force field and sampling limitations, the main advantage of the MD technique is its ability to investigate the region close to K-turn ribosomal-like geometries. This part of the conformational space is not well characterized by the solution experiments due to large-scale conformational changes seen in the experiments. We suggest that K-turns are well suited to act as flexible structural elements of ribosomal RNA. They can for example be involved in mediation of large-scale motions or they can allow a smooth assembling of the other parts of the ribosome.

Base Sequence↗
PubMed2004

Molecular dynamics simulations of Guanine quadruplex loops: advances and force field limitations.

A computational analysis of d(GGGGTTTTGGGG)(2) guanine quadruplexes containing either lateral or diagonal four-thymidine loops was carried out using molecular dynamics (MD) simulations in explicit solvent, locally enhanced sampling (LES) simulations, systematic conformational search, and free energy molecular-mechanics, Poisson Boltzmann, surface area (MM-PBSA) calculations with explicit inclusion of structural monovalent cations. The study provides, within the approximations of the applied all-atom additive force field, a qualitatively complete analysis of the available loop conformational space. The results are independent of the starting structures. Major conformational transitions not seen in conventional MD simulations are observed when LES is applied. The favored LES structures consistently provide lower free energies (as estimated by molecular-mechanics, Poisson Boltzmann, surface area) than other structures. Unfortunately, the predicted optimal structure for the diagonal loop arrangement differs substantially from the atomic resolution experiments. This result is attributed to force field deficiencies, such as the potential misbalance between solute-cation and solvent-cation terms. The MD simulations are unable to maintain the stable coordination of the monovalent cations inside the diagonal loops as reported in a recent x-ray study. The optimal diagonal and lateral loop arrangements appear to be close in energy although a proper inclusion of the loop monovalent cations could stabilize the diagonal architecture.

Computer Simulation↗
PubMed2003

Molecular dynamics simulations of RNA kissing-loop motifs reveal structural dynamics and formation of cation-binding pockets.

Explicit solvent molecular dynamics (MD) simulations were carried out for three RNA kissing-loop complexes. The theoretical structure of two base pairs (2 bp) complex of H3 stem-loop of Moloney murine leukemia virus agrees with the NMR structure with modest violations of few NMR restraints comparable to violations present in the NMR structure. In contrast to the NMR structure, however, MD shows relaxed intermolecular G-C base pairs. The core region of the kissing complex forms a cation-binding pocket with highly negative electrostatic potential. The pocket shows nanosecond-scale breathing motions coupled with oscillations of the whole molecule. Additional simulations were carried out for 6 bp kissing complexes of the DIS HIV-1 subtypes A and B. The simulated structures agree well with the X-ray data. The subtype B forms a novel four-base stack of bulged-out adenines. Both 6 bp kissing complexes have extended cation-binding pockets in their central parts. While the pocket of subtype A interacts with two hexacoordinated Mg2+ ions and one sodium ion, pocket of subtype B is filled with a string of three delocalized Na+ ions with residency times of individual cations 1-2 ns. The 6 bp complexes show breathing motions of the cation-binding pockets and loop major grooves.

Binding Sites↗
PubMed2003

Unique tertiary and neighbor interactions determine conservation patterns of Cis Watson-Crick A/G base-pairs.

X-ray, phylogenetic and quantum chemical analysis of molecular interactions and conservation patterns of cis Watson-Crick (W.C.) A/G base-pairs in 16S rRNA, 23S rRNA and other molecules was carried out. In these base-pairs, the A and G nucleotides interact with their W.C. edges with glycosidic bonds oriented cis relative to each other. The base-pair is stabilised by two hydrogen bonds, the C1'-C1' distance is enlarged and the G(N2) amino group is left unpaired. Quantum chemical calculations show that, in the absence of other interactions, the unpaired amino group is substantially non-planar due to its partial sp(3) pyramidalization, while the whole base-pair is internally propeller twisted and very flexible. The unique molecular properties of the cis W.C. A/G base-pairs make them distinct from other base-pairs. They occur mostly at the ends of canonical helices, where they serve as interfaces between the helix and other motifs. The cis W.C. A/G base-pairs play crucial roles in natural RNA structures with salient sequence conservation patterns. The key contribution to conservation is provided by the unpaired G(N2) amino group that is involved in a wide range of tertiary and neighbor contacts in the crystal structures. Many of them are oriented out of the plane of the guanine base and utilize the partial sp(3) pyramidalization of the G(N2). There is a lack of A/G to G/A covariation, which, except for the G(N2) position, would be entirely isosteric. On the contrary, there is a rather frequent occurrence of G/A to G/U covariation, as the G/U wobble base-pair has an unpaired amino group in the same position as the cis W.C. G/A base-pair. The cis W.C. A/G base-pairs are not conserved when there is no tertiary or neighbor interaction. Obtaining the proper picture of the interactions and phylogenetic patterns of the cis W.C. A/G base-pairs requires a detailed analysis of the relation between the molecular structures and the energetics of interactions at a level of single H-bonds and contacts.

Base Pairing↗
PubMed2003

Molecular dynamics simulations and thermodynamics analysis of DNA-drug complexes. Minor groove binding between 4',6-diamidino-2-phenylindole and DNA duplexes in solution.

An extended set of nanosecond-scale molecular dynamics simulations of DNA duplex sequences in explicit solvent interacting with the minor groove binding drug 4',6-diamidino-2-phenylindole (DAPI) are investigated for four different and sequence specific binding modes. Force fields for DAPI have been parametrized to properly reflect its internal nonplanarity. Sequences investigated include the binding modes observed experimentally, that is, AATT in d(CGCGAATTCGCG)(2) and ATTG in d(GGCCAATTGG)(2) and alternative shifted binding modes ATTC and AATT, respectively. In each case, stable MD simulations are obtained, well reproducing specific hydration patterns seen in the experiments. In contrast to the 2.4 A d(CGCGAATTCGCG)(2) crystal structure, the DAPI is nonplanar, consistent with its gas-phase geometry and the higher resolution crystal structure. The simulations also suggest that the DAPI molecule is able to adopt different conformational substates accompanied by specific hydration patterns that include long-residing waters. The MM_PBSA technology for estimating relative free energies was utilized. The most consistent free energy results were obtained with an approach that uses a single trajectory of the DNA-DAPI complex to estimate all free energy terms. It is demonstrated that explicit inclusion of a subset of bound water molecules shifts the calculated relative binding free energies in favor of both crystallographically observed binding modes, underlining the importance of structured hydration.

Binding Sites↗
PubMed2003

Formation pathways of a guanine-quadruplex DNA revealed by molecular dynamics and thermodynamic analysis of the substates.

The formation of a cation-stabilized guanine quadruplex (G-DNA) stem is an exceptionally slow process involving complex kinetics that has not yet been characterized at atomic resolution. Here, we investigate the formation of a parallel stranded G-DNA stem consisting of four strands of d(GGGG) using molecular dynamics simulations with explicit inclusion of counterions and solvent. Due to the limitations imposed by the nanosecond timescale of the simulations, rather than watching for the spontaneous formation of G-DNA, our approach probes the stability of possible supramolecular intermediates (including two-, three-, and four-stranded assemblies with out-of-register base pairing between guanines) on the formation pathway. The simulations suggest that "cross-like" two-stranded assemblies may serve as nucleation centers in the initial formation of parallel stranded G-DNA quadruplexes, proceeding through a series of rearrangements involving trapping of cations, association of additional strands, and progressive slippage of strands toward the full stem. To supplement the analysis, approximate free energies of the models are obtained with explicit consideration of the integral cations. The approach applied here serves as a prototype for qualitatively investigating other G-DNA molecules using molecular dynamics simulation and free-energy analysis.

Cations↗
PubMed2003

DNA basepair step deformability inferred from molecular dynamics simulations.

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.

Base Pairing↗
PubMed2003

Non-Watson-Crick basepairing and hydration in RNA motifs: molecular dynamics of 5S rRNA loop E.

Explicit solvent and counterion molecular dynamics simulations have been carried out for a total of >80 ns on the bacterial and spinach chloroplast 5S rRNA Loop E motifs. The Loop E sequences form unique duplex architectures composed of seven consecutive non-Watson-Crick basepairs. The starting structure of spinach chloroplast Loop E was modeled using isostericity principles, and the simulations refined the geometries of the three non-Watson-Crick basepairs that differ from the consensus bacterial sequence. The deep groove of Loop E motifs provides unique sites for cation binding. Binding of Mg(2+) rigidifies Loop E and stabilizes its major groove at an intermediate width. In the absence of Mg(2+), the Loop E motifs show an unprecedented degree of inner-shell binding of monovalent cations that, in contrast to Mg(2+), penetrate into the most negative regions inside the deep groove. The spinach chloroplast Loop E shows a marked tendency to compress its deep groove compared with the bacterial consensus. Structures with a narrow deep groove essentially collapse around a string of Na(+) cations with long coordination times. The Loop E non-Watson-Crick basepairing is complemented by highly specific hydration sites ranging from water bridges to hydration pockets hosting 2 to 3 long-residing waters. The ordered hydration is intimately connected with RNA local conformational variations.

Base Pair Mismatch↗
PubMed2002

Toward true DNA base-stacking energies: MP2, CCSD(T), and complete basis set calculations.

Stacking energies in low-energy geometries of pyrimidine, uracil, cytosine, and guanine homodimers were determined by the MP2 and CCSD(T) calculations utilizing a wide range of split-valence, correlation-consistent, and bond-functions basis sets. Complete basis set MP2 (CBS MP2) stacking energies extrapolated using aug-cc-pVXZ (X = D, T, and for pyrimidine dimer Q) basis sets equal to -5.3, -12.3, and -11.2 kcal/mol for the first three dimers, respectively. Higher-order correlation corrections estimated as the difference between MP2 and CCSD(T) stacking energies amount to 2.0, 0.7, and 0.9 kcal/mol and lead to final estimates of the genuine stacking energies for the three dimers of -3.4, -11.6, and -10.4 kcal/mol. The CBS MP2 stacking-energy estimate for guanine dimer (-14.8 kcal/mol) was based on the 6-31G(0.25) and aug-cc-pVDZ calculations. This simplified extrapolation can be routinely used with a meaningful accuracy around 1 kcal/mol for large aromatic stacking clusters. The final estimate of the guanine stacking energy after the CCSD(T) correction amounts to -12.9 kcal/mol. The MP2/6-31G(0.25) method previously used as the standard level to calculate aromatic stacking in hundreds of geometries of nucleobase dimers systematically underestimates the base stacking by ca. 1.0-2.5 kcal/mol per stacked dimer, covering 75-90% of the intermolecular correlation stabilization. We suggest that this correction is to be considered in calibration of force fields and other cheaper computational methods. The quality of the MP2/6-31G(0.25) predictions is nevertheless considerably better than suggested on the basis of monomer polarizability calculations. Fast and very accurate estimates of the MP2 aromatic stacking energies can be achieved using the RI-MP2 method. The CBS MP2 calculations and the CCSD(T) correction, when taken together, bring only marginal changes to the relative stability of H-bonded and stacked base pairs, with a slight shift of ca. 1 kcal/mol in favor of H-bonding. We suggest that the present values are very close to ultimate predictions of the strength of aromatic base stacking of DNA and RNA bases.

Cytosine↗
PubMed2002

Intercalators. 1. Nature of stacking interactions between intercalators (ethidium, daunomycin, ellipticine, and 4',6-diaminide-2-phenylindole) and DNA base pairs. Ab initio quantum chemical, density functional theory, and empirical potential study.

Properties of isolated intercalators (ethidium (E), daunomycin (D), ellipticine (EL), and 4,6'-diaminide-2-phenylindole (DAPI)) and their stacking interactions with adenine...thymine (AT) and guanine...cytosine (GC) nucleic acid base pairs were investigated by means of a nonempirical correlated ab initio method. All intercalators exhibit large charge delocalization, and none of them (including the DAPI dication) exhibits a site with dominant charge. All intercalators have large polarizability and are good electron acceptors, while base pairs are good electron donors. MP2/6-31G*(0.25) stabilization energies of intercalator...base pair complexes are large (E...AT, 22.4 kcal/mol; D...GC, 17.8 kcal/mol; EL...GC, 18.2 kcal/mol; DAPI...GC, 21.1 kcal/mol) and are well reproduced by modified AMBER potential (van der Waals radii of intercalator atoms are enlarged and their energy depths are increased). Standard AMBER potential underestimates binding, especially for DAPI-containing complexes. Because the DAPI dication is the best electron acceptor (among all intercalators studied), this difference is explained by the importance of the charge-transfer term, which is not included in the AMBER potential. For the neutral EL molecule, the standard AMBER force field provides correct results. The Hartree-Fock and DFT/B3LYP methods, not covering the dispersion energy, fail completely to reveal any energy minimum at the potential energy curve of the E...AT complex, and these methods thus cannot be recommended for a study of intercalation process. On the other hand, an approximate version of the DFT method, which was extended to cover London dispersion energy, yields for all complexes very good stabilization energies that are well comparable with referenced ab initio data. Besides the vertical dependence of the interaction, an energy twist dependence of the interaction energy was also investigated by a reference correlated ab initio method and empirical potentials. It is concluded that, despite the cationic (E +1, D +1, DAPI +2) or polar (EL) character of the intercalators investigated, it is the dispersion energy which predominantly contributes to the stability of intercalator...base pair complexes. Any procedure which does not cover dispersion energy is thus not suitable for studying the process of intercalation.

Adenine↗
PubMed2002

Critical effect of the N2 amino group on structure, dynamics, and elasticity of DNA polypurine tracts.

Unrestrained 5-20-ns explicit-solvent molecular dynamics simulations using the Cornell et al. force field have been carried out for d[GCG(N)11GCG]2 (N, purine base) considering guanine*cytosine (G*C), adenine*thymine (A*T), inosine*5-methyl-cytosine (I*mC), and 2-amino-adenine*thymine (D*T) basepairs. The simulations unambiguously show that the structure and elasticity of N-tracts is primarily determined by the presence of the amino group in the minor groove. Simulated A-, I-, and AI-tracts show almost identical structures, with high propeller twist and minor groove narrowing. G- and D-tracts have small propeller twisting and are partly shifted toward the A-form. The elastic properties also differ between the two groups. The sequence-dependent electrostatic component of base stacking seems to play a minor role. Our conclusions are entirely consistent with available experimental data. Nevertheless, the propeller twist and helical twist in the simulated A-tract appear to be underestimated compared to crystallographic studies. To obtain further insight into the possible force field deficiencies, additional multiple simulations have been made for d(A)10, systematically comparing four major force fields currently used in DNA simulations and utilizing B and A-DNA forms as the starting structure. This comparison shows that the conclusions of the present work are not influenced by the force field choice.

Base Pairing↗