Phosphonoacetate Modifications Enhance the Stability and Editing Yields of Guide RNAs for Cas9 Editors.

CRISPR gene editing and control systems continue to emerge and inspire novel research and clinical applications. Advances in CRISPR performance such as optimizing the duration of activity in cells, tissues, and organisms, as well as limiting off-target activities, have been extremely important for expanding the utility of CRISPR-based systems. By investigating the effects of various chemical modifications in guide RNAs (gRNAs) at defined positions and combinations, we find that 2'-O-methyl-3'-phosphonoacetate (MP) modifications can be substantially more effective than 2'-O-methyl-3'-phosphorothioate (MS) modifications at the 3' ends of single-guide RNAs (sgRNAs) to promote high editing yields, in some instances showing an order of magnitude higher editing yield in human cells. MP-modified 3' ends are especially effective at promoting the activity of guide RNAs cotransfected with Cas messenger RNA (mRNA), as the gRNA must persist in cells until the Cas protein is expressed. We demonstrate such an MP enhancement for sgRNAs cotransfected with a BE4 mRNA for cytidine base editing and also demonstrate that MP at the 3' ends of prime editing guide RNAs (pegRNAs) cotransfected with PE2 mRNA can promote maximal prime editing yields. In the presence of serum, sgRNAs with MP-modified 3' ends showed marked improvements in editing efficiency over sgRNAs with MS-modified 3' ends codelivered with Cas9 mRNA and showed more modest improvements at enhancing the activity of transfected ribonucleoprotein (RNP) complexes. Our results suggest that MP should be considered as a performance-enhancing modification for the 3' ends of synthetic gRNAs, especially in situations where the guide RNAs may be susceptible to exonuclease-mediated degradation.

Biochemical Characterization of Yeast Xrn1.

Messenger RNA degradation is an important component of overall gene expression. During the final step of eukaryotic mRNA degradation, exoribonuclease 1 (Xrn1) carries out 5' → 3' processive, hydrolytic degradation of RNA molecules using divalent metal ion catalysis. To initiate studies of the 5' → 3' RNA decay machinery in our lab, we expressed a C-terminally truncated version of Saccharomyces cerevisiae Xrn1 and explored its enzymology using a second-generation, time-resolved fluorescence RNA degradation assay. Using this system, we quantitatively explored Xrn1's preference for 5'-monophosphorylated RNA substrates, its pH dependence, and the importance of active site mutations in the molecule's conserved catalytic core. Furthermore, we explore Xrn1's preference for RNAs containing a 5' single-stranded region both in an intermolecular hairpin structure and in an RNA-DNA hybrid duplex system. These results both expand and solidify our understanding of Xrn1, a centrally important enzyme whose biochemical properties have implications in numerous RNA degradation and processing pathways.

Pore diameter effects on phase behavior of a gas condensate in graphitic one-and two-dimensional nanopores.

Molecular dynamics (MD) simulations were performed on a hydrocarbon mixture representing a typical gas condensate composed mostly of methane and other small molecules with small fractions of heavier hydrocarbons, representative of mixtures found in tight shale reservoirs. The fluid was examined both in bulk and confined to graphitic nano-scale slits and pores. Numerous widths and diameters of slits and pores respectively were examined under variable pressures at 300 K in order to find conditions in which the fluid at the center of the apertures would not be affected by capillary condensation due to the oil-wet walls. For the bulk fluid, retrograde phase behavior was verified by liquid volumes obtained from Voronoi tessellations. In cases of both one and two-dimensional confinement, for the smallest apertures, heavy molecules aggregated inside the pore space and compression of the gas outside the solid structure lead to decreases in density of the confined fluid. Normal density/pressure relationships were observed for slits having gaps of above 3 nm and pores having diameters above 6 nm. At 70 bar, the minimum gap width at which the fluid could pass through the center of slits without condensation effects was predicted to be 6 nm and the corresponding diameter in pores was predicted to be 8 nm. The models suggest that in nanoscale networks involving pores smaller than these limiting dimensions, capillary condensation should significantly impede transmission of natural gases with similar composition.

Insight into the packing pattern of ß2 fibrils: a model study of glutamic acid rich oligomers with 13C isotopic edited vibrational spectroscopy.

Polyglutamic acid at low pH forms aggregates and self-assembles into a spiral, fibril-like superstructure formed as a ß2-type sheet conformation that has a more compact intersheet packing than commonly found. This is stabilized by three-centered bifurcated hydrogen bonding of the amide carbonyl involving the protonated glutamic acid side chain. We report vibrational spectroscopic results and analyses for oligopeptides rich in glutamic acid enhanced with (13)C isotope labeling in a study modeling low pH poly-Glu self-assembly. Our results indicate bifurcated H-bonding and ß2 aggregation can be attained in these model decamers, confirming they have the same conformations as poly-Glu. We also prepared conventional ß1-sheet aggregates by rapid precipitation from the residual peptides in the higher pH supernatant. By comparing the isotope-enhanced IR and VCD spectra with theoretical predictions, we deduced that the oligo-Glu ß2 structure is based on stacked, twisted, antiparallel ß-sheets. The best fit to theoretical predictions was obtained for the strands being out of register, sequentially stepped by one residue, in a ladder-like fashion. The alternate ß1 conformer for this oligopeptide was similarly shown to be antiparallel but was less ordered and apparently had a different registry in its aggregate structure.

Structural analyses of experimental 13C edited amide I' IR and VCD for peptide ß-sheet aggregates and fibrils using DFT-based spectral simulations.

In the preceding paper, computational models based on density functional theory (DFT) were presented to characterize the sensitivity of vibrational spectroscopic methods (IR, VCD, and Raman) to structural features of ß-sheets. Isotopically edited amide I' IR for peptides labeled with (13)C in multiple different sites provides the most structurally distinct signatures of strand alignment, while VCD is sensitive to the sheet twist and intersheet stacking. In this report, we simulate the IR and VCD spectra for models approximating structures of four ß-sheet forming peptides previously experimentally studied using these methods with (13)C isotopic editing. Various register alignments are tested. Agreement with experiment is evaluated based on frequency shifts of both the (12)C and (13)C IR amide I' signals, relative intensity patterns, and VCD spectra where available. While for the simulation of IR spectra canonical planar sheets provide a sufficient model system, for VCD simulation twisted, stacked sheets are required in order to reproduce strong couplet-like amide I' VCD. Effects of the solvent (water) and amino acid side chains are also tested by using a simplified, electrostatic solvent model and atomic partial charges for the side chains. Very good agreement with experimental spectra is obtained, particularly for the relative (12)C and (13)C band frequencies. All four peptide models are shown to be antiparallel as had previously been assumed. However, in some cases our simulations are consistent with different register alignment of strands than originally proposed.

Infrared, vibrational circular dichroism, and Raman spectral simulations for ß-sheet structures with various isotopic labels, interstrand, and stacking arrangements using density functional theory.

Infrared (IR), Raman, and vibrational circular dichroism (VCD) spectral variations for different ß-sheet structures were studied using simulations based on density functional theory (DFT) force field and intensity computations. The DFT vibrational parameters were obtained for ß-sheet fragments containing nine-amides and constrained to a variety of conformations and strand arrangements. These were subsequently transferred onto corresponding larger ß-sheet models, normally consisting of five strands with ten amides each, for spectral simulations. Further extension to fibril models composed of multiple stacked ß-sheets was achieved by combining the transfer of DFT parameters for each sheet with dipole coupling methods for interactions between sheets. IR spectra of the amide I show different splitting patterns for parallel and antiparallel ß-sheets, and their VCD, in the absence of intersheet stacking, have distinct sign variations. Isotopic labeling by (13)C of selected residues yields spectral shifts and intensity changes uniquely sensitive to relative alignment of strands (registry) for antiparallel sheets. Stacking of multiple planar sheets maintains the qualitative spectral character of the single sheet but evidences some reduction in the exciton splitting of the amide I mode. Rotating sheets with respect to each other leads to a significant VCD enhancement, whose sign pattern and intensity is dependent on the handedness and degree of rotation. For twisted ß-sheets, a significant VCD enhancement is computed even for sheets stacked with either the same or opposite alignments and the inter-sheet rotation, depending on the sense, can either further increase or weaken the enhanced VCD intensity. In twisted, stacked structures (without rotation), similar VCD amide I patterns (positive couplets) are predicted for both parallel and antiparallel sheets, but different IR intensity distributions still enable their differentiation. Our simulation results prove useful for interpreting experimental vibrational spectra in terms of ß-sheet and fibril structure, as illustrated in the accompanying paper.

DFT-based simulations of amide I' IR spectra of a small protein in solution using empirical electrostatic map with a continuum solvent model.

A continuum solvent model was tested for simulations of amide I' IR spectra for a 40-residue subdomain of P22 viral coat protein in aqueous solution. Spectra obtained using DFT (BPW91/6-31G**) parameters for a reduced all-Ala representation of the protein were corrected by an electrostatic potential map obtained from the solvent cavity surface and AMBER99 side-chain atom partial charges. Various cavity sizes derived from van der Waals atomic radii with an added effective solvent radius up to 2.0 Å were tested. The interplay of the side-chain and solvent electrostatic effects was investigated by considering the side chains and solvent separately as well as together. The sensitivity to side-chain conformational fluctuations and to the parametrization of C(ß) group partial charges was also tested. Simulation results were compared to the experimental amide I' spectra of P22 subdomain, including two (13)C isotopically edited variants, as well as to the previous simulations based on the molecular dynamics trajectory in explicit solvent. For small cavity sizes, between van der Waals and that with added solvent radius of 0.5 Å, better qualitative agreement with experiment was obtained than with the explicit solvent representation, in particular for the (13)C-labeled spectra. Larger protein cavities led to progressively worse predictions due to increasingly stronger electrostatic effects of side chains, which could no longer be well compensated for by the solvent potential. Balance between side-chain and solvent electrostatic effects is important in determining the width and shape of the simulated amide I', which is also virtually unaffected by side-chain-geometry fluctuations. The continuum solvent model combined with the electrostatic map is a computationally efficient and potentially robust approach for the simulations of IR spectra of proteins in solution.