Helical Parameters Research Articles

Effect of the sulfur atom on S2 and S4 positions of the uracil ring in different DNA:RNA hybrid microhelixes with three nucleotide base pairs

The effect of the sulphur atom on the uracil ring was analyzed in different DNA:RNA microhelixes with three nucleotide base-pairs, including uridine, 2-thiouridine, 4-thiouridine, 2,4-dithiouridine, cytidine, adenosine and guanosine. Distinct backbone and helical parameters were optimized at different density functional (DFT) levels. The Watson-Crick pair with 2-thiouridine appears weaker than with uridine, but its interaction with water molecules appears easier. Two types of microhelixes were found, depending on the H-bond of H2' hydroxyl atom: A-type appears with the ribose ring in 3 E-envelope C3' -endo, and B-type in 2 E-envelope C2' -endo. B-type is less common but it is more stable and with higher dipole-moment. The sulphur atoms significantly increase the dipole-moment of the microhelix, as well as the rise and propeller twist parameters. Simulations with four Na atoms H-bonded to the phosphate groups, and further hydration with explicit water molecules were carried out. A re-definition of the numerical value calculation of several base-pair and base-stacking parameters is suggested.

Functional reconstitution of the type IVa pilus assembly system from enterohaemorrhagic Escherichia coli.

Type 4a pili (T4aP) are long, thin and dynamic fibres displayed on the surface of diverse bacteria promoting adherence, motility and transport functions. Genomes of many Enterobacteriaceae contain conserved gene clusters encoding putative T4aP assembly systems. However, their expression has been observed only in few strains including Enterohaemorrhagic Escherichia coli (EHEC) and their inducers remain unknown. Here we used EHEC genomic DNA as a template to amplify and assemble an artificial operon composed of four gene clusters encoding 13 pilus assembly proteins. Controlled expressions of this operon in nonpathogenic E. coli strains led to efficient assembly of T4aP composed of the major pilin PpdD, as shown by shearing assays and immunofluorescence microscopy. When compared with PpdD pili assembled in a heterologous Klebsiella T2SS type 2 secretion system (T2SS) by using cryo-electron microscopy (cryoEM), these pili showed indistinguishable helical parameters, emphasizing that major pilins are the principal determinants of the fibre structure. Bacterial two-hybrid analysis identified several interactions of PpdD with T4aP assembly proteins, and with components of the T2SS that allow for heterologous fibre assembly. These studies lay ground for further characterization of the T4aP structure, function and biogenesis in enterobacteria.

An In-Depth Look at DNA Crystals through the Prism of Molecular Dynamics Simulations

Summary X-ray crystallography is the primary tool for biomolecular structural determination. However, contacts formed through the crystal lattice are known to affect structures, especially for small and flexible molecules such as DNA oligomers, by introducing significant structural changes in comparison to solution. Furthermore, why molecules crystallize in certain symmetry groups, which role crystallization additives play, and whether they are just innocuous and unspecific crystallization catalysts remain unclear. By using one of the currently best-performing DNA force fields and applying significant computational effort, we described the nature of intermolecular forces that stabilize B-DNA crystals in various symmetry groups and solvent environments with an unprecedented level of detail. We showed a tight coupling between the lattice stability and the type of crystallization additives and that certain symmetry groups are stable only in the presence of a specific additive. Additives and crystal contacts induce small but non-negligible changes in the physical properties of DNA.

Characterization of the Search Complex and Recognition Mechanism of the AlkD-DNA Glycosylase.

DNA damage is a routine problem for cells, and pathways such as base excision repair have evolved to protect the genome by using DNA glycosylases to first recognize and excise lesions. The search mechanism of these enzymes is of particular interest due to the seemingly intractable problem of probing the billions of base pairs in the genome for potential damage. It has been hypothesized that glycosylases form multiple protein-DNA conformational states to efficiently search and recognize DNA lesions, ultimately only flipping out the damaged substrate into the active site. A unique DNA glycosylase, the Bacillus cereus AlkD enzyme, has been shown to excise damaged DNA without flipping the nucleobase into a protein binding pocket following lesion recognition. Here, we use microsecond-scale all-atom molecular dynamics simulations to characterize the AlkD recognition mechanism, putting it in perspective with other DNA glycosylases. We first identify and describe two distinct enzyme-DNA conformations of AlkD: the search complex (SC) and excision complex (EC). The SC is distinguished by the linearity of DNA, changes in four helical parameters in the vicinity of the lesion, and changes in distance between active site residues and the DNA. Free DNA simulations are used to demonstrate that the DNA structural deviations and increased active site interactions present in the EC are initiated by the recognition of a methylation-induced signal in the rises both 5' to the methylation and opposing this base. Our results support the hypothesis that subtle geometric distortions in DNA are recognized by AlkD and are consequently probed to initiate concerted protein and DNA conformational changes which prime excise without additional intermediate states. This mechanism is shown to be consistent among the three methylated DNA sequences that have been crystallized bound to AlkD.

Structural insight into antisense gapmer-RNA oligomer duplexes through molecular dynamics simulations

There is an extensive research carrying out on antisense technology and the molecules entering into clinical trials are increasing rapidly. Phosphorothioate (PS) is a chemical modification in which nonbridged oxygen is replaced with a sulfur, consequently providing resistance against nuclease activity. The 2'-4' conformationally restricted nucleoside has the structural features of both 2'-O-methoxy ethyl RNA (MOE), which shows good toxicity profile, and locked nucleic acid (LNA), which shows good binding affinity towards the target RNA. These modifications have been studied and suggested that they can be a potential therapeutic agents in antisense therapy. Mipomersen (ISIS 301012), which contains the novel nucleoside modification has been used to target to apolipoprotein (Apo B), which reduces LDL cholesterol by 6–41%. In this study, classical molecular dynamics (MD) simulations were performed on six different antisense gapmer/target-RNA oligomer duplexes (LNA-PS-LNA/RNA, RcMOE-PS-RcMOE/RNA, ScMOE-PS-ScMOE/RNA, MOE-PS-MOE/RNA, PS-DNA/RNA and DNA/RNA) to investigate the structural dynamics, stability and solvation properties. The LNA, MOE nucleotides present in respective duplexes are showing the structure of A-form and the PS-DNA nucleotides resemble the structure of B-form helix with respect to some of the helical parameters. Free energy calculations suggest that the oligomer, which contains LNA binds to the RNA strongly than other modifications as shown in experimental results. The MOE modified nucleotide, which although had a lower binding affinity but higher solvent accessible surface area (SASA) compared to the other modifications, may be influencing the toxicity and hence may be used it in Mipomersen, the second antisense molecule which is approved by FDA.Communicated by Ramaswamy H. Sarma

Oxidative damage to epigenetically methylated sites affects DNA stability, dynamics and enzymatic demethylation.

DNA damage can affect various regulatory elements of the genome, with the consequences for DNA structure, dynamics, and interaction with proteins remaining largely unexplored. We used solution NMR spectroscopy, restrained and free molecular dynamics to obtain the structures and investigate dominant motions for a set of DNA duplexes containing CpG sites permuted with combinations of 5-methylcytosine (mC), the primary epigenetic base, and 8-oxoguanine (oxoG), an abundant DNA lesion. Guanine oxidation significantly changed the motion in both hemimethylated and fully methylated DNA, increased base pair breathing, induced BI→BII transition in the backbone 3′ to the oxoG and reduced the variability of shift and tilt helical parameters. UV melting experiments corroborated the NMR and molecular dynamics results, showing significant destabilization of all methylated contexts by oxoG. Notably, some dynamic and thermodynamic effects were not additive in the fully methylated oxidized CpG, indicating that the introduced modifications interact with each other. Finally, we show that the presence of oxoG biases the recognition of methylated CpG dinucleotides by ROS1, a plant enzyme involved in epigenetic DNA demethylation, in favor of the oxidized DNA strand. Thus, the conformational and dynamic effects of spurious DNA oxidation in the regulatory CpG dinucleotide can have far-reaching biological consequences.

Correlation between translational and rotational kinematic abnormalities and osteoarthritis-like damage in two in vivo sheep injury models

The relations between kinematic abnormalities and post traumatic osteoarthritis have not yet been clearly elucidated. This study was conducted to determine the finite helical axes parameters and the tibiofemoral translation vector in the knee joints of two surgically induced injury sheep models: anterior cruciate ligament and medial collateral ligament transection (ACL/MCL Tx) (n = 5) and lateral meniscectomy (n = 5). We hypothesized that morphological damage in the experimental joints would be correlated to alterations in these kinematic variables. There was no strong evidence that morphological damage to the joints 20 weeks post ACL/MCL transection or meniscectomy was correlated with alterations in the finite helical axes variables. Nevertheless, significant correlations were found between the morphological damage to the joints and the magnitude of the change in the translation vectors after ACL/MCL transection (significant correlations (p = 0.005) during stance and trends (p < 0.1) at all points analyzed during swing). It can be concluded that: (1) osteoarthritic-like morphological damage after ACL/MCL transection is more critically correlated to the absolute tibiofemoral translational change and (2) alterations in analyzed kinematic variables cannot solely define osteoarthritis risk after meniscal injuries. From a clinical perspective, our results suggest that the magnitude of the change in the translation vector, which is independent of the coordinate system and combines the effects of the three translational degrees of freedom, i.e. medial–lateral, anterior-posterior and inferior-superior, would be an osteoarthritis risk factor after ligament injury, and requires validation in humans.

A Molecular Dynamics-Quantum Mechanics Theoretical Study of DNA-Mediated Charge Transport in Hydrated Ionic Liquids

Charge transport (CT) through biomolecules is of high significance in the research fields of biology, nanotechnology, and molecular devices. Inspired by our previous work that showed the binding of ionic liquid (IL) facilitated charge transport in duplex DNA, in silico simulation is a useful means to understand the microscopic mechanism of the facilitation phenomenon. Here molecular dynamics simulations (MD) of duplex DNA in water and hydrated ionic liquids were employed to explore the helical parameters. Principal component analysis was further applied to capture the subtle conformational changes of helical DNA upon different environmental impacts. Sequentially, CT rates were calculated by a QM/MM simulation of the flickering resonance model based upon MD trajectories. Herein, MD simulation illustrated that the binding of ionic liquids can restrain dynamic conformation and lower the on-site energy of the DNA base. Confined movement among the adjacent base pairs was highly related to the increase of electronic coupling among base pairs, which may lead DNA to a CT facilitated state. Sequentially combining MD and QM/MM analysis, the rational correlations among the binding modes, the conformational changes, and CT rates illustrated the facilitation effects from hydrated IL on DNA CT and supported a conformational-gating mechanism.

Unique Conformation in a Natural Interruption Sequence of Type XIX Collagen Revealed by Its High-Resolution Crystal Structure.

Naturally occurring interruptions in nonfibrillar collagen play key roles in molecular flexibility, collagen degradation, and ligand binding. The structural feature of the interruption sequences and the molecular basis for their functions have not been well studied. Here, we focused on a G5G type natural interruption sequence G-POALO-G from human type XIX collagen, a homotrimer collagen, as this sequence possesses distinct properties compared with those of a pathological similar Gly mutation sequence in collagen mimic peptides. We determined the crystal structures of the host-guest peptide (GPO)3-GPOALO-(GPO)4 to 1.03 Å resolution in two crystal forms. In these structures, the interruption zone brings localized disruptions to the triple helix and introduces a light 6-8° bend with the same directional preference to the whole molecule, which may correspond structurally to the first physiological kink site in type XIX collagen. Furthermore, at the G5G interruption site, the presence of Ala and Leu residues, both with free N-H groups, allows the formation of more direct and water-mediated interchain hydrogen bonds than in the related Gly → Ala structure. These could partly explain the difference in thermal stability between the different interruptions. In addition, our structures provide a detailed view of the dynamic property of such an interrupted zone with respect to hydrogen bonding topology, torsion angles, and helical parameters. Our results, for the first time, also identified the binding of zinc to the end of the triple helix. These findings will shed light on how the interruption sequence influences the conformation of the collagen molecule and provide a structural basis for further functional studies.

Influence Mechanisms Between Helical Motion Velocity Coefficient and Cutter Parameters and Characteristics of Tooth Surface Generated by Duplex Helical Method

Influence Mechanisms Between Helical Motion Velocity Coefficient and Cutter Parameters and Characteristics of Tooth Surface Generated by Duplex Helical Method

Structural Basis for the Influence of A1, 5A, and W51W57 Mutations on the Conductivity of the Geobacter sulfurreducens Pili

The metallic-like conductivity of the Geobacter sulfurreducens pilus and higher conductivity of its mutants reflected that biological synthesis can be utilized to improve the properties of electrically conductive pili. However, the structural basis for diverse conductivities of nanowires remains uncertain. Here, the impacts of point mutations on the flexibility and stability of pilins were investigated based on molecular dynamics simulations. Structures of the G. sulfurreducens pilus and its mutants were constructed by Rosetta. Details of the structure (i.e., electrostatic properties, helical parameters, residue interaction network, distances between amino acids, and salt bridges) were analyzed by PDB2PQR, Rosetta, RING, PyMOL, and VMD, respectively. Changes in stability, flexibility, residue interaction, and electrostatic properties of subunits directly caused wild-type pilin and its mutants assemble different structures of G. sulfurreducens pili. By comparing the structures of pili with different conductivities, the mechanism by which the G. sulfurreducens pilus transfers electron along pili was attributed, at least in part, to the density of aromatic rings, the distances between neighboring aromatic rings, and the local electrostatic environment around aromatic contacts. These results provide new insight into the potential for the biological synthesis of highly electrically conductive, nontoxic nanowires.

Overstretching of B-DNA with various pulling protocols: Appearance of structural polymorphism and S-DNA.

We report a structural polymorphism of the S-DNA when a canonical B-DNA is stretched under different pulling protocols and provide a fundamental molecular understanding of the DNA stretching mechanism. Extensive all atom molecular dynamics simulations reveal a clear formation of S-DNA when the B-DNA is stretched along the 3' directions of the opposite strands (OS3) and is characterized by the changes in the number of H-bonds, entropy, and free energy. Stretching along the 5' directions of the opposite strands (OS5) leads to force induced melting form of the DNA. Interestingly, stretching along the opposite ends of the same strand leads to a coexistence of both the S- and melted M-DNA structures. We also do the structural characterization of the S-DNA by calculating various helical parameters. We find that the S-DNA has a twist of ∼10° which corresponds to a helical repeat length of ∼36 base pairs in close agreement with the previous experimental results. Moreover, we find that the free energy barrier between the canonical and overstretched states of DNA is higher for the same termini pulling protocol in comparison to all other protocols considered in this work. Overall, our observations not only reconcile with the available experimental results qualitatively but also enhance the understanding of different overstretched DNA structures.

The Cryoelectron Microscopy Structure of the Type 1 Chaperone-Usher Pilus Rod

SummaryAdhesive chaperone-usher pili are long, supramolecular protein fibers displayed on the surface of many bacterial pathogens. The type 1 and P pili of uropathogenic Escherichia coli (UPEC) play important roles during urinary tract colonization, mediating attachment to the bladder and kidney, respectively. The biomechanical properties of the helical pilus rods allow them to reversibly uncoil in response to flow-induced forces, allowing UPEC to retain a foothold in the unique and hostile environment of the urinary tract. Here we provide the 4.2-Å resolution cryo-EM structure of the type 1 pilus rod, which together with the previous P pilus rod structure rationalizes the remarkable “spring-like” properties of chaperone-usher pili. The cryo-EM structure of the type 1 pilus rod differs in its helical parameters from the structure determined previously by a hybrid approach. We provide evidence that these structural differences originate from different quaternary structures of pili assembled in vivo and in vitro.

Molecular Dynamics of Double Stranded Xylo-Nucleic Acid.

Xylo-nucleic acid (XyloNA) is a synthetic analogue of ribo-nucleic acid (RNA), where the ribose sugar has been replaced by xylose. We present a molecular dynamics study of the conformational evolution of XyloNA double strand oligomers derived from A-RNA through the substitution of β-d-ribofuranose by β-d-xylofuranose and having lengths of 8, 16, and 29 base pairs, using a set of independent all-atom simulations performed at various time scales ranging from 55 to 100 ns, with one long 500 ns simulation of the 29-mer. In order to validate the robustness of XyloNA conformation, a set of simulations using various cutoff distances and solvation box dimensions has also been performed. These independent simulations reveal the uncoiling or elongation of the initial conformation to form an open ladder type transient state conformation and the subsequent formation of a highly flexible duplex with a tendency to coil in a left-handed fashion. The observed open ladder conformation is in line with recently obtained NMR data on the XyloNA 8-mer derived using 5'-d(GUGUACAC)-3'. The observed negative interbase pair twist leads to the observed highly flexible left-handed duplex, which is significantly less rigid than the stable left-handed dXyloNA duplex having a strong negative twist. A comparison between the xylo-analogues of DNA and RNA shows a clear distinction between the helical parameters, with implications for the pairing mechanism.

Cryoelectron Microscopy Reconstructions of the Pseudomonas aeruginosa and Neisseria gonorrhoeae Type IV Pili at Sub-nanometer Resolution

We report here cryoelectron microscopy reconstructions of type IV pili (T4P) from two important human pathogens, Pseudomonas aeruginosa and Neisseria gonorrhoeae, at ∼ 8 and 5Å resolution, respectively. The two structures reveal distinct arrangements of the pilin globular domains on the pilus surfaces, which impart different helical parameters, but similar packing of the conserved N-terminal α helices, α1, in the filament core. In contrast to the continuous α helix seen in the X-ray crystal structures of the P.aeruginosa and N.gonorrhoeae pilin subunits, α1 in the pilus filaments has a melted segment located between conserved helix-breaking residues Gly14 and Pro22, as seen for the Neisseria meningitidis T4P. Using mutagenesis we show that Pro22 is critical for pilus assembly, as are Thr2 and Glu5, which are positioned to interact in the hydrophobic filament core. These structures provide a framework for understanding T4P assembly, function, and biophysical properties.

CgHeliParm: analysis of dsDNA helical parameters for coarse-grained MARTINI molecular dynamics simulations.

We introduce cgHeliParm, a python program that provides the conformational analysis of Martini-based coarse-grained double strand DNA molecules. The software calculates the helical parameters such as base, base pair and base pair step parameters. cgHeliParm can be used for the analysis of coarse grain Martini molecular dynamics trajectories without transformation into atomistic models. This package works with Python 2.7 on MacOS and Linux. The program is freely available for download from https://github.com/ifaust83/cgheliparm. Together with the main script, the base reference files CG_X_std.lib, a number of examples and R scripts are also available from the same website. A tutorial on the use and application is also available at http://cgmartini.nl/index.php/tutorials-general-introduction/tutorial-martini-dna. i.faustino@rug.nl. Supplementary data are available at Bioinformatics online.

Electronic Structures of LNA Phosphorothioate Oligonucleotides.

Important oligonucleotides in anti-sense research have been investigated in silico and experimentally. This involves quantum mechanical (QM) calculations and chromatography experiments on locked nucleic acid (LNA) phosphorothioate (PS) oligonucleotides. iso-potential electrostatic surfaces are essential in this study and have been calculated from the wave functions derived from the QM calculations that provide binding information and other properties of these molecules. The QM calculations give details of the electronic structures in terms of e.g., energy and bonding, which make them distinguish or differentiate between the individual PS diastereoisomers determined by the position of sulfur atoms. Rules are derived from the electronic calculations of these molecules and include the effects of the phosphorothioate chirality and formation of electrostatic potential surfaces. Physical and electrochemical descriptors of the PS oligonucleotides are compared to the experiments in which chiral states on these molecules can be distinguished. The calculations demonstrate that electronic structure, electrostatic potential, and topology are highly sensitive to single PS configuration changes and can give a lead to understanding the activity of the molecules.

Effect of the sulphur atom on geometry and spectra of the biomolecule 2-thiouracil and in the WC base pair 2-thiouridine-adenosine. Influence of water in the first hydration shell

The effect of the sulphur atom on 2-thiouracil (2TU) and 2-thiouridine molecules, as compared with uracil and uridine molecules, respectively, was carried out in several environments. The predicted IR spectrum of 2TU in the isolated state was compared with that obtained for uracil molecule and with those reported experimentally in matrix isolation. Its crystal unit cell in the solid state was simulated through a tetramer form using DFT methods for the first time. The calculated Raman spectrum was compared to the experimental ones in the solid state. A linear scaling procedure was used for this task. The first hydration shell was simulated by explicit number of water molecules surrounding 2TU up to 30 and was compared with that obtained in uracil molecule. Water molecules ‘distributed’ around 2TU was preferred over that ‘clustering’, because it can better reproduce the hydration and their effects on different parameters of the molecular structure of 2TU and uracil. The total atomic charges and several calculated thermodynamic parameters were discussed. The effect of the sulphur atom on the Watson-Crick (WC) and reverse WC base pair uridine-adenosine was estimated, and the CP corrected interaction energies were calculated. 2-thiouridine has a weaker WC pair than that with uridine, although its slight higher dipole moment (μ) facilitates the interaction with the water molecules. Several helical parameters were determined.

Right-Handed Helical Foldamers Consisting of De Novo d-AApeptides

New types of foldamer scaffolds are formidably challenging to design and synthesize, yet highly desirable as structural mimics of peptides/proteins with a wide repertoire of functions. In particular, the development of peptidomimetic helical foldamers holds promise for new biomaterials, catalysts, and drug molecules. Unnatural l-sulfono-γ-AApeptides were recently developed and shown to have potential applications in both biomedical and material sciences. However, d-sulfono-γ-AApeptides, the enantiomers of l-sulfono-γ-AApeptides, have never been studied due to the lack of high-resolution three-dimensional structures to guide structure-based design. Herein, we report the first synthesis and X-ray crystal structures of a series of 2:1 l-amino acid/d-sulfono-γ-AApeptide hybrid foldamers, and elucidate their folded conformation at the atomic level. Single-crystal X-ray crystallography indicates that this class of oligomers folds into well-defined right-handed helices with unique helical parameters. The helical structures were consistent with data obtained from solution 2D NMR, CD studies, and molecular dynamics simulations. Our findings are expected to inspire the structure-based design of this type of unique folding biopolymers for biomaterials and biomedical applications.

Structure of RNA Stem Loop B from the Picornavirus Replication Platform.

The presumptive RNA cloverleaf at the start of the 5'-untranslated region of the picornavirus genome is an essential element in replication. Stem loop B (SLB) of the cloverleaf is a recognition site for the host polyC-binding protein, which initiates a switch from translation to replication. Here we present the solution structure of human rhinovirus isotype 14 SLB using nuclear magnetic resonance spectroscopy. SLB adopts a predominantly A-form helical structure. The stem contains five Watson-Crick base pairs and one wobble base pair and is capped by an eight-nucleotide loop. The wobble base pair introduces perturbations into the helical parameters but does not appear to introduce flexibility. However, the helix major groove appears to be accessible. Flexibility is seen throughout the loop and in the terminal nucleotides. The pyrimidine-rich region of the loop, the apparent recognition site for the polyC-binding protein, is the most disordered region of the structure.