Structurally Different Yet Functionally Similar: Aptamers Specific for the Ebola Virus Soluble Glycoprotein and GP1,2 and Their Application in Electrochemical Sensing.

The Ebola virus glycoprotein (GP) gene templates several mRNAs that produce either the virion-associated transmembrane protein or one of two secreted glycoproteins. Soluble glycoprotein (sGP) is the predominant product. GP1 and sGP share an amino terminal sequence of 295 amino acids but differ in quaternary structure, with GP1 being a heterohexamer with GP2 and sGP a homodimer. Two structurally different DNA aptamers were selected against sGP that also bound GP1,2. These DNA aptamers were compared with a 2'FY-RNA aptamer for their interactions with the Ebola GP gene products. The three aptamers have almost identical binding isotherms for sGP and GP1,2 in solution and on the virion. They demonstrated high affinity and selectivity for sGP and GP1,2. Furthermore, one aptamer, used as a sensing element in an electrochemical format, detected GP1,2 on pseudotyped virions and sGP with high sensitivity in the presence of serum, including from an Ebola-virus-infected monkey. Our results suggest that the aptamers interact with sGP across the interface between the monomers, which is different from the sites on the protein bound by most antibodies. The remarkable similarity in functional features of three structurally distinct aptamers suggests that aptamers, like antibodies, have preferred binding sites on proteins.

Computational Modeling of RNA Aptamers: Structure Prediction of the Apo State.

RNA aptamers are single-stranded oligonucleotides that bind to specific molecular targets with high affinity and specificity. To design aptamers for new applications, it is critical to understand the ligand binding mechanism in terms of the structure and dynamics of the ligand-bound and apo states. The problem is that most of the NMR or X-ray crystal structures available for RNA aptamers are for ligand-bound states. Available apo state structures, mostly characterized by crystallization under nonphysiological conditions or probed by low resolution techniques, might fail to represent the diverse structural variations of the apo state in solution. Here, we develop an approach to obtain a representative ensemble of apo structures that are based on in silico RNA 3D structure prediction and in vitro experiments that characterize base stacking. Using the neomycin-B aptamer as a case study, an ensemble of structures for the aptamer in the apo (unbound) state are validated and then used to investigate the ligand-binding mechanism for the aptamer in complex with neomycin-B.

Sampling Performance of Multiple Independent Molecular Dynamics Simulations of an RNA Aptamer.

Using multiple independent simulations instead of one long simulation has been shown to improve the sampling performance attained with the molecular dynamics (MD) simulation method. However, it is generally not known how long each independent simulation should be, how many independent simulations should be used, or to what extent either of these factors affects the overall sampling performance achieved for a given system. The goal of the present study was to assess the sampling performance of multiple independent MD simulations, where each independent simulation begins from a different initial molecular conformation. For this purpose, we used an RNA aptamer that is 25 nucleotides long as a case study. The initial conformations of the aptamer are derived from six de novo predicted 3D structures. Each of the six de novo predicted structures is energy minimized in solution and equilibrated with MD simulations at high temperature. Ten conformations from these six high-temperature equilibration runs are selected as initial conformations for further simulations at ambient temperature. In total, we conducted 60 independent MD simulations, each with a duration of 100 ns, to study the conformation and dynamics of the aptamer. For each group of 10 independent simulations that originated from a particular de novo predicted structure, we evaluated the potential energy distribution of the RNA and used recurrence quantification analysis to examine the sampling of RNA conformational transitions. To assess the impact of starting from different de novo predicted structures, we computed the density of structure projection on principal components to compare the regions sampled by the different groups of ten independent simulations. The recurrence rate and dependence of initial conformation among the groups were also compared. We stress the necessity of using different initial configurations as simulation starting points by showing long simulations from different initial structures suffer from being trapped in different states. Finally, we summarized the sampling efficiency for the complete set of 60 independent simulations and determined regions of under-sampling on the potential energy landscape. The results suggest that conducting multiple independent simulations using a diverse set of de novo predicted structures is a promising approach to achieve sufficient sampling. This approach avoids undesirable outcomes, such as the problem of the RNA aptamer being trapped in a local minimum. For others wishing to conduct multiple independent simulations, the analysis protocol presented in this study is a guide for examining overall sampling and determining if more simulations are necessary for sufficient sampling.

Common Secondary and Tertiary Structural Features of Aptamer-Ligand Interaction Shared by RNA Aptamers with Different Primary Sequences.

Aptamer selection can yield many oligonucleotides with different sequences and affinities for the target molecule. Here, we have combined computational and experimental approaches to understand if aptamers with different sequences but the same molecular target share structural and dynamical features. NEO1A, with a known NMR-solved structure, displays a flexible loop that interacts differently with individual aminoglycosides, its ligand affinities and specificities are responsive to ionic strength, and it possesses an adenosine in the loop that is critical for high-affinity ligand binding. NEO2A was obtained from the same selection and, although they are only 43% identical in overall sequence, NEO1A and NEO2A share similar loop sequences. Experimental analysis by 1D NMR and 2-aminopurine reporters combined with molecular dynamics modeling revealed similar structural and dynamical characteristics in both aptamers. These results are consistent with the hypothesis that the target ligand drives aptamer structure and also selects relevant dynamical characteristics for high-affinity aptamer-ligand interaction. Furthermore, they suggest that it might be possible to "migrate" structural and dynamical features between aptamer group members with different primary sequences but with the same target ligand.

Study strategies and beliefs about learning as a function of academic achievement and achievement goals.

Prior research by Hartwig and Dunlosky [(2012). Study strategies of college students: Are self-testing and scheduling related to achievement? Psychonomic Bulletin & Review, 19(1), 126-134] has demonstrated that beliefs about learning and study strategies endorsed by students are related to academic achievement: higher performing students tend to choose more effective study strategies and are more aware of the benefits of self-testing. We examined whether students' achievement goals, independent of academic achievement, predicted beliefs about learning and endorsement of study strategies. We administered Hartwig and Dunlosky's survey, along with the Achievement Goals Questionnaire [Elliot, A. J., & McGregor, H. A. (2001). A 2 × 2 achievement goal framework. Journal of Personality & Social Psychology, 80, 501-519] to a large undergraduate biology course. Similar to results by Hartwig and Dunlosky, we found that high-performing students (relative to low-performing students) were more likely to endorse self-testing, less likely to cram, and more likely to plan a study schedule ahead of time. Independent of achievement, however, achievement goals were stronger predictors of certain study behaviours. In particular, avoidance goals (e.g., fear of failure) coincided with increased use of cramming and the tendency to be driven by impending deadlines. Results suggest that individual differences in student achievement, as well as the underlying reasons for achievement, are important predictors of students' approaches to studying.

Prequestions do not enhance the benefits of retrieval in a STEM classroom.

Answering questions before a learning episode-"prequestions"-can enhance memory for that information. A number of studies have explored this effect in the laboratory; however, few studies have examined prequestions in a classroom setting. In the current study, the effects of prequestions were examined in an undergraduate course in chemical engineering. At the start of several class meetings, students were provided with a prequestion to answer about the upcoming lesson, and then were asked to provide ratings of confidence in their answers, familiarity with the content in the prequestion, and how much of the assigned reading they had completed. At the end of class, students were given the same question again (postquestion), along with a different question from the same lesson (new question). On a quiz at the end of each week, students were given the postquestions and new questions again, in addition to never-before-seen questions (quiz-only questions) from the same lessons. Performance on questions at the end of class revealed no difference in performance for postquestions vs. new questions. Although weekly quiz performance revealed an effect of retrieval practice-superior memory for material tested at the end of class (postquestions and new questions) compared to material not tested (quiz-only questions)-there was no difference in weekly quiz performance on postquestions vs. new questions. These results suggest that retrieval practice is beneficial to learning in the classroom. However, prequestions do not appear to enhance learning, nor to enhance the effects of retrieval practice.

Students' Use of Optional Online Reviews and Its Relationship to Summative Assessment Outcomes in Introductory Biology.

Retrieval practice has been shown to produce significant enhancements in student learning of course information, but the extent to which students make use of retrieval to learn information on their own is unclear. In the current study, students in a large introductory biology course were provided with optional online review questions that could be accessed as Test questions (requiring students to answer the questions before receiving feedback) or as Read questions (providing students with the question and correct answer up-front). Students more often chose to access the questions as Test compared with Read, and students who used the Test questions scored significantly higher on subsequent exams compared with students who used Read questions or did not access the questions at all. Following an in-class presentation of superior exam performance following use of the Test questions, student use of Test questions increased significantly for the remainder of the term. These results suggest that practice questions can be an effective tool for enhancing student achievement in biology and that informing students about performance-based outcomes coincides with increased use of retrieval practice.

Estimating error in diffusion coefficients derived from molecular dynamics simulations.

The computationally expensive nature of molecular dynamics simulation limits the access to length (nanometer) and time scales (nanosecond) that are orders of magnitude smaller than the experiment it models. This limitation warrants a careful estimation of statistical uncertainty associated with the properties calculated from these simulations. The assumption that a simulation is long enough so that the ergodic hypothesis applies is often invoked in the literature for the computation of properties of interest from a single molecular dynamics simulation. Here, we demonstrate that making this assumption without validation results in poor estimates of the self-diffusion coefficient from a single molecular dynamics simulation of Lennard-Jones fluid. This problem is shown to be even more severe when the diffusion coefficient of macromolecules is calculated from a single molecular dynamics simulation. We have shown that conducting multiple independent simulations is necessary to obtain reliable estimates of diffusion coefficients and their associated statistical uncertainties. We show that even a "routine" calculation of the self-diffusion coefficient for a Lennard-Jones fluid, as determined from a linear fit of the mean squared displacement of particles as a function of time, violates the key assumptions of linear regression. A rigorous approach for addressing these issues is presented.

A coarse-graining approach for molecular simulation that retains the dynamics of the all-atom reference system by implementing hydrodynamic interactions.

We report on a new approach for deriving coarse-grained intermolecular forces that retains the frictional contribution that is often discarded by conventional coarse-graining methods. The approach is tested for water and an aqueous glucose solution, and the results from the new implementation for coarse-grained molecular dynamics simulation show remarkable agreement with the dynamics obtained from reference all-atom simulations. The agreement between the structural properties observed in the coarse-grained and all-atom simulations is also preserved. We discuss how this approach may be applied broadly to any existing coarse-graining method where the coarse-grained models are rigorously derived from all-atom reference systems.

An adaptable pentaloop defines a robust neomycin-B RNA aptamer with conditional ligand-bound structures.

Aptamers can be highly specific for their targets, which implies precise molecular recognition between aptamer and target. However, as small polymers, their structures are more subject to environmental conditions than the more constrained longer RNAs such as those that constitute the ribosome. To understand the balance between structural and environmental factors in establishing ligand specificity of aptamers, we examined the RNA aptamer (NEO1A) previously reported as specific for neomycin-B. We show that NEO1A can recognize other aminoglycosides with similar affinities as for neomycin-B and its aminoglycoside specificity is strongly influenced by ionic strength and buffer composition. NMR and 2-aminopurine (2AP) fluorescence studies of the aptamer identified a flexible pentaloop and a stable binding pocket. Consistent with a well-structured binding pocket, docking analysis results correlated with experimental measures of the binding energy for most ligands. Steady state fluorescence studies of 2AP-substituted aptamers confirmed that A16 moves to a more solvent accessible position upon ligand binding while A14 moves to a less solvent accessible position, which is most likely a base stack. Analysis of binding affinities of NEO1A sequence variants showed that the base in position 16 interacts differently with each ligand and the interaction is a function of the buffer constituents. Our results show that the pentaloop provides NEO1A with the ability to adapt to external influences on its structure, with the critical base at position 16 adjusting to incorporate each ligand into a stable pocket by hydrophobic interactions and/or hydrogen bonds depending on the ligand and the ionic environment.

Molecular dynamics simulations on the Tre1 G protein-coupled receptor: exploring the role of the arginine of the NRY motif in Tre1 structure.

BACKGROUND: The arginine of the D/E/NRY motif in Rhodopsin family G protein-coupled receptors (GPCRs) is conserved in 96% of these proteins. In some GPCRs, this arginine in transmembrane 3 can form a salt bridge with an aspartic acid or glutamic acid in transmembrane 6. The Drosophila melanogaster GPCR Trapped in endoderm-1 (Tre1) is required for normal primordial germ cell migration. In a mutant form of the protein, Tre1sctt, eight amino acids RYILIACH are missing, resulting in a severe disruption of primordial germ cell development. The impact of the loss of these amino acids on Tre1 structure is unknown. Since the missing amino acids in Tre1sctt include the arginine that is part of the D/E/NRY motif in Tre1, molecular dynamics simulations were performed to explore the hypothesis that these amino acids are involved in salt bridge formation and help maintain Tre1 structure. RESULTS: Structural predictions of wild type Tre1 (Tre1+) and Tre1sctt were subjected to over 250 ns of molecular dynamics simulations. The ability of the model systems to form a salt bridge between the arginine of the D/E/NRY motif and an aspartic acid residue in transmembrane 6 was analyzed. The results indicate that a stable salt bridge can form in the Tre1+ systems and a weak salt bridge or no salt bridge, using an alternative arginine, is likely in the Tre1sctt systems. CONCLUSIONS: The weak salt bridge or lack of a salt bridge in the Tre1sctt systems could be one possible explanation for the disrupted function of Tre1sctt in primordial germ cell migration. These results provide a framework for studying the importance of the arginine of the D/E/NRY motif in the structure and function of other GPCRs that are involved in cell migration, such as CXCR4 in the mouse, zebrafish, and chicken.

Ab initio study of molecular interactions in cellulose Iα.

Biomass recalcitrance, the resistance of cellulosic biomass to degradation, is due in part to the stability of the hydrogen bond network and stacking forces between the polysaccharide chains in cellulose microfibers. The fragment molecular orbital (FMO) method at the correlated Møller-Plesset second order perturbation level of theory was used on a model of the crystalline cellulose Iα core with a total of 144 glucose units. These computations show that the intersheet chain interactions are stronger than the intrasheet chain interactions for the crystalline structure, while they are more similar to each other for a relaxed structure. An FMO chain pair interaction energy decomposition analysis for both the crystal and relaxed structures reveals an intricate interplay between electrostatic, dispersion, charge transfer, and exchange repulsion effects. The role of the primary alcohol groups in stabilizing the interchain hydrogen bond network in the inner sheet of the crystal and relaxed structures of cellulose Iα, where edge effects are absent, was analyzed. The maximum attractive intrasheet interaction is observed for the GT-TG residue pair with one intrasheet hydrogen bond, suggesting that the relative orientation of the residues is as important as the hydrogen bond network in strengthening the interaction between the residues.

Evaluation of coarse-grained mapping schemes for polysaccharide chains in cellulose.

A fundamental understanding of the intermolecular forces that bind polysaccharide chains together in cellulose is crucial for designing efficient methods to overcome the recalcitrance of lignocellulosic biomass to hydrolysis. Because the characteristic time and length scales for the degradation of cellulose by enzymatic hydrolysis or chemical pretreatment span orders of magnitude, it is important to closely integrate the molecular models used at each scale so that, ultimately, one may switch seamlessly between quantum, atomistic, and coarse-grained descriptions of the system. As a step towards that goal, four multiscale coarse-grained models for polysaccharide chains in a cellulose-Iα microfiber are considered. Using the force matching method, effective coarse-grained forces are derived from all-atom trajectories. Performance of the coarse-grained models is evaluated by comparing the intrachain radial distribution functions with those obtained using the all-atom reference data. The all-atom simulation reveals a double peak in the radial distribution function for sites within each glucose residue that arises from the distinct conformations sampled by the primary alcohol group in the glucose residues. The three-site and four-site coarse-grained models have sufficient degrees of freedom to predict this double peak while the one-site and two-site models do not. This is the first time that coarse-grained models have been shown to reproduce such subtle, yet important, molecular features in a polysaccharide chain. The relative orientations between glucose residues along the polysaccharide chain are evaluated and it is found that the four-site coarse-grained model is best at reproducing the glucose-glucose conformations observed in the all-atom simulation. The success of the four-site coarse-grained model underscores the importance of decoupling the pyranose ring from the oxygen atom in the glycosidic bond when developing all-atom to coarse-grained mapping schemes for polysaccharides.

Investigating the malleability of RNA aptamers.

Aptamers are short, single-stranded nucleic acids with structures that frequently change upon ligand binding and are sensitive to the ionic environment. To achieve facile application of aptamers in controlling cellular activities, a better understanding is needed of aptamer ligand binding parameters, structures, intramolecular mobilities and how these structures adapt to different ionic environments with consequent effects on their ligand binding characteristics. Here we discuss the integration of biochemical analysis with NMR spectroscopy and computational modeling to explore the relation between ligand binding and structural malleability of some well-studied aptamers. Several methods for determining aptamer binding affinity and specificity are discussed, including isothermal titration calorimetry, steady state fluorescence of 2-aminopurine substituted aptamers, and dye displacement assays. Also considered are aspects of molecular dynamics simulations specific to aptamers including adding ions and simulating aptamer structure in the absence of ligand when NMR spectroscopy or X-ray crystallography structures of the unoccupied aptamer are not available. We focus specifically on RNA aptamers that bind small molecule ligands as would be applied in sensors or integrated into riboswitches such as to measure the products of metabolic activity.

Dendrimer-Fullerenol Soft-Condensed Nanoassembly.

Nanoscale assembly is an area of research that has vast implications for molecular design, sensing, nanofabrication, supramolecular chemistry, catalysis, and environmental remediation. Here we show that poly(amidoamine) (PAMAM) dendrimers of both generations 1 (G1) and 4 (G4) can host 1 fullerenol per 2 dendrimer primary amines as evidenced by isothermal titration calorimetry, dynamic light scattering and spectrofluorometry. Thermodynamically, the interactions were similarly spontaneous between both generations of dendrimers and fullerenols, however, G4 formed stronger complexes with fullerenols resulting from their higher surface charge density and more internal voids, as demonstrated by spectrofluorometry. In addition to hydrogen bonding that existed between the dendrimer primary amines and the fullerenol oxygens, hydrophobic and electrostatic interactions also contributed to complex formation and dynamics. Such hybrid of soft and condensed nanoassembly may have implications for environmental remediation of discharged nanomaterials and entail new applications in drug delivery.

Reintroducing explicit solvent to a solvent-free coarse-grained model.

A unique coarse-grained modeling scheme that combines a systematic, solvent-free multiscale coarse-graining algorithm for a complex macromolecule with an existing coarse-grained solvent model is proposed. We show that this procedure efficiently and reliably describes the interactions for complex macromolecules, using the specific example of dendrimers binding phenanthrenes in water. The experimentally measured binding capacity is predicted by the unique coarse-grained modeling approach; the conditions for this simulation are beyond what could be reasonably simulated with an all-atom molecular dynamics simulation.

A biophysical perspective of understanding nanoparticles at large.

In this article we present a biophysical perspective that describes the fate of nanoparticles in both the aqueous phase and in living systems. Specifically, we show the correlations between the physicochemistry of fullerenes and their uptake, translocation, transformation, transport, and biodistribution in mammalian and plant systems, at the molecular, cellular, and whole organism level. In addition to fullerenes and their structural derivatives, we describe the biological and environmental implications and applications of the condensed matter of carbon nanotubes and quantum dots, and the soft condensed matter of plastic and dendrimers. The main purpose of this article is to demonstrate the vast opportunities and unique advantages of applying experimental and simulation biophysics to the nascent research field of understanding nanoparticles at large.

Fluorescence resonance energy transfer between phenanthrene and PAMAM dendrimers.

We describe herein an adsorption-induced energy transfer between phenanthrene, a major environmental pollutant, and a fluorescently labeled dendrimer acting as a host molecule. We find experimentally that such energy transfer is the most efficient at a solvent pH of 8 and for a phenanthrene:dendrimer molar ratio of 1:2. Using molecular dynamics simulations we show that the strongest binding interactions occur between phenanthrene and the primary amines of the dendrimer. The simulations provide evidence that at low pH, phenanthrene-phenanthrene interactions are favorable and compete with phenanthrene-dendrimer binding. This study offers a new scheme for detecting dendrimer molecular assembly and a physical basis for exploiting dendrimer nanotechnologies for water purification and environmental remediation.

Cell trafficking of carbon nanotubes based on fluorescence detection.

Cell trafficking of carbon nanotubes (CNTs) is an area of scientific inquiry that has great implications in medicine, biosensing, and environmental science and engineering. The essence of this inquiry resides in the interaction of carbon nanostructures and cell membranes, regulated by the laws of molecular cell biology and the physiochemical properties of the nanostructures. Of equal importance to this inquiry is a description of cellular responses to the integration of man-made materials; yet, how cellular responses may invoke whole-organism level reaction remains unclear. In this chapter, we show three experimental studies, which may be beneficial to obtaining such an understanding. Among the reservoir of methodologies, which have proved of merit, we focus our attention on fluorescence microscopy, one of the most powerful and yet least invasive means of probing nanoparticles in biological systems. Especially, we present the method of fluorescence energy transfer induced between a lysophospholipid molecule and a single-walled CNT upon cellular uptake, and describe coating nanotubes with RNA and suspending fullerenes with phenolic acids for facilitating their translocation across cell membranes and shuttling between cell organelles. Finally, we comment on the perspective of using molecular simulations for facilitating and guiding such experiments.

Molecular dynamics simulation of fractal aggregate diffusion.

The diffusion of fractal aggregates constructed with the method by Thouy and Jullien [J. Phys. A 27, 2953 (1994)] comprised of N(p) spherical primary particles was studied as a function of the aggregate mass and fractal dimension using molecular dynamics simulations. It is shown that finite-size effects have a strong impact on the apparent value of the diffusion coefficient (D), but these can be corrected by carrying out simulations using different simulation box sizes. Specifically, the diffusion coefficient is inversely proportional to the length of a cubic simulation box, and the constant of proportionality appears to be independent of the aggregate mass and fractal dimension. Using this result, it is possible to compute infinite dilution diffusion coefficients (D(o)) for aggregates of arbitrary size and fractal dimension, and it was found that D(o)∝N(p)(-1/df), as is often assumed by investigators simulating brownian aggregation of fractal aggregates. The ratio of hydrodynamic radius to radius of gyration is computed and shown to be independent of mass for aggregates of fixed fractal dimension, thus enabling an estimate of the diffusion coefficient for a fractal aggregate based on its radius of gyration.