Molecular dynamics simulation of the A-DNA to B-DNA transition in aqueous RbCl solution

Unrestrained molecular dynamics (MD) simulations have been carried out to characterize the stability of DNA conformations and the dynamics of A-DNA→B-DNA conformational transitions in aqueous RbCl solutions. The PARM99 force field in the AMBER8 package was used to investigate the effect of RbCl concentration on the dynamics of the A→B conformational transition in the DNA duplex d(CGCGAATTCGCG)2 . Canonical Aand B-form DNA were assumed for the initial conformation and the final conformation had a length per complete turn that matched the canonical B-DNA. The DNA structure was monitored for 3.0 ns and the distances between the C5′ atoms were obtained from the simulations. It was found that all of the double stranded DNA strands of A-DNA converged to the structure of B-form DNA within 1.0 ns during the unrestrained MD simulations. In addition, increasing the RbCl concentration in aqueous solution hindered the A→B conformational transition and the transition in aqueous RbCl solution was faster than that in aqueous NaCl solution for the same electrolyte strength. The effects of the types and concentrations of counterions on the dynamics of the A→B conformational transition can be understood in terms of the variation in water activity and the number of accumulated counterions in the major grooves of A-DNA. The rubidium ion distributions around both fixed A-DNA and B-DNA were obtained using the restrained MD simulations to help explain the effect of RbCl concentration on the dynamics of the A→B conformational transition.     

Conformationally Restricted Isoindoline‐Derived Spin Labels in Duplex DNA: Distances and Rotational Flexibility by Pulsed Electron–Electron Double Resonance Spectroscopy

Three structurally related isoindoline‐derived spin labels that have different mobilities were incorporated into duplex DNA to systematically study the effect of motion on orientation‐dependent pulsed electron–electron double resonance (PELDOR) measurements. To that end, a new nitroxide spin label, ^ExIm U , was synthesized and incorporated into DNA oligonucleotides. ^ExIm U is the first example of a conformationally unambiguous spin label for nucleic acids, in which the nitroxide N?O bond lies on the same axis as the three single bonds used to attach the otherwise rigid isoindoline‐based spin label to a uridine base. Continuous‐wave (CW) EPR measurements of ^ExIm U confirm a very high rotational mobility of the spin label in duplex DNA relative to the structurally related spin label ^Im U , which has restricted mobility due to an intramolecular hydrogen bond. The X‐band CW‐EPR spectra of ^ExIm U can be used to identify mismatches in duplex DNA. PELDOR distance measurements between pairs of the spin labels ^Im U , ^Ox U , and ^ExIm U in duplex DNA showed a strong angular dependence for ^Im U , a medium dependence for ^Ox U , and no orientation effect for ^ExIm U . Thus, precise distances can be extracted from ^ExIm U without having to take orientational effects into account.     

PELDOR on an exchange coupled nitroxide copper(II) spin pair

Transition metal ions play an important role in the design of macromolecular architectures as well as for the structure and function of proteins and oligonucleotides, which makes them interesting targets for spectroscopic investigations. In combination with site directed spin labelling, pulsed electron–electron double resonance (PELDOR or DEER) could be a well-suited method for their characterization and localization. Here, we report on the synthesis and full characterization of a copper(II) porphyrin/nitroxide model system bearing an extended π-conjugation between the spin centres and demonstrate the possibility to disentangle the dipolar through space interaction from the through bond exchange coupling contribution even in the presence of orientational selectivity and conformational flexibility. The simulations used are based on the known experimental and spin Hamiltonian parameters and on a structural model as previously employed for similar systems. The mean exchange coupling of +4(1) MHz (antiferromagnetic) is in agreement with the value of |J| = 3(1) MHz determined from room temperature continuous wave electron paramagnetic resonance (EPR). Thus, as long as the pulse excitation bandwidths are large versus the spin–spin coupling, X-band PELDOR measurements in combination with explicit time trace simulations allow for disentangling the sign and magnitude of through bond electron–electron exchange from the through space dipolar interaction D.     

Post‐synthetic Spin‐Labeling of RNA through Click Chemistry for PELDOR Measurements

Site‐directed spin labeling of RNA based on click chemistry is used in combination with pulsed electron‐electron double resonance (PELDOR) to benchmark a nitroxide spin label, called here d? . We compare this approach with another established method that employs the rigid spin label Çm for RNA labeling. By using CD spectroscopy, thermal denaturation measurements, CW‐EPR as well as PELDOR we analyzed and compared the influence of d? and Çm on a self‐complementary RNA duplex. Our results demonstrate that the conformational diversity of d? is significantly reduced near the freezing temperature of a phosphate buffer, resulting in strongly orientation‐selective PELDOR time traces of the d? ‐labeled RNA duplex.     

Dynamics of Nucleic Acids at Room Temperature Revealed by Pulsed EPR Spectroscopy

The investigation of the structure and conformational dynamics of biomolecules under physiological conditions is challenging for structural biology. Although pulsed electron paramagnetic resonance (like PELDOR) techniques provide long‐range distance and orientation information with high accuracy, such studies are usually performed at cryogenic temperatures. At room temperature (RT) PELDOR studies are seemingly impossible due to short electronic relaxation times and loss of dipolar interactions through rotational averaging. We incorporated the rigid nitroxide spin label Ç into a DNA duplex and immobilized the sample on a solid support to overcome this limitation. This enabled orientation‐selective PELDOR measurements at RT. A comparison with data recorded at 50 K revealed averaging of internal dynamics, which occur on the ns time range at RT. Thus, our approach adds a new method to study structural and dynamical processes at physiological temperature in the <10 μs time range with atomistic resolution.     

Synthesis of μ2-Oxo-Bridged Iron(III) Tetraphenylporphyrin–Spacer–Nitroxide Dimers and their Structural and Dynamics Characterization by using EPR and MD Simulations

Iron(III) porphyrins have the propensity to form μ₂-oxo-dimers, the structures of which resemble two wheels on an axle. Whereas their crystal structure is known, their solution structure and internal dynamics is not. In the present work, the structure and dynamics of such dimers were studied by means of electron paramagnetic resonance (EPR) spectroscopy and quantum chemistry based molecular dynamics (MD) simulations by using the semiempirical tight-binding method (GFN-xTB). To enable EPR investigation of the dimers, a nitroxide was attached to each of the tetraphenylporphyrin cores through a linear and a bent linker. The inter-nitroxide distance distributions within the dimers were determined by continuous-wave (cw)-EPR and pulsed electron–electron double resonance (PELDOR or DEER) experiments and, with the help of MD, interpreted in terms of the rotation of the porphyrin planes with respect to each other around the Fe–O–Fe axis. It was found that such rotation is restricted to the four registers defined by the phenyl substituents. Within the registers, the rotation angle swings between 30° and 60° in the proximal and between 125° and 145° in the distal register. With EPR, all four angles were found to be equally populated, whereas the 30° and 145° angles are strongly favored to the expense of the 60° and 125° angles in the MD simulation. In either case, the internal dynamics of these dimers thus resemble the motion of a step motor.     

Counting the monomers in nanometer-sized oligomers by pulsed electron-electron double resonance

In a lot of cases active biomolecules are complexes of higher order, thus methods capable of counting the number of building blocks and elucidating their geometric arrangement are needed. Therefore, we experimentally validate here spin-counting via 4-pulse electron-electron double resonance (PELDOR) on well-defined test samples. Two biradicals, a symmetric and an asymmetric triradical, and a tetraradical were synthesized in a convergent reaction scheme via palladium-catalyzed cross-coupling reactions. PELDOR was then used to obtain geometric information and the number of spin centers per molecule in a single experiment. The measurement yielded the expected distances (2.2-3.8 nm) and showed that different spin-spin distances in one molecule can be resolved even if the difference amounts to only 5 A. The number of spins n has been determined to be 2.1 in both biradicals, to 3.1 and 3.0 in the symmetric and asymmetric triradicals, respectively, and to 3.9 in the tetraradical. The overall error of PELDOR spin-counting was found to be 5% for up to four spins. Thus, this method is a valuable tool to determine the number of constituting spin-bearing monomers in biologically relevant homo- and heterooligomers and how their oligomerization state and geometric arrangement changes during function.     

Light‐Induced Porphyrin‐Based Spectroscopic Ruler for Nanometer Distance Measurements

We present a novel pulsed electron paramagnetic resonance (EPR) spectroscopic ruler to test the performance of a recently developed spin‐labeling method based on the photoexcited triplet state (S=1). Four‐pulse electron double resonance (PELDOR) experiments are carried out on a series of helical peptides, labeled at the N‐terminal end with the porphyrin moiety, which can be excited to the triplet state, and with the nitroxide at various sequence positions, spanning distances in the range 1.8–8 nm. The PELDOR traces provide accurate distance measurements for all the ruler series, showing deep envelope modulations at frequencies varying in a progressive way according to the increasing distance between the spin labels. The upper limit is evaluated and found to be around 8 nm. The PELDOR‐derived distances are in excellent agreement with theoretical predictions. We demonstrate that high sensitivity is acquired using the triplet state as a spin label by comparison with Cu(II)–porphyrin analogues. The new labeling approach has a high potential for measuring nanometer distances in more complex biological systems due to the properties of the porphyrin triplet state.     

Accurate Extraction of Nanometer Distances in Multimers by Pulse EPR

Pulse electron paramagnetic resonance (EPR) is gaining increasing importance in structural biology. The PELDOR (pulsed electron–electron double resonance) method allows extracting distance information on the nanometer scale. Here, we demonstrate the efficient extraction of distances from multimeric systems such as membrane‐embedded ion channels where data analysis is commonly hindered by multi‐spin effects.     

Theoretical study of a new DNA structure: the antiparallel Hoogsteen duplex

The structure of a new form of duplex DNA, the antiparallel Hoogsteen duplex, is studied in polyd(AT) sequences by means of state-of-the-art molecular dynamics simulations in aqueous solution. The structure, which was found to be stable in all of the simulations, has many similarities with the standard Watson-Crick duplex in terms of general structure, flexibility, and molecular recognition patterns. Accurate MM-PB/SA (and MM-GB/SA) analysis shows that the new structure has an effective energy similar to that of the B-type duplex, while it is slightly disfavored by intramolecular entropic considerations. Overall, MD simulations strongly suggest that the antiparallel Hoogsteen duplex is an accessible structure for a polyd(AT) sequence, which might compete under proper experimental conditions with normal B-DNA. MD simulations also suggest that chimeras containing Watson-Crick duplex and Hoogsteen antiparallel helices might coexist in a common structure, but with the differential characteristics of both type of structures preserved.     

A Bis‐Manganese(II)–DOTA Complex for Pulsed Dipolar Spectroscopy

High‐spin gadolinium(III) and manganese(II) complexes have emerged as alternatives to standard nitroxide radical spin labels for measuring nanometric distances by using pulsed electron–electron double resonance (PELDOR or DEER) at high fields/frequencies. For certain complexes, particularly those with relatively small zero‐field splitting (ZFS) and short distances between the two metal centers, the pseudosecular term of the dipolar coupling Hamiltonian is non‐negligible. However, in general, the contribution from this term during conventional data analysis is masked by the flexibility of the molecule of interest and/or the long tethers connecting them to the spin labels. The efficient synthesis of a model system consisting of two [Mn(dota)]^2? (MnDOTA; DOTA^4?=1,4,7,10‐tetraazacyclododecane‐1,4,7,10‐tetraacetate) directly connected to the ends of a central rodlike oligo(phenylene–ethynylene) (OPE) spacer is reported. The rigidity of the OPE is confirmed by Q‐band PELDOR measurements on a bis‐nitroxide analogue. The Mn^II?Mn^II distance distribution profile determined by W‐band PELDOR is in reasonable agreement with one simulated by using a simple rotamer analysis. The small degree of flexibility arising from the linking MnDOTA arm appears to outweigh the contribution from the pseudosecular term at this interspin distance. This study illustrates the potential of MnDOTA‐based spin labels for measuring fairly short nanometer distances, and also presents an interesting candidate for in‐depth studies of pulsed dipolar spectroscopy methods on Mn^II?Mn^II systems.     

PELDOR measurements on a nitroxide-labeled Cu(II) porphyrin: orientation selection, spin-density distribution, and conformational flexibility

Metal ions are functionally or structurally important centers in metalloproteins or RNAs, which makes them interesting targets for spectroscopic investigations. In combination with site-directed spin labeling, pulsed electron-electron double resonance (PELDOR or DEER) could be a well-suited method to characterize and localize them. Here, we report on the synthesis, full characterization, and PELDOR study of a copper(II) porphyrin/nitroxide model system. The X-band PELDOR time traces contain besides the distance information a convolution of orientational selectivity, conformational flexibility, exchange coupling, and spin density distribution, which can be deconvoluted by experiments with different frequency offsets and simulations. The simulations are based on the known experimental and spin Hamiltonian parameters and make use of a geometric model as employed for structurally similar bis-nitroxides and spin density parameters as obtained from density functional theory calculations. It is found that orientation selection with respect to dipolar angles is only weakly resolvable at X-band frequencies due to the large nitrogen hyperfine coupling of the copper porphyrin. On the other hand, the PELDOR time traces reveal a much faster oscillation damping than observed for structurally similar bis-nitroxides, which is mainly assigned to a small distribution in exchange couplings J. Taking the effects of orientation selectivity, distribution in J, and spin density distribution into account leads finally to a narrow distance distribution caused solely by the flexibility of the structure, which is in agreement with distributions from known bis-nitroxides of similar structure. Thus, X-band PELDOR measurements at different frequency offsets in combination with explicit time trace simulations allow for distinguishing between structural models and quantitative interpretation of copper-nitroxide PELDOR data gives access to localization of copper(II) ions.     

Free energy landscape of A-DNA to B-DNA conversion in aqueous solution

The interconversion between the well-characterized A- and B-forms of DNA is a structural transition for which the intermediate states and the free energy difference between the two endpoints are not known precisely. In the present study, the difference between the Root Mean Square Distance (RMSD) from canonical A-form and B-form DNA is used as an order parameter to characterize this free energy difference using umbrella sampling molecular dynamics (MD) simulations with explicit solvent. The constraint imposed along this order parameter allows relatively unrestricted evolution of the intermediate structures away from both canonical A- and B-forms. The free energy difference between the A- and B-forms for the hexamer DNA sequence CTCGAG in aqueous solution is conservatively estimated to be at least 2.8 kcal/mol. A continuum of intermediate structures with no well-defined local minima links the two forms. The absence of any major barriers in the free energy surface is consistent with spontaneous conversion of the A-form DNA to B-form DNA in unconstrained simulations. The extensive sampling in the MD simulations (>0.1 mus) also allowed quantitative energetic characterization of local backbone conformational variables such as sugar pseudorotation angles and BI/BII state equilibria and their dependence on base identity. The absolute minimum in the calculated free energy profile corresponds closely to the crystal structure of the hexamer sequence, indicating that the present method has the potential to identify the most stable state for an arbitrary DNA sequence in water.     

Nature of minor-groove binders-DNA complexes in the gas phase

The structure of noncovalent complexes of DNA duplex with minor groove binders (mG-binders) has been analyzed by state of the art molecular dynamics (MD) simulations. More than 3.3 micros of MD trajectories (including 4 x 0.5 micros trajectories) were collected for the Dickerson's dodecamer bound to DAPI, Hoechst 33258, and Netropsin. Comparison of these trajectories with control simulations in water allowed us to determine that the extreme dehydration and partial neutralization occurring during electrospray experiments does not produce the disruption of the DNA:mG-binder complexes or the dissociation of the two strands of the duplex. Irrespective of the drug and the simulation conditions the mG-binders remains bound to the DNA near the preferential binding position in aqueous conditions. Large distortions appear in the two DNA strands, which maintain however a memory of the original DNA duplex structure in water, and a general helical-like conformation.     

Selective High‐Resolution Detection of Membrane Protein–Ligand Interaction in Native Membranes Using Trityl–Nitroxide PELDOR

The orchestrated interaction of transmembrane proteins with other molecules mediates several crucial biological processes. Detergent solubilization may significantly alter or even abolish such hetero‐oligomeric interactions, which makes observing them at high resolution in their native environment technically challenging. Dipolar electron paramagnetic resonance (EPR) techniques such as pulsed electro–electron double resonance (PELDOR) can provide very precise distances within biomolecules. To concurrently determine the inter‐subunit interaction and the intra‐subunit conformational changes in hetero‐oligomeric complexes, a combination of different spin labels is required. Orthogonal spin labeling using a triarylmethyl (TAM) label in combination with a nitroxide label is used to detect protein–ligand interactions in native lipid bilayers. This approach provides a higher sensitivity and total selectivity and will greatly facilitate the investigation of multimeric transmembrane complexes employing different spin labels in the native lipid environment.     

Unusual long-range spin-spin coupling in fluorinated polyenes: a mechanistic analysis

Nuclear magnetic resonance (NMR) is a prospective means to realize quantum computers. The performance of a NMR quantum computer depends sensitively on the properties of the NMR-active molecule used, where one requirement is a large indirect spin-spin coupling over large distances. F-F spin-spin coupling constants (SSCCs) for fluorinated polyenes F-(CH==CH)n-F (n=1cdots, three dots, centered5) are >9 Hz across distances of more than 10 A. Analysis of the F,F spin-spin coupling mechanism with our recently developed decomposition of J into Orbital Contributions with the help of Orbital Currents and Partial Spin Polarization (J-OCOC-PSP=J-OC-PSP) method reveals that coupling is dominated by the spin-dipole (SD) term due to an interplay between the pi lone-pair orbitals at the F atoms and the pi(C2n) electron system. From our investigations we conclude that SD-dominated SSCCs should occur commonly in molecules with a contiguous pi-electron system between the two coupling nuclei and that a large SD coupling generally is the most prospective way to provide large long-range spin-spin coupling. Our results give guidelines for the design of suitable active molecules for NMR quantum computers.     

Theoretical studies of d(A:T)-based parallel-stranded DNA duplexes.

Poly d(A:T) parallel-stranded DNA duplexes based on the Hoogsteen and reverse Watson-Crick hydrogen bond pairing are studied by means of extensive molecular dynamics (MD) simulations and molecular mechanics coupled to Poisson-Boltzmann (MM-PB/SA) calculations. The structural, flexibility, and reactivity characteristics of Hoogsteen and reverse Watson-Crick parallel duplexes are described from the analysis of the trajectories. Theoretical calculations show that the two parallel duplexes are less stable than the antiparallel Watson-Crick duplex. The difference in stability between antiparallel and parallel duplexes increases steadily as the length of the duplex increases. The reverse Watson-Crick arrangement is slightly more stable than the Hoogsteen duplex, the difference being also increased linearly with the length of the duplex. A subtle balance of intramolecular and solvation terms is responsible for the preference of a given helical structure.     

NMR solution structure of an oligonucleotide hairpin with a 2'F-ANA/RNA stem: implications for RNase H specificity toward DNA/RNA hybrid duplexes.

The first structure of a 2'-deoxy-2'-fluoro-D-arabinose nucleic acid (2'F-ANA)/RNA duplex is presented. We report the structural characterization by NMR spectroscopy of a small hybrid hairpin, r(GGAC)d(TTCG)2'F-a(GTCC), containing a 2'F-ANA/RNA stem and a four-residue DNA loop. Complete (1)H, (13)C, (19)F, and (31)P resonance assignments, scalar coupling constants, and NOE constraints were obtained from homonuclear and heteronuclear 2D spectra. In the chimeric duplex, the RNA strand adopts a classic A-form structure having C3' endo sugar puckers. The 2'F-ANA strand is neither A-form nor B-form and contains O4' endo sugar puckers. This contrasts strongly with the dynamic sugar conformations previously observed in the DNA strands of DNA/RNA hybrid duplexes. Structural parameters for the duplex, such as minor groove width, x-displacement, and inclination, were intermediate between those of A-form and B-form duplexes and similar to those of DNA/RNA duplexes. These results rationalize the enhanced stability of 2'F-ANA/RNA duplexes and their ability to elicit RNase H activity. The results are relevant for the design of new antisense drugs based on sugar-modified nucleic acids.