How sequence defines structure: a crystallographic map of DNA structure and conformation

Evaluation of mutations on O6-methylguanine methyl transferase structure and its interactions: molecular dynamics simulation study

O6-methylguanine DNA methyl transferase (MGMT) is a significant vehicle for the cellular clearance of alkyl lesions, particularly the methyl group of the O-6 and O-4 positions of guanine and thymine, respectively. Many publications have studied the correlation between polymorphisms in MGMT and susceptibility to various cancers. In the present study, we investigated the consequence of L84F, common single-nucleotide polymorphism, K125E, site-specific mutagenesis, and L84F/K125E on conformation, stability, and behavior of MGMT in the free form and interaction with proliferating cell nuclear antigen (PCNA) and DNA as partners in the biochemical network by using molecular dynamics simulation method. Our results showed that all free variants of MGMT differed from the native form. However, among all free variants of MGMT, the L84F/K125E variant exhibited similar properties compared with the wild-type. In contrast, in complex modes, only amino acid residues of the L84F variant are involved in the interactions with PCNA and DNA somewhat differently relative to the wild-type. Furthermore, L84F SNP showed the highest binding free energy compared to other variants and native forms. These alterations in the amino acids and binding free energy of L84F relative to the native are the reasons for changing its region connection compared to the native form. Therefore, we propose conducting further investigations into the impact of inhibitors or chemotherapeutic agents to assess their effectiveness on MGMT variants compared to the wild-type, aiming to reduce the cost of cancer treatment that will depend on inhibiting native MGMT protein.

Molecular dynamics simulations revealed structural differences among TTHA1873-DNA interaction

This study aims to investigate the comparative binding pattern of TTHA1873 and its mutants (R55A and R138A) with DNA through molecular docking and molecular dynamics (MD) simulations. The docking results suggests that the Wild type (WT-TTHA1873), R55A, R138A and double mutant R55A/R138A having docking scores of -225.80 kcal/mol, -209.81 kcal/mol, -197.53 kcal/mol, -195.55 kcal/mol respectively and WT-TTHA1873 has more significant binding capability with DNA in comparison to mutants. The MD analysis revealed that the WT-TTHA1873 demonstrated stable interactions with DNA and exhibited a reduced conformational space compared to the mutants. By examining the atomic interactions, it was observed that significant variations in the hydrogen bonding pattern between WT-TTHA1873 and its mutants while interacting with DNA resulted in structural anomalies in the mutants and differences in DNA-binding specificity. The calculated binding free energies imply more stability of the WT-TTHA1873-DNA complex, while the mutants showed lesser binding affinity toward its interacting partner, double-stranded DNA. It is apparent that substituting single mutation R55A and R138A on TTHA1873 abolishes their DNA-binding ability. The present study portrays the critical role of R55 and R138 from TTHA1873 as likely involved in DNA binding.

Neurotransmitters and neural hormone-based probes for quadruplex DNA sequences associated with neurodegenerative diseases

The potential of neurotransmitters and neural hormones as possible G-quadruplex DNA binders was analyzed using fluorescence spectroscopy, surface-enhanced Raman spectroscopy (SERS), DNA melting analysis, and molecular docking. G-quadruplex sequences, (GGC)3 and G4C2, with roles in Fragile X syndrome and amyotrophic lateral sclerosis (ALS), respectively, were selected, and their interactions with melatonin, serotonin, and gamma-aminobutyric acid (GABA), were studied. Both melatonin and serotonin demonstrated strong interactions with the DNA sequences with hydrogen bonding being the primary mode of interaction, with some non-intercalative interactions involving the π systems. GABA demonstrated much weaker interactions and may not be a suitable candidate as a probe for low concentrations of G-quadruplex DNA.

Single-molecule structural and kinetic studies across sequence space

At the core of molecular biology lies the intricate interplay between sequence, structure, and function. Single-molecule techniques provide in-depth dynamic insights into structure and function, but laborious assays impede functional screening of large sequence libraries. We introduce high-throughput Single-molecule Parallel Analysis for Rapid eXploration of Sequence space (SPARXS), integrating single-molecule fluorescence with next-generation sequencing. We applied SPARXS to study the sequence-dependent kinetics of the Holliday junction, a critical intermediate in homologous recombination. By examining the dynamics of millions of Holliday junctions, covering thousands of distinct sequences, we demonstrated the ability of SPARXS to uncover sequence patterns, evaluate sequence motifs, and construct thermodynamic models. SPARXS emerges as a versatile tool for untangling the mechanisms that underlie sequence-specific processes at the molecular scale.

Quantum Mechanics Characterization of Non-Covalent Interaction in Nucleotide Fragments

Accurate calculation of non-covalent interaction energies in nucleotides is crucial for understanding the driving forces governing nucleic acid structure and function, as well as developing advanced molecular mechanics forcefields or machine learning potentials tailored to nucleic acids. Here, we dissect the nucleotides' structure into three main constituents: nucleobases (A, G, C, T, and U), sugar moieties (ribose and deoxyribose), and phosphate group. The interactions among these fragments and between fragments and water were analyzed. Different quantum mechanical methods were compared for their accuracy in capturing the interaction energy. The non-covalent interaction energy was decomposed into electrostatics, exchange-repulsion, dispersion, and induction using two ab initio methods: Symmetry-Adapted Perturbation Theory (SAPT) and Absolutely Localized Molecular Orbitals (ALMO). These calculations provide a benchmark for different QM methods, in addition to providing a valuable understanding of the roles of various intermolecular forces in hydrogen bonding and aromatic stacking. With SAPT, a higher theory level and/or larger basis set did not necessarily give more accuracy. It is hard to know which combination would be best for a given system. In contrast, ALMO EDA2 did not show dependence on theory level or basis set; additionally, it is faster.

Probing a Major DNA Weakness: Resolving the Groove and Sequence Selectivity of the Diimine Complex Λ-[Ru(phen)2 phi]2

The grooves of DNA provide recognition sites for many nucleic acid binding proteins and anticancer drugs such as the covalently binding cisplatin. Here we report a crystal structure showing, for the first time, groove selectivity by an intercalating ruthenium complex. The complex Λ-[Ru(phen)2 phi]2+ , where phi=9,10-phenanthrenediimine, is bound to the DNA decamer duplex d(CCGGTACCGG)2 . The structure shows that the metal complex is symmetrically bound in the major groove at the central TA/TA step, and asymmetrically bound in the minor groove at the adjacent GG/CC steps. A third type of binding links the strands, in which each terminal cytosine base stacks with one phen ligand. The overall binding stoichiometry is four Ru complexes per duplex. Complementary biophysical measurements confirm the binding preference for the Λ-enantiomer and show a high affinity for TA/TA steps and, more generally, TA-rich sequences. A striking enantiospecific elevation of melting temperatures is found for oligonucleotides which include the TATA box sequence.

Assessment of A- to B- DNA Transitions Utilizing the Drude Polarizable Force Field

In addition to the well-characterized B-form of DNA, duplex DNA can adopt various conformations, such as A or Z-DNA. Though less common, these structures can be induced biologically through protein or ligand interactions or experimentally with niche environmental conditions, such as high salt concentrations or in mixed water-ethanol. Reproducing these alternate structures through molecular dynamics simulations in recent years has been quite challenging with the currently available force fields, simulation techniques, and time scales. In this study, the Drude polarizable force field is tested for its ability to facilitate transitions between A-DNA and B-DNA or maintain A-DNA. Though transitions away from B-DNA were observed in high concentrations of ethanol, the resulting structures had hybrid properties taken from both B-DNA and A-DNA structures. This was also true for A-DNA in ethanol, which lost some of the A-DNA properties that it was expected to maintain. When B-DNA was tested in high salt environments, the resulting B-DNA structures showed no distinguishable differences with the increasing salt concentrations tested. These results with the Drude FF and recent results with additive force fields suggest that at present the current additive and polarizable force fields do not facilitate a complete transition between B- to A-DNA conformations under the conditions simulated. At present, the Drude FF favors A-B DNA hybrid structures when simulated in nonphysiological conditions.

Experimental detection of conformational transitions between forms of DNA: problems and prospects

Under different conditions, the DNA double helix can take different geometric forms. Of the large number of its conformations, in addition to the "canonical" B form, the A, C, and Z forms are widely known, and the D, Hoogsteen, and X forms are less known. DNA locally takes the A, C, and Z forms in the cell, in complexes with proteins. We compare different methods for detecting non-canonical DNA conformations: X-ray, IR, and Raman spectroscopy, linear and circular dichroism in both the infrared and ultraviolet regions, as well as NMR (measurement of chemical shifts and their anisotropy, scalar and residual dipolar couplings and inter-proton distances from NOESY (nuclear Overhauser effect spectroscopy) data). We discuss the difficulties in applying these methods, the problems of theoretical interpretation of the experimental results, and the prospects for reliable identification of non-canonical DNA conformations.

Upgrade of D+ software for hierarchical modeling of X-ray scattering data from complex structures in solution, fibers and single orientations

This article presents an upgrade of the D+ software [Ginsburg et al. (2019 ▸). J. Appl. Cryst. 52, 219-242], expanding its hierarchical solution X-ray scattering modeling capabilities for fiber diffraction and single crystallographic orientations. This upgrade was carried out using the reciprocal grid algorithm [Ginsburg et al. (2016 ▸). J. Chem. Inf. Model. 56, 1518-1527], providing D+ its computational strength. Furthermore, the extensive modifications made to the Python API of D+ are described, broadening the X-ray analysis performed with D+ to account for the effects of the instrument-resolution function and polydispersity. In addition, structure-factor and radial-distribution-function modules were added, taking into account the effects of thermal fluctuations and intermolecular interactions. Finally, numerical examples demonstrate the usage and potential of the added features.

Modular Protein-DNA Cocrystals as Precise, Programmable Assembly Scaffolds

High-precision nanomaterials to entrap DNA-binding molecules are sought after for applications such as controlled drug delivery and scaffold-assisted structural biology. Here, we engineered protein-DNA cocrystals to serve as scaffolds for DNA-binding molecules. The designed cocrystals, isoreticular cocrystals, contain DNA-binding protein and cognate DNA blocks where the DNA-DNA junctions stack end-to-end. Furthermore, the crystal symmetry allows topology preserving (isoreticular) expansion of the DNA stack without breaking protein-protein contacts, hence providing larger solvent channels for guest diffusion. Experimentally, the resulting designed isoreticular cocrystal adopted an interpenetrating I222 lattice, a phenomenon previously observed in metal-organic frameworks (MOFs). The interpenetrating lattice crystallized dependably in the same space group despite myriad modifications at the DNA-DNA junctions. Assembly was modular with respect to the DNA inserted for expansion, providing an interchangeable DNA sequence for guest-specified scaffolding. Also, the DNA-DNA junctions were tunable, accommodating varied sticky base overhang lengths and terminal phosphorylation. As a proof of concept, we used the interpenetrating scaffold crystals to separately entrap three distinct guest molecules during crystallization. Isoreticular cocrystal design offers a route to a programmable scaffold for DNA-binding molecules, and the design principles may be applied to existing cocrystals to develop scaffolding materials.

Efficient automatic construction of atom-economical QM regions with point-charge variation analysis

The setup of QM/MM calculations is not trivial since many decisions have to be made by the simulation scientist to achieve reasonable and consistent results. The main challenge to be tackled is the construction of the QM region to make sure to take into account all important parts of the adjacent environment and exclude less important ones. In our previous work [F. Brandt and Ch. R. Jacob, Systematic QM Region Construction in QM/MM Calculations Based on Uncertainty Quantification, J. Chem. Theory Comput., 2022, 18, 2584-2596.], we introduced the point charge variation analysis (PCVA) as a simple and reliable tool to systematically construct QM regions based on the sensitivity of the reaction energy with respect to variations of the MM point charges. Here, we assess several simplified variants of this PCVA approach for the example of catechol O-methyltransferase and apply PCVA for another system, the triosephosphate isomerase. Furthermore, we extend its scope by applying it to a DNA system. Our results indicate that PCVA offers an efficient and versatile approach of the automatic construction of atom-economical QM regions, but also identify possible pitfalls and limitations.

The Halogen Bond in Weakly Bonded Complexes and the Consequences for Aromaticity and Spin-Orbit Coupling

The halogen bond complexes CF3X⋯Y and C2F3X⋯Y, with Y = furan, thiophene, selenophene and X = Cl, Br, I, have been studied by using DFT and CCSD(T) in order to understand which factors govern the interaction between the halogen atom X and the aromatic ring. We found that PBE0-dDsC/QZ4P gives an adequate description of the interaction energies in these complexes, compared to CCSD(T) and experimental results. The interaction between the halogen atom X and the π-bonds in perpendicular orientation is stronger than the interaction with the in-plane lone pairs of the heteroatom of the aromatic cycle. The strength of the interaction follows the trend Cl < Br < I; the chalcogenide in the aromatic ring nor the hybridization of the C−X bond play a decisive role. The energy decomposition analysis shows that the interaction energy is dominated by all three contributions, viz., the electrostatic, orbital, and dispersion interactions: not one factor dominates the interaction energy. The aromaticity of the ring is undisturbed upon halogen bond formation: the π-ring current remains equally strong and diatropic in the complex as it is for the free aromatic ring. However, the spin-orbit coupling between the singlet and triplet π→π* states is increased upon halogen bond formation and a faster intersystem crossing between these states is therefore expected.

Polarizable Molecular Block Model: Toward the Development of an Induced Dipole Force Field for DNA

For flexible and highly ionized macromolecules such as DNA, it is important to correctly evaluate the intramolecular polarization in an induced dipole force field. In a proposed polarizable molecular block (PMB) model, a large molecule is divided into several molecular blocks. The atomic charges of the blocks are optimized by using the respective electrostatic potentials (ESPs) on the molecular surface. By using the capped hydrogen removal operation, the total charge of the blocks is controlled exactly to have an integer charge. The atomic polarizabilities of the blocks are optimized by using the respective polarized one-electron potentials that are the differences between ESPs with and without an external test charge. Induced dipole-charge interactions between the blocks are all included, but those interactions within the blocks are strictly excluded. All dipole-dipole interactions are included, but the damping functions are applied to the close dipole-dipole pairs. Several types of damping (simple scaling, exponential, linear, and Gaussian) are evaluated. The validity of the PMB model was verified by using trinucleotide duplexes which have A-, B-, and Z-DNA forms. The reference energies of trinucleotide duplexes including counterions (GGT3Na-ACC3Na, GAC3Na-GTC3Na, and GCG3Na-CGC3Na) are calculated using ωB97XD/aug-cc-pVDZ. All damping types reproduced well the reference interaction energies, dipole moments, and ESPs. Among them, the simple scaling with strong attenuation to 1-2 atomic pairs showed the highest stability against the polarization catastrophe. This study shows that it is possible to develop a high-quality polarizable force field by treating the intramolecular polarization on a block-by-block basis.

Interactions of small molecules with DNA junctions

The four natural DNA bases (A, T, G and C) associate in base pairs (A=T and G≡C), allowing the attached DNA strands to assemble into the canonical double helix of DNA (or duplex-DNA, also known as B-DNA). The intrinsic supramolecular properties of nucleobases make other associations possible (such as base triplets or quartets), which thus translates into a diversity of DNA structures beyond B-DNA. To date, the alphabet of DNA structures is ripe with approximately 20 letters (from A- to Z-DNA); however, only a few of them are being considered as key players in cell biology and, by extension, valuable targets for chemical biology intervention. In the present review, we summarise what is known about alternative DNA structures (what are they? When, where and how do they fold?) and proceed to discuss further about those considered nowadays as valuable therapeutic targets. We discuss in more detail the molecular tools (ligands) that have been recently developed to target these structures, particularly the three- and four-way DNA junctions, in order to intervene in the biological processes where they are involved. This new and stimulating chemical biology playground allows for devising innovative strategies to fight against genetic diseases.

Predicting accurate ab initio DNA electron densities with equivariant neural networks

One of the fundamental limitations of accurately modeling biomolecules like DNA is the inability to perform quantum chemistry calculations on large molecular structures. We present a machine learning model based on an equivariant Euclidean neural network framework to obtain accurate ab initio electron densities for arbitrary DNA structures that are much too large for conventional quantum methods. The model is trained on representative B-DNA basepair steps that capture both base pairing and base stacking interactions. The model produces accurate electron densities for arbitrary B-DNA structures with typical errors of less than 1%. Crucially, the error does not increase with system size, which suggests that the model can extrapolate to large DNA structures with negligible loss of accuracy. The model also generalizes reasonably to other DNA structural motifs such as the A- and Z-DNA forms, despite being trained on only B-DNA configurations. The model is used to calculate electron densities of several large-scale DNA structures, and we show that the computational scaling for this model is essentially linear. We also show that this machine learning electron density model can be used to calculate accurate electrostatic potentials for DNA. These electrostatic potentials produce more accurate results compared with classical force fields and do not show the usual deficiencies at short range.

Visualizing "Alternative Isoinformational Engineered" DNA in A- and B-Forms at High Resolution

A fundamental property of DNA built from four informational nucleotide units (GCAT) is its ability to adopt different helical forms within the context of the Watson-Crick pair. Well-characterized examples include A-, B-, and Z-DNA. For this study, we created an isoinformational biomimetic polymer, built (like standard DNA) from four informational "letters", but with the building blocks being artificial. This ALternative Isoinformational ENgineered (ALIEN) DNA was hypothesized to support two nucleobase pairs, the P:Z pair matching 2-amino-imidazo-[1,2a]-1,3,5-triazin-[8H]-4-one with 6-amino-3-5-nitro-1H-pyridin-2-one and the B:S pair matching 6-amino-4-hydroxy-5-1H-purin-2-one with 3-methyl-6-amino-pyrimidin-2-one. We report two structures of ALIEN DNA duplexes at 1.2 Å resolution and a third at 1.65 Å. All of these are built from a single self-complementary sequence (5'-CTSZZPBSBSZPPBAG) that includes 12 consecutive ALIEN nucleotides. We characterized the helical, nucleobase pair, and dinucleotide step parameters of ALIEN DNA in these structures. In addition to showing that ALIEN pairs retain basic Watson-Crick pairing geometry, two of the ALIEN DNA structures are characterized as A-form DNA and one as B-form DNA. We identified parameters that map differences effecting the transition between the two helical forms; these same parameters distinguish helical forms of isoinformational natural DNA. Collectively, our analyses suggest that ALIEN DNA retains essential structural features of natural DNA, not only its information density and Watson-Crick pairing but also its ability to adopt two canonical forms.

A unified computational view of DNA duplex, triplex, quadruplex and their donor-acceptor interactions

DNA can assume various structures as a result of interactions at atomic and molecular levels (e.g., hydrogen bonds, π–π stacking interactions, and electrostatic potentials), so understanding of the consequences of these interactions could guide development of ways to produce elaborate programmable DNA for applications in bio- and nanotechnology. We conducted advanced ab initio calculations to investigate nucleobase model structures by componentizing their donor-acceptor interactions. By unifying computational conditions, we compared the independent interactions of DNA duplexes, triplexes, and quadruplexes, which led us to evaluate a stability trend among Watson–Crick and Hoogsteen base pairing, stacking, and even ion binding. For a realistic solution-like environment, the influence of water molecules was carefully considered, and the potassium-ion preference of G-quadruplex was first analyzed at an ab initio level by considering both base-base and ion-water interactions. We devised new structure factors including hydrogen bond length, glycosidic vector angle, and twist angle, which were highly effective for comparison between computationally-predicted and experimentally-determined structures; we clarified the function of phosphate backbone during nucleobase ordering. The simulated tendency of net interaction energies agreed well with that of real world, and this agreement validates the potential of ab initio study to guide programming of complicated DNA constructs.

A Biological Take on Halogen Bonding and Other Non-Classical Non-Covalent Interactions

Classical hydrogen bonds have, for many decades, been the dominant non-covalent interaction in the toolbox that chemists and chemical engineers have used to design and control the structures of compounds and molecular assemblies as novel materials. Recently, a set of non-classical non-covalent (NC-NC) interactions have emerged that exploit the properties of the Group IV, V, VI, and VII elements of the periodic table (the tetrel, pnictogen, chalcogen, and halogen bonds, respectively). Our research group has been characterizing the prevalence, geometric constraints, and structure-function relationship specifically of the halogen bond in biological systems. We have been particularly interested in exploiting the biological halogen bonds (or BXBs) to control the structures, stabilities, and activities of biomolecules, including the DNA Holliday junction and enzymes. In this review, we first provide a set of criteria for how to determine whether BXBs or any other NC-NC interactions would have biological relevance. We then navigate the trail of studies that had led us from an initial, very biological question to our current point in the journey to establish BXBs as a tool for biomolecular engineering. Finally, we close with a perspective on future directions for this line of research.

The structural plasticity of nucleic acid duplexes revealed by WAXS and MD

Double-stranded DNA (dsDNA) and RNA (dsRNA) helices display an unusual structural diversity. Some structural variations are linked to sequence and may serve as signaling units for protein-binding partners. Therefore, elucidating the mechanisms and factors that modulate these variations is of fundamental importance. While the structural diversity of dsDNA has been extensively studied, similar studies have not been performed for dsRNA. Because of the increasing awareness of RNA's diverse biological roles, such studies are timely and increasingly important. We integrate solution x-ray scattering at wide angles (WAXS) with all-atom molecular dynamics simulations to explore the conformational ensemble of duplex topologies for different sequences and salt conditions. These tightly coordinated studies identify robust correlations between features in the WAXS profiles and duplex geometry and enable atomic-level insights into the structural diversity of DNA and RNA duplexes. Notably, dsRNA displays a marked sensitivity to the valence and identity of its associated cations.

Beyond the double helix: DNA structural diversity and the PDB

The determination of the double helical structure of DNA in 1953 remains the landmark event in the development of modern biological and biomedical science. This structure has also been the starting point for the determination of some 2000 DNA crystal structures in the subsequent 68 years. Their structural diversity has extended to the demonstration of sequence-dependent local structure in duplex DNA, to DNA bending in short and long sequences and in the DNA wound round the nucleosome, and to left-handed duplex DNAs. Beyond the double helix itself, in circumstances where DNA sequences are or can be induced to unwind from being duplex, a wide variety of topologies and forms can exist. Quadruplex structures, based on four-stranded cores of stacked G-quartets, are prevalent though not randomly distributed in the human and other genomes and can play roles in transcription, translation, and replication. Yet more complex folds can result in DNAs with extended tertiary structures and enzymatic/catalytic activity. The Protein Data Bank is the depository of all these structures, and the resource where structures can be critically examined and validated, as well as compared one with another to facilitate analysis of conformational and base morphology features. This review will briefly survey the major structural classes of DNAs and illustrate their significance, together with some examples of how the use of the Protein Data Bank by for example, data mining, has illuminated DNA structural concepts.

In silico profiling and structural insights of zinc metal ion on O6-methylguanine methyl transferase and its interactions using molecular dynamics approach

O6-methylguanine DNA methyl transferase (MGMT) is a metalloenzyme participating in the repair of alkylated DNA. In this research, we performed a comparative study for evaluating the impact of zinc metal ion on the behavior and interactions of MGMT in the both enzymatic forms of apo MGMT and holo MGMT. DNA and proliferating cell nuclear antigen (PCNA), as partners of MGMT, were utilized to evaluate molecular interactions by virtual microscopy of molecular dynamics simulation. The stability and conformational alterations of each forms (apo and holo) MGMT-PCNA, and (apo and holo) MGMT-DNA complexes were calculated by MM/PBSA method. A total of seven systems including apo MGMT, holo MGMT, free PCNA, apo MGMT-PCNA, holo MGMT-PCNA, apo MGMT-DNA, and holo MGMT-DNA complexes were simulated. In this study, we found that holo MGMT was more stable and had better folding and functional properties than that of apo MGMT. Simulation analysis of (apo and holo) MGMT-PCNA complexes displayed that the sequences of the amino acids involved in the interactions were different in the two forms of MGMT. The important amino acids of holo MGMT involved in its interaction with PCNA included E92, K101, A119, G122, N123, P124, and K125, whereas the important amino acids of apo MGMT included R128, R135, S152, N157, Y158, and L162. Virtual microscopy of molecular dynamics simulation showed that the R128 and its surrounding residues were important amino acids involved in the interaction of holo MGMT with DNA that was exactly consistent with X-ray crystallography structure. In the apo form of the protein, the N157 and its surrounding residues were important amino acids involved in the interaction with DNA. The binding free energies of - 387.976, - 396.226, - 622.227, and - 617.333 kcal/mol were obtained for holo MGMT-PCNA, apo MGMT-PCNA, holo MGMT-DNA, and apo MGMT-DNA complexes, respectively. The principle result of this research was that the area of molecular interactions differed between the two states of MGMT. Therefore, in investigations of metalloproteins, the metal ion must be preserved in their structures. Finally, it is recommended to use the holo form of metalloproteins in in vitro and in silico researches.

Origin of heat capacity increment in DNA folding: The hydration effect

BACKGROUND

Understanding DNA folding thermodynamics is crucial for prediction of DNA thermal stability. It is now well established that DNA folding is accompanied by a decrease of the heat capacity ∆c_{p, F}, however its molecular origin is not understood. In analogy to protein folding it has been assumed that this is due to dehydration of DNA constituents, however no evidence exists to support this conclusion.

METHODS

Here we analyze partial molar heat capacity of nucleic bases and nucleosides in aqueous solutions obtained from calorimetric experiments and calculate the hydration heat capacity contribution ∆c_p^hyd.

RESULTS

We present hydration heat capacity contributions of DNA constituents and show that they correlate with the solvent accessible surface area. The average contribution for nucleic base dehydration is +0.56 J mol^-1 K^-1 Å^-2 and can be used to estimate the ∆c_{p, F} contribution for DNA folding.

CONCLUSIONS

We show that dehydration is one of the major sources contributing to the observed ∆c_{p, F} increment in DNA folding. Other possible sources contributing to the overall ∆c_{p, F} should be significant but appear to compensate each other to high degree. The calculated ∆c_p^hyd for duplexes and noncanonical DNA structures agree excellently with the overall experimental ∆c_{p, F} values. By contrast, empirical parametrizations developed for proteins result in poor ∆c_p^hyd predictions and should not be applied to DNA folding.

GENERAL SIGNIFICANCE

Heat capacity is one of the main thermodynamic quantities that strongly affects thermal stability of macromolecules. At the molecular level the heat capacity in DNA folding stems from removal of water from nucleobases.

Synthesis, crystal structure and docking studies of tetracyclic 10-iodo-1,2-dihydroisoquinolino[2,1-b][1,2,4]benzothiadiazine 12,12-dioxide and its precursors

The title compound, 10-iodo-1,2-dihydroisoquinolino[2,1-b][1,2,4]benzothiadiazine 12,12-dioxide, C₁₅H₁₁IN₂O₂S (8), was synthesized via the metal-free intramolecular N-iodosuccinimide (NIS)-mediated radical oxidative sp³-C-H aminative cyclization of 2-(2'-aminobenzenesulfonyl)-1,3,4-trihydroisoquinoline, C₁₅H₁₆N₂O₂S (7). The amino adduct 7 was prepared via a two-step reaction, starting from the condensation of 2-nitrobenzenesulfonyl chloride (4) with 1,2,3,4-tetrahydroisoquinoline (5), to afford 2-(2'-nitrobenzenesulfonyl)-1,3,4-trihydroisoquinoline, C₁₅H₁₄N₂O₄S (6), in 82% yield. The catalytic hydrogenation of 6 with hydrogen gas, in the presence of 10% palladium-on-charcoal catalyst, furnished 7. Products 6-8 were characterized by their melting points, IR and NMR (¹H and ¹³C) spectroscopy, and single-crystal X-ray diffraction. The three compounds crystallized in the monoclinic space group, with 7 exhibiting classical intramolecular hydrogen bonds of 2.16 and 2.26 Å. All three crystal structures exhibit centrosymmetric pairs of intermolecular C-H...π(ring) and/or π-π stacking interactions. The docking studies of molecules 6, 7 and 8 with deoxyribonucleic acid (PDB id: 1ZEW) revealed minor-groove binding behaviours without intercalation, with 7 presenting the most favourable global energy of the three molecules. Nonetheless, molecule 8 interacted strongly with the DNA macromolecule, with an attractive van der Waals energy of -15.53 kcal mol^-1.

Structural adaptation of vertebrate endonuclease G for 5-hydroxymethylcytosine recognition and function

Modified DNA bases functionally distinguish the taxonomic forms of life—5-methylcytosine separates prokaryotes from eukaryotes and 5-hydroxymethylcytosine (^5hmC) invertebrates from vertebrates. We demonstrate here that mouse endonuclease G (mEndoG) shows specificity for both ^5hmC and Holliday junctions. The enzyme has higher affinity (>50-fold) for junctions over duplex DNAs. A ^5hmC-modification shifts the position of the cut site and increases the rate of DNA cleavage in modified versus unmodified junctions. The crystal structure of mEndoG shows that a cysteine (Cys69) is positioned to recognize ^5hmC through a thiol-hydroxyl hydrogen bond. Although this Cys is conserved from worms to mammals, a two amino acid deletion in the vertebrate relative to the invertebrate sequence unwinds an α-helix, placing the thiol of Cys69 into the mEndoG active site. Mutations of Cys69 with alanine or serine show ^5hmC-specificity that mirrors the hydrogen bonding potential of the side chain (C–H < S–H < O–H). A second orthogonal DNA binding site identified in the mEndoG structure accommodates a second arm of a junction. Thus, the specificity of mEndoG for ^5hmC and junctions derives from structural adaptations that distinguish the vertebrate from the invertebrate enzyme, thereby thereby supporting a role for ^5hmC in recombination processes.

Carboxylate-functionalized foldamer inhibitors of HIV-1 integrase and Topoisomerase 1: artificial analogues of DNA mimic proteins

Inspired by DNA mimic proteins, we have introduced aromatic foldamers bearing phosphonate groups as synthetic mimics of the charge surface of B-DNA and competitive inhibitors of some therapeutically relevant DNA-binding enzymes: the human DNA Topoisomerase 1 (Top1) and the human HIV-1 integrase (HIV-1 IN). We now report on variants of these anionic foldamers bearing carboxylates instead of phosphonates. Several new monomers have been synthesized with protecting groups suitable for solid phase synthesis (SPS). Six hexadecaamides have been prepared using SPS. Proof of their resemblance to B-DNA was brought by the first crystal structure of one of these DNA-mimic foldamers in its polyanionic form. While some of the foldamers were found to be as active as, or even more active than, the original phosphonate oligomers, others had no activity at all or could even stimulate enzyme activity in vitro. Some foldamers were found to have differential inhibitory effects on the two enzymes. These results demonstrate a strong dependence of inhibitory activity on foldamer structure and charge distribution. They open broad avenues for the development of new classes of derivatives that could inhibit the interaction of specific proteins with their DNA target thereby influencing the cellular pathways in which they are involved.

Changes in Microenvironment Modulate the B- to A-DNA Transition

B- to A-DNA transition is known to be sensitive to the macroscopic properties of the solution, such as salt and ethanol concentrations. Microenvironmental effects on DNA conformational transition have been broadly studied. Providing an intuitive picture of how DNA responds to environmental changes is, however, still needed. Analyzing the chemical equilibrium of B-to-A DNA transition at critical concentrations, employing explicit-solvent simulations, is envisioned to help understand such microenvironmental effects. In the present study, free-energy calculations characterizing the B- to A-DNA transition and the distribution of cations were carried out in solvents with different ethanol concentrations. With the addition of ethanol, the most stable structure of DNA changes from the B- to A-form, in agreement with previous experimental observation. In 60% ethanol, a chemical equilibrium is found, showing reversible transition between B- and A-DNA. Analysis of the microenvironment around DNA suggests that with the increase of ethanol concentration, the cations exhibit a significant tendency to move toward the backbone, and mobility of water molecules around the major groove and backbone decreases gradually, leading eventually to a B-to-A transition. The present results provide a free-energy view of DNA microenvironment and of the role of cation motion in the conformational transition.

Holliday Junctions Formed from Human Telomeric DNA

Cells have evolved inherent mechanisms, like homologous recombination (HR), to repair damaged DNA. However, repairs at telomeres can lead to genomic instability, often associated with cancer. While most rapidly dividing cells employ telomerase, the others maintain telomere length through HR-dependent alternative lengthening of telomeres (ALT) pathways. Here we describe the crystal structures of Holliday junction intermediates of the HR-dependent ALT mechanism. Using an extended human telomeric repeat, we also report the crystal structure of two Holliday junctions in close proximity, which associate together through strand exchange to form a hemicatenated double Holliday junction. Our combined structural results demonstrate that ACC nucleotides in the C-rich lagging strand (5'-CTAACCCTAA-3') at the telomere repeat sequence constitute a conserved structural feature that constrains crossover geometry and is a preferred site for Holliday junction formation in telomeres.

Free-Energy Profiles for A-/B-DNA Conformational Transitions in Isolated and Aggregated States from All-Atom Molecular Dynamics Simulation

In ordinary aqueous solution, B-DNA is the major structural form of DNA. After the addition of ethanol, DNA is thought to be aggregated/condensed in the A-form structure. However, there is uncertainty as to whether the B-to-A conformational change is connected to the aggregation/condensation steps. In this study, we performed all-atom molecular dynamics simulations and calculated the free-energy surface involved in the A/B conformational transition for isolated and aggregated Dickerson-Drew dodecamers (DDDs) in water and 85% ethanol environments. We found in the case of an isolated DDD, the overall free-energy profile is entirely downhill to give the B-DNA conformation in both water and 85% ethanol. However, in the aggregated state and 85% ethanol environment, there is a free-energy minimum associated with the A-DNA region in addition to the global B-DNA minimum, and there is a ∼3 kcal/mol free-energy barrier to the A-to-B conformational change. The molecular dynamics results suggest that aggregation of DNA is essential for forming A-DNA.

Single-Stranded DNA Oligonucleotides Retain Rise Coordinates Characteristic of Double Helices

The structures of single-stranded DNA oligonucleotides from dimeric to hexameric sequences have been thoroughly investigated. Computations performed at the density functional level of theory including dispersion forces and solvation show that single-stranded helices adopt conformations very close to crystallographic B-DNA, with rise coordinates amounting up to 3.3 Å. Previous results, suggesting that single strands should be shorter than double helices, largely originated from the incompleteness of the adopted basis set. Although sensible deviations with respect to standard B-DNA are predicted, computations indicate that sequences rich in stacked adenines are the most ordered ones, favoring the B-DNA pattern and inducing regular arrangements also on flanking nucleobases. Several structural properties of double helices rich in adenine are indeed already reflected by the corresponding single strands.

AMOEBA Polarizable Atomic Multipole Force Field for Nucleic Acids

The AMOEBA polarizable atomic multipole force field for nucleic acids is presented. Valence and electrostatic parameters were determined from high-level quantum mechanical data, including structures, conformational energy, and electrostatic potentials, of nucleotide model compounds. Previously derived parameters for the phosphate group and nucleobases were incorporated. A total of over 35 μs of condensed-phase molecular dynamics simulations of DNA and RNA molecules in aqueous solution and crystal lattice were performed to validate and refine the force field. The solution and/or crystal structures of DNA B-form duplexes, RNA duplexes, and hairpins were captured with an average root-mean-squared deviation from NMR structures below or around 2.0 Å. Structural details, such as base pairing and stacking, sugar puckering, backbone and χ-torsion angles, groove geometries, and crystal packing interfaces, agreed well with NMR and/or X-ray. The interconversion between A- and B-form DNAs was observed in ethanol-water mixtures at 328 K. Crystal lattices of B- and Z-form DNA and A-form RNA were examined with simulations. For the RNA tetraloop, single strand tetramers, and HIV TAR with 29 residues, the simulated conformational states, 3 J-coupling, nuclear Overhauser effect, and residual dipolar coupling data were compared with NMR results. Starting from a totally unstacked/unfolding state, the rCAAU tetranucleotide was folded into A-form-like structures during ∼1 μs molecular dynamics simulations.

A DNA structural alphabet provides new insight into DNA flexibility

DNA is a structurally plastic molecule, and its biological function is enabled by adaptation to its binding partners. To identify the DNA structural polymorphisms that are possible in such adaptations, the dinucleotide structures of 60 000 DNA steps from sequentially nonredundant crystal structures were classified and an automated protocol assigning 44 distinct structural (conformational) classes called NtC (for Nucleotide Conformers) was developed. To further facilitate understanding of the DNA structure, the NtC were assembled into the DNA structural alphabet CANA (Conformational Alphabet of Nucleic Acids) and the projection of CANA onto the graphical representation of the molecular structure was proposed. The NtC classification was used to define a validation score called confal, which quantifies the conformity between an analyzed structure and the geometries of NtC. NtC and CANA assignment were applied to analyze the structural properties of typical DNA structures such as Dickerson-Drew dodecamers, guanine quadruplexes and structural models based on fibre diffraction. NtC, CANA and confal assignment, which is accessible at the website https://dnatco.org, allows the quantitative assessment and validation of DNA structures and their subsequent analysis by means of pseudo-sequence alignment. An animated Interactive 3D Complement (I3DC) is available in Proteopedia at http://proteopedia.org/w/Journal:Acta_Cryst_D:2.

Sulfur as an Acceptor to Bromine in Biomolecular Halogen Bonds

The halogen bond (X-bond) has become an important design element in chemistry, including medicinal chemistry and biomolecular engineering. Although oxygen is the most prevalent and best characterized X-bond acceptor in biomolecules, the interaction is seen with nitrogen, sulfur, and aromatic systems as well. In this study, we characterize the structure and thermodynamics of a Br···S X-bond between a 5-bromouracil base and a phosphorothioate in a model DNA junction. The single-crystal structure of the junction shows the geometry of the Br···S to be variable, while calorimetric studies show that the anionic S acceptor is comparable to or slightly more stable than the analogous O acceptor, with a -3.5 kcal/mol difference in ΔΔH^25°C and -0.4 kcal/mol ΔΔG^25°C (including an entropic penalty ΔΔS^25°C of -10 cal/(mol K)). Thus sulfur is shown to be a favorable acceptor for bromine X-bonds, extending the application of this interaction for the design of inhibitors and biological materials.

Structure of the Holliday junction: applications beyond recombination

The Holliday junction (HJ) is an essential element in recombination and related mechanisms. The structure of this four-stranded DNA assembly, which is now well-defined alone and in complex with proteins, has led to its applications in areas well outside of molecular recombination, including nanotechnology and biophysics. This minireview explores some interesting recent research on the HJ, as it has been adapted to design regular two- or three-dimensional lattices for crystal engineering, and more complex systems through DNA origami. In addition, the sequence dependence of the structure is discussed in terms how it can be applied to characterize the geometries and energies of various noncovalent interactions, including halogen bonds in oxidatively damaged (halogenated) bases and hydrogen bonds associated with the epigenetic 5-hydroxylmethylcytosine base.

Polarizable Force Field for DNA Based on the Classical Drude Oscillator: II. Microsecond Molecular Dynamics Simulations of Duplex DNA

The structure and dynamics of DNA are governed by a sensitive balance between base stacking and pairing, hydration, and interactions with ions. Force-field models that include explicit representations of electronic polarization are capable of more accurately modeling the subtle details of these interactions versus commonly used additive force fields. In this work, we validate our recently refined polarizable force field for DNA based on the classical Drude oscillator model, in which electronic degrees of freedom are represented as negatively charged particles attached to their parent atoms via harmonic springs. The previous version of the force field, called Drude-2013, produced stable A- and B-DNA trajectories on the order of hundreds of nanoseconds, but deficiencies were identified that included weak base stacking ultimately leading to distortion of B-DNA duplexes and unstable Z-DNA. As a result of extensive refinement of base nonbonded terms and bonded parameters in the deoxyribofuranose sugar and phosphodiester backbone, we demonstrate that the new version of the Drude DNA force field is capable of simulating A- and B-forms of DNA on the microsecond time scale and the resulting conformational ensembles agree well with a broad set of experimental properties, including solution X-ray scattering profiles. In addition, simulations of Z-form duplex DNA in its crystal environment are stable on the order of 100 ns. The revised force field is to be called Drude-2017.

Relationships between hydrogen bonds and halogen bonds in biological systems

The recent recognition that halogen bonding (XB) plays important roles in the recognition and assembly of biological molecules has led to new approaches in medicinal chemistry and biomolecular engineering. When designing XBs into strategies for rational drug design or into a biomolecule to affect its structure and function, we must consider the relationship between this interaction and the more ubiquitous hydrogen bond (HB). In this review, we explore these relationships by asking whether and how XBs can replace, compete against or behave independently of HBs in various biological systems. The complex relationships between the two interactions inform us of the challenges we face in fully utilizing XBs to control the affinity and recognition of inhibitors against their therapeutic targets, and to control the structure and function of proteins, nucleic acids and other biomolecular scaffolds.

Protein/DNA interactions in complex DNA topologies: expect the unexpected

DNA supercoiling results in compacted DNA structures that can bring distal sites into close proximity. It also changes the local structure of the DNA, which can in turn influence the way it is recognised by drugs, other nucleic acids and proteins. Here, we discuss how DNA supercoiling and the formation of complex DNA topologies can affect the thermodynamics of DNA recognition. We then speculate on the implications for transcriptional control and the three-dimensional organisation of the genetic material, using examples from our own simulations and from the literature. We introduce and discuss the concept of coupling between the multiple length-scales associated with hierarchical nuclear structural organisation through DNA supercoiling and topology.

Effect of Hydroxymethylcytosine on the Structure and Stability of Holliday Junctions

5-Hydroxymethylcytosine (5hmC) is an epigenetic marker that has recently been shown to promote homologous recombination (HR). In this study, we determine the effects of 5hmC on the structure, thermodynamics, and conformational dynamics of the Holliday junction (the four-stranded DNA intermediate associated with HR) in its native stacked-X form. The hydroxymethyl and the control methyl substituents are placed in the context of an amphimorphic GxCC trinucleotide core sequence (where xC is C, 5hmC, or the methylated 5mC), which is part of a sequence also recognized by endonuclease G to promote HR. The hydroxymethyl group of the 5hmC junction adopts two distinct rotational conformations, with an in-base-plane form being dominant over the competing out-of-plane rotamer that has typically been seen in duplex structures. The in-plane rotamer is seen to be stabilized by a more stable intramolecular hydrogen bond to the junction backbone. Stabilizing hydrogen bonds (H-bonds) formed by the hydroxyl substituent in 5hmC or from a bridging water in the 5mC structure provide approximately 1.5-2 kcal/mol per interaction of stability to the junction, which is mostly offset by entropy compensation, thereby leaving the overall stability of the G5hmCC and G5mCC constructs similar to that of the GCC core. Thus, both methyl and hydroxymethyl modifications are accommodated without disrupting the structure or stability of the Holliday junction. Both 5hmC and 5mC are shown to open the structure to make the junction core more accessible. The overall consequences of incorporating 5hmC into a DNA junction are thus discussed in the context of the specificity in protein recognition of the hydroxymethyl substituent through direct and indirect readout mechanisms.

Protein/DNA interactions in complex DNA topologies: expect the unexpected

DNA supercoiling results in compacted DNA structures that can bring distal sites into close proximity. It also changes the local structure of the DNA, which can in turn influence the way it is recognised by drugs, other nucleic acids and proteins. Here, we discuss how DNA supercoiling and the formation of complex DNA topologies can affect the thermodynamics of DNA recognition. We then speculate on the implications for transcriptional control and the three-dimensional organisation of the genetic material, using examples from our own simulations and from the literature. We introduce and discuss the concept of coupling between the multiple length-scales associated with hierarchical nuclear structural organisation through DNA supercoiling and topology.

Crystal structure of a complex between β-glucopyranose and a macrocyclic receptor with dendritic multicharged water solubilizing chains

Using commercial screens for crystallization of biomolecules and taking advantage of the use of racemic crystallography allowed the production of X-ray quality single crystals and the elucidation at 1.08 Å resolution of the solid state structure of a difficult target: the complex between glucopyranose and a water soluble macrocyclic receptor equipped with dendritic multianionic solubilizing chains.