Refinement of the Sugar-Phosphate Backbone Torsion Beta for AMBER Force Fields Improves the Description of Z- and B-DNA

Exploring the interaction of biologically active compounds with DNA through the application of the SwitchSense technique, UV-Vis spectroscopy, and computational methods

DNA is a key target for anticancer and antimicrobial drugs. Assessing the bioactivity of compounds involves in silico and instrumental studies to determine their affinity for biomolecules like DNA. This study explores the potential of the switchSense technique in rapidly evaluating compound bioactivity towards DNA. By combining switchSense with computational methods and UV-Vis spectrophotometry, various bioactive compounds' interactions with DNA were analyzed. The objects of the study were: netropsin (as a model compound that binds in the helical groove), as well as derivatives of pyrazine (PTCA), sulfonamide (NbutylS), and anthraquinone (AQ-NetOH). Though no direct correlation was found between switchSense kinetics and binding modes, this research suggests the technique's broader utility in assessing new compounds' interactions with DNA. used as analytes whose interactions with DNA have not been yet fully described in the literature.

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.

HNF4A and HNF1A exhibit tissue specific target gene regulation in pancreatic beta cells and hepatocytes

HNF4A and HNF1A encode transcription factors that are important for the development and function of the pancreas and liver. Mutations in both genes have been directly linked to Maturity Onset Diabetes of the Young (MODY) and type 2 diabetes (T2D) risk. To better define the pleiotropic gene regulatory roles of HNF4A and HNF1A, we generated a comprehensive genome-wide map of their binding targets in pancreatic and hepatic cells using ChIP-Seq. HNF4A was found to bind and regulate known (ACY3, HAAO, HNF1A, MAP3K11) and previously unidentified (ABCD3, CDKN2AIP, USH1C, VIL1) loci in a tissue-dependent manner. Functional follow-up highlighted a potential role for HAAO and USH1C as regulators of beta cell function. Unlike the loss-of-function HNF4A/MODY1 variant I271fs, the T2D-associated HNF4A variant (rs1800961) was found to activate AKAP1, GAD2 and HOPX gene expression, potentially due to changes in DNA-binding affinity. We also found HNF1A to bind to and regulate GPR39 expression in beta cells. Overall, our studies provide a rich resource for uncovering downstream molecular targets of HNF4A and HNF1A that may contribute to beta cell or hepatic cell (dys)function, and set up a framework for gene discovery and functional validation.

Dynamics of terminal fraying-peeling and hydrogen bonds dictates the sequential vs. cooperative melting pathways of nanoscale DNA and PNA triplexes

Peptide nucleic acids (PNAs) are charge-neutral synthetic DNA/RNA analogues. In many aspects of biology and biotechnology, the details of DNA and PNA melting reaction coordinates are crucial, and their associative/dissociative details remain inadequately understood. In the current study, we have attempted to gain insights into comparative melting pathways and binding affinity of iso-sequences of an 18-mer PNA-DNA-PNA triplex and the analogous DNA-DNA-DNA triplex, and DNA-DNA and PNA-DNA duplexes. It is intriguing that while the DNA-DNA-DNA triplex melts in two sequential steps, the PNA-DNA-PNA triplex melts in a single step and the mechanistic aspects for this difference are still not clear. We report an all-atom molecular dynamics simulation of both complexes in the temperature range of 300 to 500 K with 20 K intervals. Based on the trajectory analysis, we provide evidence that the association and dissociation are dictated by the differences in fraying-peeling effects from either terminus to the center in a zipper pattern among the PNA-DNA-PNA triplex and DNA-DNA-DNA triplexes. These are shown to be governed by the different characteristics of H-bonding, RMSD, and Free Energy Landscape (FEL) as analyzed by PCA, leading to the DNA-DNA-DNA triplex exhibiting sequential melting, while the PNA-DNA-PNA triplex shows cooperative melting of the whole fragment in a single-step. The PNA-DNA-PNA triplex base pairs are thermodynamically more stable than the DNA-DNA-DNA triplex, with the binding affinity of PNA-TFO to the PNA : DNA duplex being higher than that of DNA-TFO to the DNA : DNA duplex. The investigation of the association/dissociation of PNA-TFO to the PNA-DNA duplex has relevance and importance in the emerging effective applications of oligonucleotide therapy.

Modeling allosteric mechanisms of eukaryotic type II topoisomerases

Type II topoisomerases (TopoIIs) are ubiquitous enzymes that are involved in crucial nuclear processes such as genome organization, chromosome segregation, and other DNA metabolic processes. These enzymes function as large, homodimeric complexes that undergo a complex cycle of binding and hydrolysis of two ATP molecules in their ATPase domains, which regulates the capture and passage of one DNA double-helix through a second, cleaved DNA molecule. This process requires the transmission of information about the state of the bound nucleotide over vast ranges in the TopoII complex. How this information is transmitted at the molecular level to regulate TopoII functions and how protein substitutions disrupt these mechanisms remains largely unknown. Here, we employed extensive microsecond-scale molecular dynamics simulations of the yeast TopoII enzyme in multiple nucleotide-bound states and with amino acid substitutions near both the N and C termini of the complex. Simulation results indicate that the ATPase domains are remarkably flexible on the sub-microsecond timescale and that these dynamics are modulated by the identity of the bound nucleotides and both local and distant amino acid substitutions. Network analyses point toward specific allosteric networks that transmit information about the hydrolysis cycle throughout the complex, which include residues in both the protein and the bound DNA molecule. Amino acid substitutions weaken many of these pathways. Together, our results provide molecular level details on how the TopoII catalytic cycle is controlled through nucleotide binding and hydrolysis and how mutations may disrupt this process.

Capture and Storage of Cell-Free DNA via Bio-Informational Hydrogel Microspheres

Excessive cell-free DNA (cfDNA) can induce chronic inflammation by activating intracellular nucleic acid sensors. Intervention in cfDNA-mediated "pro-inflammatory signaling transduction" could be a potential alleviating strategy for chronic inflammation, such as in diabetic wounds. However, effectively and specifically downgrading cfDNA concentration in the pathological microenvironment remains a challenge. Therefore, this work prepares free-standing polydopamine nanosheets through DNA-guided assembly and loaded them into microfluidic hydrogel microspheres. The π─π stacking/hydrogen bonding interactions between polydopamine nanosheets and the π-rich bases of cfDNA, along with the cage-like spatial confinement created by the hydrogel polymer network, achieved cfDNA capture and storage, respectively. Catechol in polydopamine nanosheets can also assist in reducing reactive oxygen species (ROS) levels. Efficient cfDNA binding independent of serum proteins, specific interdiction of abnormal activation of cfDNA-associated toll-like receptor 9, as well as down-regulation of inflammatory cytokines and ROS levels are shown in this system. The chronic inflammation alleviating and the pro-healing effects on the mice model with diabetic wounds are also investigated. This work presents a new strategy for capturing and storing cfDNA to intervene in cell signaling transduction. It also offers new insights into the regulatory mechanisms between inflammatory mediators and biomaterials in inflammation-related diseases.

Transient interactions modulate the affinity of NF-κB transcription factors for DNA

The dimeric nuclear factor kappa B (NF-κB) transcription factors (TFs) regulate gene expression by binding to a variety of κB DNA elements with conserved G:C-rich flanking sequences enclosing a degenerate central region. Toward defining mechanistic principles of affinity regulated by degeneracy, we observed an unusual dependence of the affinity of RelA on the identity of the central base pair, which appears to be noncontacted in the complex crystal structures. The affinity of κB sites with A or T at the central position is ~10-fold higher than with G or C. The crystal structures of neither the complexes nor the free κB DNAs could explain the differences in affinity. Interestingly, differential dynamics of several residues were revealed in molecular dynamics simulation studies, where simulation replicates totaling 148 μs were performed on NF-κB:DNA complexes and free κB DNAs. Notably, Arg187 and Arg124 exhibited selectivity in transient interactions that orchestrated a complex interplay among several DNA-interacting residues in the central region. Binding and simulation studies with mutants supported these observations of transient interactions dictating specificity. In combination with published reports, this work provides insights into the nuanced mechanisms governing the discriminatory binding of NF-κB family TFs to κB DNA elements and sheds light on cancer pathogenesis of cRel, a close homolog of RelA.

CGeNArate: a sequence-dependent coarse-grained model of DNA for accurate atomistic MD simulations of kb-long duplexes

We present CGeNArate, a new model for molecular dynamics simulations of very long segments of B-DNA in the context of biotechnological or chromatin studies. The developed method uses a coarse-grained Hamiltonian with trajectories that are back-mapped to the atomistic resolution level with extreme accuracy by means of Machine Learning Approaches. The method is sequence-dependent and reproduces very well not only local, but also global physical properties of DNA. The efficiency of the method allows us to recover with a reduced computational effort high-quality atomic-resolution ensembles of segments containing many kilobases of DNA, entering into the gene range or even the entire DNA of certain cellular organelles.

ANNaMo: Coarse-grained modeling for folding and assembly of RNA and DNA systems

The folding of RNA and DNA strands plays crucial roles in biological systems and bionanotechnology. However, studying these processes with high-resolution numerical models is beyond current computational capabilities due to the timescales and system sizes involved. In this article, we present a new coarse-grained model for investigating the folding dynamics of nucleic acids. Our model represents three nucleotides with a patchy particle and is parameterized using well-established nearest-neighbor models. Thanks to the reduction of degrees of freedom and to a bond-swapping mechanism, our model allows for simulations at timescales and length scales that are currently inaccessible to more detailed models. To validate the performance of our model, we conducted extensive simulations of various systems: We examined the thermodynamics of DNA hairpins, capturing their stability and structural transitions, the folding of an MMTV pseudoknot, which is a complex RNA structure involved in viral replication, and also explored the folding of an RNA tile containing a k-type pseudoknot. Finally, we evaluated the performance of the new model in reproducing the melting temperatures of oligomers and the dependence on the toehold length of the displacement rate in toehold-mediated displacement processes, a key reaction used in molecular computing. All in all, the successful reproduction of experimental data and favorable comparisons with existing coarse-grained models validate the effectiveness of the new model.

Effects of the Y432S Cancer-Associated Variant on the Reaction Mechanism of Human DNA Polymerase κ

Human DNA polymerases are vital for genetic information management. Their function involves catalyzing the synthesis of DNA strands with unparalleled accuracy, which ensures the fidelity and stability of the human genomic blueprint. Several disease-associated mutations and their functional impact on DNA polymerases have been reported. One particular polymerase, human DNA polymerase kappa (Pol κ), has been reported to be susceptible to several cancer-associated mutations. The Y432S mutation in Pol κ, associated with various cancers, is of interest due to its impact on polymerization activity and markedly reduced thermal stability. Here, we have used computational simulations to investigate the functional consequences of the Y432S using classical molecular dynamics (MD) and coupled quantum mechanics/molecular mechanics (QM/MM) methods. Our findings suggest that Y432S induces structural alterations in domains responsible for nucleotide addition and ternary complex stabilization while retaining structural features consistent with possible catalysis in the active site. Calculations of the minimum energy path associated with the reaction mechanism of the wild type (WT) and Y432S Pol κ indicate that, while both enzymes are catalytically competent (in terms of energetics and the active site's geometries), the cancer mutation results in an endoergic reaction and an increase in the catalytic barrier. Interactions with a third magnesium ion and environmental effects on nonbonded interactions, particularly involving key residues, contribute to the kinetic and thermodynamic distinctions between the WT and mutant during the catalytic reaction. The energetics and electronic findings suggest that active site residues favor the catalytic reaction with dCTP3- over dCTP4-.

A pH-Responsive Topological Switch Based on a DNA Quadruplex-Duplex Hybrid

A guanine-rich oligonucleotide based on a human telomeric sequence but with the first three-nucleotide intervening stretch replaced by a putative 15-nucleotide hairpin-forming sequence shows a pH-dependent folding into different quadruplex-duplex hybrids in a potassium containing buffer. At slightly acidic pH, the quadruplex domain adopts a chair-type conformation. Upon increasing the pH, a transition with a midpoint close to neutral pH to a major and minor (3+1) hybrid topology with either a coaxially stacked or orthogonally oriented duplex stem-loop occurs. NMR-derived high-resolution structures reveal that an adenine protonation is prerequisite for the formation of a non-canonical base quartet, capping the outer G-tetrad at the quadruplex-duplex interface and stabilizing the antiparallel chair conformation in an acidic environment. Being directly associated with interactions at the quadruplex-duplex interface, this unique pH-dependent topological transition is fully reversible. Coupled with a conformation-sensitive optical readout demonstrated as a proof of concept using the fluorescent dye thiazole orange, the present quadruplex-duplex hybrid architecture represents a potentially valuable pH-sensing system responsive in a physiological pH range of 7±1.

Mechanical Stability and Unfolding Pathways of Parallel Tetrameric G-Quadruplexes Probed by Pulling Simulations

Guanine quadruplex (GQ) is a noncanonical nucleic acid structure formed by guanine-rich DNA and RNA sequences. Folding of GQs is a complex process, where several aspects remain elusive, despite being important for understanding structure formation and biological functions of GQs. Pulling experiments are a common tool for acquiring insights into the folding landscape of GQs. Herein, we applied a computational pulling strategy─steered molecular dynamics (SMD) simulations─in combination with standard molecular dynamics (MD) simulations to explore the unfolding landscapes of tetrameric parallel GQs. We identified anisotropic properties of elastic conformational changes, unfolding transitions, and GQ mechanical stabilities. Using a special set of structural parameters, we found that the vertical component of pulling force (perpendicular to the average G-quartet plane) plays a significant role in disrupting GQ structures and weakening their mechanical stabilities. We demonstrated that the magnitude of the vertical force component depends on the pulling anchor positions and the number of G-quartets. Typical unfolding transitions for tetrameric parallel GQs involve base unzipping, opening of the G-stem, strand slippage, and rotation to cross-like structures. The unzipping was detected as the first and dominant unfolding event, and it usually started at the 3'-end. Furthermore, results from both SMD and standard MD simulations indicate that partial spiral conformations serve as a transient ensemble during the (un)folding of GQs.

SNP-Associated Substitutions of Amino Acid Residues in the dNTP Selection Subdomain Decrease Polβ Polymerase Activity

In the cell, DNA polymerase β (Polβ) is involved in many processes aimed at maintaining genome stability and is considered the main repair DNA polymerase participating in base excision repair (BER). Polβ can fill DNA gaps formed by other DNA repair enzymes. Single-nucleotide polymorphisms (SNPs) in the POLB gene can affect the enzymatic properties of the resulting protein, owing to possible amino acid substitutions. For many SNP-associated Polβ variants, an association with cancer, owing to changes in polymerase activity and fidelity, has been shown. In this work, kinetic analyses and molecular dynamics simulations were used to examine the activity of naturally occurring polymorphic variants G274R, G290C, and R333W. Previously, the amino acid substitutions at these positions have been found in various types of tumors, implying a specific role of Gly-274, Gly-290, and Arg-333 in Polβ functioning. All three polymorphic variants had reduced polymerase activity. Two substitutions-G274R and R333W-led to the almost complete disappearance of gap-filling and primer elongation activities, a decrease in the deoxynucleotide triphosphate-binding ability, and a lower polymerization constant, due to alterations of local contacts near the replaced amino acid residues. Thus, variants G274R, G290C, and R333W may be implicated in an elevated level of unrepaired DNA damage.

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.

Dynamic Profiling and Binding Affinity Prediction of NBTI Antibacterials against DNA Gyrase Enzyme by Multidimensional Machine Learning and Molecular Dynamics Simulations

Bacterial type II topoisomerases are well-characterized and clinically important targets for antibacterial chemotherapy. Novel bacterial topoisomerase inhibitors (NBTIs) are a newly disclosed class of antibacterials. Prediction of their binding affinity to these enzymes would be beneficial for de novo design/optimization of new NBTIs. Utilizing in vitro NBTI experimental data, we constructed two comprehensive multidimensional DNA gyrase surrogate models for Staphylococcus aureus (q2 = 0.791) and Escherichia coli (q2 = 0.806). Both models accurately predicted the IC50s of 26 NBTIs from our recent studies. To investigate the NBTI's dynamic profile and binding to both targets, 10 selected NBTIs underwent molecular dynamics (MD) simulations. The analysis of MD production trajectories confirmed key hydrogen-bonding and hydrophobic contacts that NBTIs establish in both enzymes. Moreover, the binding free energies of selected NBTIs were computed by the linear interaction energy (LIE) method employing an in-house derived set of fitting parameters (α = 0.16, β = 0.029, γ = 0.0, and intercept = -1.72), which are successfully applicable to DNA gyrase of Gram-positive/Gram-negative pathogens. Both methods offer accurate predictions of the binding free energies of NBTIs against S. aureus and E. coli DNA gyrase. We are confident that this integrated modeling approach could be valuable in the de novo design and optimization of efficient NBTIs for combating resistant bacterial pathogens.

Structural insights into the morpholino nucleic acid/RNA duplex using the new XNA builder Ducque in a molecular modeling pipeline

The field of synthetic nucleic acids with novel backbone structures [xenobiotic nucleic acids (XNAs)] has flourished due to the increased importance of XNA antisense oligonucleotides and aptamers in medicine, as well as the development of XNA processing enzymes and new XNA genetic materials. Molecular modeling on XNA structures can accelerate rational design in the field of XNAs as it contributes in understanding and predicting how changes in the sugar-phosphate backbone impact on the complementation properties of the nucleic acids. To support the development of novel XNA polymers, we present a first-in-class open-source program (Ducque) to build duplexes of nucleic acid analogs with customizable chemistry. A detailed procedure is described to extend the Ducque library with new user-defined XNA fragments using quantum mechanics (QM) and to generate QM-based force field parameters for molecular dynamics simulations within standard packages such as AMBER. The tool was used within a molecular modeling workflow to accurately reproduce a selection of experimental structures for nucleic acid duplexes with ribose-based as well as non-ribose-based nucleosides. Additionally, it was challenged to build duplexes of morpholino nucleic acids bound to complementary RNA sequences.

The Impact of SNP-Induced Amino Acid Substitutions L19P and G66R in the dRP-Lyase Domain of Human DNA Polymerase β on Enzyme Activities

Base excision repair (BER), which involves the sequential activity of DNA glycosylases, apurinic/apyrimidinic endonucleases, DNA polymerases, and DNA ligases, is one of the enzymatic systems that preserve the integrity of the genome. Normal BER is effective, but due to single-nucleotide polymorphisms (SNPs), the enzymes themselves-whose main function is to identify and eliminate damaged bases-can undergo amino acid changes. One of the enzymes in BER is DNA polymerase β (Polβ), whose function is to fill gaps in DNA. SNPs can significantly affect the catalytic activity of an enzyme by causing an amino acid substitution. In this work, pre-steady-state kinetic analyses and molecular dynamics simulations were used to examine the activity of naturally occurring variants of Polβ that have the substitutions L19P and G66R in the dRP-lyase domain. Despite the substantial distance between the dRP-lyase domain and the nucleotidyltransferase active site, it was found that the capacity to form a complex with DNA and with an incoming dNTP is significantly altered by these substitutions. Therefore, the lower activity of the tested polymorphic variants may be associated with a greater number of unrepaired DNA lesions.

Sequence Dependence in Nucleosome Dynamics

The basic packaging unit of eukaryotic chromatin is the nucleosome that contains 145-147 base pair duplex DNA wrapped around an octameric histone protein. While the DNA sequence plays a crucial role in controlling the positioning of the nucleosome, the molecular details behind the interplay between DNA sequence and nucleosome dynamics remain relatively unexplored. This study analyzes this interplay in detail by performing all-atom molecular dynamics simulations of nucleosomes, comparing the human α-satellite palindromic (ASP) and the strong positioning "Widom-601" DNA sequence at time scales of 12 μs. The simulations are performed at salt concentrations 10-20 times higher than physiological salt concentrations to screen the electrostatic interactions and promote unwrapping. These microsecond-long simulations give insight into the molecular-level sequence-dependent events that dictate the pathway of DNA unwrapping. We find that the "ASP" sequence forms a loop around SHL ± 5 for three sets of simulations. Coincident with loop formation is a cooperative increase in contacts with the neighboring N-terminal H2B tail and C-terminal H2A tail and the release of neighboring counterions. We find that the Widom-601 sequence exhibits a strong breathing motion of the nucleic acid ends. Coincident with the breathing motion is the collapse of the full N-terminal H3 tail and formation of an α-helix that interacts with the H3 histone core. We postulate that the dynamics of these histone tails and their modification with post-translational modifications (PTMs) may play a key role in governing this dynamics.

Atomistic simulations of RNA duplex thermal denaturation: Sequence- and forcefield-dependence

Double-stranded RNA is the end-product of template-based replication, and is also the functional state of some biological RNAs. Similarly to proteins and DNA, they can be denatured by temperature, with important physiological and technological implications. Here, we use an in silico strategy to probe the thermal denaturation of RNA duplexes. Following previous results that were obtained on a few different duplexes, and which nuanced the canonical 2-state picture of nucleic acid denaturation, we here specifically address three different aspects that greatly improve our description of the temperature-induced dsRNA separation. First, we investigate the effect of the spatial distribution of weak and strong base-pairs among the duplex sequence. We show that the deviations from the two-state dehybridization mechanism are more pronounced when a strong core is flanked with weak extremities, while duplexes with a weak core but strong extremities exhibit a two-state behavior, which can be explained by the key role played by base fraying. This was later verified by generating artificial hairpin or circular states containing one or two locked duplex extremities, which results in an important reinforcement of the entire HB structure of the duplex and higher melting temperatures. Finally, we demonstrate that our results are little sensitive to the employed combination of RNA and water forcefields. The trends in thermal stability among the different sequences as well as the observed unfolding mechanisms (and the deviations from a two-state scenario) remain the same regardless of the employed atomistic models. However, our study points to possible limitations of recent reparametrizations of the Amber RNA forcefield, which sometimes results in duplexes that readily denature under ambient conditions, in contradiction with available experimental results.

Modeling Allosteric Mechanisms of Eukaryotic Type II Topoisomerases

Type II topoisomerases (TopoIIs) are ubiquitous enzymes that are involved in crucial nuclear processes such as genome organization, chromosome segregation, and other DNA metabolic processes. These enzymes function as large, homodimeric complexes that undergo a complex cycle of binding and hydrolysis of two ATP molecules in their ATPase domains, which regulates the capture and passage of one DNA double-helix through a second, cleaved DNA molecule. This process requires the transmission of information about the state of the bound nucleotide over vast ranges in the TopoII complex. How this information is transmitted at the molecular level to regulate TopoII functions and how protein substitutions disrupt these mechanisms remains largely unknown. Here, we employed extensive microsecond scale molecular dynamics simulations of the yeast TopoII enzyme in multiple nucleotide-bound states and with amino acid substitutions near both the N- and C-terminals of the complex. Simulation results indicate that the ATPase domains are remarkably flexible on the sub-microsecond timescale and that these dynamics are modulated by the identity of the bound nucleotides and both local and distant amino acid substitutions. Network analyses point towards specific allosteric networks that transmit information about the hydrolysis cycle throughout the complex, which include residues in both the protein and the bound DNA molecule. Amino acid substitutions weaken many of these pathways. Together, our results provide molecular-level details on how the TopoII catalytic cycle is controlled through nucleotide binding and hydrolysis and how mutations may disrupt this process.

Assessment of Amino Acid Electrostatic Parametrizations of the Polarizable Gaussian Multipole Model

Accurate parametrization of amino acids is pivotal for the development of reliable force fields for molecular modeling of biomolecules such as proteins. This study aims to assess amino acid electrostatic parametrizations with the polarizable Gaussian Multipole (pGM) model by evaluating the performance of the pGM-perm (with atomic permanent dipoles) and pGM-ind (without atomic permanent dipoles) variants compared to the traditional RESP model. The 100-conf-combterm fitting strategy on tetrapeptides was adopted, in which (1) all peptide bond atoms (-CO-NH-) share identical set of parameters and (2) the total charges of the two terminal N-acetyl (ACE) and N-methylamide (NME) groups were set to neutral. The accuracy and transferability of electrostatic parameters across peptides with varying lengths and real-world examples were examined. The results demonstrate the enhanced performance of the pGM-perm model in accurately representing the electrostatic properties of amino acids. This insight underscores the potential of the pGM-perm model and the 100-conf-combterm strategy for the future development of the pGM force field.

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.

Surface-Accelerated String Method for Locating Minimum Free Energy Paths

We present a surface-accelerated string method (SASM) to efficiently optimize low-dimensional reaction pathways from the sampling performed with expensive quantum mechanical/molecular mechanical (QM/MM) Hamiltonians. The SASM accelerates the convergence of the path using the aggregate sampling obtained from the current and previous string iterations, whereas approaches like the string method in collective variables (SMCV) or the modified string method in collective variables (MSMCV) update the path only from the sampling obtained from the current iteration. Furthermore, the SASM decouples the number of images used to perform sampling from the number of synthetic images used to represent the path. The path is optimized on the current best estimate of the free energy surface obtained from all available sampling, and the proposed set of new simulations is not restricted to being located along the optimized path. Instead, the umbrella potential placement is chosen to extend the range of the free energy surface and improve the quality of the free energy estimates near the path. In this manner, the SASM is shown to improve the exploration for a minimum free energy pathway in regions where the free energy surface is relatively flat. Furthermore, it improves the quality of the free energy profile when the string is discretized with too few images. We compare the SASM, SMCV, and MSMCV using 3 QM/MM applications: a ribozyme methyltransferase reaction using 2 reaction coordinates, the 2'-O-transphosphorylation reaction of Hammerhead ribozyme using 3 reaction coordinates, and a tautomeric reaction in B-DNA using 5 reaction coordinates. We show that SASM converges the paths using roughly 3 times less sampling than the SMCV and MSMCV methods. All three algorithms have been implemented in the FE-ToolKit package made freely available.

A comprehensive evolutionary scenario for the origin and neofunctionalization of the Drosophila speciation gene Odysseus (OdsH)

Odysseus (OdsH) was the first speciation gene described in Drosophila related to hybrid sterility in offspring of mating between Drosophila mauritiana and Drosophila simulans. Its origin is attributed to the duplication of the gene unc-4 in the subgenus Sophophora. By using a much larger sample of Drosophilidae species, we showed that contrary to what has been previously proposed, OdsH origin occurred 62 MYA. Evolutionary rates, expression, and transcription factor-binding sites of OdsH evidence that it may have rapidly experienced neofunctionalization in male sexual functions. Furthermore, the analysis of the OdsH peptide allowed the identification of mutations of D. mauritiana that could result in incompatibility in hybrids. In order to find if OdsH could be related to hybrid sterility, beyond Sophophora, we explored the expression of OdsH in Drosophila arizonae and Drosophila mojavensis, a pair of sister species with incomplete reproductive isolation. Our data indicated that OdsH expression is not atypical in their male-sterile hybrids. In conclusion, we have proposed that the origin of OdsH occurred earlier than previously proposed, followed by neofunctionalization. Our results also suggested that its role as a speciation gene might be restricted to D. mauritiana and D. simulans.

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 stability of D5SIC-DNAM-incorporated DNA duplex in Taq polymerase binary system: a systematic classical MD approach

DNA polymerases are fundamental enzymes that play a crucial role in processing DNA with high fidelity and accuracy ensuring the faithful transmission of genetic information. The recognition of unnatural base pairs (UBPs) by polymerases, enabling their replication, represents a significant and groundbreaking discovery with profound implications for genetic expansion. Romesberg et al. examined the impact of DNA containing 2,6-dimethyl-2H-isoquiniline-1-thione: D5SIC (DS) and 2-methoxy-3-methylnaphthalene: DNAM (DN) UBPs bound to T. aquaticus DNA polymerase (Taq) through crystal structure analysis. Here, we have used polarizable and nonpolarizable classical molecular dynamics (MD) simulations to investigate the structural aspects and stability of Taq in complex with a DNA duplex including a DS-DN pair in the terminal 3' and 5' positions. Our results suggest that the flexibility of UBP-incorporated DNA in the terminal position is arrested by the polymerase, thus preventing fraying and mispairing. Our investigation also reveals that the UBP remains in an intercalated conformation inside the active site, exhibiting two distinct orientations in agreement with experimental findings. Our analysis pinpoints particular residues responsible for favorable interactions with the UBP, with some relying on van der Waals interactions while other on Coulombic forces.

Investigation of dynamical flexibility of D5SIC-DNAM inside DNA duplex in aqueous solution: a systematic classical MD approach

Incorporation of artificial 3rd base pairs (unnatural base pairs, UBPs) has emerged as a fundamental technique in pursuit of expanding the genetic alphabet. 2,6-Dimethyl-2H-isoquiniline-1-thione: D5SIC (DS) and 2-methoxy-3-methylnaphthalene: DNAM (DN), a potential unnatural base pair (UBP) developed by Romesberg and colleagues, has been shown to have remarkable capability for replication within DNA. Crystal structures of a Taq polymerase/double-stranded DNA (ds-DNA) complex containing a DS-DN pair in the 3' terminus showed a parallelly stacked geometry for the pre-insertion, and an intercalated geometry for the post-insertion structure. Unconventional orientations of DS-DN inside a DNA duplex have inspired scientists to investigate the conformational orientations and structural properties of UBP-incorporated DNA. In recent years, computational simulations have been used to investigate the geometry of DS-DN within the DNA duplex; nevertheless, unresolved questions persist owing to inconclusive findings. In this work, we investigate the structural and dynamical properties of DS and DN inside a ds-DNA strand in aqueous solution considering both short and long DNA templates using polarizable, and non-polarizable classical MD simulations. Flexible conformational change of UBP with major populations of Watson-Crick-Franklin (WCF) and three distinct non-Watson-Crick-Franklin (nWCFP1, nWCFP2, nWCFO) conformations through intra and inter-strand flipping have been observed. Our results suggest that a dynamical conformational change leads to the production of diffierent conformational distribution for the systems. Simulations with a short ds-DNA duplex suggest nWCF (P1 and O) as the predominant structures, whereas long ds-DNA duplex simulations indicate almost equal populations of WCF, nWCFP1, nWCFO. DS-DN in the terminal position is found to be more flexible with occasional mispairing and fraying. Overall, these results suggest flexibility and dynamical conformational change of the UBP as well as indicate varied conformational distribution irrespective of starting orientation of the UBP and length og DNA strand.

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.

Structure-activity relationships for the G-quadruplex-targeting experimental drug QN-302 and two analogues probed with comparative transcriptome profiling and molecular modeling

The tetrasubstituted naphthalene diimide compound QN-302 binds to G-quadruplex (G4) DNA structures. It shows high potency in pancreatic ductal adenocarcinoma (PDAC) cells and inhibits the transcription of cancer-related genes in these cells and in PDAC animal models. It is currently in Phase 1a clinical evaluation as an anticancer drug. A study of structure-activity relationships of QN-302 and two related analogues (CM03 and SOP1247) is reported here. These have been probed using comparisons of transcriptional profiles from whole-genome RNA-seq analyses, together with molecular modelling and molecular dynamics simulations. Compounds CM03 and SOP1247 differ by the presence of a methoxy substituent in the latter: these two compounds have closely similar transcriptional profiles. Whereas QN-302 (with an additional benzyl-pyrrolidine group), although also showing down-regulatory effects in the same cancer-related pathways, has effects on distinct genes, for example in the hedgehog pathway. This distinctive pattern of genes affected by QN-302 is hypothesized to contribute to its superior potency compared to CM03 and SOP1247. Its enhanced ability to stabilize G4 structures has been attributed to its benzyl-pyrrolidine substituent fitting into and filling most of the space in a G4 groove compared to the hydrogen atom in CM03 or the methoxy group substituent in SOP1247.

Are kuravirus capsid diameters quantized? The first all-atom genome tracing method for double-stranded DNA viruses

The revolution in cryo-electron microscopy has resulted in unprecedented power to resolve large macromolecular complexes including viruses. Many methods exist to explain density corresponding to proteins and thus entire protein capsids have been solved at the all-atom level. However methods for nucleic acids lag behind, and no all-atom viral double-stranded DNA genomes have been published at all. We here present a method which exploits the spiral winding patterns of DNA in icosahedral capsids. The method quickly generates shells of DNA wound in user-specified, idealized spherical or cylindrical spirals. For transition regions, the method allows guided semiflexible fitting. For the kuravirus SU10, our method explains most of the density in a semiautomated fashion. The results suggest rules for DNA turns in the end caps under which two discrete parameters determine the capsid inner diameter. We suggest that other kuraviruses viruses may follow the same winding scheme, producing a discrete rather than continuous spectrum of capsid inner diameters. Our software may be used to explain the published density maps of other double-stranded DNA viruses and uncover their genome packaging principles.

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.

Enhanced molecular dynamic simulation studies unravel long-range effects caused by sequence variations and partner binding in RNA aptamers

Intrinsic flexibility and structural modularity are two common features of RNA molecules. Although functionally crucial, RNA plasticity often represents a major complication in high-resolution structural studies. To overcome this problem, RNAs may be rigidified through the complexation with high-affinity partners such as Fab molecules. This approach has been previously used to characterize the DIR2-aptamer. However, possible perturbations induced by the insertion of the Fab binding site on the DIR2-aptamer conformational properties were not investigated. Here, using enhanced molecular dynamics simulations, we compared the dynamics of the DIR2 aptamer holding the Fab binding site with that of the parental sequence. Our results suggest that the L2-loop modification for the Fab recognition leads to a significant increase in local flexibility that also affects the mobility of distant regions. The trajectories provide clear indications of the groups and the interactions mediating the dynamics transfer in DIR2. The effectiveness of our approach in addressing RNA flexibility was further corroborated by showing its ability to reproduce the most important events affecting the NF-κB RNA aptamer upon dissociation from the partner. Therefore, REMD analyses, a rarely adopted technique to unravel the structural/dynamical properties of aptamers, could efficiently complement experimental data guiding the rational design of nucleic acid therapeutics.

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.

Click and Detect: Versatile Ampicillin Aptasensor Enabled by Click Chemistry on a Graphene-Alkyne Derivative

Tackling the current problem of antimicrobial resistance (AMR) requires fast, inexpensive, and effective methods for controlling and detecting antibiotics in diverse samples at the point of interest. Cost-effective, disposable, point-of-care electrochemical biosensors are a particularly attractive option. However, there is a need for conductive and versatile carbon-based materials and inks that enable effective bioconjugation under mild conditions for the development of robust, sensitive, and selective devices. This work describes a simple and fast methodology to construct an aptasensor based on a novel graphene derivative equipped with alkyne groups prepared via fluorographene chemistry. Using click chemistry, an aptamer is immobilized and used as a successful platform for the selective determination of ampicillin in real samples in the presence of interfering molecules. The electrochemical aptasensor displayed a detection limit of 1.36 nM, high selectivity among other antibiotics, the storage stability of 4 weeks, and is effective in real samples. Additionally, structural and docking simulations of the aptamer shed light on the ampicillin binding mechanism. The versatility of this platform opens up wide possibilities for constructing a new class of aptasensor based on disposable screen-printed carbon electrodes usable in point-of-care devices.

Simple Adjustment of Intranucleotide Base-Phosphate Interaction in the OL3 AMBER Force Field Improves RNA Simulations

Molecular dynamics (MD) simulations represent an established tool to study RNA molecules. The outcome of MD studies depends, however, on the quality of the force field (ff). Here we suggest a correction for the widely used AMBER OL3 ff by adding a simple adjustment of the nonbonded parameters. The reparameterization of the Lennard-Jones potential for the -H8···O5'- and -H6···O5'- atom pairs addresses an intranucleotide steric clash occurring in the type 0 base-phosphate interaction (0BPh). The nonbonded fix (NBfix) modification of 0BPh interactions (NBfix0BPh modification) was tuned via a reweighting approach and subsequently tested using an extensive set of standard and enhanced sampling simulations of both unstructured and folded RNA motifs. The modification corrects minor but visible intranucleotide clash for the anti nucleobase conformation. We observed that structural ensembles of small RNA benchmark motifs simulated with the NBfix0BPh modification provide better agreement with experiments. No side effects of the modification were observed in standard simulations of larger structured RNA motifs. We suggest that the combination of OL3 RNA ff and NBfix0BPh modification is a viable option to improve RNA MD simulations.

Understanding ATP Binding to DosS Catalytic Domain with a Short ATP-Lid

DosS is a heme-containing histidine kinase that triggers dormancy transformation inMycobacterium tuberculosis. Sequence comparison of the catalytic ATP-binding (CA) domain of DosS to other well-studied histidine kinases reveals a short ATP-lid. This feature has been thought to block binding of ATP to DosS's CA domain in the absence of interactions with DosS's dimerization and histidine phospho-transfer (DHp) domain. Here, we use a combination of computational modeling, structural biology, and biophysical studies to re-examine ATP-binding modalities in DosS. We show that the closed-lid conformation observed in crystal structures of DosS CA is caused by the presence of Zn2+ in the ATP binding pocket that coordinates with Glu537 on the ATP-lid. Furthermore, circular dichroism studies and comparisons of DosS CA's crystal structure with its AlphaFold model and homologous DesK reveal that residues 503-507 that appear as a random coil in the Zn2+-coordinated crystal structure are in fact part of the N-box α helix needed for efficient ATP binding. Such random-coil transformation of an N-box α helix turn and the closed-lid conformation are both artifacts arising from large millimolar Zn2+ concentrations used in DosS CA crystallization buffers. In contrast, in the absence of Zn2+, the short ATP-lid of DosS CA has significant conformational flexibility and can effectively bind AMP-PNP (Kd = 53 ± 13 μM), a non-hydrolyzable ATP analog. Furthermore, the nucleotide affinity remains unchanged when CA is conjugated to the DHp domain (Kd = 51 ± 6 μM). In all, our findings reveal that the short ATP-lid of DosS CA does not hinder ATP binding and provide insights that extend to 2988 homologous bacterial proteins containing such ATP-lids.

AutoParams: An Automated Web-Based Tool To Generate Force Field Parameters for Molecular Dynamics Simulations

Many research questions benefit from molecular dynamics simulations to observe the motions and conformations of molecules over time, which rely on force fields that describe sets of common molecules by category. With the increase of importance for large data sets used in machine learning and growing computational efficiency, the ability to rapidly create large numbers of force field inputs is of high importance. Unusual molecules, such as nucleotide analogues, functionalized carbohydrates, and modified amino acids, are difficult to describe consistently using standard force fields, requiring the development of custom parameters for each unique molecule. While these parameters may be created by individual users, the process can become time-consuming or may introduce errors that may not be immediately apparent. We present an open-source automated parameter generation service, AutoParams, which requires minimal input from the user and creates useful Amber force field parameter sets for most molecules, particularly those that combine molecular types (e.g., a carbohydrate functionalized with a benzene). We include hierarchical atom-typing logic that makes it straightforward to expand with additional force fields and settings, and options for creating monomers in polymers, such as functionalized amino acids. It can be straightforwardly linked to any charge generation program and currently has interfaces to Psi4, PsiRESP, and TeraChem. It is open source and is available via GitHub. It includes error checking and testing protocols to ensure the parameters will be sufficient for subsequent molecular dynamics simulations and streamlines the creation of force field databases.

Transcriptional repression by a secondary DNA binding surface of DNA topoisomerase I safeguards against hypertranscription

Regulation of global transcription output is important for normal development and disease, but little is known about the mechanisms involved. DNA topoisomerase I (TOP1) is an enzyme well-known for its role in relieving DNA supercoils for enabling transcription. Here, we report a non-enzymatic function of TOP1 that downregulates RNA synthesis. This function is dependent on specific DNA-interacting residues located on a conserved protein surface. A loss-of-function knock-in mutation on this surface, R548Q, is sufficient to cause hypertranscription and alter differentiation outcomes in mouse embryonic stem cells (mESCs). Hypertranscription in mESCs is accompanied by reduced TOP1 chromatin binding and change in genomic supercoiling. Notably, the mutation does not impact TOP1 enzymatic activity; rather, it diminishes TOP1-DNA binding and formation of compact protein-DNA structures. Thus, TOP1 exhibits opposing influences on transcription through distinct activities which are likely to be coordinated. This highlights TOP1 as a safeguard of appropriate total transcription levels in cells.

Showcasing Different G-Quadruplex Folds of a G-Rich Sequence: Between Rule-Based Prediction and Butterfly Effect

In better understanding the interactions of G-quadruplexes in a cellular or noncellular environment, a reliable sequence-based prediction of their three-dimensional fold would be extremely useful, yet is often limited by their remarkable structural diversity. A G-rich sequence related to a promoter sequence of the PDGFR-β nuclease hypersensitivity element (NHE) comprises a G3-G3-G2-G4-G3 pattern of five G-runs with two to four G residues. Although the predominant formation of three-layered canonical G-quadruplexes with uninterrupted G-columns can be expected, minimal base substitutions in a non-G-tract domain were shown to guide folding into either a basket-type antiparallel quadruplex, a parallel-stranded quadruplex with an interrupted G-column, a quadruplex with a V-shaped loop, or a (3+1) hybrid quadruplex. A 3D NMR structure for each of the different folds was determined. Supported by thermodynamic profiling on additional sequence variants, formed topologies were rationalized by the identification and assessment of specific critical interactions of loop and overhang residues, giving valuable insights into their contribution to favor a particular conformer. The variability of such tertiary interactions, together with only small differences in quadruplex free energies, emphasizes current limits for a reliable sequence-dependent prediction of favored topologies from sequences with multiple irregularly positioned G-tracts.

Mechanical codes of chemical-scale specificity in DNA motifs

In gene transcription, certain sequences of double-stranded (ds)DNA play a vital role in nucleosome positioning and expression initiation. That dsDNA is deformed to various extents in these processes leads us to ask: Could the genomic DNA also have sequence specificity in its chemical-scale mechanical properties? We approach this question using statistical machine learning to determine the rigidity between DNA chemical moieties. What emerges for the polyA, polyG, TpA, and CpG sequences studied here is a unique trigram that contains the quantitative mechanical strengths between bases and along the backbone. In a way, such a sequence-dependent trigram could be viewed as a DNA mechanical code. Interestingly, we discover a compensatory competition between the axial base-stacking interaction and the transverse base-pairing interaction, and such a reciprocal relationship constitutes the most discriminating feature of the mechanical code. Our results also provide chemical-scale understanding for experimental observables. For example, the long polyA persistence length is shown to have strong base stacking while its complement (polyAc) exhibits high backbone rigidity. The mechanical code concept enables a direct reading of the physical interactions encoded in the sequence which, with further development, is expected to shed new light on DNA allostery and DNA-binding drugs.

Streamlining and Optimizing Strategies of Electrostatic Parameterization

Accurate characterization of electrostatic interactions is crucial in molecular simulation. Various methods and programs have been developed to obtain electrostatic parameters for additive or polarizable models to replicate electrostatic properties obtained from experimental measurements or theoretical calculations. Electrostatic potentials (ESPs), a set of physically well-defined observables from quantum mechanical (QM) calculations, are well suited for optimization efforts due to the ease of collecting a large amount of conformation-dependent data. However, a reliable set of QM ESP computed at an appropriate level of theory and atomic basis set is necessary. In addition, despite the recent development of the PyRESP program for electrostatic parameterizations of induced dipole-polarizable models, the time-consuming and error-prone input file preparation process has limited the widespread use of these protocols. This work aims to comprehensively evaluate the quality of QM ESPs derived by eight methods, including wave function methods such as Hartree-Fock (HF), second-order Møller-Plesset (MP2), and coupled cluster-singles and doubles (CCSD), as well as five hybrid density functional theory (DFT) methods, used in conjunction with 13 different basis sets. The highest theory levels CCSD/aug-cc-pV5Z (a5z) and MP2/aug-cc-pV5Z (a5z) were selected as benchmark data over two homemade data sets. The results show that the hybrid DFT method, ωB97X-D, combined with the aug-cc-pVTZ (a3z) basis set, performs well in reproducing ESPs while taking both accuracy and efficiency into consideration. Moreover, a flexible and user-friendly program called PyRESP_GEN was developed to streamline input file preparation. The restraining strengths, along with strategies for polarizable Gaussian multipole (pGM) model parameterizations, were also optimized. These findings and the program presented in this work facilitate the development and application of induced dipole-polarizable models, such as pGM models, for molecular simulations of both chemical and biological significance.

Structure and Dynamics of DNA and RNA Double Helices Formed by d(CTG), d(GTC), r(CUG), and r(GUC) Trinucleotide Repeats and Associated DNA-RNA Hybrids

Myotonic dystrophy type 1 is the most frequent form of muscular dystrophy in adults caused by an abnormal expansion of the CTG trinucleotide. Both the expanded DNA and the expanded CUG RNA transcript can fold into hairpins. Co-transcriptional formation of stable RNA·DNA hybrids can also enhance the instability of repeat tracts. We performed molecular dynamics simulations of homoduplexes associated with the disease, d(CTG)n and r(CUG)n, and their corresponding r(CAG)n:d(CTG)n and r(CUG)n:d(CAG)n hybrids that can form under bidirectional transcription and of non-pathological d(GTC)n and d(GUC)n homoduplexes. We characterized their conformations, stability, and dynamics and found that the U·U and T·T mismatches are dynamic, favoring anti-anti conformations inside the helical core, followed by anti-syn and syn-syn conformations. For DNA, the secondary minima in the non-expanding d(GTC)n helices are deeper, wider, and longer-lived than those in d(CTG)n, which constitutes another biophysical factor further differentiating the expanding and non-expanding sequences. The hybrid helices are closer to A-RNA, with the A-T and A-U pairs forming two stable Watson-Crick hydrogen bonds. The neutralizing ion distribution around the non-canonical pairs is also described.

Terahertz waves regulate the mechanical unfolding of tau pre-mRNA hairpins

Intermolecular interactions, including hydrogen bonds, dominate the pairing and unpairing of nucleic acid chains in the transfer process of genetic information. The energy of THz waves just matches with the weak interactions, so THz waves may interact with biomolecules. Here, the dynamic effects of THz electromagnetic (EM) waves on the mechanical unfolding process of RNA hairpins (WT-30nt and its mutants, rHP, SARS-CoV-2, and SRV-1 SF206) are investigated using steered molecular dynamics (SMD) simulations. The results show that THz waves can either promote the unfolding of the double helix of the RNA hairpin during the initial unfolding phase (4-21.8 THz) or significantly enhance (23.8 and 25.5 THz) or weaken (37.4 and 41.2 THz) its structural stability during unfolding. Our findings have important implications for applying THz waves to regulate dynamic deconvolution processes, such as gene replication, transcription, and translation.

Charge Engineering of the Nucleic Acid Binding Cleft of a Thermostable HIV-1 Reverse Transcriptase Reveals Key Interactions and a Novel Mechanism of RNase H Inactivation

Coupled with PCR, reverse transcriptases (RTs) have been widely used for RNA detection and gene expression analysis. Increased thermostability and nucleic acid binding affinity are desirable RT properties to improve yields and sensitivity of these applications. The effects of amino acid substitutions in the RT RNase H domain were tested in an engineered HIV-1 group O RT, containing mutations K358R/A359G/S360A and devoid of RNase H activity due to the presence of E478Q (O3MQ RT). Twenty mutant RTs with Lys or Arg at positions interacting with the template-primer (i.e., at positions 473-477, 499-502 and 505) were obtained and characterized. Most of them produced significant amounts of cDNA at 37, 50 and 65 °C, as determined in RT-PCR reactions. However, a big loss of activity was observed with mutants A477K/R, S499K/R, V502K/R and Y505K/R, particularly at 65 °C. Binding affinity experiments confirmed that residues 477, 502 and 505 were less tolerant to mutations. Amino acid substitutions Q500K and Q500R produced a slight increase of cDNA synthesis efficiency at 50 and 65 °C, without altering the KD for model DNA/DNA and RNA/DNA heteroduplexes. Interestingly, molecular dynamics simulations predicted that those mutations inactivate the RNase H activity by altering the geometry of the catalytic site. Proof of this unexpected effect was obtained after introducing Q500K or Q500R in the wild-type HIV-1BH10 RT and mutant K358R/A359G/S360A RT. Our results reveal a novel mechanism of RNase H inactivation that preserves RT DNA binding and polymerization efficiency without substituting RNase H active site residues.

SELEX based aptamers with diagnostic and entry inhibitor therapeutic potential for SARS-CoV-2

Frequent mutation and variable immunological protection against vaccination is a common feature for COVID-19 pandemic. Early detection and confinement remain key to controlling further spread of infection. In response, we have developed an aptamer-based system that possesses both diagnostic and therapeutic potential towards the virus. A random aptamer library (~ 1017 molecules) was screened using systematic evolution of ligands by exponential enrichment (SELEX) and aptamer R was identified as a potent binder for the SARS-CoV-2 spike receptor binding domain (RBD) using in vitro binding assay. Using a pseudotyped viral entry assay we have shown that aptamer R specifically inhibited the entry of a SARS-CoV-2 pseudotyped virus in HEK293T-ACE2 cells but did not inhibit the entry of a Vesicular Stomatitis Virus (VSV) glycoprotein (G) pseudotyped virus, hence establishing its specificity towards SARS-CoV-2 spike protein. The antiviral potential of aptamers R and J (same central sequence as R but lacking flanked primer regions) was tested and showed 95.4% and 82.5% inhibition, respectively, against the SARS-CoV-2 virus. Finally, intermolecular interactions between the aptamers and the RBD domain were analyzed using in silico docking and molecular dynamics simulations that provided additional insight into the binding and inhibitory action of aptamers R and J.

Comprehensive analysis unveils altered binding kinetics of 5-/6-methylCytosine/adenine modifications in R2R3-DNA system

Recent studies have shown that DNA methylation is an important epigenetic marker. Two prominent forms are methylation of the C5 position of cytosine and methylation of the C6 position of adenine. Given the vital significance of DNA methylation, investigating the mechanisms that influence protein binding remains a compelling pursuit. This study used molecular dynamics simulations to investigate the binding patterns of R2R3 protein and four differentially methylated DNAs. The alanine scanning combined with interaction entropy method was used to identify key residues that respond to different methylation patterns. The order of protein binding ability to DNA is as follows: unmethylated DNA > A11 methylation (5'-A6mAC-3') (6m2A system) > A10 methylation (5'-6mAAC-3') (6m1A system) > both A10 and A11 methylation (5'-6mA6mAC-3') (6mAA system) > C12 methylation (5'-AA5mC-3') (5mC system). All methylation systems lead to the sixth α helix (H6) (residues D105 to L116) moving away from the binding interface, and in the 5mC and 6m1A systems, the third α helix (H3) (residues G54 to L65) exhibits a similar trend. When the positively charged amino acids in H3 and H6 move away from the binding interface, their electrostatic and van der Waals interactions with the negatively charged DNA are weakened. Structural changes induced by methylation contributed to the destabilization of the hydrogen bond network near the original binding site, except for the 6m2A system. Moreover, there is a positive correlation between the number of methylated sites and the probability of distorting the DNA structure. Our study explores how different methylation patterns affect binding and structural adaptability, and have implications for drug discovery and understanding diseases related to abnormal methylation.

Quantitative investigation of the effects of DNA modifications and protein mutations on MeCP2-MBD-DNA interactions

DNA methylation as an important epigenetic marker, has gained attention for the significance of three oxidative modifications (hydroxymethyl-C (hmC), formyl-C (fC), and carboxyl-C (caC)). Mutations occurring in the methyl-CpG-binding domain (MBD) of MeCP2 result in Rett. However, uncertainties persist regarding DNA modification and MBD mutation-induced interaction changes. Here, molecular dynamics simulations were used to investigate the underlying mechanisms behind changes due to different modifications of DNA and MBD mutations. Alanine scanning combined with the interaction entropy method was employed to accurately evaluate the binding free energy. The results show that MBD has the strongest binding ability for mCDNA, followed by caC, hmC, and fCDNA, with the weakest binding ability observed for CDNA. Further analysis revealed that mC modification induces DNA bending, causing residues R91 and R162 closer to the DNA. This proximity enhances van der Waals and electrostatic interactions. Conversely, the caC/hmC and fC modifications lead to two loop regions (near K112 and K130) closer to DNA, respectively. Furthermore, DNA modifications promote the formation of stable hydrogen bond networks, however mutations in the MBD significantly reduce the binding free energy. This study provides detailed insight into the effects of DNA modifications and MBD mutations on binding ability. It emphasizes the necessity for research and development of targeted Rett compounds that induce conformational compatibility between MBD and DNA, enhancing the stability and strength of their interactions.

Complexity of Guanine Quadruplex Unfolding Pathways Revealed by Atomistic Pulling Simulations

Guanine quadruplexes (GQs) are non-canonical nucleic acid structures involved in many biological processes. GQs formed in single-stranded regions often need to be unwound by cellular machinery, so their mechanochemical properties are important. Here, we performed steered molecular dynamics simulations of human telomeric GQs to study their unfolding. We examined four pulling regimes, including a very slow setup with pulling velocity and force load accessible to high-speed atomic force microscopy. We identified multiple factors affecting the unfolding mechanism, i.e.,: (i) the more the direction of force was perpendicular to the GQ channel axis (determined by GQ topology), the more the base unzipping mechanism happened, (ii) the more parallel the direction of force was, GQ opening and cross-like GQs were more likely to occur, (iii) strand slippage mechanism was possible for GQs with an all-anti pattern in a strand, and (iv) slower pulling velocity led to richer structural dynamics with sampling of more intermediates and partial refolding events. We also identified that a GQ may eventually unfold after a force drop under forces smaller than those that the GQ withstood before the drop. Finally, we found out that different unfolding intermediates could have very similar chain end-to-end distances, which reveals some limitations of structural interpretations of single-molecule spectroscopic data.