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.

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.

Simulated pathogenic conformational switch regions matched well with the biochemical findings

Pathogenic conformational conversion is a general causation of many disease, such as transmissible spon- giform encephalopathy (TSE) caused by misfolding of prion, sickle cell anemia, and etc. In such structural changes, misfolding occurs in regions important for the stability of native structure firstly. This destabi- lizes the normal conformation and leads to subsequent errors in folding pathway. Sites involved in the first stage can be deemed switch regions of the protein, and are vital for conformational conversion. Namely it could be a switch of disease at residue level. Here we report an algorithm that can identify such sites computationally with an accuracy of 93%, by calculating the probability of the native structure of a short segment jumping to a mistake one. Knowledge of such switch sites could be used to target clinical therapy, study physiological and pathologic mechanism of protein, and etc.

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.

Generating artificial homologous proteins according to the representative family character in molecular mechanics properties - an attempt in validating an underlying rule of protein evolution

The molecular mechanics property is the foundation of many characters of proteins. Based on intramolecular hydrophobic force network, the representative family character underlying a protein’s mechanics property is described by a simple two-letter scheme. The tendency of a sequence to become a member of a protein family is scored according to this mathematical representation. Remote homologs of the WW-domain family could be easily designed using such a mechanistic signature of protein homology. Experimental validation showed that nearly all artificial homologs have the representative folding and bioactivity of their assigned family. Since the molecular mechanics property is the only consideration in this study, the results indicate its possible role in the generation of new members of a protein family during evolution.

Optical two-dimensional SW-banyan network: Optical implementation, routing control, and determination of permissible permutations

A 2-D SW-banyan network is introduced by properly folding the 1-D SW-banyan network, and its corresponding optical setup is proposed by means of polarizing beamsplitters and 2-D phase spatial light modulators. Then, based on the characteristics and the proposed optical setup, the control for the routing path between any source-destination pair is given, and the method to determine whether a given permutation is permissible or not is discussed. Because the proposed optical setup consists of only optical polarization elements, it is compact in structure, its corresponding energy loss and crosstalk are low, and its corresponding available number of channels is high. (C) 1996 Society of Photo-Optical Instrumentation Engineers.

Perspectives on the origin of microfilaments, microtubules, the relevant chaperonin system and cytoskeletal motors - a commentary on the spirochaete origin of flagella

The origin of cytoskeleton and the origin of relevant intracellular transportation system are big problems for understanding the emergence of eukaryotic cells. The present article summarized relevant information of evidences and molecular traces on the origin of actin, tubulin, the chaperonin system for folding them, myosins, kinesins, axonemal dyneins and cytoplasmic dyneins. On this basis the authors proposed a series of works, which should be done in the future, and indicated the ways for reaching the targets. These targets are mainly: 1) the reconstruction of evolutionary path from MreB protein of archaeal ancestor of eukaryotic cells to typical actin; 2) the finding of the MreB or MreB-related proteins in crenarchaea and using them to examine J. A. Lake's hypothesis on the origin of eukaryote from "eocytes" (crenarchaea); 3) the examinations of the existence and distribution of cytoskeleton made of MreB-related protein within coccoid archaea, especially in amoeboid archaeon Thermoplasm acidophilum; 4) using Thermoplasma as a model of archaeal ancestor of eukaryotic cells; 5) the searching for the homolog of ancestral dynein in present-day living archaea. During the writing of this article, Margulis' famous spirochaete hypothesis on the origin of flagella and cilia was unexpectedly involved and analyzed from aspects of tubulins, dyneins and spirochaetes. Actually, spirochaete cannot be reasonably assumed as the ectosymbiotic ancestor of eukaryotic flagella and cilia, since their swing depends upon large amount of bacterial flagella beneath the flexible outer wall, but not depends upon their intracellular tubules and the assumed dyneins. In this case, if they had "evolved" into cilia and lost their bacterial flagella, they would immediately become immobile! In fact, tubulin and dynein-like proteins have not been found in any spirochaete.

Secondary structural analysis of the mRNA regions encoding the hemagglutinin cleavage site basic amino acids of the avian influenza virus H5N1 subtype samples

Here we report the codon bias and the mRNA secondary structural features of the hemagglutinin (HA) cleavage site basic amino acid regions of avian influenza virus H5N1 subtypes. We have developed a dynamic extended folding strategy to predict RNA secondar

大肠杆菌mRNA 编码区长度、形成二级结构倾向与密码子偏好性的关系

从GenBank获得大肠杆菌K-12MG1655株的全基因组序列,计算了与基因密码子偏好性相关的多个参数(Nc、CAI、GC、GC3s),对其mRNA编码区长度、形成二级结构倾向与密码子偏好性之间的关系进行了统计学分析,发现虽然翻译效率(包括翻译速度和翻译精度)是制约大肠杆菌高表达基因的密码子偏好性的主要因素,同时,mRNA编码区长度及其形成二级结构的倾向也是形成这种偏好性的不可忽略的原因,而且对偏好性有一定程度的削弱。另外对mRNA编码区形成二级结构倾向的生物学意义进行了讨论分析。

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.

Secondary structures of mRNAs and splicing of introns in eukaryotes.

After analyzing the secondary structures of 68 exon-intron-exon and the corresponding exon-exon sequence segments, it is found that about 90% of 5' and 3' terminal bases G (splicing sites) of introns are situated in the loops of secondary structures or at the ends of stems near the loops, and most of "G" s in loops are closed to the ends of loops. Approximately 92% of the connecting sites of the adjoining exons also show the similar features. About 82% of the branch point "A" s are situated in loops or at the ends of stems near the loops. Splicing sites and branch points approach each other in space because of the folding.

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.

Stable RNA hairpins in 88 coding regions of human mRNA

RNA hairpins containing UNCG, GNRA, CUUG (N = A, U, C or G, R = G or A) loops are unusually thermodynamic stable and conserved structures. The structural features of these hairpin loops are very special, and they play very important roles in vivo. They are prevalent in rRNA, catalytic RNA and non-coding mRNA. However, the 5' C(UUCG)G 3' hairpin is not found in the folding structure of 88 human mRNA coding regions. It is also different from rRNA in that there is no preference for certain sequences among tetraloops in these 88 mRNA folding structures.

Mutations stabilize small subunit ribosomal RNA in desiccation-tolerant cyanobacteria Nostoc

The ribosomal RNA molecule is an ideal model for evaluating the stability of a gene product under desiccation stress. We isolated 8 Nostoc strains that had the capacity to withstand desiccation in habitats and sequenced their 16S rRNA genes. The stabilities of 16S rRNAs secondary structures, indicated by free energy change of folding, were compared among Nostoc and other related species. The results suggested that 163 rRNA secondary structures of the desiccation-tolerant Nostoc strains were more stable than that of planktonic Nostocaceae species. The stabilizing mutations were divided into two categories: (1) those causing GC to replace other types of base pairs in stems and (2) those causing extension of stems. By mapping stabilizing mutations onto the Nostoc phylogenetic tree based on 16S rRNA gene, it was shown that most of stabilizing mutations had evolved during adaptive radiation among Nostoc spp. The evolution of 16S rRNA along the Nostoc lineage is suggested to be selectively advantageous under desiccation stress.

A full spd tight-binding treatment for electronic bands of graphitic tubes

A tight-binding (TB) treatment with the inclusion of d orbitals is applied to the electronic structures of graphitic tubes. The results show that the high angular moment bases in TB scheme are necessary to account the severe curvature effect in ultra-thin single wall carbon nanotubes, especially for properly reproducing the band edge overlap behavior in (5, 0) tube, predicted by the existing ab initio calculations. In the large diameter limit, the participation of two synnmetry-allowed d bases provides a natural replication to the recent measured electronic dispersions of valence band of graphene when the strong anisotropy due to the two-dimensional planar hexagonal sheet structure is dealt with properly. In addition, the detailed relation between the two sets of quantum numbers of screw symmetry and that of zone folding is formulated in appendix. (C) 2008 Elsevier Ltd. All rights reserved.