Structural characterization of (3-mercaptopropyl)sulfonate monolayer on gold surfaces

We have investigated the structure of (3-mercaptopropyl)sulfonate (MPS) monolayer self-assembled onto gold surfaces by quartz crystal microbalance with energy dissipation monitoring (QCM-D) and various electrochemical methods. QCM-D experiments show that the MPS monolayer behaves as a thin rigid film with a surfacic mass of 166 ng cm(-2). Interfacial capacitance measurements demonstrate that the MPS monolayer is a rather open structure that can be penetrated by the ionic species of the phosphate buffer electrolyte. From MPS reductive desorption experiments, MPS surface concentration corresponds to 4.6 x 10(-10) mol cm(-2), which represents 60% of the coverage reported for a densely packed thiol monolayer. Despite this low packing density, oxidation of catechol is strongly inhibited leading to voltammograms that are free ofdiffusional contribution. This unique behavior has been exploited to show that the MPS monolayer covers the entire gold surface with a surface coverage at least equal to theta = 0.9981, which means a very low number of MPS-free pinholes and/or defects. Kinetics of electron transfer toward soluble redox species has been studied using catechol as a neutral hydrophilic probe, but also ferrocyanide as hydrophilic anion and ferrocenemethanol as neutral hydrophobic molecule. It is proposed that the MPS monolayer provides a high kinetic barrier toward permeation of these species and that electron transfer mainly occurs by electron tunneling through the MPS monolayer.     

Correlations Between Amino Acids at Different Sites in Local Sequences of Protein Fragments with Given Structural Patterns

Ample evidence suggests that the local structures of peptide fragments in native proteins are to some extent encoded by their local sequences. Detecting such local correlations is important but it is still an open question what would be the most appropriate method. This is partly because conventional sequence analyses treat amino acid preferences at each site of a protein sequence independently, while it is often the inter-site interactions that bring about local sequence-structure correlations. Here a new scheme is introduced to capture the correlation between amino acid preferences at different sites for different local structure types. A library of nine-residue fragments is constructed, and the fragments are divided into clusters based on their local structures. For each local structure cluster or type, chi-square tests are used to identify correlated preferences of amino acid combinations at pairs of sites. A score function is constructed including both the single site amino acid preferences and the dual-site amino acid combination preferences, which can be used to identify whether a sequence fragment would have a strong tendency to form a particular local structure in native proteins. The results show that, given a local structure pattern, dual-site amino acid combinations contain different information from single site amino acid preferences. Representative examples show that many of the statistically identified correlations agree with previously-proposed heuristic rules about local sequence-structure correlations, or are consistent with physical-chemical interactions required to stabilize particular local structures. Results also show that such dual-site correlations in the score function significantly improves the Z-score matching a sequence fragment to its native local structure relative to nonnative local structures, and certain local structure types are highly predictable from the local sequence alone if inter-site correlations are considered.     

Design of amino acid sequences to fold into C(alpha)-model proteins

In order to extend the results obtained with minimal lattice models to more realistic systems, we study a model where proteins are described as a chain of 20 kinds of structureless amino acids moving in a continuum space and interacting through a contact potential controlled by a 20x20 quenched random matrix. The goal of the present work is to design and characterize amino acid sequences folding to the SH3 conformation, a 60-residue recognition domain common to many regulatory proteins. We show that a number of sequences can fold, starting from a random conformation, to within a distance root-mean-square deviation between 2.6 and 4.0 A from the native state. Good folders are those sequences displaying in the native conformation an energy lower than a sequence-independent threshold energy.     

A computational study of the correlations between structure and dynamics in free and surface-immobilized single polymer chains

The correlations between structure and dynamics in free and surface-immobilized polymers were investigated via Langevin dynamics simulations of a free-jointed homopolymer. A detailed analysis was performed for a polymer in free solution and a polymer attached to a surface. The cases of repulsive and attractive surfaces, as well as poor and good solvents, were considered. The analysis focuses on properties that are particularly relevant to single molecule measurements, namely: (1) the distribution of end-to-end distance, (2) the correlations between the conformational structure and the time scale of its motion, (3) the correlations, at equilibrium, between the end-to-end distance and its displacement, and (4) the correlation between the initial coil conformation and the collapse pathway into the globular state. The differences and similarities between this model and a previously considered model of a protein, with two-state folding kinetics and a well-defined native state, are also discussed.     

Mechanical properties of the domains of titin in a Go-like model

Comparison of properties of three domains of titin, I1, I27, and I28, in a simple geometry-based model shows that despite a high structural homology between their native states different domains show similar but distinguishable mechanical properties. Folding properties of the separate domains are predicted to be diversified which reflects sensitivity of the kinetics to the details of native structures. The Go-like model corresponding to the experimentally resolved native structure of the I1 domain is found to provide the biggest thermodynamic and mechanical stability compared to the other domains studied here. We analyze elastic, thermodynamic, and kinetic properties of several structures corresponding to the I28 domain as obtained through homology-based modeling. We discuss the ability of the models of the I28 domain to reproduce experimental results qualitatively. A strengthening of contacts that involve hydrophobic amino acids does not affect theoretical comparisons of the domains. Tandem linkages of up to five identical or different domains unravel in a serial fashion at low temperatures. We study the nature of the intermediate state that arises in the early stages of the serial unraveling and find it to qualitatively agree with the results of Marszalek et al.     

Frustration and hydrophobicity interplay in protein folding and protein evolution

A lattice model is used to study mutations and compacting effects on protein folding rates and folding temperature. In the context of protein evolution, we address the question regarding the best scenario for a polypeptide chain to fold: either a fast nonspecific collapse followed by a slow rearrangement to form the native structure or a specific collapse from the unfolded state with the simultaneous formation of the native state. This question is investigated for optimized sequences, whose native state has no frustrated contacts between monomers, and also for mutated sequences, whose native state has some degree of frustration. It is found that the best scenario for folding may depend on the amount of frustration of the native structure. The implication of this result on protein evolution is discussed.     

二维HP格点模型中的紧密高分子链构象及其热力学性质的研究

采用二维HP模型用精确计数法和MonteCarlo方法研究了链长为N(≤ 2 2 )的紧密高分子链的构象和热力学性质 .发现不同HP序列的紧密高分子链的平均自由能和平均配分函数与链长N存在关系 :〈F〉=aN+b ,　ln〈Z〉=a′N +b′ .同时发现对于可折叠成基态且简并度为 1的紧密高分子链 ,其平均自由能和平均配分函数与链长N也存在相似的关系 .在HP模型中对于链长为N的紧密高分子链 ,存在着 2 N + 1 个不同的HP序列 .我们发现可以折叠成基态且简并度为 1的蛋白质分子的HP序列数目NS 为NS =a× 2 N+ 1 　 (a =0 0 2 5 ) ,对应的HP序列中 ,疏水基团 (H)数目的含量为 4 0 %～ 6 0 %的序列出现的几率最大 .同时在这些紧密高分子链中有些具有相同的结构 ,发现结构的‘简并度’为 3 3～ 4 0 (10≤N≤ 16 ) .在紧密高分子链折叠过程中 ,折叠的初期能量下降比较快 ,折叠的中期能量下降比较缓慢 ,折叠的后期能量下降也是比较快     

Understanding the folding rates and folding nuclei of globular proteins

The first part of this paper contains an overview of protein structures, their spontaneous formation ("folding"), and the thermodynamic and kinetic aspects of this phenomenon, as revealed by in vitro experiments. It is stressed that universal features of folding are observed near the point of thermodynamic equilibrium between the native and denatured states of the protein. Here the "two-state" ("denatured state" <--> "native state") transition proceeds without accumulation of metastable intermediates, but includes only the unstable "transition state". This state, which is the most unstable in the folding pathway, and its structured core (a "nucleus") are distinguished by their essential influence on the folding/unfolding kinetics. In the second part of the paper, a theory of protein folding rates and related phenomena is presented. First, it is shown that the protein size determines the range of a protein's folding rates in the vicinity of the point of thermodynamic equilibrium between the native and denatured states of the protein. Then, we present methods for calculating folding and unfolding rates of globular proteins from their sizes, stabilities and either 3D structures or amino acid sequences. Finally, we show that the same theory outlines the location of the protein folding nucleus (i.e., the structured part of the transition state) in reasonable agreement with experimental data.     

Effective interactions in lysozyme aqueous solutions: a small-angle neutron scattering and computer simulation study

We report protein-protein structure factors of aqueous lysozyme solutions at different pH and ionic strengths, as determined by small-angle neutron scattering experiments. The observed upturn of the structure factor at small wavevectors, as the pH increases, marks a crossover between two different regimes, one dominated by repulsive forces, and another one where attractive interactions become prominent, with the ensuing development of enhanced density fluctuations. In order to rationalize such experimental outcome from a microscopic viewpoint, we have carried out extensive simulations of different coarse-grained models. We have first studied a model in which macromolecules are described as soft spheres interacting through an attractive r(-6) potential, plus embedded pH-dependent discrete charges; we show that the uprise undergone by the structure factor is qualitatively predicted. We have then studied a Derjaguin-Landau-Verwey-Overbeek (DLVO) model, in which only central interactions are advocated; we demonstrate that this model leads to a protein-rich/protein-poor coexistence curve that agrees quite well with the experimental counterpart; experimental correlations are instead reproduced only at low pH and ionic strengths. We have finally investigated a third, "mixed" model in which the central attractive term of the DLVO potential is imported within the distributed-charge approach; it turns out that the different balance of interactions, with a much shorter-range attractive contribution, leads in this latter case to an improved agreement with the experimental crossover. We discuss the relationship between experimental correlations, phase coexistence, and features of effective interactions, as well as possible paths toward a quantitative prediction of structural properties of real lysozyme solutions.     

Efficient quantification of the importance of contacts for the dynamical stability of proteins

Understanding the stability of the native state and the dynamics of a protein is of great importance for all areas of biomolecular design. The efficient estimation of the influence of individual contacts between amino acids in a protein structure is a first step in the reengineering of a particular protein for technological or pharmacological purposes. At the same time, the functional annotation of molecular evolution can be facilitated by such insight. Here, we use a recently suggested, information theoretical measure in biomolecular design - the Kullback-Leibler-divergence - to quantify and therefore rank residue-residue contacts within proteins according to their overall contribution to the molecular mechanics. We implement this protocol on the basis of a reduced molecular model, which allows us to use a well-known lemma of linear algebra to speed up the computation. The increase in computational performance is around 10(1)- to 10(4)-fold. We applied the method to two proteins to illustrate the protocol and its results. We found that our method can reliably identify key residues in the molecular mechanics and the protein fold in comparison to well-known properties in the serine protease inhibitor. We found significant correlations to experimental results, e.g., dissociation constants and Φ values.     

A designed well-folded monomeric four-helix bundle protein prepared by Fmoc solid-phase peptide synthesis and native chemical ligation

The design and total chemical synthesis of a monomeric native-like four-helix bundle protein is presented. The designed protein, GTD-Lig, consists of 90 amino acids and is based on the dimeric structure of the de novo designed helix-loop-helix GTD-43. GTD-Lig was prepared by the native chemical ligation strategy and the fragments (45 residues long) were synthesized by applying standard fluorenylmethoxycarbonyl (Fmoc) chemistry. The required peptide-thioester fragment was prepared by anchoring the free gamma-carboxy group of Fmoc-Glu-allyl to the solid phase. After chain elongation the allyl moiety was orthogonally removed and the resulting carboxy group was functionalized with a glycine-thioester followed by standard trifluoroacetic acid (TFA) cleavage to produce the unprotected peptide-thioester. The structure of the synthetic protein was examined by far- and near-UV circular dichroism (CD), sedimentation equilibrium ultracentrifugation, and NMR and fluorescence spectroscopy. The spectroscopic methods show a highly helical and native-like monomeric protein consistent with the design. Heat-induced unfolding was studied by tryptophan absorbance and far-UV CD. The thermal unfolding of GTD-Lig occurs in two steps; a cooperative transition from the native state to an intermediate state and thereafter by noncooperative melting to the unfolded state. The intermediate exhibits the properties of a molten globule such as a retained native secondary structure and a compact hydrophobic core. The thermodynamics of GuHCl-induced unfolding were evaluated by far-UV CD monitoring and the unfolding exhibited a cooperative transition that is well-fitted by a two-state mechanism from the native to the unfolded state. GTD-Lig clearly shows the characteristics of a native protein with a well-defined structure and typical unfolding transitions. The design and synthesis presented herein is of general applicability for the construction of large monomeric proteins.     

A lattice protein with an amyloidogenic latent state: stability and folding kinetics

We have designed a model lattice protein that has two stable folded states, the lower free energy native state and a latent state of somewhat higher energy. The two states have a sizable part of their structures in common (two "alpha-helices") and differ in the content of "alpha-helices" and "beta-strands" in the rest of their structures; i.e. for the native state, this part is alpha-helical, and for the latent state it is composed of beta-strands. Thus, the lattice protein free energy surface mimics that of amyloidogenic proteins that form well organized fibrils under appropriate conditions. A Go-like potential was used and the folding process was simulated with a Monte Carlo method. To gain insight into the equilibrium free energy surface and the folding kinetics, we have combined standard approaches (reduced free energy surfaces, contact maps, time-dependent populations of the characteristic states, and folding time distributions) with a new approach. The latter is based on a principal coordinate analysis of the entire set of contacts, which makes possible the introduction of unbiased reaction coordinates and the construction of a kinetic network for the folding process. The system is found to have four characteristic basins, namely a semicompact globule, an on-pathway intermediate (the bifurcation basin), and the native and latent states. The bifurcation basin is shallow and consists of the structure common to the native and latent states, with the rest disorganized. On the basis of the simulation results, a simple kinetic model describing the transitions between the characteristic states was developed, and the rate constants for the essential transitions were estimated. During the folding process the system dwells in the bifurcation basin for a relatively short time before it proceeds to the native or latent state. We suggest that such a bifurcation may occur generally for proteins in which native and latent states have a sizable part of their structures in common. Moreover, there is the possibility of introducing changes in the system (e.g., mutations), which guide the system toward the native or misfolded state.     

Evidence of a two-stage thermal denaturation process in lysozyme: a Raman scattering and differential scanning calorimetry investigation

Raman spectroscopy (in the low-frequency range and the amide I band region) and modulated differential scanning calorimetry investigations have been used to analyze temperature-induced structural changes in lysozyme dissolved in 1H2O and 2H2O in the thermal denaturation process. Low-frequency Raman data reveal a change in tertiary structure without concomitant unfolding of the secondary structure. Calorimetric data show that this structural change is responsible for the configurational entropy change associated with the strong-to-fragile liquid transition and correspond to about 1/3 of the native-denaturated transition enthalpy. This is the first stage of the thermal denaturation which is a precursor of the secondary structure change and is determined to be strongly dependent on the stability of the hydrogen-bond network in water. Low-frequency Raman spectroscopy provides information on the flexibility of the tertiary structure (in the native state and the transient folding state) in relation to the fragility of the mixture. The unfolding of the secondary structure appears as a consequence of the change in the tertiary structure and independent of the solvent. Protein conformational stability is directly dependent on the stability of the native tertiary structure. The structural transformation of tertiary structure can be detected through the enhanced 1H/2H exchange inhibited in native proteins. Taking into account similar features reported in the literature observed for different proteins it can be considered that the two-stage transformation observed in lysozyme dissolved in water is a general mechanism for the thermal denaturation of proteins.     

Glassy phases in random heteropolymers with correlated sequences

We develop an analytic approach for the study of lattice heteropolymers and apply it to copolymers with correlated Markovian sequences. According to our analysis, heteropolymers present three different dense phases depending upon the temperature, the nature of the monomer interactions, and the sequence correlations: (i) a liquid phase, (ii) a "soft glass" phase, and (iii) a "frozen glass" phase. The presence of the intermediate "soft glass" phase is predicted, for instance, in the case of polyampholytes with sequences that favor the alternation of monomers. Our approach is based on the cavity method, a refined Bethe-Peierls approximation adapted to frustrated systems. It amounts to a mean-field treatment in which the nearest-neighbor correlations, which are crucial in the dense phases of heteropolymers, are handled exactly. This approach is powerful and versatile; it can be improved systematically and generalized to other polymeric systems.     

Morphological phase separation in unstable thin films: pattern formation and growth

We present results from a comprehensive numerical study of morphological phase separation (MPS) in unstable thin liquid films on a 2-dimensional substrate. We study the quantitative properties of the evolution morphology via several experimentally relevant markers, e.g., correlation function, structure factor, domain-size and defect-size probability distributions, and growth laws. Our results suggest that the late-stage morphologies exhibit dynamical scaling, and their evolution is self-similar in time. We emphasize the analogies and differences between MPS in films and segregation kinetics in unstable binary mixtures.     

Enhanced Wang Landau sampling of adsorbed protein conformations

Using computer simulations to model the folding of proteins into their native states is computationally expensive due to the extraordinarily low degeneracy of the ground state. In this paper, we develop an efficient way to sample these folded conformations using Wang Landau sampling coupled with the configurational bias method (which uses an unphysical "temperature" that lies between the collapse and folding transition temperatures of the protein). This method speeds up the folding process by roughly an order of magnitude over existing algorithms for the sequences studied. We apply this method to study the adsorption of intrinsically disordered hydrophobic polar protein fragments on a hydrophobic surface. We find that these fragments, which are unstructured in the bulk, acquire secondary structure upon adsorption onto a strong hydrophobic surface. Apparently, the presence of a hydrophobic surface allows these random coil fragments to fold by providing hydrophobic contacts that were lost in protein fragmentation.     

Sensing Picomolar Concentrations of RNA Using Switchable Plasmonic Chirality

Detecting small sequences of RNA in biological samples such as microRNA or viral RNA demands highly sensitive and specific methods. Here, a reconfigurable DNA origami template has been used where a chiral arrangement of gold nanorods on the structure can lead to the generation of strong circular dichroism (CD). Switching of the cross‐like DNA structure is achieved by the addition of nucleic acid sequences, which arrests the structure in one of the possible chiral states by specific molecular recognition. A specific sequence can thus be detected through the resulting changes in the plasmonic CD spectrum. We show the sensitive and selective detection of a target RNA sequence from the hepatitis C virus genome. The RNA binds to a complementary sequence that is part of the lock mechanism, which leads to the formation of a defined state of the plasmonic system with a distinct optical response. With this approach, we were able to detect this specific RNA sequence at concentrations as low as 100 pm .     

Lagrange multiplier based transport theory for quantum wires

We discuss how the Lagrange multiplier method of nonequilibrium steady state statistical mechanics can be applied to describe the electronic transport in a quantum wire. We describe the theoretical scheme using a tight-binding model. The Hamiltonian of the wire is extended via a Lagrange multiplier to "open" the quantum system and to drive current through it. The diagonalization of the extended Hamiltonian yields the transport properties of wire. We show that the Lagrange multiplier method is equivalent to the Landauer approach within the considered model.     

Computing motif correlations in proteins

Protein motifs, which are specific regions and conserved regions, are found by comparing multiple protein sequences. These conserved regions in general play an important role in protein functions and protein folds, for example, for their binding properties or enzymatic activities. The aim here is to find the existence correlations of protein motifs. The knowledge of protein motif/domain sharing should be important in shedding new light on the biologic functions of proteins and offering a basis in analyzing the evolution in the human genome or other genomes. The protein sequences used here are obtained from the PIR-NREF database and the protein motifs are retrieved from the PROSITE database. We apply data mining approach to discover the occurrence correlations of motif in protein sequences. The correlation of motifs mined can be used in evolution analyses and protein structure prediction. We discuss the latter, i.e., protein structure prediction in this study. The correlations mined are stored and maintained in a database system. The database is now available at http://bioinfo.csie.ncu.edu.tw/ProMotif/.