Conformations produced by interactions of the side chains in poly(silylenemethylene) with the repeating sequence [Si(CH3)R?CH2]x [R = ?O(CH2)3OC6H4C6H5]

A coarse-grained model has been developed for asymmetrically substituted poly(silylenemethylene)s in which the side chain is a flexible spacer terminated by a biphenyl unit. Each monomer unit is represented by four coarse-grained beads that interact via a Lennard–Jones potential and are subject to the first- and second-order interactions deduced from the atomistically detailed model. Metropolis Monte Carlo simulations were performed for isolated syndiotactic, isotactic, and atactic chains. Snapshots from the equilibrated coarse-grained chain on the discrete space of a high coordination lattice were reverse-mapped to atomistically detailed structures in continuous space. At 373 K, the chains were disordered independent of the stereochemical composition. The occupancy of bond pairs depended on the stereochemical composition, with the trans-gauche (tg) sequence being favored by the isotactic chain. When the simulation was performed with the backbone constrained to specific periodic structures, the g helix was the lowest energy structure for either the atactic or isotactic chains. For the syndiotactic chain, the g and gt helices were favored. The appearance of the g helix as the favored periodic structure of the isolated chain was consistent with the chain conformation reported previously for the smectic phase of this polymer in the bulk state. The g helix was disrupted when the backbone was allowed to access nonhelical conformations, even though these conformations may have been slightly higher in energy. © 2005 Wiley Periodicals, Inc. J Polym Sci Part B: Polym Phys 43: 886–896, 2005     

Statistics of time-limited ensembles of bent DNA conformations

The paper considers statistical properties of ensembles of chain conformations obtained by short-time Brownian dynamics (BD) of a coarse-grained DNA model in order to find out if the conditions necessary for accurate evaluation of the polymer elasticity are attainable in atom-level molecular dynamics (MD) simulations. To measure the bending persistence length (PL) with a 10% error using data accumulated in a single trajectory of a double helix of 15 base pairs, dynamics should be continued for a few microseconds. However, these estimates should be scaled down by about 2 orders of magnitude because the bending dynamics of short double helices in MD features much smaller relaxation times. As a result, good qualitative agreement with the worm-like chain (WLC) theory is reached in MD after tens of nanoseconds. The presently accessible durations of MD trajectories provide reasonably accurate evaluation of DNA elasticity and allow modeling of its mesoscopic properties. The surprisingly fast bending dynamics of short double helices in MD suggests that the microscopic mechanisms of DNA flexibility differ from a simple harmonic model.     

Duplex structure of a minimal nucleic acid

The crystal structure of an 8-mer (S)-GNA duplex is presented. As a tool for phasing, the anomalous diffraction of two copper(II) ions within two artificial metallo-base pairs was employed. The duplex structure confirms a canonical Watson-Crick base pairing scheme of GNA with antiparallel strands. The duplex secondary structure is distinct from canonical A- and B-form nucleic acids and can be described as a right-handed helical ribbon wrapped around the helix axis, resulting in a large hollow core. Most intriguingly, neighboring base pairs slide strongly against each other, resulting in extensive interstrand base-base hydrophobic interactions along with unusual hydrophobic intrastrand interactions of nucleobases with their backbone. These results reveal how a minimal nucleic acid backbone can support highly stable Watson-Crick-like duplex formation.     

MD simulations of homomorphous PNA, DNA, and RNA single strands: characterization and comparison of conformations and dynamics

MD simulations of homomorphous single-stranded PNA, DNA, and RNA with the same base sequence have been performed in aqueous solvent. For each strand two separate simulations were performed starting from a (i) helical conformation and (ii) random coiled state. Comparisons of the simulations with the single-stranded helices (case i) show that the differences in the covalent nature of the backbones cause significant differences in the structural and dynamical properties of the strands. It is found that the PNA strand maintains its nice base-stacked initial helical structure throughout the 1.5-ns MD simulation at 300 K, while DNA/RNA show relatively larger fluctuations in the structures with a few local unstacking events during -ns MD simulation each. It seems that the weak physical coupling between the bases and the backbone in PNA causes a loss of correlation between the dynamics of the bases and the backbone compared to the DNA/RNA and helps maintain the base-stacked helical conformation. The global flexibility of a single-stranded PNA helix was also found to be lowest, while RNA appears to be the most flexible single-stranded helix. The sugar pucker of several nucleotides in single-stranded DNA and RNA was found to adopt both C2'-endo and C3'-endo conformations for significant times. This effect is more pronounced for single strands in completely coiled states. The simulations with single-stranded coils as the initial structure also indicate that a PNA can adopt a more compact globular structure, while DNA/RNA of the same size adopts a more extended coil structure. This allows even a short PNA in the coiled state to form a significantly stable nonsequentially base-stacked globular structure in solution. Due to the hydrophobic nature of the PNA backbone, it interacts with surrounding water rather weakly compared to DNA/RNA.     

A Janus-Wedge DNA triplex with A-W1-T and G-W2-C base triplets

A new type of double-stranded DNA targeting format by formation of a Janus-Wedge (J-W) triple helix is described. The "wedge" residue W1 is used for A-T and T-A base pairs while W2 is used for G-C and C-G base pairs. Both wedge residues are attached to a PNA backbone that is designed to insert the probe strand into double-stranded DNA and base pair with both Watson-Crick faces. To study the stability of such an assembly, we have examined the formation of the J-W triplex with various sequences.     

Thermodynamics of twisted DNA with solvent interaction

The imaginary time path integral formalism is applied to a nonlinear Hamiltonian for a short fragment of heterogeneous DNA with a stabilizing solvent interaction term. Torsional effects are modeled by a twist angle between neighboring base pairs stacked along the molecule backbone. The base pair displacements are described by an ensemble of temperature dependent paths thus incorporating those fluctuational effects which shape the multisteps thermal denaturation. By summing over ~10(7)-10(8) base pair paths, a large number of double helix configurations is taken into account consistently with the physical requirements of the model potential. The partition function is computed as a function of the twist. It is found that the equilibrium twist angle, peculiar of B-DNA at room temperature, yields the stablest helicoidal geometry against thermal disruption of the base pair hydrogen bonds. This result is corroborated by the computation of thermodynamical properties such as fractions of open base pairs and specific heat.     

DNA and estrogen receptor interaction revealed by fragment molecular orbital calculations

Molecular orbital calculations of the complex between DNA-ERE (estrogen response element) and ER (estrogen receptor)-DBD (DNA-binding domain) were performed using the fragment molecular orbital (FMO) method, which enables large-scale MO (molecular orbital) calculations by reducing the computational cost and by significantly increasing efficiency for parallel computation. Such a large system, which contains 3354 atoms, is impractical via conventional MO methods due to the immense computational cost. Details of the interaction between DNA-ERE and ER-DBD were revealed in this study as follows by using the FMO calculations to analyze the interfragment interaction energies (IFIEs) and the electrostatic potentials (ESPs). An area with a high positive ESP is identified on the DNA-binding side of ER-DBD and is the main driving force behind access to the DNA. The position of the ER-DBD monomer can be fixed on a phosphate group of DNA-ERE by the strong electrostatic interactions, whereas the rotation cannot be fixed. In contrast, both the position and rotation of the ER-DBD dimer can be fixed and can therefore form the stable (ER-DBD)2...DNA-ERE complex. Dimerization of the ER-DBD monomers, each of which have a charge of +5 , is mainly due to large attractive interaction energies of the second Zn fragments. The base pairs in the consensus sequence of DNA-ERE interact only with the recognition helix located in the major groove due to the large shielding effect of the phosphate groups of DNA. The recognition helix has weaker interactions with the base pairs than the electrostatic interactions with the phosphate groups. Thus, the DNA-binding machinery of the ER-DBD dimer, which can secure the recognition helix in the major groove of DNA, is crucial for interactions between the recognition helix and base pairs.     

On the stability of single- and double-stranded DNA helices: The application of the PPT-MCF method on large fragments of DNA

The pseudo-polarization tensor mutually consistent field (PPT -MCF ) method recently introduced [1] has been applied to study the stacking interactions between the nucleotide bases in large periodic B-DNA fragments. The effects on the global and local binding properties caused by replacing one base in the periodic sequence by another base are investigated. The increase in the stability for comparable fragments owing to this base substitution is further enforced in the case of periodic alternating helices. The most important results are that the stacking interaction between two bases is slowly converging with the interbase distance and that the average contribution per base to the binding energy is repulsive. Furthermore, the energetical properties of double helix models in B- and Z-DNA configurations, respectively, consisting of up to five base pairs have been compared. It turns out that the G C G C sequence in Z-DNA is significantly more stable than either in periodic or periodic alternating B-DNA. In these cases the average energy contribution of a single Watson–Crick-type base pair is predicted also to be positive. From the calculations it follows that the double helix is not stabilized owing to the hydrogen bonding between the bases belonging to both strands, in contradiction to most other investigations.     

Melting temperature measurement and mesoscopic evaluation of single,double and triple DNA mismatches

Unlike the canonical base pairs AT and GC, the molecular properties of mismatches such as hydrogen bonding and stacking interactions are strongly dependent on the identity of the neighbouring base pairs. As a result, due to the sheer number of possible combinations of mismatches and flanking base pairs, only a fraction of these have been studied in varying experiments or theoretical models. Here, we report on the melting temperature measurement and mesoscopic analysis of contiguous DNA mismatches in nearest-neighbours and next-nearest neighbour contexts. A total of 4032 different mismatch combinations, including single, double and triple mismatches were covered. These were compared with 64 sequences containing all combinations of canonical base pairs in the same location under the same conditions. For a substantial number of single mismatch configurations, 15%, the measured melting temperatures were higher than the least stable AT base pair. The mesoscopic calculation, using the Peyrard–Bishop model, was performed on the set of 4096 sequences, and resulted in estimates of on-site and nearest-neighbour interactions that can be correlated to hydrogen bonding and base stacking. Our results confirm many of the known properties of mismatches, including the peculiar sheared stacking of tandem GA mismatches. More intriguingly, it also reveals that a number of mismatches present strong hydrogen bonding when flanked on both sites by other mismatches. To highlight the applicability of our results, we discuss a number of practical situations such as enzyme binding affinities, thymine DNA glycosylase repair activity, and trinucleotide repeat expansions.

A comprehensive experimental and theoretical evaluation of all DNA mismatch contexts, providing an insight into the intra-molecular interactions.     

Solution Structure of a Hexitol Nucleic Acid Duplex with Four Consecutive T⋅T Base Pairs

The structure of the hexitol nucleic acid (HNA) h(GCGCTTTTGCGC) was determined by NMR spectroscopy. This unnatural nucleic acid was developed as a mimic for A‐RNA. In solution, the studied sequence is forming a symmetric double‐stranded structure with four central consecutive T⋅T wobble pairs flanked by G⋅C Watson‐Crick base pairs. The stem regions adopt an A‐type helical structure. Discrete changes in backbone angles are altering the course of the helix axis in the internal loop region. Two H‐bonds are formed in each wobble pair, and base stacking is preserved in the duplex, explaining the stability of the duplex. This structure elucidation provides information about the influence of a (T)₄ fragment on local helix geometries as well as on the nature of the T⋅T mismatch base pairing in a TTTT tract.     

Effective thermostat induced by coarse graining of simple point charge water

We investigate how the transport properties of a united atom fluid with a dissipative particle dynamics thermostat depend on the functional form and magnitude of both the conservative and the stochastic interactions. We demonstrate how the thermostat strongly affects the hydrodynamics, especially diffusion, viscosity, and local escape times. As model system we use simple point charge (SPC) water, from which projected trajectories are used to determine the effective interactions in the united atom model. The simulation results support our argument that the thermostat should be viewed as an integral part of the coarse-grained dynamics rather than a tool for approaching thermal equilibrium. As our main result we show that the united atom model with the adjusted effective interactions approximately reproduces the diffusion constant and the viscosity of the underlying detailed SPC water model.     

MD simulation of the transitions between B-DNA and A-DNA in the framework of a coarse-grained model

Transitions between the B and A forms of a short DNA double helix (12 base pairs) at different salt concentrations in an aqueous solution have been studied by the molecular dynamics method in the framework of a coarse-grained model with explicit ions but without friction. It has been shown that the A-DNA, stable at high salt concentrations, is a dynamic conglomerate of the molecule and the ions coming from the solution into the deep major groove and then leaving it. In such a short helix, in the model without friction, even at low salt concentrations, transitions from B-DNA to A-DNA and back are frequent and fast. Stable ADNA (without transitions to B-DNA) forms at salt concentrations greater than 0.45 mol/L.     

Conserved patterns in backbone torsional changes allow for single base flipping from duplex DNA with minimal distortion of the double helix

Base flipping is a structural mechanism common to many DNA processing and repair enzymes. Changes in the local backbone torsions that occur during base flipping and the effect of environment on their behavior are of particular interest in understanding different base flipping mechanisms. In the present study, structures sampled during umbrella sampling molecular dynamics (MD) simulations of base flipping in aqueous and protein-bound environments, carried out with two different MD simulation strategies, are analyzed to find the most significant phosphodiester backbone distortions in the vicinity of the flipping base. Torsional sampling on the 5' side of the flipping base during flipping through the major groove shows similarities to the torsional sampling on the 3' side during flipping through the minor groove and vice versa. In differing environments, this behavior varies only marginally. These compensating torsional changes in the DNA backbone on 5' and 3' sides of the flipping base limit overall distortion of the DNA double helix during single base flipping. Rotameric intermediate states observed during base flipping are identified and postulated to be metastable states implicated in both large-scale structural changes and functional effects of chemical modifications in DNA.     

Coarse-grained computations for a micellar system

We establish, through coarse-grained computation, a connection between traditional, continuum numerical algorithms (initial value problems as well as fixed point algorithms), and atomistic simulations of the Larson model of micelle formation. The procedure hinges on the (expected) evolution of a few slow, coarse-grained mesoscopic observables of the Monte Carlo simulation, and on (computational) time scale separation between these and the remaining "slaved," fast variables. Short bursts of appropriately initialized atomistic simulation are used to estimate the (coarse grained, deterministic) local dynamics of the evolution of the observables. These estimates are then in turn used to accelerate the evolution to computational stationarity through traditional continuum algorithms (forward Euler integration, Newton-Raphson fixed point computation). This "equation-free" framework, bypassing the derivation of explicit, closed equations for the observables (e.g., equations of state), may provide a computational bridge between direct atomistic/stochastic simulation and the analysis of its macroscopic, system-level consequences.     

Analysis of stacking overlap in nucleic acid structures: algorithm and application

RNA contains different secondary structural motifs like pseudo-helices, hairpin loops, internal loops, etc. in addition to anti-parallel double helices and random coils. The secondary structures are mainly stabilized by base-pairing and stacking interactions between the planar aromatic bases. The hydrogen bonding strength and geometries of base pairs are characterized by six intra-base pair parameters. Similarly, stacking can be represented by six local doublet parameters. These dinucleotide step parameters can describe the quality of stacking between Watson–Crick base pairs very effectively. However, it is quite difficult to understand the stacking pattern for dinucleotides consisting of non canonical base pairs from these parameters. Stacking interaction is a manifestation of the interaction between two aromatic bases or base pairs and thus can be estimated best by the overlap area between the planar aromatic moieties. We have calculated base pair overlap between two consecutive base pairs as the buried van der Waals surface between them. In general, overlap values show normal distribution for the Watson–Crick base pairs in most double helices within a range from 45 to 50 Å² irrespective of base sequence. The dinucleotide steps with non-canonical base pairs also are seen to have high overlap value, although their twist and few other parameters are rather unusual. We have analyzed hairpin loops of different length, bulges within double helical structures and pseudo-continuous helices using our algorithm. The overlap area analyses indicate good stacking between few looped out bases especially in GNRA tetraloop, which was difficult to quantitatively characterise from analysis of the base pair or dinucleotide step parameters. This parameter is also seen to be capable to distinguish pseudo-continuous helices from kinked helix junctions.     

α-Synuclein induces both positive mean curvature and negative Gaussian curvature in membranes

Using a combination of X-ray scattering, fluorescence correlation spectroscopy, coarse-grained molecular dynamics (MD) simulations and potential of mean force calculations, we have explored the membrane remodeling effects of monomeric α-synuclein (αS). Our initial findings from multiple approaches are that αS (1) causes a significant thinning of the bilayer and (2) stabilizes positive mean curvature, such that the maximum principle curvature matches that of synaptic vesicles, αS-induced tubules, and the synthetic lipid vesicles to which the protein binds most tightly. This suggests that αS binding to synaptic vesicles likely stabilizes their intrinsic curvature. We then show that αS induces local negative Gaussian curvature, an effect that occurs in regions of αS shown previously via NMR and corroborated by MD simulation to have significant conformational flexibility. The induction of negative Gaussian curvature, which has implications for all curvature-sensing and curvature-generating amphipathic α-helices, supports a hypothesis that connects helix insertion to fusion and fission of vesicles, processes that have recently been linked to αS function. Then, in an effort to explain these biophysical properties of αS, we promote an intrinsic curvature-field model that recasts long-range protein-protein interactions in terms of the interactions between the local curvature fields generated by lipid-protein complexes.