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.

Sensitivity of the RNA Structure to Ion Conditions as Probed by Molecular Dynamics Simulations of Common Canonical RNA Duplexes.

RNA molecules play a key role in countless biochemical processes. RNA interactions, which are of highly diverse nature, are determined by the fact that RNA is a highly negatively charged polyelectrolyte, which leads to intimate interactions with an ion atmosphere. Although RNA molecules are formally single-stranded, canonical (Watson-Crick) duplexes are key components of folded RNAs. A double-stranded (ds) RNA is also important for the design of RNA-based nanostructures and assemblies. Despite the fact that the description of canonical dsRNA is considered the least problematic part of RNA modeling, the imperfect shape and flexibility of dsRNA can lead to imbalances in the simulations of larger RNAs and RNA-containing assemblies. We present a comprehensive set of molecular dynamics (MD) simulations of four canonical A-RNA duplexes. Our focus was directed toward the characterization of the influence of varying ion concentrations and of the size of the solvation box. We compared several water models and four RNA force fields. The simulations showed that the A-RNA shape was most sensitive to the RNA force field, with some force fields leading to a reduced inclination of the A-RNA duplexes. The ions and water models played a minor role. The effect of the box size was negligible, and even boxes with a small fraction of the bulk solvent outside the RNA hydration sphere were sufficient for the simulation of the dsRNA.

Computer Aided Development of Nucleic Acid Applications in Nanotechnologies.

Utilization of nucleic acids (NAs) in nanotechnologies and nanotechnology-related applications is a growing field with broad application potential, ranging from biosensing up to targeted cell delivery. Computer simulations are useful techniques that can aid design and speed up development in this field. This review focuses on computer simulations of hybrid nanomaterials composed of NAs and other components. Current state-of-the-art molecular dynamics simulations, empirical force fields (FFs), and coarse-grained approaches for the description of deoxyribonucleic acid and ribonucleic acid are critically discussed. Challenges in combining biomacromolecular and nanomaterial FFs are emphasized. Recent applications of simulations for modeling NAs and their interactions with nano- and biomaterials are overviewed in the fields of sensing applications, targeted delivery, and NA templated materials. Future perspectives of development are also highlighted.

Parallel G-triplexes and G-hairpins as potential transitory ensembles in the folding of parallel-stranded DNA G-Quadruplexes.

Guanine quadruplexes (G4s) are non-canonical nucleic acids structures common in important genomic regions. Parallel-stranded G4 folds are the most abundant, but their folding mechanism is not fully understood. Recent research highlighted that G4 DNA molecules fold via kinetic partitioning mechanism dominated by competition amongst diverse long-living G4 folds. The role of other intermediate species such as parallel G-triplexes and G-hairpins in the folding process has been a matter of debate. Here, we use standard and enhanced-sampling molecular dynamics simulations (total length of ∼0.9 ms) to study these potential folding intermediates. We suggest that parallel G-triplex per se is rather an unstable species that is in local equilibrium with a broad ensemble of triplex-like structures. The equilibrium is shifted to well-structured G-triplex by stacked aromatic ligand and to a lesser extent by flanking duplexes or nucleotides. Next, we study propeller loop formation in GGGAGGGAGGG, GGGAGGG and GGGTTAGGG sequences. We identify multiple folding pathways from different unfolded and misfolded structures leading towards an ensemble of intermediates called cross-like structures (cross-hairpins), thus providing atomistic level of description of the single-molecule folding events. In summary, the parallel G-triplex is a possible, but not mandatory short-living (transitory) intermediate in the folding of parallel-stranded G4.

RNA Structural Dynamics As Captured by Molecular Simulations: A Comprehensive Overview.

With both catalytic and genetic functions, ribonucleic acid (RNA) is perhaps the most pluripotent chemical species in molecular biology, and its functions are intimately linked to its structure and dynamics. Computer simulations, and in particular atomistic molecular dynamics (MD), allow structural dynamics of biomolecular systems to be investigated with unprecedented temporal and spatial resolution. We here provide a comprehensive overview of the fast-developing field of MD simulations of RNA molecules. We begin with an in-depth, evaluatory coverage of the most fundamental methodological challenges that set the basis for the future development of the field, in particular, the current developments and inherent physical limitations of the atomistic force fields and the recent advances in a broad spectrum of enhanced sampling methods. We also survey the closely related field of coarse-grained modeling of RNA systems. After dealing with the methodological aspects, we provide an exhaustive overview of the available RNA simulation literature, ranging from studies of the smallest RNA oligonucleotides to investigations of the entire ribosome. Our review encompasses tetranucleotides, tetraloops, a number of small RNA motifs, A-helix RNA, kissing-loop complexes, the TAR RNA element, the decoding center and other important regions of the ribosome, as well as assorted others systems. Extended sections are devoted to RNA-ion interactions, ribozymes, riboswitches, and protein/RNA complexes. Our overview is written for as broad of an audience as possible, aiming to provide a much-needed interdisciplinary bridge between computation and experiment, together with a perspective on the future of the field.

Structural dynamics of propeller loop: towards folding of RNA G-quadruplex.

We have carried out an extended set of standard and enhanced-sampling MD simulations (for a cumulative simulation time of 620 µs) with the aim to study folding landscapes of the rGGGUUAGGG and rGGGAGGG parallel G-hairpins (PH) with propeller loop. We identify folding and unfolding pathways of the PH, which is bridged with the unfolded state via an ensemble of cross-like structures (CS) possessing mutually tilted or perpendicular G-strands interacting via guanine-guanine H-bonding. The oligonucleotides reach the PH conformation from the unfolded state via a conformational diffusion through the folding landscape, i.e. as a series of rearrangements of the H-bond interactions starting from compacted anti-parallel hairpin-like structures. Although isolated PHs do not appear to be thermodynamically stable we suggest that CS and PH-types of structures are sufficiently populated during RNA guanine quadruplex (GQ) folding within the context of complete GQ-forming sequences. These structures may participate in compact coil-like ensembles that involve all four G-strands and already some bound ions. Such ensembles can then rearrange into the fully folded parallel GQs via conformational diffusion. We propose that the basic atomistic folding mechanism of propeller loops suggested in this work may be common for their formation in RNA and DNA GQs.

Folding of guanine quadruplex molecules-funnel-like mechanism or kinetic partitioning? An overview from MD simulation studies.

BACKGROUND: Guanine quadruplexes (GQs) play vital roles in many cellular processes and are of much interest as drug targets. In contrast to the availability of many structural studies, there is still limited knowledge on GQ folding. SCOPE OF REVIEW: We review recent molecular dynamics (MD) simulation studies of the folding of GQs, with an emphasis paid to the human telomeric DNA GQ. We explain the basic principles and limitations of all types of MD methods used to study unfolding and folding in a way accessible to non-specialists. We discuss the potential role of G-hairpin, G-triplex and alternative GQ intermediates in the folding process. We argue that, in general, folding of GQs is fundamentally different from funneled folding of small fast-folding proteins, and can be best described by a kinetic partitioning (KP) mechanism. KP is a competition between at least two (but often many) well-separated and structurally different conformational ensembles. MAJOR CONCLUSIONS: The KP mechanism is the only plausible way to explain experiments reporting long time-scales of GQ folding and the existence of long-lived sub-states. A significant part of the natural partitioning of the free energy landscape of GQs comes from the ability of the GQ-forming sequences to populate a large number of syn-anti patterns in their G-tracts. The extreme complexity of the KP of GQs typically prevents an appropriate description of the folding landscape using just a few order parameters or collective variables. GENERAL SIGNIFICANCE: We reconcile available computational and experimental studies of GQ folding and formulate basic principles characterizing GQ folding landscapes. This article is part of a Special Issue entitled "G-quadruplex" Guest Editor: Dr. Concetta Giancola and Dr. Daniela Montesarchio.

Hairpins participating in folding of human telomeric sequence quadruplexes studied by standard and T-REMD simulations.

DNA G-hairpins are potential key structures participating in folding of human telomeric guanine quadruplexes (GQ). We examined their properties by standard MD simulations starting from the folded state and long T-REMD starting from the unfolded state, accumulating ∼130 µs of atomistic simulations. Antiparallel G-hairpins should spontaneously form in all stages of the folding to support lateral and diagonal loops, with sub-µs scale rearrangements between them. We found no clear predisposition for direct folding into specific GQ topologies with specific syn/anti patterns. Our key prediction stemming from the T-REMD is that an ideal unfolded ensemble of the full GQ sequence populates all 4096 syn/anti combinations of its four G-stretches. The simulations can propose idealized folding pathways but we explain that such few-state pathways may be misleading. In the context of the available experimental data, the simulations strongly suggest that the GQ folding could be best understood by the kinetic partitioning mechanism with a set of deep competing minima on the folding landscape, with only a small fraction of molecules directly folding to the native fold. The landscape should further include non-specific collapse processes where the molecules move via diffusion and consecutive random rare transitions, which could, e.g. structure the propeller loops.

Molecular dynamic simulations of protein/RNA complexes: CRISPR/Csy4 endoribonuclease.

BACKGROUND: Many prokaryotic genomes comprise Clustered Regularly Interspaced Short Palindromic Repeats (CRISPRs) offering defense against foreign nucleic acids. These immune systems are conditioned by the production of small CRISPR-derived RNAs matured from long RNA precursors. This often requires a Csy4 endoribonuclease cleaving the RNA 3'-end. METHODS: We report extended explicit solvent molecular dynamic (MD) simulations of Csy4/RNA complex in precursor and product states, based on X-ray structures of product and inactivated precursor (55 simulations; ~3.7µs in total). RESULTS: The simulations identify double-protonated His29 and deprotonated terminal phosphate as the likely dominant protonation states consistent with the product structure. We revealed potential substates consistent with Ser148 and His29 acting as the general base and acid, respectively. The Ser148 could be straightforwardly deprotonated through solvent and could without further structural rearrangements deprotonate the nucleophile, contrasting similar studies investigating the general base role of nucleobases in ribozymes. We could not locate geometries consistent with His29 acting as general base. However, we caution that the X-ray structures do not always capture the catalytically active geometries and then the reactive structures may be unreachable by the simulation technique. CONCLUSIONS: We identified potential catalytic arrangement of the Csy4/RNA complex but we also report limitations of the simulation technique. Even for the dominant protonation state we could not achieve full agreement between the simulations and the structural data. GENERAL SIGNIFICANCE: Potential catalytic arrangement of the Csy4/RNA complex is found. Further, we provide unique insights into limitations of simulations of protein/RNA complexes, namely, the influence of the starting experimental structures and force field limitations. This article is part of a Special Issue entitled Recent developments of molecular dynamics.

Exponential repulsion improves structural predictability of molecular docking.

Molecular docking is a powerful tool for theoretical prediction of the preferred conformation and orientation of small molecules within protein active sites. The obtained poses can be used for estimation of binding energies, which indicate the inhibition effect of designed inhibitors, and therefore might be used for in silico drug design. However, the evaluation of ligand binding affinity critically depends on successful prediction of the native binding mode. Contemporary docking methods are often based on scoring functions derived from molecular mechanical potentials. In such potentials, nonbonded interactions are typically represented by electrostatic interactions between atom-centered partial charges and standard 6-12 Lennard-Jones potential. Here, we present implementation and testing of a scoring function based on more physically justified exponential repulsion instead of the standard Lennard-Jones potential. We found that this scoring function significantly improved prediction of the native binding modes in proteins bearing narrow active sites such as serine proteases and kinases. © 2016 Wiley Periodicals, Inc.

Disparate HDV ribozyme crystal structures represent intermediates on a rugged free-energy landscape.

The hepatitis delta virus (HDV) ribozyme is a member of the class of small, self-cleaving catalytic RNAs found in a wide range of genomes from HDV to human. Both pre- and post-catalysis (precursor and product) crystal structures of the cis-acting genomic HDV ribozyme have been determined. These structures, together with extensive solution probing, have suggested that a significant conformational change accompanies catalysis. A recent crystal structure of a trans-acting precursor, obtained at low pH and by molecular replacement from the previous product conformation, conforms to the product, raising the possibility that it represents an activated conformer past the conformational change. Here, using fluorescence resonance energy transfer (FRET), we discovered that cleavage of this ribozyme at physiological pH is accompanied by a structural lengthening in magnitude comparable to previous trans-acting HDV ribozymes. Conformational heterogeneity observed by FRET in solution appears to have been removed upon crystallization. Analysis of a total of 1.8 µsec of molecular dynamics (MD) simulations showed that the crystallographically unresolved cleavage site conformation is likely correctly modeled after the hammerhead ribozyme, but that crystal contacts and the removal of several 2'-oxygens near the scissile phosphate compromise catalytic in-line fitness. A cis-acting version of the ribozyme exhibits a more dynamic active site, while a G-1 residue upstream of the scissile phosphate favors poor fitness, allowing us to rationalize corresponding changes in catalytic activity. Based on these data, we propose that the available crystal structures of the HDV ribozyme represent intermediates on an overall rugged RNA folding free-energy landscape.

Chemical feasibility of the general acid/base mechanism of glmS ribozyme self-cleavage.

In numerous Gram-positive bacteria, the glmS ribozyme or catalytic riboswitch regulates the expression of glucosamine-6-phosphate (GlcN6P) synthase via site-specific cleavage of its sugar-phosphate backbone in response to GlcN6P ligand binding. Biochemical data have suggested a crucial catalytic role for an active site guanine (G40 in Thermoanaerobacter tengcongensis, G33 in Bacillus anthracis). We used hybrid quantum chemical/molecular mechanical (QM/MM) calculations to probe the mechanism where G40 is deprotonated and acts as a general base. The calculations suggest that the deprotonated guanine G40(-) is sufficiently reactive to overcome the thermodynamic penalty arising from its rare protonation state, and thus is able to activate the A-1(2'-OH) group toward nucleophilic attack on the adjacent backbone. Furthermore, deprotonation of A-1(2'-OH) and nucleophilic attack are predicted to occur as separate steps, where activation of A-1(2'-OH) precedes nucleophilic attack. Conversely, the transition state associated with the rate-determining step corresponds to concurrent nucleophilic attack and protonation of the G1(O5') leaving group by the ammonium moiety of the GlcN6P cofactor. Overall, our calculations help to explain the crucial roles of G40 (as a general base) and GlcN6P (as a general acid) during glmS ribozyme self-cleavage. In addition, we show that the QM/MM description of the glmS ribozyme self-cleavage reaction is significantly more sensitive to the size of the QM region and the quality of the QM-MM coupling than that of other small ribozymes.

The role of an active site Mg(2+) in HDV ribozyme self-cleavage: insights from QM/MM calculations.

The hepatitis delta virus (HDV) ribozyme is a catalytic RNA motif embedded in the human pathogenic HDV RNA. It catalyzes self-cleavage of its sugar-phosphate backbone with direct participation of the active site cytosine C75. Biochemical and structural data support a general acid role of C75. Here, we used hybrid quantum mechanical/molecular mechanical (QM/MM) calculations to probe the reaction mechanism and changes in Gibbs energy along the ribozyme's reaction pathway with an N3-protonated C75H(+) in the active site, which acts as the general acid, and a partially hydrated Mg(2+) ion with one deprotonated, inner-shell coordinated water molecule that acts as the general base. We followed eight reaction paths with a distinct position and coordination of the catalytically important active site Mg(2+) ion. For six of them, we observed feasible activation barriers ranging from 14.2 to 21.9 kcal mol(-1), indicating that the specific position of the Mg(2+) ion in the active site is predicted to strongly affect the kinetics of self-cleavage. The deprotonation of the U-1(2'-OH) nucleophile and the nucleophilic attack of the resulting U-1(2'-O(-)) on the scissile phosphodiester are found to be separate steps, as deprotonation precedes the nucleophilic attack. This sequential mechanism of the HDV ribozyme differs from the concerted nucleophilic activation and attack suggested for the hairpin ribozyme. We estimate the pKa of the U-1(2'-OH) group to range from 8.8 to 11.2, suggesting that it is lowered by several units from that of a free ribose, comparable to and most likely smaller than the pKa of the solvated active site Mg(2+) ion. Our results thus support the notion that the structure of the HDV ribozyme, and particularly the positioning of the active site Mg(2+) ion, facilitate deprotonation and activation of the 2'-OH nucleophile.

Wobble pairs of the HDV ribozyme play specific roles in stabilization of active site dynamics.

The hepatitis delta virus (HDV) is the only known human pathogen whose genome contains a catalytic RNA motif (ribozyme). The overall architecture of the HDV ribozyme is that of a double-nested pseudoknot, with two GU pairs flanking the active site. Although extensive studies have shown that mutation of either wobble results in decreased catalytic activity, little work has focused on linking these mutations to specific structural effects on catalytic fitness. Here we use molecular dynamics simulations based on an activated structure to probe the active site dynamics as a result of wobble pair mutations. In both wild-type and mutant ribozymes, the in-line fitness of the active site (as a measure of catalytic proficiency) strongly depends on the presence of a C75(N3H3+)N1(O5') hydrogen bond, which positions C75 as the general acid for the reaction. Our mutational analyses show that each GU wobble supports catalytically fit conformations in distinct ways; the reverse G25U20 wobble promotes high in-line fitness, high occupancy of the C75(N3H3+)G1(O5') general-acid hydrogen bond and stabilization of the G1U37 wobble, while the G1U37 wobble acts more locally by stabilizing high in-line fitness and the C75(N3H3+)G1(O5') hydrogen bond. We also find that stable type I A-minor and P1.1 hydrogen bonding above and below the active site, respectively, prevent local structural disorder from spreading and disrupting global conformation. Taken together, our results define specific, often redundant architectural roles for several structural motifs of the HDV ribozyme active site, expanding the known roles of these motifs within all HDV-like ribozymes and other structured RNAs.

Anatomy of enzyme channels.

BACKGROUND: Enzyme active sites can be connected to the exterior environment by one or more channels passing through the protein. Despite our current knowledge of enzyme structure and function, surprisingly little is known about how often channels are present or about any structural features such channels may have in common. RESULTS: Here, we analyze the long channels (i.e. >15 Å) leading to the active sites of 4,306 enzyme structures. We find that over 64% of enzymes contain two or more long channels, their typical length being 28 Å. We show that amino acid compositions of the channel significantly differ both to the composition of the active site, surface and interior of the protein. CONCLUSIONS: The majority of enzymes have buried active sites accessible via a network of access channels. This indicates that enzymes tend to have buried active sites, with channels controlling access to, and egress from, them, and that suggests channels may play a key role in helping determine enzyme substrate.

How to understand quantum chemical computations on DNA and RNA systems? A practical guide for non-specialists.

In this review primarily written for non-experts we explain basic methodological aspects and interpretation of modern quantum chemical (QM) computations applied to nucleic acids. We introduce current reference QM computations on small model systems consisting of dozens of atoms. Then we comment on recent advance of fast and accurate dispersion-corrected density functional theory methods, which will allow computations of small but complete nucleic acids building blocks in the near future. The qualitative difference between QM and molecular mechanics (MM, force field) computations is discussed. We also explain relation of QM and molecular simulation computations to experiments.

MOLEonline 2.0: interactive web-based analysis of biomacromolecular channels.

Biomolecular channels play important roles in many biological systems, e.g. enzymes, ribosomes and ion channels. This article introduces a web-based interactive MOLEonline 2.0 application for the analysis of access/egress paths to interior molecular voids. MOLEonline 2.0 enables platform-independent, easy-to-use and interactive analyses of (bio)macromolecular channels, tunnels and pores. Results are presented in a clear manner, making their interpretation easy. For each channel, MOLEonline displays a 3D graphical representation of the channel, its profile accompanied by a list of lining residues and also its basic physicochemical properties. The users can tune advanced parameters when performing a channel search to direct the search according to their needs. The MOLEonline 2.0 application is freely available via the Internet at http://ncbr.muni.cz/mole or http://mole.upol.cz.

Nature and magnitude of aromatic base stacking in DNA and RNA: Quantum chemistry, molecular mechanics, and experiment.

Base stacking is a major interaction shaping up and stabilizing nucleic acids. During the last decades, base stacking has been extensively studied by experimental and theoretical methods. Advanced quantum-chemical calculations clarified that base stacking is a common interaction, which in the first approximation can be described as combination of the three most basic contributions to molecular interactions, namely, electrostatic interaction, London dispersion attraction and short-range repulsion. There is not any specific π-π energy term associated with the delocalized π electrons of the aromatic rings that cannot be described by the mentioned contributions. The base stacking can be rather reasonably approximated by simple molecular simulation methods based on well-calibrated common force fields although the force fields do not include nonadditivity of stacking, anisotropy of dispersion interactions, and some other effects. However, description of stacking association in condensed phase and understanding of the stacking role in biomolecules remain a difficult problem, as the net base stacking forces always act in a complex and context-specific environment. Moreover, the stacking forces are balanced with many other energy contributions. Differences in definition of stacking in experimental and theoretical studies are explained.

Benchmark quantum-chemical calculations on a complete set of rotameric families of the DNA sugar-phosphate backbone and their comparison with modern density functional theory.

The DNA sugar-phosphate backbone has a substantial influence on the DNA structural dynamics. Structural biology and bioinformatics studies revealed that the DNA backbone in experimental structures samples a wide range of distinct conformational substates, known as rotameric DNA backbone conformational families. Their correct description is essential for methods used to model nucleic acids and is known to be the Achilles heel of force field computations. In this study we report the benchmark database of MP2 calculations extrapolated to the complete basis set of atomic orbitals with aug-cc-pVTZ and aug-cc-pVQZ basis sets, MP2(T,Q), augmented by ΔCCSD(T)/aug-cc-pVDZ corrections. The calculations are performed in the gas phase as well as using a COSMO solvent model. This study includes a complete set of 18 established and biochemically most important families of DNA backbone conformations and several other salient conformations that we identified in experimental structures. We utilize an electronically sufficiently complete DNA sugar-phosphate-sugar (SPS) backbone model system truncated to prevent undesired intramolecular interactions. The calculations are then compared with other QM methods. The BLYP and TPSS functionals supplemented with Grimme's D3(BJ) dispersion term provide the best tradeoff between computational demands and accuracy and can be recommended for preliminary conformational searches as well as calculations on large model systems. Among the tested methods, the best agreement with the benchmark database has been obtained for the double-hybrid DSD-BLYP functional in combination with a quadruple-ζ basis set, which is, however, computationally very demanding. The new hybrid density functionals PW6B95-D3 and MPW1B95-D3 yield outstanding results and even slightly outperform the computationally more demanding PWPB95 double-hybrid functional. B3LYP-D3 is somewhat less accurate compared to the other hybrids. Extrapolated MP2(D,T) calculations are not as accurate as the less demanding DFT-D3 methods. Preliminary force field tests using several charge sets reveal an almost order of magnitude larger deviations from the reference QM data compared to modern DFT-D3, underlining the challenges facing force field simulations of nucleic acids. As expected, inclusion of the solvent environment approximated by a continuum approach has a large impact on the relative stabilities of different backbone substates and is important when comparing the QM data with structural bioinformatics and other experimental data.