Motions in proteins

Studies on the location of aromatic amino acids in alpha-crystallin

The locations of tryptophan residues in alpha-crystallin and homopolymers constructed from the alpha A- and alpha B-chains were examined by comparing their fluorescence emission properties and their accessibilities to quenchers. Two classes of tryptophan could be distinguished on the basis of differences in their spectral characteristics, fluorescence decay lifetimes, quenching with acrylamide and exposure by increasing concentrations of urea. Polarization measurements indicated that the tryptophan residues were associated with flexible segments of the polypeptide chains. The two classes could be assigned, one to Trp-9 (in both A- and B-chains) which is in an hydrophobic environment, and one to Trp-60 (B-chain) which appeared to be nearer the surface of the aggregate. No evidence was found for residues inaccessible to the quenchers. An apparent partition coefficient of 40 was obtained for the association of acrylamide with the protein. From temperature-dependence studies, it was concluded that there was a significant energy barrier to the penetration of acrylamide into the protein matrix (Ea = 5.8 kcal/mol) and that entry of the quencher was through channels produced by temporary disruption of the matrix (delta s = 1.5 eu). The phenolic side-chains of tyrosine residues in several different alpha-crystallins were found to ionize with pK values above pH 11, indicative of residues highly shielded from the solvent. Comparisons of polypeptide sequences, together with tyrosine fluorescence emission data and the pK values, permitted a tentative assignment of residue locations. All of the data are consistent with a possible micelle-like structure for alpha-crystallin but not with a layered structure.

Halophilic class I aldolase and glyceraldehyde-3-phosphate dehydrogenase: some salt-dependent structural features

Aldolase and glyceraldehyde-3-phosphate dehydrogenase from the extremely halophilic archaebacterium Haloarcula vallismortis are stable only in high concentrations of KCl present within the physiological environment. Data concerning the structural changes in the two enzymes as a result of lowering of salt concentration and changes in pH were obtained by monitoring the intrinsic protein fluorescence in the presence of quenchers. When the KCl concentrations were lowered below 2 M or in the presence of 6 M guanidine hydrochloride, the emission maximum shifted to a longer wavelength, indicating enhanced exposure of tryptophyl residues to the solvent. The spectral characteristics of the two proteins in guanidine hydrochloride and 0.4 M KCl were identical. However, these denatured states appear to be different than those observed after acid denaturation. Further perturbation of fluorescence was observed due to I-, and application of the Stern-Volmer law showed that the total fluorescence was available to the quenchers only in 0.4 M KCl solutions. The unfolding of proteins in 0.4 M KCl was a gradual process which was accompanied by a time-dependent loss in enzyme activity. The activity loss was complete within 30 min for aldolase whereas in the case of GAPDH nearly 3 h was required for the destruction of activity. For both enzymes, inactivation and protein denaturation were strongly correlated. The data on activity and thermostability measurements of the two enzymes in varying concentrations of KCl and potassium phosphate revealed that though both proteins are halophilic, the forces in the maintenance of their stability could be different.(ABSTRACT TRUNCATED AT 250 WORDS)

Backbone dynamics of the Bacillus subtilis glucose permease IIA domain determined from 15N NMR relaxation measurements

The backbone dynamics of the uniformly 15N-labeled IIA domain of the glucose permease of Bacillus subtilis have been characterized using inverse-detected two-dimensional 1H-15N NMR spectroscopy. Longitudinal (T1) and transverse (T2) 15N relaxation time constants and steady-state (1H)-15N NOEs were measured, at a spectrometer proton frequency of 500 MHz, for 137 (91%) of the 151 protonated backbone nitrogens. These data were analyzed by using a model-free dynamics formalism to determine the generalized order parameter (S2), the effective correlation time for internal motions (tau e), and 15N exchange broadening contributions (Rex) for each residue, as well as the overall molecular rotational correlation time (tau m). The T1 and T2 values for most residues were in the ranges 0.45-0.55 and 0.11-0.15 s, respectively; however, a small number of residues exhibited significantly slower relaxation. Similarly, (1H)-15N NOE values for most residues were in the range 0.72-0.80, but a few residues had much smaller positive NOEs and some exhibited negative NOEs. The molecular rotational correlation time was 6.24 +/- 0.01 ns; most residues had order parameters in the range 0.75-0.90 and tau e values of less than ca. 25 ps. Residues found to be more mobile than the average were concentrated in three areas: the N-terminal residues (1-13), which were observed to be highly disordered; the loop from P25 to D41, the apex of which is situated adjacent to the active site and may have a role in binding to other proteins; and the region from A146 to S149. All mobile residues occurred in regions close to termini, in loops, or in irregular secondary structure.

Internal motional averaging and three-dimensional structure determination by nuclear magnetic resonance

Dynamic averaging effects from internal motions on interproton distances estimated from nuclear Overhauser effects (NOE) are determined by using a molecular dynamics simulation of lysozyme. Generalized order parameters measuring angular averaging and radial averaging parameters are calculated. The product of these two parameters describes the full averaging effects on cross-relaxation. Analysis of 2778 non-methyl NOE interactions from the protein interior and surface indicates that distances estimated by assuming a rigid molecule have less than 10% error for 89% of the NOE interactions. However, analysis of 1854 methyl interactions found that only 68% of the distances estimated from cross-relaxation rates would have less than 10% error. Qualitative evaluation of distances according to strong, medium and weak NOE intensities, when used to define only the upper bound for interproton separation, would misassign less than 1% of the distance constraints because of motional averaging. Internal motions do not obscure the identification of secondary structure, although some instances of significant averaging effects were found for interactions in alpha-helical regions. Interresidue NOEs for amino acids more than three residues apart in the primary sequence are more extensively averaged than intraresidue or short-range interresidue NOEs. Intraresidue interactions exhibit a greater degree of angular averaging than those involving interresidue proton pairs. An internal motion does not equally affect all NOE interactions for a particular proton. Thus, incorporation of averaging parameters in nuclear magnetic resonance structure determination procedures must be made on a proton-pair-wise basis. On the basis of the motional averaging results, particular fixed-distance proton pairs in proteins are suggested for use as distance references. A small percentage of NOE pairs localized to three regions of the protein exhibit extreme averaging effects from internal motions. The regions and types of motions involved are described.

Conformational dynamics and solvent viscosity effects in carboxypeptidase-A-catalyzed benzoylglycylphenyllactate hydrolysis

We have used a new approach to the dynamics of hydrolytic metalloenzyme catalysis based on investigations of both external solvent viscosity effects and kinetic 2H isotope effects. The former reflects solvent and protein dynamics, and the nuclear reorganization distribution among damped protein motion and intramolecular friction-free nuclear motion. The isotope effect represents proton tunnelling and reorganization in the hydrogen bond network around the active site. We illustrate the approach by new spectrophotometric and pH-titration data for carboxypeptidase-A-catalyzed benzoylglycyl-L-phenyllactate hydrolysis. This substrate exhibits both a significant inverse fractional power law viscosity dependence over wide ranges controlled by glycerol and sucrose, and a kinetic 2H isotope effect of 1.65. The analogous benzoylglycylphenylalanine hydrolysis has a smaller isotope effect (1.3) and no viscosity dependence. Viscosity variation has no effect on the CD spectra in the 180-240-nm range. In terms of stochastic chemical rate theory, the data correspond to an enzyme-peptide substrate complex with a 'tight' structure protected from the solvent. In comparison, the enzyme-ester substrate complex is 'softer', strongly coupled to the solvent, and the rate-determining step is accompanied by proton transfer or by substantial reorganization in the hydrogen bonds near the active site.

The effects of the molecular environment on the kinetics of catalase reaction and its relevance to cell aging

Catalase-mediated oxygen evolution from H(2)O(2) was measured with a Clark electrode in a closed, buffered aqueous medium in a properly designed experimental model. Addition of various molecules (bovine serum albumin, egg albumin, dextrane, starch and PEG-8000) in increasing concentrations up to greater than 20% by volume decreased the rate of oxygen evolution in a linear fashion, although there were some quantitative differences between the effects of various components. The rotational correlation time of the spin label TMPN in BSA and egg albumin increased considerably in both protein solutions indicating a very significant increase of the microviscosity of these media. The results are interpreted in terms of the molecular enzyme kinetic model (MEKM) predicting that apart from microviscosity, the average mass distribution and the distance of the nearest two lattice points formed by components of the medium are also inversely proportional to both the forward and backward rate constants of the formation of the enzyme-substrate complex. These observations, on one hand, lend further experimental support to the MEKM and on the other one, are consistent with the membrane hypothesis of aging which explains the loss of cell functions by an age-dependent increase of intracellular density.

Effect of salt and membrane fluidity on fluorophore motions of a gramicidin C derivative

We have used fluorescence spectroscopy to investigate the effect of salt and membrane fluidity on the rotational motion of a 5-(dimethylamino)naphthalene-1-sulfonyl (dansyl) derivative of gramicidin C (dansyl-gC) in dimyristoylphosphatidylcholine bilayers, under conditions where the peptide is a formyl-NH to formyl-NH terminal dimer, and in octyl glucoside micelles, where the peptide is an intertwined helical dimer. Energy-transfer experiments showed no changes in either conformation or dimer aggregation under the conditions explored here (15-40 degrees C, 1-350 bar, 0-0.33 M Mg2+, and 0-1 M Na+). The addition of permeable (Na+) or nonpermeable (Mg2+) ions did not affect the temperature or pressure behavior of dansyl rotation. However, fluorescence lifetime measurements indicated an increase in solvent accessibility in the presence of sodium. In bilayers, the temperature dependence of the fluorescence polarization and lifetime shows strong interactions between the dansyl residue and the peptide, and at no time did the dansyl motions become solvent controlled as has been observed for aqueous solvent peptides [Scarlata, S. F., Rholam, M., & Weber, G. (1984) Biochemistry 23, 6789]. In micelles, the change in rotational motion with temperature followed solvent expansion, showing that in this case the dansyl residue does not associate extensively with the peptide. Our results indicate that because of the extensive coupling between the dansyl residue and the rest of the peptide, membrane fluidity does not play a major role in controlling side-chain motions.

A nuclear magnetic resonance study of the hydrogen-exchange behaviour of lysozyme in crystals and solution

Amide hydrogen/deuterium exchange behaviour has been studied for all of the peptide amides of hen lysozyme by means of two-dimensional n.m.r. spectroscopy. The amides have been grouped into four categories on the basis of their rates of exchange in solution at pH 4.2 and 7.5. The distribution of the amides into the different categories has been examined in the light of the crystallographic structural information, considering the type of secondary structure, the nature of hydrogen bonding and the distance from the protein surface. None of these features was found to determine uniquely the pattern of hydrogen exchange rates within the protein. The exchange behaviour of the individual amides could, however, in general be rationalized by a combination of these features. Hydrogen exchange was also monitored in both tetragonal and triclinic crystals of lysozyme, by allowing exchange to take place in the crystals prior to dissolution and recording of n.m.r. spectra under conditions where further exchange was minimized. This enabled direct comparison to be made of the exchange behaviour in the crystals and solution. A reduction in exchange rate was observed in the crystalline state relative to solution for a substantial number of amides and distinct differences between exchange in the different crystals could be observed. These differences between the solution and the different crystal states do not, however, correlate in a simple manner with proximity to intermolecular contacts in the crystals. However, the existence of these contacts, which are on the surface of the protein molecule, have a profound effect on the exchange of amides in the interior of the protein. The results indicate that the spectrum of fluctuations giving rise to hydrogen exchange may be significantly altered by the intermolecular interactions present within the crystalline state.

Substrate specificity of solvent viscosity effects in carboxypeptidase A catalyzed peptide hydrolysis

We have investigated the viscosity of carboxypeptidase A catalyzed Bz-Gly-Phe hydrolysis at pH 7.5 (Tris) and 0.5 mol.l-1 NaCl over the range 10-100 mp, varied by addition of glycerol or sucrose. In contrast to previous reports of strong viscosity effects on the corresponding Cbz-Ala-Ala-Ala hydrolysis, both the catalytic constant and the Michaelis constant are virtually independent of viscosity over the 10-fold range investigated. Furthermore, the CD spectra of carboxypeptidase A in the high-viscosity media point to no change in the alpha-helix and beta-sheet structure in these media. The data are compatible either with a compacter, more rigid enzyme-substrate structure or with a more prominent role of intramolecular nuclear reorganization compared to protein reorganization for Bz-Gly-Phe than for Cbz-Ala-Ala-Ala. These views can be given a preciser frame in terms of stochastic chemical rate theory.

Anisotropy decays of single tryptophan proteins measured by GHz frequency-domain fluorometry with collisional quenching

We used harmonic-content frequency-domain fluorometry to determine the anisotropy decays of a variety of single tryptophan peptides and proteins. Resolution of the rapid and complex anisotropy decays was enhanced by global analysis of the data measured in the presence of quenching by either oxygen or acrylamide. For each protein, and for each quencher, data were obtained at four to six quencher concentrations, and the data analyzed globally to recover the anisotropy decay. The decrease in decay times produced by quenching allows measurements to an upper frequency limit of 2 GHz. The chosen proteins provided a range of exposures of the tryptophan residues to the aqueous phase, these being ACTH, monellin, Staphylococcus nuclease and ribonuclease T1, in order of decreasing exposure. Examination of indole and several small peptides demonstrates the resolution limitations of the measurements; a correlation time of 12 ps was measured for indole in methanol at 40 degrees C. Comparison of the anisotropy decays of gly-trp-gly with leu-trp-leu revealed stearic effects of the larger leucine side chains on the indole ring. The anisotropy decay of gly-trp-gly revealed a 40 ps component for the indole side chain, which was resolved from the overall 150 ps correlation time of the tripeptide. Only the longer correlation time was observed for leu-trp-leu. With the exception of ribonuclease T1, each of the proteins displayed a subnanosecond component in the anisotropy decay which we assign to independent motions of the tryptophan residues. For example, Staphylococcus nuclease and monellin displayed segmental tryptophan motions with correlation times of 80 and 275 ps, respectively. The amplitudes of the rapid components increased with increasing exposure to the aqueous phase. These highly resolved anisotropy decays for proteins of known structure are suitable for comparison with molecular dynamic simulations.

Structure and fluctuations of bacteriophage T4 glutaredoxin modelled by molecular dynamics

A 120 ps molecular dynamics (MD) trajectory for bacteriophage T4 glutaredoxin was calculated including non-inertial solvent effects. The potential energy attains an equilibrated regime after the first 20 ps. The r.m.s. difference of all non-hydrogen atoms between X-ray and average MD structures for the regular secondary structure is 0.99A which shows that the MD simulation reproduces the essentials of the structure with high accuracy. Loop displacements are detected, shown by the larger full structure all non-hydrogen atom r.m.s. difference of 1.2A. The fluctuation pattern derived from MD agrees fairly well with that derived from X-ray isotropic temperature factors. The active site is a stable structural region in this MD modellization. Structural changes are put in context with the protein's function.

Calculation of conformational ensembles from potentials of mean force. An approach to the knowledge-based prediction of local structures in globular proteins

We present a prototype of a new approach to the folding problem of polypeptide chains. This approach is based on the analysis of known protein structures. It derives the energy potentials for the atomic interactions of all amino acid residue pairs as a function of the distance between the involved atoms. These potentials are then used to calculate the energies of all conformations that exist in the data base with respect to a given sequence. Then, by using only the most stable conformations, clusters of the most probable conformations for the given sequence are obtained. To discuss the results properly we introduce a new classification of segments based on their conformational stability. Special care is taken to allow for sparse data sets. The use of the method is demonstrated in the discussion of the identical oligopeptide sequences found in different conformations in unrelated proteins. VNTFV, for example, adopts a beta-strand in ribonuclease but it is found in an alpha-helical conformation in erythrocruorin. In the case of VNTFV the ensemble obtained consists of a single cluster of beta-strand conformations, indicating that this may be the preferred conformation for the pentapeptide. When the flanking residues are included in the calculation the hepapeptide P-VNTFV-H (ribonuclease) again yields an ensemble of beta-strands. However, in the ensemble of D-VNTFV-A (erythrocruorin) the major cluster is of alpha-helical type. In the present study we concentrate on the local aspects of protein conformations. However, the theory presented is quite general and not restricted to oligopeptides. We indicate extensions of the approach to the calculation of global conformations of proteins as well as conceivable applications to a number of molecular systems.

On the role of intracellular physicochemistry in quantitative gene expression during aging and the effect of centrophenoxine. A review

The turnover of proteins in biological systems is due mostly to an ever-occurring damaging (cross-linking) effect of the OH. free radicals. The replacement of the damaged proteins requires a continuous gene expression. A key issue of experimental gerontology is why the gene expression maintains the fidelity but loses the speed during aging. The membrane hypothesis of aging (MHA) proposes a cellular mechanism based on the fact that the more compact cellular structures (e.g., membranes) are damaged faster than the more diluted ones (e.g., cytosol). In addition, the cell membrane is exposed also to the residual heat-induced damage deriving from a frequent discharge of its electric polarity. Therefore, one can assume that even an extremely small incompleteness of the replacement of the damaged membrane components per turnover cycle may result in an error accumulation, which may be responsible first for the inhibition of growth, then for aging of cells. In agreement with this hypothesis, neurons display a life-long, gradual loss of the passive potassium permeability of the cell membrane, resulting in a continuous dehydration of the intracellular mass, i.e., an increase of physical density. Theory and experimental models show that this latter process causes a slowing down of all enzyme functions including those realizing the gene expression and the elimination of the damaged components. Increase of the dry mass content of cells and tissues is an obligatory process for maturation; however, later on it leads to aging. The known nootropic effects of centrophenoxine (CPH) can be interpreted on the basis of the OH. radical scavenger properties of dimethylaminoethanol (DMAE) which is incorporated in the neuronal membranes of the brain in form of phosphatidyl-DMAE. The protective effects of CPH (and of similar, newly synthesized other drugs) on the membrane components can slow down the age-dependent deteriorations of the intracellular physicochemistry, in agreement with the predictions of the MHA.

How the ear's works work

The senses of hearing and equilibrium depend on sensory receptors called hair cells which can detect motions of atomic dimensions and respond more than 100,000 times a second. Biophysical studies suggest that mechanical forces control the opening and closing of transduction channels by acting through elastic components in each hair cell's mechanoreceptive hair bundle. Other ion channels, as well as the mechanical and hydrodynamic properties of hair bundles, tune individual hair cells to particular frequencies of stimulation.

Applications of ultrafast laser spectroscopy for the study of biological systems

The discovery of mode-locked laser operation now nearly two decades ago has started a development which enables researchers to probe the dynamics of ultrafast physical and chemical processes at the molecular level on shorter and shorter time scales. Naturally the first applications were in the fields of photophysics and photochemistry where it was then possible for the first time to probe electronic and vibrational relaxation processes on a sub-nanosecond timescale. The development went from lasers producing pulses of many picoseconds to the shortest pulses which are at present just a few femtoseconds long. Soon after their discovery ultrashort pulses were applied also to biological systems which has revealed a wealth of information contributing to our understanding of a broadrange of biological processes on the molecular level.It is the aim of this review to discuss the recent advances and point out some future trends in the study of ultrafast processes in biological systems using laser techniques. The emphasis will be mainly on new results obtained during the last 5 or 6 years. The term ultrafast means that I shall restrict myself to sub-nanosecond processes with a few exceptions.

Studies on the biotin-binding sites of avidin and streptavidin. A chemically induced dynamic nuclear polarization investigation of the status of tyrosine residues

We applied the protein photochemically induced dynamic nuclear polarization (photo-c.i.d.n.p.) method to explore the conformation of the side chains of tyrosine, tryptophan and histidine residues in three biotin-binding proteins. The c.i.d.n.p. spectra of avidin, streptavidin and 'core' streptavidin were compared with those of their complexes with biotin and its derivatives. The data indicate that the single tyrosine residue (Tyr-33) of avidin is clearly inaccessible to the triplet flavin photo-c.i.d.n.p. probe. The same holds for all tryptophan and histidine side chains. Although the analogous Tyr-43 residue of streptavidin is also buried, at least three of the other tyrosine residues of this protein are exposed. The same conclusions apply to the truncated form of the protein, core streptavidin. As judged by the photo-c.i.d.n.p. results, complexing of avidin and streptavidin with biotin, N-epsilon-biotinyl-L-lysine (biocytin) or biotinyltyrosine has little or no effect on tyrosine accessibility in these proteins. Biotinyltyrosine can be used to probe the depth of the corresponding binding site. The accessibility of the tyrosine side chain of biotinyltyrosine in the complex demonstrates the exquisite fit of the biotin-binding cleft of avidin: only the biotin moiety appears to be accommodated, leaving the tyrosine side chain exposed.

Application of a molecular enzyme kinetic model for aging cells and tissues

According to the membrane hypothesis of aging (MHA), cellular senescence is attributable to a life-long, cross-linking action of oxygen-free radicals in the cell plasma membrane, resulting in a continuous decrease of the passive ion permeabilities. The consequent increase in the intracellular potassium content is accompanied by a considerable condensation of the intracellular mass (i.e., by loss of water). MHA suggested that an age-dependent increase in the physical density of the intracellular mass can underly the well-known age-dependent decreases of the macromolecular synthetic processes, the enzymic turnover rates, etc. MHA was partly based on a molecular enzyme kinetic model (MEKM) suggesting that environmental factors can substantially influence the enzyme catalysis and regulation through collisional coupling. However, the possible quantitative ranges of alterations in enzyme activities have not been estimated. This paper concludes, using principal features of the two models, that known age-dependent changes in the membrane lipid fluidity and intracellular density may result in even a 10-fold overall decrease in the enzyme activities (characterized by kcat and k-1) during the life.

Viscosity and transient solvent accessibility of Trp-63 in the native conformation of lysozyme

We have measured the rates of isotope exchange at the nitrogen of the indole ring of Trp-63 of lysozyme and of L-tryptophan as a function of solution viscosity. We have used two cosolvents, glycerol and ethylene glycol, to modify the relative viscosity. We have derived the appropriate kinetic equations for the alternative possibilities that the exchange takes place either in solution or in the intact protein matrix. Because we chose to study the proton-catalyzed exchange reaction, the rate of it is not expected to be diffusion-limited. We confirmed this by measuring the exchange from tryptophan. These results and the known effects of glycerol and ethylene glycol on the solvation of indole allow us to predict that if the exchange reaction takes place in a protein matrix the effects of the two cosolvents when compared under isoviscous conditions should be identical. This is what we find for Trp-63 in lysozyme at 15, 20 and 26 degrees C. The slope of the linear plot of log k vs. log relative viscosity is 0.6. This strongly supports a model for conformational fluctuations where transient solvation takes place without major changes in protein folding. The most interesting feature of our findings is the fact that a slow reaction admittedly not diffusion-limited shows, when taking place in a protein matrix, a linear dependence on solution viscosity. We suggest that what we observe is the effect of damping of movement of the side chain expressed as a change in the friction along the reaction coordinate in the corresponding phase space. The presence of such effects stresses the validity and usefulness of Kramers model of rate processes for reactions taking place in a protein matrix. Such behavior is predicted by several of the recently proposed general mechanisms of enzyme catalysis.

Strategy for obtaining non-native protein structures using antibody cross-reactions

Antibodies made to short peptides or to unfolded forms of proteins are often found to cross-react with intact proteins. These cross-reactions can be used to populate non-native protein conformations, possibly including protein folding intermediates, and the structures of the non-native conformations can be characterized using amide proton exchange and two-dimensional NMR.

Tryptophan sidechain dynamics in hydrophobic oligopeptides determined by use of 13C nuclear magnetic resonance spectroscopy

Two oligopeptides, t-boc-LAWAL-OMe and t-boc-LALALW-OMe, were synthesized for the purpose of examining the sidechain dynamics of the tryptophan residue in hydrophobic environments by 13C nuclear magnetic resonance and fluorescence spectroscopy. In both peptides, the tryptophan sidechain was greater than 95% enriched with 13C at the C delta 1 position. Spin-lattice relaxation time (T1) and steady-state nuclear Overhauser effect (NOE) data were obtained at 50.3 and 75.4 MHz for both peptides in CD3OD, and at 75.4 MHz for t-boc-LALALW-OMe in lysolecithin-D2O micelles. We have adapted the model-free approach of G. Lipari and A. Szabo (1982, J. Am. Chem. Soc. 104:4546) to interpret the 13C-NMR data. Computer-generated curves based on experimental data obtained at a single frequency demonstrate relationships between an effective correlation time for tryptophan sidechain motion (tau e), a generalized order parameter (sigma) describing the extent of motional restriction, and an overall correlation time for the peptide (tau m). Assuming predominantly dipolar relaxation, least-squares fits of the dual frequency relaxation data provide values for these parameters for both peptides. The contribution of chemical shift anisotropy (CSA), however, is also explicitly assessed in the data analysis, and is shown to perturb the predicted sigma, tau e, and tau m values and to decrease chi(2) values observed in nonlinear least-squares analysis of the data. Because of uncertainty in the contribution of CSA to the relaxation of the indole ring 13C delta 1 atom, nonlinear least-squares analysis of the relaxation data were performed with and without inclusion of a CSA term in the appropriate relaxation equations. Neglecting CSA, an overall peptide correlation time of 0.69 ns is predicted for t-boc-LAWAL-OMe in CD3OD at 20 degrees C compared with 1.28 ns for t-boc-LALALW-OMe. Given these tau m values and taking into account the effect of measurement error in the T1 and NOE data, the internal dynamics of the tryptophan residue of t-boc-LAWAL-OMe in this isotropic environment are described by a range of tau e values from 70 to 112 ps and sigma values between 0.22 and 0.36. Similarly, for t-boc-LALALW-OMe, 68 less than or equal to tau e less than or equal to 93 ps and 0.09 less than or equal to sigma less than or equal to 0.17. The Ch-terminal position of the tryptophan residue in the hexapeptide may account for its lower order parameter. In lysolecithin micelles, the model-free approach applied tot-boc-LALALW-OMe predicts a Te between 0.87 and 1.08 ns, and an order parameter range of 0.72-0.80, assuming an average Tm of 14 ns (Saunders, L. 1966. Biochim. Biophys. Acta. 125:70) for a typical peptide-micelle complex. In this case, measurement of only two 13C relaxation parameters at a single frequency yields sufficient information, plotted in the form of a composite T1-NOE solution curve, to constrain the allowed values of the model-free motional parameters within a relatively narrow range. The predicted range of eV and Te values for the peptide-micelle complex demonstrate that both the rate and spatial mobility of the indole moiety are markedly restrained in the anisotropic micelle environment relative to free methanol solution. Steady-state fluorescence anisotropy measurements made on the peptides dissolved in methanol or with synthetic lysolecithins in water were used to calculate apparent order parameters for tryptophan motion; these values agree well with order parameters calculated from 13C NMR data. The reported results are relevant to the issue of protein dynamic events occurring on the picosecond time scale predicted by molecular dynamics simulations.

Glucose-6-phosphate dehydrogenase from Leuconostoc mesenteroides: ligand-induced conformational changes

Glucose-6-phosphate dehydrogenase from Leuconostoc mesenteroides is inactivated by trypsin, chymotrypsin, pronase E, thermolysin, 4.0 M urea, and by heating to 49 degrees C. It is protected, to varying degrees, against all these forms of inactivation by glucose 6-phosphate, NAD+, and NADP+. When these ligands are present at 10 times their respective KD concentrations, protection by NAD+ or glucose 6-phosphate is substantially greater than protection by NADP+. A detailed analysis was undertaken of the protective effects of these ligands, at varying concentrations, on proteolysis of glucose-6-phosphate dehydrogenase by thermolysin. This study confirmed the above conclusion and permitted calculation of KD values for NAD+, NADP+, and glucose 6-phosphate that agree with such values determined by independent means. For NADP+, two KD values, 6.1 microM and 8.0 mM, can be derived, associated with protection against thermolysin by low and high NADP+ concentrations, respectively. The former value is in agreement with other determinations of KD and the latter value appears to represent binding of NADP+ to a second site which causes inhibition of catalysis. A Ki value of 10.5 mM for NADP+ was derived from inhibition studies. The principal conclusion from these studies is that NAD+ binding to L. mesenteroides glucose-6-phosphate dehydrogenase results in a larger global conformational change of the enzyme than does NADP+ binding. Presumably, a substantially larger proportion of the free energy of binding of NAD+, compared to NADP+, is used to alter the enzyme's conformation, as reflected in a much higher KD value. This may play an important role in enabling this dual nucleotide-specific dehydrogenase to accommodate either NAD+ or NADP+ at the same binding site.

Dynamics of protein conformational fluctuation in enzyme catalysis with special attention to proton transfers in serine proteinases

We have provided a quantum mechanical model for proteinase-catalyzed peptide, amide and ester hydrolysis. The model rests on electron and atom transfer theory, but incorporates the dynamics of conformational nuclear modes as a new element. The model is applied to acylation, but can straightaway be extended to deacylation, and is substantiated by recent structural and kinetic data for proteinase enzyme catalysis. The role of the conformational modes is found to be two-fold. First, the crystallographic distances for the proton transfers involved are far too large for direct transfer. His-57 mobility, handled stochastically, to bring the donor and acceptor groups within suitable reach, is therefore a crucial element of the theory. Secondly, the charge alignment in the Asp-102/His-57/tetrahedral intermediate system implies that the curvature of the potential surface along the conformational coordinates in this state is much lower than in the initial enzyme-substrate and final acyl states. A consequence of this is that the activation energy liberated after the first proton transfer is not dissipated, but stored in the conformational system and used in the second proton transfer step.

Internal motions in proteins and gating kinetics of ionic channels

Single-channel current recordings have revealed a complex kinetic behavior of ionic channels. Many channels exhibit closed-time distributions in which long waiting times occur with a much higher frequency than predicted by a simple exponential decay function. In this paper a model for opening-closing transitions that accounts for internal motions in the protein matrix is discussed. The model is based on the notion that the transition between a conductive and a nonconductive state of the channel represents a local process in the protein, such as the movement of a small segment of a peptide chain or the rotation of a single amino-acid residue. When the blocking group moves into the ion pathway, a structural defect is created consisting in a region of loose packing and/or poor hydrogen bonding. By rearrangements of neighboring groups, the defect may migrate within the protein matrix, carrying out a kind of random walk. Once the defect has moved away from the site where it was formed, a transition back to the open state of the channel is possible only when the defect has returned by chance to the original position. The kinetic properties of this model are analyzed by stochastic simulation of defect diffusion in a small domain of the protein. With a suitable choice of domain size and diffusion rate, the model is found to predict closed-time distributions that agree with experimental observations.

Tryptophan phosphorescence and the conformation of liver alcohol dehydrogenase in solution and in the crystalline state

Information on the effects of crystallization upon the structure of liver alcohol dehydrogenase from horse is obtained from a comparison of the phosphorescence properties of its tryptophan residues in solution and in the crystalline state. In the crystalline state the red shift in the phosphorescence spectrum of the solvent-exposed Trp-15 attests to a decreased polarity of its environment consistent with its shielding away from the aqueous solvent probably through its involvement in an intermolecular contact. On the other hand, the triplet-state lifetime of Trp-314 which is buried deeply in the coenzyme-binding domain demonstrates that the flexibility of this region of the macromolecule is unaffected by crystallization; a conclusion supported also by the similarity in the rate of oxygen quenching of its phosphorescence. Given that lattice constraints strongly inhibit large-scale conformational changes these results allow us to identify the average solution structure with the 'open' conformer determined crystallographically.

Structure and stability of thermophilic enzymes. Studies on thermolysin

The molecular mechanisms responsible for the unusual stability of enzymes isolated from thermophilic microorganisms are much more complex and subtle than was originally thought. In particular, a general mechanism cannot be proposed, since individual enzymes can be stabilized by specific molecular interactions and forces. The results of studies on thermophilic enzymes obtained in recent years in our laboratory will be summarized, with particular emphasis being placed on those obtained with thermolysin, a stable metalloendopeptidase isolated from Bacillus thermoproteolyticus. Fragmentation of thermolysin by limited proteolysis by added protease (subtilisin) or autolysis mediated by heat or the ion-chelating agent EDTA leads to quite selective peptide bond fissions, allowing isolation of 'nicked' thermolysin species. Correlation of the sites of proteolytic cleavage with the known three-dimensional structure of thermolysin allowed us to infer some of the key characteristics of the structure, folding, dynamics and stability of the thermolysin molecule. The potential utility of these and other studies on thermophilic enzymes in devising strategies for enhancing the stability of mesophilic enzymes using genetic engineering techniques is discussed.

Are spontaneous conformational interconversions a molecular basis for long-period oscillations in enzyme activity?

An unconventional hypothesis to the molecular basis of enzyme rhythms is that the intrinsic physical instability of the protein molecules which, in an aqueous medium, tend to move continuously from one conformational state to another could lead, in the population of enzyme molecules, to sizeable long-period oscillations in affinity for substrate and sensitivity to ligands and regulatory effects. To investigate this hypothesis, malate dehydrogenase was extracted and purified from leaves of the plant Kalanchoe blossfeldiana. The enzyme solutions were maintained under constant conditions and sampled at regular intervals for up to 40 or 70 h for measurements of activity as a function of substrate concentration, Km for oxaloacetic acid and sensitivity to the action of 2,3-butanedione, a modifier of active site arginyl residues. The results show that continuous slow oscillations in the catalytic capacity of the enzyme occur in all the extracts checked, together with fluctuations in Km. Apparent circadian periodicities were observed in accordance with previous data established during long run (100 h) experiments. The saturation curves for substrate showed multiple kinetic functions, with various pronounced intermediary plateaus and "bumps" depending on the time of sampling. Variation in the response to the effect of butanedione indicated fluctuation in the accessibility to the active site. Taken together, the results suggest that, under constant conditions, the enzyme in solution shifts continuously and reversibly between different configurations. This was confirmed by parallel studies on the proton-NMR spectrum of water aggregates in the enzyme solution and proton exchange rates.(ABSTRACT TRUNCATED AT 250 WORDS)

Fractal analysis of a voltage-dependent potassium channel from cultured mouse hippocampal neurons

The kinetics of ion channels have been widely modeled as a Markov process. In these models it is assumed that the channel protein has a small number of discrete conformational states and the kinetic rate constants connecting these states are constant. In the alternative fractal model the spontaneous fluctuations of the channel protein at many different time scales are represented by a kinetic rate constant k = At1-D, where A is the kinetic setpoint and D the fractal dimension. Single-channel currents were recorded at 146 mM external K+ from an inwardly rectifying, 120 pS, K+ selective, voltage-sensitive channel in cultured mouse hippocampal neurons. The kinetics of these channels were found to be statistically self-similar at different time scales as predicted by the fractal model. The fractal dimensions were approximately 2 for the closed times and approximately 1 for the open times and did not depend on voltage. For both the open and closed times the logarithm of the kinetic setpoint was found to be proportional to the applied voltage, which indicates that the gating of this channel involves the net inward movement of approximately one negative charge when this channel opens. Thus, the open and closed times and the voltage dependence of the gating of this channel are well described by the fractal model.

Time-domain reflectrometry studies of water binding and structural flexibility in chymotrypsin

Time-domain dielectric spectroscopy has been employed to probe the hydration properties and structural flexibility of chymotrypsin (EC 3.4.21.1). The dielectric properties of the hydrated protein above 100 MHz have been used to identify two categories of protein-bound water, the first being irrotationally bound to the protein with a second, relatively weakly bound, having a rotational freedom comparable with that of normal bulk water. A dielectric dispersion observed, centred at 12 MHz, has been attributed to the relaxation of the polar components of the protein structure. This dielectric loss became increasingly significant above a transition in the hydration dependence, where water is relatively weakly bound to the chymotrypsin. This is discussed in terms of the formation of water clusters on the protein surface which screen electrostatic interactions between protein-charged groups.

Monte Carlo-minimization approach to the multiple-minima problem in protein folding

A Monte Carlo-minimization method has been developed to overcome the multiple-minima problem. The Metropolis Monte Carlo sampling, assisted by energy minimization, surmounts intervening barriers in moving through successive discrete local minima in the multidimensional energy surface. The method has located the lowest-energy minimum thus far reported for the brain pentapeptide [Met5]enkephalin in the absence of water. Presumably it is the global minimum-energy structure. This supports the concept that protein folding may be a Markov process. In the presence of water, the molecules appear to exist as an ensemble of different conformations.

13C-NMR study of labeled vinyl groups in paramagnetic myoglobin derivatives

The 13C-NMR spectra of high-spin met-aquo myoglobin, spin-equilibrium met-azido myoglobin, low-spin met-cyano myoglobin, deoxy myoglobin and carbonmonoxy myoglobin from sperm whale reconstituted with hemin 13C enriched at both vinyl alpha or beta positions have been recorded. In all cases the labeled vinyl 13C signals are clearly resolved and useful spectra could be obtained within approx. 15 minutes. The decoupling of multiplet structure due to attached proton(s) has led to the specific assignment of vinyl 13C alpha signals in all paramagnetic derivatives and the 13C beta signals in met-cyano myoglobin. In all other cases, the collapse of the proton multiplet structure as a function of 1H decoupling frequency has located, but not assigned, the attached 1H resonance positions which are obscured by the intense diamagnetic envelope in the 1H-NMR spectrum. The resulting vinyl 13C hyperfine shifts follow Curie behavior, and the patterns closely resemble those in the appropriate model complexes in the same oxidation/spin/ligation state, except that the protein exhibits more in-plane asymmetry. The hyperfine shift patterns are indicative of dominant pi contact shifts for all ferric complexes. Deoxy myoglobin vinyl 13C and 1H contact shifts provide little evidence for pi bonding.

Protein structure probed by polarization spectroscopy. II. A time-resolved fluorescence study of human fibrinogen

Human fibrinogen in solution was studied by monitoring the time-resolved depolarization of the fluorescence emitted by two spectroscopic labels of which the fluorescence lifetimes differ by an order of magnitude. Contrary to a long-held view, no evidence of molecular flexibility was found in the 10-1000 ns range. In addition, from the rate of the overall rotation, it is proposed that a prolate and symmetric ellipsoid of 47 X 10.5 nm may represent the time-averaged hydrodynamic size and shape of the protein in solution. This rigid and highly hydrated structure (4 g water/g protein) accommodates the latest nodular models obtained from electron microscopy, explains the singular hydrodynamics of fibrinogen and, apparently, it would perform the two main functions of the protein in haemostasis, blood coagulation and platelet aggregation, more efficiently than the flexible molecule.

Fluorescence quenching dynamics of tryptophan in proteins. Effect of internal rotation under potential barrier

In many proteins fluorescence from single tryptophan exhibits a nonexponential decay function. To elucidate the origin of this nonexponential decay, we have examined the fluorescence decay function and time-resolved fluorescence anisotropy of a fluorophore covalently bound to a macromolecule by solving a rotational analogue of the Smoluchowski equation. An angular-dependent quenching constant and potential energy for the fluorophore undergoing internal rotation were introduced into the equation of motion for fluorophore. Results of numerical calculations using the equations thus obtained predict that both the fluorescence decay function and time-resolved anisotropy are dependent on rotational diffusion coefficients of fluorophore and potential energy for the internal rotation. The method was applied to the observed fluorescence decay curve of the single tryptophan in apocytochrome c from horse heart. The calculated decay curves fit the observed ones well.

Protein structure by solid-state NMR spectroscopy

The three-dimensional structures of proteins are among the most valuable contributions of biophysics to the understanding of biological systems (Dickerson & Geis, 1969; Creighton, 1983). Protein structures are utilized in the description and interpretation of a wide variety of biological phenomena, including genetic regulation, enzyme mechanisms, antibody recognition, cellular energetics, and macroscopic mechanical and structural properties of molecular assemblies. Virtually all of the information currently available about the structures of proteins at atomic resolution has been obtained from diffraction studies of single crystals of proteins (Wyckoffet al, 1985). However, recently developed NMR methods are capable of determining the structures of proteins and are now being applied to a variety of systems, including proteins in solution and other non-crystalline environments that are not amenable for X-ray diffraction studies. Solid-state NMR methods are useful for proteins that undergo limited overall reorientation by virtue of their being in the crystalline solid state or integral parts of supramolecular structures that do not reorient rapidly in solution. For reviews of applications of solid-state NMR spectroscopy to biological systems see Torchia and VanderHart (1979), Griffin (1981), Oldfieldet al. (1982), Opella (1982), Torchia (1982), Gauesh (1984), Torchia (1984) and Opella (1986). This review describes how solid-state NMR can be used to obtain structural information about proteins. Methods applicable to samples with macroscopic orientation are emphasized.

Anisotropy and anharmonicity of atomic fluctuations in proteins: analysis of a molecular dynamics simulation

Positional probability density functions (pdf) for the atomic fluctuations are determined from a molecular dynamics simulation for hen egg-white lysozyme. Most atoms are found to have motions that are highly anisotropic but only slightly anharmonic. The largest deviations from harmonic motion are in the direction of the largest rms fluctuations in the local principal axis frame. Backbone atoms tend to be more nearly harmonic than sidechain atoms. The atoms with the largest anharmonicities tend to have pdfs with multiple peaks, each of which is close to harmonic. Several model pdfs are evaluated on the basis of how well they fit probability densities from the dynamics simulations when parameterized in terms of the moments of the distribution. Gram-Charlier and Edgeworth perturbation expansions, which have been successful in describing the motions of small molecules in crystals, are shown to be inadequate for the distributions found in the dynamics of proteins. Multipeaked distribution functions are found to be more appropriate.