Configuration-dependent diffusion can shift the kinetic transition state and barrier height of protein folding

We show that diffusion can play an important role in protein-folding kinetics. We explicitly calculate the diffusion coefficient of protein folding in a lattice model. We found that diffusion typically is configuration- or reaction coordinate-dependent. The diffusion coefficient is found to be decreasing with respect to the progression of folding toward the native state, which is caused by the collapse to a compact state constraining the configurational space for exploration. The configuration- or position-dependent diffusion coefficient has a significant contribution to the kinetics in addition to the thermodynamic free-energy barrier. It effectively changes (increases in this case) the kinetic barrier height as well as the position of the corresponding transition state and therefore modifies the folding kinetic rates as well as the kinetic routes. The resulting folding time, by considering both kinetic diffusion and the thermodynamic folding free-energy profile, thus is slower than the estimation from the thermodynamic free-energy barrier with constant diffusion but is consistent with the results from kinetic simulations. The configuration- or coordinate-dependent diffusion is especially important with respect to fast folding, when there is a small or no free-energy barrier and kinetics is controlled by diffusion.Including the configurational dependence will challenge the transition state theory of protein folding.

The complex kinetics of protein folding in wide temperature ranges

The complex protein folding kinetics in wide temperature ranges is studied through diffusive dynamics on the underlying energy landscape. The well-known kinetic chevron rollover behavior is recovered from the mean first passage time, with the U-shape dependence on temperature. The fastest folding temperature T-0 is found to be smaller than the folding transition temperature T-f. We found that the fluctuations of the kinetics through the distribution of first passage time show rather universal behavior, from high-temperature exponential Poissonian kinetics to the relatively low-temperature highly nonexponential kinetics. The transition temperature is at T-k and T-0, T-k, T-f. In certain low-temperature regimes, a power law behavior at long time emerges. At very low temperatures ( lower than trapping transition temperature T< T-0/(4&SIM;6)), the kinetics is an exponential Poissonian process again.

Probing the kinetics of single molecule protein folding

We propose an approach to integrate the theory, simulations, and experiments in protein-folding kinetics. This is realized by measuring the mean and high-order moments of the first-passage time and its associated distribution. The full kinetics is revealed in the current theoretical framework through these measurements. In the experiments, information about the statistical properties of first-passage times can be obtained from the kinetic folding trajectories of single molecule experiments ( for example, fluorescence). Theoretical/simulation and experimental approaches can be directly related. We study in particular the temperature-varying kinetics to probe the underlying structure of the folding energy landscape. At high temperatures, exponential kinetics is observed; there are multiple parallel kinetic paths leading to the native state. At intermediate temperatures, nonexponential kinetics appears, revealing the nature of the distribution of local traps on the landscape and, as a result, discrete kinetic paths emerge. At very low temperatures, exponential kinetics is again observed; the dynamics on the underlying landscape is dominated by a single barrier.

Diffusion dynamics, moments, and distribution of first-passage time on the protein-folding energy landscape, with applications to single molecules

We study the dynamics of protein folding via statistical energy-landscape theory. In particular, we concentrate on the local-connectivity case with the folding progress described by the fraction of native conformations. We found that the first passage-time (FPT) distribution undergoes a dynamic transition at a temperature below which the FPT distribution develops a power-law tail, a signature of the intermittent nonexponential kinetic phenomena for the folding dynamics. Possible applications to single-molecule dynamics experiments are discussed.

The topological configuration and conformational analysis of mRNA in translation

The theoretical model construction of mRNA hairpin structure and single-stranded structure as well as the simulation studies on RNA structure determined by the X-ray crystal diffraction and nuclear magnetic resonance revealed that in translation, after mRNA being unfolded into single-stranded structure, its topological configuration was closely correlative with the original hairpin structure. The conformational features of single-stranded mRNA appeared as helical regions alternating with curly regions to different extents, which might exert the influence on the folding of nascent polypeptide by various regulating effects including different translational rates.

Exploring consensus mRNA secondary (folding) structure units by stochastic sampling and folding simulation

It is extremely difficult to explore mRNA folding structure by biological experiments. In this report, we use stochastic sampling and folding simulation to test the existence of the stable secondary structural units of-mRNA, look for the folding units, and explore the probabilistic stabilization of the units. Using this method, We made simulations for all possible local optimum secondary structures of a single strand mRNA within a certain range, and searched for the common parts of the secondary structures. The consensus secondary structure units (CSSUs) extracted from the above method are mainly hairpins, with a few single strands. These CSSUs suggest that the mRNA folding units could be relatively stable and could perform specific biological function. The significance of these observations for the mRNA folding problem in general is also discussed. (c) 2004 Elsevier B.V. All rights reserved.

First-passage time distribution and non-Markovian diffusion dynamics of protein folding

We study the kinetics of protein folding via statistical energy landscape theory. We concentrate on the local-connectivity case, where the configurational changes can only occur among neighboring states, with the folding progress described in terms of an order parameter given by the fraction of native conformations. The non-Markovian diffusion dynamics is analyzed in detail and an expression for the mean first-passage time (MFPT) from non-native unfolded states to native folded state is obtained. It was found that the MFPT has a V-shaped dependence on the temperature. We also find that the MFPT is shortened as one increases the gap between the energy of the native and average non-native folded states relative to the fluctuations of the energy landscape. The second- and higher-order moments are studied to infer the first-passage time distribution. At high temperature, the distribution becomes close to a Poisson distribution, while at low temperatures the distribution becomes a Levy-type distribution with power-law tails, indicating a nonself-averaging intermittent behavior of folding dynamics. We note the likely relevance of this result to single-molecule dynamics experiments, where a power law (Levy) distribution of the relaxation time of the underlined protein energy landscape is observed.

Temperature dependence of the distribution of the first passage time: Results from discontinuous molecular dynamics simulations of an all-atom model of the second beta-hairpin fragment of protein G

More than 22 000 folding kinetic simulations were performed to study the temperature dependence of the distribution of first passage time (FPT) for the folding of an all-atom Go-like model of the second beta-hairpin fragment of protein G. We find that the mean FPT (MFPT) for folding has a U (or V)-shaped dependence on the temperature with a minimum at a characteristic optimal folding temperature T-opt*. The optimal folding temperature T-opt* is located between the thermodynamic folding transition temperature and the solidification temperature based on the Lindemann criterion for the solid. Both the T-opt* and the MFPT decrease when the energy bias gap against nonnative contacts increases. The high-order moments are nearly constant when the temperature is higher than T-opt* and start to diverge when the temperature is lower than T-opt*. The distribution of FPT is close to a log-normal-like distribution at T* greater than or equal to T-opt*. At even lower temperatures, the distribution starts to develop long power-law-like tails, indicating the non-self-averaging intermittent behavior of the folding dynamics. It is demonstrated that the distribution of FPT can also be calculated reliably from the derivative of the fraction not folded (or fraction folded), a measurable quantity by routine ensemble-averaged experimental techniques at dilute protein concentrations.

The involvement of HSP22 from bay scallop Argopecten irradians in response to heavy metal stress

Heat shock protein 22 (HSP22) is an important member of small heat shock protein (sHSP) subfamily which plays a key role in the process of protecting cells, facilitating the folding of nascent peptides, and responding to stress. In the present study, the cDNA of HSP22 was cloned from Argopecten irradians (designated as AiHSP22) by rapid amplification cDNA end (RACE) based on the expressed sequence tags (ESTs). The full-length cDNA of AiHSP22 was of 1,112 bp, with an open reading frame of 588 bp encoding a polypeptide of 195 amino acids. The deduced amino acid sequence of AiHSP22 showed high similarity to previously identified HSP22s. The expression patterns of AiHSP22 mRNA in different tissues and in haemocytes of scallops exposed to Cd2+, Pb2+ or Cu2+ were investigated by real-time quantitative RT-PCR. The mRNA of AiHSP22 was constitutively expressed in all examined tissues, including haemocyte, muscle, kidney, gonad, gill and heart. The expression level in heart and muscle was higher than that in other tissues. The mRNA level of AiHSP22 in haemocytes was up-regulated after a 10 days exposure of scallops to Cu2+, Pb2+ and Cd2+. However, the expression of AiHSP22 did not increase linearly along with the rise of heavy metal concentration. Different concentrations of the same metal resulted in different effects on AiHSP22 expression. The sensitive response of AiHSP22 to Cu2+, Pb2+ and Cd2+ stress indicated that it could be developed as an indicator of exposure to heavy metals for the pollution monitoring programs in aquatic environment.

Study on the structure of C-phycocyanin in Spirulina platensis with scanning tunneling microscope

The C-phycocyanin (C-PC) trimmer was isolated from the blue-green alga Spirulina platensis, and scanning tunnelling microscope (STM) was used to investigate its structure, High resolution STM images of C-PC were obtained. From the STM images, it could be observed that the C-PC molecules were disk-like in shape and the subunits of C-PC arranged in ring-like pattern with a channel in the center. After filter treatment, the folding of the polypeptide chains could be. seen clearly. This is the first time to observe directly the topography of phycobiliprotein, and the results showed STM to be a powerful tool for the structural study of phycobiliproteins.

cDNA cloning and mRNA expression of heat shock protein 90 gene in the haemocytes of Zhikong scallop Chlamys farreri

Heat shock protein 90 (HSP90) is a highly conserved molecular chaperone contributing to the folding, maintenance of structural integrity and proper regulation of a subset of cytosolic proteins. The full-length cDNA of Zhikong scallop Chlamysfarreri HSP90 (designated CfHSP90) was cloned by EST and rapid RACE techniques. It was of 2710 bp, including an open reading frame (ORF) of 2181 bp encoding a polypeptide of 726 amino acids with all the five HSP90 family signatures. BLAST analysis revealed that the CfHSP90 gene shared high similarity with other known HSP90 genes. Fluorescent real-time quantitative RT-PCR was used to examine the expression pattern of CfHSP90 mRNA in haemocytes of scallops exposed to Cd2+, Pb2+ and Cu2+ for 10 and 20 days, respectively. All the three heavy metals could induce CfHSP90 expression. There was a clear dose-dependent expression pattern of CfHSP90 after heavy metals exposure for 10 days or 20 days. Different concentrations of the same metal resulted in different effects on CfHSP90 expression. The results indicated that CfHSP90 responded to various heavy metal stresses with a dose-dependent expression pattern as well as exposure time effect, and could be used as a molecular biomarker in a heavy metal polluted environment. (c) 2007 Elsevier Inc. All rights reserved.

Molecular cloning, characterization and expression of heat shock protein 90 gene in the haemocytes of bay scallop Argopecten irradians

Heat shock protein 90 (HSP90) is a highly conserved molecular chaperone that plays key roles in the folding, maintenance of structural integrity and regulation of a subset of cytosolic proteins. In the present study, the cDNA of Argopecten irradians HSP90 (designated AiHSP90) was cloned by the combination of homology cloning and rapid amplification of cDNA ends (RACE) approaches. The full-length cDNA of AiHSP90 was of 2669 bp, including an open reading frame (ORF) of 2175 bp encoding a polypeptide of 724 amino acids with predicted molecular weight of 83.08 kDa and theoretical isoelectric point of 4.81. BLAST analysis revealed that AiHSP90 shared high similarity with other known HSP90s, and the five conserved amino acid blocks defined as HSP90 protein family signatures were also identified in AiHSP90, which indicated that AiHSP90 should be a cytosolic member of the HSP90 family. Fluorescent real-time quantitative PCR was employed to examine the expression pattern of AiHSP90 mRNA in haemocytes of scallops challenged by Gram-negative bacteria Vibrio anguillarum and Gram-positive bacteria Micrococcus luteus. In both bacterial challenged groups, the relative expression level of AiHSP90 transcript was up-regulated and reached maximal. level at 9 h after injection, and then dropped progressively to the original level at about 48 h post challenge. The results indicated that AiHSP90 was potentially involved in the immune responses against bacteria challenge in scallop A. irradian. (c) 2007 Elsevier Ltd. All rights reserved.

The cDNA cloning and mRNA expression of heat shock protein 70 gene in the haemocytes of bay scallop (Argopecten irradians, Lamarck 1819) responding to bacteria challenge and naphthalin stress

Heat shock protein 70 (HSP70) is an important member of the heat shock protein superfamily, and it plays a key role in the process of protecting cells, facilitating the folding of nascent peptides and responding to stress. The cDNA of bay scallop Argopecten irradians HSP70 (designated AIHSP70) was cloned by the techniques of homological cloning and rapid amplification of cDNA end (RACE). The full length of AIHSP70 cDNA was 2651 bp in length, having a 5' untranslated region (UTR) of 96 bp, a 3' UTR of 575 bp, and an open reading frame (ORF) of 1980 bp encoding a polypeptide of 659 amino acids with an estimated molecular mass of 71.80 kDa and an estimated isoelectric point of 5.26. BLAST analysis revealed that the AIHSP70 gene shared high identity with other known HSP70 genes. Three classical HSP signature motifs were detected in AIHSP70 by InterPro, analysis. 3-D structural prediction of AIHSP70 showed that its N terminal ATPase activity domain and,C terminal substrate-binding domain shared high similarity with that in human heat shock protein 70. The results indicated that the AIHSP70 was a member of the heat shock protein 70 family. A semi-quantitive RT-PCR method was used to analyse the expression of AIHSP70 gene after the treatment of naphthalin which is one kind of polycyclic aromatic hydrocarbon (PAH) and the challenge of bacteria. mRNA expression of AIHSP70 in scallop was up-regulated significantly after the stimulation of naphthalin and increased with increasing naphthalin concentration. A clearly time-dependent expression pattern of AIHSP70 was observed after the scallops were infected by Vibrio anguillarum, and the mRNA expression reached a maximum level at 8 h and lasted to 16 h, and then dropped progressively. The results indicated that AIHSP70 could play an important role in mediating the environmental stress and immune response in scallop. (c) 2006 Elsevier Ltd. All rights reserved.

Switch Region for Pathogenic Structural Change in Conformational Disease and Its Prediction

Many diseases are believed to be related to abnormal protein folding. In the first step of such pathogenic structural changes, misfolding occurs in regions important for the stability of the native structure. This destabilizes the normal protein conformation, while exposing the previously hidden aggregation-prone regions, leading to subsequent errors in the folding pathway. Sites involved in this first stage can be deemed switch regions of the protein, and can represent perfect binding targets for drugs to block the abnormal folding pathway and prevent pathogenic conformational changes. In this study, a prediction algorithm for the switch regions responsible for the start of pathogenic structural changes is introduced. With an accuracy of 94%, this algorithm can successfully find short segments covering sites significant in triggering conformational diseases (CDs) and is the first that can predict switch regions for various CDs. To illustrate its effectiveness in dealing with urgent public health problems, the reason of the increased pathogenicity of H5N1 influenza virus is analyzed; the mechanisms of the pandemic swine-origin 2009 A(H1N1) influenza virus in overcoming species barriers and in infecting large number of potential patients are also suggested. It is shown that the algorithm is a potential tool useful in the study of the pathology of CDs because: (1) it can identify the origin of pathogenic structural conversion with high sensitivity and specificity, and (2) it provides an ideal target for clinical treatment.

mRNA 很可能携带三维遗传信息

The forming mechanism of the three - dimensional structures of proteins，i.e.the mechanism of protein folding，is a basic problem in molecular biology which is still unsolved unitl now. In which a core problem is whether there is the three – dimensional genetic information that decide the three - dimensional structures of proteins. However, the research on this field has mot yet been reported. Recently，we made a comparative study on the folded structures of more than 70 mature messeneger RNAs (mRNAs) and the three - dimensional structures of the proteins encoded by them，it has been found that there exist marked correspondences between their featured structures in the following aspects: 1.The number of the structural units. An RNA molecule can form a secondary structure(stem and loop structure) by the folding and the base pairing of itself. The elementary structural unit of an RNA secondary structure is hairpin(or compound hair pin).The regular structural unit in the secondary structure of a protein is # alpha # - helix or #beta# - sheet . We have found that the hairpin number in the secondary structure of each mature mRNA is equal or approximately equal to the number of the regular secondary structural unis of the encoded protein. 2 .Turning region. Turn is a main structrual element in the secondary structure of a protein, which decides the backbone orientation of a protein molecule to some extent .Our analysis shows that the nucleotide sequence segments in an mRNA which encode the turns of the corresponding protein are overall situated in the turning regions of the mRNA secondary structure such as haipin，bulge loop or multibaranch loops. 3 .The arrangement of structural elements in space. In order to understand the backbone orientation of an RNA molecule and the arangement of its structural elements in space，we have modeled the three一dimensional structure of the mRNA molecule on SGI workstation based on its secondary structure.The result shows that the spatial arrangement of most of the nucleotide sequence segments encoding the structural elements of a protein is consistent with that of these stretural exements in the protein. For instance，the nucleotide sequences corresponding to each pleated sheet of a # beta # - sheet structure are close to each other in the mRNA secondary stucture and in the three - dimensional structure，although some of the nucleotide segments are far apart from each other in the one - dimensional sequence. For another instance，the two triplet codons of cysteines which form a disulphide bridge geneal1y are very close to each other in the mRNA folded structure. In addition，we also analyzed the locations of the codons proline - coding and the distrbution of the nucleotide sequences #alpha# - helix - coding in the folded structures of mRNAs . Some distribution laws have been found. All of these results suggest that the transfer of the genetic information from mRNA to protein not only is one – dimensional but also is three - dime ns ional. That is，there exists the genetic information that decide the three - dimensional structures of proteins. To a certain extent，we could say that the mRNA folding detemines the protein folding. Based on these results，it would be possible to predict the three - dimensional structures of proteins from the primary，secondary and tertiary structures of the m RNAs at a higher accuracy.And more important is that a new clue has been provided to uncover the“spatial coding" of the genetic information.