Search PubMed⌕ Search

Biomedical subjects

Martin Gruebele

Publications and source records attributed to Martin Gruebele.

At least 19 recordsLinked to original sources

Structure-function-folding relationship in a WW domain.

Protein folding barriers result from a combination of factors including unavoidable energetic frustration from nonnative interactions, natural variation and selection of the amino acid sequence for function, and/or selection pressure against aggregation. The rate-limiting step for human Pin1 WW domain folding is the formation of the loop 1 substructure. The native conformation of this six-residue loop positions side chains that are important for mediating protein-protein interactions through the binding of Pro-rich sequences. Replacement of the wild-type loop 1 primary structure by shorter sequences with a high propensity to fold into a type-I' beta-turn conformation or the statistically preferred type-I G1 bulge conformation accelerates WW domain folding by almost an order of magnitude and increases thermodynamic stability. However, loop engineering to optimize folding energetics has a significant downside: it effectively eliminates WW domain function according to ligand-binding studies. The energetic contribution of loop 1 to ligand binding appears to have evolved at the expense of fast folding and additional protein stability. Thus, the two-state barrier exhibited by the wild-type human Pin1 WW domain principally results from functional requirements, rather than from physical constraints inherent to even the most efficient loop formation process.

Amino Acid Sequence↗

Exploring the energy landscape of a small RNA hairpin.

The energy landscape of a small RNA tetraloop hairpin is explored by temperature jump kinetics and base-substitution. The folding kinetics are single-exponential near the folding transition midpoint T(m). An additional fast phase appears below the midpoint, and an additional slow phase appears above the midpoint. Stem mutation affects the high-temperature phase, while loop mutation affects the low-temperature phase. An adjusted 2-D lattice model reproduces the temperature-dependent phases, although it oversimplifies the structural interpretation. A four-state free energy landscape model is generated based on the lattice model. This model explains the thermodynamics and multiphase kinetics over the full temperature range of the experiments. An analysis of three variants shows that one of the intermediate RNA structures is a stacking-related trap affected by stem but not loop modification, while the other is an early intermediate that forms some stem and loop structure. Even a very fast-folding 8-mer RNA with an ideal tetraloop sequence has a rugged energy landscape, ideal for testing analytical and computational models.

Kinetics↗

Low barrier kinetics: dependence on observables and free energy surface.

Dynamics on free energy surfaces with high activation barriers are usually treated by few-state kinetics models, yielding characteristic rate coefficients and amplitudes depending on the connectivity of the states. When the barriers are low (< 3 kT), the assumption of instantaneous equilibration of the transition state, and hence, the few-state kinetics treatment, break down. Langevin dynamics is used here to explore the characteristic trends that occur in such cases, as a function of barrier height, number of barriers, dimensionality of the free energy surface, and switching functions that describe how spectroscopic probes vary from reactant to product. The result is a systematic phenomenological description of low barrier kinetics and dynamics.

Computer Simulation↗

Beta-sheet folding mechanisms from perturbation energetics.

Amide backbone and sidechain mutagenesis data can be used in combination with kinetic and thermodynamic measurements to understand the energetic contributions of backbone hydrogen bonding and the hydrophobic effect to the acquisition of beta-sheet structure. For example, it has been revealed that loop 1 of the WW domain forms in the transition state, consistent with the emerging theme that reverse turn formation is rate limiting in beta-sheet folding. A distinct subset of WW domain residues principally influences thermodynamic stability by forming hydrogen bonds and hydrophobic interactions that stabilize the native state. Energetic data and sequence mining reveal that only a small subset of the molecular information contained in sequences or observed in high-resolution structures is required to generate folded functional beta-sheets, consistent with evolutionary robustness.

Animals↗

Binary and ternary aggregation within tethered protein constructs.

The free energy per monomer of a protein aggregate varies with the number of participating monomers n. The change of this free energy with aggregate size, DeltaDeltaG(n), is difficult to determine by sedimentation or concentration studies. We introduce a kinetic approach to quantitate the free energy of aggregates in the presence of tethers. By linking the protein U1A into dimers and trimers, a high effective concentration of the monomers is achieved, together with exact size control of the aggregates. We found that the free energy of the aggregate relative to the native monomer reached a maximum for n = 2, and decreased by DeltaDeltaG(2) = -3.1 kT between dimer and trimer.

Dimerization↗

Molecular vibrational energy flow and dilution factors in an anharmonic state space.

A fourth-order resonance Hamiltonian is derived from the experimental normal-mode Hamiltonian of SCCl2. The anharmonic vibrational state space constructed from the effective Hamiltonian provides a realistic model for vibrational energy flow from bright states accessible by pulsed laser excitation. We study the experimentally derived distribution PE(sigma) of dilution factors sigma as a function of energy. This distribution characterizes the dynamics in the long-time limit. State space models predict that PE(sigma) should be bimodal, with some states undergoing facile intramolecular vibrational energy redistribution (small sigma), while others at the same total energy remain "protected" (sigma approximately 1). The bimodal distribution is in qualitative agreement with analytical and numerical local density of states models. However, there are fewer states protected from energy flow, and the protected states begin to fragment at higher energy, shifting from sigma approximately 1 to sigma approximately 0.5. We also examine how dilution factors are distributed in the vibrational state space of SCCl2 and how the power law specifying the survival probability of harmonic initial states correlates with the dilution factor distribution of anharmonic initial states.

Chemistry, Physical↗

Laser absorption scanning tunneling microscopy of carbon nanotubes.

We report single molecule laser absorption by carbon nanotubes stamped under ultrahigh vacuum onto Si(100)2x1:H surfaces. Absorption is detected by scanning tunneling microscopy. Images are obtained with and without modulated laser excitation using lock-in amplification and a rear-illumination geometry to reduce thermal effects. Noise appears at topographic edges and is analyzed by a quantitative model in terms of scan speed, mechanical instabilities, and feedback current fluctuations at the edge of the nanotubes. Noise due to mechanical instabilities is shown to persist even in the limit of slow scan speed.

Journal Article↗

Solvent-tuning the collapse and helix formation time scales of lambda(6-85)*.

The lambda(6-85)(*) pseudo-wild type of lambda repressor fragment is a fast two-state folder (k(f) approximately 35 microsec(-1) at 58 degrees C). Previously, highly stable lambda(6-85)(*) mutants with k(f) > 30 microsec(-1) have been engineered to fold nearly or fully downhill. Stabilization of the native state by solvent tuning might also tune lambda(6-85)(*) away from two-state folding. We test this prediction by examining the folding thermodynamics and kinetics of lambda(6-85)(*) in a stabilizing solvent, 45% by weight aqueous ethylene glycol at -28 degrees C. Detection of kinetics by circular dichroism at 222 nm (sensitive to helix content) and small angle X-ray scattering (measuring the radius of gyration) shows that refolding from guanidine hydrochloride denatured conditions exhibits very different time scales for collapse and secondary structure formation: the two processes become decoupled. Collapse remains a low-barrier activated process, while the fastest of several secondary structure formation time scales approaches the downhill folding limit. Two-state folding of lambda(6-85)(*) is not a robust process.

Circular Dichroism↗

Engineering a beta-sheet protein toward the folding speed limit.

Recent experimental studies have shown that alpha-helical proteins can approach the folding "speed limit", where folding switches from an activated to a downhill process in free energy. beta-sheet proteins are generally thought to fold more slowly than helix bundles. However, based on studies of hairpins, folding should still be able to approach the microsecond time scale. Here we demonstrate how the hPin1 WW domain, a triple-stranded beta-sheet protein with a sharp thermodynamic melting transition, can be engineered toward the folding "speed limit" without a significant loss in thermal denaturation cooperativity.

Amino Acid Sequence↗

Downhill protein folding: evolution meets physics.

Proteins can be redesigned to fold downhill on a free energy surface characterized by only a few coordinates, confirming a principal prediction of the 'energy-landscape' model. Nonetheless, natural proteins have small but significant barriers. Spectroscopy and kinetics reveal potential biological causes for activation barriers during protein folding: evolution against protein aggregation and for protein function.

Animals↗

Kinetic equivalence of the heat and cold structural transitions of lambda6-85.

Heat- and cold-denatured proteins are considered separate thermodynamic states because temperature tuning requires the protein to pass through two 'soft' first-order phase transitions. When both pressure and temperature changes are allowed, the heat- and cold-denatured states of proteins can be interconverted without going through the native state. This raises the question of whether these states are distinguished from one another by their folding kinetics. For the Tyr22Trp/Ala37Gly/Ala49Gly mutant of the 80 residue five-helix bundle protein lambda(6-85), we show that viscosity-corrected folding rates do not distinguish the cold- and heat-denatured states. We attribute this to a folding mechanism dominated by hydrophobic collapse. Our 'temperature-symmetric' approach offers an alternative to viscosity tuning with solvent additives in such cases.

Cold Temperature↗

Kinetics are probe-dependent during downhill folding of an engineered lambda6-85 protein.

The Y22W/Q33Y/G46,48A mutant of the protein lambda6-85 folds in a few microseconds at room temperature. We find that its folding kinetics are probe-dependent under a strong bias toward the native state, a new signature for downhill folding. The IR- and fluorescence-detected relaxation time scales converge when the native bias is removed by raising the temperature, recovering activated two-state folding. Langevin dynamics simulations on one- and 2D free energy surfaces tunable from two-state to downhill folding reproduce the difference between the IR and fluorescence experiments, as well as the temperature and viscosity trends. In addition, the 2D surface reproduces the stretched exponential dynamics that we fit to the glucose solution experimental data at short times. Nonexponential dynamics at <10 micros is a signature either for local free energy minima along the reaction coordinate ("longitudinal roughness"), or for folding on a higher-dimensional free energy surface ("transverse roughness").

Kinetics↗

Folding thermodynamics and kinetics of YNMG RNA hairpins: specific incorporation of 8-bromoguanosine leads to stabilization by enhancement of the folding rate.

Modified nucleotides allow fundamental energetic and kinetic properties of nucleic acids to be probed. Here, we demonstrate that an RNA hairpin containing the nucleotide analogue 8-bromoguanosine (8BrG or G), gcUUCGgc, has enhanced stability relative to the unmodified hairpin, with DeltaDeltaG(37)(degrees)= -0.69 +/- 0.15 kcal mol(-1) and DeltaT(M) = +6.8 +/- 1.4 degrees C. NMR spectroscopic data suggest that the enhanced stability of gcUUCGgc does not arise from the native state; laser temperature-jump experiments support this notion, as gcUUCGgc and gcUUCGgc have similar unfolding rate constants, but the folding rate constant of gcUUCGgc is 4.1-fold faster at 37.5 degrees C and 2.8-fold faster under isoenergetic conditions. On the basis of these findings, we propose that 8BrG reduces the conformational entropy of the denatured state, resulting in an accelerated conformational search for the native state and enhanced stability.

Base Pairing↗

Rate-temperature relationships in lambda-repressor fragment lambda 6-85 folding.

Two classes of lambda(6-85) mutants (those richer in alanine, and those richer in glycine) have very similar slopes in an Arrhenius plot of the unfolding rates but very different temperature dependencies of the folding rates. Temperature-dependent interactions (e.g., hydrophobicity) play a large role in the initial stages of folding but not in the initial stages of unfolding of lambda(6-85). Placement of the transition state in terms of its surface exposure and entropy shows that at least two reaction coordinates are required to describe folding of all mutants over the full temperature range. The unusual Arrhenius plots of the very fastest mutant provide an additional kinetic signature for downhill folding.

Alanine↗

Context-dependent contributions of backbone hydrogen bonding to beta-sheet folding energetics.

Backbone hydrogen bonds (H-bonds) are prominent features of protein structures; however, their role in protein folding remains controversial because they cannot be selectively perturbed by traditional methods of protein mutagenesis. Here we have assessed the contribution of backbone H-bonds to the folding kinetics and thermodynamics of the PIN WW domain, a small beta-sheet protein, by individually replacing its backbone amides with esters. Amide-to-ester mutations site-specifically perturb backbone H-bonds in two ways: a H-bond donor is eliminated by replacing an amide NH with an ester oxygen, and a H-bond acceptor is weakened by replacing an amide carbonyl with an ester carbonyl. We perturbed the 11 backbone H-bonds of the PIN WW domain by synthesizing 19 amide-to-ester mutants. Thermodynamic studies on these variants show that the protein is most destabilized when H-bonds that are enveloped by a hydrophobic cluster are perturbed. Kinetic studies indicate that native-like secondary structure forms in one of the protein's loops in the folding transition state, but the backbone is less ordered elsewhere in the sequence. Collectively, our results provide an unusually detailed picture of the folding of a beta-sheet protein.

Amides↗

Detection-dependent kinetics as a probe of folding landscape microstructure.

The folding landscapes of polypeptides and proteins exhibit a hierarchy of local minima. The causes range from proline isomerization all the way down to microstructure in the free energy caused by residual frustration inherent in even the best 20 amino acid design. The corresponding time scales range from hours to submicroseconds. The smallest microstructures are difficult to detect. We have measured the folding/unfolding kinetics of the engineered trpzip2 peptide at different tryptophan fluorescence wavelengths, each yielding a different rate. Wavelength-dependent folding kinetics on 0.1-2 mus time scales show that different microstructures with a range of solvent exposure and local dynamics are populated. We estimate a lower limit for the roughness of the free energy surface based on the range of rates observed.

Animals↗

Heterogeneous folding of the trpzip hairpin: full atom simulation and experiment.

The beta-hairpin trpzip2 can be tuned continuously from a two-state folder to folding on a rough energy landscape without a dominant refolding barrier. At high denaturant concentration, this extremely stable peptide exhibits a single apparent "two-state" transition temperature when monitored by different spectroscopic techniques. However, under optimal folding conditions the hairpin undergoes an unusual folding process with three clusters of melting transitions ranging from 15 degrees C to 160 degrees C, as monitored by 12 different experimental and computational observables. We explain this behavior in terms of a rough free energy landscape of the unfolded peptide caused by multiple tryptophan interactions and alternative backbone conformations. The landscape is mapped out by potentials of mean force derived from replica-exchange molecular dynamics simulations. Implications for deducing cooperativity from denaturant titrations, for the origin of folding cooperativity, and for the folding of thermophilic proteins are pointed out. trpzip is an excellent small tunable model system for the glass-like folding transitions predicted by landscape theory.

Circular Dichroism↗

Folding lambda-repressor at its speed limit.

We show that the five-helix bundle lambda(6-85) can be engineered and solvent-tuned to make the transition from activated two-state folding to downhill folding. The transition manifests itself as the appearance of additional dynamics faster than the activated kinetics, followed by the disappearance of the activated kinetics when the bias toward the native state is increased. Our fastest value of 1 micros for the "speed" limit of lambda(6-85) is measured at low concentrations of a denaturant that smooths the free-energy surface. Complete disappearance of the activated phase is obtained in stabilizing glucose buffer. Langevin dynamics on a rough free-energy surface with variable bias toward the native state provides a robust and quantitative description of the transition from activated to downhill folding. Based on our simulation, we estimate the residual energetic frustration of lambda(6-85) to be delta(2) G approximately 0.64 k(2)T(2). We show that lambda(6-86), as well as very fast folding proteins or folding intermediates estimated to lie near the speed limit, provide a better rate-topology correlation than proteins with larger energetic frustration. A limit of beta > or = 0.7 on any stretching of lambda(6-85) barrier-free dynamics suggests that a low-dimensional free-energy surface is sufficient to describe folding.

Circular Dichroism↗