Molecular dynamics simulations of N-terminal peptides from a nucleotide binding protein

Molecular dynamics (MD) simulations of N-terminal peptides from lactate dehydrogenase (LDH) with increasing length and individual secondary structure elements were used to study their stability in relation to folding. Ten simulations of 1–2 ns of different peptides in water starting from the coordinates of the crystal structure were performed. The stability of the peptides was compared qualitatively by analyzing the root mean square deviation (RMSD) from the crystal structure, radius of gyration, secondary and tertiary structure, and solvent accessible surface area. In agreement with earlier MD studies, relatively short (< 15 amino acids) peptides containing individual secondary structure elements were generally found to be unstable; the hydrophobic α₁-helix of the nucleotide binding fold displayed a significantly higher stability, however. Our simulations further showed that the first βαβ supersecondary unit of the characteristic dinucleotide binding fold (Rossmann fold) of LDH is somewhat more stable than other units of similar length and that the α₂-helix, which unfolds by itself, is stabilized by binding to this unit. This finding suggests that the first βαβ unit could function as an N-terminal folding nucleus, upon which the remainder of the polypeptide chain can be assembled. Indeed, simulations with longer units (βαβα and βαβαββ) showed that all structural elements of these units are rather stable. The outcome of our studies is in line with suggestions that folding of the N-terminal portion of LDH in vivo can be a cotranslational process that takes place during the ribosomal peptide synthesis.     

Bacterial IscU is a well folded and functional single domain protein.

Iron-sulfur clusters are widely represented in most organisms, but the mechanism of their formation is not fully understood. Of the two main proteins involved in cluster formation, NifS/IscS and NifU/IscU, only the former has been well studied from a structural point of view. Here we report an extensive structural characterization of Escherichia coli IscU. We show by a variety of physico-chemical techniques that E. coli IscU construct can be expressed to high purity as a monomeric protein, characterized by an alphabeta fold with high alpha-helix content. The high melting temperature and the reversibility of the thermal unfolding curve (as measured by CD spectroscopy) hint at a well ordered stable fold. The excellent dispersion of cross peaks in the (1)H-(15)N correlation spectrum is consistent with these observations. Monomeric E. coli IscU is able to provide a scaffold for Iron-sulfur cluster assembly, but has no direct interaction with either Fe(II) or Fe(III) ions, suggesting the need of further partners to achieve a stable interaction.     

The crystal structure of NlpI. A prokaryotic tetratricopeptide repeat protein with a globular fold

There are several different families of repeat proteins. In each, a distinct structural motif is repeated in tandem to generate an elongated structure. The nonglobular, extended structures that result are particularly well suited to present a large surface area and to function as interaction domains. Many repeat proteins have been demonstrated experimentally to fold and function as independent domains. In tetratricopeptide (TPR) repeats, the repeat unit is a helix-turn-helix motif. The majority of TPR motifs occur as three to over 12 tandem repeats in different proteins. The majority of TPR structures in the Protein Data Bank are of isolated domains. Here we present the high-resolution structure of NlpI, the first structure of a complete TPR-containing protein. We show that in this instance the TPR motifs do not fold and function as an independent domain, but are fully integrated into the three-dimensional structure of a globular protein. The NlpI structure is also the first TPR structure from a prokaryote. It is of particular interest because it is a membrane-associated protein, and mutations in it alter septation and virulence.     

Evolution of protein sequences and structures.

The relationship between sequence similarity and structural similarity has been examined in 36 protein families with five or more diverse members whose structures are known. The structural similarity within a family (as determined with the DALI structure comparison program) is linearly related to sequence similarity (as determined by a Smith-Waterman search of the protein sequences in the structure database). The correlation between structural similarity and sequence similarity is very high; 18 of the 36 families had linear correlation coefficients r>/=0.878, and only nine had correlation coefficients r相似文献   

Secondary structure of the 3' terminus of hepatitis C virus minus-strand RNA

The 3'-terminal ends of both the positive and negative strands of the hepatitis C virus (HCV) RNA, the latter being the replicative intermediate, are most likely the initiation sites for replication by the viral RNA-dependent RNA polymerase, NS5B. The structural features of the very conserved 3' plus [(+)] strand untranslated region [3' (+) UTR] are well established (K. J. Blight and C. M. Rice, J. Virol. 71:7345-7352, 1997). However, little information is available concerning the 3' end of the minus [(-)] strand RNA. In the present work, we used chemical and enzymatic probing to investigate the conformation of that region, which is complementary to the 5' (+) UTR and the first 74 nucleotides of the HCV polyprotein coding sequence. By combining our experimental data with computer predictions, we have derived a secondary-structure model of this region. In our model, the last 220 nucleotides, where initiation of the (+) strand RNA synthesis presumably takes place, fold into five stable stem-loops, forming domain I. Domain I is linked to an overall less stable structure, named domain II, containing the sequences complementary to the pseudoknot of the internal ribosomal entry site in the 5' (+) UTR. Our results show that, even though the (-) strand 3'-terminal region has the antisense sequence of the 5' (+) UTR, it does not fold into its mirror image. Interestingly, comparison of the replication initiation sites on both strands reveals common structural features that may play key functions in the replication process.     

Crystal structure of Escherichia coli alkanesulfonate monooxygenase SsuD

The FMNH(2)-dependent alkanesulfonate monooxygenase SsuD catalyzes the conversion of alkanesulfonates to the corresponding aldehyde and sulfite. The enzyme allows Escherichia coli to use a wide range of alkanesulfonates as sulfur sources for growth when sulfate or cysteine are not available. The structure of SsuD was solved using the multiwavelength anomalous dispersion method from only four ordered selenium sites per asymmetric unit (one site per 20,800 Da). The final model includes 328 of 380 amino acid residues and was refined to an R-factor of 23.5% (R(free)=27.5%) at 2.3A resolution. The X-ray crystal structure of SsuD shows a homotetrameric state for the enzyme, each subunit being composed of a TIM-barrel fold enlarged by four insertion regions that contribute to intersubunit interactions. SsuD is structurally related to a bacterial luciferase and an archaeal coenzyme F(420)-dependent reductase in spite of a low level of sequence identity with these enzymes. The structural relationship is not limited to the beta-barrel region; it includes most but not all extension regions and shows distinct properties for the SsuD TIM-barrel. A likely substrate-binding site is postulated on the basis of the SsuD structure presented here, results from earlier biochemical studies, and structure relatedness to bacterial luciferase. SsuD is related to other FMNH(2)-dependent monooxygenases that show distant sequence relationship to luciferase. Thus, the structure reported here provides a model for enzymes belonging to this family and suggests that they might all fold as TIM-barrel proteins.     

A parallel stranded G‐quadruplex composed of threose nucleic acid (TNA)

G‐rich sequences can adopt four‐stranded helical structures, called G‐quadruplexes, that self‐assemble around monovalent cations like sodium (Na⁺) and potassium (K⁺). Whether similar structures can be formed from xeno‐nucleic acid (XNA) polymers with a shorter backbone repeat unit is an unanswered question with significant implications on the fold space of functional XNA polymers. Here, we examine the potential for TNA (α‐l ‐threofuranosyl nucleic acid) to adopt a four‐stranded helical structure based on a planar G‐quartet motif. Using native polyacrylamide gel electrophoresis (PAGE), circular dichroism (CD) and solution‐state nuclear magnetic resonance (NMR) spectroscopy, we show that despite a backbone repeat unit that is one atom shorter than the backbone repeat unit found in DNA and RNA, TNA can self‐assemble into stable G‐quadruplex structures that are similar in thermal stability to equivalent DNA structures. However, unlike DNA, TNA does not appear to discriminate between Na⁺ and K⁺ ions, as G‐quadruplex structures form equally well in the presence of either ion. Together, these findings demonstrate that despite a shorter backbone repeat unit, TNA is capable of self‐assembling into stable G‐quadruplex structures.     

Analysis of the sequence and structural features of the left-handed beta-helical fold

The left-handed parallel beta-helix (LbetaH) is a structurally repetitive, highly regular, and symmetrical fold formed by coiling of elongated beta-sheets into helical "rungs." This canonical fold has recently received interest as a possible solution to the fibril structure of amyloid and as a building block of self-assembled nanotubular structures. In light of this interest, we aimed to understand the structural requirements of the LbetaH fold. We first sought to determine the sequence characteristics of the repeats by analyzing known structures to identify positional preferences of specific residues types. We then used molecular dynamics simulations to demonstrate the stabilizing effect of successive rungs and the hydrophobic core of the LbetaH. We show that a two-rung structure is the minimally stable LbetaH structure. In addition, we defined the structure-based sequence preference of the LbetaH and undertook a genome-wide sequence search to determine the prevalence of this unique protein fold. This profile-based LbetaH search algorithm predicted a large fraction of LbetaH proteins from microbial origins. However, the relative number of predicted LbetaH proteins per specie was approximately equal across the genomes from prokaryotes to eukaryotes.     

Physical origins of protein superfamilies

In this work, we discovered a fundamental connection between selection for protein stability and emergence of preferred structures of proteins. Using a standard exact three-dimensional lattice model we evolve sequences starting from random ones and determine the exact native structure after each mutation. Acceptance of mutations is biased to select for stable proteins. We found that certain structures, "wonderfolds", are independently discovered numerous times as native states of stable proteins in many unrelated runs of selection. The strong dependence of lattice fold usage on the structural determinant of designability quantitatively reproduces uneven fold usage in natural proteins. Diversity of sequences that fold into wonderfold structures gives rise to superfamilies, i.e. sets of dissimilar sequences that fold into the same or very similar structures. The present work establishes a model of pre-biotic structure selection, which identifies dominant structural patterns emerging upon optimization of proteins for survival in a hot environment. Convergently discovered pre-biotic initial superfamilies with wonderfold structures could have served as a seed for subsequent biological evolution involving gene duplications and divergence.     

The exception that confirms the rule: a higher-order telomeric G-quadruplex structure more stable in sodium than in potassium

DNA and RNA guanine-quadruplexes (G4s) are stabilized by several cations, in particular by potassium and sodium ions. Generally, potassium stabilizes guanine-quartet assemblies to a larger extent than sodium; in this article we report about a higher-order G4 structure more stable in sodium than in potassium. Repeats of the DNA GGGTTA telomeric motif fold into contiguous G4 units. Using three independent approaches (thermal denaturation experiments, isothermal molecular-beacon and protein-binding assays), we show that the (GGGTTA)₇GGG sequence, folding into two contiguous G4 units, exhibits an unusual feature among G4 motifs: despite a lower thermal stability, its sodium conformation is more stable than its potassium counterpart at physiological temperature. Using differential scanning calorimetry and mutated sequences, we show that this switch in the relative stability of the sodium and potassium conformations (occurring around 45°C in 100 mM cation concentration) is the result of a more favorable enthalpy change upon folding in sodium, generated by stabilizing interactions between the two G4 units in the sodium conformation. Our work demonstrates that interactions between G4 structural domains can make a higher-order structure more stable in sodium than in potassium, even though its G4 structural domains are individually more stable in potassium than in sodium.     

Improving sequence-based fold recognition by using 3D model quality assessment

MOTIVATION: The ability of a simple method (MODCHECK) to determine the sequence-structure compatibility of a set of structural models generated by fold recognition is tested in a thorough benchmark analysis. Four Model Quality Assessment Programs (MQAPs) were tested on 188 targets from the latest LiveBench-9 automated structure evaluation experiment. We systematically test and evaluate whether the MQAP methods can successfully detect native-like models. RESULTS: We show that compared with the other three methods tested MODCHECK is the most reliable method for consistently performing the best top model selection and for ranking the models. In addition, we show that the choice of model similarity score used to assess a model's similarity to the experimental structure can influence the overall performance of these tools. Although these MQAP methods fail to improve the model selection performance for methods that already incorporate protein three dimension (3D) structural information, an improvement is observed for methods that are purely sequence-based, including the best profile-profile methods. This suggests that even the best sequence-based fold recognition methods can still be improved by taking into account the 3D structural information. CONTACT: d.jones@cs.ucl.ac.uk     

Fold independent structural comparisons of protein-ligand binding sites for exploring functional relationships

The rapid growth in protein structural data and the emergence of structural genomics projects have increased the need for automatic structure analysis and tools for function prediction. Small molecule recognition is critical to the function of many proteins; therefore, determination of ligand binding site similarity is important for understanding ligand interactions and may allow their functional classification. Here, we present a binding sites database (SitesBase) that given a known protein-ligand binding site allows rapid retrieval of other binding sites with similar structure independent of overall sequence or fold similarity. However, each match is also annotated with sequence similarity and fold information to aid interpretation of structure and functional similarity. Similarity in ligand binding sites can indicate common binding modes and recognition of similar molecules, allowing potential inference of function for an uncharacterised protein or providing additional evidence of common function where sequence or fold similarity is already known. Alternatively, the resource can provide valuable information for detailed studies of molecular recognition including structure-based ligand design and in understanding ligand cross-reactivity. Here, we show examples of atomic similarity between superfamily or more distant fold relatives as well as between seemingly unrelated proteins. Assignment of unclassified proteins to structural superfamiles is also undertaken and in most cases substantiates assignments made using sequence similarity. Correct assignment is also possible where sequence similarity fails to find significant matches, illustrating the potential use of binding site comparisons for newly determined proteins.     

All-atom ab initio folding of a diverse set of proteins

Natural proteins fold to a unique, thermodynamically dominant state. Modeling of the folding process and prediction of the native fold of proteins are two major unsolved problems in biophysics. Here, we show successful all-atom ab initio folding of a representative diverse set of proteins by using a minimalist transferable-energy model that consists of two-body atom-atom interactions, hydrogen bonding, and a local sequence-energy term that models sequence-specific chain stiffness. Starting from a random coil, the native-like structure was observed during replica exchange Monte Carlo (REMC) simulation for most proteins regardless of their structural classes; the lowest energy structure was close to native-in the range of 2-6 A root-mean-square deviation (rmsd). Our results demonstrate that the successful folding of a protein chain to its native state is governed by only a few crucial energetic terms.     

Divergent functional properties of the ribosome-associated molecular chaperone Ssb compared with other Hsp70s.

Ssbs of Saccharomyces cerevisiae are ribosome-associated molecular chaperones, which can be cross-linked to nascent polypeptide chains. Because Ssbs are members of a divergent subclass of Hsp70s found thus far only in fungi, we asked if the structural requirements for in vivo function were similar to those of "classic" Hsp70s. An intact peptide-binding domain is essential and an alteration of a conserved residue in the peptide-binding cleft (V442) affects function. However, Ssb tolerates a number of alterations in the peptide-binding cleft, revealing a high degree of flexibility in its functional requirements. Because binding of Ssb to peptide substrates in vitro was undetectable, we assessed the importance of substrate binding using the chimera BAB, in which the peptide binding domain of Ssb is exchanged for the analogous domain of the more "classical" Hsp70, Ssa. BAB, which binds peptide substrates in vitro, can substitute for Ssb in vivo. Alteration of a residue in the peptide-binding cleft of BAB creates a protein with a reduced affinity for peptide and altered ribosome binding that is unable to substitute for Ssb in vivo. These results indicate that Ssb's ability to bind unfolded polypeptides is likely critical for its function. This binding accounts, in part, for its stable interaction with translating ribosomes, even although it has a low affinity for peptides that detectably bind to other Hsp70s in vitro. These unusual properties may allow Ssb to function efficiently as a chaperone for ribosome-bound nascent chains.     

The nsp9 replicase protein of SARS-coronavirus, structure and functional insights

As part of a high-throughput structural analysis of SARS-coronavirus (SARS-CoV) proteins, we have solved the structure of the non-structural protein 9 (nsp9). This protein, encoded by ORF1a, has no designated function but is most likely involved with viral RNA synthesis. The protein comprises a single beta-barrel with a fold previously unseen in single domain proteins. The fold superficially resembles an OB-fold with a C-terminal extension and is related to both of the two subdomains of the SARS-CoV 3C-like protease (which belongs to the serine protease superfamily). nsp9 has, presumably, evolved from a protease. The crystal structure suggests that the protein is dimeric. This is confirmed by analytical ultracentrifugation and dynamic light scattering. We show that nsp9 binds RNA and interacts with nsp8, activities that may be essential for its function(s).     

Structural and Functional Characterization of the Recombinant Death Domain from Death-Associated Protein Kinase

Death-associated protein kinase (DAPk) is a calcium/calmodulin-regulated Ser/Thr-protein kinase that functions at an important point of integration for cell death signaling pathways. DAPk has a structurally unique multi-domain architecture, including a C-terminally positioned death domain (DD) that is a positive regulator of DAPk activity. In this study, recombinant DAPk-DD was observed to aggregate readily and could not be prepared in sufficient yield for structural analysis. However, DAPk-DD could be obtained as a soluble protein in the form of a translational fusion protein with the B1 domain of streptococcal protein G. In contrast to other DDs that adopt the canonical six amphipathic α-helices arranged in a compact fold, the DAPk-DD was found to possess surprisingly low regular secondary structure content and an absence of a stable globular fold, as determined by circular dichroism (CD), NMR spectroscopy and a temperature-dependent fluorescence assay. Furthermore, we measured the in vitro interaction between extracellular-regulated kinase-2 (ERK2) and various recombinant DAPk-DD constructs. Despite the low level of structural order, the recombinant DAPk-DD retained the ability to interact with ERK2 in a 1∶1 ratio with a K _d in the low micromolar range. Only the full-length DAPk-DD could bind ERK2, indicating that the apparent ‘D-motif’ located in the putative sixth helix of DAPk-DD is not sufficient for ERK2 recognition. CD analysis revealed that binding of DAPk-DD to ERK2 is not accompanied by a significant change in secondary structure. Taken together our data argue that the DAPk-DD, when expressed in isolation, does not adopt a classical DD fold, yet in this state retains the capacity to interact with at least one of its binding partners. The lack of a stable globular structure for the DAPk-DD may reflect either that its folding would be supported by interactions absent in our experimental set-up, or a limitation in the structural bioinformatics assignment of the three-dimensional structure.     

The emergence of catalytic and structural diversity within the beta-clip fold

The beta-clip fold includes a diverse group of protein domains that are unified by the presence of two characteristic waist-like constrictions, which bound a central extended region. Members of this fold include enzymes like deoxyuridine triphosphatase and the SET methylase, carbohydrate-binding domains like the fish antifreeze proteins/Sialate synthase C-terminal domains, and functionally enigmatic accessory subunits of urease and molybdopterin biosynthesis protein MoeA. In this study, we reconstruct the evolutionary history of this fold using sensitive sequence and structure comparisons methods. Using sequence profile searches, we identified novel versions of the beta-clip fold in the bacterial flagellar chaperone FlgA and the related pilus protein CpaB, the StrU-like dehydrogenases, and the UxaA/GarD-like hexuronate dehydratases (SAF superfamily). We present evidence that these versions of the beta-clip domain, like the related type III anti-freeze proteins and C-terminal domains of sialic acid synthases, are involved in interactions with carbohydrates. We propose that the FlgA and CpaB-like proteins mediate the assembly of bacterial flagella and Flp pili by means of their interactions with the carbohydrate moieties of peptidoglycan. The N-terminal beta-clip domain of the hexuronate dehydratases appears to have evolved a novel metal-binding site, while their C-terminal domain is likely to adopt a metal-binding TIM barrel-like fold. Using structural comparisons, we show that the beta-clip fold can be further classified into two major groups, one that includes the SAF, SET, dUTPase superfamilies, and the other that includes the phage lambda head decoration protein, the beta subunit of urease and the C-terminal domain of the molybdenum cofactor biosynthesis protein MoeA. Structural comparisons also suggest the beta-clip fold was assembled through the duplication of a three-stranded unit. Though the three-stranded units are likely to have had a common origin, we present evidence that complete beta-clip domains were assembled through such duplications, independently on multiple occasions. There is also evidence for circular permutation of the basic three-stranded unit on different occasions in the evolution of the beta-clip unit. We also describe how assembly of this fold from a basic three-stranded unit has been utilized to accommodate a variety of activities in its different versions.     

Structure-approximating inverse protein folding problem in the 2D HP model.

The inverse protein folding problem is that of designing an amino acid sequence which has a particular native protein fold. This problem arises in drug design where a particular structure is necessary to ensure proper protein-protein interactions. In this paper, we show that in the 2D HP model of Dill it is possible to solve this problem for a broad class of structures. These structures can be used to closely approximate any given structure. One of the most important properties of a good protein (in drug design) is its stability--the aptitude not to fold simultaneously into other structures. We show that for a number of basic structures, our sequences have a unique fold.