Hydrated myoglobin's anharmonic fluctuations are not primarily due to dihedral transitions.

To characterize the functionally important anharmonic motions of proteins, simulations of carboxymyoglobin (MbCO) dynamics have been performed during which dihedral transitions were prohibited. Comparison of torsionally restrained and unrestrained protein dynamics simulated at three levels of hydration and at temperatures ranging from 100 to 400 K suggests that hydration "catalyzes" protein mobility by facilitating collective anharmonic motions that do not require dihedral transitions. When dihedral transitions were prohibited, dehydrated MbCO, to a good approximation, exhibited only harmonic fluctuations, whereas hydrated MbCO exhibited both harmonic and anharmonic motions. The fluctuation of helix centers of mass also remained highly anharmonic in the torsionally restrained hydrated system. Atomic mean-square fluctuation at 300 K was reduced upon prohibition of dihedral transitions by only 28% and 10% for MbCO hydrated by 350 and 3830 water molecules, respectively.

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.

Region-dependent dynamics of cAMP response element-binding protein phosphorylation in the basal ganglia

The cAMP response element-binding protein (CREB) is an activity-dependent transcription factor that is involved in neural plasticity. The kinetics of CREB phosphorylation have been suggested to be important for gene activation, with sustained phosphorylation being associated with downstream gene expression. If so, the duration of CREB phosphorylation might serve as an indicator for time-sensitive plastic changes in neurons. To screen for regions potentially involved in dopamine-mediated plasticity in the basal ganglia, we used organotypic slice cultures to study the patterns of dopamine- and calcium-mediated CREB phosphorylation in the major subdivisions of the striatum. Different durations of CREB phosphorylation were evoked in the dorsal and ventral striatum by activation of dopamine D1-class receptors. The same D1 stimulus elicited (i) transient phosphorylation (≤15 min) in the matrix of the dorsal striatum; (ii) sustained phosphorylation (≤2 hr) in limbic-related structures including striosomes, the nucleus accumbens, the fundus striati, and the bed nucleus of the stria terminalis; and (iii) prolonged phosphorylation (up to 4 hr or more) in cellular islands in the olfactory tubercle. Elevation of Ca2+ influx by stimulation of L-type Ca2+ channels, NMDA, or KCl induced strong CREB phosphorylation in the dorsal striatum but not in the olfactory tubercle. These findings differentiate the response of CREB to dopamine and calcium signals in different striatal regions and suggest that dopamine-mediated CREB phosphorylation is persistent in limbic-related regions of the neonatal basal ganglia. The downstream effects activated by persistent CREB phosphorylation may include time-sensitive neuroplasticity modulated by dopamine.

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.

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.

Identification of two distinct human SMC protein complexes involved in mitotic chromosome dynamics

The structural maintenance of chromosomes (SMC) family member proteins previously were shown to play a critical role in mitotic chromosome condensation and segregation in yeast and Xenopus. Other family members were demonstrated to be required for DNA repair in yeast and mammals. Although several different SMC proteins were identified in different organisms, little is known about the SMC proteins in humans. Here, we report the identification of four human SMC proteins that form two distinct heterodimeric complexes in the cell, the human chromosome-associated protein (hCAP)-C and hCAP-E protein complex (hCAP-C/hCAP-E), and the human SMC1 (hSMC1) and hSMC3 protein complex (hSMC1/hSMC3). The hCAP-C/hCAP-E complex is the human ortholog of the Xenopus chromosome-associated protein (XCAP)-C/XCAP-E complex required for mitotic chromosome condensation. We found that a second complex, hSMC1/hSMC3, is required for metaphase progression in mitotic cells. Punctate vs. diffuse distribution patterns of the hCAP-C/hCAP-E and hSMC1/hSMC3 complexes in the interphase nucleus indicate independent behaviors of the two complexes during the cell cycle. These results suggest that two distinct classes of SMC protein complexes are involved in different aspects of mitotic chromosome organization in human cells.

Dynamics of fluorescence fluctuations in green fluorescent protein observed by fluorescence correlation spectroscopy

We have investigated the pH dependence of the dynamics of conformational fluctuations of green fluorescent protein mutants EGFP (F64L/S65T) and GFP-S65T in small ensembles of molecules in solution by using fluorescence correlation spectroscopy (FCS). FCS utilizes time-resolved measurements of fluctuations in the molecular fluorescence emission for determination of the intrinsic dynamics and thermodynamics of all processes that affect the fluorescence. Fluorescence excitation of a bulk solution of EGFP decreases to zero at low pH (pKa = 5.8) paralleled by a decrease of the absorption at 488 nm and an increase at 400 nm. Protonation of the hydroxyl group of Tyr-66, which is part of the chromophore, induces these changes. When FCS is used the fluctuations in the protonation state of the chromophore are time resolved. The autocorrelation function of fluorescence emission shows contributions from two chemical relaxation processes as well as diffusional concentration fluctuations. The time constant of the fast, pH-dependent chemical process decreases with pH from 300 μs at pH 7 to 45 μs at pH 5, while the time-average fraction of molecules in a nonfluorescent state increases to 80% in the same range. A second, pH-independent, process with a time constant of 340 μs and an associated fraction of 13% nonfluorescent molecules is observed between pH 8 and 11, possibly representing an internal proton transfer process and associated conformational rearrangements. The FCS data provide direct measures of the dynamics and the equilibrium properties of the protonation processes. Thus FCS is a convenient, intrinsically calibrated method for pH measurements in subfemtoliter volumes with nanomolar concentrations of EGFP.

The Dynamics of Golgi Protein Traffic Visualized in Living Yeast Cells

We describe for the first time the visualization of Golgi membranes in living yeast cells, using green fluorescent protein (GFP) chimeras. Late and early Golgi markers are present in distinct sets of scattered, moving cisternae. The immediate effects of temperature-sensitive mutations on the distribution of these markers give clues to the transport processes occurring. We show that the late Golgi marker GFP-Sft2p and the glycosyltransferases, Anp1p and Mnn1p, disperse into vesicle-like structures within minutes of a temperature shift in sec18, sft1, and sed5 cells, but not in sec14 cells. This is consistent with retrograde vesicular traffic, mediated by the vesicle SNARE Sft1p, to early cisternae containing the target SNARE Sed5p. Strikingly, Sed5p itself moves rapidly to the endoplasmic reticulum (ER) in sec12 cells, implying that it cycles through the ER. Electron microscopy shows that Golgi membranes vesiculate in sec18 cells within 10 min of a temperature shift. These results emphasize the dynamic nature of Golgi cisternae and satisfy the kinetic requirements of a cisternal maturation model in which all resident proteins must undergo retrograde vesicular transport, either within the Golgi complex or from there to the ER, as anterograde cargo advances.

Dynamics and elasticity of the fibronectin matrix in living cell culture visualized by fibronectin–green fluorescent protein

Fibronectin (FN) forms the primitive fibrillar matrix in both embryos and healing wounds. To study the matrix in living cell cultures, we have constructed a cell line that secretes FN molecules chimeric with green fluorescent protein. These FN–green fluorescent protein molecules were assembled into a typical matrix that was easily visualized by fluorescence over periods of several hours. FN fibrils remained mostly straight, and they were seen to extend and contract to accommodate movements of the cells, indicating that they are elastic. When fibrils were broken or detached from cells, they contracted to less than one-fourth of their extended length, demonstrating that they are highly stretched in the living culture. Previous work from other laboratories has suggested that cryptic sites for FN assembly may be exposed by tension on FN. Our results show directly that FN matrix fibrils are not only under tension but are also highly stretched. This stretched state of FN is an obvious candidate for exposing the cryptic assembly sites.

Cell Cycle-Dependent Changes in Microtubule Dynamics in Living Cells Expressing Green Fluorescent Protein-α TubulinV⃞

LLCPK-1 cells were transfected with a green fluorescent protein (GFP)-α tubulin construct and a cell line permanently expressing GFP-α tubulin was established (LLCPK-1α). The mitotic index and doubling time for LLCPK-1α were not significantly different from parental cells. Quantitative immunoblotting showed that 17% of the tubulin in LLCPK-1α cells was GFP-tubulin; the level of unlabeled tubulin was reduced to 82% of that in parental cells. The parameters of microtubule dynamic instability were compared for interphase LLCPK-1α and parental cells injected with rhodamine-labeled tubulin. Dynamic instability was very similar in the two cases, demonstrating that LLCPK-1α cells are a useful tool for analysis of microtubule dynamics throughout the cell cycle. Comparison of astral microtubule behavior in mitosis with microtubule behavior in interphase demonstrated that the frequency of catastrophe increased twofold and that the frequency of rescue decreased nearly fourfold in mitotic compared with interphase cells. The percentage of time that microtubules spent in an attenuated state, or pause, was also dramatically reduced, from 73.5% in interphase to 11.4% in mitosis. The rates of microtubule elongation and rapid shortening were not changed; overall dynamicity increased 3.6-fold in mitosis. Microtubule release from the centrosome and a subset of differentially stable astral microtubules were also observed. The results provide the first quantitative measurements of mitotic microtubule dynamics in mammalian cells.

Symmetry, stability, and dynamics of multidomain and multicomponent protein systems

Symmetry is commonly observed in many biological systems. Here we discuss representative examples of the role of symmetry in structural molecular biology. Point group symmetries are observed in many protein oligomers whose three-dimensional atomic structures have been elucidated by x-ray crystallography. Approximate symmetry also occurs in multidomain proteins. Symmetry often confers stability on the molecular system and results in economical usage of basic components to build the macromolecular structure. Symmetry is also associated with cooperativity. Mild perturbation from perfect symmetry may be essential in some systems for dynamic functions.

Quantitative ER ↔ Golgi Transport Kinetics and Protein Separation upon Golgi Exit Revealed by Vesicular Integral Membrane Protein 36 Dynamics in Live CellsV⃞

To quantitatively investigate the trafficking of the transmembrane lectin VIP36 and its relation to cargo-containing transport carriers (TCs), we analyzed a C-terminal fluorescent-protein (FP) fusion, VIP36-SP-FP. When expressed at moderate levels, VIP36-SP-FP localized to the endoplasmic reticulum, Golgi apparatus, and intermediate transport structures, and colocalized with epitope-tagged VIP36. Temperature shift and pharmacological experiments indicated VIP36-SP-FP recycled in the early secretory pathway, exhibiting trafficking representative of a class of transmembrane cargo receptors, including the closely related lectin ERGIC53. VIP36-SP-FP trafficking structures comprised tubules and globular elements, which translocated in a saltatory manner. Simultaneous visualization of anterograde secretory cargo and VIP36-SP-FP indicated that the globular structures were pre-Golgi carriers, and that VIP36-SP-FP segregated from cargo within the Golgi and was not included in post-Golgi TCs. Organelle-specific bleach experiments directly measured the exchange of VIP36-SP-FP between the Golgi and endoplasmic reticulum (ER). Fitting a two-compartment model to the recovery data predicted first order rate constants of 1.22 ± 0.44%/min for ER → Golgi, and 7.68 ± 1.94%/min for Golgi → ER transport, revealing a half-time of 113 ± 70 min for leaving the ER and 1.67 ± 0.45 min for leaving the Golgi, and accounting for the measured steady-state distribution of VIP36-SP-FP (13% Golgi/87% ER). Perturbing transport with AlF4− treatment altered VIP36-SP-GFP distribution and changed the rate constants. The parameters of the model suggest that relatively small differences in the first order rate constants, perhaps manifested in subtle differences in the tendency to enter distinct TCs, result in large differences in the steady-state localization of secretory components.

Ultra-fast excited state dynamics in green fluorescent protein: multiple states and proton transfer.

The green fluorescent protein (GFP) of the jellyfish Aequorea Victoria has attracted widespread interest since the discovery that its chromophore is generated by the autocatalytic, posttranslational cyclization and oxidation of a hexapeptide unit. This permits fusion of the DNA sequence of GFP with that of any protein whose expression or transport can then be readily monitored by sensitive fluorescence methods without the need to add exogenous fluorescent dyes. The excited state dynamics of GFP were studied following photo-excitation of each of its two strong absorption bands in the visible using fluorescence upconversion spectroscopy (about 100 fs time resolution). It is shown that excitation of the higher energy feature leads very rapidly to a form of the lower energy species, and that the excited state interconversion rate can be markedly slowed by replacing exchangeable protons with deuterons. This observation and others lead to a model in which the two visible absorption bands correspond to GFP in two ground-state conformations. These conformations can be slowly interconverted in the ground state, but the process is much faster in the excited state. The observed isotope effect suggests that the initial excited state process involves a proton transfer reaction that is followed by additional structural changes. These observations may help to rationalize and motivate mutations that alter the absorption properties and improve the photo stability of GFP.

From residue matching patterns to protein folding topographies: General model and bovine pancreatic trypsin inhibitor

A coarse-grained model for protein-folding dynamics is introduced based on a discretized representation of torsional modes. The model, based on the Ramachandran map of the local torsional potential surface and the class (hydrophobic/polar/neutral) of each residue, recognizes patterns of both torsional conformations and hydrophobic-polar contacts, with tolerance for imperfect patterns. It incorporates empirical rates for formation of secondary and tertiary structure. The method yields a topological representation of the evolving local torsional configuration of the folding protein, modulo the basins of the Ramachandran map. The folding process is modeled as a sequence of transitions from one contact pattern to another, as the torsional patterns evolve. We test the model by applying it to the folding process of bovine pancreatic trypsin inhibitor, obtaining a kinetic description of the transitions between the contact patterns visited by the protein along the dominant folding pathway. The kinetics and detailed balance make it possible to invert the result to obtain a coarse topographic description of the potential energy surface along the dominant folding pathway, in effect to go backward or forward between a topological representation of the chain conformation and a topographical description of the potential energy surface governing the folding process. As a result, the strong structure-seeking character of bovine pancreatic trypsin inhibitor and the principal features of its folding pathway are reproduced in a reasonably quantitative way.