Search

While CRISPR/Cas9 is a powerful tool in genome engineering, the ontarget activity and off-target effects of the system widely vary because of the differences in guide RNA (gRNA) sequences and genomic environments. Traditional approaches rely on separate models and parameters to treat on- and off-target cleavage activities. Here, we demonstrate that a free-energy scheme dominates the Cas9 editing efficacy and delineate a method that simultaneously considers on-target activities and off-target effects. While data-driven machinelearning approaches learn rules to model particular datasets, they may not be as transferrable to new systems or capable of producing new mechanistic insights as principled physical approaches. By integrating the energetics of R-loop formation under Cas9 binding, the effect of the protospacer adjacent motif sequence, and the folding stability of the whole single guide RNA, we devised a unified, physical model that can apply to any cleavage-activity dataset. This unified framework improves predictions for both on-target activities and offtarget efficiencies of spCas9 and may be readily transferred to other systems with different guide RNAs or Cas9 ortholog proteins. [ABSTRACT FROM AUTHOR]

RNA functions are intrinsically tied tofolding kinetics. The most elementary step in RNA folding is the closing and opening of a base pair. Understanding this elementary rate process is the basis for RNA folding kinetics studies. Previous studies mostly focused on the unfolding of base pairs. Here, based on a hybrid approach, we investigate the folding process at level of single base pairing/stacking. The study, which integrates molecular dynamics simulation, kinetic Monte Carlo simulation, and master equation methods, uncovers two alternative dominant pathways: Starting from the unfolded state, the nucleotide backbone first folds to the native conformation, followed by subsequent adjustment of the base conformation. During the base conformational rearrangement, the backbone either retains the native conformation or switches to nonnative conformations in order to lower the kinetic barrier for base rearrangement. The method enables quantification of kinetic partitioning among the different pathways. Moreover, the simulation reveals several intriguing ion binding/ dissociation signatures for the conformational changes. Our approach may be useful for developing a base pair opening/closing rate model. [ABSTRACT FROM AUTHOR]

Examines the RNA hairpin-folding kinetics. Creation of statistical mechanics model to combine eigenvalue solution; Visualization of distinct time separation and rate of limiting step; Involvement of trapping and detrapping in folding of kinetic intermediates.