PubMed Health⌕ Search

Biomedical subjects

Tobin R Sosnick

Publications and source records attributed to Tobin R Sosnick.

At least 19 recordsLinked to original sources
PubMed2007

Polypeptide motions are dominated by peptide group oscillations resulting from dihedral angle correlations between nearest neighbors.

To identify basic local backbone motions in unfolded chains, simulations are performed for a variety of peptide systems using three popular force fields and for implicit and explicit solvent models. A dominant "crankshaft-like" motion is found that involves only a localized oscillation of the plane of the peptide group. This motion results in a strong anticorrelated motion of the phi angle of the ith residue (phi(i)) and the psi angle of the residue i - 1 (psi(i-1)) on the 0.1 ps time scale. Only a slight correlation is found between the motions of the two backbone dihedral angles of the same residue. Aside from the special cases of glycine and proline, no correlations are found between backbone dihedral angles that are separated by more than one torsion angle. These short time, correlated motions are found both in equilibrium fluctuations and during the transit process between Ramachandran basins, e.g., from the beta to the alpha region. A residue's complete transit from one Ramachandran basin to another, however, occurs in a manner independent of its neighbors' conformational transitions. These properties appear to be intrinsic because they are robust across different force fields, solvent models, nonbonded interaction routines, and most amino acids.

Alanine↗
PubMed2007

Characterizing protein folding transition States using Psi-analysis.

We discuss the implementation of Psi-analysis for the structural characterization of protein folding transition states. In Psi-analysis, engineered bi-histidine metal ion binding sites are introduced at surface positions to stabilize secondary and tertiary structures. The addition of metal ions stabilizes the interaction between the two known histidines in a continuous fashion. Measuring the ratio of transition state stabilization to that of the native state provides information about the presence of the metal binding site in the transition state. Psi-Analysis uses noninvasive surface mutations and does not require specialized equipment, so it can be readily applied to characterize the folding of many proteins. As a result, this method can provide a wealth of high-resolution quantitative data for comparison with theoretical folding simulations. Additionally, investigations of other biological processes also may utilize metal binding sites and Psi-analysis to detect conformational events during catalysis, assembly, and function.

Binding Sites↗
PubMed2006

Minimalist representations and the importance of nearest neighbor effects in protein folding simulations.

In order to investigate the level of representation required to simulate folding and predict structure, we test the ability of a variety of reduced representations to identify native states in decoy libraries and to recover the native structure given the advanced knowledge of the very broad native Ramachandran basin assignments. Simplifications include the removal of the entire side-chain or the retention of only the Cbeta atoms. Scoring functions are derived from an all-atom statistical potential that distinguishes between atoms and different residue types. Structures are obtained by minimizing the scoring function with a computationally rapid simulated annealing algorithm. Results are compared for simulations in which backbone conformations are sampled from a Protein Data Bank-based backbone rotamer library generated by either ignoring or including a dependence on the identity and conformation of the neighboring residues. Only when the Cbeta atoms and nearest neighbor effects are included do the lowest energy structures generally fall within 4 A of the native backbone root-mean square deviation (RMSD), despite the initial configuration being highly expanded with an average RMSD > or = 10 A. The side-chains are reinserted into the Cbeta models with minimal steric clash. Therefore, the detailed, all-atom information lost in descending to a Cbeta-level representation is recaptured to a large measure using backbone dihedral angle sampling that includes nearest neighbor effects and an appropriate scoring function.

Amino Acids↗
PubMed2006

Small proteins fold through transition states with native-like topologies.

The folding pathway of common-type acyl phosphatase (ctAcP) is characterized using psi-analysis, which identifies specific chain-chain contacts using bi-histidine (biHis) metal-ion binding sites. In the transition state ensemble (TSE), the majority of the protein is structured with a near-native topology, only lacking one beta-strand and an alpha-helix. psi-Values are zero or unity for all sites except one at the amino terminus of helix H2. This fractional psi-value remains unchanged when three metal ions of differing coordination geometries are used, indicating this end of the helix experiences microscopic heterogeneity through fraying in the TSE. Ubiquitin, the other globular protein characterized using psi-analysis, also exhibits a single consensus TSE structure. Hence, the TSE of both proteins have converged to a single configuration, albeit one that contains some fraying at the periphery. Models of the TSE of both proteins are created using all-atom Langevin dynamics simulations using distance constraints derived from the experimental psi-values. For both proteins, the relative contact order of the TS models is approximately 80% of the native value. This shared value viewed in the context of the known correlation between contact order and folding rates, suggests that other proteins will have a similarly high fraction of the native contact order. This constraint greatly limits the range of possible configurations at the rate-limiting step.

Computer Simulation↗
PubMed2006

PII structure in the model peptides for unfolded proteins: studies on ubiquitin fragments and several alanine-rich peptides containing QQQ, SSS, FFF, and VVV.

A great deal of attention has been paid lately to the structures in unfolded proteins due to the recent discovery of many biologically functional but natively unfolded proteins and the far-reaching implications of order in unfolded states for protein folding. Recently, studies on oligo-Ala, oligo-Lys, oligo-Asp, and oligo-Glu, as well as oligo-Pro, have indicated that the left-handed polyproline II (PII) is the major local structure in these short peptides. Here, we show by NMR and CD studies that ubiquitin fragments, model unfolded peptides composed of nonrepeating amino acids, and four alanine-rich peptides containing QQQ, SSS, FFF, and VVV sequences are all present in aqueous solution predominantly in the extended PII or beta conformation. The results from this and related studies indicate that PII might be a major backbone conformation in unfolded proteins. The presence of defined local backbone structure in unfolded proteins is inconsistent with predictions from random coil models.

Alanine↗
PubMed2006

On the precision of experimentally determined protein folding rates and phi-values.

Phi-values, a relatively direct probe of transition-state structure, are an important benchmark in both experimental and theoretical studies of protein folding. Recently, however, significant controversy has emerged regarding the reliability with which phi-values can be determined experimentally: Because phi is a ratio of differences between experimental observables it is extremely sensitive to errors in those observations when the differences are small. Here we address this issue directly by performing blind, replicate measurements in three laboratories. By monitoring within- and between-laboratory variability, we have determined the precision with which folding rates and phi-values are measured using generally accepted laboratory practices and under conditions typical of our laboratories. We find that, unless the change in free energy associated with the probing mutation is quite large, the precision of phi-values is relatively poor when determined using rates extrapolated to the absence of denaturant. In contrast, when we employ rates estimated at nonzero denaturant concentrations or assume that the slopes of the chevron arms (mf and mu) are invariant upon mutation, the precision of our estimates of phi is significantly improved. Nevertheless, the reproducibility we thus obtain still compares poorly with the confidence intervals typically reported in the literature. This discrepancy appears to arise due to differences in how precision is calculated, the dependence of precision on the number of data points employed in defining a chevron, and interlaboratory sources of variability that may have been largely ignored in the prior literature.

Fluorometry↗
PubMed2006

Methods for the accurate estimation of confidence intervals on protein folding phi-values.

Phi-values provide an important benchmark for the comparison of experimental protein folding studies to computer simulations and theories of the folding process. Despite the growing importance of phi measurements, however, formulas to quantify the precision with which phi is measured have seen little significant discussion. Moreover, a commonly employed method for the determination of standard errors on phi estimates assumes that estimates of the changes in free energy of the transition and folded states are independent. Here we demonstrate that this assumption is usually incorrect and that this typically leads to the underestimation of phi precision. We derive an analytical expression for the precision of phi estimates (assuming linear chevron behavior) that explicitly takes this dependence into account. We also describe an alternative method that implicitly corrects for the effect. By simulating experimental chevron data, we show that both methods accurately estimate phi confidence intervals. We also explore the effects of the commonly employed techniques of calculating phi from kinetics estimated at non-zero denaturant concentrations and via the assumption of parallel chevron arms. We find that these approaches can produce significantly different estimates for phi (again, even for truly linear chevron behavior), indicating that they are not equivalent, interchangeable measures of transition state structure. Lastly, we describe a Web-based implementation of the above algorithms for general use by the protein folding community.

Algorithms↗
PubMed2006

Structural basis for altering the stability of homologous RNAs from a mesophilic and a thermophilic bacterium.

Tertiary RNA structures from thermophilic bacteria generally are more stable than their mesophilic homologs. To understand the structural basis of the increase in stability, we investigated equilibrium folding of the specificity domain (S-domain) of RNase P RNA from a mesophilic (Escherichia coli) and a thermophilic (Thermus thermophilus) bacterium. Equilibrium folding of both S-domains is described by a minimal, three-state folding scheme, U-to-I-to-N. In the I-to-N transition of the thermophilic S-domain, more structure forms and protections are stronger against T1 nuclease and hydroxyl radical reactions. Phylogenetic comparison in the context of the native structure reveals that among 39 nucleotide differences between these S-domains, 12 likely contribute to higher stability. These residues participate in extensive networks of hydrogen bonding, stacking, and metal ion coordination throughout the molecule. The thermophilic S-domain achieves higher stability by mutating strategic base pairs to G-C, decreasing surface accessibility of the native state, and increasing the amount of structure formation in the native folding transition. An E. coli S-domain mutant containing these 12 nt has the same stability and folding cooperativity as the T. thermophilus S-domain. E. coli S-domain mutants containing a subset of 4 or 6 nt have the same stability as the T. thermophilus S-domain but the same folding cooperativity as the E. coli S-domain. These results show that increasing stability can be accomplished by mutations within a local structure, but increasing folding cooperativity needs concerted changes among multiple structural units.

Base Sequence↗
PubMed2005

Structure of a folding intermediate reveals the interplay between core and peripheral elements in RNA folding.

Though the molecular architecture of many native RNA structures has been characterized, the structures of folding intermediates are poorly defined. Here, we present a nucleotide-level model of a highly structured equilibrium folding intermediate of the specificity domain of the Bacillus subtilis RNase P RNA, obtained using chemical and nuclease mapping, circular dichroism spectroscopy, small-angle X-ray scattering and molecular modeling. The crystal structure indicates that the 154 nucleotide specificity domain is composed of several secondary and tertiary structural modules. The structure of the intermediate contains modules composed of secondary structures and short-range tertiary interactions, implying a sequential order of tertiary structure formation during folding. The intermediate lacks the native core and several long-range interactions among peripheral regions, such as a GAAA tetraloop and its receptor. Folding to the native structure requires the local rearrangement of a T-loop in the core in concert with the formation of the GAAA tetraloop-receptor interaction. The interplay of core and peripheral structure formation rationalizes the high degree of cooperativity observed in the folding transition leading to the native structure.

Bacillus subtilis↗
PubMed2005

Statistical coil model of the unfolded state: resolving the reconciliation problem.

An unfolded state ensemble is generated by using a self-avoiding statistical coil model that is based on backbone conformational frequencies in a coil library, a subset of the Protein Data Bank. The model reproduces two apparently contradicting behaviors observed in the chemically denatured state for a variety of proteins, random coil scaling of the radius of gyration and the presence of significant amounts of local backbone structure (NMR residual dipolar couplings). The most stretched members of our unfolded ensemble dominate the residual dipolar coupling signal, whereas the uniformity of the sign of the couplings follows from the preponderance of polyproline II and beta conformers in the coil library. Agreement with the NMR data substantially improves when the backbone conformational preferences include correlations arising from the chemical and conformational identity of neighboring residues. Although the unfolded ensembles match the experimental observables, they do not display evidence of native-like topology. By providing an accurate representation of the unfolded state, our statistical coil model can be used to improve thermodynamic and kinetic modeling of protein folding.

Computational Biology↗
PubMed2005

Helix, sheet, and polyproline II frequencies and strong nearest neighbor effects in a restricted coil library.

A central issue in protein folding is the degree to which each residue's backbone conformational preferences stabilize the native state. We have studied the conformational preferences of each amino acid when the amino acid is not constrained to be in a regular secondary structure. In this large but highly restricted coil library, the backbone preferentially adopts dihedral angles consistent with the polyproline II conformation rather than alpha or beta conformations. The preference for the polyproline II conformation is independent of the degree of solvation. In conjunction with a new masking procedure, the frequencies in our coil library accurately recapitulate both helix and sheet frequencies for the amino acids in structured regions, as well as polyproline II propensities. Therefore, structural propensities for alpha-helices and beta-sheets and for polyproline II conformations in unfolded peptides can be rationalized solely by local effects. In addition, these propensities are often strongly affected by both the chemical nature and the conformation of neighboring residues, contrary to the Flory isolated residue hypothesis.

Computational Biology↗
PubMed2005

Mechanistic insights on the folding of a large ribozyme during transcription.

RNA folding during transcription resembles folding in a cellular environment. We previously investigated the folding of a large ribozyme derived from a bacterial RNase P RNA during its transcription by the Escherichia coli RNA polymerase and the effect of the elongation factor NusA. We found that transcriptional pausing at a specific site induced by NusA significantly altered the folding pathway. In this work, we compare folding during transcription by the E. coli RNA polymerase of circularly permuted variants and site-specific mutants of the RNase P ribozyme to elucidate the molecular mechanism of transcriptional pausing and RNA folding. The effect of NusA-induced pausing depends on the order of RNA synthesis and only affects local folding of the RNA. Pausing likely prevents a misfolded structure between the 5' strand of a helix and its adjacent junction located in the specificity domain and a region known to bind single-stranded RNA located in the catalytic domain. These results lead to a structural model on how transcriptional pausing affects folding of RNase P RNA. Structural rearrangements of a nascent RNA transcript enhanced by transcriptional pausing may be a general feature of RNA folding during transcription.

Base Sequence↗
PubMed2005

Protein folding: defining a "standard" set of experimental conditions and a preliminary kinetic data set of two-state proteins.

Recent years have seen the publication of both empirical and theoretical relationships predicting the rates with which proteins fold. Our ability to test and refine these relationships has been limited, however, by a variety of difficulties associated with the comparison of folding and unfolding rates, thermodynamics, and structure across diverse sets of proteins. These difficulties include the wide, potentially confounding range of experimental conditions and methods employed to date and the difficulty of obtaining correct and complete sequence and structural details for the characterized constructs. The lack of a single approach to data analysis and error estimation, or even of a common set of units and reporting standards, further hinders comparative studies of folding. In an effort to overcome these problems, we define here a "consensus" set of experimental conditions (25 degrees C at pH 7.0, 50 mM buffer), data analysis methods, and data reporting standards that we hope will provide a benchmark for experimental studies. We take the first step in this initiative by describing the folding kinetics of 30 apparently two-state proteins or protein domains under the consensus conditions. The goal of our efforts is to set uniform standards for the experimental community and to initiate an accumulating, self-consistent data set that will aid ongoing efforts to understand the folding process.

Biochemistry↗
PubMed2004

Efficient fluorescence labeling of a large RNA through oligonucleotide hybridization.

We present an efficient method of introducing fluorophore labels at selected locations in a large RNA. The method is based on specific and highly efficient hybridization between a fluorophore-containing DNA oligonucleotide and a modular hairpin loop replacing a functionally unimportant hairpin loop in the RNA. We demonstrate its feasibility using a 255-nucleotide RNA derived from the catalytic domain of RNase P from Bacillus subtilis. Hybridization of the DNA oligonucleotide to the modular hairpin loop minimally perturbs the structure and function of this RNA. This labeling scheme should be applicable in studies of RNA conformational dynamics by ensemble and single molecule fluorescence methods.

Bacillus subtilis↗
PubMed2004

Differences in the folding transition state of ubiquitin indicated by phi and psi analyses.

We compare the folding transition state (TS) of ubiquitin previously identified by using psi analysis to that determined by using analysis. Both methods attempt to identify interactions and their relative populations at the rate-limiting step for folding. The TS ensemble derived from psi analysis has an extensive native-like chain topology, with a four-stranded beta-sheet network and a portion of the major helix. According to analysis, however, the TS is much smaller and more polarized, with only a local helix/hairpin motif. We find that structured regions can have values far from unity, the canonical value for such sites, because of structural relaxation of the TS. Consequently, these sites may be incorrectly interpreted as contributing little to the structure of the TS. These results stress the need for caution when interpreting and drawing conclusions from analysis alone and highlight the need for more specific tools for examining the structure and energetics of the TS ensemble.

Computer Simulation↗
PubMed2004

Barrier-limited, microsecond folding of a stable protein measured with hydrogen exchange: Implications for downhill folding.

Folding experiments are conducted to test whether a covalently cross-linked coiled-coil folds so quickly that the process is no longer limited by a free-energy barrier. This protein is very stable and topologically simple, needing merely to "zipper up," while having an extrapolated folding rate of k(f) = 2 x 10(5) s(-1). These properties make it likely to attain the elusive "downhill folding" limit, at which a series of intermediates can be characterized. To measure the ultra-fast kinetics in the absence of denaturant, we apply NMR and hydrogen-exchange methods. The stability and its denaturant dependence for the hydrogen bonds in the central part of protein equal the values calculated for whole-molecule unfolding. Like-wise, their closing and opening rates indicate that these hydrogen bonds are broken and reformed in a single cooperative event representing the folding transition from the fully unfolded state to the native state. Additionally, closing rates for these hydrogen bonds agree with the extrapolated barrier-limited folding rate observed near the melting transition. Therefore, even in the absence of denaturant, where DeltaG(eq) approximately -6 kcal.mol(-1) (1 cal = 4.18 J) and tau(f) approximately 6 mus, folding remains cooperative and barrier-limited. Given that this prime candidate for downhill folding fails to do so, we propose that protein folding will remain barrier-limited for proteins that fold cooperatively.

Amino Acid Sequence↗
PubMed2004

Fast folding of a helical protein initiated by the collision of unstructured chains.

To examine whether helix formation necessarily precedes chain collision, we have measured the folding of a fully helical coiled coil that has been specially engineered to have negligible intrinsic helical propensity but high overall stability. The folding rate approaches the diffusion-limited value and is much faster than possible if folding is contingent on precollision helix formation. Therefore, the collision of two unstructured chains is the initial step of the dominant kinetic pathway, whereas helicity exerts its influence only at a later step. Folding from an unstructured encounter complex may be efficient and robust, which has implications for any biological process that couples folding to binding.

Amino Acid Sequence↗