Assessing the Current State of Amber Force Field Modifications for DNA

Helicobacter pylori East Asian type CagA hijacks more SHIP2 by its EPIYA-D motif to potentiate the oncogenicity

CagA is a significant oncogenic factor injected into host cells by Helicobacter pylori, which is divided into two subtypes: East Asian type (CagAE), characterized by the EPIYA-D motif, and western type (CagAW), harboring the EPIYA-C motif. CagAE has been reported to have higher carcinogenicity than CagAW, although the underlying reason is not fully understood. SHIP2 is an intracellular phosphatase that can be recruited by CagA to perturb the homeostasis of intracellular signaling pathways. In this study, we found that SHIP2 contributes to the higher oncogenicity of CagAE. Co-Immunoprecipitation and Pull-down assays showed that CagAE bind more SHIP2 than CagAW. Immunofluorescence staining showed that a higher amount of SHIP2 recruited by CagAE to the plasma membrane catalyzes the conversion of PI(3,4,5)P3 into PI(3,4)P2. This alteration causes higher activation of Akt signaling, which results in enhanced IL-8 secretion, migration, and invasion of the infected cells. SPR analysis showed that this stronger interaction between CagAE and SHIP2 stems from the higher affinity between the EPIYA-D motif of CagAE and the SH2 domain of SHIP2. Structural analysis revealed the crucial role of the Phe residue at the Y + 5 position in EPIYA-D. After mutating Phe of CagAE into Asp (the corresponding residue in the EPIYA-C motif) or Ala, the activation of downstream Akt signaling was reduced and the malignant transformation of infected cells was alleviated. These findings revealed that CagAE hijacks SHIP2 through its EPIYA-D motif to enhance its carcinogenicity, which provides a better understanding of the higher oncogenic risk of H. pylori CagAE.

Molecular architecture and functional dynamics of the pre-incision complex in nucleotide excision repair

Nucleotide excision repair (NER) is vital for genome integrity. Yet, our understanding of the complex NER protein machinery remains incomplete. Combining cryo-EM and XL-MS data with AlphaFold2 predictions, we build an integrative model of the NER pre-incision complex(PInC). Here TFIIH serves as a molecular ruler, defining the DNA bubble size and precisely positioning the XPG and XPF nucleases for incision. Using simulations and graph theoretical analyses, we unveil PInC's assembly, global motions, and partitioning into dynamic communities. Remarkably, XPG caps XPD's DNA-binding groove and bridges both junctions of the DNA bubble, suggesting a novel coordination mechanism of PInC's dual incision. XPA rigging interlaces XPF/ERCC1 with RPA, XPD, XPB, and 5' ssDNA, exposing XPA's crucial role in licensing the XPF/ERCC1 incision. Mapping disease mutations onto our models reveals clustering into distinct mechanistic classes, elucidating xeroderma pigmentosum and Cockayne syndrome disease etiology.

One-pot assay using a target-driven split aptamer recognition and assembly strategy for convenient and rapid detection of gliotoxin

An aptamer targeting gliotoxin (GTX) was optimized to increase the binding affinity by approximately 20 times and achieve higher structural stability and targeting specificity. Molecular dynamics simulations were used to explore the molecular mechanism and key action sites underlying the recognition of GTX by the optimized aptamer. Subsequently, the optimized aptamer was split into two fragments and a convenient and rapid one-pot assay for GTX detection was successfully established using a target-driven split aptamer recognition and assembly strategy. The method exhibited a good linear range of 0.128 nM to 2 μM, a low detection limit of 0.07 nM, and excellent selectivity for GTX. Furthermore, the method had good accuracy and stability in real sample analysis. Therefore, the developed one-pot method provides a reliable, convenient, and cost-effective approach for the widespread application of GTX detection.

Mechanism of Dual-Site Recognition in a Classic DNA Aptamer

Nucleic acid aptamers possess unique advantages in specific recognition. However, the lack of in-depth investigation into their dynamic recognition mechanisms has restricted their rational design and potential applications in fields such as biosensing and targeted therapy. We herein utilized enhanced sampling molecular dynamics to address affinities of adenosine monophosphate (AMP) to the dual binding sites in the DNA aptamer, focusing on the dynamic recognition mechanism and pathways. The present results indicate that in addition to the widely known intermolecular interactions, inequivalence of chemical environments of the two binding sites leads to slightly higher stability of AMP binding to the site proximal to the aptamer terminus. In the presence of two AMPs captured by the two sites, each binding free energy is enhanced. In particular, an additional hydrogen bond of AMP to A10 is introduced in the dual-site binding complex, which increases the binding energy from -4.25 ± 0.47 to -9.48 ± 0.33 kcal mol-1 in the site close to the loop. For the dual-site recognition process, the free energy landscape and minimum free energy pathway calculations elucidate the crucial role of electrostatic interactions between the AMP phosphate groups and Na+ ions in positively cooperative binding mechanisms.

The molecular mechanism for TERRA recruitment and annealing to telomeres

Telomeric repeat containing RNA (TERRA) is a noncoding RNA that is transcribed from telomeres. Previous study showed that TERRA trans anneals by invading into the telomeric duplex to form an R-loop in mammalian cells. Here, we elucidate the molecular mechanism underlying TERRA recruitment and invasion into telomeres in the context of shelterin proteins, RAD51 and RNase H using single molecule (sm) assays. We demonstrate that TERRA trans annealing into telomeric DNA exhibits dynamic movement that is stabilized by TRF2. TERRA annealing to the telomeric duplex results in the formation of a stable triplex structure which differs from a conventional R-loop. We identified that the presence of a sub-telomeric DNA and a telomeric overhang in the form of a G-quadruplex significantly enhances TERRA annealing to telomeric duplex. We also demonstrate that RAD51-TERRA complex invades telomere duplex more efficiently than TERRA alone. Additionally, TRF2 increases TERRA affinity to telomeric duplex and protects it from RNase H digestion. In contrast, TRF1 represses TERRA annealing to telomeric duplex and fails to provide protection against RNase H digestion. Our findings provide an in-depth molecular mechanism underpinning TERRA recruitment and annealing to the telomere.

STAG2-RAD21 complex: A unidirectional DNA ratchet mechanism in loop extrusion

DNA loop extrusion plays a key role in the regulation of gene expression and the structural arrangement of chromatin. Most existing mechanistic models of loop extrusion depend on some type of ratchet mechanism, which should permit the elongation of loops while preventing their collapse, by enabling DNA to move in only one direction. STAG2 is already known to exert a role as DNA anchor, but the available structural data suggest a possible role in unidirectional DNA motion. In this work, a computational simulation framework was constructed to evaluate whether STAG2 could enforce such unidirectional displacement of a DNA double helix. The results reveal that STAG2 V-shape allows DNA sliding in one direction, but blocks opposite DNA movement via a linear ratchet mechanism. Furthermore, these results suggest that RAD21 binding to STAG2 controls its flexibility by narrowing the opening of its V-shape, which otherwise remains widely open in absence of RAD21. Therefore, in the proposed model, in addition to its already described role as a DNA anchor, the STAG2-RAD21 complex would be part of a ratchet mechanism capable of exerting directional selectivity on DNA sliding during loop extrusion. The identification of the molecular basis of the ratchet mechanism of loop extrusion is a critical step in unraveling new insights into a broad spectrum of chromatin activities and their implications for the mechanisms of chromatin-related diseases.

Aptamer biosensor design for the detection of endocrine-disrupting chemicals small organic molecules using novel bioinformatics methods

Endocrine-disrupting chemicals (EDCs) are substances that can disrupt the normal functioning of hormones.Using aptamers, which are biological recognition elements, biosensors can quickly and accurately detect EDCs in environmental samples. However, the elucidation of aptamer structures by conventional methods is highly challenging due to their complexity. This has led to the development of three-dimensional aptamer structures based on different models and techniques. To do this, we developed a way to predict the 3D structures of the SS DNA needed for this sequence by starting with an aptamer sequence that has biosensor properties specific to bisphenol-A (BPA), one of the chemicals found in water samples that can interfere with hormones. In addition, we will elucidate the intermolecular mechanisms and binding affinity between aptamers and endocrine disruptors using bioinformatics techniques such as molecular docking, molecular dynamics simulation, and binding energies. The outcomes of our study are to compare modeling programs and force fields to see how reliable they are and how well they agree with results found in the existing literature, to understand the intermolecular mechanisms and affinity of aptamer-based biosensors, and to find a new way to make aptamers that takes less time and costs less.

Mechanical response of double-stranded DNA: Bend, twist, and overwind

We employed all-atom molecular dynamics simulations to explore the mechanical response of bending, twisting, and overwinding for double-stranded DNA (dsDNA). We analyzed the bending and twisting deformations, as well as their stiffnesses, using the tilt, roll, and twist modes under stretching force. Findings indicate that the roll and twist angles vary linearly with the stretching force but show opposite trends. The tilt, roll, and twist elastic moduli are considered constants, while the coupling between roll and twist modes slightly decreases under stretching force. The effect of the stretching force on the roll and twist modes, including both their deformations and elasticities, exhibits sequence-dependence, with symmetry around the base pair step. Furthermore, we examined the overwinding path and mechanism of dsDNA from the perspective of the stiffness matrix, based on the tilt, roll, and twist modes. The correlations among tilt, roll, and twist angles imply an alternative overwinding pathway via twist-roll coupling when dsDNA is stretched, wherein entropic contribution prevails.

A Matter of Charge: Electrostatically Tuned Coassembly of Amphiphilic Peptides

Coassembly of peptide biomaterials offers a compelling avenue to broaden the spectrum of hierarchically ordered supramolecular nanoscale structures that may be relevant for biomedical and biotechnological applications. In this work coassemblies of amphiphilic and oppositely charged, anionic and cationic, β-sheet peptides are studied, which may give rise to a diverse range of coassembled forms. Mixtures of the peptides show significantly lower critical coassembly concentration (CCC) values compared to those of the individual pure peptides. Intriguingly, the highest formation of coassembled fibrils is found to require excess of the cationic peptide whereas equimolar mixtures of the peptides exhibited the maximum folding into β-sheet structures. Mixtures of the peptides coassembled sequentially from solutions at concentrations surpassing each peptide's intrinsic critical assembly concentration (CAC), are also found to require a higher portion of the cationic peptide to stabilize hydrogels. This study illuminates a systematic investigation of oppositely charged β-sheet peptides over a range of concentrations, in solutions and in hydrogels. The results may be relevant to the fundamental understanding of such intricate charge-driven assembly systems and to the formulation of peptide-based nanostructures with diverse functionalities.

Molecular basis of facilitated target search and sequence discrimination of TALE homeodomain transcription factor Meis1

Transcription factors specifically bind to their consensus sequence motifs and regulate transcription efficiency. Transcription factors are also able to non-specifically contact the phosphate backbone of DNA through electrostatic interaction. The homeodomain of Meis1 TALE human transcription factor (Meis1-HD) recognizes its target DNA sequences via two DNA contact regions, the L1-α1 region and the α3 helix (specific binding mode). This study demonstrates that the non-specific binding mode of Meis1-HD is the energetically favored process during DNA binding, achieved by the interaction of the L1-α1 region with the phosphate backbone. An NMR dynamics study suggests that non-specific binding might set up an intermediate structure which can then rapidly and easily find the consensus region on a long section of genomic DNA in a facilitated binding process. Structural analysis using NMR and molecular dynamics shows that key structural distortions in the Meis1-HD-DNA complex are induced by various single nucleotide mutations in the consensus sequence, resulting in decreased DNA binding affinity. Collectively, our results elucidate the detailed molecular mechanism of how Meis1-HD recognizes single nucleotide mutations within its consensus sequence: (i) through the conformational features of the α3 helix; and (ii) by the dynamic features (rigid or flexible) of the L1 loop and the α3 helix. These findings enhance our understanding of how single nucleotide mutations in transcription factor consensus sequences lead to dysfunctional transcription and, ultimately, human disease.

Machine-learned molecular mechanics force fields from large-scale quantum chemical data

The development of reliable and extensible molecular mechanics (MM) force fields-fast, empirical models characterizing the potential energy surface of molecular systems-is indispensable for biomolecular simulation and computer-aided drug design. Here, we introduce a generalized and extensible machine-learned MM force field, espaloma-0.3, and an end-to-end differentiable framework using graph neural networks to overcome the limitations of traditional rule-based methods. Trained in a single GPU-day to fit a large and diverse quantum chemical dataset of over 1.1 M energy and force calculations, espaloma-0.3 reproduces quantum chemical energetic properties of chemical domains highly relevant to drug discovery, including small molecules, peptides, and nucleic acids. Moreover, this force field maintains the quantum chemical energy-minimized geometries of small molecules and preserves the condensed phase properties of peptides and folded proteins, self-consistently parametrizing proteins and ligands to produce stable simulations leading to highly accurate predictions of binding free energies. This methodology demonstrates significant promise as a path forward for systematically building more accurate force fields that are easily extensible to new chemical domains of interest.

Computational Modeling Study of the Molecular Basis of dNTP Selectivity in Human Terminal Deoxynucleotidyltransferase

Human terminal deoxynucleotidyl transferase (TdT) can catalyze template-independent DNA synthesis during the V(D)J recombination and DNA repair through nonhomologous end joining. The capacity for template-independent random addition of nucleotides to single-stranded DNA makes this polymerase useful in various molecular biological applications involving sequential stepwise synthesis of oligonucleotides using modified dNTP. Nonetheless, a serious limitation to the applications of this enzyme is strong selectivity of human TdT toward dNTPs in the order dGTP > dTTP ≈ dATP > dCTP. This study involved molecular dynamics to simulate a potential impact of amino acid substitutions on the enzyme's selectivity toward dNTPs. It was found that the formation of stable hydrogen bonds between a nitrogenous base and amino acid residues at positions 395 and 456 is crucial for the preferences for dNTPs. A set of single-substitution and double-substitution mutants at these positions was analyzed by molecular dynamics simulations. The data revealed two TdT mutants-containing either substitution D395N or substitutions D395N+E456N-that possess substantially equalized selectivity toward various dNTPs as compared to the wild-type enzyme. These results will enable rational design of TdT-like enzymes with equalized dNTP selectivity for biotechnological applications.

Quadruplex-duplex junction in LTR-III: A molecular insight into the complexes with BMH-21, namitecan and doxorubicin

Quadruplex-Duplex (Q-D) junctions are unique structural motifs garnering increasing interest as drug targets, due to their frequent occurrence in genomic sequences. The viral HIV LTR-III sequence was chosen as a Q-D junction model to study the affinity of the selected compounds BMH-21, namitecan (ST-1968), and doxorubicin (DOXO), all containing a planar polycyclic aromatic moiety, linked to either one short aminoalkyl or an aminoglycosyl group. A multidisciplinary approach that combines NMR spectroscopy, molecular modelling, circular dichroism (CD) and fluorescence spectroscopy was employed. The studied ligands induced moderate but clear stabilization to the Q-D junction by interacting with the interfacial tetrad. DOXO was found to be the best Q-D junction binder. Interestingly, the removal of the aminoglycosyl group significantly changed the pattern of the interactions, indicating that highly polar substituents have a stronger affinity with the exposed regions of the Q-D junction, particularly at the level of the interfacial tetrad.

DeltaGzip: Computing Biopolymer-Ligand Binding Affinity via Kolmogorov Complexity and Lossless Compression

The design of biosequences for biosensing and therapeutics is a challenging multistep search and optimization task. In principle, computational modeling may speed up the design process by virtual screening of sequences based on their binding affinities to target molecules. However, in practice, existing machine-learned models trained to predict binding affinities lack the flexibility with respect to reaction conditions, and molecular dynamics simulations that can incorporate reaction conditions suffer from high computational costs. Here, we describe a computational approach called DeltaGzip that evaluates the free energy of binding in biopolymer-ligand complexes from ultrashort equilibrium molecular dynamics simulations. The entropy of binding is evaluated using the Kolmogorov complexity definition of entropy and approximated using a lossless compression algorithm, Gzip. We benchmark the method on a well-studied data set of protein-ligand complexes comparing the predictions of DeltaGzip to the free energies of binding obtained using Jarzynski equality and experimental measurements.

Amodiaquine Nonspecifically Binds Double Stranded and Three-Way Junction DNA Structures

The 4-aminoquinoline class of compounds includes the important antimalarial compounds amodiaquine and chloroquine. Despite their medicinal importance, the mode of action of these compounds is poorly understood. In a previous study we observed these compounds, as well as quinine and mefloquine, tightly bind the DNA cocaine-binding aptamer. Here, we further explore the range of nucleic acid structures bound by these compounds. To gauge a wide range of binding affinities, we used isothermal titration calorimetry to explore high affinity binding (nM to tens of μM) and NMR spectroscopy to assay weak binding biding in the hundreds of micromolar range. We find that amodiaquine tightly binds all double stranded DNA structures explored. Mefloquine binds double stranded DNA duplex molecules tightly and weakly associates with a three-way junction DNA construct. Quinine and chloroquine only weakly bind duplex DNA but do not tightly bind any of the DNA constructs explored. A simulation of the free energy of binding of these ligands to the Dickerson-Drew dodecamer resulted in an excellent agreement between the simulated and experimental free energy. These results provide new insight into the DNA binding of clinically important antimalarial compounds and may play a role in future development of new antimalarials.

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.

Simulating the Lyotropic Phase Behavior of a Partially Self-Complementary DNA Tetramer

DNA oligomers in solution have been found to develop liquid crystal phases via a hierarchical process that involves Watson-Crick base pairing, supramolecular assembly into columns of duplexes, and long-range ordering. The multiscale nature of this phenomenon makes it difficult to quantitatively describe and assess the importance of the various contributions, particularly for very short strands. We performed molecular dynamics simulations based on the coarse-grained oxDNA model, aiming to depict all of the assembly processes involved and the phase behavior of solutions of the DNA GCCG tetramers. We find good quantitative matching to experimental data at both levels of molecular association (thermal melting) and collective ordering (phase diagram). We characterize the isotropic state and the low-density nematic and high-density columnar liquid crystal phases in terms of molecular order, size of aggregates, and structure, together with their effects on diffusivity processes. We observe a cooperative aggregation mechanism in which the formation of dimers is less thermodynamically favored than the formation of longer aggregates.

The photoactivated dynamics of dGpdC and dCpdG sequences in DNA: a comprehensive quantum mechanical study

Study of alternating DNA GC sequences by different time-resolved spectroscopies has provided fundamental information on the interaction between UV light and DNA, a process of great biological importance. Multiple decay paths have been identified, but their interplay is still poorly understood. Here, we characterize the photophysics of GC-DNA by integrating different computational approaches, to study molecular models including up to 6 bases described at a full quantum mechanical level. Quantum dynamical simulations, exploiting a nonadiabatic linear vibronic coupling (LVC) model, coupled with molecular dynamics sampling of the initial structures of a (GC)5 DNA duplex, provide new insights into the photophysics in the sub-picosecond time-regime. They indicate a substantial population transfer, within 50 fs, from the spectroscopic states towards G → C charge transfer states involving two stacked bases (CTintra), thus explaining the ultrafast disappearance of fluorescence. This picture is consistent with that provided by quantum mechanical geometry optimizations, using time dependent-density functional theory and a polarizable continuum model, which we use to parametrize the LVC model and to map the main excited state deactivation pathways. For the first time, the infrared and excited state absorption signatures of the various states along these pathways are comprehensively mapped. The computational models suggest that the main deactivation pathways, which, according to experiment, lead to ground state recovery on the 10-50 ps time scale, involve CTintra followed by interstrand proton transfer from the neutral G to C-. Our calculations indicate that CTintra is populated to a larger extent and more rapidly in GC than in CG steps and suggest the likely involvement of monomer-like and interstrand charge transfer decay routes for isolated and less stacked CG steps. These findings underscore the importance of the DNA sequence and thermal fluctuations for the dynamics. They will also aid the interpretation of experimental results on other sequences.

Conformational Dynamics of the Nucleosomal Histone H2B Tails Revealed by Molecular Dynamics Simulations

Epigenetic modifications of histone N-terminal tails play a critical role in regulating the chromatin structure and biological processes such as transcription and DNA repair. One of the key post-translational modifications (PTMs) is the acetylation of lysine residues on histone tails. Epigenetic modifications are ubiquitous in the development of diseases, such as cancer and neurological disorders. Histone H2B tails are critical regulators of nucleosome dynamics, biological processes, and certain diseases. Here, we report all-atomistic molecular dynamics (MD) simulations of the nucleosome to demonstrate that acetylation of the histone tails changes their conformational space and interaction with DNA. We perform simulations of H2B tails, critical regulators of gene regulation, in both the lysine-acetylated (ACK) and unacetylated wild type (WT) states. To explore the effects of salt concentration, we use two different NaCl concentrations to perform simulations at microsecond time scales. Salt can modulate the effects of electrostatic interactions between the DNA phosphate backbone and histone tails. Upon acetylation, H2B tails shift their secondary structure helical propensity. The number of contacts between the DNA and the H2B tail decreases. We characterize the conformational dynamics of the H2B tails by principal component analysis (PCA). The ACK tails become more compact at increased salt concentrations, but conformations from the WT tails display the most contacts with DNA at both salt concentrations. Mainly, H2B acetylation may increase the DNA accessibility for regulatory proteins to bind, which can aid in gene regulation and NCP stability.

Amber free energy tools: Interoperable software for free energy simulations using generalized quantum mechanical/molecular mechanical and machine learning potentials

We report the development and testing of new integrated cyberinfrastructure for performing free energy simulations with generalized hybrid quantum mechanical/molecular mechanical (QM/MM) and machine learning potentials (MLPs) in Amber. The Sander molecular dynamics program has been extended to leverage fast, density-functional tight-binding models implemented in the DFTB+ and xTB packages, and an interface to the DeePMD-kit software enables the use of MLPs. The software is integrated through application program interfaces that circumvent the need to perform "system calls" and enable the incorporation of long-range Ewald electrostatics into the external software's self-consistent field procedure. The infrastructure provides access to QM/MM models that may serve as the foundation for QM/MM-ΔMLP potentials, which supplement the semiempirical QM/MM model with a MLP correction trained to reproduce ab initio QM/MM energies and forces. Efficient optimization of minimum free energy pathways is enabled through a new surface-accelerated finite-temperature string method implemented in the FE-ToolKit package. Furthermore, we interfaced Sander with the i-PI software by implementing the socket communication protocol used in the i-PI client-server model. The new interface with i-PI allows for the treatment of nuclear quantum effects with semiempirical QM/MM-ΔMLP models. The modular interoperable software is demonstrated on proton transfer reactions in guanine-thymine mispairs in a B-form deoxyribonucleic acid helix. The current work represents a considerable advance in the development of modular software for performing free energy simulations of chemical reactions that are important in a wide range of applications.

Electronic and Nuclear Quantum Effects on Proton Transfer Reactions of Guanine-Thymine (G-T) Mispairs Using Combined Quantum Mechanical/Molecular Mechanical and Machine Learning Potentials

Rare tautomeric forms of nucleobases can lead to Watson-Crick-like (WC-like) mispairs in DNA, but the process of proton transfer is fast and difficult to detect experimentally. NMR studies show evidence for the existence of short-time WC-like guanine-thymine (G-T) mispairs; however, the mechanism of proton transfer and the degree to which nuclear quantum effects play a role are unclear. We use a B-DNA helix exhibiting a wGT mispair as a model system to study tautomerization reactions. We perform ab initio (PBE0/6-31G*) quantum mechanical/molecular mechanical (QM/MM) simulations to examine the free energy surface for tautomerization. We demonstrate that while the ab initio QM/MM simulations are accurate, considerable sampling is required to achieve high precision in the free energy barriers. To address this problem, we develop a QM/MM machine learning potential correction (QM/MM-ΔMLP) that is able to improve the computational efficiency, greatly extend the accessible time scales of the simulations, and enable practical application of path integral molecular dynamics to examine nuclear quantum effects. We find that the inclusion of nuclear quantum effects has only a modest effect on the mechanistic pathway but leads to a considerable lowering of the free energy barrier for the GT*⇌G*T equilibrium. Our results enable a rationalization of observed experimental data and the prediction of populations of rare tautomeric forms of nucleobases and rates of their interconversion in B-DNA.

Oligo-Adenine and Cyanuric Acid Supramolecular DNA-Based Hydrogels Exhibiting Acid-Resistance and Physiological pH-Responsiveness

Expanding the functions and applications of DNA by integrating noncanonical bases and structures into biopolymers is a continuous scientific effort. An adenine-rich strand (A-strand) is introduced as functional scaffold revealing, in the presence of the low-molecular-weight cofactor cyanuric acid (CA, pKa 6.9), supramolecular hydrogel-forming efficacies demonstrating multiple pH-responsiveness. At pH 1.2, the A-strand transforms into a parallel A-motif duplex hydrogel cross-linked by AH+-H+A units due to the protonation of adenine (pKa 3.5). At pH 5.2, and in the presence of coadded CA, a helicene-like configuration is formed between adenine and protonated CA, generating a parallel A-CA triplex cross-linked hydrogel. At pH 8.0, the hydrogel undergoes transition into a liquid state by deprotonation of CA cofactor units and disassembly of A-CA triplex into its constituent components. Density functional theory calculations and molecular dynamics simulations, supporting the structural reconfigurations of A-strand in the presence of CA, are performed. The sequential pH-stimulated hydrogel states are rheometrically characterized. The hydrogel framework is loaded with fluorescein-labeled insulin, and the pH-stimulated release of insulin from the hydrogel across the pH barriers present in the gastrointestinal tract is demonstrated. The results provide principles for future application of the hydrogel for oral insulin administration for diabetes.

Multiscale modeling reveals the ion-mediated phase separation of nucleosome core particles

Due to the vast length scale inside the cell nucleus, multiscale models are required to understand chromatin folding, structure, and dynamics and how they regulate genomic activities such as DNA transcription, replication, and repair. We study the interactions and structure of condensed phases formed by the universal building block of chromatin, the nucleosome core particle (NCP), using bottom-up multiscale coarse-grained (CG) simulations with a model extracted from all-atom MD simulations. In the presence of the multivalent cations Mg(H2O)62+ or CoHex3+, we analyze the internal structures of the NCP aggregates and the contributions of histone tails and ions to the aggregation patterns. We then derive a "super" coarse-grained (SCG) NCP model to study the macroscopic scale phase separation of NCPs. The SCG simulations show the formation of NCP aggregates with Mg(H2O)62+ concentration-dependent densities and sizes. Variation of the CoHex3+ concentrations results in highly ordered lamellocolumnar and hexagonal columnar phases in agreement with experimental data. The results give detailed insights into nucleosome interactions and for understanding chromatin folding in the cell nucleus.

An interaction network in the polymerase active site is a prerequisite for Watson-Crick base pairing in Pol γ

The replication accuracy of DNA polymerase gamma (Pol γ) is essential for mitochondrial genome integrity. Mutation of human Pol γ arginine-853 has been linked to neurological diseases. Although not a catalytic residue, Pol γ arginine-853 mutants are void of polymerase activity. To identify the structural basis for the disease, we determined a crystal structure of the Pol γ mutant ternary complex with correct incoming nucleotide 2'-deoxycytidine 5'-triphosphate (dCTP). Opposite to the wild type that undergoes open-to-closed conformational changes when bound to a correct nucleotide that is essential for forming a catalytically competent active site, the mutant complex failed to undergo the conformational change, and the dCTP did not base pair with its Watson-Crick complementary templating residue. Our studies revealed that arginine-853 coordinates an interaction network that aligns the 3'-end of primer and dCTP with the catalytic residues. Disruption of the network precludes the formation of Watson-Crick base pairing and closing of the active site, resulting in an inactive polymerase.

Effect of temperature on anisotropic bending elasticity of dsRNA: an all-atom molecular dynamics simulation

Employing all-atom molecular dynamics simulations, we examined the temperature-dependent behavior of bending elasticity in double-stranded RNA (dsRNA). Specifically, we focused on the bending persistence length and its constituent components, namely, the tilt and roll stiffness. Our results revealed a near-linear decrease in these stiffness components as a function of temperature, thereby highlighting the increased flexibility of dsRNA at elevated temperatures. Furthermore, our data revealed a significant anisotropy in dsRNA bending elasticity, which diminished with increasing temperature, attributable to marked disparities in tilt and roll stiffness components. We delineated the underlying biophysical mechanisms and corroborated our findings with extant literature. These observations offer salient implications for advancing our understanding of nucleic acid elasticity, and are pertinent to potential medical applications.

XMECP: Reaching State-of-the-Art MECP Optimization in Multiscale Complex Systems

The Python-based program, XMECP, is developed for realizing robust, efficient, and state-of-the-art minimum energy crossing point (MECP) optimization in multiscale complex systems. This article introduces the basic capabilities of the XMECP program by theoretically investigating the MECP mechanism of several example systems including (1) the photosensitization mechanism of benzophenone, (2) photoinduced proton-coupled electron transfer in the cytosine-guanine base pair in DNA, (3) the spin-flip process in oxygen activation catalyzed by an iron-containing 2-oxoglutarate-dependent oxygenase (Fe/2OGX), and (4) the photochemical pathway of flavoprotein adjusted by the intensity of an external electric field. MECPs related to multistate reaction and multistate reactivity in large-scale complex biochemical systems can be well-treated by workflows suggested by the XMECP program. The branching plane updating the MECP optimization algorithm is strongly recommended as it provides derivative coupling vector (DCV) with explicit calculation and can equivalently evaluate contributions from non-QM residues to DCV, which can be nonadiabatic coupling or spin-orbit coupling in different cases. In the discussed QM/MM examples, we also found that the influence on the QM region by DCV can occur through noncovalent interactions and decay with distance. In the example of DNA base pairs, the nonadiabatic coupling occurs across the π-π stacking structure formed in the double-helix system. In contrast to general intuition, in the example of Fe/2OGX, the central ferrous and oxygen part contribute little to the spin-orbit coupling; however, a nearby arginine residue, which is treated by molecular mechanics in the QM/MM method, contributes significantly via two hydrogen bonds formed with α-ketoglutarate (α-KG). This indicates that the arginine residue plays a significant role in oxygen activation, driving the initial triplet state toward the productive quintet state, which is more than the previous knowledge that the arginine residue can bind α-KG at the reaction site by hydrogen bonds.

Directional Pumpless Transport of Biomolecules through Self-Propelled Nanodroplets on Patterned Heterostructures

The surface patterning in natural systems has exhibited appreciable functional advantages for life activities, which serve as inspiration for the design of artificial counterparts to achieve functions such as directional liquid transport at the nanoscale. Here, we propose a patterned two-dimensional (2D) in-plane heterostructure with a triangle-shaped hexagonal boron nitride (hBN) track embedded in graphene nanosheets, which can achieve unidirectional and self-propelled transport of nanodroplets carrying various biomolecules such as DNA, RNA, and peptides. Our extensive MD simulations show that the wettability gradient on the patterned heterostructure can drive the motion of nanodroplet with an instantaneous acceleration, which also permits long-distance transport (>100 nm) at the microsecond time scale. The different behaviors of various types of biomolecules have been further studied systematically within the transporting nanodroplets. These findings suggest that these specially designed, patterned heterostructures have the potential for spontaneous, directional transport of important biomolecules, which might be useful in biosensing, drug delivery, and biomedical nanodevices.

Systematic assessment of the flexibility of uracil damaged DNA

Uracil is a common DNA lesion which is recognized and removed by uracil DNA-glycosylase (UDG) as a part of the base excision repair pathway. Excision proceeds by base flipping, and UDG efficiency is thought to depend on the ease of deformability of the bases neighboring the lesion. We used molecular dynamics simulations to assess the flexibility of a large library of dsDNA strands, containing all tetranucleotide motifs with U:A, U:G, T:A or C:G base pairs. Our study demonstrates that uracil damaged DNA largely follows trends in flexibility of undamaged DNA. Measured bending persistence lengths, groove widths, step parameters and base flipping propensities demonstrate that uracil increases the flexibility of DNA, and that U:G base paired strands are more flexible than U:A strands. Certain sequence contexts are more deformable than others, with a key role for the 3' base next to uracil. Flexibilities are large when this base is an A or G, and repressed for a C or T. A 5' T adjacent to the uracil strongly promotes flexibility, but other 5' bases are less influential. DNA bending is correlated to step deformations and base flipping, and bending aids flipping. Our study implies that the link between substrate flexibility and UDG efficiency is widely valid, helps explain why UDG prefers to bind U:G base paired strands, and suggests that the DNA bending angle of the UDG-substrate complex is optimal for base flipping.Communicated by Ramaswamy H. Sarma.

Development of the 12-Base Short Dimeric Myogenetic Oligodeoxynucleotide That Induces Myogenic Differentiation

A myogenetic oligodeoxynucleotide (myoDN), iSN04 (5'-AGA TTA GGG TGA GGG TGA-3'), is a single-stranded 18-base telomeric DNA that serves as an anti-nucleolin aptamer and induces myogenic differentiation, which is expected to be a nucleic acid drug for the prevention of disease-associated muscle wasting. To improve the drug efficacy and synthesis cost of myoDN, shortening the sequence while maintaining its structure-based function is a major challenge. Here, we report the novel 12-base non-telomeric myoDN, iMyo01 (5'-TTG GGT GGG GAA-3'), which has comparable myogenic activity to iSN04. iMyo01 as well as iSN04 promoted myotube formation of primary-cultured human myoblasts with upregulation of myogenic gene expression. Both iMyo01 and iSN04 interacted with nucleolin, but iMyo01 did not bind to berberine, the isoquinoline alkaloid that stabilizes iSN04. Nuclear magnetic resonance revealed that iMyo01 forms a G-quadruplex structure despite its short sequence. Native polyacrylamide gel electrophoresis and a computational molecular dynamics simulation indicated that iMyo01 forms a homodimer to generate a G-quadruplex. These results provide new insights into the aptamer truncation technology that preserves aptamer conformation and bioactivity for the development of efficient nucleic acid drugs.

Identification of a novel FOXO3 agonist that protects against alcohol induced liver injury

Alcohol-related liver disease (ALD) is a global healthcare concern which caused by excessive alcohol consumption with limited treatment options. The pathogenesis of ALD is complex and involves in hepatocyte damage, hepatic inflammation, increased gut permeability and microbiome dysbiosis. FOXO3 is a well-recognized transcription factor which associated with longevity via promoting antioxidant stress response, preventing senescence and cell death, and inhibiting inflammation. We and many others have reported that FOXO3-/- mice develop more severe liver injury in response to alcohol. In the present study, we aimed to develop compounds that activate FOXO3 and further investigate their effects in alcohol induced liver injury. Through virtual screening, we discovered series of small molecular compounds that showed high affinity to FOXO3. We confirmed effects of compounds on FOXO3 target gene expression, as well as antioxidant and anti-apoptotic effects in vitro. Subsequently we evaluated the protective efficacy of compounds in alcohol induced liver injury in vivo. As a result, the leading compound we identified, 214991, activated downstream target genes expression of FOXO3, inhibited intracellular ROS accumulation and cell apoptosis induced by H2O2 and sorafenib. By using Lieber-DeCarli alcohol feeding mouse model, 214991 showed protective effects against alcohol-induced liver inflammation, macrophage and neutrophil infiltration, and steatosis. These findings not only reinforce the potential of FOXO3 as a valuable target for therapeutic intervention of ALD, but also suggested that compound 214991 as a promising candidate for the development of innovative therapeutic strategies of ALD.

An Ensemble Docking Approach for Analyzing and Designing Aptamer Heterodimers Targeting VEGF165

Vascular endothelial growth factor 165 (VEGF165) is a prominent isoform of the VEGF-A protein that plays a crucial role in various angiogenesis-related diseases. It is homodimeric, and each of its monomers is composed of two domains connected by a flexible linker. DNA aptamers, which have emerged as potent therapeutic molecules for many proteins with high specificity and affinity, can also work for VEGF165. A DNA aptamer heterodimer composed of monomers of V7t1 and del5-1 connected by a flexible linker (V7t1:del5-1) exhibits a greater binding affinity with VEGF165 compared to either of the two monomers alone. Although the structure of the complex formed between the aptamer heterodimer and VEGF165 is unknown due to the highly flexible linkers, gaining structural information will still be valuable for future developments. Toward this end of accessing structural information, we adopt an ensemble docking approach here. We first obtain an ensemble of structures for both VEGF165 and the aptamer heterodimer by considering both small- and large-scale motions. We then proceed through an extraction process based on ensemble docking, molecular dynamics simulations, and binding free energy calculations to predict the structures of the VEGF165/V7t1:del5-1 complex. Through the same procedures, we reach a new aptamer heterodimer that bears a locked nucleic acid-modified counterpart of V7t1, namely RNV66:del5-1, which also binds well with VEGF165. We apply the same protocol to the monomeric units V7t1, RNV66, and del5-1 to target VEGF165. We observe that V7t1:del5-1 and RNV66:del5-1 show higher binding affinities with VEGF165 than any of the monomers, consistent with experiments that support the notion that aptamer heterodimers are more effective anti-VEGF165 aptamers than monomeric aptamers. Among the five different aptamers studied here, the newly designed RNV66:del5-1 shows the highest binding affinity with VEGF165. We expect that our ensemble docking approach can help in de novo designs of homo/heterodimeric anti-angiogenic drugs to target the homodimeric VEGF165.

Deciphering the Lexicon of Protein Targets: A Review on Multifaceted Drug Discovery in the Era of Artificial Intelligence

Understanding protein sequence and structure is essential for understanding protein-protein interactions (PPIs), which are essential for many biological processes and diseases. Targeting protein binding hot spots, which regulate signaling and growth, with rational drug design is promising. Rational drug design uses structural data and computational tools to study protein binding sites and protein interfaces to design inhibitors that can change these interactions, thereby potentially leading to therapeutic approaches. Artificial intelligence (AI), such as machine learning (ML) and deep learning (DL), has advanced drug discovery and design by providing computational resources and methods. Quantum chemistry is essential for drug reactivity, toxicology, drug screening, and quantitative structure-activity relationship (QSAR) properties. This review discusses the methodologies and challenges of identifying and characterizing hot spots and binding sites. It also explores the strategies and applications of artificial-intelligence-based rational drug design technologies that target proteins and protein-protein interaction (PPI) binding hot spots. It provides valuable insights for drug design with therapeutic implications. We have also demonstrated the pathological conditions of heat shock protein 27 (HSP27) and matrix metallopoproteinases (MMP2 and MMP9) and designed inhibitors of these proteins using the drug discovery paradigm in a case study on the discovery of drug molecules for cancer treatment. Additionally, the implications of benzothiazole derivatives for anticancer drug design and discovery are deliberated.

Why Are Left-Handed G-Quadruplexes Scarce?

G-quadruplexes (G4s) are nucleic acid structures crucial for the regulation of gene expression and genome maintenance. While they hold promise as nanodevice components, achieving desired G4 folds requires understanding the interplay between stability and structural properties, like helicity. Although right-handed G4 structures dominate the experimental data, the molecular basis for this preference over left-handed helicity is unclear. To address this, we employ all-atom molecular dynamics simulations and quantum chemical methods. Our results reveal that right-handed G4s exhibit greater thermodynamic and kinetic stability as a result of favorable sugar-phosphate backbone conformations in guanine tracts. Moreover, while hydrogen-bonding patterns influence helicity-specific G4 loop conformations, they minimally affect stability differences. We also elucidate the strong correlation between helicity and the strand progression direction, essential for G4 structures. These findings deepen our understanding of G4s, providing molecular-level insights into their structural and energetic preferences, which could inform the design of novel nanodevices.

Elucidation of the catalytic mechanism of a single-metal dependent homing endonuclease using QM and QM/MM approaches: the case study of I-PpoI

Homing endonucleases (HEs) are highly specific DNA cleaving enzymes, with I-PpoI having been suggested to use a single metal to accelerate phosphodiester bond cleavage. Although an I-PpoI mechanism has been proposed based on experimental structural data, no consensus has been reached regarding the roles of the metal or key active site amino acids. This study uses QM cluster and QM/MM calculations to provide atomic-level details of the I-PpoI catalytic mechanism. Minimal QM cluster and large-scale QM/MM models demonstrate that the experimentally-proposed pathway involving direct Mg2+ coordination to the substrate coupled with leaving group protonation through a metal-activated water is not feasible due to an inconducive I-PpoI active site alignment. Despite QM cluster models of varying size uncovering a pathway involving leaving group protonation by a metal-activated water, indirect (water-mediated) metal coordination to the substrate is required to afford this pathway, which renders this mechanism energetically infeasible. Instead, QM cluster models reveal that the preferred pathway involves direct Mg2+-O3' coordination to stabilize the charged substrate and assist leaving group departure, while H98 activates the water nucleophile. These calculations also underscore that both catalytic residues that directly interact with the substrate and secondary amino acids that position or stabilize these residues are required for efficient catalysis. QM/MM calculations on the solvated enzyme-DNA complex verify the preferred mechanism, which is fully consistent with experimental kinetic, structural, and mutational data. The fundamental understanding of the I-PpoI mechanism of action, gained from the present work can be used to further explore potential uses of this enzyme in biotechnology and medicine, and direct future computational investigations of other members of the understudied HE family.

Systematic Comparison of Atomistic Force Fields for the Mechanical Properties of Double-Stranded DNA

The response of double-stranded DNA to external mechanical stress plays a central role in its interactions with the protein machinery in the cell. Modern atomistic force fields have been shown to provide highly accurate predictions for the fine structural features of the duplex. In contrast, and despite their pivotal function, less attention has been devoted to the accuracy of the prediction of the elastic parameters. Several reports have addressed the flexibility of double-stranded DNA via all-atom molecular dynamics, yet the collected information is insufficient to have a clear understanding of the relative performance of the various force fields. In this work, we fill this gap by performing a systematic study in which several systems, characterized by different sequence contexts, are simulated with the most popular force fields within the AMBER family, bcs1 and OL15, as well as with CHARMM36. Analysis of our results, together with their comparison with previous work focused on bsc0, allows us to unveil the differences in the predicted rigidity between the newest force fields and suggests a roadmap to test their performance against experiments. In the case of the stretch modulus, we reconcile these differences, showing that a single mapping between sequence-dependent conformation and elasticity via the crookedness parameter captures simultaneously the results of all force fields, supporting the key role of crookedness in the mechanical response of double-stranded DNA.

Melting and Bubble Formation in a Double-Stranded DNA: Microscopic Aspects of Early Base-Pair Opening Events and the Role of Water

Despite its rigid structure, DNA is a remarkably flexible molecule. Flexibility is essential for biological functions (such as transcription and gene repair), which require large-amplitude structural changes such as bubble formation. The bubbles thus formed are required to have a certain stability of their own and survive long on the time scale of molecular motions. A molecular understanding of fluctuations leading to quasi-stable structures is not available. Through extensive atomistic molecular dynamics simulations, we identify a sequence of microscopic events that culminate in local bubble formation, which is initiated by base-pair (BP) opening, resulting from the cleavage of native BP hydrogen bonds (HBs). This is followed by the formation of mismatched BPs with non-native contacts. These metastable structures can either revert to their original forms or undergo a flipping transition to form a local bubble that can span across 3-4 BPs. A substantial distortion of the DNA backbone and a disruption of BP stacking are observed because of the structural changes induced by these local perturbations. We also explored how water helps in the entire process. A small number of water molecules undergo rearrangement to stabilize the intermediate states by forming HBs with DNA bases. Water thus acts as a lubricant that counteracts the enthalpic penalty suffered from the loss of native BP contacts. Although the process of bubble formation is reversible, the sequence of steps involved poses an entropic barrier, preventing it from easily retracing the path to the native state.

Influence of temperature on bend, twist and twist-bend coupling of dsDNA

The temperature-dependent bend and twist elasticities of dsDNA, as well as their couplings, were explored through all-atom molecular dynamics simulations. Three rotational parameters, tilt, roll, and twist, were employed to assess the bend and twist elasticities through their stiffness matrix. Our analysis indicates that the bend and twist stiffnesses decrease as the temperature rises, primarily owing to entropic influences stemming from thermodynamic fluctuations. Furthermore, the couplings between these rotational parameters also exhibit a decline with increasing temperature, although the roll-twist coupling displays greater strength than the tilt-roll and tilt-twist couplings, attributed to its more robust correction component. We elucidated the influence of temperature on bend and twist elasticities based on the comparisons between various models and existing data.

Temperature-dependent elasticity of DNA, RNA, and hybrid double helices

Nucleic acid double helices in their DNA, RNA, and DNA-RNA hybrid form play a fundamental role in biology and are main building blocks of artificial nanostructures, but how their properties depend on temperature remains poorly understood. Here, we report thermal dependence of dynamic bending persistence length, twist rigidity, stretch modulus, and twist-stretch coupling for DNA, RNA, and hybrid duplexes between 7°C and 47°C. The results are based on all-atom molecular dynamics simulations using different force field parameterizations. We first demonstrate that unrestrained molecular dynamics can reproduce experimentally known mechanical properties of the duplexes at room temperature. Beyond experimentally known features, we also infer the twist rigidity and twist-stretch coupling of the hybrid duplex. As for the temperature dependence, we found that increasing temperature softens all the duplexes with respect to bending, twisting, and stretching. The relative decrease of the stretch moduli is 0.003-0.004/°C, similar for all the duplex variants despite their very different stretching stiffness, whereas RNA twist stiffness decreases by 0.003/°C, and smaller values are found for the other elastic moduli. The twist-stretch couplings are nearly unaffected by temperature. The stretching, bending, and twisting stiffness all include an important entropic component. Relation of our results to the two-state model of DNA flexibility is discussed. Our work provides temperature-dependent elasticity of nucleic acid duplexes at the microsecond scale relevant for initial stages of protein binding.

Investigation of the Flipping Dynamics of 1, N6-Ethenoadenine in Alkyladenine DNA Glycosylase

Alkyladenine DNA glycosylase (AAG) is an essential enzyme responsible for maintaining genome integrity by repairing several DNA lesions damaged by alkylation or deamination. Understanding how it can recognize and excise the lesions thus lays the foundation for therapeutic treatment against lesion-associated diseases or cancers. However, the molecular details of how the lesion can be distinguished from the matched base by AAG and how it enters the cleavage site, ready for excision, are not fully elucidated. In this study, we have revealed the molecular details of the flipping dynamics of 1, N6-ethenoadenine (εA) not only in the form of free double-stranded DNA (dsDNA) but also in the form of the AAG-dsDNA complex. Our MD simulations and PMF calculations have shown that the flipping of εA and dA is thermodynamically disfavored in the free dsDNA, even though εA has a lower flipping energy barrier than dA. By sharp contrast, the flipping of εA is thermodynamically favored in AAG with an obvious free energy drop, while dA is equally stabilized before and after the flipping. Moreover, a comparison of the PMFs in the forms of free dsDNA and the AAG-dsDNA complex has pinpointed the role of AAG in discriminating εA against dA and facilitating the flipping of εA. Besides, the flipping process is simulated along the major and minor grooves, and our results have additionally demonstrated that the flipping is not directional in the free dsDNA while flipping along the major groove is kinetically more favorable than the minor groove in the AAG-dsDNA complex. Overall, our study has offered molecular insights into the flipping dynamics of εA and revealed its discrimination mechanism by AAG, which is expected to guide further enzyme engineering for therapeutic applications.

Estrone-Based Derivatives Stabilize the c-MYC and c-KIT G-Quadruplex DNA Structures

G-rich sequences are present across the genome and can fold to form dynamic secondary structures, namely, G-quadruplexes (G4). These structures play a pivotal role in regulating numerous biological processes including replication, transcription, and translation. Therefore, targeting these structures using molecular scaffolds is an attractive approach to modulating their functions. Herein, we report the synthesis of three estrone-based derivatives (Est-1, Est-2, and Est-3) with a nonplanar core and a cationic alkyl side chain as G4 stabilizers. CD melting and polymerase stop assay results indicate that these ligands preferentially stabilize parallel c-MYC and c-KIT1 G4s over the other G4s and duplex DNAs. The ligand Est-3 shows cytotoxicity against cancer cell lines and effectively downregulates the c-KIT gene in HepG2 cell lines. Molecular modeling and dynamics studies showed that the ligand prefers stacking over the 5'-quartet of c-MYC G4 using the aromatic ring of the ligand. Overall, the findings of this study demonstrate that even G4 ligands can accommodate nonplanar scaffolds, which opens up new avenues for ligand design.

Copper(II) coordination compounds based on bis-hydrazones of 2,6-diacetylpyridine: synthesis, structure, and cytotoxic activity

By reacting a series of 2,6-diacetylpyridine bis-hydrazones containing pyrimidine (H2L1), benzimidazole (H2L2) and phthalazine (H2L3) heterocyclic fragments with copper(II) chloride and bromide, a variety of pentacoordinated complexes of the composition [Cu(H2L1)X]X, [Cu(HL2)X] and [Cu(HL3)X], where X = Cl-, Br-, are formed. The properties and structure of the compounds were studied by means of NMR, IR, UV-vis, ESR, and X-ray absorption spectroscopy, cyclic voltammetry and X-Ray single crystal diffraction methods. It was shown that complexes of the cationic type [Cu(H2L1)X]X have an asymmetric structure with a distorted square-pyramidal geometry of the coordination unit. The coordination polyhedron of metal chelates [Cu(HL2)X] and [Cu(HL3)X] is an almost ideal square pyramid. Investigations of the cytotoxic activity of the obtained compounds in vitro on human hepatocellular carcinoma (HepG2) and non-tumor human lung fibroblast (MRC-5) cell lines demonstrated that complexes show higher activity compared with the well-known anticancer agent cisplatin. In addition, metal chelates [Cu(H2L1)Cl]Cl, [Cu(HL2)Cl], [Cu(HL2)Br] and [Cu(HL3)Cl] were less toxic to non-tumor cells MRC-5. A study of the binding of complexes to bovine serum albumin (BSA) protein using fluorescence spectroscopy showed that copper complexes are strongly bound to BSA. To study the mechanism of interaction of the complexes with the DNA of cancer cells, molecular dynamics simulation of the compound [Cu(HL3)Cl] was carried out. It was shown that the complex enters into π-stacking interactions predominantly with adenine and thymine bases.

Is Metal Stabilization of the Leaving Group Required or Can Lysine Facilitate Phosphodiester Bond Cleavage in Nucleic Acids? A Computational Study of EndoV

Endonuclease V (EndoV) is a single-metal-dependent enzyme that repairs deaminated DNA nucleobases in cells by cleaving the phosphodiester bond, and this enzyme has proven to be a powerful tool in biotechnology and medicine. The catalytic mechanism used by EndoV must be understood to design new disease detection and therapeutic solutions and further exploit the enzyme in interdisciplinary applications. This study has used a mixed molecular dynamics (MD) and quantum mechanics/molecular mechanics (QM/MM) approach to compare eight distinct catalytic pathways and provides the first proposed mechanism for bacterial EndoV. The calculations demonstrate that mechanisms involving either direct or indirect metal coordination to the leaving group of the substrate previously proposed for other nucleases are unlikely for EndoV, regardless of the general base (histidine, aspartate, and substrate phosphate moiety). Instead, distinct catalytic pathways are characterized for EndoV that involve K139 stabilizing the leaving group, a metal-coordinated water stabilizing the transition structure, and either H214 or a substrate phosphate group activating the water nucleophile. In silico K139A and H214A mutational results support the newly proposed roles of these residues. Although this is a previously unseen combination of general base, general acid, and metal-binding architecture for a one-metal-dependent endonuclease, our proposed catalytic mechanisms are fully consistent with experimental kinetic, structural, and mutational data. In addition to substantiating a growing body of literature, suggesting that one metal is enough to catalyze P-O bond cleavage in nucleic acids, this new fundamental understanding of the catalytic function will promote the exploration of new and improved applications of EndoV.

DNA Quadruplex Structure with a Unique Cation Dependency

DNA quadruplex structures provide an additional layer of regulatory control in genome maintenance and gene expression and are widely used in nanotechnology. We report the discovery of an unprecedented tetrastranded structure formed from a native G-rich DNA sequence originating from the telomeric region of Caenorhabditis elegans. The structure is defined by multiple properties that distinguish it from all other known DNA quadruplexes. Most notably, the formation of a stable so-called KNa-quadruplex (KNaQ) requires concurrent coordination of K+ and Na+ ions at two distinct binding sites. This structure provides novel insight into G-rich DNA folding under ionic conditions relevant to eukaryotic cell physiology and the structural evolution of telomeric DNA. It highlights the differences between the structural organization of human and nematode telomeric DNA, which should be considered when using C. elegans as a model in telomere biology, particularly in drug screening applications. Additionally, the absence/presence of KNaQ motifs in the host/parasite introduces an intriguing possibility of exploiting the KNaQ fold as a plausible antiparasitic drug target. The structure's unique shape and ion dependency and the possibility of controlling its folding by using low-molecular-weight ligands can be used for the design or discovery of novel recognition DNA elements and sensors.

Design and Characterization of Compact, Programmable, Multistranded Nonimmunostimulatory Nucleic Acid Nanoparticles Suitable for Biomedical Applications

We report a thorough investigation of the role of single-stranded thymidine (ssT) linkers in the stability and flexibility of minimal, multistranded DNA nanostructures. We systematically explore the impact of varying the number of ssTs in three-way junction motifs (3WJs) on their formation and properties. Through various UV melting experiments and molecular dynamics simulations, we demonstrate that while the number of ssTs minimally affects thermodynamic stability, the increasing ssT regions significantly enhance the structural flexibility of 3WJs. Utilizing this knowledge, we design triangular DNA nanoparticles with varying ssTs, all showing exceptional assembly efficiency except for the 0T triangle. All triangles demonstrate enhanced stability in blood serum and are nonimmunostimulatory and nontoxic in mammalian cell lines. The 4T 3WJ is chosen as the building block for constructing other polygons due to its enhanced flexibility and favorable physicochemical characteristics, making it a versatile choice for creating cost-effective, stable, and functional DNA nanostructures that can be stored in the dehydrated forms while retaining their structures. Our study provides valuable insights into the design and application of nucleic acid nanostructures, emphasizing the importance of understanding stability and flexibility in the realm of nucleic acid nanotechnology. Our findings suggest the intricate connection between these ssTs and the structural adaptability of DNA 3WJs, paving the way for more precise design and engineering of nucleic acid nanosystems suitable for broad biomedical applications.

Non-covalent binding interaction of bioactive coumarin esculetin with calf thymus DNA and yeast transfer RNA: A detailed investigation to decipher the binding affinities, binding location, interacting forces and structural alterations at a molecular level

Esculetin is a well-known coumarin derivative found abundantly in nature possessing an extensive array of pharmacological and therapeutic properties. Consequently, to comprehend its molecular recognition mechanism, our objective is to conduct a complete investigation of its interactions with the nucleic acid, specifically ct-DNA, and t-RNA, using spectroscopic and computational techniques. The intrinsic fluorescence of esculetin is quenched when it interacts with ct-DNA and t-RNA, and this occurs through a static quenching mechanism. The thermodynamic parameters demonstrated that the interaction is influenced by hydrogen bonding and weak van der Waals forces. CD and FT-IR results revealed no conformational changes in ct-DNA and t-RNA structure on binding with esculetin. Furthermore, competitive displacement assay with ethidium bromide, melting temperature, viscosity measurement, and potassium iodide quenching experiments, reflected that esculetin probably binds to the minor groove of ct-DNA. The molecular docking results provided further confirmation for the spectroscopic findings, including the binding location of esculetin and binding energies of esculetin complexes with ct-DNA and t-RNA. Molecular dynamics simulation studies demonstrated the conformational stability and flexibility of nucleic acids.

Computed interactions of berenil with restricted foldamers of c-MYC DNA G-quadruplexes

G-quadruplexes (G4s) are secondary four-stranded DNA helical structures made up of guanine-rich nucleic acids that can assemble in the promoter regions of multiple genes under the appropriate conditions. Stabilization of G4 structures by small molecules can regulate transcription in non-telomeric regions, including in proto-oncogenes and promoter regions, contributing to anti-proliferative and anti-tumor activities. Because G4s are detectable in cancer cells but not in normal cells, they make excellent drug discovery targets. Diminazene, DMZ (or berenil), has been shown to be an efficient G-quadruplex binder. Due to the stability of the folding topology, G-quadruplex structures are frequently found in the promotor regions of oncogenes and may play a regulatory role in gene activation. Using molecular docking and molecular dynamics simulations on several different binding poses, we have studied DMZ binding toward multiple G4 topologies of the c-MYC G-quadruplex. DMZ binds preferentially to G4s that have extended loops and flanking bases. This preference arises from its interactions with the loops and the flanking nucleotides, which were not found in the structure lacking extended regions. The binding to the G4s with no extended regions instead occurred mostly through end stacking. All binding sites for DMZ were confirmed by 100 ns molecular dynamics simulations and through binding enthalpies calculated using the MM-PBSA method. The primary driving forces were electrostatic, as the cationic DMZ interacts with the anionic phosphate backbone, and through van der Waals interactions, which primarily contributed in end stacking interactions.Communicated by Ramaswamy H. Sarma.

Temperature-Dependent Twist of Double-Stranded RNA Probed by Magnetic Tweezer Experiments and Molecular Dynamics Simulations

RNA plays critical roles in the transmission and regulation of genetic information and is increasingly used in biomedical and biotechnological applications. Functional RNAs contain extended double-stranded regions, and the structure of double-stranded RNA (dsRNA) has been revealed at high resolution. However, the dependence of the properties of the RNA double helix on environmental effects, notably temperature, is still poorly understood. Here, we use single-molecule magnetic tweezer measurements to determine the dependence of the dsRNA twist on temperature. We find that dsRNA unwinds with increasing temperature, even more than DNA, with ΔTwRNA = -14.4 ± 0.7°/(°C·kbp), compared to ΔTwDNA = -11.0 ± 1.2°/(°C·kbp). All-atom molecular dynamics (MD) simulations using a range of nucleic acid force fields, ion parameters, and water models correctly predict that dsRNA unwinds with rising temperature but significantly underestimate the magnitude of the effect. These MD data, together with additional MD simulations involving DNA and DNA-RNA hybrid duplexes, reveal a linear correlation between the twist temperature decrease and the helical rise, in line with DNA but at variance with RNA experimental data. We speculate that this discrepancy might be caused by some unknown bias in the RNA force fields tested or by as yet undiscovered transient alternative structures in the RNA duplex. Our results provide a baseline to model more complex RNA assemblies and to test and develop new parametrizations for RNA simulations. They may also inspire physical models of the temperature-dependent dsRNA structure.

van der Waals Parameter Scanning with Amber Nucleic Acid Force Fields: Revisiting Means to Better Capture the RNA/DNA Structure through MD

Molecular dynamics simulations can be used in combination with experimental techniques to uncover the intricacies of biomolecular structure, dynamics, and the resulting interactions. However, many noncanonical nucleic acid structures have proven to be challenging to replicate in accurate agreement with experimental data, often attributed to known force field deficiencies. A common force field criticism is the handling of van der Waals (vdW) parameters, which have not been updated since the regular use of Ewald's methods became routine. This work dives into the effects of minute vdW radii shifts on RNA tetranucleotide, B-DNA, and Z-DNA model systems described by commonly used Amber force fields. Using multidimensional replica exchange molecular dynamics (M-REMD), the GACC RNA tetranucleotide demonstrated changes in the structural distribution between the NMR minor and anomalous structure populations based on the O2' vdW radii scanning. However, no significant change in the NMR Major conformation population was observed. There were minimal changes in the B-DNA structure but there were more substantial improvements in Z-DNA structural descriptions, specifically with the Tumuc1 force field. This occurred with both LJbb vdW radii adjustments and incorporation of the CUFIX nonbonded parameter modifications. Though the limited vdW modifications tested did not provide a universal fix to the challenge of simulating the various known nucleic acid structures, they do provide direction and a greater understanding for future force field development efforts.

AMBERff at Scale: Multimillion-Atom Simulations with AMBER Force Fields in NAMD

All-atom molecular dynamics (MD) simulations are an essential structural biology technique with increasing application to multimillion-atom systems, including viruses and cellular machinery. Classical MD simulations rely on parameter sets, such as the AMBER family of force fields (AMBERff), to accurately describe molecular motion. Here, we present an implementation of AMBERff for use in NAMD that overcomes previous limitations to enable high-performance, massively parallel simulations encompassing up to two billion atoms. Single-point potential energy comparisons and case studies on model systems demonstrate that the implementation produces results that are as accurate as running AMBERff in its native engine.

Structural basis of antiphage immunity generated by a prokaryotic Argonaute-associated SPARSA system

Argonaute (Ago) proteins are ubiquitous across all kingdoms of life. Eukaryotic Agos (eAgos) use small RNAs to recognize transcripts for RNA silencing in eukaryotes. In contrast, the functions of prokaryotic counterparts (pAgo) are less well known. Recently, short pAgos in conjunction with the associated TIR or Sir2 (SPARTA or SPARSA) were found to serve as antiviral systems to combat phage infections. Herein, we present the cryo-EM structures of nicotinamide adenine dinucleotide (NAD+)-bound SPARSA with and without nucleic acids at resolutions of 3.1 Å and 3.6 Å, respectively. Our results reveal that the APAZ (Analogue of PAZ) domain and the short pAgo form a featured architecture similar to the long pAgo to accommodate nucleic acids. We further identified the key residues for NAD+ binding and elucidated the structural basis for guide RNA and target DNA recognition. Using structural comparisons, molecular dynamics simulations, and biochemical experiments, we proposed a putative mechanism for NAD+ hydrolysis in which an H186 loop mediates nucleophilic attack by catalytic water molecules. Overall, our study provides mechanistic insight into the antiphage role of the SPARSA system.