The role of structure, energy landscape, dynamics, and allostery in the enzymatic function of myoglobin

The grail of protein science is the connection between structure and function. For myoglobin (Mb) this goal is close. Described as only a passive dioxygen storage protein in texts, we argue here that Mb is actually an allosteric enzyme that can catalyze reactions among small molecules. Studies of the structural, spectroscopic, and kinetic properties of Mb lead to a model that relates structure, energy landscape, dynamics, and function. Mb functions as a miniature chemical reactor, concentrating and orienting diatomic molecules such as NO, CO, O2, and H2O2 in highly conserved internal cavities. Reactions can be controlled because Mb exists in distinct taxonomic substates with different catalytic properties and connectivities of internal cavities.

The use of site-directed fluorophore labeling and donor–donor energy migration to investigate solution structure and dynamics in proteins

The use of molecular genetics for introducing fluorescent molecules enables the use of donor–donor energy migration to determine intramolecular distances in a variety of proteins. This approach can be applied to examine the overall molecular dimensions of proteins and to investigate structural changes upon interactions with specific target molecules. In this report, the donor–donor energy migration method is demonstrated by experiments with the latent form of plasminogen activator inhibitor type 1. Based on the known x-ray structure of plasminogen activator inhibitor type 1, three positions forming the corners of a triangle were chosen. Double Cys substitution mutants (V106C-H185C, H185C-M266C, and M266C-V106C) and corresponding single substitution mutants (V106C, H185C, and M266C) were created and labeled with a sulfhydryl specific derivative of BODIPY (=the D molecule). The side lengths of this triangle were obtained from analyses of the experimental data. The analyses account for the local anisotropic order and rotational motions of the D molecules, as well as for the influence of a partial DD-labeling. The distances, as determined from x-ray diffraction, between the Cα-atoms of the positions V106C–H185C, H185C–M266C, and M266C–V106C were 60.9, 30.8, and 55.1 Å, respectively. These are in good agreement with the distances of 54 ± 4, 38 ± 3, and 55 ± 3 Å, as determined between the BODIPY groups attached via linkers to the same residues. Although the positions of the D-molecules and the Cα-atoms physically cannot coincide, there is a reasonable agreement between the methods.

Molecular dynamics and free-energy calculations applied to affinity maturation in antibody 48G7

We investigated the relative free energies of hapten binding to the germ line and mature forms of the 48G7 antibody Fab fragments by applying a continuum model to structures sampled from molecular dynamics simulations in explicit solvent. Reasonable absolute and very good relative free energies were obtained. As a result of nine somatic mutations that do not contact the hapten, the affinity-matured antibody binds the hapten >104 tighter than the germ line antibody. Energetic analysis reveals that van der Waals interactions and nonpolar contributions to solvation are similar and drive the formations of both the germ line and mature antibody–hapten complexes. Affinity maturation of the 48G7 antibody therefore appears to occur through reorganization of the combining site geometry in a manner that optimizes the balance of gaining favorable electrostatic interactions with the hapten and losing those with solvent during the binding process. As reflected by lower rms fluctuations in the antibody–hapten complex, the mature complex undergoes more restricted fluctuations than the germ line complex. The dramatically increased affinity of the 48G7 antibody over its germ line precursor is thus made possible by electrostatic optimization.

Dynamics and folding of single two-stranded coiled-coil peptides studied by fluorescent energy transfer confocal microscopy

We report single-molecule measurements on the folding and unfolding conformational equilibrium distributions and dynamics of a disulfide crosslinked version of the two-stranded coiled coil from GCN4. The peptide has a fluorescent donor and acceptor at the N termini of its two chains and a Cys disulfide near its C terminus. Thus, folding brings the two N termini of the two chains close together, resulting in an enhancement of fluorescent resonant energy transfer. End-to-end distance distributions have thus been characterized under conditions where the peptide is nearly fully folded (0 M urea), unfolded (7.4 M urea), and in dynamic exchange between folded and unfolded states (3.0 M urea). The distributions have been compared for the peptide freely diffusing in solution and deposited onto aminopropyl silanized glass. As the urea concentration is increased, the mean end-to-end distance shifts to longer distances both in free solution and on the modified surface. The widths of these distributions indicate that the molecules are undergoing millisecond conformational fluctuations. Under all three conditions, these fluctuations gave nonexponential correlations on 1- to 100-ms time scale. A component of the correlation decay that was sensitive to the concentration of urea corresponded to that measured by bulk relaxation kinetics. The trajectories provided effective intramolecular diffusion coefficients as a function of the end-to-end distances for the folded and unfolded states. Single-molecule folding studies provide information concerning the distributions of conformational states in the folded, unfolded, and dynamically interconverting states.

Single-pair fluorescence resonance energy transfer on freely diffusing molecules: Observation of Förster distance dependence and subpopulations

Photon bursts from single diffusing donor-acceptor labeled macromolecules were used to measure intramolecular distances and identify subpopulations of freely diffusing macromolecules in a heterogeneous ensemble. By using DNA as a rigid spacer, a series of constructs with varying intramolecular donor-acceptor spacings were used to measure the mean and distribution width of fluorescence resonance energy transfer (FRET) efficiencies as a function of distance. The mean single-pair FRET efficiencies qualitatively follow the distance dependence predicted by Förster theory. Possible contributions to the widths of the FRET efficiency distributions are discussed, and potential applications in the study of biopolymer conformational dynamics are suggested. The ability to measure intramolecular (and intermolecular) distances for single molecules implies the ability to distinguish and monitor subpopulations of molecules in a mixture with different distances or conformational states. This is demonstrated by monitoring substrate and product subpopulations before and after a restriction endonuclease cleavage reaction. Distance measurements at single-molecule resolution also should facilitate the study of complex reactions such as biopolymer folding. To this end, the denaturation of a DNA hairpin was examined by using single-pair FRET.

Brownian dynamics simulations of protein folding: Access to milliseconds time scale and beyond

Protein folding occurs on a time scale ranging from milliseconds to minutes for a majority of proteins. Computer simulation of protein folding, from a random configuration to the native structure, is nontrivial owing to the large disparity between the simulation and folding time scales. As an effort to overcome this limitation, simple models with idealized protein subdomains, e.g., the diffusion–collision model of Karplus and Weaver, have gained some popularity. We present here new results for the folding of a four-helix bundle within the framework of the diffusion–collision model. Even with such simplifying assumptions, a direct application of standard Brownian dynamics methods would consume 10,000 processor-years on current supercomputers. We circumvent this difficulty by invoking a special Brownian dynamics simulation. The method features the calculation of the mean passage time of an event from the flux overpopulation method and the sampling of events that lead to productive collisions even if their probability is extremely small (because of large free-energy barriers that separate them from the higher probability events). Using these developments, we demonstrate that a coarse-grained model of the four-helix bundle can be simulated in several days on current supercomputers. Furthermore, such simulations yield folding times that are in the range of time scales observed in experiments.

Linking topography of its potential surface with the dynamics of folding of a protein model

The “3-color, 46-bead” model of a folding polypeptide is the vehicle for adapting to proteins a mode of analysis used heretofore for atomic clusters, to relate the topography of the potential surface to the dynamics that lead to formation of selected structures. The analysis is based on sequences of stationary points—successive minima, joined by saddles—that rise monotonically in energy from basin bottoms. Like structure-seeking clusters, the potential surface of the model studied here is staircase-like, rather than sawtooth-like, with highly collective motions required for passage from one minimum to the next. The surface has several deep basins whose minima correspond to very similar structures, but which are separated by high energy barriers.

Fluorescence and photobleaching dynamics of single light-harvesting complexes

Single light-harvesting complexes LH-2 from Rhodopseudomonas acidophila were immobilized on various charged surfaces under physiological conditions. Polarized light experiments showed that the complexes were situated on the surface as nearly upright cylinders. Their fluorescence lifetimes and photobleaching properties were obtained by using a confocal fluorescence microscope with picosecond time resolution. Initially all molecules fluoresced with a lifetime of 1 ± 0.2 ns, similar to the bulk value. The photobleaching of one bacteriochlorophyll molecule from the 18-member assembly caused the fluorescence to switch off completely, because of trapping of the mobile excitations by energy transfer. This process was linear in light intensity. On continued irradiation the fluorescence often reappeared, but all molecules did not show the same behavior. Some LH-2 complexes displayed a variation of their quantum yields that was attributed to photoinduced confinement of the excited states and thereby a diminution of the superradiance. Others showed much shorter lifetimes caused by excitation energy traps that are only ≈3% efficient. On repeated excitation some molecules entered a noisy state where the fluorescence switched on and off with a correlation time of ≈0.1 s. About 490 molecules were examined.

Stability and dynamics in a hyperthermophilic protein with melting temperature close to 200°C

The rubredoxin protein from the hyperthermophilic archaebacterium Pyrococcus furiosus was examined by a hydrogen exchange method. Even though the protein does not exhibit reversible thermal unfolding, one can determine its stability parameters—free energy, enthalpy, entropy, and melting temperature—and also the distribution of stability throughout the protein, by using hydrogen exchange to measure the reversible cycling of the protein between native and unfolded states that occurs even under native conditions.

Vibrational population relaxation of carbon monoxide in the heme pocket of photolyzed carbonmonoxy myoglobin: Comparison of time-resolved mid-IR absorbance experiments and molecular dynamics simulations

The vibrational energy relaxation of carbon monoxide in the heme pocket of sperm whale myoglobin was studied by using molecular dynamics simulation and normal mode analysis methods. Molecular dynamics trajectories of solvated myoglobin were run at 300 K for both the δ- and ɛ-tautomers of the distal His-64. Vibrational population relaxation times of 335 ± 115 ps for the δ-tautomer and 640 ± 185 ps for the ɛ-tautomer were estimated by using the Landau–Teller model. Normal mode analysis was used to identify those protein residues that act as the primary “doorway” modes in the vibrational relaxation of the oscillator. Although the CO relaxation rates in both the ɛ- and δ-tautomers are similar in magnitude, the simulations predict that the vibrational relaxation of the CO is faster in the δ-tautomer with the distal His playing an important role in the energy relaxation mechanism. Time-resolved mid-IR absorbance measurements were performed on photolyzed carbonmonoxy hemoglobin (Hb13CO). From these measurements, a T1 time of 600 ± 150 ps was determined. The simulation and experimental estimates are compared and discussed.

Femtosecond dynamics of the forbidden carotenoid S1 state in light-harvesting complexes of purple bacteria observed after two-photon excitation

Time-resolved excited-state absorption intensities after direct two-photon excitation of the carotenoid S1 state are reported for light-harvesting complexes of purple bacteria. Direct excitation of the carotenoid S1 state enables the measurement of subsequent dynamics on a fs time scale without interference from higher excited states, such as the optically allowed S2 state or the recently discovered dark state situated between S1 and S2. The lifetimes of the carotenoid S1 states in the B800-B850 complex and B800-B820 complex of Rhodopseudomonas acidophila are 7 ± 0.5 ps and 6 ± 0.5 ps, respectively, and in the light-harvesting complex 2 of Rhodobacter sphaeroides ≈1.9 ± 0.5 ps. These results explain the differences in the carotenoid to bacteriochlorophyll energy transfer efficiency after S2 excitation. In Rps. acidophila the carotenoid S1 to bacteriochlorophyll energy transfer is found to be quite inefficient (φET1 <28%) whereas in Rb. sphaeroides this energy transfer is very efficient (φET1 ≈80%). The results are rationalized by calculations of the ensemble averaged time constants. We find that the Car S1 → B800 electronic energy transfer (EET) pathway (≈85%) dominates over Car S1 → B850 EET (≈15%) in Rb. sphaeroides, whereas in Rps. acidophila the Car S1 → B850 EET (≈60%) is more efficient than the Car S1 → B800 EET (≈40%). The individual electronic couplings for the Car S1 → BChl energy transfer are estimated to be approximately 5–26 cm−1. A major contribution to the difference between the energy transfer efficiencies can be explained by different Car S1 energy gaps in the two species.

The vibrational energy flow transition in organic molecules: Theory meets experiment

Most large dynamical systems are thought to have ergodic dynamics, whereas small systems may not have free interchange of energy between degrees of freedom. This assumption is made in many areas of chemistry and physics, ranging from nuclei to reacting molecules and on to quantum dots. We examine the transition to facile vibrational energy flow in a large set of organic molecules as molecular size is increased. Both analytical and computational results based on local random matrix models describe the transition to unrestricted vibrational energy flow in these molecules. In particular, the models connect the number of states participating in intramolecular energy flow to simple molecular properties such as the molecular size and the distribution of vibrational frequencies. The transition itself is governed by a local anharmonic coupling strength and a local state density. The theoretical results for the transition characteristics compare well with those implied by experimental measurements using IR fluorescence spectroscopy of dilution factors reported by Stewart and McDonald [Stewart, G. M. & McDonald, J. D. (1983) J. Chem. Phys. 78, 3907–3915].

Spatiotemporal dynamics of guanosine 3′,5′-cyclic monophosphate revealed by a genetically encoded, fluorescent indicator

To investigate the dynamics of guanosine 3′,5′-cyclic monophosphate (cGMP) in single living cells, we constructed genetically encoded, fluorescent cGMP indicators by bracketing cGMP-dependent protein kinase (cGPK), minus residues 1–77, between cyan and yellow mutants of green fluorescent protein. cGMP decreased fluorescence resonance energy transfer (FRET) and increased the ratio of cyan to yellow emissions by up to 1.5-fold with apparent dissociation constants of ≈2 μM and >100:1 selectivity for cGMP over cAMP. To eliminate constitutive kinase activity, Thr516 of cGPK was mutated to Ala. Emission ratio imaging of the indicators transfected into rat fetal lung fibroblast (RFL)-6 showed cGMP transients resulting from activation of soluble and particulate guanylyl cyclase, respectively, by nitric oxide (NO) and C-type natriuretic peptide (CNP). Whereas all naive cells tested responded to CNP, only 68% responded to NO. Both sets of signals showed large and variable (0.5–4 min) latencies. The phosphodiesterase (PDE) inhibitor 3-isobutyl-1-methylxanthine (IBMX) did not elevate cGMP on its own but consistently amplified responses to NO or CNP, suggesting that basal activity of guanylate cyclase is very low and emphasizing the importance of PDEs in cGMP recycling. A fraction of RFL cells showed slowly propagating tides of cGMP spreading across the cell in response to delocalized application of NO. Biolistically transfected Purkinje neurons showed cGMP responses to parallel fiber activity and NO donors, confirming that single-cell increases in cGMP occur under conditions appropriate to cause synaptic plasticity.

Dynamics of Fos-Jun-NFAT1 complexes

Transcription initiation in eukaryotes is controlled by nucleoprotein complexes formed through cooperative interactions among multiple transcription regulatory proteins. These complexes may be assembled via stochastic collisions or defined pathways. We investigated the dynamics of Fos-Jun-NFAT1 complexes by using a multicolor fluorescence resonance energy transfer assay. Fos-Jun heterodimers can bind to AP-1 sites in two opposite orientations, only one of which is populated in mature Fos-Jun-NFAT1 complexes. We studied the reversal of Fos-Jun binding orientation in response to NFAT1 by measuring the efficiencies of energy transfer from donor fluorophores linked to opposite ends of an oligonucleotide to an acceptor fluorophore linked to one subunit of the heterodimer. The reorientation of Fos-Jun by NFAT1 was not inhibited by competitor oligonucleotides or heterodimers. The rate of Fos-Jun reorientation was faster than the rate of heterodimer dissociation at some binding sites. The facilitated reorientation of Fos-Jun heterodimers therefore can enhance the efficiency of Fos-Jun-NFAT1 complex formation. We also examined the influence of the preferred orientation of Fos-Jun binding on the stability and transcriptional activity of Fos-Jun-NFAT1 complexes. Complexes formed at sites where Fos-Jun favored the same binding orientation in the presence and absence of NFAT1 exhibited an 8-fold slower dissociation rate than complexes formed at sites where Fos-Jun favored the opposite binding orientation. Fos-Jun-NFAT1 complexes also exhibited greater transcription activation at promoter elements that favored the same orientation of Fos-Jun binding in the presence and absence of NFAT1. Thus, the orientation of heterodimer binding can influence both the dynamics and promoter selectivity of multiprotein transcription regulatory complexes.

Ultrafast diffraction and structural dynamics: The nature of complex molecules far from equilibrium

Studies of molecular structures at or near their equilibrium configurations have long provided information on their geometry in terms of bond distances and angles. Far-from-equilibrium structures are relatively unknown—especially for complex systems—and generally, neither their dynamics nor their average geometries can be extrapolated from equilibrium values. For such nonequilibrium structures, vibrational amplitudes and bond distances play a central role in phenomena such as energy redistribution and chemical reactivity. Ultrafast electron diffraction, which was developed to study transient molecular structures, provides a direct method for probing the nature of complex molecules far from equilibrium. Here we present our ultrafast electron diffraction observations of transient structures for two cyclic hydrocarbons. At high internal energies of ≈4 eV, these molecules display markedly different behavior. For 1,3,5-cycloheptatriene, excitation results in the formation of hot ground-state structures with bond distances similar to those of the initial structure, but with nearly three times the average vibrational amplitude. Energy is redistributed within 5 ps, but with a negative temperature characterizing the nonequilibrium population. In contrast, the ring-opening reaction of 1,3-cyclohexadiene is shown to result in hot structures with a C—C bond distance of over 1.7 Å, which is 0.2 Å away from any expected equilibrium value. Even up to 400 ps, energy remains trapped in large-amplitude motions comprised of torsion and asymmetric stretching. These studies promise a new direction for studying structural dynamics in nonequilibrium complex systems.