Entropy and Mg2+ control ligand affinity and specificity in the malachite green binding RNA aptamer

Effect of montmorillonite K10 clay on RNA structure and function

One of the earliest living systems was likely based on RNA ("the RNA world"). Mineral surfaces have been postulated to be an important environment for the prebiotic chemistry of RNA. In addition to adsorbing RNA and thus potentially reducing the chance of parasitic takeover through limited diffusion, minerals have been shown to promote a range of processes related to the emergence of life, including RNA polymerization, peptide bond formation, and self-assembly of vesicles. In addition, self-cleaving ribozymes have been shown to retain activity when adsorbed to the clay mineral montmorillonite. However, simulation studies suggest that adsorption to minerals is likely to interfere with RNA folding and, thus, function. To further evaluate the plausibility of a mineral-adsorbed RNA world, here we studied the effect of the synthetic clay montmorillonite K10 on the malachite green RNA aptamer, including binding of the clay to malachite green and RNA, as well as on the formation of secondary structures in model RNA and DNA oligonucleotides. We evaluated the fluorescence of the aptamer complex, adsorption to the mineral, melting curves, Förster resonance energy transfer interactions, and 1H-NMR signals to study the folding and functionality of these nucleic acids. Our results indicate that while some base pairings are unperturbed, the overall folding and binding of the malachite green aptamer are substantially disrupted by montmorillonite. These findings suggest that minerals would constrain the structures, and possibly the functions, available to an adsorbed RNA world.

Docking-aided rational tailoring of a fluorescence- and affinity-enhancing aptamer for a label-free ratiometric malachite green point-of-care aptasensor

Although nucleic acid aptasensors are increasingly applied in the detection of environmentally hazardous biomolecules, several formidable challenges remain with this technique because of their vulnerability, high cost and suboptimal sensitivity. Here, a docking-aided rational tailoring (DART) strategy was established at three levels and in two dimensions for the refinement of malachite green (MG) DNA aptamers. Guided by in silico molecular docking, coarse and fine tailoring were conducted at three levels each, to significantly enhance fluorescence activation intensity and binding affinity in two dimensions. Empowered by the results of the rational tailoring, a mechanistic view of the MG DNA aptamer-target interaction was thoroughly analyzed via four types of interactions. To meet the demand for point-of-care testing (POCT), a label-free and ratiometric fluorescent aptasensor was developed leveraging the tailored MG aptamer, based on the binding site competition-equilibrium effect via the introduction of a reference dye. This sensitive, specific, low-cost and rapid aptasensor subsequently demonstrated outstanding detection performance, achieving an ideal signal response range of 5 nmol·L^-1 - 6 μmol·L^-1 and a low limit of detection (LOD) of 1.49 nmol·L^-1. The DART strategy and systematic exploration of the MG DNA luminescent aptamers herein will provide a valuable reference in the field of aptamer tailoring, biosensing and bioimaging. The proposed label-free ratiometric aptasensor also provides a highly generalizable strategy for hazardous biomolecular detection.

In silico Druggability Assessment of the NUDIX Hydrolase Protein Family as a Workflow for Target Prioritization

Computational chemistry has now been widely accepted as a useful tool for shortening lead times in early drug discovery. When selecting new potential drug targets, it is important to assess the likelihood of finding suitable starting points for lead generation before pursuing costly high-throughput screening campaigns. By exploiting available high-resolution crystal structures, an in silico druggability assessment can facilitate the decision of whether, and in cases where several protein family members exist, which of these to pursue experimentally. Many of the algorithms and software suites commonly applied for in silico druggability assessment are complex, technically challenging and not always user-friendly. Here we applied the intuitive open access servers of DoGSite, FTMap and CryptoSite to comprehensively predict ligand binding pockets, druggability scores and conformationally active regions of the NUDIX protein family. In parallel we analyzed potential ligand binding sites, their druggability and pocket parameter using Schrödinger's SiteMap. Then an in silico docking cascade of a subset of the ZINC FragNow library using the Glide docking program was performed to assess identified pockets for large-scale small-molecule binding. Subsequently, this initial dual ranking of druggable sites within the NUDIX protein family was benchmarked against experimental hit rates obtained both in-house and by others from traditional biochemical and fragment screening campaigns. The observed correlation suggests that the presented user-friendly workflow of a dual parallel in silico druggability assessment is applicable as a standalone method for decision on target prioritization and exclusion in future screening campaigns.

Investigations on the interface of nucleic acid aptamers and binding targets

Nucleic acid aptamers are single-stranded DNA or RNA of 20-100 nucleotides in length that have attracted substantial scientific interest due to their ability to specifically bind to target molecules via the formation of three-dimensional structures. Compared to traditional protein antibodies, aptamers have several advantages, such as their small size, high binding affinity, specificity, flexible structure, being chemical synthesizable and modifiable, good biocompatibility, high stability and low immunogenicity, which all contribute to their widely applications in the biomedical field. To date, much progress has been made in the study and applications of aptamers, however, detailed information on how aptamers bind to their targets is still scarce. Over the past few decades, many methods have been introduced to investigate the aptamer-target binding process, such as measuring the main kinetic or thermodynamic parameters, detecting the structural changes of the binding complexes, etc. Apart from traditional physicochemical methods, various types of molecular docking programs have been applied to simulate the aptamer-target interactions, while these simulations also have limitations. To facilitate the further research on the interactions, herein, we provide a brief review to illustrate the recent advances in the study of aptamer-target interactions. We summarize the binding targets of aptamers, such as small molecules, macromolecules, and even cells. Their binding constants (KD) are also summarized. Methods to probe the aptamer-target binding process, such as surface plasmon resonance (SPR), circular dichroism spectroscopy (CD), isothermal titration calorimetry (ITC), footprinting assay, truncation and mutation assay, nuclear magnetic resonance spectroscopy (NMR), X-ray crystallography and molecular docking simulation are indicated. The binding forces mediating the aptamer-target interactions, such as hydrogen bonding, electrostatic interaction, the hydrophobic effect, π-π stacking and van der Waals forces are summarized. The challenges and future perspectives are also discussed.

New insights into a classic aptamer: binding sites, cooperativity and more sensitive adenosine detection

The DNA aptamer for adenosine (also for AMP and ATP) is a highly conserved sequence that has recurred in a few selections. It it a widely used model aptamer for biosensor development, and its nuclear magnetic resonance structure shows that each aptamer binds two AMP molecules. In this work, each binding site was individually removed by rational sequence design, while the remaining site still retained a similar binding affinity and specificity as confirmed by isothermal titration calorimetry. The thermodynamic parameters of binding are presented, and its biochemical implications are discussed. The number of binding sites can also be increased, and up to four sites are introduced in a single DNA sequence. Finally, the different sequences are made into fluorescent biosensors based on the structure-switching signaling aptamer design. The one-site aptamer has 3.8-fold higher sensitivity at lower adenosine concentration with a limit of detection of 9.1 μM adenosine, but weaker fluorescence signal at higher adenosine concentrations, consistent with a moderate cooperativity in the original aptamer. This work has offered insights into a classic aptamer for the relationship between the number of binding sites and sensitivity, and a shorter aptamer for improved biosensor design.

Organic additives stabilize RNA aptamer binding of malachite green

Aptamer-ligand binding has been utilized for biological applications due to its specific binding and synthetic nature. However, the applications will be limited if the binding or the ligand is unstable. Malachite green aptamer (MGA) and its labile ligand malachite green (MG) were found to have increasing apparent dissociation constants (Kd) as determined through the first order rate loss of emission intensity of the MGA-MG fluorescent complex. The fluorescent intensity loss was hypothesized to be from the hydrolysis of MG into malachite green carbinol base (MGOH). Random screening organic additives were found to reduce or retain the fluorescence emission and the calculated apparent Kd of MGA-MG binding. The protective effect became more apparent as the percentage of organic additives increased up to 10% v/v. The mechanism behind the organic additive protective effects was primarily from a ~5X increase in first order rate kinetics of MGOH→MG (kMGOH→MG), which significantly changed the equilibrium constant (Keq), favoring the generation of MG, versus MGOH without organic additives. A simple way has been developed to stabilize the apparent Kd of MGA-MG binding over 24h, which may be beneficial in stabilizing other triphenylmethane or carbocation ligand-aptamer interactions that are susceptible to SN1 hydrolysis.

Molecular Dynamics Simulation Analysis of Anti-MUC1 Aptamer and Mucin 1 Peptide Binding

Aptasensors utilize aptamers as bioreceptors. Aptamers are highly efficient, have a high specificity and are reusable. Within the biosensor the aptamers are immobilized to maximize their access to target molecules. Knowledge of the orientation and location of the aptamer and peptide during binding could be gained through computational modeling. Experimentally, the aptamer (anti-MUC1 S2.2) has been identified as a bioreceptor for breast cancer biomarker mucin 1 (MUC1) protein. However, within this protein lie several peptide variants with the common sequence APDTRPAP that are targeted by the aptamer. Understanding orientation and location of the binding region for a peptide-aptamer complex is critical in their biosensor applicability. In this study, we investigate through computational modeling how this peptide sequence and its minor variants affect the peptide-aptamer complex binding. We use molecular dynamics simulations to study multiple peptide-aptamer systems consisting of MUC1 (APDTRPAP) and MUC1-G (APDTRPAPG) peptides with the anti-MUC1 aptamer under similar physiological conditions reported experimentally. Multiple simulations of the MUC1 peptide and aptamer reveal that the peptide interacts between 3' and 5' ends of the aptamer but does not fully bind. Multiple simulations of the MUC1-G peptide indicate consistent binding with the thymine loop of the aptamer, initiated by the arginine residue of the peptide. We find that the binding event induces structural changes in the aptamer by altering the number of hydrogen bonds within the aptamer and establishes a stable peptide-aptamer complex. In all MUC1-G cases the occurrence of binding was confirmed by systematically studying the distance distributions between peptide and aptamers. These results are found to corroborate well with experimental study reported in the literature that indicated a strong binding in the case of MUC1-G peptide and anti-MUC1 aptamer. Present MD simulations highlight the role of the arginine residue of MUC1-G peptide in initiating the binding. The addition of the glycine residue to the peptide, as in the case of MUC1-G, is shown to yield a stable binding. Our study clearly demonstrates the ability of MD simulations to obtain molecular insights for peptide-aptamer binding, and to provide details on the orientation and location of binding between the peptide-aptamer that can be instrumental in biosensor development.

Malachite green interacts with the membrane skeletal protein, spectrin

Energy minimized complex of MG with the self association domain of spectrin.

Thermodynamics and kinetics of adaptive binding in the malachite green RNA aptamer

Adaptive binding, the ability of molecules to fold themselves around the structure of a ligand and thereby incorporating it into their three-dimensional fold, is a key feature of most RNA aptamers. The malachite green aptamer (MGA) has been shown to bind several closely related triphenyl dyes with planar and nonplanar structures in this manner. Competitive binding studies using isothermal titration calorimetry and stopped flow kinetics have been conducted with the aim of understanding the adaptive nature of RNA-ligand interaction. The results of these studies reveal that binding of one ligand can reduce the ability of the aptamer pocket to adapt to another ligand, even if this second ligand has a significantly higher affinity to the free aptamer. A similar effect is observed in the presence of Mg(2+) ions which stabilize the binding pocket in a more ligand bound-like conformation.

Potential toxicity and affinity of triphenylmethane dye malachite green to lysozyme

Malachite green is a triphenylmethane dye that is used extensively in many industrial and aquacultural processes, generating environmental concerns and health problems to human being. In this contribution, the complexation between lysozyme and malachite green was verified by means of computer-aided molecular modeling, steady state and time-resolved fluorescence, and circular dichroism (CD) approaches. The precise binding patch of malachite green in lysozyme has been identified from molecular modeling and ANS displacement, Trp-62, Trp-63, and Trp-108 residues of lysozyme were earmarked to possess high-affinity for this dye, the principal forces in the lysozyme-malachite green adduct are hydrophobic and π-π interactions. Steady state fluorescence proclaimed the complex of malachite green with lysozyme yields quenching through static type, which substantiates time-resolved fluorescence measurements that lysozyme-malachite green conjugation formation has an affinity of 10(3)M(-1). Moreover, via molecular modeling and also CD data, we can safely arrive at a conclusion that the polypeptide chain of lysozyme partially destabilized upon complexation with malachite green. The data emerged here will help to further understand the toxicological action of malachite green in human body.

Thermodynamics of ligand binding to a heterogeneous RNA population in the malachite green aptamer

The malachite green aptamer binds two closely related ligands, malachite green (MG) and tetramethylrosamine (TMR), with nearly equal affinity. The MG ligand consists of three phenyl rings emanating from a central carbon, while TMR has two of the three rings connected by an ether linkage. The binding pockets for MG and TMR in the aptamer, known from high-resolution structures, differ only in the conformation of a few nucleotides. Herein, we applied isothermal titration calorimetry (ITC) to compare the thermodynamics of binding of MG and TMR to the aptamer. Binding heat capacities were obtained from ITC titrations over the temperature range of 15-60 °C. Two temperature regimes were found for MG binding: one from 15 to 45 °C where MG bound with a large negative heat capacity and an apparent stoichiometry (n) of ~0.4 and another from 50 to 60 °C where MG bound with a positive heat capacity and an n of ~1.1. The binding of TMR, on the other hand, revealed only one temperature regime for binding, with a more modest negative heat capacity and an n of ~1.2. The large difference in heat capacity between the two ligands suggests that significantly more conformational rearrangement occurs upon the binding of MG than that of TMR, which is consistent with differences in solvent accessible surface area calculated for available ligand-bound structures. Lastly, we note that the binding stoichiometry of MG was improved not only by an increase in the temperature but also by a decrease in the concentration of Mg(2+) or an increase in the time between ITC injections. These studies suggest that binding of a dynamical ligand to a functional RNA requires the RNA itself to have significant dynamics.