Hydrodynamic study of flexibility in immunoglobulin IgG1 using Brownian dynamics and the Monte Carlo simulations of a simple model

A simple bead model is proposed for the antibody molecule immunoglobulin IgG1. The partial flexibility of the hinge is represented by a quadratic potential associated to the angles between arms. Conformational and hydrodynamic properties are calculated using Monte Carlo (rigid-body) and Brownian dynamics simulations. Comparison of experimental and calculated values for some overall properties allows the assignment of dimensions and other model parameters. The Brownian dynamics technique is used next to simulate a rotational correlation function that is comparable with the decay of fluorescence emission anisotropy. This is done with varying flexibility at the hinge. The longest relaxation time shows a threefold decrease when going from the rigid Y-shaped conformation to the completely flexible case. The calculations are in good agreement with the decay times observed for IgG1. A flexibility analysis of the latter indicates that a variability of +/- 55 degrees (standard deviation) in the angle between the Fab arms.     

Backbone dynamics of proteins derived from carbonyl carbon relaxation times at 500, 600 and 800 MHz: Application to ribonuclease T1

The backbone dynamics of uniformly 13C/15N-enriched ribonuclease T1 have beeninvestigated using carbonyl carbon relaxation times recorded at three different spectrometerfrequencies. Pulse sequences for the determination of the longitudinal (T1) and transverse (T2)relaxation times are presented. The relaxation behaviour was analysed in terms of a multispinsystem. Although the chemical shift anisotropy relaxation mechanism dominates at highmagnetic field strength, the contributions of the dipole–dipole interactions and thecross-correlation between these two relaxation mechanisms have also been considered.Information about internal motions has been extracted from the relaxation data using themodel-free approach of Lipari and Szabo in order to determine order parameters (S2) andeffective internal correlation times (i). Using a relatively simple relation between themeasured relaxation rates and the spectral density function, an analytical expression for themicrodynamical parameters in dependence of T1 and T2 has been derived. The spectraldensity mapping technique has been applied in order to study the behaviour of the carbonylcarbon resonances in more detail.     

1H- and 2H-NMR relaxation in serum albumin solutions

The proton and deuteron relaxation times T1 and T2 were investigated in water and heavy water solutions of fatless human serum albumin. The temperature, concentration and Larmor precession frequency dependences can be well described by the conception of fast exchange in a simple biphasic model of water molecules rotation in the first hydration layer with slight anisotropy of motion. In the protonated systems the intermolecular dipole-dipole relaxation mechanism must be taken into account.     

Sampling of protein dynamics in nanosecond time scale by 15N NMR relaxation and self-diffusion measurements.

This paper presents a procedure for detection of intermediate nanosecond internal dynamics in globular proteins. The procedure uses 1H-15N relaxation measurements at several spectrometer frequencies and hydrodynamic calculations based on experimental self-diffusion coefficients. New heteronuclear experiments, using pulse field gradients, are introduced for the measurement of translation diffusion coefficients of 15N labeled proteins. An advanced interpretation of recently published (Luginbühl et al., Biochemistry, 36, 7305-7312 (1997)) backbone amide 15N relaxation data, measured at two spectrometers (400 and 750 MHz for 1H) for N-terminal DNA-binding domain (1-63) of 434 repressor, is presented. Non-applicability of commonly used fast (picosecond) dynamics model (FD) was justified by (i) poor fit of relaxation data by the FD model-free spectral density function both for isotropic and anisotropic models of the overall molecular tumbling; (ii) specific dependence of the overall rotation correlation times calculated from T1/T2 ratio on the spectrometer frequency; (iii) mismatch of the ratio of longitudinal 15N relaxation times T1, measured at different spectrometer frequencies, in comparison with that anticipated for the FD model; (iv) significantly underestimated overall rotation correlation time provided by the FD model (5.50+/-0.15 and 5.80+/-0.15 ns for 750 and 400 MHz spectrometer frequency respectively) in comparison with correlation time obtained from hydrodynamics. On the other hand, all relaxation and hydrodynamics data are in good correspondence with the model of intermediate (nanoseconds) dynamics. Overall rotation correlation time of 7.5+/-0.7 ns was calculated from experimental translation self-diffusion rate using hydrodynamics formalism (Garcia de la Torre, J. and Bloomfield, V.A. Quart. Rev. Biophys., 14, 81-139 (1981)). The statistical analysis of 15N relaxation data along with the hydrodynamic consideration clearly revealed that most of the residues in 434(1-63) repressor are involved in the nanosecond internal dynamics characterized by the the mean order parameters of 0.59+/-0.06 and the correlation times of ca. 5 ns.     

Rotation of lipids in membranes: molecular dynamics simulation, 31P spin-lattice relaxation, and rigid-body dynamics

Molecular dynamics simulations and ³¹P-NMR spin-lattice (R₁) relaxation rates from 0.022 to 21.1 T of fluid phase dipalmitoylphosphatidylcholine bilayers are compared. Agreement between experiment and direct prediction from simulation indicates that the dominant slow relaxation (correlation) times of the dipolar and chemical shift anisotropy spin-lattice relaxation are ∼10 ns and 3 ns, respectively. Overall reorientation of the lipid body, consisting of the phosphorus, glycerol, and acyl chains, is well described within a rigid-body model. Wobble, with D_⊥ = 1-2 × 10⁸ s⁻¹, is the primary component of the 10 ns relaxation; this timescale is consistent with the tumbling of a lipid-sized cylinder in a medium with the viscosity of liquid hexadecane. The value for D_|| the diffusion constant for rotation about the long axis of the lipid body, is difficult to determine precisely because of averaging by fast motions and wobble; it is tentatively estimated to be 1 × 10⁷ s⁻¹. The resulting D_||/D_⊥ ≈ 0.1 implies that axial rotation is strongly modulated by interactions at the lipid/water interface. Rigid-body modeling and potential of mean force evaluations show that the choline group is relatively uncoupled from the rest of the lipid. This is consistent with the ratio of chemical shift anisotropy and dipolar correlation times reported here and the previous observations that ³¹P-NMR lineshapes are axially symmetric even in the gel phase of dipalmitoylphosphatidylcholine.     

Time-resolved fluorescence of the single tryptophan residue in rat alpha-fetoprotein and rat serum albumin: analysis by the maximum-entropy method.

The time-resolved fluorescence emissions of the lone tryptophan residues in rat alpha-fetoprotein (RFP) and rat serum albumin (RSA) were studied. The total fluorescence intensity decays in both proteins were multiexponential. Analysis of the data by nonlinear least squares as a sum of discrete exponentials showed that four exponentials were needed for a satisfactory fit for both proteins. Analysis by the maximum entropy method using 150 logarithmically equally spaced exponentials yielded four well-resolved excited-state lifetime classes with barycenters and relative amplitudes values (ci) that corresponded to those obtained from the nonlinear least-squares method. Changing the temperature affected the relative amplitudes of the lifetime classes but had little effect on the lifetime values themselves. This suggests that the four classes reflect local conformational substates that exchange slowly with respect to the time window of observation defined by the longest lifetime. The internal rotational dynamics of the tryptophan in each protein was monitored by fluorescence anisotropy decay measurements. The mobility of the tryptophan appeared to be larger and faster in RFP than in RSA. The nonlinear least-squares analysis suggests the existence of three rotational correlation times of 0.1, 3, and 55 ns for this protein. As a function of temperature, the long correlation time did not follow the Perrin's law expected for a rigid rotating body. This suggests that this correlation time may reflect not only the Brownian rotation of the whole protein but also the flexibilities of domains in the protein. For RSA a two-component model with correlation times of 0.4 and 31 ns was sufficient to describe the data.(ABSTRACT TRUNCATED AT 250 WORDS)     

Polarization of the intrinsic fluorescence of proteins. I. Factors determining the form of the polarization spectrum

The dependence of the slopes of normalized Perrin plots on the excitation wavelength was established. It was shown that the cause of this effect is the anisotropy of the Brownian rotation of proteins, which must be regarded as asymmetrical particles with specific, but not random, orientation of tryptophan with respect to the main axes of the macromolecule. These findings were analysed on the basis of the rotational depolarization theory of such systems, applied for the case when bands of two absorption oscillators overlap as it is for oscillators 1La and 1Lb in the longest wavelength absorption band of tryptophan. It was shown that anisotropy of Brownian molecular rotation is one of the factors that determines the form of the polarization spectrum. The difference of the polarization spectrum of proteins from that of tryptophan, extrapolated to the infinite viscosity, is determined by energy transfer processes in proteins.     

Rotational dynamics of the Ca-ATPase in sarcoplasmic reticulum studied by time-resolved phosphorescence anisotropy

We have investigated the microsecond rotational motions of the Ca-ATPase in rabbit skeletal sarcoplasmic reticulum (SR), by measuring the time-resolved phosphorescence anisotropy of erythrosin 5-isothiocyanate (ERITC) covalently and specifically attached to the enzyme. Over a wide range of solvent conditions and temperatures, the phosphorescence anisotropy decay was best fit by a sum of three exponentials plus a constant term. At 4 degrees C, the rotational correlation times were phi 1 = 13 +/- 3 microseconds, phi 2 = 77 +/- 11 microseconds, and phi 3 = 314 +/- 23 microseconds. Increasing the solution viscosity with glycerol caused very little effect on the correlation times, while decreasing the lipid viscosity with diethyl ether decreased the correlation times substantially, indicating that the decay corresponds to rotation of the protein within the membrane, not to vesicle tumbling. The normalized residual anisotropy (A infinity) is insensitive to viscosity and temperature changes, supporting the model of uniaxial rotation of the protein about the membrane normal. The value of A infinity (0.20 +/- .02) indicates that each of the three decay components can be analyzed as a separate rotational species, with the preexponential factor Ai equal to 1.25X the mole fraction. An empirically accurate measurement of the membrane lipid viscosity was obtained, permitting a theoretical analysis of the correlation times in terms of the sizes of the rotating species. At 4 degrees C, the dominant correlation time (phi 3) is too large for a Ca-ATPase monomer, strongly suggesting that the enzyme is primarily aggregated (oligomeric).(ABSTRACT TRUNCATED AT 250 WORDS)     

Quenching-resolved emission anisotropy studies with single and multitryptophan-containing proteins

Measurements of the anisotropy of protein fluorescence as a function of an added collisional quencher, such as acrylamide, are used to construct Perrin plots. For single tryptophan containing proteins, such plots yield an apparent rotational correlation time for the depolarization process, which, in most cases, is approximately the value expected for Brownian rotation of the entire protein. Apparent limiting fluorescence anisotropy values, which range from 0.20 to 0.32 for the proteins studied, are also obtained from the Perrin plots. The lower values for the limiting anisotropy found for some proteins are interpreted as indicating the existence of relatively rapid, limited (within a cone of angle 0 degrees--30 degrees) motion of the tryptophan side chains that is independent of the overall rotation of the protein. Examples of the use of this fluorescence technique to study protein conformational changes are presented, including the monomer in equilibrium dimer equilibrium of beta-lactoglobulin, the monomer in equilibrium tetramer equilibrium of melittin, the N in equilibrium F transition of human serum albumin, and the induced change in the conformation of cod parvalbumin caused by the removal of Ca+2. Because multitryptophan-containing proteins have certain tryptophans that are accessible to solute quencher and others that are inaccessible, this method can be used to determine the steady state anisotropy of each class of tryptophan residues.     

(13)C-(1)H NMR relaxation and fluorescence anisotropy decay study of tyrosine dynamics in motilin

Tyrosine ring dynamics of the gastrointestinal hormone motilin was studied using two independent physical methods: fluorescence polarization anisotropy decay and NMR relaxation. Motilin, a 22-residue peptide, was selectively (13)C labeled in the ring epsilon-carbons of the single tyrosine residue. To eliminate effects of differences in peptide concentration, the same motilin sample was used in both experiments. NMR relaxation rates of the tyrosine ring C(epsilon)-H(epsilon) vectors, measured at four magnetic field strengths (9.4, 11.7, 14.1, and 18.8 Tesla) were used to map the spectral density function. When the data were analyzed using dynamic models with the same number of components, the dynamic parameters from NMR and fluorescence are in excellent agreement. However, the estimated rotational correlation times depend on the choice of dynamic model. The correlation times estimated from the two-component model-free approach and the three-component models were significantly different (1.7 ns and 2.2 ns, respectively). Various earlier studies of protein dynamics by NMR and fluorescence were compared. The rotational correlation times estimated by NMR for samples with high protein concentration were on average 18% longer for folded monomeric proteins than the corresponding times estimated by fluorescence polarization anisotropy decay, after correction for differences in viscosity due to temperature and D(2)O/H(2)O ratio.     

Speed of intersubunit communication in proteins.

To determine the speed of communication between protein subunits, time-resolved absorption spectra were measured following partial photodissociation of the carbon monoxide complex of hemoglobin. The experiments were carried out using linearly polarized, 10-ns laser pulses, with the polarization of the excitation pulse both parallel and perpendicular to the polarization of the probe pulse. The substantial contribution to the observed spectra from photoselection effects was eliminated by isotropically averaging the polarized spectra, allowing a detailed comparison of the kinetics as a function of the degree of photolysis. These results show that prior to 1 microsecond both geminate ligand rebinding and conformational relaxation are independent of the number of ligands dissociated from the hemoglobin tetramer, as expected for a two-state allosteric model. After this time the kinetics depend on the ligation state of the tetramer. The conformational relaxation at 10 microseconds can be interpreted in terms of the two-state allosteric model as arising from the R to T quaternary conformational change of both unliganded and singly liganded molecules. These results suggest that communication between subunits requires about 1 microsecond and that the mechanism of the communication which occurs after this time is via the R to T conformational change. The optical anisotropy provides a novel means of accurately determining the extinction coefficients of the transient photoproduct. The decay in the optical anisotropy, moreover, provides an accurate determination of the rotational correlation time of 36 +/- 3 ns.     

Interactions of pyridoxal kinase and aspartate aminotransferase emission anisotropy and compartmentation studies

Physical interactions between pyridoxal kinase and aspartate aminotransferase were detected by means of emission anisotropy and affinity chromatography techniques. Binding of aspartate aminotransferase (apoenzymes) to pyridoxal kinase tagged with a fluorescent probe was detected by emission anisotropy measurements at pH 6.8 (150 mM KCl). Upon saturation of the kinase with the aminotransferase, the emission anisotropy increases 22%. The protein complex is characterized by a dissociation constant of 3 microM. Time-dependent emission anisotropy measurements conducted with the mixture 5-naphthylamine-1-sulfonic acid-kinase aspartate aminotransferase (apoenzyme), revealed the presence of two rotational correlation times of phi 1 = 36 and phi 2 = 62 ns. The longer correlation time is attributed to the stable protein complex. By immobilizing one enzyme (pyridoxal kinase) through interactions with pyridoxal-Sepharose, it was possible to demonstrate that aspartate aminotransferase releases pyridoxal kinase. A test of compartmentation of pyridoxal-5-phosphate within the protein complex using alkaline phosphatase as trapping agent, indicates that the cofactor generated by the catalytic action of the kinase is channeled to the apotransaminase. The main function of the stable complex formed by the kinase and the aminotransferase is to hinder the release of free pyridoxal-5-phosphate into the bulk solvent.     

Time-resolved rotational dynamics of phosphorescent-labeled myosin heads in contracting muscle fibers.

We have measured the microsecond rotational motions of myosin heads in contracting rabbit psoas muscle fibers by detecting the transient phosphorescence anisotropy of eosin-5-maleimide attached specifically to the myosin head. Experiments were performed on small bundles (10-20 fibers) of glycerinated rabbit psoas muscle fibers at 4 degrees C. The isometric tension and physiological ATPase activity of activated fibers were unaffected by labeling 60-80% of the heads. Following excitation of the probes by a 10-ns laser pulse polarized parallel to the fiber axis, the time-resolved emission anisotropy of muscle fibers in rigor (no ATP) showed no decay from 1 microsecond to 1 ms (r infinity = 0.095), indicating that all heads are rigidly attached to actin on this time scale. In relaxation (5 mM MgATP but no Ca2+), the anisotropy decayed substantially over the microsecond time range, from an initial anisotropy (r0) of 0.066 to a final anisotropy (r infinity) of 0.034, indicating large-amplitude rotational motions with correlation times of about 10 and 150 microseconds and an overall angular range of 40-50 degrees. In isometric contraction (MgATP plus saturating Ca2+), the amplitude of the anisotropy decay (and thus the amplitude of the microsecond motion) is slightly less than in relaxation, and the rotational correlation times are about twice as long, indicating slower motions than those observed in relaxation. While the residual anisotropy (at 1 ms) in contraction is much closer to that in relaxation than in rigor, the initial anisotropy (at 1 microsecond) is approximately equidistant between those of rigor and relaxation.(ABSTRACT TRUNCATED AT 250 WORDS)     

Calcium-dependent changes in the flexibility of the regulatory domain of troponin C in the troponin complex

With the recent advances in structure determination of the troponin complex, it becomes even more important to understand the dynamics of its components and how they are affected by the presence or absence of Ca(2+). We used NMR techniques to study the backbone dynamics of skeletal troponin C (TnC) in the complex. Transverse relaxation-optimized spectroscopy pulse sequences and deuteration of TnC were essential to assign most of the TnC residues in the complex. Backbone amide (15)N relaxation times were measured in the presence of Ca(2+) or EGTA/Mg(2+). T(1) relaxation times could not be interpreted precisely, because for a molecule of this size, the longitudinal backbone amide (15)N relaxation rate due to chemical shift anisotropy and dipole-dipole interactions becomes too small, and other relaxation mechanisms become relevant. T(2) relaxation times were of the expected magnitude for a complex of this size, and most of the variation of T(2) times in the presence of Ca(2+) could be explained by the anisotropy of the complex, suggesting a relatively rigid molecule. The only exception was EF-hand site III and helix F immediately after, which are more flexible than the rest of the molecule. In the presence of EGTA/Mg(2+), relaxation times for residues in the C-domain of TnC are very similar to values in the presence of Ca(2+), whereas the N-domain becomes more flexible. Taken together with the high flexibility of the linker between the two domains, we concluded that in the absence of Ca(2+), the N-domain of TnC moves independently from the rest of the complex.     

Saturation transfer electron spin resonance study on the rotational diffusion of calcium- and magnesium-dependent adenosine triphosphatase in sarcoplasmic reticulum membranes.

The technique of saturation transfer electron spin resonance has been applied to study the rotational diffusion of spin-labeled Ca2+, Mg2+-dependent ATPase molecules in the membranes of sarcoplasmic reticulum vesicles. Comparison of the present data with those for spin-labeled hemoglobin undergoing isotropic rotation leads to a value of 2 X 10(-4) s for the apparent rotational correlation time at 20 degrees C for the membrane-bound protein. Consideration of the anisotropy of the Brownian rotation of the membrane-bound ATPase suggests that the true correlation time for the expected axial rotation may be somewhat smaller than the apparent value. An Arrhenius plot of the rotational motion shows a break, which is interpreted as indicating the occurrence of a conformational change of the ATPase molecule at about 15 degrees C.     

Proton magnetic relaxation dispersion in solutions of the cuproprotein diamine oxidase.

Proton nuclear magnetic resonance relaxation measurements were made over the range 4.7--220 MHz for aqueous solutions of hog kidney diamine oxidase. The values of 1/T1 give rise in two distinct dispersions, at 16 and 75 MHz, whereas 1/T2 displays a minimum at 20 MHz. The temperature dependence of relaxation rates in all cases yield apparent activation energies less than 0.6 kcal/mol. These data indicate to us that the two Cu(II) ions of diamine oxidase are intrinsically different in terms of their electronic relaxation characteristics and hence, chemical environments. Low field limits of the two electronic relaxation times are 2 and 10 ns, with one of these correlation times being frequency dependent. The value of the frequency-dependent electronic relaxation time is governed by interactions that are modulated by a process having a correlation time of 5 ps.     

Reversible unfolding of cytochrome c upon interaction with cardiolipin bilayers. 2. Evidence from phosphorus-31 NMR measurements.

31P NMR measurements were conducted to determine the structural and chemical environment of beef heart cardiolipin when bound to cytochrome c. 31P NMR line shapes infer that the majority of lipid remains in the bilayer state and that the average conformation of the lipid phosphate is not greatly affected by binding to the protein. An analysis of the spin-lattice (T1) relaxation times of hydrated cardiolipin as a function of temperature describes a T1 minimum at around 25 degrees C which leads to a correlation time for the phosphates in the lipid headgroup of 0.71 ns. The relaxation behavior of the protein-lipid complex was markedly different, showing a pronounced enhancement in the phosphorus spin-lattice relaxation rate. This effect of the protein increased progressively with increasing temperature, giving no indication of a minimum in T1 up to 75 degrees C. The enhancement in lipid phosphorus T1 relaxation was observed with protein in both oxidation states, being somewhat less marked for the reduced form. The characteristics of the T1 effects and the influence of the protein on other relaxation processes determined for the lipid phosphorus (spin-spin relaxation and longitudinal relaxation in the rotating frame) point to a strong paramagnetic interaction from the protein. A comparison with the relaxation behavior of samples spinning at the "magic angle" was also consistent with this mechanism. The results suggest that cytochrome c reversibly denatures on complexation with cardiolipin bilayers, such that the electronic ground state prevailing in the native structure of both oxidized and reduced protein can convert to high-spin states with greater magnetic susceptibility.