Solution NMR structure of the cold-shock protein from the hyperthermophilic bacterium Thermotoga maritima

Balanced Three-Point Water Model OPC3-B for Intrinsically Disordered and Ordered Proteins

Intrinsically disordered proteins (IDPs) play a critical role in many biological processes. Due to the inherent structural flexibility of IDPs, experimental methods present significant challenges for sampling their conformational information at the atomic level. Therefore, molecular dynamics (MD) simulations have emerged as the primary tools for modeling IDPs whose accuracy depend on force field and water model. To enhance the accuracy of physical modeling of IDPs, several force fields have been developed. However, current water models lack precision and underestimate the interaction between water molecules and proteins. Here, we used Monte-Carlo re-weighting method to re-parameterize a three-point water model based on OPC3 for IDPs (named OPC3-B). We benchmarked the performance of OPC3-B compared with nine different water models for 10 IDPs and three ordered proteins. The results indicate that the performance of OPC3-B is better than other water models for both IDPs and ordered proteins. At the same time, OPC3-B possess the power of transferability with other force field to simulate IDPs. This newly developed water model can be used to insight into the research of sequence-disordered-function paradigm for IDPs.

Thermodynamics of globular protein native structure

The temperature dependence of the partial heat capacity of the native protein structure in an aqueous solution has been analyzed. It is shown that the strictly linear temperature dependence is due to the contributions of the vibrational and conformational components, which indicates volume consistensy and the absence of conformational transitions up to the main two-state transition. The two-level structural and functional organization of the protein three-dimensional structure are considered in relation to the energy and conformational entropy properties in accordance with the principles of the organization of the protein macromolecule.Communicated by Ramaswamy H. Sarma.

NMR Structure and Biophysical Characterization of Thermophilic Single-Stranded DNA Binding Protein from Sulfolobus Solfataricus

Proteins from Sulfolobus solfataricus (S. solfataricus), an extremophile, are active even at high temperatures. The single-stranded DNA (ssDNA) binding protein of S. solfataricus (SsoSSB) is overexpressed to protect ssDNA during DNA metabolism. Although SsoSSB has the potential to be applied in various areas, its structural and ssDNA binding properties at high temperatures have not been studied. We present the solution structure, backbone dynamics, and ssDNA binding properties of SsoSSB at 50 °C. The overall structure is consistent with the structures previously studied at room temperature. However, the loop between the first two β sheets, which is flexible and is expected to undergo conformational change upon ssDNA binding, shows a difference from the ssDNA bound structure. The ssDNA binding ability was maintained at high temperature, but different interactions were observed depending on the temperature. Backbone dynamics at high temperature showed that the rigidity of the structured region was well maintained. The investigation of an N-terminal deletion mutant revealed that it is important for maintaining thermostability, structure, and ssDNA binding ability. The structural and dynamic properties of SsoSSB observed at high temperature can provide information on the behavior of proteins in thermophiles at the molecular level and guide the development of new experimental techniques.

Small Proteins in Archaea, a Mainly Unexplored World

In recent years, increasing numbers of small proteins have moved into the focus of science. Small proteins have been identified and characterized in all three domains of life, but the majority remains functionally uncharacterized, lack secondary structure, and exhibit limited evolutionary conservation. While quite a few have already been described for bacteria and eukaryotic organisms, the amount of known and functionally analyzed archaeal small proteins is still very limited. In this review, we compile the current state of research, show strategies for systematic approaches for global identification of small archaeal proteins, and address selected functionally characterized examples. Besides, we document exemplarily for one archaeon the tool development and optimization to identify small proteins using genome-wide approaches.

Balanced Solvent Model for Intrinsically Disordered and Ordered Proteins

Intrinsically disordered proteins (IDPs) have no fixed three-dimensional (3D) structures under physiological conditions, with the content being about 51% in human proteomics. IDPs are associated with many human diseases, such as cancer, diabetes, and neurodegenerative diseases. Because IDPs do not crystallize and have diverse conformers, traditional experimental methods such as crystallization and NMR can hardly capture their conformation ensemble and just provide average structural characters of IDPs. Therefore, molecular dynamics (MD) simulations become a valuable complement to the experimental data. However, the accuracy of molecular dynamics simulation for IDPs depends on the combination of force fields and solvent models. Recently, we released an environment-specific force field (ESFF1) for IDPs, which can well reproduce the local structural properties (such as J-coupling and secondary chemical shifts). However, there is still a large deviation for the radius of gyration (R_g). Therefore, a solvent model combined with ESFF1 is necessary to capture the local and global characters for IDPs and ordered proteins. Here, we investigated the underestimation or overestimation of the solvent interaction for four solvent models (TIP3P, TIP4P-Ew, TIP4P-D, OPC) under ESFF1 and found the important ε parameter of the solvent model to play a key role in scaling R_g. A near-linear relationship between the simulation R_g and the ε parameter was used to develop the new solvent model, named TIP4P-B. The results indicate that the simulated R_g with TIP4P-B is in better agreement with the experimental observations than the other four solvent models. Simultaneously, TIP4P-B can also maintain the advantages of the ESFF1 force field for the local structural properties. Additionally, TIP4P-B can successfully sample the conformation of ordered proteins. These findings confirm that TIP4P-B is a balanced solvent model and can improve sampling R_g performance for folded proteins and IDPs.

Force Field Benchmark of Amino Acids. 3. Hydration with Scaled Lennard-Jones Interactions

Classical protein force fields were reported with too weak protein-water interactions relative to protein-protein interactions, leading to more compact structures and artificial protein aggregation. Here we investigated the impacts of scaled Lennard-Jones (LJ) interactions on the hydration of amino acids and the simulation of folded and intrinsically disordered proteins (IDPs). The obtained optimal scaling parameters reproduce accurately hydration free energies of neutral amino acid side chain analogues and do not affect the compactness and structural stability of folded proteins significantly. The scaling leads to less compact IDPs and varies from case to case. Strengthening the interactions between protein and water oxygen or hydrogen atoms by increasing the interacting LJ well depth (ε) appears more effective than weakening protein-protein interactions by reducing the interacting dispersion coefficients (C₆). We demonstrate that weakening water-water interactions is a solution as well to obtaining more favorable protein-water interactions in an indirect way, although modern force fields like Amber ff19SB and a99SB-disp tend to use water models with strong water-water interactions. This is likely a compromise between strong protein-protein interactions and strong water-water interactions. Independent optimization of protein force fields and water models is therefore needed to make both interactions more close to reality, leading to good accuracy without bias or scaling.

Cold-Shock Domains-Abundance, Structure, Properties, and Nucleic-Acid Binding

Simple Summary

Proteins are composed of compact domains, often of known three-dimensional structure, and natively unstructured polypeptide regions. The abundant cold-shock domain is among the set of canonical nucleic acid-binding domains and conserved from bacteria to man. Proteins containing cold-shock domains serve a large variety of biological functions, which are mostly linked to DNA or RNA binding. These functions include the regulation of transcription, RNA splicing, translation, stability and sequestration. Cold-shock domains have a simple architecture with a conserved surface ideally suited to bind single-stranded nucleic acids. Because the binding is mostly by non-specific molecular interactions which do not involve the sugar-phosphate backbone, cold-shock domains are not strictly sequence-specific and do not discriminate reliably between DNA and RNA. Many, but not all functions of cold shock-domain proteins in health and disease can be understood based of the physical and structural properties of their cold-shock domains.

Abstract

The cold-shock domain has a deceptively simple architecture but supports a complex biology. It is conserved from bacteria to man and has representatives in all kingdoms of life. Bacterial cold-shock proteins consist of a single cold-shock domain and some, but not all are induced by cold shock. Cold-shock domains in human proteins are often associated with natively unfolded protein segments and more rarely with other folded domains. Cold-shock proteins and domains share a five-stranded all-antiparallel β-barrel structure and a conserved surface that binds single-stranded nucleic acids, predominantly by stacking interactions between nucleobases and aromatic protein sidechains. This conserved binding mode explains the cold-shock domains’ ability to associate with both DNA and RNA strands and their limited sequence selectivity. The promiscuous DNA and RNA binding provides a rationale for the ability of cold-shock domain-containing proteins to function in transcription regulation and DNA-damage repair as well as in regulating splicing, translation, mRNA stability and RNA sequestration.

A suite of 19F based relaxation dispersion experiments to assess biomolecular motions

Proteins and nucleic acids are highly dynamic bio-molecules that can populate a variety of conformational states. NMR relaxation dispersion (RD) methods are uniquely suited to quantify the associated kinetic and thermodynamic parameters. Here, we present a consistent suite of ¹⁹F-based CPMG, on-resonance R_1ρ and off-resonance R_1ρ RD experiments. We validate these experiments by studying the unfolding transition of a 7.5 kDa cold shock protein. Furthermore we show that the ¹⁹F RD experiments are applicable to very large molecular machines by quantifying dynamics in the 360 kDa half-proteasome. Our approach significantly extends the timescale of chemical exchange that can be studied with ¹⁹F RD, adds robustness to the extraction of exchange parameters and can determine the absolute chemical shifts of excited states. Importantly, due to the simplicity of ¹⁹F NMR spectra, it is possible to record complete datasets within hours on samples that are of very low costs. This makes the presented experiments ideally suited to complement static structural information from cryo-EM and X-ray crystallography with insights into functionally relevant motions.

Pleiotropic roles of cold shock proteins with special emphasis on unexplored cold shock protein member of Plasmodium falciparum

The cold shock domain (CSD) forms the hallmark of the cold shock protein family that provides the characteristic feature of binding with nucleic acids. While much of the information is available on bacterial, plants and human cold shock proteins, their existence and functions in the malaria parasite remains undefined. In the present review, the available information on functions of well-characterized cold shock protein members in different organisms has been collected and an attempt was made to identify the presence and role of cold shock proteins in malaria parasite. A single Plasmodium falciparum cold shock protein (PfCoSP) was found in P. falciparum which is reported to be essential for parasite survival. Essentiality of PfCoSP underscores its importance in malaria parasite life cycle. In silico tools were used to predict the features of PfCoSP and to identify its homologues in bacteria, plants, humans, and other Plasmodium species. Modelled structures of PfCoSP and its homologues in Plasmodium species were compared with human cold shock protein 'YBOX-1' (Y-box binding protein 1) that provide important insights into their functioning. PfCoSP model was subjected to docking with B-form DNA and RNA to reveal a number of residues crucial for their interaction. Transcriptome analysis and motifs identified in PfCoSP implicate its role in controlling gene expression at gametocyte, ookinete and asexual blood stages of malaria parasite. Overall, this review emphasizes the functional diversity of the cold shock protein family by discussing their known roles in gene expression regulation, cold acclimation, developmental processes like flowering transition, and flower and seed development, and probable function in gametocytogenesis in case of malaria parasite. This enables readers to view the cold shock protein family comprehensively.

RNA and DNA Binding Epitopes of the Cold Shock Protein TmCsp from the Hyperthermophile Thermotoga maritima

Prokaryotic cold shock proteins (CSPs) are considered to play an important role in the transcriptional and translational regulation of gene expression, possibly by acting as transcription anti-terminators and “RNA chaperones”. They bind with high affinity to single-stranded nucleic acids. Here we report the binding epitope of TmCsp from Thermotoga maritima for both single-stranded DNA and RNA, using heteronuclear 2D NMR spectroscopy. At “physiological” growth temperatures of TmCsp (≥ 343 K), all oligonucleotides studied have dissociation constants between 1.6 ((dT)₇) and 25.2 ((dA)₇) μM as determined by tryptophan fluorescence quenching. Reduction of the temperature to 303 K leads to a pronounced increase of affinity for thymidylate (dT)₇ and uridylate (rU)₇ heptamers with dissociation constants of 4.0 and 10.8 nM, respectively, whereas the weak binding of TmCsp to cytidylate, adenylate, and guanylate heptamers (dC)₇, (dA)₇, and (dT)₇ is almost unaffected by temperature. The change of affinities of TmCsp for (dT)₇ and (rU)₇ by approximately 3 orders of magnitude shows that it represents a cold chock sensor that switches on the cold shock reaction of the cell. A temperature dependent conformational switch of the protein is required for this action. The binding epitope on TmCsp for the ssDNA and RNA heptamers is very similar and comprises β-strands 1 and 2, the loop β1–β2 as well as the loops connecting β3 with β4 and β4 with β5. Besides the loop regions, surprisingly, mainly the RNA-binding motif RNP1 is involved in ssDNA and RNA binding, while only two amino acids, H28 and W29, of the postulated RNA-binding motif RNP2 interact with the uridylate and thymidylate homonucleotides, although a high affinity in the nanomolar range is achieved. This is in contrast to the binding properties of other CSPs or cold shock domains, where RNP1 as well as RNP2 are involved in binding. TmCsp takes up a unique position since it is the only one which possesses a tryptophan residue instead of a usually highly conserved phenylalanine or tyrosine residue at the end of RNP2. NMR titrations suggest that neither (dT)₇ nor (rU)₇ represent the full binding motif and that non-optimal intercalation of W29 into these oligonucleotides blocks the access of the RNP2 site to the DNA or RNA. NMR-experiments with (dA)₇ suggest an interaction of W29 with the adenine ring. Full binding seems to require at least one single purine base well-positioned within a thymine- or uracil-rich stretch of nucleic acids.

Comparison and Evaluation of Force Fields for Intrinsically Disordered Proteins

Molecular dynamics (MD) simulations of six upgraded empirical force fields were compared and evaluated with short peptides, intrinsically disordered proteins, and folded proteins using trajectories of 1, 1.5, 5, or 10 μs (five replicates of 200 ns, 300 ns, 1 μs, or 2 μs) for each system. Previous studies have shown that different force fields, water models, simulation methods, and parameters can affect simulation outcomes. Here, the MD simulations were done in an explicit solvent with RS-peptide, HEWL19, HIV-rev, β amyloid (Aβ)-40, Aβ-42, phosphodiesterase-γ, CspTm, and ubiquitin using ff99IDPs, ff14IDPs, ff14IDPSFF, ff03w, CHARMM36m, and CHARMM22* force fields. The IDP ensembles generated by six all-atom empirical force fields were compared against NMR data. Despite using identical starting structures and simulation parameters, ensembles obtained with different force fields exhibit significant differences in NMR RMDs, secondary structure contents, and global properties such as the radius of gyration. The intrinsically disordered protein (IDP)-specific force fields could substantially reproduce the experimental observables in force field comparison: they have the lowest error in chemical shifts and J-couplings for short peptides/proteins, reasonably well for large IDPs and reasonably well with the radius of gyration. A high population of disorderness was observed in the IDP-specific force field for the IDP ensemble with a fraction of β sheets for β-amyloids. CHARMM22* performs better for many observables; however, it still has a preference toward the helicity for short peptides. The results of β-amyloid 42 starting from two different initial structures (Aβ42^1Z0Q and Aβ42^model) were also compared with DSSP and NMR data. The results obtained with IDP-specific force fields within 2 μs simulation time are similar, even though starting from different structures. The current force fields perform equally well for folded proteins. The results of currently developed or modified force fields for IDPs are capable of enlightening the overall performance of the force field for disordered as well as folded proteins, thereby contributing to force field development.

The Role of Electrostatics and Folding Kinetics on the Thermostability of Homologous Cold Shock Proteins

Understanding which aspects contribute to the thermostability of proteins is a challenge that has persisted for decades, and it is of great relevance for protein engineering. Several types of interactions can influence the thermostability of a protein. Among them, the electrostatic interactions have been a target of particular attention. Aiming to explore how this type of interaction can affect protein thermostability, this paper investigated four homologous cold shock proteins from psychrophilic, mesophilic, thermophilic, and hyperthermophilic organisms using a set of theoretical methodologies. It is well-known that electrostatics as well as hydrophobicity are key-elements for the stabilization of these proteins. Therefore, both interactions were initially analyzed in the native structure of each protein. Electrostatic interactions present in the native structures were calculated with the Tanford-Kirkwood model with solvent accessibility, and the amount of hydrophobic surface area buried upon folding was estimated by measuring both folded and extended structures. On the basis of Energy Landscape Theory, the local frustration and the simplified alpha-carbon structure-based model were modeled with a Debye-Hückel potential to take into account the electrostatics and the effects of an implicit solvent. Thermodynamic data for the structure-based model simulations were collected and analyzed using the Weighted Histogram Analysis and Stochastic Diffusion methods. Kinetic quantities including folding times, transition path times, folding routes, and Φ values were also obtained. As a result, we found that the methods are able to qualitatively infer that electrostatic interactions play an important role on the stabilization of the most stable thermophilic cold shock proteins, showing agreement with the experimental data.

Tyr51: Key Determinant of the Low Thermostability of the Colwellia psychrerythraea Cold-Shock Protein

Cold-shock proteins (Csps) are expressed at lower-than-optimum temperatures, and they function as RNA chaperones; however, no structural studies on psychrophilic Csps have been reported. Here, we aimed to investigate the structure and dynamics of the Csp of psychrophile Colwellia psychrerythraea 34H, ( Cp-Csp). Although Cp-Csp shares sequence homology, common folding patterns, and motifs, including a five β-stranded barrel, with its thermophilic counterparts, its thermostability (37 °C) was markedly lower than those of other Csps. Cp-Csp binds heptathymidine with an affinity of 10^-7 M, thereby increasing its thermostability to 50 °C. Nuclear magnetic resonance spectroscopic analysis of the Cp-Csp structure and backbone dynamics revealed a flexible structure with only one salt bridge and 10 residues in the hydrophobic cavity. Notably, Cp-Csp contains Tyr51 instead of the conserved Phe in the hydrophobic core, and its phenolic hydroxyl group projects toward the surface. The Y51F mutation increased the stability of hydrophobic packing and may have allowed for the formation of a K3-E21 salt bridge, thereby increasing its thermostability to 43 °C. Cp-Csp exhibited conformational exchanges in its ribonucleoprotein motifs 1 and 2 (754 and 642 s^-1), and heptathymidine binding markedly decreased these motions. Cp-Csp lacks salt bridges and has longer flexible loops and a less compact hydrophobic cavity resulting from Tyr51 compared to mesophilic and thermophilic Csps. These might explain the low thermostability of Cp-Csp. The conformational flexibility of Cp-Csp facilitates its accommodation of nucleic acids at low temperatures in polar oceans and its function as an RNA chaperone for cold adaptation.

Structure and interaction of Corynebacterium pseudotuberculosis cold shock protein A with Y-box single-stranded DNA fragment

Cold shock proteins (Csps) function to preserve cell viability at low temperatures by binding to nucleic acids and consequently control gene expression. The mesophilic bacterium Corynebacterium pseudotuberculosis is the causative agent of caseous lymphadenitis in animals, and infection in livestock is a considerable economic burden worldwide. In this report, the structure of cold shock protein A from Cp (Cp-CspA) and biochemical analysis of its temperature-dependent interaction with a Y-box ssDNA motif is presented. The Cp-CspA structure contains five β-strands making up a β-barrel fold with 11 hydrophobic core residues and two salt bridges that confers it with a melting temperature of ~ 54 °C that is similar to mesophilic Bs-CspB. Chemical shift perturbations analysis revealed that residues in the nucleic acid-binding motifs (RNP 1 and 2) and loop 3 are involved in binding to the Y-box fragment either by direct interaction or by conformational rearrangements remote from the binding region. Fluorescence quenching experiments of Cp-CspA showed that the dissociation constants for Y-box ssDNA binding is nanomolar and the binding affinity decreased as the temperature increased, indicating that the interaction is enthalpically driven and the hydrogen bonds and van der Waals forces are important contributions for complex stabilization. The Y31 of Cp-CspA is a particular occurrence among Csps from mesophilic bacteria that provide a possible explanation for the higher binding affinity to ssDNA than that observed for Bs-CspB. Anisotropy measurements indicated that the reduction in molecular mobility of Cp-CspA upon Y-box binding is characterized by a cooperative process.

DATABASE

Resonance assignment and structural data are available in the Biological Magnetic Resonance Data Bank and Protein Data Bank under accession number 26802 and 5O6F, respectively.

Mechanics of Protein Adaptation to High Temperatures

Inspired by Somero's corresponding state principle that relates protein enhanced thermal stability with mechanical rigidity, we deployed state of the art computational techniques (based on atomistic steered molecular dynamics and Hamiltonian-replica exchange simulations) to study the in silico realization of mechanical and thermal unfolding of two homologous Csp proteins that have evolved to thrive in different thermal environments. By complementing recent single-molecule experiments, we unambiguously show that, for these homologues whose structures are very similar, the increased thermal resistance of the thermophilic variant is not associated with an increased mechanical stability. Our approach provides microscopic insights that are otherwise inaccessible to experimental techniques, and explains why the protein weak spots for thermal and mechanical denaturation are distinct.

RNA binding and chaperone activity of the E. coli cold-shock protein CspA

Ensuring the correct folding of RNA molecules in the cell is of major importance for a large variety of biological functions. Therefore, chaperone proteins that assist RNA in adopting their functionally active states are abundant in all living organisms. An important feature of RNA chaperone proteins is that they do not require an external energy source to perform their activity, and that they interact transiently and non-specifically with their RNA targets. So far, little is known about the mechanistic details of the RNA chaperone activity of these proteins. Prominent examples of RNA chaperones are bacterial cold shock proteins (Csp) that have been reported to bind single-stranded RNA and DNA. Here, we have used advanced NMR spectroscopy techniques to investigate at atomic resolution the RNA-melting activity of CspA, the major cold shock protein of Escherichia coli, upon binding to different RNA hairpins. Real-time NMR provides detailed information on the folding kinetics and folding pathways. Finally, comparison of wild-type CspA with single-point mutants and small peptides yields insights into the complementary roles of aromatic and positively charged amino-acid side chains for the RNA chaperone activity of the protein.

Unusual dimerization of a BcCsp mutant leads to reduced conformational dynamics

Cold shock proteins (Csp) constitute a family of ubiquitous small proteins that act as RNA-chaperones to avoid cold-induced termination of translation. All members contain two subdomains composed of 2 and 3 β-strands, respectively, which are connected by a hinge loop and fold into a β-barrel. Bacillus caldolyticus Csp (BcCsp) is one of the most studied members of the family in terms of its folding, function, and structure. This protein has been described as a monomer in solution, although a recent crystal structure showed dimerization via domain swapping (DS). In contrast, other cold shock proteins of the same fold are known to dimerize in a nonswapped arrangement. Hypothesizing that reducing the size of the hinge loop may promote swapping as in several other DS proteins with different folds we deleted two residues from these region (BcCsp∆36-37), leading to a protein in monomer-dimer equilibrium with similar folding stability to that of the wild-type. Strikingly, the crystal structure of BcCsp∆36-37 revealed a nonswapped dimer with its interface located at the nucleic acid-binding surface, showing that the deletion led to structural consequences far from the perturbation site. Concomitantly, circular dichroism experiments on BcCsp∆36-37 demonstrated that binding of the oligonucleotide hexathymidine disrupts the dimer. Additionally, HDXMS shows a protective effect on the protein structure upon dimerization, where the resulting interactions between ligand-binding surfaces in the dimer reduced the extent of exchange throughout the whole protein. Our work provides evidence of the complex interplay between conformational dynamics, deletions, and oligomerization within the Csp protein family.

DATABASES

Structural data are available in the Protein Data Bank under accession number 5JX4.

Terahertz Frequency Spectroscopy to Determine Cold Shock Protein Stability upon Solvation and Evaporation - A Molecular Dynamics Study

Infrared (IR) and Terahertz (THz) spectroscopy simulations were carried out using CHARMM35b2 to determine protein stability. The stabilities of three bacterial cold shock proteins (Csps) originating from mesophiles, thermophiles and hyper- thermophiles respectively were investigated in this study. The three different Csps were investigated by Normal-Mode analysis and Molecular Dynamics simulation of THz spectra using the Hessian matrix for solvated systems, interpreted in the harmonic approximation at optimum near-melting temperatures of each homologue, by incorporating differences in the hydrous and anhydrous states of the Csps. The results show slight variations in the large scale protein motion. However, the IR spectra of Csps observed at the low frequency saddle surface region, clearly distinguishes the thermophilic and mesophilic proteins based on their stability. Further studies on protein stability employing low-frequency collective modes have the potential to reveal functionally important conformational changes that are biologically significant.

Reconciling Intermediates in Mechanical Unfolding Experiments with Two-State Protein Folding in Bulk

Most experimentally well-characterized single domain proteins of less than 100 residues have been found to be two-state folders. That is, only two distinct populations can explain both equilibrium and kinetic measurements. Results from single molecule force spectroscopy, where a protein is unfolded by applying a mechanical pulling force to its ends, have largely confirmed this description for proteins found to be two-state in ensemble experiments. Recently, however, stable intermediates have been reported in mechanical unfolding experiments on a cold-shock protein previously found to be a prototypical two-state folder. Here, we tackle this discrepancy using free energy landscapes and Markov state models derived from coarse-grained molecular simulations. We show that protein folding intermediates can be selectively stabilized by the pulling force and that the populations of these intermediates vary in a force-dependent manner. Our model qualitatively captures the experimental results and suggests a possible origin of the apparent discrepancy.

The physics of pulling polyproteins: a review of single molecule force spectroscopy using the AFM to study protein unfolding

One of the most exciting developments in the field of biological physics in recent years is the ability to manipulate single molecules and probe their properties and function. Since its emergence over two decades ago, single molecule force spectroscopy has become a powerful tool to explore the response of biological molecules, including proteins, DNA, RNA and their complexes, to the application of an applied force. The force versus extension response of molecules can provide valuable insight into its mechanical stability, as well as details of the underlying energy landscape. In this review we will introduce the technique of single molecule force spectroscopy using the atomic force microscope (AFM), with particular focus on its application to study proteins. We will review the models which have been developed and employed to extract information from single molecule force spectroscopy experiments. Finally, we will end with a discussion of future directions in this field.

Modest influence of FRET chromophores on the properties of unfolded proteins

Single-molecule Förster resonance energy transfer (FRET) experiments are often used to study the properties of unfolded and intrinsically disordered proteins. Because of their large extinction coefficients and quantum yields, synthetic heteroaromatic chromophores covalently linked to the protein are often used as donor and acceptor fluorophores. A key issue in the interpretation of such experiments is the extent to which the properties of the unfolded chain may be affected by the presence of these chromophores. In this article, we investigate this question using all-atom explicit solvent replica exchange molecular dynamics simulations of three different unfolded or intrinsically disordered proteins. We find that the secondary structure and long-range contacts are largely the same in the presence or absence of the fluorophores, and that the dimensions of the chain with and without chromophores are similar. This suggests that, at least in the cases studied, extrinsic fluorophores have little effect on the structural properties of unfolded or disordered proteins. We also find that the critical FRET orientational factor κ(2), has an average value and equilibrium distribution very close to that expected for isotropic orientations, which supports one of the assumptions frequently made when interpreting FRET efficiency in terms of distances.

Accounting for dye diffusion and orientation when relating FRET measurements to distances: three simple computational methods

Förster resonance energy transfer (FRET) allows in principal for the structural changes of biological systems to be revealed by monitoring distributions and distance fluctuations between parts of individual molecules. However, because flexible probes usually have to be attached to the macromolecule to conduct these experiments, they suffer from uncertainty in probe positions and orientations. One of the way to address this issue is to use molecular dynamics simulations to explicitly model the likely positions of the probes, but, this is still not widely accessible because of the large computational effort required. Here we compare three simpler methods that can potentially replace MD simulations in FRET data interpretation. In the first, the volume accessible for dye movement is calculated using a fast, geometrical algorithm. The next method, adapted from the analysis of electron paramagnetic studies, utilises a library of rotamers describing probe conformations. The last method uses preliminary MD simulations of fluorescent dyes in solution, to identify all conformational states of dyes and overlays this on the macromolecular system. A comparison of these methods in the simple system of dye-labelled polyproline, shows that in the case of lack of interaction between the dye and host, all give results comparable with MD simulations but require much less time. Differences between these three methods and their ability to compete with MD simulations in the analysis of real experiment are demonstrated and discussed using the examples of cold shock protein and leucine transporter systems.

Water dispersion interactions strongly influence simulated structural properties of disordered protein states

Many proteins can be partially or completely disordered under physiological conditions. Structural characterization of these disordered states using experimental methods can be challenging, since they are composed of a structurally heterogeneous ensemble of conformations rather than a single dominant conformation. Molecular dynamics (MD) simulations should in principle provide an ideal tool for elucidating the composition and behavior of disordered states at an atomic level of detail. Unfortunately, MD simulations using current physics-based models tend to produce disordered-state ensembles that are structurally too compact relative to experiments. We find that the water models typically used in MD simulations significantly underestimate London dispersion interactions, and speculate that this may be a possible reason for these erroneous results. To test this hypothesis, we create a new water model, TIP4P-D, that approximately corrects for these deficiencies in modeling water dispersion interactions while maintaining compatibility with existing physics-based models. We show that simulations of solvated proteins using this new water model typically result in disordered states that are substantially more expanded and in better agreement with experiment. These results represent a significant step toward extending the range of applicability of MD simulations to include the study of (partially or fully) disordered protein states.

Life in extreme environments: single molecule force spectroscopy as a tool to explore proteins from extremophilic organisms

Extremophiles are organisms which survive and thrive in extreme environments. The proteins from extremophilic single-celled organisms have received considerable attention as they are structurally stable and functionally active under extreme physical and chemical conditions. In this short article, we provide an introduction to extremophiles, the structural adaptations of proteins from extremophilic organisms and the exploitation of these proteins in industrial applications. We provide a review of recent developments which have utilized single molecule force spectroscopy to mechanically manipulate proteins from extremophilic organisms and the information which has been gained about their stability, flexibility and underlying energy landscapes.

Structure and flexibility of the thermophilic cold-shock protein of Thermus aquaticus

The thermophilic bacterium Thermus aquaticus is a well-known source of Taq polymerase. Here, we studied the structure and dynamics of the T. aquaticus cold-shock protein (Ta-Csp) to better understand its thermostability using NMR spectroscopy. We found that Ta-Csp has a five-stranded β-barrel structure with five salt bridges which are important for more rigid structure and a higher melting temperature (76 °C) of Ta-Csp compared to mesophilic and psychrophilic Csps. Microsecond to millisecond time scale exchange processes occur only at the β1-β2 surface region of the nucleic acid binding site with an average conformational exchange rate constant of 674 s(-1). The results imply that thermophilic Ta-Csp has a more rigid structure and may not need high structural flexibility to accommodate nucleic acids upon cold shock compared to its mesophile and psychrophile counterparts.

Structural and dynamic features of cold-shock proteins of Listeria monocytogenes, a psychrophilic bacterium

Cold-shock proteins (Csps), proteins expressed when the ambient temperature drops below the growth-supporting temperature, bind to single-stranded nucleic acids and act as RNA chaperones to regulate translation. Listeria monocytogenes is a psychrophilic food-borne pathogen that is problematic for the food industry. Structures of Csps from psychrophilic bacteria have not yet been studied. Despite dramatic differences in the thermostability of Csps of various thermophilic microorganisms, these proteins share a high degree of primary sequence homology and a high degree of three-dimensional structural similarity. Here, we investigated the structural and dynamic features as well as the thermostability of L. monocytogenes CspA (Lm-CspA). Lm-CspA has a five-stranded β-barrel structure with hydrophobic core packing and two salt bridges. When heptathymidine (dT(7)) binds, values for the heteronuclear nuclear Overhauser effect and order parameters of residues in surface loop regions near nucleic acid binding sites increase dramatically. Moreover, Carr-Purcell-Meiboom-Gill experiments showed that slow motions observed for the nucleic acid binding residues K7, W8, F15, F27, and R56 disappeared in Lm-CspA-dT(7). Lm-CspA is less thermostable than mesophilic and thermophilic Csps, with a lower melting temperature (40 °C). The structural flexibility that accompanies longer surface loops and less hydrophobic core packing and a number of salt bridges and unfavorable electrostatic repulsion are likely key factors in the low thermostability of Lm-CspA. This implies that the large conformational flexibility of psychrophilic Lm-CspA, which more easily accommodates nucleic acids at low temperature, is required for RNA chaperone function under cold-shock conditions and for the cold adaptation of L. monocytogenes.

Single-molecule force spectroscopy identifies a small cold shock protein as being mechanically robust

Single-molecule force spectroscopy has emerged as a powerful approach to examine the stability and dynamics of single proteins. We have completed force extension experiments on the small cold shock protein B from Thermotoga maritima, using a specially constructed chimeric polyprotein. The protein's simple topology, which is distinct from the mechanically well-characterized β-grasp and immunoglobulin (Ig)-like folds, in addition to the wide range of structural homologues resulting from its ancient origin, provides an attractive model protein for single-molecule force spectroscopy studies. We have determined that the protein has mechanical stability, unfolding at greater than 70 pN at a pulling velocity of 100 nm s(-1). We reveal features of the unfolding energy landscape by measuring the dependence of the mechanical stability on pulling velocity, in combination with Monte Carlo simulations. We show that the cold shock protein has mechanically robust, yet malleable, features that may be important in providing the protein with stability and flexibility to function over a range of environmental conditions. These results provide insights into the relationship between the secondary structure and topology of a protein and its mechanical strength. This lays the foundation for the investigation of the effects of changes in environmental conditions on the mechanical and dynamic properties of cold shock proteins.

Solution structure of the cold-shock-like protein from Rickettsia rickettsii

The solution structure of the cold-shock-like protein from R. rickettsii, the causative agent of Rocky Mountain spotted fever, is reported.

Rocky Mountain spotted fever is caused by Rickettsia rickettsii infection. R. rickettsii can be transmitted to mammals, including humans, through the bite of an infected hard-bodied tick of the family Ixodidae. Since the R. rickettsii genome contains only one cold-shock-like protein and given the essential nature of cold-shock proteins in other bacteria, the structure of the cold-shock-like protein from R. rickettsii was investigated. With the exception of a short α-helix found between β-strands 3 and 4, the solution structure of the R. rickettsii cold-shock-like protein has the typical Greek-key five-stranded β-barrel structure found in most cold-shock domains. Additionally, the R. rickettsii cold-shock-like protein, with a ΔG of unfolding of 18.4 kJ mol⁻¹, has a similar stability when compared with other bacterial cold-shock proteins.

Configuration-dependent diffusion dynamics of downhill and two-state protein folding

Configuration-dependent diffusion (CDD) is important for protein folding kinetics with small thermodynamic barriers. CDD can be even more crucial in downhill folding without thermodynamic barriers. We explored the CDD of a downhill protein (BBL), and a two-state protein (CI2). The hidden kinetic barriers due to CDD were revealed. The increased ~1 k(B)T kinetic barrier is in line with experimental value based on other fast folding proteins. Compared to that of CI2, the effective free-energy profile of BBL is found to be significantly influenced by CDD, and the kinetics are totally determined by diffusion. These findings are consistent with both earlier bulk and single-molecule fluorescence measurements. In addition, we found the temperature dependence of CDD. We also found that the ratio of folding transition temperature against optimal kinetic folding temperature can provide both a quantitative measure for the underlying landscape topography and an indicator for the possible appearance of downhill folding. Our study can help for a better understanding of the role of diffusion in protein folding dynamics.

A non-cold-inducible cold shock protein homolog mainly contributes to translational control under optimal growth conditions

Cold shock proteins (Csps) include both cold-induced and non-cold-induced proteins, contrary to their name. Cold-induced Csps are well studied; they function in cold acclimation by controlling transcription and translation. Some Csps have been reported to contribute to other cellular processes. However, the functions of non-cold-induced Csps under optimal growth conditions remain unknown. To elucidate these functions, we used transcriptome and proteome analyses as comprehensive approaches and have compared the outputs of wild-type and non-cold-induced Csp-deletion mutant cells. As a model organism, we selected Thermus thermophilus HB8 because it has only two csp genes (ttcsp1 and ttcsp2); ttCsp1 is the only non-cold-induced Csp. Surprisingly, the amount of transcripts and proteins upon deletion of the ttcsp1 gene was quite different. DNA microarray analysis revealed that the deletion of ttcsp1 did not affect the amount of transcripts, although the ttcsp1 gene was constantly expressed in the wild-type cell. Nonetheless, proteomic analysis revealed that the expression levels of many proteins were significantly altered when ttcsp1 was deleted. These results suggest that ttCsp1 functions in translation independent of transcription. Furthermore, ttCsp1 is involved in both the stimulation and inhibition of translation of specific proteins. Here, we have determined the crystal structure of ttCsp1 at 1.65 Å. This is the first report to present the structure of a non-cold-inducible cold shock protein. We also report the nucleotide binding affinity of ttCsp1. Finally, we discuss the functions of non-cold-induced Csps and propose how they modulate the levels of specific proteins to suit the prevailing environmental conditions.

RNA single strands bind to a conserved surface of the major cold shock protein in crystals and solution

Bacterial cold shock proteins (CSPs) regulate the cellular response to temperature downshift. Their general principle of function involves RNA chaperoning and transcriptional antitermination. Here we present two crystal structures of cold shock protein B from Bacillus subtilis (Bs-CspB) in complex with either a hexanucleotide (5'-UUUUUU-3') or heptanucleotide (5'-GUCUUUA-3') single-stranded RNA (ssRNA). Hydrogen bonds and stacking interactions between RNA bases and aromatic sidechains characterize individual binding subsites. Additional binding subsites which are not occupied by the ligand in the crystal structure were revealed by NMR spectroscopy in solution on Bs-CspB·RNA complexes. Binding studies demonstrate that Bs-CspB associates with ssDNA as well as ssRNA with moderate sequence specificity. Varying affinities of oligonucleotides are reflected mainly in changes of the dissociation rates. The generally lower binding affinity of ssRNA compared to its ssDNA analog is attributed solely to the substitution of thymine by uracil bases in RNA.

Two-intermediate model to characterize the structure of fast-folding proteins

This paper introduces a new model that enables researchers to conduct protein folding simulations. A two-step in silico process is used in the course of structural analysis of a set of fast-folding proteins. The model assumes an early stage (ES) that depends solely on the backbone conformation, as described by its geometrical properties--specifically, by the V-angle between two sequential peptide bond planes (which determines the radius of curvature, also called R-radius, according to a second-degree polynomial form). The agreement between the structure under consideration and the assumed model is measured in terms of the magnitude of dispersion of both parameters with respect to idealized values. The second step, called late-stage folding (LS), is based on the "fuzzy oil drop" model, which involves an external hydrophobic force field described by a three-dimensional Gauss function. The degree of conformance between the structure under consideration and its idealized model is expressed quantitatively by means of the Kullback-Leibler entropy, which is a measure of disparity between the observed and expected hydrophobicity distributions. A set of proteins, representative of the fast-folding group - specifically, cold shock proteins - is shown to agree with the proposed model.

Structure and function of a cold shock domain fold protein, CspD, in Janthinobacterium sp. Ant5-2 from East Antarctica

A cold shock domain (CSD)-containing protein, CspD, of molecular mass ~7.28 kDa in a psychrotolerant Antarctic Janthinobacterium sp. Ant5-2 (ATCC BAA-2154) exhibited constitutive expression at 37, 22, 15, 4 and -1°C. The cspD gene encoding the CspD protein of Ant5-2 was cloned, sequenced and analyzed. The deduced protein sequence was highly similar to the conserved domains of the cold shock proteins (Csps) from bacteria belonging to the class Betaproteobacteria. Its expression was both time- and growth phase-dependent and increased when exposed to 37°C and UV radiation (UVC, dose: 1.8 and 2.8 mJ cm(-2)). The results from the electrophoretic mobility shift and subcellular localization study confirmed its single-stranded DNA-binding property. In silico analysis of the deduced tertiary structure of CspD from Ant5-2 showed a highly stable domain-swapped dimer, forming two similar monomeric Csp folds. This study established an overall framework of the structure, function and phylogenetic analysis of CspD from an Antarctic Janthinobacterium sp. Ant5-2, which may facilitate and stimulate the study of CSD fold proteins in the class Betaproteobacteria.

The origin of nonmonotonic complex behavior and the effects of nonnative interactions on the diffusive properties of protein folding

We present a method for calculating the configurational-dependent diffusion coefficient of a globular protein as a function of the global folding process. Using a coarse-grained structure-based model, we determined the diffusion coefficient, in reaction coordinate space, as a function of the fraction of native contacts formed Q for the cold shock protein (TmCSP). We find nonmonotonic behavior for the diffusion coefficient, with high values for the folded and unfolded ensembles and a lower range of values in the transition state ensemble. We also characterized the folding landscape associated with an energetically frustrated variant of the model. We find that a low-level of frustration can actually stabilize the native ensemble and increase the associated diffusion coefficient. These findings can be understood from a mechanistic standpoint, in that the transition state ensemble has a more homogeneous structural content when frustration is present. Additionally, these findings are consistent with earlier calculations based on lattice models of protein folding and more recent single-molecule fluorescence measurements.

Coordinate and time-dependent diffusion dynamics in protein folding

We developed both analytical and simulation methods to explore the diffusion dynamics in protein folding. We found the diffusion as a quantitative measure of escape from local traps along the protein folding funnel with chosen reaction coordinates has two remarkable effects on kinetics. At a fixed coordinate, local escape time depends on the distribution of barriers around it, therefore the diffusion is often time distributed. On the other hand, the environments (local escape barriers) change along the coordinates, therefore diffusion is coordinate dependent. The effects of time-dependent diffusion on folding can lead to non-exponential kinetics and non-Poisson statistics of folding time distribution. The effects of coordinate dependent diffusion on folding can lead to the change of the kinetic barrier height as well as the position of the corresponding transition state and therefore modify the folding kinetic rates as well as the kinetic routes. Our analytical models for folding are based on a generalized Fokker-Planck diffusion equation with diffusion coefficient both dependent on coordinate and time. Our simulation for folding are based on structure-based folding models with a specific fast folding protein CspTm studied experimentally on diffusion and folding with single molecules. The coordinate and time-dependent diffusion are especially important to be considered in fast folding and single molecule studies, when there is a small or no free energy barrier and kinetics is controlled by diffusion while underlying statistics of kinetics become important. Including the coordinate dependence of diffusion will challenge the transition state theory of protein folding. The classical transition state theory will have to be modified to be consistent. The more detailed folding mechanistic studies involving phi value analysis based on the classical transition state theory will also have to be quantitatively modified. Complex kinetics with multiple time scales may allow us not only to explore the folding kinetics but also probe the local landscape and barrier height distribution with single-molecule experiments.

Integrase-directed recovery of functional genes from genomic libraries

Large population sizes, rapid growth and 3.8 billion years of evolution firmly establish microorganisms as a major source of the planet's biological and genetic diversity. However, up to 99% of the microorganisms in a given environment cannot be cultured. Culture-independent methods that directly access the genetic potential of an environmental sample can unveil new proteins with diverse functions, but the sequencing of random DNA can generate enormous amounts of extraneous data. Integrons are recombination systems that accumulate open reading frames (gene cassettes), many of which code for functional proteins with enormous adaptive potential. Some integrons harbor hundreds of gene cassettes and evidence suggests that the gene cassette pool may be limitless in size. Accessing this genetic pool has been hampered since sequence-based techniques, such as hybridization or PCR, often recover only partial genes or a small subset of those present in the sample. Here, a three-plasmid genetic strategy for the sequence-independent recovery of gene cassettes from genomic libraries is described and its use by retrieving functional gene cassettes from the chromosomal integron of Vibrio vulnificus ATCC 27562 is demonstrated. By manipulating the natural activity of integrons, we can gain access to the caches of functional genes amassed by these structures.

Structure of the cold-shock domain protein from Neisseria meningitidis reveals a strand-exchanged dimer

The structure of the cold-shock domain protein from Neisseria meningitidis has been solved to 2.6 A resolution and shown to comprise a dimer formed by the exchange of two beta-strands between protein monomers. The overall fold of the monomer closely resembles those of other bacterial cold-shock proteins. The neisserial protein behaved as a monomer in solution and was shown to bind to a hexathymidine oligonucleotide with a stoichiometry of 1:1 and a K(d) of 1.25 microM.

Characterization of cold-shock protein A of Antarctic Streptomyces sp. AA8321

Polar organisms should have mechanisms to survive the extremely cold environment. Four genes encoding cold-shock proteins, which are small, cold-induced bacterial proteins, have been cloned from the Antarctic bacterium Streptomyces sp. AA8321. Since the specific functions of any polar bacterial or Streptomyces cold-shock proteins have not yet been determined, we examined the role of cold-shock protein A from Streptomyces sp. AA8321 (CspA(St)). Gel filtration chromatography showed that purified CspA(St) exists as a homodimer under physiological conditions, and gel shift assays showed that it binds to single-stranded, but not double-stranded, DNA. Overexpression of CspA(St) in Escherichia coli severely impaired the ability of the host cells to form colonies, and the cells developed an elongated morphology. Incorporation of a deoxynucleoside analogue, 5-bromo-2'-deoxyuridine, into newly synthesized DNA was also drastically diminished in CspA(St)-overexpressing cells. These results suggest that CspA(St) play a role in inhibition of DNA replication during cold-adaptation.