Mechanical and chemical unfolding of a single protein: A comparison

Is the mechanical unraveling of protein domains by atomic force microscopy (AFM) just a technological feat or a true measurement of their unfolding? By engineering a protein made of tandem repeats of identical Ig modules, we were able to get explicit AFM data on the unfolding rate of a single protein domain that can be accurately extrapolated to zero force. We compare this with chemical unfolding rates for untethered modules extrapolated to 0 M denaturant. The unfolding rates obtained by the two methods are the same. Furthermore, the transition state for unfolding appears at the same position on the folding pathway when assessed by either method. These results indicate that mechanical unfolding of a single protein by AFM does indeed reflect the same event that is observed in traditional unfolding experiments. The way is now open for the extensive use of AFM to measure folding reactions at the single-molecule level. Single-molecule AFM recordings have the added advantage that they define the reaction coordinate and expose rare unfolding events that cannot be observed in the absence of chemical denaturants.

Force-mediated kinetics of single P-selectin/ligand complexes observed by atomic force microscopy

Leukocytes roll along the endothelium of postcapillary venules in response to inflammatory signals. Rolling under the hydrodynamic drag forces of blood flow is mediated by the interaction between selectins and their ligands across the leukocyte and endothelial cell surfaces. Here we present force-spectroscopy experiments on single complexes of P-selectin and P-selectin glycoprotein ligand-1 by atomic force microscopy to determine the intrinsic molecular properties of this dynamic adhesion process. By modeling intermolecular and intramolecular forces as well as the adhesion probability in atomic force microscopy experiments we gain information on rupture forces, elasticity, and kinetics of the P-selectin/P-selectin glycoprotein ligand-1 interaction. The complexes are able to withstand forces up to 165 pN and show a chain-like elasticity with a molecular spring constant of 5.3 pN nm−1 and a persistence length of 0.35 nm. The dissociation constant (off-rate) varies over three orders of magnitude from 0.02 s−1 under zero force up to 15 s−1 under external applied forces. Rupture force and lifetime of the complexes are not constant, but directly depend on the applied force per unit time, which is a product of the intrinsic molecular elasticity and the external pulling velocity. The high strength of binding combined with force-dependent rate constants and high molecular elasticity are tailored to support physiological leukocyte rolling.

Conformational fluctuations in single DNA molecules

Measurement of fluorescent lifetimes of dye-tagged DNA molecules reveal the existence of different conformations. Conformational fluctuations observed by fluorescence correlation spectroscopy give rise to a relaxation behavior that is described by “stretched” exponentials and indicates the presence of a distribution of transition rates between two conformations. Whether this is an inhomogeneous distribution, where each molecule contributes with its own reaction rate to the overall distribution, or a homogeneous distribution, where the reaction rate of each molecule is time-dependent, is not yet known. We used a tetramethylrhodamine-linked 217-bp DNA oligonucleotide as a probe for conformational fluctuations. Fluorescence fluctuations from single DNA molecules attached to a streptavidin-coated surface directly show the transitions between two conformational states. The conformational fluctuations typical for single molecules are similar to those seen in single ion channels in cell membranes.

Evolution of plastid gene rps2 in a lineage of hemiparasitic and holoparasitic plants: Many losses of photosynthesis and complex patterns of rate variation

The plastid genomes of some nonphotosynthetic parasitic plants have experienced an extreme reduction in gene content and an increase in evolutionary rate of remaining genes. Nothing is known of the dynamics of these events or whether either is a direct outcome of the loss of photosynthesis. The parasitic Scrophulariaceae and Orobanchaceae, representing a continuum of heterotrophic ability ranging from photosynthetic hemiparasites to nonphotosynthetic holoparasites, are used to investigate these issues. We present a phylogenetic hypothesis for parasitic Scrophulariaceae and Orobanchaceae based on sequences of the plastid gene rps2, encoding the S2 subunit of the plastid ribosome. Parasitic Scrophulariaceae and Orobanchaceae form a monophyletic group in which parasitism can be inferred to have evolved once. Holoparasitism has evolved independently at least five times, with certain holoparasitic lineages representing single species, genera, and collections of nonphotosynthetic genera. Evolutionary loss of the photosynthetic gene rbcL is limited to a subset of holoparasitic lineages, with several holoparasites retaining a full length rbcL sequence. In contrast, the translational gene rps2 is retained in all plants investigated but has experienced rate accelerations in several hemi- as well as holoparasitic lineages, suggesting that there may be substantial molecular evolutionary changes to the plastid genome of parasites before the loss of photosynthesis. Independent patterns of synonymous and nonsynonymous rate acceleration in rps2 point to distinct mechanisms underlying rate variation in different lineages. Parasitic Scrophulariaceae (including the traditional Orobanchaceae) provide a rich platform for the investigation of molecular evolutionary process, gene function, and the evolution of parasitism.

Limited internal friction in the rate-limiting step of a two-state protein folding reaction

Small, single-domain proteins typically fold via a compact transition-state ensemble in a process well fitted by a simple, two-state model. To characterize the rate-limiting conformational changes that underlie two-state folding, we have investigated experimentally the effects of changing solvent viscosity on the refolding of the IgG binding domain of protein L. In conjunction with numerical simulations, our results indicate that the rate-limiting conformational changes of the folding of this domain are strongly coupled to solvent viscosity and lack any significant “internal friction” arising from intrachain collisions. When compared with the previously determined solvent viscosity dependencies of other, more restricted conformational changes, our results suggest that the rate-limiting folding transition involves conformational fluctuations that displace considerable amounts of solvent. Reconciling evidence that the folding transition state ensemble is comprised of highly collapsed species with these and similar, previously reported results should provide a significant constraint for theoretical models of the folding process.

Inverse radiation dose-rate effects on somatic and germ-line mutations and DNA damage rates

The mutagenic effect of low linear energy transfer ionizing radiation is reduced for a given dose as the dose rate (DR) is reduced to a low level, a phenomenon known as the direct DR effect. Our reanalysis of published data shows that for both somatic and germ-line mutations there is an opposite, inverse DR effect, with reduction from low to very low DR, the overall dependence of induced mutations being parabolically related to DR, with a minimum in the range of 0.1 to 1.0 cGy/min (rule 1). This general pattern can be attributed to an optimal induction of error-free DNA repair in a DR region of minimal mutability (MMDR region). The diminished activation of repair at very low DRs may reflect a low ratio of induced (“signal”) to spontaneous background DNA damage (“noise”). Because two common DNA lesions, 8-oxoguanine and thymine glycol, were already known to activate repair in irradiated mammalian cells, we estimated how their rates of production are altered upon radiation exposure in the MMDR region. For these and other abundant lesions (abasic sites and single-strand breaks), the DNA damage rate increment in the MMDR region is in the range of 10% to 100% (rule 2). These estimates suggest a genetically programmed optimatization of response to radiation in the MMDR region.

Determination of the lifetime and force dependence of interactions of single bonds between surface-attached CD2 and CD48 adhesion molecules

We studied single molecular interactions between surface-attached rat CD2, a T-lymphocyte adhesion receptor, and CD48, a CD2 ligand found on antigen-presenting cells. Spherical particles were coated with decreasing densities of CD48–CD4 chimeric molecules then driven along CD2-derivatized glass surfaces under a low hydrodynamic shear rate. Particles exhibited multiple arrests of varying duration. By analyzing the dependence of arrest frequency and duration on the surface density of CD48 sites, it was concluded that (i) arrests were generated by single molecular bonds and (ii) the initial bond dissociation rate was about 7.8 s−1. The force exerted on bonds was increased from about 11 to 22 pN; the detachment rate exhibited a twofold increase. These results agree with and extend studies on the CD2–CD48 interaction by surface plasmon resonance technology, which yielded an affinity constant of ≈104 M−1 and a dissociation rate of ≥6 s−1. It is concluded that the flow chamber technology can be an useful complement to atomic force microscopy for studying interactions between isolated biomolecules, with a resolution of about 20 ms and sensitivity of a few piconewtons. Further, this technology might be extended to actual cells.

Mutation rate varies among alleles at a microsatellite locus: Phylogenetic evidence

The understanding of the mutational mechanism that generates high levels of variation at microsatellite loci lags far behind the application of these genetic markers. A phylogenetic approach was developed to study the pattern and rate of mutations at a dinucleotide microsatellite locus tightly linked to HLA-DQB1 (DQCAR). A random Japanese population (n = 129) and a collection of multiethnic samples (n = 941) were typed at the DQB1 and DQCAR loci. The phylogeny of DQB1 alleles was then reconstructed and DQCAR alleles were superimposed onto the phylogeny. This approach allowed us to group DQCAR alleles that share a common ancestor. The results indicated that the DQCAR mutation rate varies drastically among alleles within this single microsatellite locus. Some DQCAR alleles never mutated during a long period of evolutionary time. Sequencing of representative DQCAR alleles showed that these alleles lost their ability to mutate because of nucleotide substitutions that shorten the length of uninterrupted CA repeat arrays; in contrast, all mutating alleles had relatively longer perfect CA repeat sequences.

Stretch-activated single K+ channels account for whole-cell currents elicited by swelling

Functionally significant stretch-activated ion channels have been clearly identified in excitable cells. Although single-channel studies suggest their expression in other cell types, their activity in the whole-cell configuration has not been shown. This discrepancy makes their physiological significance doubtful and suggests that their mechanical activation is artifactual. Possible roles for these molecules in nonexcitable cells are acute cell-volume regulation and, in epithelial cells, the complex adjustment of ion fluxes across individual cell membranes when the rate of transepithelial transport changes. We report the results of experiments on isolated epithelial cells expressing in the basolateral membrane stretch-activated K+ channels demonstrable by the cell-attached patch-clamp technique. In these cells, reversible whole-cell currents were elicited by both isosmotic and hyposmotic cell swelling. Cation selectivity and block by inorganic agents were the same for single-channel and whole-cell currents, indicating that the same entity underlies single-channel and whole-cell currents and that the single-channel events are not artifactual. In these cells, when the rate of apical-membrane NaCl entry increases, the cell Na+ content and volume also increase, stimulating the Na+,K+-ATPase at the basolateral membrane, i.e., both Na+ extrusion and K+ uptake increase. We speculate that, under these conditions, the parallel activation of basolateral K+ channels (by the swelling) elevates conductive K+ loss, tending to maintain the cell K+ content constant (“pump-leak parallelism”). This study describes a physiologically relevant stretch-activated channel, at both the single-channel and whole-cell levels, in a nonneural cell type.

Unbinding forces of single antibody-antigen complexes correlate with their thermal dissociation rates

Point mutants of three unrelated antifluorescein antibodies were constructed to obtain nine different single-chain Fv fragments, whose on-rates, off-rates, and equilibrium binding affinities were determined in solution. Additionally, activation energies for unbinding were estimated from the temperature dependence of the off-rate in solution. Loading rate-dependent unbinding forces were determined for single molecules by atomic force microscopy, which extrapolated at zero force to a value close to the off-rate measured in solution, without any indication for multiple transition states. The measured unbinding forces of all nine mutants correlated well with the off-rate in solution, but not with the temperature dependence of the reaction, indicating that the same transition state must be crossed in spontaneous and forced unbinding and that the unbinding path under load cannot be too different from the one at zero force. The distance of the transition state from the ground state along the unbinding pathway is directly proportional to the barrier height, regardless of the details of the binding site, which most likely reflects the elasticity of the protein in the unbinding process. Atomic force microscopy thus can be a valuable tool for the characterization of solution properties of protein-ligand systems at the single molecule level, predicting relative off-rates, potentially of great value for combinatorial chemistry and biology.

Involvement of a cytosine side chain in proton transfer in the rate-determining step of ribozyme self-cleavage

Ribozymes of hepatitis delta virus have been proposed to use an active-site cytosine as an acid-base catalyst in the self-cleavage reaction. In this study, we have examined the role of cytosine in more detail with the antigenomic ribozyme. Evidence that proton transfer in the rate-determining step involved cytosine 76 (C76) was obtained from examining cleavage activity of the wild-type and imidazole buffer-rescued C76-deleted (C76Δ) ribozymes in D2O and H2O. In both reactions, a similar kinetic isotope effect and shift in the apparent pKa indicate that the buffer is functionally substituting for the side chain in proton transfer. Proton inventory of the wild-type reaction supported a mechanism of a single proton transfer at the transition state. This proton transfer step was further characterized by exogenous base rescue of a C76Δ mutant with cytosine and imidazole analogues. For the imidazole analogues that rescued activity, the apparent pKa of the rescue reaction, measured under kcat/KM conditions, correlated with the pKa of the base. From these data a Brønsted coefficient (β) of 0.51 was determined for the base-rescued reaction of C76Δ. This value is consistent with that expected for proton transfer in the transition state. Together, these data provide strong support for a mechanism where an RNA side chain participates directly in general acid or general base catalysis of the wild-type ribozyme to facilitate RNA cleavage.

Chemically distinct transition states govern rapid dissociation of single L-selectin bonds under force

Carbohydrate–protein bonds interrupt the rapid flow of leukocytes in the circulation by initiation of rolling and tethering at vessel walls. The cell surface carbohydrate ligands are glycosylated proteins like the mucin P-selectin glycoprotein ligand-1 (PSGL-1), which bind ubiquitously to the family of E-, P-, and L-selectin proteins in membranes of leukocytes and endothelium. The current view is that carbohydrate–selectin bonds dissociate a few times per second, and the unbinding rate increases weakly with force. However, such studies have provided little insight into how numerous hydrogen bonds, a Ca2+ metal ion bond, and other interactions contribute to the mechanical strength of these attachments. Decorating a force probe with very dilute ligands and controlling touch to achieve rare single-bond events, we have varied the unbinding rates of carbohydrate–selectin bonds by detachment with ramps of force/time from 10 to 100,000 pN/sec. Testing PSGL-1, its outer 19 aa (19FT), and sialyl LewisX (sLeX) against L-selectin in vitro on glass microspheres and in situ on neutrophils, we found that the unbinding rates followed the same dependence on force and increased by nearly 1,000-fold as rupture forces rose from a few to ≈200 pN. Plotted on a logarithmic scale of loading rate, the rupture forces reveal two prominent energy barriers along the unbinding pathway. Strengths above 75 pN arise from rapid detachment (<0.01 sec) impeded by an inner barrier that requires a Ca2+ bond between a single sLeX and the lectin domain. Strengths below 75 pN occur under slow detachment (>0.01 sec) impeded by the outer barrier, which appears to involve an array of weak (putatively hydrogen) bonds.

Tuning DNA “strings”: Modulating the rate of DNA replication with mechanical tension

Recent experiments have measured the rate of replication of DNA catalyzed by a single enzyme moving along a stretched template strand. The dependence on tension was interpreted as evidence that T7 and related DNA polymerases convert two (n = 2) or more single-stranded template bases to double helix geometry in the polymerization site during each catalytic cycle. However, we find structural data on the T7 enzyme–template complex indicate n = 1. We also present a model for the “tuning” of replication rate by mechanical tension. This model considers only local interactions in the neighborhood of the enzyme, unlike previous models that use stretching curves for the entire polymer chain. Our results, with n = 1, reconcile force-dependent replication rate studies with structural data on DNA polymerase complexes.

Ribonuclease P (RNase P) RNA is converted to a Cd(2+)-ribozyme by a single Rp-phosphorothioate modification in the precursor tRNA at the RNase P cleavage site.

To study the cleavage mechanism of bacterial Nase P RNA, we have synthesized precursor tRNA substrates carrying a single Rp- or Sp-phosphorothioate modification at the RNase P cleavage site. Both the Sp- and the Rp-diastereomer reduced the rate of processing by Escherichia coli RNase P RNA at least 1000-fold under conditions where the chemical step is rate-limiting. The Rp-modification had no effect and the Sp-modification had a moderate effect on precursor tRNA ground state binding to RNase P RNA. Processing of the Rp-diastereomeric substrate was largely restored in the presence of the "thiophilic" Cd2+ as the only divalent metal ion, demonstrating direct metal ion coordination to the (pro)-Rp substituent at the cleavage site and arguing against a specific role for Mg(2+)-ions at the pro-Sp oxygen. For the Rp-diastereomeric substrate, Hill plot analysis revealed a cooperative dependence upon [Cd2+] of nH = 1.8, consistent with a two-metal ion mechanism. In the presence of the Sp-modification, neither Mn2+ nor Cd2+ was able to restore detectable cleavage at the canonical site. Instead, the ribozyme promotes cleavage at the neighboring unmodified phosphodiester with low efficiency. Dramatic inhibition of the chemical step by both the Rp- and Sp-phosphorothioate modification is unprecedented among known ribozymes and points to unique features of transition state geometry in the RNase P RNA-catalyzed reaction.

Replication infidelity during a single cycle of Ty1 retrotransposition.

Retroviruses undergo a high frequency of genetic alterations during the process of copying their RNA genomes. However, little is known about the replication fidelity of other elements that transpose via reverse transcription of an RNA intermediate. The complete sequence of 29 independently integrated copies of the yeast retrotransposon Ty1 (173,043 nt) was determined, and the mutation rate during a single cycle of replication was calculated. The observed base substitution rate of 2.5 x 10(-5) bp per replication cycle suggests that this intracellular element can mutate as rapidly as retroviruses. The pattern and distribution of errors in the Ty1 genome is nonrandom and provides clues to potential in vivo molecular mechanisms of reverse transcriptase-mediated error generation, including heterogeneous RNase H cleavage of Ty1 RNA, addition of terminal nontemplated bases, and transient dislocation and realignment of primer-templates. Overall, analysis of errors generated during Ty1 replication underscores the utility of a genetically tractable model system for the study of reverse transcriptase fidelity.