Monitoring aromatic picosecond to nanosecond dynamics in proteins via 13C relaxation: expanding perturbation mapping of the rigidifying core mutation, V54A, in eglin c

Long-range effects, such as allostery, have evolved in proteins as a means of regulating function via communication between distal sites. An NMR-based perturbation mapping approach was used to more completely probe the dynamic response of the core mutation V54A in the protein eglin c by monitoring changes in picosecond to nanosecond aromatic side-chain dynamics and H/D exchange stabilities. Previous side-chain dynamics studies on this mutant were limited to methyl-bearing residues, most of which were found to rigidify on the picosecond to nanosecond time scale in the form of a contiguous "network". Here, high precision (13)C relaxation data from 13 aromatic side chains were acquired by applying canonical relaxation experiments to a newly developed carbon labeling scheme [Teilum et al. (2006) J. Am. Chem. Soc. 128, 2506-2507]. The fitting of model-free parameters yielded S (2) variability which is intermediate with respect to backbone and methyl-bearing side-chain variability and tau e values that are approximately 1 ns. Inclusion of the aromatic dynamic response results in an expanded network of dynamically coupled residues, with some aromatics showing increases in flexibility, which partially offsets the rigidification in methyl side chains. Using amide hydrogen exchange, dynamic propagation on a slower time scale was probed in response to the V54A perturbation. Surprisingly, regional stabilization (slowed exchange) 10-12 A from the site of mutation was observed despite a global destabilization of 1.5 kcal x mol (-1). Furthermore, this unlikely pocket of stabilized residues colocalizes with increases in aromatic flexibility on the faster time scale. Because the converse is also true (destabilized residues colocalize with rigidification on the fast time scale), a plausible entropy-driven mechanism is discussed for relating colocalization of opposing dynamic trends on vastly different time scales.     

Monoclonal antibodies demonstrate limited structural homology between myosin isozymes from Acanthamoeba

We used a library of 31 monoclonal and six polyclonal antibodies to compare the structures of the two classes of cytoplasmic myosin isozymes isolated from Acanthamoeba: myosin-I, a 150,000-mol-wt, globular molecule; and myosin-II, a 400,000-mol-wt molecule with two heads and a 90-nm tail. This analysis confirms that myosin-I and -II are unique gene products and provides the first evidence that these isozymes have at least one structurally homologous region functionally important for myosin's role in contractility. Characterization of the 23 myosin-II monoclonal antibody binding sites by antibody staining of one-dimensional peptide maps and solid phase, competitive binding assays demonstrate that they bind to at least 15 unique sites on the myosin-II heavy chain. The antibodies can be grouped into six families, whose members bind close to one another. None of the monoclonal antibodies bind to myosin-II light chains and polyclonal antibodies against myosin-II light or heavy chain bind only to myosin-II light or heavy chains, respectively: no antibody binds both heavy and light chains. Six of eight monoclonal antibodies and one of two polyclonal sera that react with the myosin-I heavy chain also bind to determinants on the myosin-II heavy chain. The cross-reactive monoclonal antibodies bind to the region of myosin-II recognized by the largest family of myosin-II monoclonal antibodies. In the two papers that immediately follow, we show that this family of monoclonal antibodies to myosin-II binds to the myosin-II tail near the junction with the heads and inhibits both the actin-activated ATPase of myosin-II and contraction of gelled cytoplasmic extracts of Acanthamoeba cytoplasm. Further, this structurally homologous region may play a key role in energy transduction by cytoplasmic myosins.     

Large-amplitude picosecond anisotropy decay of the intrinsic fluorescence of double-stranded DNA.

The conformational flexibility of the DNA double helix is of great interest because of its potential role in protein recognition, packaging into chromosomes, formation of photodefects, and interaction with drugs. Theory finds that DNA is very flexible; however, there is a scarcity of experimental results that examine intrinsic properties of the DNA bases for the inherent flexibility in solution. We have studied the dynamics of poly(dA).poly(dT) and (dA)20.(dT)20 in a 50 mM cacodylate, 0.1 M NaCl, pH 7 buffer by using the time-correlated picosecond fluorescence anisotropy of thymine selectively excited at 293 nm. For both nucleic acids, a large-amplitude biphasic decrease in the anisotropy is observed that has a very fast, large-amplitude component on the picosecond time scale and a slower, smaller-amplitude component on the nanosecond time scale. These modes are sensitive to sucrose concentration, and are greatly attenuated at 77% sucrose by volume. This observation suggests that motions of the bases make a significant contribution to the observed fluorescence depolarization (in the absence of sucrose). Measurements on the single-stranded systems poly(dT) and (dT)20 reveal a much smaller amplitude of the very fast depolarization mode. These observations are consistent with a mechanism that involves concerted motions in the interior of the double-stranded systems.     

Interpretation of NMR relaxation properties of Pin1, a two-domain protein, based on Brownian dynamic simulations

Many important proteins contain multiple domains connected by flexible linkers. Inter-domain motion is suggested to play a key role in many processes involving molecular recognition. Heteronuclear NMR relaxation is sensitive to motions in the relevant time scales and could provide valuable information on the dynamics of multi-domain proteins. However, the standard analysis based on the separation of global tumbling and fast local motions is no longer valid for multi-domain proteins undergoing internal motions involving complete domains and that take place on the same time scale than the overall motion.The complexity of the motions experienced even for the simplest two-domain proteins are difficult to capture with simple extensions of the classical Lipari-Szabo approach. Hydrodynamic effects are expected to dominate the motion of the individual globular domains, as well as that of the complete protein. Using Pin1 as a test case, we have simulated its motion at the microsecond time scale, at a reasonable computational expense, using Brownian Dynamic simulations on simplified models. The resulting trajectories provide insight on the interplay between global and inter-domain motion and can be analyzed using the recently published method of isotropic Reorientational Mode Dynamics which offer a way of calculating their contribution to heteronuclear relaxation rates. The analysis of trajectories computed with Pin1 models of different flexibility provides a general framework to understand the dynamics of multi-domain proteins and explains some of the observed features in the relaxation rate profile of free Pin1.     

Characterization of azido-NAD to assess its potential as a two-dimensional infrared probe of enzyme dynamics

Enzyme active-site dynamics at femtosecond to picosecond time scales are of great biochemical importance, but remain relatively unexplored due to the lack of appropriate analytical methods. Two-dimensional infrared (2D IR) spectroscopy is one of the few methods that can examine chemical biological motions at this time scale, but all the IR probes used so far were specific to a few unique enzymes. The lack of IR probes of broader specificity is a major limitation to further 2D IR studies of enzyme dynamics. Here we describe the synthesis of a general IR probe for nicotinamide-dependent enzymes. This azido analog of the ubiquitous cofactor nicotinamide adenine dinucleotide is found to be stable and bind to several dehydrogenases with dissociation constants similar to that for the native cofactor. The infrared absorption spectra of this probe bound to several enzymes indicate that it has significant potential as a 2D IR probe to investigate femtosecond dynamics of nicotinamide-dependent enzymes.     

Multi‐constraint computational design suggests that native sequences of germline antibody H3 loops are nearly optimal for conformational flexibility

The limited size of the germline antibody repertoire has to recognize a far larger number of potential antigens. The ability of a single antibody to bind multiple ligands due to conformational flexibility in the antigen‐binding site can significantly enlarge the repertoire. Among the six complementarity determining regions (CDRs) that generally comprise the binding site, the CDR H3 loop is particularly variable. Computational protein design studies showed that predicted low energy sequences compatible with a given backbone structure often have considerable similarity to the corresponding native sequences of naturally occurring proteins, indicating that native protein sequences are close to optimal for their structures. Here, we take a step forward to determine whether conformational flexibility, believed to play a key functional role in germline antibodies, is also central in shaping their native sequence. In particular, we use a multi‐constraint computational design strategy, along with the Rosetta scoring function, to propose that the native sequences of CDR H3 loops from germline antibodies are nearly optimal for conformational flexibility. Moreover, we find that antibody maturation may lead to sequences with a higher degree of optimization for a single conformation, while disfavoring sequences that are intrinsically flexible. In addition, this computational strategy allows us to predict mutations in the CDR H3 loop to stabilize the antigen‐bound conformation, a computational mimic of affinity maturation, that may increase antigen binding affinity by preorganizing the antigen binding loop. In vivo affinity maturation data are consistent with our predictions. The method described here can be useful to design antibodies with higher selectivity and affinity by reducing conformational diversity. Proteins 2009. © 2008 Wiley‐Liss, Inc.     

Backbone dynamics of Tet repressor alpha8intersectionalpha9 loop

A set of single Trp mutants of class B Tet repressor (TetR), in which Trp residues are located from positions 159 to 167, has been engineered to investigate the dynamics of the loop joining the alpha-helices 8 and 9. The fluorescence anisotropy decay of most mutants can be described by the sum of three exponential components. The longest rotational correlation time, 30 ns at 10 degrees C, corresponds to the overall rotation of the protein. The shortest two components, on the subnanosecond and nanosecond time scale, are related to internal motions of the protein. The initial anisotropy, in the 0.16-0.22 range, indicates the existence of an additional ultrafast motion on the picosecond time scale. Examination of physical models for underlying motions indicates that librational motions of the Trp side chain within the rotameric chi(1) x chi(2) potential wells contribute to the picosecond depolarization process, whereas the subnanosecond and nanosecond depolarization processes are related to backbone dynamics. In the absence of inducer, the order parameters of these motions, about 0.90 and 0.80 for most positions, indicate limited flexibility of the loop backbone. Anhydrotetracycline binding to TetR induces an increased mobility of the loop on the nanosecond time scale. This suggests that entropic factors might play a role in the mechanism of allosteric transition.     

Human Germline Antibody Gene Segments Encode Polyspecific Antibodies

Structural flexibility in germline gene-encoded antibodies allows promiscuous binding to diverse antigens. The binding affinity and specificity for a particular epitope typically increase as antibody genes acquire somatic mutations in antigen-stimulated B cells. In this work, we investigated whether germline gene-encoded antibodies are optimal for polyspecificity by determining the basis for recognition of diverse antigens by antibodies encoded by three V_H gene segments. Panels of somatically mutated antibodies encoded by a common V_H gene, but each binding to a different antigen, were computationally redesigned to predict antibodies that could engage multiple antigens at once. The Rosetta multi-state design process predicted antibody sequences for the entire heavy chain variable region, including framework, CDR1, and CDR2 mutations. The predicted sequences matched the germline gene sequences to a remarkable degree, revealing by computational design the residues that are predicted to enable polyspecificity, i.e., binding of many unrelated antigens with a common sequence. The process thereby reverses antibody maturation in silico. In contrast, when designing antibodies to bind a single antigen, a sequence similar to that of the mature antibody sequence was returned, mimicking natural antibody maturation in silico. We demonstrated that the Rosetta computational design algorithm captures important aspects of antibody/antigen recognition. While the hypervariable region CDR3 often mediates much of the specificity of mature antibodies, we identified key positions in the V_H gene encoding CDR1, CDR2, and the immunoglobulin framework that are critical contributors for polyspecificity in germline antibodies. Computational design of antibodies capable of binding multiple antigens may allow the rational design of antibodies that retain polyspecificity for diverse epitope binding.     

13C relaxation studies of the DNA target sequence for hhai methyltransferase reveal unique motional properties

The goal of this work was to examine if sequence-dependent conformational flexibility in DNA plays a role in base extrusion, a common conformational change induced by many DNA-modifying enzymes. We studied the dynamics of the double-stranded DNA target of the HhaI methyltransferase by recording an extensive set of (13)C NMR relaxation parameters. We observe that the cytidine furanose rings experience fast (picosecond to nanosecond) motions that are not present in other nucleotides; the methylation site experiences particularly high mobility. We also observe that the bases of guanosine and cytidine residues within the HhaI recognition sequence GCGC experience motions on a much slower (1-100 micros) time scale. We compare these observations with previous solution and solid-state NMR studies of the EcoRI nuclease target sequence, and solid-state NMR studies of a similar HhaI target construct. While an increased mobility of cytidine furanose rings compared to those of other nucleotides is observed for both sequences, the slower motions are only observed in the HhaI target DNA. We propose that this inherent flexibility lowers the energetic barriers that must occur when the DNA binds to the HhaI methyltransferase and for extrusion of the cytidine prior to its methylation.     

Influence of heme post-translational modification and distal ligation on the backbone dynamics of a monomeric hemoglobin

The cyanobacterium Synechococcus sp. PCC 7002 uses a hemoglobin of the truncated lineage (GlbN) in the detoxification of reactive species generated in the assimilation of nitrate. In view of a sensing or enzymatic role, several states of GlbN are of interest with respect to its structure-activity relationship. Nuclear magnetic resonance spectroscopy was applied to compare the structure and backbone dynamics of six GlbN forms differing in their oxidation state [Fe(II) or Fe(III)], distal ligand to the iron (histidine, carbon monoxide, or cyanide), or heme post-translational modification (b heme or covalently attached heme). Structural properties were assessed with pseudocontact shift calculations. (15)N relaxation data were analyzed by reduced spectral density mapping (picosecond to nanosecond motions) and by inspection of elevated R(2) values (microsecond to millisecond motions). On the picosecond to nanosecond time scale, GlbN exhibited little flexibility and was unresponsive to the differences among the various forms. Regions of slightly higher mobility were the CE turn, the EF loop, and the H-H' kink. In contrast, fluctuations on the microsecond to millisecond time scale depended on the form. Cyanide binding to the ferric state did not enhance motions, whereas reduction to the ferrous bis-histidine state resulted in elevated R(2) values for several amides. This response was attributed, at least in part, to a weakening of the distal histidine coordination. Carbon monoxide binding quenched some of these fluctuations. The results emphasized the role of the distal ligand in dictating backbone flexibility and illustrated the multiple ways in which motions are controlled by the hemoglobin fold.     

Backbone dynamics of the ribonuclease binase active site area using multinuclear ((15)N and (13)CO) NMR relaxation and computational molecular dynamics

The nano-pico second backbone dynamics of the ribonuclease binase, homologous to barnase, is investigated with (15)N, (13)C NMR relaxation at 11.74 and 18.78 T and with a 1.1 ns molecular dynamics simulation. The data are compared with the temperature factors reported for the X-ray structure of this enzyme. The molecular dynamics and X-ray data correspond well and predict motions in the loops 56-61 and 99-104 that contain residues that specifically recognize substrate and are catalytic (His101), respectively. In contrast, the (15)N relaxation data indicate that these loops are mostly ordered at the nano-pico second time scale. Nano-pico second motions in the recognition loop 56-61 are evident from (13)CO-(13)C cross relaxation data, but the mobility of the catalytic loop 99-104 is not detected by (13)CO cross relaxation either. From the results of this and previous work [Wang, L., Pang, Y., Holder, T., Brender, J. R., Kurochkin, A., and Zuiderweg, E. R. P. (2001) Proc. Natl. Acad. Sci. U.S.A., 98, 7684-7689], the following dynamical characterization of the active site area of binase emerges: a beta sheet, rigid at all probed time scales, supports the catalytic residue Glu 72. Both substrate-encapsulating loops are mobile on both fast and slow time scales, but the fast motions of the loop which contains the other catalytic residue, His 101, as predicted by B-factors and computational molecular dynamics is not detected by NMR relaxation. This work strongly argues for the use of several measures in the study of protein dynamics.     

Exploration of specificity in germline monoclonal antibody recognition of a range of natural and synthetic epitopes

To explore the molecular basis of antigen recognition by germline antibodies, we have determined to high resolution the structures of the near-germline monoclonal antibody S25-2 in complex with seven distinct carbohydrate antigens based on the bacterial sugar 3-deoxy-α-d-manno-oct-2-ulosonic acid (Kdo). In contrast to previous findings, the inherited germline Kdo monosaccharide binding site is not restricted to this bacterial sugar but is able to accommodate an array of substitutions and chemical modifications of Kdo, including naturally occurring antigens containing the related monosaccharide d-glycero-α-d-talo-oct-2-ulosonic acid as well as nonterminal Kdo residues. However, we show by surface plasmon resonance and ELISA how antibody S25-2 specificity is so dependent on the context in which the antigen is presented that a free disaccharide displays strong binding while the same lipid-A-bound disaccharide does not bind. These structures provide insight into how inherited germline genes code for immunoglobulins of limited flexibility that are capable of binding a range of epitopes from which affinity-matured antibodies are generated.     

Dynamics of cellular retinoic acid binding protein I on multiple time scales with implications for ligand binding

Cellular retinoic acid binding protein I (CRABPI) belongs to the family of intracellular lipid binding proteins (iLBPs), all of which bind a hydrophobic ligand within an internal cavity. The structures of several iLBPs reveal minimal structural differences between the apo (ligand-free) and holo (ligand-bound) forms, suggesting that dynamics must play an important role in the ligand recognition and binding processes. Here, a variety of nuclear magnetic resonance (NMR) spectroscopy methods were used to systematically study the dynamics of both apo and holo CRABPI at various time scales. Translational and rotational diffusion constant measurements were used to study the overall motions of the proteins. Both apo and holo forms of CRABPI tend to self-associate at high (1.2 mM) concentrations, while at low concentrations (0.2 mM), they are predominantly monomeric. Rapid amide exchange rate and laboratory frame relaxation rate measurements at two spectrometer field strengths (500 and 600 MHz) were used to probe the internal motions of the individual residues. Several residues in the apo form, notably within the ligand recognition region, exhibit millisecond time scale motions that are significantly arrested in the holo form. In contrast, no significant differences in the high-frequency motions were observed between the two forms. These results provide direct experimental evidence for dynamics-induced ligand recognition and binding at a specifically defined time scale. They also exemplify the importance of dynamics in providing a more comprehensive understanding of how a protein functions.     

Characterization of the internal motions of a chimeric protein by 13C NMR highlights the important dynamic consequences of the engineering on a millisecond time scale.

By transferring the central curaremimetic beta hairpin of the snake toxin alpha into the scaffold of the scorpion charybdotoxin, a chimeric protein was constructed that reproduced the three-dimensional structure and partially reproduced the function of the parent beta hairpin, without perturbing the three-dimensional structure of the scaffold [1]. Picosecond to hour time scale motions of charybdotoxin and the engineered protein were observed, in order to evaluate the dynamic consequences of the six deletions and eight mutations differentiating the two molecules. The chimeric protein dynamics were also compared to that of toxin alpha, in order to examine the beta hairpin motions in both structural contexts. Thus, 13C R1, R1rho and 1H-->13C nOe were measured for all the CalphaHalpha and threonine CbetaHbeta vectors. As the proteins were not labeled, accordion techniques combined to coherence selection by pulsed field gradients and preservation of magnetization following equivalent pathways were used to considerably reduce the spectrometer time needed. On one hand, we observed that the chimeric protein and charybdotoxin are subjected to similar picosecond to nanosecond time scale motions except around the modified beta sheet region. The chimeric protein also exhibits an additional millisecond time scale motion on its whole sequence, and its beta structure is less stable on a minute to hour time scale. On the other hand, when the beta hairpin dynamics is compared in two different structural contexts, i.e. in the chimeric protein and the curaremimetic toxin alpha, the picosecond to nanosecond time scale motions are fairly conserved. However, the microsecond to millisecond time scale motions are different on most of the beta hairpin sequence, and the beta sheet seems more stable in toxin alpha than in the chimera. The slower microsecond to hour time scale motions seem to be extremely sensitive to the structural context, and thus poorly transferred from one protein to another.     

Identification of peptide mimotopes for the fluorescein hapten binding of monoclonal antibody B13‐DE1

Using 6mer and 12mer phage peptide libraries three unique phage clones were identified which specifically bind to a monoclonal anti‐FITC antibody, B13‐DE1. The two 6mer and one 12mer peptide insert sequences are clearly related to each other and contain a high proportion of hydrophobic amino acids. The peptides are bound by the antibody combining site of B13‐DE1 probably in a similar manner to FITC and represent therefore true peptidic mimics of the fluorescein hapten. No reactivity of the peptides could be demonstrated with another monoclonal anti‐fluorescein antibody or with polyclonal anti‐fluorescein antibodies. Immunization of mice with the peptides resulted in the production of antibodies cross‐reacting with all peptides but not with fluorescein. The results show that phage peptide libraries can be used to isolate mimotope peptides which can mimic low molecular weight structures seen by a specific antibody and probably other recognition molecules. Copyright © 1999 John Wiley & Sons, Ltd.     

Coarse graining the dynamics of nano-confined solutes: the case of ions in clays

We investigate the possibility of describing by a continuous solvent model the dynamics of solutes confined down to the molecular scale. We derive a generalised Langevin equation (GLE) for the generic motion of a solute in an external potential using the Mori–Zwanzig formalism. We then compute the corresponding memory function from molecular simulations, in the case of cesium ions confined in the interlayer porosity of montmorillonite clays, with a very low water content (only six solvent molecules per ion). Previous attempts to describe the dynamics of cesium in this system by a simple Langevin equation were unsuccessful. The purpose of this work is not to carry out GLE simulations using the memory function from molecular simulations, but rather to analyse the separation of time scales between the confined ions and solvent. We show that such a separation is not achieved and discuss the relative contribution of the ion–surface, ion–solvent and ion–ion interactions to the dynamics. On the picosecond time scale, the ion oscillates in a surface-and-solvent cage, which relaxes on much longer time scales extending to several nanoseconds. The resulting overall dynamics resembles that of glasses or diffusion inside a solid by site-to-site hopping.     

Structural equilibrium fluctuations in mesophilic and thermophilic alpha-amylase

By comparing a mesophilic alpha-amylase with its thermophilic homolog, we investigated the relationship between thermal stability and internal equilibrium fluctuations. Fourier transform infrared spectroscopy monitoring hydrogen/deuterium (H/D) exchange kinetics and incoherent neutron scattering measuring picosecond dynamics were used to study dynamic features of the folded state at room temperature. Fairly similar rates of slowly exchanging amide protons indicate about the same free energy of stabilization DeltaG(stab) for both enzymes at room temperature. With respect to motions on shorter time scales, the thermophilic enzyme is characterized by an unexpected higher structural flexibility as compared to the mesophilic counterpart. In particular, the picosecond dynamics revealed a higher degree of conformational freedom for the thermophilic alpha-amylase. The mechanism proposed for increasing thermal stability in the present case is characterized by entropic stabilization and by flattening of the curvature of DeltaG(stab) as a function of temperature.     

Dynamics of a de novo designed three-helix bundle protein studied by 15N, 13C, and 2H NMR relaxation methods

Understanding how the amino acid sequence of a polypeptide chain specifies a unique, functional three-dimensional structure remains an important goal, especially in the context of the emerging discipline of de novo protein design. Alpha3D is a single chain protein of 73 amino acids resulting from a de novo design effort. Previous solution nuclear magnetic resonance studies of alpha3D confirm that the protein adopts the designed structure of a three-helix bundle. Furthermore, alpha3D has been previously shown to possess all of the major thermodynamic and structural characteristics of natural proteins, though it shares no sequence homology to any protein sequence in the database. In this work, the backbone and side-chain dynamics of alpha3D were investigated using 15N, 13C, and 2H nuclear magnetic resonance relaxation methods with the aim of assessing the character of the internal motions of this native-like protein of de novo design. At the backbone level, both 15N and 13C(alpha) relaxation studies indicate highly restrictive motion on the picosecond to nanosecond time scale in the alpha-helical regions of alpha3D, with increasing mobility at the ends of the alpha-helices and in the two loop regions. This is largely consistent with what is seen in proteins of natural origin. Overall, the view provided by both 2H and 13C methyl relaxation methods suggest that the side chains of alpha3D are more dynamic compared to natural proteins. Regions of relative flexibility bound clusters of rigid methyl-bearing side-chain groups that are interspersed with aromatic and beta-branched amino acids. The time scale of motions associated with methyl-bearing side chains of alpha3D are significantly longer than that seen in natural proteins. These results indicate that the strategies underlying the design of alpha3D have largely, but not completely, captured both the structural and dynamic character of natural proteins.     

A view of dynamics changes in the molten globule-native folding step by quasielastic neutron scattering

In order to understand the changes in protein dynamics that occur in the final stages of protein folding, we have used neutron scattering to probe the differences between a protein in its folded state and the molten globule states. The internal dynamics of bovine alpha-lactalbumin (BLA) and its molten globules (MBLA) have been examined using incoherent, quasielastic neutron scattering (IQNS). The IQNS results show length scale dependent, pico-second dynamics changes on length scales from 3.3 to 60 A studied. On shorter-length scales, the non-exchangeable protons undergo jump motions over potential barriers, as those involved in side-chain rotamer changes. The mean potential barrier to local jump motions is higher in BLA than in MBLA, as might be expected. On longer length scales, the protons undergo spatially restricted diffusive motions with the diffusive motions being more restricted in BLA than in MBLA. Both BLA and MBLA have similar mean square amplitudes of high frequency motions comparable to the chemical bond vibrational motions. Bond vibrational motions thus do not change significantly upon folding. Interestingly, the quasielastic scattering intensities show pronounced maxima for both BLA and MBLA, suggesting that "clusters" of atoms are moving collectively within the proteins on picosecond time scales. The correlation length, or "the cluster size", of such atom clusters moving collectively is dramatically reduced in the molten globules with the correlation length being 6.9 A in MBLA shorter than that of 18 A in BLA. Such collective motions may be important for the stability of the folded state, and may influence the protein folding pathways from the molten globules.