Search

We describe a method for direct tRNA sequencing using the Oxford Nanopore MinION. The principal technical advance is custom adapters that facilitate end-to-end sequencing of individual transfer RNA (tRNA) molecules at subnanometer precision. A second advance is a nanopore sequencing pipeline optimized for tRNA. We tested this method using purified E. coli tRNA fMet , tRNA Lys , and tRNA Phe samples. 76-92% of individual aligned tRNA sequence reads were full length. As a proof of concept, we showed that nanopore sequencing detected all 43 expected isoacceptors in total E. coli MRE600 tRNA as well as isodecoders that further define that tRNA population. Alignment-based comparisons between the three purified tRNAs and their synthetic controls revealed systematic nucleotide miscalls that were diagnostic of known modifications. Systematic miscalls were also observed proximal to known modifications in total E. coli tRNA alignments, including a highly conserved pseudouridine in the T loop. This work highlights the potential of nanopore direct tRNA sequencing as well as improvements needed to implement tRNA sequencing for human healthcare applications.

Previously we showed that the protein unfoldase ClpX could facilitate translocation of individual proteins through the α-hemolysin nanopore. This results in ionic current fluctuations that correlate with unfolding and passage of intact protein strands through the pore lumen. It is plausible that this technology could be used to identify protein domains and structural modifications at the single-molecule level that arise from subtle changes in primary amino acid sequence (e.g., point mutations). As a test, we engineered proteins bearing well-characterized domains connected in series along an ∼700 amino acid strand. Point mutations in a titin immunoglobulin domain (titin I27) and point mutations, proteolytic cleavage, and rearrangement of beta-strands in green fluorescent protein (GFP), caused ionic current pattern changes for single strands predicted by bulk phase and force spectroscopy experiments. Among these variants, individual proteins could be classified at 86-99% accuracy using standard machine learning tools. We conclude that a ClpXP-nanopore device can discriminate among distinct protein domains, and that sequence-dependent variations within those domains are detectable.

Individual DNA molecules can be read at single nucleotide precision using nanopores coupled to processive enzymes. Discrimination among the four canonical bases has been achieved, as has discrimination among cytosine, 5-methylcytosine (mC), and 5-hydroxymethylcytosine (hmC). Two additional modified cytosine bases, 5-carboxylcytosine (caC) and 5-formylcytosine (fC), are produced during enzymatic conversion of hmC to cytosine in mammalian cells. Thus, an accurate picture of the cytosine epigenetic status in target cells should also include these C5-cytosine variants. In the present study, we used a patch clamp amplifier to acquire ionic current traces caused by phi29 DNA polymerase-controlled translocation of DNA templates through the M2MspA pore. Decision boundaries based on three consecutive ionic current states were implemented to call mC, hmC, caC, fC, or cytosine at CG dinucleotides in ∼4400 individual DNA molecules. We found that the percentage of correct base calls for single pass reads ranged from 91.6% to 98.3%. This accuracy depended upon the identity of nearest neighbor bases surrounding the CG dinucleotide.

Complexes formed between phi29 DNA polymerase (DNAP) and DNA fluctuate discretely between the pre-translocation and post-translocation states on the millisecond time scale. The translocation fluctuations can be observed in ionic current traces when individual complexes are captured atop the α-hemolysin nanopore in an electric field. The presence of complementary 2'-deoxynucleoside triphosphate (dNTP) shifts the equilibrium across the translocation step toward the post-translocation state. Here we have determined quantitatively the kinetic relationship between the phi29 DNAP translocation step and dNTP binding. We demonstrate that dNTP binds to phi29 DNAP-DNA complexes only after the transition from the pre-translocation state to the post-translocation state; dNTP binding rectifies the translocation but it does not directly drive the translocation. Based on the measured time traces of current amplitude, we developed a method for determining the forward and reverse translocation rates and the dNTP association and dissociation rates, individually at each dNTP concentration and each voltage. The translocation rates, and their response to force, match those determined for phi29 DNAP-DNA binary complexes and are unaffected by dNTP. The dNTP association and dissociation rates do not vary as a function of voltage, indicating that force does not distort the polymerase active site and that dNTP binding does not directly involve a displacement in the translocation direction. This combined experimental and theoretical approach and the results obtained provide a framework for separately evaluating the effects of biological variables on the translocation transitions and their effects on dNTP binding.

Complexes formed between the bacteriophage phi29 DNA polymerase (DNAP) and DNA fluctuate between the pre-translocation and post-translocation states on the millisecond time scale. These fluctuations can be directly observed with single-nucleotide precision in real-time ionic current traces when individual complexes are captured atop the α-hemolysin nanopore in an applied electric field. We recently quantified the equilibrium across the translocation step as a function of applied force (voltage), active-site proximal DNA sequences, and the binding of complementary dNTP. To gain insight into the mechanism of this step in the DNAP catalytic cycle, in this study, we have examined the stochastic dynamics of the translocation step. The survival probability of complexes in each of the two states decayed at a single exponential rate, indicating that the observed fluctuations are between two discrete states. We used a robust mathematical formulation based on the autocorrelation function to extract the forward and reverse rates of the transitions between the pre-translocation state and the post-translocation state from ionic current traces of captured phi29 DNAP-DNA binary complexes. We evaluated each transition rate as a function of applied voltage to examine the energy landscape of the phi29 DNAP translocation step. The analysis reveals that active-site proximal DNA sequences influence the depth of the pre-translocation and post-translocation state energy wells and affect the location of the transition state along the direction of the translocation.

Coupling nucleic acid processing enzymes to nanoscale pores allows controlled movement of individual DNA or RNA strands that is reported as an ionic current/time series. Hundreds of individual enzyme complexes can be examined in single-file order at high bandwidth and spatial resolution. The bacteriophage phi29 DNA polymerase (phi29 DNAP) is an attractive candidate for this technology, due to its remarkable processivity and high affinity for DNA substrates. Here we show that phi29 DNAP-DNA complexes are stable when captured in an electric field across the α-hemolysin nanopore. DNA substrates were activated for replication at the nanopore orifice by exploiting the 3'-5' exonuclease activity of wild-type phi29 DNAP to excise a 3'-H terminal residue, yielding a primer strand 3'-OH. In the presence of deoxynucleoside triphosphates, DNA synthesis was initiated, allowing real-time detection of numerous sequential nucleotide additions that was limited only by DNA template length. Translocation of phi29 DNAP along DNA substrates was observed in real time at Ångstrom-scale precision as the template strand was drawn through the nanopore lumen during replication.

By controlling the microarchitecture of bioengineered scaffolds for artificial tissues, their material and cell-interaction properties can be designed to mimic native correspondents. Current understanding of this relationship is sparse and based on microscopy requiring harsh sample preparation and labeling, leaving it open to which extent the natural morphology is studied. This work introduces multimodal nonlinear microscopy for label-free imaging of tissue scaffolds, exemplified by bacterial cellulose. Unique three-dimensional images visualizing the formation of nanofiber networks throughout the biosynthesis, revealing that supra-structures (layered structures, cavities) are formed. Cell integration in compact scaffolds was visualized and compared with porous scaffolds. While the former showed distinct boundaries to the native tissue, gradual cell integration was observed for the porous material. Thus, the degree of cell integration can be controlled through scaffold supra-structures. This illustrates the potential of nonlinear microscopy for noninvasive imaging of the intriguing interaction mechanisms between scaffolds and cells.

DNA polymerases are molecular motors that catalyze template-dependent DNA replication, advancing along template DNA by one nucleotide with each catalytic cycle. Nanopore-based measurements have emerged as a single molecule technique for the study of these enzymes. Using the alpha-hemolysin nanopore, we determined the position of DNA templates bearing inserts of abasic (1',2'-dideoxy) residues, bound to the Klenow fragment of Escherichia coli DNA polymerase I (KF) or to bacteriophage T7 DNA polymerase. Hundreds of individual polymerase complexes were analyzed at 5 A precision within minutes. We generated a map of current amplitudes for DNA-KF-deoxynucleoside triphosphate (dNTP) ternary complexes, using a series of templates bearing blocks of three abasic residues that were displaced by approximately 5 A in the nanopore lumen. Plotted as a function of the distance of the abasic insert from n = 0 in the active site of the enzyme held atop the pore, this map has a single peak. The map is similar when the primer length, the DNA sequences flanking the abasic insert, and the DNA sequences in the vicinity of the KF active site are varied. Primer extension catalyzed by KF using a three abasic template in the presence of a mixture of dNTPs and 2',3'-dideoxynucleoside triphosphates resulted in a ladder of ternary complexes with discrete amplitudes that closely corresponded to this map. An ionic current map measured in the presence of 0.15 M KCl mirrored the map obtained with 0.3 M KCl, permitting experiments with a broader range of mesophilic DNA and RNA processing enzymes. We used the abasic templates to show that capture of complexes with the KF homologue, T7 DNA polymerase, yields an amplitude map nearly indistinguishable from the KF map.

DNA polymerases catalyze template-dependent genome replication. The assembly of a high affinity ternary complex between these enzymes, the double strand-single strand junction of their DNA substrate, and the deoxynucleoside triphosphate (dNTP) complementary to the first template base in the polymerase active site is essential to this process. We present a single molecule method for iterative measurements of DNA-polymerase complex assembly with high temporal resolution, using active voltage control of individual DNA substrate molecules tethered noncovalently in an alpha-hemolysin nanopore. DNA binding states of the Klenow fragment of Escherichia coli DNA polymerase I (KF) were diagnosed based upon their ionic current signature, and reacted to with submillisecond precision to execute voltage changes that controlled exposure of the DNA substrate to KF and dNTP. Precise control of exposure times allowed measurements of DNA-KF complex assembly on a time scale that superimposed with the rate of KF binding. Hundreds of measurements were made with a single tethered DNA molecule within seconds, and dozens of molecules can be tethered within a single experiment. This approach allows statistically robust analysis of the assembly of complexes between DNA and RNA processing enzymes and their substrates at the single molecule level.

Nanoscale pores are a tool for single molecule analysis of DNA or RNA processing enzymes. Monitoring catalytic activity in real time using this technique requires that these enzymes retain function while held atop a nanopore in an applied electric field. Using an alpha-hemolysin nanopore, we measured the dwell time for complexes of DNA with the Klenow fragment of Escherichia coli DNA polymerase I (KF) as a function of the concentration of deoxynucleoside triphosphate (dNTP) substrate. We analyzed these dwell time measurements in the framework of a two-state model for captured complexes (DNA-KF binary and DNA-KF-dNTP ternary states). Average nanopore dwell time increased without saturating as a function of correct dNTP concentration across 4 orders of magnitude. This arises from two factors that are proportional to dNTP concentration: (1) The fraction of complexes that are in the ternary state when initially captured predominantly affects dwell time at low dNTP concentrations. (2) The rate of binding and rebinding of dNTP to captured complexes affects dwell time at higher dNTP concentrations. Thus there are two regimes that display a linear relationship between average dwell time and dNTP concentration. The transition from one linear regime to the other occurs near the equilibrium dissociation constant (K(d)) for dNTP binding to KF-DNA complexes in solution. We conclude from the combination of titration experiments and modeling that DNA-KF complexes captured atop the nanopore retain iterative, sequence-specific dNTP binding, as required for catalysis and fidelity in DNA synthesis.

The mammalian bombesin receptor subfamily of G protein-coupled receptors currently consists of the gastrin-releasing peptide receptor (GRP-R), neuromedin B receptor, and bombesin receptor subtype 3. All three receptors contain a conserved aspartate residue (D98) at the extracellular boundary of transmembrane domain II and a conserved arginine residue (R309) near the extracellular boundary of transmembrane domain VII. To evaluate the functional role of these residues, site-directed GRP-R mutants were expressed in fibroblasts and assayed for their ability to both bind agonist and catalyze exchange of guanine nucleotides. Alanine substitution at GRP-R position 98 or 309 reduced agonist binding affinity by 24- and 56-fold, respectively, compared to wild-type GRP-R. Single swap GRP-R mutations either resulted in no receptor expression in the membrane (D98R) or the protein was not able to bind agonist (R309D). In contrast, the double swap mutation (D98R/R309D) had high-affinity agonist binding, reduced from wild-type GRP-R by only 6-fold. In situ reconstitution of urea-extracted membranes expressing either wild-type or mutant (D98A or R309A) GRP-R with G(q) indicated that alanine substitution greatly reduced G protein catalytic exchange compared to wild-type GRP-R. The D98R/R309D GRP-R had both a higher intrinsic basal activity and a higher overall catalytic exchange activity compared to wild-type; however, the wild-type GRP-R produced a larger agonist-stimulated response relative to the double swap mutant. Taken together, these data show that GRP-R residues D98 and R309 are critical for efficient coupling of GRP-R to G(q). Furthermore, our findings are consistent with a salt bridge interaction between these two polar and oppositely charged amino acids that maintains the proper receptor conformation necessary to interact with G proteins.

Recently, a fourth member of the bombesin (Bn) receptor family (fBB4-R) was isolated from a cDNA library from the brain of the frog, Bombina orientalis. Its pharmacology and cell biology are largely unknown, and no known natural cell lines or tissues possess sufficient numbers of fBB4-R's to allow either of these to be determined. To address these issues, we have used three different strategies. fBB4-R expression in cells widely used for other Bn receptor subtypes was unsuccessful as was expression in two frog cell lines. However, stable fBB4-R cell lines were obtained in CHO-K1 cells which were shown to faithfully demonstrate the correct pharmacology of the related Bn receptor, the GRP receptor, when expressed in these cells. [DPhe6,betaAla11,Phe13,Nle14]Bn(6-14) was found to have high affinity (Ki = 0.4 nM) for the fBB4 receptor and 125I-[DTyr6,betaala11,Phe13,Nle14]Bn(6-14) to be an excellent ligand for this receptor. The fBB4-R had a unique pharmacology for naturally occurring Bn-related agonists, with the presence of a penultimate phenylalanine being critical for high-affinity interaction. It also had a unique profile for six classes of Bn antagonists. The fBB4-R was coupled to phospholipase C with activation increasing [3H]inositol phosphates and mobilizing Ca2+ almost entirely from cellular sources. There was a close correlation between agonist the receptor occupation and the receptor activation. Three of the five classes of Bn receptor antagonists that interacted with higher affinity with the fBB4-R functioned as fBB4-R antagonists and two as partial agonists. fBB4-R activation stimulated increases in phospholipase D (PLD) over the same range of concentrations at which it activated phospholipase C. These results demonstrate that the fBB4 receptor has a unique pharmacology for agonists and antagonists and is coupled to phospholipase C and D. The availability of these cell lines, this novel ligand, and the identification of three classes of antagonists that can be used as lead compounds should facilitate the further investigation of the pharmacology and cell biology of the BB4 receptor.

The molecular mechanism of general anesthesia is not understood. Possible modes of action include binding at a protein site, such as a receptor or channel, or physical effects on membrane lipid properties. The pump-leak hypothesis suggests that anesthetics perturb the bilayer of synaptic vesicles, thereby increasing ionic permeability. This results in decay of proton gradients required for transport and accumulation of neurotransmitters. The subsequent loss of neurotransmitters from synaptic vesicles reduces the efficiency of synaptic transmission and results in the anesthetized state. We have determined the effects of general anesthetics on certain parameters of enzyme activity and membrane permeability relevant to the pump-leak hypothesis. We used chromaffin granules as a convenient model system and focused on clinically relevant anesthetic concentrations (ED50), quantitative measurements of permeability changes, and the kinetics of gradient decay. General anesthetics at ED50 have little or no effect on the proton-transport ATPase activity, but do cause modest increments in proton permeability that change the catecholamine distribution in actively pumping chromaffin granule preparations. We found that pH gradients do not collapse entirely under these conditions and that only a fraction of total catecholamine is lost from the chromaffin granules. When total collapse is induced by other means, efflux of catecholamines occurs with a half-time near 30 min. These results suggest that if the pump-leak hypothesis is valid, then very small losses of catecholamines must be sufficient to induce anesthesia. We conclude that the weight of evidence favors other mechanisms, notably direct binding of anesthetics to sensitive proteins.