Design of an adenosine analogue that selectively improves the affinity of a mutant U1A protein for RNA

The RNA recognition motif (RRM), one of the most common RNA binding domains, contains three highly conserved aromatic amino acids that participate in stacking interactions with RNA bases. We have investigated the contribution of these highly conserved aromatic amino acids to the affinity of the complex formed between the N-terminal RRM of the U1A protein and stem loop 2 of U1 snRNA. Previously, we found that substitution of one of these conserved aromatic amino acids, Phe56, with Ala resulted in a large destabilization of the complex. Here, we have modified A6, the base in stem loop 2 RNA that stacks with Phe56, to compensate for a portion of the destabilization caused by the Phe56Ala mutation. We have designed two modified adenosines, A-3CPh and A-4CPh, in which a phenyl group is linked to the adenosine such that it may replace the phenyl group that is eliminated by the Phe56Ala mutation in the complex. We have found that incorporation of A-3CPh into stem loop 2 RNA stabilizes the complex formed with Phe56Ala by 0.6 kcal/mol, while incorporation of A-4CPh into stem loop 2 RNA stabilizes this complex by 1.8 kcal/mol. Either base modification destabilizes the wild-type complex by 0.8-0.9 kcal/mol. Experiments with other U1A mutant proteins suggest that the stabilization of the complex between the Phe56Ala U1A protein and stem loop 2 RNA is due to a specific interaction between the Phe56Ala U1A protein and A6-4CPh stem loop 2 RNA.     

Nanosecond folding dynamics of a three-stranded beta-sheet

We studied conformational stability and folding kinetics of a three-stranded beta-sheet containing two rigid turns. Static infrared measurements indicate that this beta-sheet undergoes a broad but cooperative thermal unfolding transition with a midpoint at approximately 53 degrees C. Interestingly, time-resolved infrared experiments show that its relaxation kinetics in response to a temperature-jump (T-jump) occur on the nanosecond time scale (e.g., the relaxation time is approximately 140 ns at 35.0 degrees C), thereby suggesting that the conformational relaxation encounters only a small free energy barrier or even proceeds in a downhill manner. Further Langevin dynamics simulations suggest that the observed T-jump relaxation kinetics could be modeled by a conformational diffusion process along a single-well free energy profile, which allowed us to determine the effective diffusion constant and also the roughness of the folding energy landscape.     

Molecular dynamics and thermodynamics of protein-RNA interactions: mutation of a conserved aromatic residue modifies stacking interactions and structural adaptation in the U1A-stem loop 2 RNA complex.

Molecular dynamics (MD) simulations and free energy component analysis have been performed to evaluate the molecular origins of the 5.5 kcal/mol destabilization of the complex formed between the N-terminal RNP domain of U1A and stem loop 2 of U1 snRNA upon mutation of a conserved aromatic residue, Phe56, to Ala. MD simulations, including counterions and water, have been carried out on the wild type and Phe56Ala peptide-stem loop 2 RNA complexes, the free wild type and Phe56Ala peptides, and the free stem loop 2 RNA. The MD structure of the Phe56Ala-stem loop 2 complex is similar to that of the wild type complex except the stacking interaction between Phe56 and A6 of stem loop 2 is absent and loop 3 of the peptide is more dynamic. However, the MD simulations predict large changes in the structure and dynamics of helix C and increased dynamic range of loop 3 for the free Phe56Ala peptide compared to the wild type peptide. Since helix C and loop 3 are highly variable regions of RNP domains, this indicates that a significant contribution to the reduced affinity of the Phe56Ala peptide for RNA results from cooperation between highly conserved and highly variable regions of the RNP domain of U1A. Surprisingly, these structural effects, which are manifested as cooperative free energy changes, occur in the free peptide, rather than in the complex, and are revealed only by study of both the initial and final states of the complexation process. Free energy component analysis correctly accounts for the destabilization of the Phe56Ala-stem loop 2 complex, and indicates that approximately 80% of the destabilization is due to the loss of the stacking interaction and approximately 20% is due to differences in U1A adaptation.     

The large scale conformational change of the human DPP III-substrate prefers the "closed" form

Human dipeptidyl peptidase III (DPP III) is a two domain metallo-peptidase from the M49 family. The wide interdomain cleft and broad substrate specificity suggest that this enzyme could experience significant conformational change. Long (>100 ns) molecular dynamics (MD) simulations of DPP III revealed large range conformational changes of the protein, suggesting the pre-existing equilibrium model for a substrate binding. The binding free energy calculations revealed tighter binding of the preferred synthetic substrate Arg-Arg-2-naphtylamide to the "closed" than to the "open" DPP III conformation. Our assumption that Asp372 plays a crucial role in the large scale interdomain closure was proved by the MD simulations of the Asp372Ala variant. During the same simulation time, the variant remained more "open" than the wild type protein. Apparently, Ala was not as efficient as Asp in establishing the interdomain interactions. According to the MM-PBSA calculations, the electrostatic component of the free energy of solvation turned out to be higher for the "closed" protein than for its less compact form. However, the gain in entropy due to water released from the interdomain cleft nicely balanced this negative effect.     

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.     

Structural analysis of alanine tripeptide with antiparallel and parallel beta-sheet structures in relation to the analysis of mixed beta-sheet structures in Samia cynthia ricini silk protein fiber using solid-state NMR spectroscopy

The structural analysis of natural protein fibers with mixed parallel and antiparallel beta-sheet structures by solid-state NMR is reported. To obtain NMR parameters that can characterize these beta-sheet structures, (13)C solid-state NMR experiments were performed on two alanine tripeptide samples: one with 100% parallel beta-sheet structure and the other with 100% antiparallel beta-sheet structure. All (13)C resonances of the tripeptides could be assigned by a comparison of the methyl (13)C resonances of Ala(3) with different [3-(13)C]Ala labeling schemes and also by a series of RFDR (radio frequency driven recoupling) spectra observed by changing mixing times. Two (13)C resonances observed for each Ala residue could be assigned to two nonequivalent molecules per unit cell. Differences in the (13)C chemical shifts and (13)C spin-lattice relaxation times (T(1)) were observed between the two beta-sheet structures. Especially, about 3 times longer T(1) values were obtained for parallel beta-sheet structure as compared to those of antiparallel beta-sheet structure, which could be explicable by the difference in the hydrogen-bond networks of both structures. This very large difference in T(1) becomes a good measure to differentiate between parallel or antiparallel beta-sheet structures. These differences in the NMR parameters found for the tripeptides may be applied to assign the parallel and antiparallel beta-sheet (13)C resonances in the asymmetric and broad methyl spectra of [3-(13)C]Ala silk protein fiber of a wild silkworm, Samia cynthia ricini.     

Examination of laser microbeam cell lysis in a PDMS microfluidic channel using time-resolved imaging

We use time-resolved imaging to examine the lysis dynamics of non-adherent BAF-3 cells within a microfluidic channel produced by the delivery of single highly-focused 540 ps duration laser pulses at lambda = 532 nm. Time-resolved bright-field images reveal that the delivery of the pulsed laser microbeam results in the formation of a laser-induced plasma followed by shock wave emission and cavitation bubble formation. The confinement offered by the microfluidic channel constrains substantially the cavitation bubble expansion and results in significant deformation of the PDMS channel walls. To examine the cell lysis and dispersal of the cellular contents, we acquire time-resolved fluorescence images of the process in which the cells were loaded with a fluorescent dye. These fluorescence images reveal cell lysis to occur on the nanosecond to microsecond time scale by the plasma formation and cavitation bubble dynamics. Moreover, the time-resolved fluorescence images show that while the cellular contents are dispersed by the expansion of the laser-induced cavitation bubble, the flow associated with the bubble collapse subsequently re-localizes the cellular contents to a small region. This capacity of pulsed laser microbeam irradiation to achieve rapid cell lysis in microfluidic channels with minimal dilution of the cellular contents has important implications for their use in lab-on-a-chip applications.     

Probing polar solvation dynamics in proteins: a molecular dynamics simulation analysis

Measurements of time-resolved Stokes shifts on picosecond to nanosecond time scales have been used to probe the polar solvation dynamics of biological systems. Since it is difficult to decompose the measurements into protein and solvent contributions, computer simulations are useful to aid in understanding the details of the molecular behavior. Here we report the analysis of simulations of the electrostatic interactions of the rest of the protein and the solvent with 11 residues of the immunoglobulin binding domain B1 of protein G. It is shown that the polar solvation dynamics are position-dependent and highly heterogeneous. The contributions due to interactions with the protein and with the solvent are determined. The solvent contributions are found to vary from negligible after a few picoseconds to dominant on a scale of hundreds of picoseconds. The origin for the latter is found to involve coupled hydration and protein conformational dynamics. The resulting microscopic picture demonstrates that a wide range of possibilities have to be considered in the interpretation of time-resolved Stokes shift measurements.     

Time-resolved fluorescence study of the single tryptophan in thiocyanate and azide derivatives of horseradish peroxidase: Implication for a p H-induced conformational change in the heme cavity

Detailed pH-dependent steady state and picosecond time-resolved tryptophan fluorescence studies on thiocyanate and azide complexes of horseradish peroxidase have been carried out. The fluorescence decay of the single tryptophan in these species was fitted to a discrete three exponential model. Maximum entropy method analysis also gave three distinct regions of lifetime distributions. The fast subnanosecond lifetime component was found to have > 97% amplitude contribution while other two longer lifetime components have small contributions. Small contributions from the nanosecond lifetime components possibly arise from apoprotein impurity or some small amount of disordered heme conformer of the protein. pH dependence of the fast picosecond lifetime components was found to show a systematic behavior which has been interpreted in the light of obligatory conformation change associated with activation of the enzyme at low pH.     

The energetic and structural effects of single amino acid substitutions upon capped alpha-helical peptides containing 17 amino acid residues. An ONIOM DFT/AM1 study

We evaluate the effect of the amino acid mutations of glycine, leucine, valine, phenylalanine, serine, and proline for the 10th alanine in the capped peptide, acetly(Ala)(17)NH(2), upon the energies of the alpha-helices and beta-strands using ONIOM DFT/AM1 molecular orbital calculations. The relative stabilities of the alpha-helix (to the beta-strand) derive from the differences between the effects upon not only the helix but the strand as well. Thus, Ala --> Pro significantly destabilizes both but destabilizes the alpha-helix more, while Ala --> Gly stabilizes both but stabilizes the beta-strand more. The theoretical results are discussed in the context of the known experimental reports. We suggest that the solvation of the unfolded state drives the helix/coil equilibrium in solution.     

Quenching and blinking of fluorescence of a single dye molecule bound to gold nanoparticles

A fluorescein derivative (SAMSA) bound to gold nanoparticles of different diameters is investigated by time-resolved fluorescence at the single molecule level in a wide dynamic range, from nanosecond to second time scale. The significant decrease of both SAMSA excited state lifetime and fluorescence quantum yield observed upon binding to gold nanoparticles can be essentially traced back to an increase of the nonradiative deactivation rate, probably due to energy transfer, that depends on the nanoparticle size. A slow single molecule fluorescence blinking, in the ms time scale, has a marked dependence on the excitation intensity both under single and under two photon excitation. The blinking dynamics is limited by a low probability nonlinear excitation to a high energy state from which a transition to a dark state occurs. The results point out a strong coupling between the vibro-electronic configuration of the dye and the plasmonic features of the metal nanoparticles that provide dye radiationless deactivation channels on a wide dynamic range.     

Controlling the DNA binding specificity of bHLH proteins through intramolecular interactions

Reversible control of the conformation of proteins was employed to probe the relationship between flexibility and specificity of the basic helix-loop-helix protein MyoD. A fusion protein (apaMyoD) was designed where the basic DNA binding helix of MyoD was stablized by an amino-terminal extension with a sequence derived from the bee venom peptide apamin. The disulfide-stabilized helix from apamin served as a nucleus for a helix that extended for a further ten residues, thereby holding apaMyoD's DNA recognition helix in a predominantly alpha-helical conformation. The thermal stability of the DNA complexes of apaMyoD was increased by 13 degrees C relative to MyoD-bHLH. Measurements of the fluorescence anisotropy change on DNA binding indicated that apaMyoD bound to E-box-containing DNA sequences with enhanced affinity relative to MyoD-bHLH. Consequently, the DNA binding specificity of apaMyoD was increased 10-fold.     

Measurement of solvation responses at multiple sites in a globular protein

Proteins respond to electrostatic perturbations through complex reorganizations of their charged and polar groups, as well as those of the surrounding media. These solvation responses occur both in the protein interior and on its surface, though the exact mechanisms of solvation are not well understood, in part because of limited data on the solvation responses for any given protein. Here, we characterize the solvation kinetics at sites throughout the sequence of a small globular protein, the B1 domain of streptococcal protein G (GB1), using the synthetic fluorescent amino acid Aladan. Aladan was incorporated into seven different GB1 sites, and the time-dependent Stokes shift was measured over the femtosecond to nanosecond time scales by fluorescence upconversion and time-correlated single photon counting. The seven sites range from buried within the protein core to fully solvent-exposed on the protein surface, and are located on different protein secondary structures including beta-sheets, helices, and loops. The dynamics in the protein sites were compared against the free fluorophore in buffer. All protein sites exhibited an initial, ultrafast Stokes shift on the subpicosecond time scale similar to that observed for the free fluorophore, but smaller in magnitude. As the probe is moved from the surface to more buried sites, the dynamics of the solvation response become slower, while no clear correlation between dynamics and secondary structure is observed. We suggest that restricted movements of the surrounding protein residues give rise to the observed long time dynamics and that such movements comprise a large portion of the protein's solvation response. The proper treatment of dynamic Stokes shift data when the time scale for solvation is comparable to the fluorescence lifetime is discussed.     

Intramolecular disulfide bridges as a phototrigger to monitor the dynamics of small cyclic peptides

Two cyclic disulfide-bridged tetrapeptides [cyclo(Boc-Cys-Pro-Aib-Cys-OMe) (1) and cyclo(Boc-Cys-Pro-Phe-Cys-OMe) (2)] have been monitored by time-resolved mid-IR spectroscopy in the C=O vibrational range. A conformational change is induced by cleavage of the intramolecular disulfide bridge upon UV excitation (lambda(exc) = 260 nm), giving rise to a pair of cysteinyl radicals (thiyl radicals), which diffuse apart allowing the peptide to change conformation before they undergo quenching. The amide I band reports on the dynamics of the peptide backbone, which evolves on a 100 ps time scale and then stays constant up to 10 micros at low enough concentrations ( approximately 100 mM). To probe specifically the lifetime of the free cysteinyl radicals, time-resolved UV laser flash photolysis has been applied. The concentration of the cysteinyl radical decays nonexponentially, but about 50% are still present after 1 ms. The photocleavable disulfide bridge hence may serve as an intrinsic, naturally occurring phototrigger to study peptide dynamics that opens a wide time-window from a few picoseconds to many hundreds of microseconds.     

Identifying residual structure in intrinsically disordered systems: a 2D IR spectroscopic study of the GVGXPGVG peptide

The peptide amide-I vibration of a proline turn encodes information on the turn structure. In this study, FTIR, two-dimensional IR spectroscopy and molecular dynamics simulations were employed to characterize the varying turn conformations that exist in the GVGX(L)PGVG family of disordered peptides. This analysis revealed that changing the size of the side chain at the X amino acid site from Gly to Ala to Val substantially alters the conformation of the peptide. To quantify this effect, proline peak shifts and intensity changes were compared to a structure-based spectroscopic model. These simulated spectra were used to assign the population of type-II β turns, bulged turns, and irregular β turns for each peptide. Of particular interest was the Val variant commonly found in the protein elastin, which contained a 25% population of irregular β turns containing two peptide hydrogen bonds to the proline C═O.     

Backbone motions of free and pheromone-bound major urinary protein I studied by molecular dynamics simulation

Molecular motions of free and pheromone-bound mouse major urinary protein I, previously investigated by NMR relaxation, were simulated in 30 ns molecular dynamics (MD) runs. The backbone flexibility was described in terms of order parameters and correlation times, commonly used in the NMR relaxation analysis. Special attention was paid to the effect of conformational changes on the nanosecond time scale. Time-dependent order parameters were determined in order to separate motions occurring on different time scales. As an alternative approach, slow conformational changes were identified from the backbone torsion angle variances, and "conformationally filtered" order parameters were calculated for well-defined conformation states. A comparison of the data obtained for the free and pheromone-bound protein showed that some residues are more rigid in the bound form, but a larger portion of the protein becomes more flexible upon the pheromone binding. This finding is in general agreement with the NMR results. The higher flexibility observed on the fast (fs-ps) time scale was typically observed for the residues exhibiting higher conformational freedom on the ns time scale. An inspection of the hydrogen bond network provided a structural explanation for the flexibility differences between the free and pheromone-bound proteins in the simulations.     

A suite of 2H NMR spin relaxation experiments for the measurement of RNA dynamics

A suite of (2)H-based spin relaxation NMR experiments is presented for the measurement of molecular dynamics in a site-specific manner in uniformly (13)C, randomly fractionally deuterated ( approximately 50%) RNA molecules. The experiments quantify (2)H R(1) and R(2) relaxation rates that can subsequently be analyzed to obtain information about dynamics on a pico- to nanosecond time scale. Sensitivity permitting, the consistency of the data can be evaluated by measuring all five rates that are accessible for a spin 1 particle and establishing that the rates obey relations that are predicted from theory. The utility of the methodology is demonstrated with studies of the dynamics of a 14-mer RNA containing the UUCG tetraloop at temperatures of 25 and 5 degrees C. The high quality of the data, even at 5 degrees C, suggests that the experiments will be of use for the study of RNA molecules that are as large as 30 nucleotides.     

Changes in dynamics of SRE-RNA on binding to the VTS1p-SAM domain studied by 13C NMR relaxation

RNA recognition by proteins is often accompanied by significant changes in RNA dynamics in addition to conformational changes. However, there are very few studies which characterize the changes in molecular motions in RNA on protein binding. We present a quantitative (13)C NMR relaxation study of the changes in RNA dynamics in the pico-nanosecond time scale and micro-millisecond time scale resulting from interaction of the stem-loop SRE-RNA with the VTS1p-SAM domain. (13)C relaxation rates of the protonated carbons of the nucleotide base and anomeric carbons have been analyzed by employing the model-free formalism, for a fully (13)C/(15)N-labeled sample of the SRE-RNA in the free and protein-bound forms. In the free RNA, the nature of molecular motions are found to be distinctly different in the stem and the loop region. On binding to the protein, the nature of motions becomes more homogeneous throughout the RNA, with many residues showing increased flexibility at the aromatic carbon sites, while the anomeric carbon sites become more rigid. Surprisingly, we also observe indications of a slow collective motion of the RNA in the binding pocket of the protein. The observation of increased motions on binding is interesting in the context of growing evidence that binding does not always lead to motional restrictions and the resulting entropy gain could favor the free energy of association.