Coordinate and time-dependent diffusion dynamics in protein folding

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Configuration-Dependent Diffusion Dynamics of Downhill and Two-State Protein Folding

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Effect of the thermostat in the molecular dynamics simulation on the folding of the model protein chignolin

Molecular dynamics simulations of the model protein chignolin with explicit solvent were carried out, in order to analyze the influence of the Berendsen thermostat on the evolution and folding of the peptide. The dependence of the peptide behavior on temperature was tested with the commonly employed thermostat scheme consisting of one thermostat for the protein and another for the solvent. The thermostat coupling time of the protein was increased to infinity, when the protein is not in direct contact with the thermal bath, a situation known as minimally invasive thermostat. In agreement with other works, it was observed that only in the last situation the instantaneous temperature of the model protein obeys a canonical distribution. As for the folding studies, it was shown that, in the applications of the commonly utilized thermostat schemes, the systems are trapped in local minima regions from which it has difficulty escaping. With the minimally invasive thermostat the time that the protein needs to fold was reduced by two to three times. These results show that the obstacles to the evolution of the extended peptide to the folded structure can be overcome when the temperature of the peptide is not directly controlled.

Optical projection tomography reveals dynamics of HEV growth after immunization with protein plus CFA and features shared with HEVs in acute autoinflammatory lymphadenopathy

The vascular-stromal compartment of lymph nodes is important for lymph node function, and high endothelial venules (HEVs) play a critical role in controlling the entry of recirculating lymphocytes. In autoimmune and autoinflammatory diseases, lymph node swelling is often accompanied by apparent HEV expansion and, potentially, targeting HEV expansion could be used therapeutically to limit autoimmunity. In previous studies using mostly flow cytometry analysis, we defined three differentially regulated phases of lymph node vascular-stromal growth: initiation, expansion, and the re-establishment of vascular quiescence and stabilization. In this study, we use optical projection tomography to better understand the morphologic aspects of HEV growth upon immunization with ovalbumin/CFA (OVA/CFA). We find HEV elongation as well as modest arborization during the initiation phase, increased arborization during the expansion phase, and, finally, vessel narrowing during the re-establishment of vascular quiescence and stabilization. We also examine acutely enlarged autoinflammatory lymph nodes induced by regulatory T cell depletion and show that HEVs are expanded and morphologically similar to the expanded HEVs in OVA/CFA-stimulated lymph nodes. These results reinforce the idea of differentially regulated, distinct phases of vascular-stromal growth after immunization and suggest that insights gained from studying immunization-induced lymph node vascular growth may help to understand how the lymph node vascular-stromal compartment could be therapeutically targeted in autoimmune and autoinflammatory diseases.

Characterization of the Effect of the Mitochondrial Protein Hint2 on Intracellular Ca2+ dynamics

Hint2, one of the five members of the superfamily of the histidine triad AMP-lysine hydrolase proteins, is expressed in mitochondria of various cell types. In human adrenocarcinoma cells, Hint2 modulates Ca2+ handling by mitochondria. As Hint2 is highly expressed in hepatocytes, we investigated if this protein affects Ca2+ dynamics in this cell type. We found that in hepatocytes isolated from Hint2−/− mice, the frequency of Ca2+ oscillations induced by 1 μM noradrenaline was 150% higher than in the wild-type. Using spectrophotometry, we analyzed the rates of Ca2+ pumping in suspensions of mitochondria prepared from hepatocytes of either wild-type or Hint2−/− mice; we found that Hint2 accelerates Ca2+ pumping into mitochondria. We then resorted to computational modeling to elucidate the possible molecular target of Hint2 that could explain both observations. On the basis of a detailed model for mitochondrial metabolism proposed in another study, we identified the respiratory chain as the most probable target of Hint2. We then used the model to predict that the absence of Hint2 leads to a premature opening of the mitochondrial permeability transition pore in response to repetitive additions of Ca2+ in suspensions of mitochondria. This prediction was then confirmed experimentally.

Unveiling Protein Functions through the Dynamics of the Interaction Network

Protein interaction networks have become a tool to study biological processes, either for predicting molecular functions or for designing proper new drugs to regulate the main biological interactions. Furthermore, such networks are known to be organized in sub-networks of proteins contributing to the same cellular function. However, the protein function prediction is not accurate and each protein has traditionally been assigned to only one function by the network formalism. By considering the network of the physical interactions between proteins of the yeast together with a manual and single functional classification scheme, we introduce a method able to reveal important information on protein function, at both micro- and macro-scale. In particular, the inspection of the properties of oscillatory dynamics on top of the protein interaction network leads to the identification of misclassification problems in protein function assignments, as well as to unveil correct identification of protein functions. We also demonstrate that our approach can give a network representation of the meta-organization of biological processes by unraveling the interactions between different functional classes

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.