A perspective on the molecular simulation of DNA from structural and functional aspects

Spermine and spermidine SI-PPCs: Molecular dynamics reveals enhanced biomolecular interactions

In this paper the effects on the interaction of highly positively charged substitution-inert platinum polynuclear complexes (SI-PPCs) with negatively charged DNA and heparin are examined and compared by theoretical chemistry methods. Electrostatic and hydrogen bonding interactions contribute to the overall effects on the biomolecule. Root Mean Square (RMS) deviation, Solvent Accessible Surface, RMS fluctuation, and interaction analysis all confirm similar effects on both biomolecules, dictated predominantly by the total positive charge and total number of hydrogen bonds formed. Especially, changes in structural parameters suggesting condensation and reduction of available surface area will reduce or prevent normal protein recognition and may thus potentially inhibit biological mechanisms related to apoptosis (DNA) or reduced vascularization viability (HEP). Thermodynamic analyses supported these findings with favourable interaction energies. The comparison of DNA and heparin confirms the general intersectionality between the two biomolecules and confirms the intrinsic dual-nature function of this chemotype. The distinction between the two-limiting mode of actions (HS or DNA-centred) could reflect an intriguing balance between extracellular (GAG) and intracellular (DNA) binding and affinities. The results underline the need to fully understand GAG-small molecule interactions and their contribution to drug pharmacology and related therapeutic modalities. This report contributes to that understanding.

Generating High-Precision Force Fields for Molecular Dynamics Simulations to Study Chemical Reaction Mechanisms Using Molecular Configuration Transformer

Theoretical studies on chemical reaction mechanisms have been crucial in organic chemistry. Traditionally, calculating the manually constructed molecular conformations of transition states for chemical reactions using quantum chemical calculations is the most commonly used method. However, this way is heavily dependent on individual experience and chemical intuition. In our previous study, we proposed a research paradigm that used enhanced sampling in molecular dynamics simulations to study chemical reactions. This approach can directly simulate the entire process of a chemical reaction. However, the computational speed limited the use of high-precision potential energy functions for simulations. To address this issue, we presented a scheme for training high-precision force fields for molecular modeling using a previously developed graph-neural-network-based molecular model, molecular configuration transformer. This potential energy function allowed for highly accurate simulations at a low computational cost, leading to more precise calculations of the mechanism of chemical reactions. We applied this approach to study a Claisen rearrangement reaction and a carbonyl insertion reaction catalyzed by manganese.

Structural and dynamical aspect of DNA motif sequence specific binding of AP-1 transcription factor

Activator protein-1 (AP-1) comprises one of the largest and most evolutionary conserved families of ubiquitous eukaryotic transcription factors that act as a pioneer factor. Diversity in DNA binding interaction of AP-1 through a conserved basic-zipper (bZIP) domain directs in-depth understanding of how AP-1 achieves its DNA binding selectivity and consequently gene regulation specificity. Here, we address the structural and dynamical aspects of the DNA target recognition process of AP-1 using microsecond-long atomistic simulations based on the structure of the human AP-1 FosB/JunD bZIP-DNA complex. Our results show the unique role of DNA shape features in selective base specific interactions, characteristic ion population, and solvation properties of DNA grooves to form the motif sequence specific AP-1-DNA complex. The TpG step at the two terminals of the AP-1 site plays an important role in the structural adjustment of DNA by modifying the helical twist in the AP-1 bound state. We addressed the role of intrinsic motion of the bZIP domain in terms of opening and closing gripper motions of DNA binding helices, in target site recognition and binding of AP-1 factors. Our observations suggest that binding to the cognate motif in DNA is mainly accompanied with the precise adjustment of closing gripper motion of DNA binding helices of the bZIP domain.

Multi-scale gene regulation mechanism: Spatiotemporal transmission of genetic information

Gene expression is regulated by many factors, including transcription factors, chromatin three-dimensional topology, modifications of DNA and histone proteins, and non-coding RNAs. The execution of these complex mechanisms requires an effectively coordinated regulation system. In this review, we emphasize that the multi-scale heterogeneous DNA sequence plays a fundamental and important role for gene expression activity and usage of different means of epigenetic regulation. We illustrate here that the chromatin structure organization provides a stage for spatiotemporal regulation between different genes or gene modules and to realize their downstream functional cooperation. Such a perspective expands our understanding of the central dogma: In addition to one-dimensional sequence information, inter-gene interactions can also be transferred from DNA and RNA to protein levels.

Markov state models elucidate the stability of DNA influenced by the chiral 5S-Tg base

The static and dynamic structures of DNA duplexes affected by 5S-Tg (Tg, Thymine glycol) epimers were studied using MD simulations and Markov State Models (MSMs) analysis. The results show that the 5S,6S-Tg base caused little perturbation to the helix, and the base-flipping barrier was determined to be 4.4 kcal mol-1 through the use of enhanced sampling meta-eABF calculations, comparable to 5.4 kcal mol-1 of the corresponding thymine flipping. Two conformations with the different hydrogen bond structures between 5S,6R-Tg and A19 were identified in several independent MD trajectories. The 5S,6R-Tg:O6HO6•••N1:A19 hydrogen bond is present in the high-energy conformation displaying a clear helical distortion, and near barrier-free Tg base flipping. The low-energy conformation always maintains Watson-Crick base pairing between 5S,6R-Tg and A19, and 5S-Tg base flipping is accompanied by a small barrier of ca. 2.0 KBT (T = 298 K). The same conformations are observed in the MSMs analysis. Moreover, the transition path and metastable structures of the damaged base flipping are for the first time verified through MSMs analysis. The data clearly show that the epimers have completely different influence on the stability of the DNA duplex, thus implying different enzymatic mechanisms for DNA repair.

Staggered intercalation of DNA duplexes with base-pair modulation by two distinct drug molecules induces asymmetric backbone twisting and structure polymorphism

The use of multiple drugs simultaneously targeting DNA is a promising strategy in cancer therapy for potentially overcoming single drug resistance. In support of this concept, we report that a combination of actinomycin D (ActD) and echinomycin (Echi), can interact in novel ways with native and mismatched DNA sequences, distinct from the structural effects produced by either drug alone. Changes in the former with GpC and CpG steps separated by a A:G or G:A mismatch or in a native DNA with canonical G:C and C:G base pairs, result in significant asymmetric backbone twists through staggered intercalation and base pair modulations. A wobble or Watson–Crick base pair at the two drug-binding interfaces can result in a single-stranded ‘chair-shaped’ DNA duplex with a straight helical axis. However, a novel sugar-edged hydrogen bonding geometry in the G:A mismatch leads to a ‘curved-shaped’ duplex. Two non-canonical G:C Hoogsteen base pairings produce a sharply kinked duplex in different forms and a four-way junction-like superstructure, respectively. Therefore, single base pair modulations on the two drug-binding interfaces could significantly affect global DNA structure. These structures thus provide a rationale for atypical DNA recognition via multiple DNA intercalators and a structural basis for the drugs’ potential synergetic use.

Microscopic Insight into pH-Dependent Conformational Dynamics and Noncanonical Base Pairing in Telomeric i-Motif DNA

Gene regulatory functions of noncanonical i-motif DNA are associated with dynamic i-motif formation in the cellular environment and pH variation. With atomistic simulations, we show the dramatic influence of solvent pH on the conformational dynamics of biologically relevant telomeric i-motif DNA coupled with protonation of cytosine bases in different conformations. We rationalized the pH-dependent dynamics and conformational variability of the i-motif in terms of base pairing and specific loop motions. The human telomeric i-motif is found to acquire various metastable folded conformations at pH values near the pKa of cytosine with the formation of a noncanonical C:C W:W trans base pair along with the hemiprotonated C:C+ pairs in the i-motif core. pH-dependent dynamics and the local solvent structure of i-motif DNA imply that the presence of a cosolvent or molecular crowding can promote i-motif formation in vivo by changing the conformational fluctuations and hydration state of the structure.

Dynamics of 5R-Tg Base Flipping in DNA Duplexes Based on Simulations─Agreement with Experiments and Beyond

Damaged or mismatched DNA bases are normally thought to be able to flip out of the helical stack, providing enzymes with access to the faulty genetic information otherwise hidden inside the helix. Thymine glycol (Tg) is one of the most common products of nucleic acid damage. However, the static and dynamic structures of DNA duplexes affected by 5R-Tg epimers are still not clearly understood, including the ability of these to undergo spontaneous base flipping. Structural effects of the 5R-Tg epimers on the duplex DNA are herein studied using molecular dynamics together with reliable DFT based calculations. In comparison with the corresponding intact DNA, the cis-5R,6S-Tg epimer base causes little perturbation to the duplex DNA, and a barrier of 4.9 kcal mol^–1 is obtained by meta-eABF for cis-5R,6S-Tg base flipping out of the duplex DNA, comparable to the 5.4 kcal mol^–1 obtained for the corresponding thymine flipping in intact DNA. For the trans-5R,6R-Tg epimer, three stable local structures were identified, of which the most stable disrupts the Watson–Crick hydrogen-bonded G5/C20 base pair, leading to conformational distortion of the duplex. Interestingly, the relative barrier height of the 5R-Tg flipping is only 1.0 kcal mol^–1 for one of these trans-5R,6R-Tg epimers. Water bridge interactions were identified to be essential for 5R-Tg flipping. The study clearly demonstrates the occurrence of partial trans-5R,6R-Tg epimer flipping in solution.