Ligand-Induced Protein Responses and Mechanical Signal Propagation Described by Linear Response Theories

In this study, a general linear response theory (LRT) is formulated to describe time-dependent and -independent protein conformational changes upon CO binding with myoglobin. Using the theory, we are able to monitor protein relaxation in two stages. The slower relaxation is found to occur from 4.4 to 81.2 picoseconds and the time constants characterized for a couple of aromatic residues agree with those observed by UV Resonance Raman (UVRR) spectrometry and time resolved x-ray crystallography. The faster “early responses”, triggered as early as 400 femtoseconds, can be best described by the theory when impulse forces are used. The newly formulated theory describes the mechanical propagation following ligand-binding as a function of time, space and types of the perturbation forces. The “disseminators”, defined as the residues that propagate signals throughout the molecule the fastest among all the residues in protein when perturbed, are found evolutionarily conserved and the mutations of which have been shown to largely change the CO rebinding kinetics in myoglobin.     

A study of vibrational relaxation of B-state carbon monoxide in the heme pocket of photolyzed carboxymyoglobin.

The vibrational energy relaxation of dissociated carbon monoxide in the heme pocket of sperm whale myoglobin has been studied using equilibrium molecular dynamics simulation and normal mode analysis methods. Molecular dynamics trajectories of solvated myoglobin were run at 300 K for both the delta- and epsilon-tautomers of the distal histidine, His64. Vibrational population relaxation times were estimated using the Landau-Teller model. For carbon monoxide (CO) in the myoglobin epsilon-tautomer, for a frequency of omega0 = 2131 cm-1 corresponding to the B1 state, T1epsilon(B1) = 640 +/- 185 ps, and for a frequency of omega0 = 2119 cm-1 corresponding to the B2 state, T1epsilon(B2) = 590 +/- 175 ps. Although the CO relaxation rates in both the epsilon- and delta-tautomers are similar in magnitude, the simulations predict that the vibrational relaxation of the CO is faster in the delta-tautomer. For CO in the myoglobin delta-tautomer, it was found that the relaxation times were identical within error for the two CO substate frequencies, T1delta(B1) = 335 +/- 115 ps and T1delta(B2) = 330 +/- 145 ps. These simulation results are in reasonable agreement with experimental results of Anfinrud and coworkers (unpublished results). Normal mode calculations were used to identify the dominant coupling between the protein and CO molecules. The calculations suggest that the residues of the myoglobin pocket, acting as a first solvation shell to the CO molecule, contribute the primary "doorway" modes in the vibrational relaxation of the oscillator.     

Myoglobin and cytochrome b5: a nuclear magnetic resonance study of a highly dynamic protein complex

The transient complex of bovine myoglobin and cytochrome b(5) has been investigated using a combination of NMR chemical shift mapping, (15)N relaxation data, and protein docking simulations. Chemical shift perturbations observed for cytochrome b(5) amide resonances upon complex formation with either metmyoglobin (Fe(III)) or carbon monoxide-bound myoglobin (Fe(II)) are more than 10-fold smaller than in other transient redox protein complexes. From (15)N relaxation experiments, an increase in the overall correlation time of cytochrome b(5) in the presence of myoglobin is observed, confirming that complex formation is occurring. The chemical shift perturbations of proton and nitrogen amide nuclei as well as heme protons of cytochrome b(5) titrate with increasing myoglobin concentrations, also demonstrating the formation of a weak complex with a K(a) in the inverse millimolar range. The perturbed residues map over a wide surface area of cytochrome b(5), with patches of residues located around the exposed heme 6-propionate as well as at the back of the protein. The nature of the affected residues is mostly negatively charged contrary to perturbed residues in other transient complexes, which are mainly hydrophobic or polar. Protein docking simulations using the NMR data as constraints show several docking geometries both close to and far away from the exposed heme propionates of myoglobin. Overall, the data support the emerging view that this complex consists of a dynamic ensemble of orientations in which each protein constantly diffuses over the surface of the other. The characteristic NMR features may serve as a structural tool for the identification of such dynamic complexes.     

Influence of histidine tag attachment on picosecond protein dynamics

Polyhistidine affinity tags are routinely employed as a convenient means of purifying recombinantly expressed proteins. A tacit assumption is commonly made that His tags have little influence on protein structure and function. Attachment of a His tag to the N-terminus of the robust globular protein myoglobin leads to only minor changes to the electrostatic environment of the heme pocket, as evinced by the nearly unchanged Fourier transform infrared spectrum of CO bound to the heme of His-tagged myoglobin. Experiments employing two-dimensional infrared vibrational echo spectroscopy of the heme-bound CO, however, find that significant changes occur to the short time scale (picoseconds) dynamics of myoglobin as a result of His tag incorporation. The His tag mainly reduces the dynamics on the 1.4 ps time scale and also alters protein motions of myoglobin on the slower, >100 ps time scale, as demonstrated by the His tag's influence on the fluctuations of the CO vibrational frequency, which reports on protein structural dynamics. The results suggest that affinity tags may have effects on protein function and indicate that investigators of affinity-tagged proteins should take this into consideration when investigating the dynamics and other properties of such proteins.     

The transition from inhomogeneous to homogeneous kinetics in CO binding to myoglobin.

Heme proteins react inhomogeneously with ligands at cryogenic temperatures and homogeneously at room temperature. We have identified and characterized a transition from inhomogeneous to homogeneous behavior at intermediate temperatures in the time dependence of CO binding to horse myoglobin. The turnover is attributed to a functionally important tertiary protein relaxation process during which the barrier increases dynamically. This is verified by a combination of theory and multipulse measurements. A likely biological significance of this effect is in the autocatalysis of the ligand release process.     

From metmyoglobin to deoxy myoglobin: relaxations of an intermediate state

Metmyoglobin has been reduced at low temperature (below 100 K) using x-rays or by excitation of tris(2,2′,bipyridine)ruthenium(II) chloride with visible light. Upon reduction, an intermediate state is formed where the structure of the protein is very similar to that of metmyoglobin with the water molecule still bound to the heme iron, but the iron is II low spin. The nature of the intermediate state has been investigated with optical spectroscopy. The Q_o and Q_v bands of the intermediate state are split, suggesting that the protoporphyrin is distorted. The intermediate state undergoes a relaxation observed by a shifting of the Soret band at temperatures above 80 K. Above 140 K, the protein begins to relax to the deoxy conformation. The relaxation kinetics of the protein have been monitored optically as a function of time and temperature from minutes to several hours and from 150 K to 190 K. By measuring the entire visible spectrum, we are able to distinguish between electron transfer processes and the protein relaxation from the intermediate state to deoxy myoglobin. The relaxation has been measured in both horse myoglobin and sperm whale myoglobin with the relaxation occurring on faster time scales in horse myoglobin. Both the reduction kinetics and the relaxation show non-exponential behavior. The reduction kinetics can be fit well to a stretched exponential. The structural relaxation from the intermediate state to the deoxy conformation shows a more complex, dynamical behavior and the reaction is most likely affected by the relaxation of the protein within the intermediate state. Received: 30 June 1997 / Accepted: 6 November 1997     

Selective lateral muscle activation in moderate medial knee osteoarthritis subjects does not unload medial knee condyle

There is some debate in the literature regarding the role of quadriceps-hamstrings co-contraction in the onset and progression of knee osteoarthritis. Does co-contraction during walking increase knee contact loads, thereby causing knee osteoarthritis, or might it be a compensatory mechanism to unload the medial tibial condyle? We used a detailed musculoskeletal model of the lower limb to test the hypothesis that selective activation of lateral hamstrings and quadriceps, in conjunction with inhibited medial gastrocnemius, can actually reduce the joint contact force on the medial compartment of the knee, independent of changes in kinematics or external forces. “Baseline” joint loads were computed for eight subjects with moderate medial knee osteoarthritis (OA) during level walking, using static optimization to resolve the system of muscle forces for each subject?s scaled model. Holding all external loads and kinematics constant, each subject?s model was then perturbed to represent non-optimal “OA-type” activation based on mean differences detected between electromyograms (EMG) of control and osteoarthritis subjects. Knee joint contact forces were greater for the “OA-type” than the “Baseline” distribution of muscle forces, particularly during early stance. The early-stance increase in medial contact load due to the “OA-type” perturbation could implicate this selective activation strategy as a cause of knee osteoarthritis. However, the largest increase in the contact load was found at the lateral condyle, and the “OA-type” lateral activation strategy did not increase the overall (greater of the first or second) medial peak contact load. While “OA-type” selective activation of lateral muscles does not appear to reduce the medial knee contact load, it could allow subjects to increase knee joint stiffness without any further increase to the peak medial contact load.     

Protein conformational relaxation following photodissociation of CO from carbonmonoxymyoglobin: picosecond circular dichroism and absorption studies

Picosecond time-resolved polarization spectroscopy is used to study relaxation dynamics in myoglobin following photoelimination of CO from carbonmonoxymyoglobin. Evolution of the transient circular dichroism signal of the N band of myoglobin (probed at 355 nm) to that characteristic of equilibrium myoglobin requires approximately 300 ps. This time scale is significantly longer than that corresponding to the photoinitiated bond cleavage. Transient linear dichroism of the Soret band and picosecond time-resolved magnetic circular dichroism measurements of the Q band demonstrate that the circular dichroism kinetics do not result from either time-dependent changes in the orientation of the transition moments of the heme ring or the doming of the heme that accompanies the out-of-plane motion of the iron. Finally, transient absorption data of the near-IR optical transition of photogenerated myoglobin suggest that the circular dichroism data are not a measure of the tilting of the proximal histidine. The circular dichroism data are discussed in terms of a relaxation in the tertiary structure of the protein following dissociation.     

Relaxation of ultralarge VWF bundles in a microfluidic-AFM hybrid reactor

The crucial role of the biopolymer “Von Willebrand factor” (VWF) in blood platelet binding is tightly regulated by the shear forces to which the protein is exposed in the blood flow. Under high-shear conditions, VWFs ability to immobilize blood platelets is strongly increased due to a change in conformation which at sufficient concentration is accompanied by the formation of ultra large VWF bundles (ULVWF). However, little is known about the dynamic and mechanical properties of such bundles.Combining a surface acoustic wave (SAW) based microfluidic reactor with an atomic force microscope (AFM) we were able to study the relaxation of stretched VWF bundles formed by hydrodynamic stress. We found that the dynamical response of the network is well characterized by stretched exponentials, indicating that the relaxation process proceeds through hopping events between a multitude of minima. This finding is in accordance with current ideas of VWF self-association. The longest relaxation time does not show a clear dependence on the length of the bundle, and is dominated by the internal conformations and effective friction within the bundle.     

Protein Collective Motions Coupled to Ligand Migration in Myoglobin

Ligand migration processes inside myoglobin and protein dynamics coupled to the migration were theoretically investigated with molecular dynamics simulations. Based on a linear response theory, we identified protein motions coupled to the transient migration of ligand, carbon monoxide (CO), through channels. The result indicates that the coupled protein motions involve collective motions extended over the entire protein correlated with local gating motions at the channels. Protein motions, coupled to opening of a channel from the distal pocket to a neighboring xenon site, were found to share the collective motion with experimentally observed protein motions coupled to a doming motion of the heme Fe atom upon photodissociation of the ligand. Analysis based on generalized Langevin dynamics elucidated slow and diffusive features of the protein response motions. Remarkably small transmission coefficients for rates of the CO migrations through myoglobin were found, suggesting that the CO migration dynamics are characterized as motions governed by the protein dynamics involving the collective motions, rather than as thermally activated transitions across energy barriers of well-structured channels.     

Ligand dynamics in an electron transfer protein. Picosecond geminate recombination of carbon monoxide to heme in mutant forms of cytochrome c

Substitution of the heme coordination residue Met-80 of the electron transport protein yeast iso-1-cytochrome c allows external ligands like CO to bind and thus increase the effective redox potential. This mutation, in principle, turns the protein into a quasi-native photoactivable electron donor. We have studied the kinetic and spectral characteristics of geminate recombination of heme and CO in a series of single M80X (X = Ala, Ser, Asp, Arg) mutants, using femtosecond transient absorption spectroscopy. In these proteins, all geminate recombination occurs on the picosecond and early nanosecond time scale, in a multiphasic manner, in which heme relaxation takes place on the same time scale. The extent of geminate recombination varies from >99% (Ala, Ser) to approximately 70% (Arg), the latter value being in principle low enough for electron injection studies. The rates and extent of the CO geminate recombination phases are much higher than in functional ligand-binding proteins like myoglobin, presumably reflecting the rigid and hydrophobic properties of the heme environment, which are optimized for electron transfer. Thus, the dynamics of CO recombination in cytochrome c are a tool for studying the heme pocket, in a similar way as NO in myoglobin. We discuss the differences in the CO kinetics between the mutants in terms of the properties of the heme environment and strategies to enhance the CO escape yield. Experiments on double mutants in which Phe-82 is replaced by Asp or Gly as well as the M80D substitution indicate that such steric changes substantially increase the motional freedom-dissociated CO.     

Key residue-dominated protein folding dynamics

A “key-residue” hypothesis that a few residues’ characteristics contain the essential dynamics of the whole protein is proposed for the study of side-chain relaxation near native states. Molecular dynamics simulation is performed on the folding of Trp-cage, and four key residues are discovered and shown to be highly sensitive to the change of state of the protein away from the native state. Order parameters that characterize the geometrical properties of key residues are shown to form valuable phase plane on which one distinguishes different reaction pathways. Furthermore, one of the key residues, Trp6, is observed to display two reconfiguration processes, in which one is induced by an unconstrained torsion of the side-chain of Trp6, with a rate faster by almost an order of magnitude than the other one described by Kussell’s model. The faster process seems to occur more frequently in our simulation and thus represent a significant mechanism in folding dynamics.     

Amino-terminal domain interactions of lambda integrase on arm-type DNA

In contrast to the other tyrosine recombinase family members, integrase protein (Int) of bacteriophage λ has an additional amino-terminal domain that binds to “arm-type” DNA sequences distant from those involved in strand exchange. The homomeric interaction between neighboring amino-terminal domains of Int is contributed by R30-D71 salt-bridge in a non-equivalent manner on Holliday-junction intermediates. In this report, R30 and D71 residues were investigated in regard to Int’s cooperative binding to “arm-type” DNA and the attenuating function of “arm-type” DNA. The results suggest the electrostatic interaction between residues 30 and 71 is dependent on “arm-type” DNA and contributes the “selective” inhibition of catalytic activity of λ Int by “arm-type” DNA.     

Configural processing of different topologically structured figures: an ERP Study

According to Chen’s theory, topological differences are perceived faster than feature differences in early visual perception. We hypothesized that topological perception is caused by the sensitivity in discriminating figures with and without “holes”. An ERP experiment was conducted utilizing a passive paradigm to investigate the differences in perceiving figures with “hole” and with “no-hole”. The results showed differences in N170 components between figures with “holes” and with “no-holes”. The inversion of the “hole” could influence the latency of N170, but the inversion of the “no-hole” could not, which indicated that global features are processed first in the “hole” perception whilst local features are given priority to the “no-hole” perception. This result was similar to studies concerning face and non-face objects, suggesting a configural processing of the “hole”.     

Time-resolved resonance Raman study on ultrafast structural relaxation and vibrational cooling of photodissociated carbonmonoxy myoglobin

A localized small structural change is converted to a higher order conformational change of protein and extends to a mesoscopic scale to induce a physiological function. To understand such features of protein, ultrafast dynamics of myoglobin (Mb) following photolysis of carbon monoxide were investigated. Recent results are summarized here with a stress on structural and vibrational energy relaxation. The core expansion of heme takes place within 2 ps but the out of plane displacement of the heme iron and the accompanying protein conformational change occur in 10 and 100 s of the picosecond regimes, respectively. Unexpectedly, it was found from UV resonance Raman spectra that Trp7 in the N-terminal region and Tyr151 in the C-terminal region undergo appreciable structural changes upon ligand binding-dissociation while Tyr104, Tyr146, and Trp14 do not. Because of the communication between the movements of these surface residues and the heme iron, the rate of spectral change of the iron-histidine (Fe- His) stretching band after CO photodissociation is influenced by the viscosity of solvent. Temporal changes of the anti-Stokes Raman intensity demonstrated immediate generation of vibrationally excited heme upon photodissociation and its decay with a time constant of 1-2 ps.     

The role of RGS protein in agonist-dependent relaxation of GIRK currents in Xenopus oocytes

G protein coupled inward rectifier K⁺ channels (GIRK) are activated by the G_βγ subunits of G proteins upon activation of G protein coupled receptors (GPCRs). Receptor-stimulated GIRK currents are known to possess a curious property, termed “agonist-dependent relaxation,” denoting a slow current increase upon stepping the membrane voltage from positive to negative potentials. Regulators of G protein signaling (RGS) proteins have earlier been implicated in this phenomenon since RGS coexpression was required for relaxation to be observed in heterologous expression systems. However, a recent study presented contrasting evidence that GIRK current relaxation reflects voltage sensitive agonist binding to the GPCR. The present study re-examined the role of RGS protein in agonist-dependent relaxation and found that RGS coexpression is not necessary for the relaxation phenomenon. However, RGS4 speeds up relaxation kinetics, allowing the phenomenon to be observed using shorter voltage steps. These findings resolve the controversy regarding the role of RGS protein vs. GPCR voltage sensitivity in mediating agonist-dependent relaxation of GIRK currents.     

Microwave-assisted specific chemical digestion for rapid protein identification

We have developed a rapid microwave-assisted protein digestion technique based on classic acid hydrolysis reaction with 2% formic acid solution. In this mild chemical environment, proteins are hydrolyzed to peptides, which can be directly analyzed by MALDI-MS or ESI-MS without prior sample purification. Dilute formic acid cleaves proteins specifically at the C-terminal of aspartyl (Asp) residues within 10 min of exposure to microwave irradiation. By adjusting the irradiation time, we found that the extent of protein fragmentation can be controlled, as shown by the single fragmentation of myoglobin at the C-terminal of any of the Asp residues. The efficacy and simplicity of this technique for protein identification are demonstrated by the peptide mass maps of in-gel digested myoglobin and BSA, as well as proteins isolated from Escherichia coli K12 cells.     

Ligand binding and protein relaxation in heme proteins: a room temperature analysis of NO geminate recombination

Ultrafast absorption spectroscopy is used to study heme-NO recombination at room temperature in aqueous buffer on time scales where the ligand cannot leave its cage environment. While a single barrier is observed for the cage recombination of NO with heme in the absence of globin, recombination in hemoglobin and myoglobin is nonexponential. Examination of hemoglobin with and without inositol hexaphosphate points to proximal constraints as important determinants of the geminate rebinding kinetics. Molecular dynamics simulations of myoglobin and heme-imidazole subsequent to ligand dissociation were used to investigate the transient behavior of the Fe-proximal histidine coordinate and its possible involvement in geminate recombination. The calculations, in the context of the absorption measurements, are used to formulate a distinction between nonexponential rebinding that results from multiple protein conformations (substates) present at equilibrium or from nonequilibrium relaxation of the protein triggered by a perturbation such as ligand dissociation. The importance of these two processes is expected to depend on the time scale of rebinding relative to equilibrium fluctuations and nonequilibrium relaxation. Since NO rebinding occurs on the picosecond time scale of the calculated myoglobin relaxation, a time-dependent barrier is likely to be an important factor in the observed nonexponential kinetics. The general implications of the present results for ligand binding in heme proteins and its time and temperature dependence are discussed. It appears likely that, at low temperatures, inhomogeneous protein populations play an important role and that as the temperature is raised, relaxation effects become significant as well.