Search

Despite the popularity of computer-aided study and design of RNA molecules, little is known about the accuracy of commonly used structure modeling packages in tasks sensitive to ensemble properties of RNA. Here, we demonstrate that the EternaBench dataset, a set of more than 20,000 synthetic RNA constructs designed on the RNA design platform Eterna, provides incisive discriminative power in evaluating current packages in ensemble-oriented structure prediction tasks. We find that CONTRAfold and RNAsoft, packages with parameters derived through statistical learning, achieve consistently higher accuracy than more widely used packages in their standard settings, which derive parameters primarily from thermodynamic experiments. We hypothesized that training a multitask model with the varied data types in EternaBench might improve inference on ensemble-based prediction tasks. Indeed, the resulting model, named EternaFold, demonstrated improved performance that generalizes to diverse external datasets including complete messenger RNAs, viral genomes probed in human cells and synthetic designs modeling mRNA vaccines.

AlphaFold2 (AF2) has revolutionized structural biology by accurately predicting single structures of proteins and protein-protein complexes. However, biological function is rooted in a protein’s ability to sample different conformational substates, and disease-causing point mutations are often due to population changes of these substates. This has sparked immense interest in expanding AF2’s capability to predict conformational substates. We demonstrate that clustering an input multiple sequence alignment (MSA) by sequence similarity enables AF2 to sample alternate states of known metamorphic proteins, including the circadian rhythm protein KaiB, the transcription factor RfaH, and the spindle checkpoint protein Mad2, and score these states with high confidence. Moreover, we use AF2 to identify a minimal set of two point mutations predicted to switch KaiB between its two states. Finally, we used our clustering method, AF-cluster, to screen for alternate states in protein families without known fold-switching, and identified a putative alternate state for the oxidoreductase DsbE. Similarly to KaiB, DsbE is predicted to switch between a thioredoxin-like fold and a novel fold. This prediction is the subject of future experimental testing. Further development of such bioinformatic methods in tandem with experiments will likely have profound impact on predicting protein energy landscapes, essential for shedding light into biological function.

RNA hydrolysis presents problems in manufacturing, long-term storage, world-wide delivery, and in vivo stability of messenger RNA (mRNA)-based vaccines and therapeutics. A largely unexplored strategy to reduce mRNA hydrolysis is to redesign RNAs to form double-stranded regions, which are protected from in-line cleavage and enzymatic degradation, while coding for the same proteins. The amount of stabilization that this strategy can deliver and the most effective algorithmic approach to achieve stabilization remain poorly understood. Here, we present simple calculations for estimating RNA stability against hydrolysis, and a model that links the average unpaired probability of an mRNA, or AUP, to its overall hydrolysis rate. To characterize the stabilization achievable through structure design, we compare AUP optimization by conventional mRNA design methods to results from more computationally sophisticated algorithms and crowdsourcing through the OpenVaccine challenge on the Eterna platform. These computational tests were carried out on both model mRNAs and COVID-19 mRNA vaccine candidates. We find that rational design on Eterna and the more sophisticated algorithms lead to constructs with low AUP, which we term ‘superfolder’ mRNAs. These designs exhibit wide diversity of sequence and structure features that may be desirable for translation, biophysical size, and immunogenicity, and their folding is robust to temperature, choice of flanking untranslated regions, and changes in target protein sequence, as illustrated by rapid redesign of superfolder mRNAs for B.1.351, P.1, and B.1.1.7 variants of the prefusion-stabilized SARS-CoV-2 spike protein. Increases in in vitro mRNA half-life by at least two-fold appear immediately achievable.Significance statementMessenger RNA (mRNA) medicines that encode and promote translation of a target protein have shown promising use as vaccines in the current SARS-CoV-2 pandemic as well as infectious diseases due to their speed of design and manufacturing. However, these molecules are intrinsically prone to hydrolysis, leading to poor stability in aqueous buffer and major challenges in distribution. Here, we present a principled biophysical model for predicting RNA degradation, and demonstrate that the stability of any mRNA can be increased at least two-fold over conventional design techniques. Furthermore, the predicted stabilization is robust to post-design modifications. This conceptual framework and accompanying algorithm can be immediately deployed to guide re-design of mRNA vaccines and therapeutics to increase in vitro stability.

As the COVID-19 outbreak spreads, there is a growing need for a compilation of conserved RNA genome regions in the SARS-CoV-2 virus along with their structural propensities to guide development of antivirals and diagnostics. Here we present a first look at RNA sequence conservation and structural propensities in the SARS-CoV-2 genome. Using sequence alignments spanning a range of betacoronaviruses, we rank genomic regions by RNA sequence conservation, identifying 79 regions of length at least 15 nt as exactly conserved over SARS-related complete genome sequences available near the beginning of the COVID-19 outbreak. We then confirm the conservation of the majority of these genome regions across 739 SARS-CoV-2 sequences subsequently reported from the COVID-19 outbreak, and we present a curated list of 30 “SARS-related-conserved” regions. We find that known RNA structured elements curated as Rfam families and in prior literature are enriched in these conserved genome regions, and we predict additional conserved, stable secondary structures across the viral genome. We provide 106 “SARS-CoV-2-conserved-structured” regions as potential targets for antivirals that bind to structured RNA. We further provide detailed secondary structure models for the extended 5′ UTR, frameshifting stimulation element, and 3′ UTR. Lastly, we predict regions of the SARS-CoV-2 viral genome that have low propensity for RNA secondary structure and are conserved within SARS-CoV-2 strains. These 59 “SARS-CoV-2-conserved-unstructured” genomic regions may be most easily accessible by hybridization in primer-based diagnostic strategies.

Internet-based scientific communities promise a means to apply distributed, diverse human intelligence toward previously intractable scientific problems. However, current implementations have not allowed communities to propose experiments to test all emerging hypotheses at scale or to modify hypotheses in response to experiments. We report high-throughput methods for molecular characterization of nucleic acids that enable the large-scale video game–based crowdsourcing of RNA sensor design, followed by high-throughput functional characterization. Iterative design testing of thousands of crowdsourced RNA sensor designs produced near–thermodynamically optimal and reversible RNA switches that act as self-contained molecular sensors and couple five distinct small molecule inputs to three distinct protein binding and fluorogenic outputs. This work suggests a paradigm for widely distributed experimental bioscience.

Medicines based on messenger RNA (mRNA) hold immense potential, as evidenced by their rapid deployment as COVID-19 vaccines. However, worldwide distribution of mRNA molecules has been limited by their thermostability, which is fundamentally limited by the intrinsic instability of RNA molecules to a chemical degradation reaction called in-line hydrolysis. Predicting the degradation of an RNA molecule is a key task in designing more stable RNA-based therapeutics. Here, we describe a crowdsourced machine learning competition (‘Stanford OpenVaccine’) on Kaggle, involving single-nucleotide resolution measurements on 6,043 diverse 102–130-nucleotide RNA constructs that were themselves solicited through crowdsourcing on the RNA design platform Eterna. The entire experiment was completed in less than 6 months, and 41% of nucleotide-level predictions from the winning model were within experimental error of the ground truth measurement. Furthermore, these models generalized to blindly predicting orthogonal degradation data on much longer mRNA molecules (504–1,588 nucleotides) with improved accuracy compared with previously published models. These results indicate that such models can represent in-line hydrolysis with excellent accuracy, supporting their use for designing stabilized messenger RNAs. The integration of two crowdsourcing platforms, one for dataset creation and another for machine learning, may be fruitful for other urgent problems that demand scientific discovery on rapid timescales.

International audience; The residence time of a drug on its target has been suggested as a more pertinent metric of therapeutic efficacy than the traditionally used affinity constant. Here we introduce junctured-DNA (J-DNA) tweezers as a generic platform which enables real-time observation, at the single-molecule level, of biomolecular interactions. This tool corresponds to a double-strand DNA scaffold which can be nanomanipulated and on which proteins of interest can be engrafted thanks to widely-used genetic tagging strategies. Thus, J-DNA tweezers allow straightforward and robust access to single molecule force spectroscopy in drug discovery, and more generally in biophysics. Proof-of-principle experiments are provided for the rapamycin-mediated association between FKBP12 and FRB, a system relevant in both medicine and chemical biology. Individual interactions were monitored under a range of applied forces and temperatures, yielding after analysis the characteristic features of the energy profile along the dissociation landscape.

The rational design of RNA is becoming important for rapidly developing technologies in medicine and biochemistry. Recent work has led to the development of several RNA secondary structure design algorithms and corresponding benchmarks to evaluate their performance. However, the performance of these algorithms is linked to the nature of the underlying algorithms for predicting secondary structure from sequences. Here, we show that an online community of RNA design experts is capable of modifying an existing RNA secondary structure design benchmark (Eterna100) with minimal alterations to address changes in the folding engine used (Vienna 1.8 updated to Vienna 2.4). We tested this new Eterna100-V2 benchmark with five RNA design algorithms, and found that neural network-based methods exhibited reduced performance in the folding engine they were evaluated on in their respective papers. We investigated this discrepancy, and determined that structural features, previously classified as difficult, may be dependent on parameters inherent to the RNA energy function itself. These findings suggest that for optimal performance, future algorithms should focus on finding strategies capable of solving RNA secondary structure design benchmarks independently of the free energy benchmark used. Eterna100-V1 and Eterna100-V2 benchmarks and example solutions are freely available at https://github.com/eternagame/eterna100-benchmarking.

Therapeutic mRNAs and vaccines are being developed for a broad range of human diseases, including COVID-19. However, their optimization is hindered by mRNA instability and inefficient protein expression. Here, we describe design principles that overcome these barriers. We develop an RNA sequencing-based platform called PERSIST-seq to systematically delineate in-cell mRNA stability, ribosome load, as well as in-solution stability of a library of diverse mRNAs. We find that, surprisingly, in-cell stability is a greater driver of protein output than high ribosome load. We further introduce a method called In-line-seq, applied to thousands of diverse RNAs, that reveals sequence and structure-based rules for mitigating hydrolytic degradation. Our findings show that highly structured “superfolder” mRNAs can be designed to improve both stability and expression with further enhancement through pseudouridine nucleoside modification. Together, our study demonstrates simultaneous improvement of mRNA stability and protein expression and provides a computational-experimental platform for the enhancement of mRNA medicines.

The computer-aided study and design of RNA molecules is increasingly prevalent across a range of disciplines, yet little is known about the accuracy of commonly used structure modeling packages in tasks sensitive to ensemble properties of RNA. Here, we demonstrate that the EternaBench dataset, a set of over 20,000 synthetic RNA constructs designed in iterative cycles on the RNA design platform Eterna, provides incisive discriminative power in evaluating current packages in ensemble-oriented structure prediction tasks. We find that CONTRAfold and RNAsoft, packages with parameters derived through statistical learning, achieve consistently higher accuracy than more widely used packages in their standard settings, which derive parameters primarily from thermodynamic experiments. Motivated by these results, we develop a multitask-learning-based model, EternaFold, which demonstrates improved performance that generalizes to diverse external datasets, including complete mRNAs and viral genomes probed in human cells and synthetic designs modeling mRNA vaccines.

Internet-based scientific communities promise a means to apply distributed, diverse human intelligence towards previously intractable scientific problems. However, current implementations have not allowed communities to propose experiments to test all emerging hypotheses at scale or to modify hypotheses in response to experiments. We report high-throughput methods for molecular characterization of nucleic acids that enable the large-scale videogame-based crowdsourcing of functional RNA sensor design, followed by high-throughput functional characterization. Iterative design testing of thousands of crowdsourced RNA sensor designs produced near-thermodynamically optimal and reversible RNA switches that act as self-contained molecular sensors and couple five distinct small molecule inputs to three distinct protein binding and fluorogenic outputs—results that surpass computational and expert-based design. This work represents a new paradigm for widely distributed experimental bioscience.One Sentence SummaryOnline community discovers standalone RNA sensors.

Riboswitches are RNA elements that recognize diverse chemical and biomolecular inputs, and transduce this recognition process to genetic, fluorescent, and other engineered outputs using RNA conformational changes. These systems are pervasive in cellular biology and are a promising biotechnology with applications in genetic regulation and biosensing. Here, we derive a simple expression bounding the activation ratio—the proportion of RNA in the active vs. inactive states—for both ON and OFF riboswitches that operate near thermodynamic equilibrium: [Formula: see text] , where [I] is the input ligand concentration and [Formula: see text] is the intrinsic dissociation constant of the aptamer module toward the input ligand. A survey of published studies of natural and synthetic riboswitches confirms that the vast majority of empirically measured activation ratios have remained well below this thermodynamic limit. A few natural and synthetic riboswitches achieve activation ratios close to the limit, and these molecules highlight important principles for achieving high riboswitch performance. For several applications, including “light-up” fluorescent sensors and chemically-controlled CRISPR/Cas complexes, the thermodynamic limit has not yet been achieved, suggesting that current tools are operating at suboptimal efficiencies. Future riboswitch studies will benefit from comparing observed activation ratios to this simple expression for the optimal activation ratio. We present experimental and computational suggestions for how to make these quantitative comparisons and suggest new molecular mechanisms that may allow non-equilibrium riboswitches to surpass the derived limit.

RNA hydrolysis presents problems in manufacturing, long-term storage, world-wide delivery and in vivo stability of messenger RNA (mRNA)-based vaccines and therapeutics. A largely unexplored strategy to reduce mRNA hydrolysis is to redesign RNAs to form double-stranded regions, which are protected from in-line cleavage and enzymatic degradation, while coding for the same proteins. The amount of stabilization that this strategy can deliver and the most effective algorithmic approach to achieve stabilization remain poorly understood. Here, we present simple calculations for estimating RNA stability against hydrolysis, and a model that links the average unpaired probability of an mRNA, or AUP, to its overall hydrolysis rate. To characterize the stabilization achievable through structure design, we compare AUP optimization by conventional mRNA design methods to results from more computationally sophisticated algorithms and crowdsourcing through the OpenVaccine challenge on the Eterna platform. We find that rational design on Eterna and the more sophisticated algorithms lead to constructs with low AUP, which we term 'superfolder' mRNAs. These designs exhibit a wide diversity of sequence and structure features that may be desirable for translation, biophysical size, and immunogenicity. Furthermore, their folding is robust to temperature, computer modeling method, choice of flanking untranslated regions, and changes in target protein sequence, as illustrated by rapid redesign of superfolder mRNAs for B.1.351, P.1 and B.1.1.7 variants of the prefusion-stabilized SARS-CoV-2 spike protein. Increases in in vitro mRNA half-life by at least two-fold appear immediately achievable.

Describing the dynamics and conformational landscapes of Intrinsically Disordered Proteins (IDPs) is of paramount importance to understanding their functions. Markov State Models (MSMs) are often used to characterize the dynamics of more structured proteins, but models of IDPs built using conventional MSM modelling protocols can be difficult to interpret due to the inherent nature of IDPs, which exhibit fast transitions between disordered microstates. We propose a new method of determining MSM states from all-atom molecular dynamics simulation data of IDPs by using per-residue secondary structure assignments as input features in a MSM model. Because such secondary structure algorithms use a select set of features for assignment (dihedral angles, contact distances, etc.), they represent a knowledge-based refinement of feature sets used for model-building. This method adds interpretability to IDP conformational landscapes, which are increasingly viewed as composed of transient secondary structure, and allows us to readily use MSM analysis tools in this paradigm. We demonstrate the use of our method with the transcription factor p53 c-terminal domain (p53-CTD), a commonly-studied IDP. We are able to characterize the full secondary structure phase space observed for p53-CTD, and describe characteristics of p53-CTD as a network of transient helical and beta-hairpin structures with different network behaviors in different domains of secondary structure. This analysis provides a novel example of how IDPs can be studied and how researchers might better understand a disordered protein conformational landscape.

Aluminum has attracted great attention recently as it has been suggested by several studies to be associated with increased risks for Alzheimer's and Parkinson's disease. The toxicity of the trivalent ion is assumed to derive from structural changes induced in lipid bilayers upon binding, though the mechanism of this process is still not well understood. In the present study we elucidate the effect of Al(3+) on supported lipid bilayers (SLBs) using fluorescence microscopy, the quartz crystal microbalance with dissipation (QCM-D) technique, dual-polarization interferometry (DPI), and molecular dynamics (MD) simulations. Results from these techniques show that binding of Al(3+) to SLBs containing negatively charged and neutral phospholipids induces irreversible changes such as domain formation. The measured variations in SLB thickness, birefringence, and density indicate a phase transition from a disordered to a densely packed ordered phase.

Variational auto-encoder frameworks have demonstrated success in reducing complex nonlinear dynamics in molecular simulation to a single non-linear embedding. In this work, we illustrate how this non-linear latent embedding can be used as a collective variable for enhanced sampling, and present a simple modification that allows us to rapidly perform sampling in multiple related systems. We first demonstrate our method is able to describe the effects of force field changes in capped alanine dipeptide after learning a model using AMBER99. We further provide a simple extension to variational dynamics encoders that allows the model to be trained in a more efficient manner on larger systems by encoding the outputs of a linear transformation using time-structure based independent component analysis (tICA). Using this technique, we show how such a model trained for one protein, the WW domain, can efficiently be transferred to perform enhanced sampling on a related mutant protein, the GTT mutation. This method shows promise for its ability to rapidly sample related systems using a single transferable collective variable and is generally applicable to sets of related simulations, enabling us to probe the effects of variation in increasingly large systems of biophysical interest., Comment: 20 pages, 10 figures

As deep Variational Auto-Encoder (VAE) frameworks become more widely used for modeling biomolecular simulation data, we emphasize the capability of the VAE architecture to concurrently maximize the time scale of the latent space while inferring a reduced coordinate, which assists in finding slow processes as according to the variational approach to conformational dynamics. We provide evidence that the VDE framework [Hernandez et al., Phys. Rev. E 97, 062412 (2018)], which uses this autocorrelation loss along with a time-lagged reconstruction loss, obtains a variationally optimized latent coordinate in comparison with related loss functions. We thus recommend leveraging the autocorrelation of the latent space while training neural network models of biomolecular simulation data to better represent slow processes.

Markov state models (MSMs) are a powerful framework for the analysis of molecular dynamics data sets, such as protein folding simulations, because of their straightforward construction and statistical rigor. The coarse-graining of MSMs into an interpretable number of macrostates is a crucial step for connecting theoretical results with experimental observables. Here we present the minimum variance clustering approach (MVCA) for the coarse-graining of MSMs into macrostate models. The method utilizes agglomerative clustering with Ward's minimum variance objective function, and the similarity of the microstate dynamics is determined using the Jensen-Shannon divergence between the corresponding rows in the MSM transition probability matrix. We first show that MVCA produces intuitive results for a simple tripeptide system and is robust toward long-duration statistical artifacts. MVCA is then applied to two protein folding simulations of the same protein in different force fields to demonstrate that a different number of macrostates is appropriate for each model, revealing a misfolded state present in only one of the simulations. Finally, we show that the same method can be used to analyze a data set containing many MSMs from simulations in different force fields by aggregating them into groups and quantifying their dynamical similarity in the context of force field parameter choices. The minimum variance clustering approach with the Jensen-Shannon divergence provides a powerful tool to group dynamics by similarity, both among model states and among dynamical models themselves.

Zeolitic imidazolate frameworks (ZIFs) have been widely investigated for numerous applications including energy storage, heterogeneous catalysis, and greenhouse gas adsorption. Much of the early work has focused on the bulk properties of microcrystalline ZIFs. Herein, we focus on identifying the nature of the surface of ZIF-8 by studying a supported ZIF-8 nanoparticle film using surface characterization techniques. We have experimentally identified the presence of a zinc-rich surface terminated by carbonates and water/hydroxyl groups (in addition to the expected methylimidazole terminations) using X-ray photoelectron spectroscopy (XPS). The thermal stability of ZIF-8 thin films was also investigated using scanning electron microscopy (SEM) and temperature-programmed reaction spectroscopy (TPRS). We determined the onset of decomposition of ZIF-8 thin films to be approximately 630 K using TPRS in an ultrahigh vacuum (UHV) environment. This work presents the first characterization steps needed to study the e...

In the standard DNA brick set-up, distinct 32-nucleotide strands of single-stranded DNA are each designed to bind specifically to four other such molecules. Experimentally, it has been demonstrated that the overall yield is increased if certain bricks which occur on the outer faces of target structures are merged with adjacent bricks. However, it is not well understood by what mechanism such 'boundary bricks' increase the yield, as they likely influence both the nucleation process and the final stability of the target structure. Here, we use Monte Carlo simulations with a patchy particle model of DNA bricks to investigate the role of boundary bricks in the self-assembly of complex multicomponent target structures. We demonstrate that boundary bricks lower the free-energy barrier to nucleation and that boundary bricks on edges stabilize the final structure. However, boundary bricks are also more prone to aggregation, as they can stabilize partially assembled intermediates. We explore some design strategies that permit us to benefit from the stabilizing role of boundary bricks whilst minimizing their ability to hinder assembly; in particular, we show that maximizing the total number of boundary bricks is not an optimal strategy.

Often the analysis of time-dependent chemical and biophysical systems produces high-dimensional time-series data for which it can be difficult to interpret which individual features are most salient. While recent work from our group and others has demonstrated the utility of time-lagged co-variate models to study such systems, linearity assumptions can limit the compression of inherently nonlinear dynamics into just a few characteristic components. Recent work in the field of deep learning has led to the development of variational autoencoders (VAE), which are able to compress complex datasets into simpler manifolds. We present the use of a time-lagged VAE, or variational dynamics encoder (VDE), to reduce complex, nonlinear processes to a single embedding with high fidelity to the underlying dynamics. We demonstrate how the VDE is able to capture nontrivial dynamics in a variety of examples, including Brownian dynamics and atomistic protein folding. Additionally, we demonstrate a method for analyzing the VDE model, inspired by saliency mapping, to determine what features are selected by the VDE model to describe dynamics. The VDE presents an important step in applying techniques from deep learning to more accurately model and interpret complex biophysics., Comment: Fixed typos and added references

Formation and quality of single solid supported lipid membranes and double lipid membranes were investigated with single vesicle resolution using label-free evanescence light scattering microscopy (EvSM). For the formation of double lipid membranes we made use of electrostatic interaction between charged lipids and oppositely charged cations.

Understanding the kinetics of dye adsorption and desorption on semiconductors is crucial for optimizing the performance of dye-sensitized solar cells (DSSCs). Quartz crystal microbalance with dissipation monitoring (QCM-D) measures adsorbed mass in real time, allowing determination of binding kinetics. In this work, we characterize adsorption of the common RuBipy dye N3 to the native oxide layer of a planar, sputter-coated titanium surface, simulating the TiO2 substrate of a DSSC. We report adsorption equilibrium constants consistent with prior optical measurements of N3 adsorption. Dye binding and surface integrity were also verified by scanning electron microscopy, energy-dispersive X-ray spectroscopy, and X-ray photoelectron spectroscopy (XPS). We further study desorption of the dye from the native oxide layer on the QCM sensors using tetrabutylammonium hydroxide (TBAOH), a commonly used industrial desorbant. We find that using TBAOH as a desorbant does not fully regenerate the surface, though little ruthenium or nitrogen is observed by XPS after desorption, suggesting that carboxyl moieties of N3 remain bound. We demonstrate the native oxide layer of a titanium sensor as a valid and readily available planar TiO2 morphology to study dye adsorption and desorption and begin to investigate the mechanism of dye desorption in DSSCs, a system that requires further study.

Understanding the kinetics of dye adsorption on semiconductors is crucial for designing dye-sensitized solar cells (DSSCs) with enhanced efficiency. Harms et al. recently applied the Quartz-Crystal Microbalance with Dissipation Monitoring (QCM-D) to study in situ dye adsorption on flat TiO2 surfaces. QCM-D measures adsorption in real time and therefore allows one to determine the kinetics of the process. In this work, we characterize the adsorption of N3, a commercial RuBipy dye, using the native oxide layer of a titanium sensor to simulate the TiO2 substrate of a DSSC. We report equilibrium constants that are in agreement with previous absorbance studies of N3 adsorption, and therefore demonstrate the native oxide layer of a titanium sensor as a valid and readily available planar TiO2 morphology to study dye adsorption.