PubMed Health⌕ Search

Biomedical subjects

J Sponer

Publications and source records attributed to J Sponer.

At least 19 recordsLinked to original sources
PubMed2001

Molecular dynamics of the frame-shifting pseudoknot from beet western yellows virus: the role of non-Watson-Crick base-pairing, ordered hydration, cation binding and base mutations on stability and unfolding.

Molecular dynamics simulations of the frame-shifting pseudoknot from beet western yellows virus (BWYV, NDB file UR0004) were performed with explicit inclusion of solvent and counterions. In all, 33 ns of simulation were carried out, including 10 ns of the native structure with protonation of the crucial cytosine residue, C8(N3+). The native structure exhibited stable trajectories retaining all Watson-Crick and tertiary base-pairs, except for fluctuations or transient disruptions at specific sites. The most significant fluctuations involved the change or disruption of hydrogen-bonding between C8(N3+) and bases G12, A25, and C26, as well as disruption of the water bridges linking C8(N3+) with A25 and C26. To increase sampling of rare events, the native simulation was continued at 400 K. A partial, irreversible unfolding of the molecule was initiated by slippage of C8(N3+) relative to G12 and continued by sudden concerted changes in hydrogen-bonding involving A23, A24, and A25. These events were followed by a gradual loss of stacking interactions in loop 2. Of the Watson-Crick base-pairs, only the 5'-terminal pair of stem 1 dissociated at 400 K, while the trans sugar-edge/sugar-edge A20.G4 interaction remained surprisingly stable. Four additional room-temperature simulations were carried out to obtain insights into the structural and dynamic effects of selected mutations. In two of these, C8 was left unprotonated. Considerable local rearrangements occurred that were not observed in the crystal structure, thus confirming N3-protonation of C8 in the native molecule. We also investigated the effect of mutating C8(N3+) to U8, to correlate with experimental and phylogenetic studies, and of changing the G4 x C17 base-pair to A4 x U17 to weaken the trans sugar-edge interaction between positions 4 and 20 and to test models of unfolding. The simulations indicate that the C8 x G12 x C26 base-triple at the junction is the most labile region of the frame-shifting pseudoknot. They provide insights into the roles of the other non-Watson-Crick base-pairs in the early stages of unfolding of the pseudoknot, which must occur to allow readthrough of the message by the ribosome. The simulations revealed several critical, highly ordered hydration sites with close to 100 % occupancies and residency times of individual water molecules of up to 5 ns. Sodium cation coordination sites with occupancies above 50 % were also observed.

Base Pairing↗
PubMed2001

Protonation of platinated adenine nucleobases. Gas phase vs condensed phase picture.

Protonation of adenine carrying a Pt(II) moiety either at N7, N3, or N1 is possible in solution, but the site of protonation is influenced by the location of the Pt(II) electrophile and to some extent also by the overall charge of the metal entity (+2, +1, 0, -1), hence the other ligands (NH(3), Cl(-), OH(-)) bound to Pt(II). Quantum chemical calculations based on density functional theory (DFT) have been carried out for intrinsic protonation energies of adenine complexes carrying the following Pt(II) species at either of the three ring N atoms: [Pt(NH(3))(3)](2+) (1), trans- [Pt(NH(3))(2)Cl](+) (2a), cis-[Pt(NH(3))(2)Cl](+) (2b), trans-[Pt(NH(3))(2)Cl(2)] (3a), cis-[Pt(NH(3))Cl(2)] (3b), [PtCl(3)](-) (4), trans-[Pt(NH(3))(2)OH](+) (5a), cis-[Pt(NH(3))(2)(OH)](+) (5b), trans-[Pt(NH(3))(OH)(2)] (6a), cis-[Pt(NH(3))(OH)(2)] (6b), and [Pt(OH)(3)](-) (7). The data have been compared with results derived from solution studies (water) and X-ray crystallography, whenever available. The electrostatic effects associated with the charge of the metal entity have the major influence on the calculated intrinsic (gas phase) proton affinities, unlike the condensed phase data. Nevertheless, the relative gas phase trends correlate surprisingly well with condensed phase data; i.e., variation of the pK(a) values measured in solution is consistent with the calculated gas phase protonation energies. In addition to a systematic study of the ring proton affinities, proton transfer processes within the platinated adenine species were often observed when investigating Pt adducts with OH(-) ligands, and they are discussed in more detail. To the best of our knowledge, this is the first study attempting to find a systematic correlation between gas phase and condensed phase data on protonation of metalated nucleobases. The gas phase data provide a very useful complement to the condensed phase and X-ray experiments, showing that the gas phase studies are capable of valuable predictions and contribute to our understanding of the solvent and counterion effects on metal-assisted proton shift processes.

Adenine Nucleotides↗
PubMed2001

Structural dynamics and cation interactions of DNA quadruplex molecules containing mixed guanine/cytosine quartets revealed by large-scale MD simulations.

Large-scale molecular dynamics (MD) simulations have been utilized to study G-DNA quadruplex molecules containing mixed GCGC and all-guanine GGGG quartet layers. Incorporation of mixed GCGC quartets into G-DNA stems substantially enhances their sequence variability. The mixed quadruplexes form rigid assemblies that require integral monovalent cations for their stabilization. The interaction of cations with the all-guanine quartets is the leading contribution for the stability of the four-stranded assemblies, while the mixed quartets are rather tolerated within the structure. The simulations predict that two cations are preferred to stabilize a four-layer quadruplex stem composed of two GCGC and two all-guanine quartets. The distribution of cations in the structure is influenced by the position of the GCGC quartets within the quadruplex, the presence and arrangement of thymidine loops connecting the guanine/cytosine stretches forming the stems, and the cation type present (Na(+) or K(+)). The simulations identify multiple nanosecond-scale stable arrangements of the thymidine loops present in the molecules investigated. In these thymidine loops, several structured pockets are identified capable of temporarily coordinating cations. However, no stable association of cations to a loop has been observed. The simulations reveal several paths through the thymidine loop regions that can be followed by the cations when exchanging between the central ion channel in the quadruplex stem and the surrounding solvent. We have carried out 20 independent simulations while the length of simulations reaches a total of 90 ns, rendering this study one of the most extensive MD investigations carried out on nucleic acids so far. The trajectories provide a largely converged characterization of the structural dynamics of these four-stranded G-DNA molecules.

Base Sequence↗
PubMed2001

Molecular dynamics of DNA quadruplex molecules containing inosine, 6-thioguanine and 6-thiopurine.

The ability of the four-stranded guanine (G)-DNA motif to incorporate nonstandard guanine analogue bases 6-oxopurine (inosine, I), 6-thioguanine (tG), and 6-thiopurine (tI) has been investigated using large-scale molecular dynamics simulations. The simulations suggest that a G-DNA stem can incorporate inosines without any marked effect on its structure and dynamics. The all-inosine quadruplex stem d(IIII)(4) shows identical dynamical properties as d(GGGG)(4) on the nanosecond time scale, with both molecular assemblies being stabilized by monovalent cations residing in the channel of the stem. However, simulations carried out in the absence of these cations show dramatic differences in the behavior of d(GGGG)(4) and d(IIII)(4). Whereas vacant d(GGGG)(4) shows large fluctuations but does not disintegrate, vacant d(IIII)(4) is completely disrupted within the first nanosecond. This is a consequence of the lack of the H-bonds involving the N2 amino group that is not present in inosine. This indicates that formation of the inosine quadruplex could involve entirely different intermediate structures than formation of the guanosine quadruplex, and early association of cations in this process appears to be inevitable. In the simulations, the incorporation of 6-thioguanine and 6-thiopurine sharply destabilizes four-stranded G-DNA structures, in close agreement with experimental data. The main reason is the size of the thiogroup leading to considerable steric conflicts and expelling the cations out of the channel of the quadruplex stem. The G-DNA stem can accommodate a single thioguanine base with minor perturbations. Incorporation of a thioguanine quartet layer is associated with a large destabilization of the G-DNA stem whereas the all-thioguanine quadruplex immediately collapses.

Biophysical Phenomena↗
PubMed2001

Solution structure of the dodecamer d-(CATGGGCC-CATG)2 is B-DNA. Experimental and molecular dynamics study.

The DNA duplex d-(CATGGGCCCATG)2 has been studied in solution by FTIR, NMR and CD. The experimental approaches have been complemented by series of large-scale unrestrained molecular dynamics simulation with explicit inclusion of solvent and counterions. Typical proton-proton distances extracted from the NMR spectra and the CD spectra are completely in agreement with slightly modified B-DNA. By molecular dynamics simulation, starting from A-type sugar pucker, a spontaneous repuckering to B-type sugar pucker was observed. Both experimental and theoretical approaches suggest for the dodecamer d-(CATGGGCCCATG)2 under solution conditions puckering of all 2'-deoxyribose residues in the south conformation (mostly C2'-endo) and can exclude significant population of sugars in the north conformation (C3'-endo). NMR, FTIR and CD data are in agreement with a B-form of the dodecamer in solution. Furthermore, the duplex shows a cooperative B-A transition in solution induced by addition of trifluorethanol. This contrasts a recently published crystal structure of the same oligonucleotide found as an intermediate between B- and A-DNA where 23 out of 24 sugar residues were reported to adopt the north (N-type) conformation (C3'-endo) like in A-DNA (Ng, H. L., Kopka, M. L. and Dickerson, R. E., Proc. Natl. Acad. Sci. U S A 97, 2035-2039 (2000)). The simulated structures resemble standard B-DNA. They nevertheless show a moderate shift towards A-type stacking similar to that seen in the crystal, despite the striking difference in sugar puckers between the MD and X-ray structures. This is in agreement with preceding MD reports noticing special stacking features of G-tracts exhibiting a tendency towards the A-type stacking supported by the CD spectra also reflecting the G-tract stacking. MD simulations reveal several noticeable local conformational variations, such as redistribution of helical twist and base pair roll between the central GpC steps and the adjacent G-tract segments, as well as a substantial helical twist variability in the CpA(TpG) steps combined with a large positive base pair roll. These local variations are rather different from those seen in the crystal.

Base Sequence↗
PubMed2000

Sequence-dependent elastic properties of DNA.

Harmonic elastic constants of 3-11 bp duplex DNA fragments were evaluated using four 5 ns unrestrained molecular dynamics simulation trajectories of 17 bp duplexes with explicit inclusion of solvent and counterions. All simulations were carried out with the Cornell et al. force-field and particle mesh Ewald method for long-range electrostatic interactions. The elastic constants including anisotropic bending and all coupling terms were derived by analyzing the correlations of fluctuations of structural properties along the trajectories. The following sequences have been considered: homopolymer d(ApA)(n) and d(GpG)(n), and alternating d(GPC)(n) and d(APT)(n). The calculated values of elastic constants are in very good overall agreement with experimental values for random sequences. The atomic-resolution molecular dynamics approach, however, reveals a pronounced sequence-dependence of the stretching and torsional rigidity of DNA, while sequence-dependence of the bending rigidity is smaller for the sequences considered. The earlier predicted twist-bend coupling emerged as the most important cross-term for fragments shorter than one helical turn. The calculated hydrodynamic relaxation times suggest that damping of bending motions may play a role in molecular dynamics simulations of long DNA fragments. A comparison of elasticity calculations using global and local helicoidal analyses is reported. The calculations reveal the importance of the fragment length definition. The present work shows that large-scale molecular dynamics simulations represent a unique source of data to study various aspects of DNA elasticity including its sequence-dependence.

Algorithms↗
PubMed2000

The interactions of square platinum(II) complexes with guanine and adenine: a quantum-chemical ab initio study of metalated tautomeric forms.

The influence of binding of square planar platinum complexes on tautomeric equilibria of the DNA bases guanine and adenine was investigated using the density functional B3LYP method. Neutral trans-dichloro(amine)-, +1 charged chloro(diamine)-, and +2 charged triamine-platinum(II) species were chosen for coordination to bases. Only the N7 interaction site of the bases was considered. The calculations demonstrate that the neutral platinum adduct does not change the tautomeric equilibria of the bases. Furthermore, N7 binding of the neutral Pt adduct moderately reduces the probability of protonation of the N1 position of adenine. Larger effects can be observed for +1 and mainly +2 adducts, but these can be rationalized by electrostatic effects. Since the electrostatic effects are expected to be efficiently compensated for by a charged backbone of DNA and counterions in a polar solvent, no dramatic increase in mispair formation is predicted for Pt(II) adducts, which is in agreement with experiment. The interaction energies between Pt adducts and the nucleobases were also evaluated. These interaction energies range from ca. 210 kJ/mol for neutral adducts, interacting with both bases and their tautomers, up to 500 kJ/mol for the +2 charged adducts, interacting with the keto-guanine tautomer and the anti-imino-adenine tautomer. The surprisingly large interaction energy for the latter structure is due to the strong H-bond between the NH3 ligand group of the metal adduct and the N6 nitrogen atom of the base.

Adenine↗
PubMed2000

Structure, energetics, vibrational frequencies and charge transfer of base pairs, nucleoside pairs, nucleotide pairs and B-DNA pairs of trinucleotides: ab initio HF/MINI-1 and empirical force field study.

Geometries, interaction energies and vibrational frequencies of base pairs, nucleoside pairs and nucleotide pairs were studied by ab initio Hartree-Fock (HF) method using MINI-1 basis set and empirical Cornell et al. force field (AMBER 4.1). A good agreement was found between HF/MINI-1 and AMBER results. In addition, both methods provide reasonable agreement with available high-level ab initio data. Finally, AMBER potential was used to determine the structure, energetics and vibrational frequencies of B-DNA pairs of trinucleotides. Stabilization energies of clusters are lowered when passing from base pairs to nucleoside pairs, nucleotide pairs and to pairs of trinucleotides. The lowest vibrations of base pairs and nucleoside pairs correspond to intermolecular motions of bases, specifically to buckle and propeller motions. In the case of pairs of larger subunits the lowest vibrations are of intramolecular nature (rotation around glycosidic bond, sugar and phosphate vibration). The spectra of these clusters became more complicated and quasi-degenerate. Intermolecular charge transfer between bases in H-bonded and stacked pairs is negligible, while a significant intramolecular charge transfer was observed.

Base Pairing↗
PubMed2000

Cation-pi and amino-acceptor interactions between hydrated metal cations and DNA bases. A quantum-chemical view.

Cation-pi interactions between cytosine and hexahydrated cations have been characterized using ab initio method with inclusion of electron correlation effects, assuming idealized and crystal geometries of the interacting species. Hydrated metal cations can interact with nucleobases in a cation-pi manner. The stabilization energy of such complexes would be large and comparable to the one for cation-pi complex with benzene. Further, polarized water molecules belonging to the hydration shell of the cation are capable to form a strong hydrogen bond interaction with the nitrogen lone electron pair of the amino groups of bases and enforce a pronounced sp3 pyramidalization of the nucleobase amino groups. However, in contrast to the benzene-cation complexes, the cation-pi configurations are highly unstable for a nucleobase since the conventional in plane binding of hydrated cations to the acceptor sites on the nucleobase is strongly preferred. Thus, a cation-pi interaction with a nucleobase can occur only if the position of the cation is locked above the nucleobase plane by another strong interaction. This indeed can occur in biopolymers and may have an effect on the local DNA architecture. Nevertheless, nucleobases have no intrinsic propensity to form cation-pi interactions.

Cations↗
PubMed1999

Crystal structure of d(GGCCAATTGG) complexed with DAPI reveals novel binding mode.

The single-crystal X-ray structure of the complex between the minor groove binder 4',6-diamidino-2-phenylindole (DAPI) and d(GGCCAATTGG) reveals a novel way of off-centered binding, with an unique hydrogen bond between the minor groove binder and a CG base pair. Application of crystal engineering and cryocooling techniques helped to extend the resolution to 1.9 A, resulting in an unambiguous determination of drug conformation and orientation. The structure was refined to completion using SHELXL-93, resulting in a residual factor R of 18. 0% for 3562 reflections with F(o) > 4sigma(F(o)) including 81 water molecules. As the bulky NH(2)-group on guanine is believed to prevent drug binding in the minor groove, the nature and stability of the CG-DAPI contact was further addressed in full detail using ab initio quantum chemical methods. The amino groups involved in the guanine-drug interaction are substantially nonplanar, resulting in an energy gain of about 5 kcal/mol. The combined structural and theoretical data suggest that the guanine NH(2)-group does not destabilize the drug binding to an extent that it prevents complexation.

Base Pairing↗
PubMed1999

Metal ions in non-complementary DNA base pairs: an ab initio study of Cu(I), Ag(I), and Au(I) complexes with the cytosine-adenine base pair.

Ab initio calculations have been carried out to characterize the structure and energetics of a silver(I) complex with the cytosine-adenine DNA base pair and an aqua ligand in the coordination sphere of Ag. In addition, we have also studied analogous complexes with Cu(I) and Au(I), and structures in which adenine has been replaced by purine in order to investigate the structural role of the adenine amino group. The calculations revealed that all metal-modified structures are dominated by the metal-base interactions, while the water-metal ion interaction and many-body interligand repulsion are less important contributions. Nevertheless, the structural role of the water molecule in the complex is quite apparent and in agreement with an earlier crystallographic study. The metal-modified base pairs exhibit large conformational flexibility toward out-of-plane motions (propeller twist and buckle), comparable or, in some cases, even larger than that observed in the base pairs without metal ions. All structures have been optimized within the Hartree-Fock approximation, while interaction energies were evaluated with the inclusion of electron correlation.

Adenine↗
PubMed1999

Interactions of hydrated IIa and IIb group metal cations with thioguanine-cytosine DNA base pair: Ab initio and density functional theory investigation of polarization effects, differences among cations, and flexibility of the cation hydration shell.

The structures and energies of the thioguanine-cytosine Watson-Crick (thioGC WC) base pair interacting with hydrated IIa (Mg2+, Ca2+, Ba2+) and IIb group (Zn2+, Cd2+, Hg2+) cations have been studied using ab initio techniques. Furthermore, complexes between guanine and thioguanine with hydrated cations have been characterized assuming various structures of the hydration shells. The complexes of the thioGC WC base pair with hydrated cations have similar properties as the previously studied GC WC base pair. There is substantial polarization stabilization of the base pairing due to cation binding which amounts to 7 - 11 kcal/mol. Soft Cd2+ and Hg2+ cations have a uniquely strong interaction with the thiogroup and induce substantial nonplanarity of the pairing. The thiogroup tends to reduce the number of water molecules in the first hydration shell of the cation. All complexes were optimized within the Hartree-Fock (HF) approximation while their energetics has been evaluated using the second-order Moller-Plesset perturbational method (MP2). All interaction energy evaluations and a substantial portion of the optimizations of the hydrated cation-(thio)guanine complexes have been repeated using Becke-3LYP Density Functional Theory method. All three approximations used (HF, Becke-3LYP, and MP2) give qualitatively the same results for the present cationic complexes. The results demonstrate specific differences among the cations and provide a set of reference structures and energies for verification and/or parametrization of empirical potentials and other theoretical methods.

Barium↗
PubMed1998

On the potential role of the amino nitrogen atom as a hydrogen bond acceptor in macromolecules.

Crystallographic studies of duplex DNA have indicated that opposing exocyclic amino groups may form close NH⋯:N contacts. To study the nature of such interactions, we have examined the database of small molecule, high-resolution crystal structures for more accurate examples of this type of unconventional interaction. We have found cases where the amino groups in guanine and adenine bases accept hydrogen bonds from conventional donors, such as amino or hydroxyl groups. More frequently, the purine amino group was found to contact closely electropositive C-H groups. Searches of the nucleic acid structural databases also yielded several examples where the purine amino group is contacted by hydrogen bond donors in macromolecules. Ab initio calculations indicate that the hydrogen-amino contact is improved energetically when the amino group moves from the conventional geometry, where all atoms are co-planar with the base, to one in which the hydrogen atoms lie out of the plane and the nitrogen is at the apex of a pyramid, resulting in polarization of the amino group. The combined structural and theoretical data suggest that the amino group is flexible, and can accommodate close contacts, because the resulting polarization permits electropositive atoms to approach the amino group nitrogen more closely than expected for their conventional van der Waals radii. The flexibility of the amino group may permit particular DNA conformations that enforce hydrogen-amino contacts to optimize favorable stacking interactions, and it may play a role in the recognition of nucleosides. We speculate that the amino group can accept hydrogen bonds under special circumstances in macromolecules, and that this ability might play a mechanistic role in catalytic processes such as deamination or amino transfer.

Crystallography↗
PubMed1997

Base-base and deoxyribose-base stacking interactions in B-DNA and Z-DNA: a quantum-chemical study.

Base-stacking interactions in canonical and crystal B-DNA and in Z-DNA steps are studied using the ab initio quantum-chemical method with inclusion of electron correlation. The stacking energies in canonical B-DNA base-pair steps vary from -9.5 kcal/mol (GG) to -13.2 kcal/mol (GC). The many-body nonadditivity term, although rather small in absolute value, influences the sequence dependence of stacking energy. The base-stacking energies calculated for CGC and a hypothetical TAT sequence in Z-configuration are similar to those in B-DNA. Comparison with older quantum-chemical studies shows that they do not provide even a qualitatively correct description of base stacking. We also evaluate the base-(deoxy)ribose stacking geometry that occurs in Z-DNA and in nucleotides linked by 2',5'-phosphodiester bonds. Although the molecular orbital analysis does not rule out the charge-transfer n-pi* interaction of the sugar 04' with the aromatic base, the base-sugar contact is stabilized by dispersion energy similar to that of stacked bases. The stabilization amounts to almost 4 kcal/mol and is thus comparable to that afforded by normal base-base stacking. This enhancement of the total stacking interaction could contribute to the propensity of short d(CG)n sequences to adopt the Z-conformation.

Base Composition↗
PubMed1997

Hydrogen-bonded trimers of DNA bases and their interaction with metal cations: ab initio quantum-chemical and empirical potential study.

Neutral (G.GC, A.AT, G.AT, T.AT, and C(imino).GC) and protonated (CH+.GC and AH+.GC) hydrogen-bonded trimers of nucleic acid bases were characterized by ab initio methods with the inclusion of electron correlation. In addition, the influence of metal cations on the third-strand binding in Purine-Purine-Pyrimidine (Pu.PuPy) reverse-Hoogsteen triplets has been studied. The ab initio calculations were compared with those from recently introduced force fields (AMBER4.1, CHARMM23, and CFF95). The three-body term in neutral trimers is mostly negligible, and the use of empirical potentials is justified. The only exception is the neutral G.GC Hoogsteen trimer with a three-body term of -4 kcal/mol. Protonated trimers are stabilized by molecular ion-molecular dipole attraction and the interaction within the complex is nonadditive, with the three-body term on the order of -3 kcal/mol. There is a significant induction interaction between the third-strand protonated base and guanine. The calculations indicate an enhancement of the third-strand binding in the G.GC reverse-Hoogsteen trimer due to-metal cation coordination to the N7/O6 position of the third-strand guanine. Interactions between metal cations and complexes of DNA bases are in general highly non-additive; the three-body term is above-10 kcal/mol in a complex of a divalent cation (Ca2+) with the GG reverse-Hoogsteen pair. The pairwise additive empirical potentials qualitatively underestimate the binding energy between cation and base.

Adenine↗
PubMed1996

Base stacking and hydrogen bonding in protonated cytosine dimer: the role of molecular ion-dipole and induction interactions.

An ab initio quantum-chemical study of stacked and hydrogen-bonded protonated cytosine dimer has been carried out. The calculations were made using the second-order Moller-Plesset perturbational method (MP2) with a medium-sized polarized set of atomic orbitals. H-bonded as well as stacked protonated base pairs are more stable than the neutral base pairs. Two energy contributions not present in the neutral base pairs stabilize the protonated base pairs: the molecular ion - dipole interaction, and the induction interaction. The molecular ion - dipole stabilization dominates in base pairs with highly polar neutral monomers, such as the C...CH+ base pair. The induction interaction is not included in the commonly used empirical potentials, which do not reproduce the changes in intermolecular stabilization due to protonation. We demonstrate that the base stacking of several consecutive C...CH + pairs, as proposed for polycytidylic acid and i-DNA, is strongly repulsive. We also show that the intermolecular interactions strongly prefer protonation of adenine in protonated adenine-cytosine pairs.

Base Composition↗