Open AccessE2EDNA: Simulation Protocol for DNA Aptamers with Ligands
Author(s) -
Michael Kilgour,
Tao Liu,
Brandon Walker,
Pengyu Ren,
Lena Simine
Publication year - 2021
Publication title -
journal of chemical information and modeling
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 1.24
H-Index - 160
eISSN - 1549-960X
pISSN - 1549-9596
DOI - 10.1021/acs.jcim.1c00696
Subject(s) - aptamer , computer science , flexibility (engineering) , dna , systematic evolution of ligands by exponential enrichment , molecular dynamics , folding (dsp implementation) , protocol (science) , computational biology , biological system , chemistry , engineering , biology , rna , computational chemistry , mathematics , biochemistry , alternative medicine , pathology , gene , medicine , statistics , genetics , electrical engineering
We present E2EDNA, a simulation protocol and accompanying code for the molecular biophysics and materials science communities. This protocol is both easy to use and sufficiently efficient to simulate single-stranded (ss)DNA and small analyte systems that are important to cellular processes and nanotechnologies such as DNA aptamer-based sensors. Existing computational tools used for aptamer design focus on cost-effective secondary structure prediction and motif analysis in the large data sets produced by SELEX experiments. As a rule, they do not offer flexibility with respect to the choice of the theoretical engine or direct access to the simulation platform. Practical aptamer optimization often requires higher accuracy predictions for only a small subset of sequences suggested, e.g., by SELEX experiments, but in the absence of a streamlined procedure, this task is extremely time and expertise intensive. We address this gap by introducing E2EDNA, a computational framework that accepts a DNA sequence in the FASTA format and the structures of the desired ligands and performs approximate folding followed by a refining step, analyte complexation, and molecular dynamics sampling at the desired level of accuracy. As a case study, we simulate a DNA-UTP (uridine triphosphate) complex in water using the state-of-the-art AMOEBA polarizable force field. The code is available at https://github.com/InfluenceFunctional/E2EDNA.