Biomedical subjects

Jan Kubelka

Publications and source records attributed to Jan Kubelka.

15 recordsLinked to original sources

PubMed2006

Simulation of infrared spectra for beta-hairpin peptides stabilized by an Aib-Gly turn sequence: correlation between conformational fluctuation and vibrational coupling.

Vibrational spectra of a 12-residue beta-hairpin peptide, RYVEVBGKKILQ (HBG), stabilized by an Aib-Gly turn sequence (B = Aib) were investigated theoretically using a combination of molecular dynamics (MD) and density functional theory (DFT) calculations. Selected conformations of HBG were extracted from a classical MD trajectory and used for spectral simulations. DFT calculations, based on the Cartesian coordinate spectral property transfer protocol, were carried out for peptide structures in which all residues are replaced with Ala, except for the Aib and Gly residues, but the backbone (phi, psi, omega) structure of the original configuration is retained. The simulations provide a basis for interpretation of the HBG amide I infrared spectra in terms of structural variables such as detailed secondary structure and thermal conformational fluctuation as well as vibrational coupling as indicated by spectra of 13C isotope-labeled variants. The characteristic amide I band shape of such small beta-hairpin peptides appears to arise from the structure of the short antiparallel beta-sheet strands. The role of structural parameter fluctuation in vibrational coupling is evaluated by comparison of DFT-derived amide coupling constants for selected configurations and from transition dipole coupling calculations of coupling parameters between (13)C isotopically labeled residues for a MD-derived ensemble of configurations. Calculated results were compared with the experimentally obtained spectra for several (13)C isotope-labeled peptides of this sequence.

Journal Article↗

PubMed2006

Specificity of the initial collapse in the folding of the cold shock protein.

The two-state folding reaction of the cold shock protein from Bacillus caldolyticus (Bc-Csp) is preceded by a rapid chain collapse. A fast shortening of intra-protein distances was revealed by Förster resonance energy transfer (FRET) measurements with protein variants that carried individual pairs of donor and acceptor chromophores at various positions along the polypeptide chain. Here we investigated the specificity of this rapid compaction. Energy transfer experiments that probed the stretching of strand beta2 and the close approach between the strands beta1 and beta2 revealed that the beta1-beta2 hairpin is barely formed in the collapsed form, although it is native-like in the folding transition state of Bc-Csp. The time course of the collapse could not be resolved by pressure or temperature jump experiments, indicating that the collapsed and extended forms are not separated by an energy barrier. The co-solute (NH4)2SO4 stabilizes both native Bc-Csp and the collapsed form, which suggests that the large hydrated SO4(2-) ions are excluded from the surface of the collapsed form in a similar fashion as they are excluded from folded Bc-Csp. Ethylene glycol increases the stability of proteins because it is excluded preferentially from the backbone, which is accessible in the unfolded state. The collapsed form of Bc-Csp resembles the unfolded form in its interaction with ethylene glycol, suggesting that in the collapsed form the backbone is still accessible to water and small molecules. Our results thus rule out that the collapsed form is a folding intermediate with native-like chain topology. It is better described as a mixture of compact conformations that belong to the unfolded state ensemble. However, some of its structural elements are reminiscent of the native protein.

Ammonium Sulfate↗

PubMed2006

Sub-microsecond protein folding.

We have investigated the structure, equilibria, and folding kinetics of an engineered 35-residue subdomain of the chicken villin headpiece, an ultrafast-folding protein. Substitution of two buried lysine residues by norleucine residues stabilizes the protein by 1 kcal/mol and increases the folding rate sixfold, as measured by nanosecond laser T-jump. The folding rate at 300 K is (0.7 micros)(-1) with little or no temperature dependence, making this protein the first sub-microsecond folder, with a rate only twofold slower than the theoretically predicted speed limit. Using the 70 ns process to obtain the effective diffusion coefficient, the free energy barrier height is estimated from Kramers theory to be less than approximately 1 kcal/mol. X-ray crystallographic determination at 1A resolution shows no significant change in structure compared to the single-norleucine-substituted molecule and suggests that the increased stability is electrostatic in origin. The ultrafast folding rate, very accurate X-ray structure, and small size make this engineered villin subdomain an ideal system for simulation by atomistic molecular dynamics with explicit solvent.

Amino Acid Substitution↗

PubMed2005

High-resolution x-ray crystal structures of the villin headpiece subdomain, an ultrafast folding protein.

The 35-residue subdomain of the villin headpiece (HP35) is a small ultrafast folding protein that is being intensely studied by experiments, theory, and simulations. We have solved the x-ray structures of HP35 and its fastest folding mutant [K24 norleucine (nL)] to atomic resolution and compared their experimentally measured folding kinetics by using laser temperature jump. The structures, which are in different space groups, are almost identical to each other but differ significantly from previously solved NMR structures. Hence, the differences between the x-ray and NMR structures are probably not caused by lattice contacts or crystal/solution differences, but reflect the higher accuracy of the x-ray structures. The x-ray structures reveal important details of packing of the hydrophobic core and some additional features, such as cross-helical H bonds. Comparison of the x-ray structures indicates that the nL substitution produces only local perturbations. Consequently, the finding that the small stabilization by the mutation is completely reflected in an increased folding rate suggests that this region of the protein is as structured in the transition state as in the folded structure. It is therefore a target for engineering to increase the folding rate of the subdomain from approximately 0.5 micros(-1) for the nL mutant to the estimated theoretical speed limit of approximately 3 micros(-1).

Crystallography, X-Ray↗

PubMed2005

Solvent effects on IR and VCD spectra of helical peptides: DFT-based static spectral simulations with explicit water.

Simulations of IR and VCD spectra are carried out for model alpha-helical, 3(10)-helical, and 3(1)-helical (polyProII-like) oligopeptides, with up to 21 amide groups, and including explicit consideration of effects of directly hydrogen-bonded solvent (water). Parameters used were obtained from ab initio density functional theory (DFT) computations of force field, atomic polar and axial tensors for oligopeptides of 5 to 7 amides, whose structures were constrained in (phi,psi) to target the secondary structure type but otherwise fully optimized. By comparison with experimental data as well as with calculations for identical but isolated (gas phase) peptides, the computed effects of an inner shell of aqueous solvent on the vibrational spectra of helical oligopeptides are illustrated. The interaction with solvent causes significant frequency shifts of the amide bands, but only minor changes in the characteristic IR intensity distributions and splittings and the VCD band shapes. Better agreement with experimental band shapes is achieved for the alpha-helical amide I' (N-deuterated) VCD by inclusion of explicit solvent in the calculations. Some improvements are also observed in theoretical VCD predictions for 13C labeled alpha-helical peptides when solvated models are used. However, the qualitative isotopic splitting patterns are preserved and just shifted in frequency due to consistent, solvent independent interamide coupling constants. The critical match of experiment and theory for relative positions of transitions in peptides with specifically separated 13C=O labels, including and neglecting solvent, confirms the stability of the coupling interactions. Despite these solvation effects, the calculated VCD band shape of the amide I mode is shown to be a reliable conformational probe, since it remains basically insensitive to frequency shifts caused by environment. Thus theoretical VCD simulations, even vacuum calculations, are shown to provide useful spectral predictions for solution-phase peptides.

Circular Dichroism↗

PubMed2005

IR study of cross-strand coupling in a beta-hairpin peptide using isotopic labels.

Model beta-hairpin peptides can be used to develop understanding of fundamental elements of beta-sheet secondary structure formation and stability. We have studied two 13C-labeled variants of a beta-hairpin peptide modified from a design originally proposed by Gellman: Arg-Tyr-Val-Glu-Val-Aib-Gly-Lys-Lys-Ile-Leu-Gln. (In this peptide, the two italicized residues form a beta-turn, while 13C-labels are on the amide C=O of Val3, Lys8 in HBG-L and Val3, Ile10 in HBG-S.) Both these peptides are labeled on opposite strands of the hairpin, but differ in the labeling pattern. One (HBG-L) forms a large (14-atom) H-bonded ring of labeled C=Os, while the other (HBG-S) forms a small (10-atom) H-bonded ring. These impact the amide I infrared spectra, with HBG-L having a 13C frequency and intensity higher than that of HBG-S, in good agreement with our spectral simulations based on quantum mechanically derived force fields. The thermal behavior of both peptides yields a broad thermal transition and lacks an isosbestic point. The 13C band for HBG-L has the largest intensity change with temperature, distinct from the 12C change and the HBG-S 13C change.

Amino Acid Sequence↗

PubMed2005

Determination of ultrafast protein folding rates from loop formation dynamics.

Quenching of the triplet state of tryptophan by contact with cysteine can be used to measure the kinetics of loop formation in unfolded proteins. Here we show that cysteine quenching dynamics also provide a novel method for measuring folding rates when the exchange between folded and unfolded states is faster than the unquenched triplet lifetime (approximately 100 micros). We use this technique to investigate folding/unfolding kinetics of the 35 residue headpiece subdomain of the protein villin, which contains a single tryptophan residue and was engineered to contain a cysteine residue at the N terminus. At intermediate concentrations of denaturant the time-course of the triplet decay consists of two relaxations, the rates and amplitudes of which reveal the fast kinetics for folding and unfolding of this protein. The folding rates extracted using a simple kinetic model are close to those reported previously from laser-induced temperature-jump experiments that employ the change in tryptophan fluorescence as a probe. However, the results differ significantly from those reported from dynamic NMR line shape analysis on a variant with methionine at the N terminus, an issue that remains to be resolved. The analysis of the triplet quenching kinetics also shows that the quenching rates in the unfolded state increase with decreasing denaturant concentration, indicating a compaction of the unfolded protein.

Antioxidants↗

PubMed2004

Nature of vibrational coupling in helical peptides: an isotopic labeling study.

Infrared (IR) and vibrational circular dichroism (VCD) spectra were measured for a series of isotopically ((13)C on two or more amide Cdouble bond]O) labeled, 25 residue, alpha-helical peptides of the sequence Ac-(AAAAK)(4)AAAAY-NH(2) that were also studied in the previous paper. Theoretical IR and VCD simulations were performed for correspondingly isotopically labeled Ac-A(24)-NHCH(3) constrained to an alpha-helical conformation by use of property tensor transfer from density functional theory (DFT) calculations on Ac-A(10)-NHCH(3). The simulations predicted and experiments confirmed that the vibrational coupling constants between i, i + 1 and i, i + 2 residues differ in sign, thus leading to a reversal of the (13)C VCD pattern and explaining the large shift in the (13)C amide I frequency as reported in the previous paper. The sign of the coupling constant remained consistent for larger label separation (with the exception of i, i + 4) and for more labels with uniform separation. Such effects confirm that the isotopically labeled group vibrations are essentially only coupled to each other and are effectively uncoupled from those of the unlabeled groups. This development confirms the utility of isotopic labels for site-specific structural studies with vibrational spectra. Observed spectral effects cannot be explained by considering only transition dipole coupling (TDC) between amide oscillators, particularly for smaller label separations, but the TDC and ab initio predicted couplings roughly converge at large separation.

Alanine↗

PubMed2004

The protein folding 'speed limit'.

How fast can a protein possibly fold? This question has stimulated experimentalists to seek fast folding proteins and to engineer them to fold even faster. Proteins folding at or near the speed limit are prime candidates for all-atom molecular dynamics simulations. They may also have no free energy barrier, allowing the direct observation of intermediate structures on the pathways from the unfolded to the folded state. Both experimental and theoretical approaches predict a speed limit of approximately N/100micros for a generic N-residue single-domain protein, with alpha proteins folding faster than beta or alphabeta. The predicted limits suggest that most known ultrafast folding proteins can be engineered to fold more than ten times faster.

Computer Simulation↗

PubMed2003

Optical spectroscopic investigations of model beta-sheet hairpins in aqueous solution.

In this contribution we report optical spectroscopic data on a series of designed beta hairpins previously shown by NMR to contain a substantial population of beta-sheet structure. These models contain a designed hydrophobic cluster and a (D)Pro-Gly sequence to promote formation of a turn geometry. FTIR, electronic and vibrational CD (ECD and VCD) spectra for these small peptides are comparable to expected bandshapes for peptides of high beta-sheet content. The (D)Pro-Gly sequence provides a better turn motif than Asn-Gly as measured by its beta-sheet spectral characteristics. IR and VCD spectra are in qualitative agreement with theoretical simulations based on transfer of parameters from ab initio quantum mechanical force field and intensity computations for the turn and strands. These calculations provide assignments for some distinguishing modes in both IR and VCD spectra. Increased sheet structure can be induced in these hairpins by use of mixed solvent conditions. Thermal denaturation studies reveal that these hairpins undergo very broad unfolding transitions. Guanidine hydrochloride unfolding transitions for the selected hairpin models are similarly broad. However, the "end-states" of temperature and chaotropic denaturation are spectroscopically differentiable.

Circular Dichroism↗

PubMed2003

Experimental tests of villin subdomain folding simulations.

We have used laser temperature-jump to investigate the kinetics and mechanism of folding the 35 residue subdomain of the villin headpiece. The relaxation kinetics are biphasic with a sub-microsecond phase corresponding to a helix-coil transition and a slower microsecond phase corresponding to overall unfolding/refolding. At 300 K, the folding time is 4.3(+/-0.6) micros, making it the fastest folding, naturally occurring protein, with a rate close to the theoretical speed limit. This time is in remarkable agreement with the prediction of 5 (+11,-3) micros by Zagrovic et al. from atomistic molecular dynamics simulations using an implicit solvent model. We test their prediction that replacement of the C-terminal phenylalanine residue with alanine will increase the folding rate by removing a transient non-native interaction. We find that the alanine substitution has no effect on the folding rate or on the equilibrium constant. Implications of this result for the validity of the simulated folding mechanism are discussed.

Alanine↗

PubMed2002

Discriminating 3(10)- from alpha-helices: vibrational and electronic CD and IR absorption study of related Aib-containing oligopeptides.

Model peptides based on -(Aib-Ala)(n)-, and (Aib)(n)-Leu-(Aib)(2) sequences, which have varying amounts of 3(10)-helical character, were studied by use of vibrational and electronic circular dichroism (VCD and ECD) and Fourier transform infrared (FTIR) absorption spectroscopies to test the correlation of spectral response and conformation. The data indicate that these peptides, starting from a length of about four to six residues, predominantly adopt a 3(10)-helical conformation at room temperature. The longest model peptides, depending on the series, may evidence some alpha-helical contribution to the spectra, while the shorter ones, with less than six residues, have much less order. The IR absorption spectra (as supported by theory) showed only small frequency changes between 3(10)- and alpha-helices. By contrast, solvent effects are a source of much bigger perturbations. The ECD results show that the intensity ratio for the approximately 222-nm to approximately 208-nm bands, while useful for distinguishing between these two helical types in some sequences, may have a narrower range of application than VCD. However, the VCD data presented here continue to support the proposed discrimination between alpha- and 3(10)-helices based on qualitative amide I and II bandshape differences. The present study shows the intensities of the 3(10)-helical amide I (peak-to-peak) to its amide II VCD to be of the same order and useful for discriminating them from alpha-helices, whose amide I dominates the amide II in intensity. This qualitative result is experimentally independent of the amount of alphaMe-substituted residues in the sequence. These experimental VCD results are consistent in detail with theoretical spectral simulations for Ac-(Ala)(8)-NH(2), Ac-(Aib-Ala)(4)-NH(2), and Ac-(Aib)(8)-NH(2) in 3(10)- and alpha-helical conformations.

Amino Acid Sequence↗

PubMed2002

Ab initio quantum mechanical models of peptide helices and their vibrational spectra.

Structural parameters for standard peptide helices (alpha, 3(10), 3(1) left-handed) were fully ab initio optimized for Ac-(L-Ala)(9)-NHMe and for Ac-(L-Pro)(9)-NHMe (poly-L-proline-PLP I and PLP II-forms), in order to better understand the relative stability and minimum energy geometries of these conformers and the dependence of the ir absorption and vibrational CD (VCD) spectra on detailed variation in these conformations. Only the 3(10)-helical Ala-based conformation was stable in vacuum for this decaamide structure, but both Pro-based conformers minimized successfully. Inclusion of solvent effects, by use of the conductor-like screening solvent model (COSMO), enabled ab initio optimizations [at the DFT/B3LYP/SV(P) level] without any constraints for the alpha- and 3(10)-helical Ala-based peptides as well as the two Pro-based peptides. The geometries obtained compare well with peptide chain torsion angles and hydrogen-bond distances found for these secondary structure types in x-ray structures of peptides and proteins. For the simulation of VCD spectra, force field and intensity response tensors were obtained ab initio for the complete Ala-based peptides in vacuum, but constrained to the COSMO optimized torsional angles, due to limitations of the solvent model. Resultant spectral patterns reproduce well many aspects of the experimental spectra and capture the differences observed for these various helical types.

Crystallography, X-Ray↗

PubMed2002

Discrimination between peptide 3(10)- and alpha-helices. Theoretical analysis of the impact of alpha-methyl substitution on experimental spectra.

Detailed spectral simulations based on ab initio density functional theory computations of the amide I and II infrared (IR) and vibrational circular dichroism (VCD) spectra for Ac-(Ala)(4)-NH(2), Ac-(Aib-Ala)(2)-NH(2), and Ac-(Aib)(4)-NH(2) constrained to 3(10)- and alpha-helical conformations are presented. Parameters from these ab initio calculations are transferred onto corresponding larger oligopeptides to simulate the spectra for dodecamers. The differences between conformations and for different Aib substitution patterns within a conformation are reflected in observable spectral patterns where data are available. Simulated IR spectra show small frequency shifts in the amide I maxima between 3(10)- and alpha-helices, but the same magnitude shifts occur within one conformation upon Aib substitution. Thus, from a computational basis, the frequency of the amide I maximum does not discriminate between the 3(10)- or alpha-helical conformations. Calculated VCD band shapes for 3(10)-helices showed more significant changes in amplitude, with change in the fraction of Aib, than those for alpha-helices. Generally, with increasing Aib content, the overall amide I VCD intensity becomes weaker and the amide I couplet becomes more conservative, while the amide II VCD is less affected. Although the detailed band shape is shown to be sensitive to alpha-Me substitution, the basic pattern of amide I and II relative VCD intensities still differs between alpha- and 3(10)-helices and, as a consequence, successfully discriminates between them. These predictions are all borne out in experimental spectra of Aib, mixed Aib-Ala, and Ala-based helical peptides, where available.

Alanine↗

PubMed2002

Spectroscopic characterization of selected beta- sheet hairpin models.

The IR and vibrational circular dichroism (VCD) spectra of a model two-stranded beta hairpin are compared to those of a related cyclic two-stranded model, which are both stabilized by DPro- Gly turns. The spectra are compared to ab initio based simulations to support specific assignments of the dominant features and suggest a revised interpretation of the IR and VCD spectra for beta sheet containing proteins.

Circular Dichroism↗