Globular Fold Research Articles

Mutational Analysis of the Binding-Induced Folding Reaction of the Mixed-Lineage Leukemia Protein to the KIX Domain.

Intrinsically disordered proteins represent a large class of proteins that lack a well-defined three-dimensional structure in isolation but can undergo a disorder to order transition upon binding to their physiological ligands. Understanding the mechanism by which these proteins fold upon binding represents a challenge. Here we present a detailed mutational study of the kinetics of the binding reaction between the transactivation domain of the mixed-lineage leukemia protein, an intrinsically disordered protein, and the KIX domain, performed under different experimental conditions. The experimental data allow us to infer the mechanism of folding upon binding and to pinpoint the key interactions present in the transition state. Furthermore, we identify a peculiar malleability of the observed mechanism upon changes in reaction conditions. This finding, which is in opposition to the robustness typically observed in the folding of globular proteins, is discussed in the context of previous work on intrinsically disordered proteins.

Globular Protein Folding In Vitro and In Vivo.

In vitro, computational, and theoretical studies of protein folding have converged to paint a rich and complex energy landscape. This landscape is sensitively modulated by environmental conditions and subject to evolutionary pressure on protein function. Of these environments, none is more complex than the cell itself, where proteins function in the cytosol, in membranes, and in different compartments. A wide variety of kinetic and thermodynamics experiments, ranging from single-molecule studies to jump kinetics and from nuclear magnetic resonance to imaging on the microscope, have elucidated how protein energy landscapes facilitate folding and how they are subject to evolutionary constraints and environmental perturbation. Here we review some recent developments in the field and refer the reader to some original work and additional reviews that cover this broad topic in protein science.

The first crystal structure of human RNase 6 reveals a novel substrate-binding and cleavage site arrangement.

Human RNase 6 is a cationic secreted protein that belongs to the RNase A superfamily. Its expression is induced in neutrophils and monocytes upon bacterial infection, suggesting a role in host defence. We present here the crystal structure of RNase 6 obtained at 1.72 Å (1 Å=0.1nm) resolution, which is the first report for the protein 3D structure and thereby setting the basis for functional studies. The structure shows an overall kidney-shaped globular fold shared with the other known family members. Three sulfate anions bound to RNase 6 were found, interacting with residues at the main active site (His(15), His(122) and Gln(14)) and cationic surface-exposed residues (His(36), His(39), Arg(66) and His(67)). Kinetic characterization, together with prediction of protein-nucleotide complexes by molecular dynamics, was applied to analyse the RNase 6 substrate nitrogenous base and phosphate selectivity. Our results reveal that, although RNase 6 is a moderate catalyst in comparison with the pancreatic RNase type, its structure includes lineage-specific features that facilitate its activity towards polymeric nucleotide substrates. In particular, enzyme interactions at the substrate 5' end can provide an endonuclease-typecleavage pattern. Interestingly, the RNase 6 crystal structure revealed a novel secondary active site conformed by the His(36)-His(39) dyad that facilitates the polynucleotide substrate catalysis.

Fast Closuer of Long Loops at the Initiation of a Protein Folding Pathway

What are the initial steps in the folding of globular proteins? In order to test the hypothesis that long loop closure transitions are part of the initiation of the folding pathways of globular proteins, our groups are studying the initiation of the folding transition of a model protein by time resolved excitation energy transfer (trFRET) detected fast kinetics experiments. Site specific double labeling is used to study the timing of conformational transitions of individual loop forming chain segments at the µs time regime. Previously, it was shown that at least three long loops in the Escherichia coli adenylate kinase (AK) molecule close within the first 5 ms of the 6 seconds folding kinetics of AK. We applied a continuous flow based double kinetics experiment. These measurements enabled us to obtain µs series of transient time dependent distributions of distances between the ends of the labeled loops. Analysis of the trFRET experiments shows that the N terminal loop (loop I) is closed within less than 60 µs after the initiation of refolding. Loop II is also mostly closed within that time step but shows an additional small reduction of the mean end to end distance in a second phase at a rate of 0.005µs−1. This second phase can either reflect tightening of a loosely closed loop in the ensemble of conformers or may reflect two subpopulations in the ensemble, which differ in the rate of closure of loop II, but not in the rate of closure of loop I. This study shows the very fast closure of long loops in the otherwise disordered backbone and fine details of the very early hidden pre-transition state stepsthat are essential for the fast and efficient folding of the protein molecule.

Nuclear Magnetic Resonance Solution Structures of Lacticin Q and Aureocin A53 Reveal a Structural Motif Conserved among Leaderless Bacteriocins with Broad-Spectrum Activity.

Lacticin Q (LnqQ) and aureocin A53 (AucA) are leaderless bacteriocins from Lactococcus lactis QU5 and Staphylococcus aureus A53, respectively. These bacteriocins are characterized by the absence of an N-terminal leader sequence and are active against a broad range of Gram-positive bacteria. LnqQ and AucA consist of 53 and 51 amino acids, respectively, and have 47% identical sequences. In this study, their three-dimensional structures were elucidated using solution nuclear magnetic resonance and were shown to consist of four α-helices that assume a very similar compact, globular overall fold (root-mean-square deviation of 1.7 Å) with a highly cationic surface and a hydrophobic core. The structures of LnqQ and AucA resemble the shorter two-component leaderless bacteriocins, enterocins 7A and 7B, despite having low levels of sequence identity. Homology modeling revealed that the observed structural motif may be shared among leaderless bacteriocins with broad-spectrum activity against Gram-positive organisms. The elucidated structures of LnqQ and AucA also exhibit some resemblance to circular bacteriocins. Despite their similar overall fold, inhibition studies showed that LnqQ and AucA have different antimicrobial potency against the Gram-positive strains tested, suggesting that sequence disparities play a crucial role in their mechanisms of action.

A composite approach towards a complete model of the myosin rod.

Sarcomeric myosins have the remarkable ability to form regular bipolar thick filaments that, together with actin thin filaments, constitute the fundamental contractile unit of skeletal and cardiac muscle. This has been established for over 50 years and yet a molecular model for the thick filament has not been attained. In part this is due to the lack of a detailed molecular model for the coiled-coil that constitutes the myosin rod. The ability to self-assemble resides in the C-terminal section of myosin known as light meromyosin (LMM) which exhibits strong salt-dependent aggregation that has inhibited structural studies. Here we evaluate the feasibility of generating a complete model for the myosin rod by combining overlapping structures of five sections of coiled-coil covering 164 amino acid residues which constitute 20% of LMM. Each section contains ∼ 7-9 heptads of myosin. The problem of aggregation was overcome by incorporating the globular folding domains, Gp7 and Xrcc4 which enhance crystallization. The effect of these domains on the stability and conformation of the myosin rod was examined through biophysical studies and overlapping structures. In addition, a computational approach was developed to combine the sections into a contiguous model. The structures were aligned, trimmed to form a contiguous model, and simulated for >700 ns to remove the discontinuities and achieve an equilibrated conformation that represents the native state. This experimental and computational strategy lays the foundation for building a model for the entire myosin rod.

Crystal Structure of the Human Cytomegalovirus pUL50-pUL53 Core Nuclear Egress Complex Provides Insight into a Unique Assembly Scaffold for Virus-Host Protein Interactions

Nuclear replication of cytomegalovirus relies on elaborate mechanisms of nucleocytoplasmic egress of viral particles. Thus, the role of two essential and conserved viral nuclear egress proteins, pUL50 and pUL53, is pivotal. pUL50 and pUL53 heterodimerize and form a core nuclear egress complex (NEC), which is anchored to the inner nuclear membrane and provides a scaffold for the assembly of a multimeric viral-cellular NEC. Here, we report the crystal structure of the pUL50-pUL53 heterodimer (amino acids 1-175 and 50-292, respectively) at 2.44 Å resolution. Both proteins adopt a globular fold with mixed α and β secondary structure elements. pUL53-specific features include a zinc-binding site and a hook-like N-terminal extension, the latter representing a hallmark element of the pUL50-pUL53 interaction. The hook-like extension (amino acids 59-87) embraces pUL50 and contributes 1510 Å(2) to the total interface area (1880 Å(2)). The pUL50 structure overall resembles the recently published NMR structure of the murine cytomegalovirus homolog pM50 but reveals a considerable repositioning of the very C-terminal α-helix of pUL50 upon pUL53 binding. pUL53 shows structural resemblance with the GHKL domain of bacterial sensory histidine kinases. A close examination of the crystal structure indicates partial assembly of pUL50-pUL53 heterodimers to hexameric ring-like structures possibly providing additional scaffolding opportunities for NEC. In combination, the structural information on pUL50-pUL53 considerably improves our understanding of the mechanism of HCMV nuclear egress. It may also accelerate the validation of the NEC as a unique target for developing a novel type of antiviral drug and improved options of broad-spectrum antiherpesviral therapy.

Dissecting and reprogramming the folding and assembly of tandem-repeat proteins.

Studying protein folding and protein design in globular proteins presents significant challenges because of the two related features, topological complexity and co-operativity. In contrast, tandem-repeat proteins have regular and modular structures composed of linearly arrayed motifs. This means that the biophysics of even giant repeat proteins is highly amenable to dissection and to rational design. Here we discuss what has been learnt about the folding mechanisms of tandem-repeat proteins. The defining features that have emerged are: (i) accessibility of multiple distinct routes between denatured and native states, both at equilibrium and under kinetic conditions; (ii) different routes are favoured for folding compared with unfolding; (iii) unfolding energy barriers are broad, reflecting stepwise unravelling of an array repeat by repeat; (iv) highly co-operative unfolding at equilibrium and the potential for exceptionally high thermodynamic stabilities by introducing consensus residues; (v) under force, helical-repeat structures are very weak with non-co-operative unfolding leading to elasticity and buffering effects. This level of understanding should enable us to create repeat proteins with made-to-measure folding mechanisms, in which one can dial into the sequence the order of repeat folding, number of pathways taken, step size (co-operativity) and fine-structure of the kinetic energy barriers.

Binding interaction of a novel fluorophore with serum albumins: steady state fluorescence perturbation and molecular modeling analysis.

Fluorescence emission and anisotropy are widely used to measure the binding parameters and kinetic behavior of reactions that cause a change in the rotational time of a fluorescent molecule. We report here fluorescence emission and anisotropy behavior of a newly synthesized novel naphthalene base fluorophore (methyl 3-[(6-{[2-(tert-butoxy)-2-oxoethyl] (4-methoxyphenyl)amino}naphthalen-2-yl)formamido]propanoate) in several solution conditions including its binding to human and bovine serum albumin proteins both in their native and denatured states. The fluorescence yield of the compound substantially increased inside hydrophobic protein surface and ~30 nm decrease in Stokes’ shift, compared to aqueous solution, was observed. Shift in fluorescence excitation peak position from the absorption peak of the molecule was ~8 nm in protein solution. This indicated possible alteration of excited state geometry of the compound by the globular fold of albumins. In addition, we measured the steady state fluorescence anisotropy of the molecule to evaluate several thermodynamic parameters and the results suggested the binding was energetically favorable. The measured ΔG° was ~−30 kJ mol−1 and the derived dissociation constant was ~10−6 M. The molecular docking analysis further highlighted the nonspecific association of the compound with the proteins and hydrophobic forces may have a significant role in the binding processes. Under the denatured condition of the protein, the compound lost its binding efficacy and reduction in fluorescence intensity was observed. Thus, the molecule appears as a new fluorescence probe to report the nature of its binding site in terms of increased fluorescence quantum yield and decreased Stokes’ shift. It can also report the changes in the binding site due to global change in protein structure such as unfolding/misfolding often linked to several human disorder. Further it could be useful to detect and study the drug binding site of specific protein of interest.Electronic supplementary materialThe online version of this article (doi:10.1186/s40064-015-1333-8) contains supplementary material, which is available to authorized users.

TRNA acceptor stem and anticodon bases form independent codes related to protein folding

Aminoacyl-tRNA synthetases recognize tRNA anticodon and 3' acceptor stem bases. Synthetase Urzymes acylate cognate tRNAs even without anticodon-binding domains, in keeping with the possibility that acceptor stem recognition preceded anticodon recognition. Representing tRNA identity elements with two bits per base, we show that the anticodon encodes the hydrophobicity of each amino acid side-chain as represented by its water-to-cyclohexane distribution coefficient, and this relationship holds true over the entire temperature range of liquid water. The acceptor stem codes preferentially for the surface area or size of each side-chain, as represented by its vapor-to-cyclohexane distribution coefficient. These orthogonal experimental properties are both necessary to account satisfactorily for the exposed surface area of amino acids in folded proteins. Moreover, the acceptor stem codes correctly for β-branched and carboxylic acid side-chains, whereas the anticodon codes for a wider range of such properties, but not for size or β-branching. These and other results suggest that genetic coding of 3D protein structures evolved in distinct stages, based initially on the size of the amino acid and later on its compatibility with globular folding in water.

Universal distribution of mutational effects on protein stability, uncoupling of protein robustness from sequence evolution and distinct evolutionary modes of prokaryotic and eukaryotic proteins

Robustness to destabilizing effects of mutations is thought of as a key factor of protein evolution. The connections between two measures of robustness, the relative core size and the computationally estimated effect of mutations on protein stability (ΔΔG), protein abundance and the selection pressure on protein-coding genes (dN/dS) were analyzed for the organisms with a large number of available protein structures including four eukaryotes, two bacteria and one archaeon. The distribution of the effects of mutations in the core on protein stability is universal and indistinguishable in eukaryotes and bacteria, centered at slightly destabilizing amino acid replacements, and with a heavy tail of more strongly destabilizing replacements. The distribution of mutational effects in the hyperthermophilic archaeon Thermococcus gammatolerans is significantly shifted toward strongly destabilizing replacements which is indicative of stronger constraints that are imposed on proteins in hyperthermophiles. The median effect of mutations is strongly, positively correlated with the relative core size, in evidence of the congruence between the two measures of protein robustness. However, both measures show only limited correlations to the expression level and selection pressure on protein-coding genes. Thus, the degree of robustness reflected in the universal distribution of mutational effects appears to be a fundamental, ancient feature of globular protein folds whereas the observed variations are largely neutral and uncoupled from short term protein evolution. A weak anticorrelation between protein core size and selection pressure is observed only for surface residues in prokaryotes but a stronger anticorrelation is observed for all residues in eukaryotic proteins. This substantial difference between proteins of prokaryotes and eukaryotes is likely to stem from the demonstrable higher compactness of prokaryotic proteins.

Structure and Assembly of Functional Amyloids

The ability of proteins to self‐assemble into beta‐sheet‐rich aggregates called amyloid fibrils is considered to be universal, although certain polypeptide sequences have a particularly high propensity to adopt these conformations. In many cases the formation of amyloid fibrils is deleterious and associated with the progression of disease, but there are also examples of proteins for which the cross‐beta structure represents the functional conformation. Ure2 is the protein determinant of the yeast prion [URE3]. Ure2 consists of an N‐terminal prion‐inducing domain that is disordered in the native state, whereas the C‐terminal functional domain has a globular fold with structural similarity to glutathione transferase enzymes. The C‐terminal domain shows enzymatic activity in both the soluble and fibrillar forms of Ure2. We have used a variety of biophysical approaches to investigate the structure of Ure2 fibrils and their mechanism of assembly. We have also created chimeric constructs where the prion domain is genetically fused to other enzymes of different sizes and architectures. These chimeric polypeptide constructs spontaneously self‐assemble into nanofibrils with fused active enzyme subunits displayed on the amyloid fibril surface. We can measure steady‐state kinetic parameters for the appended enzymes in situ within fibrils, and compare these for the identical protein constructs in solution. We are currently investigating the use of microfluidic techniques to form enzymatically‐active microgel particles from the chimeric self‐assembling protein nanofibrils. These approaches allow amyloid fibrils to be used as a scaffold for enzyme immobilization, facilitating recycling and reuse of the appended enzymes in functional assays, both in batch processes and for microflow chemistry.

An Efficiently-folding Purely-symmetric de novo Designed Protein

Many globular protein folds exhibit some form of rotational symmetry – the most common example being the TIM barrel (having 8-fold rotational symmetry of a repeating beta-strand/turn/alpha- helix/turn motif). This type of common symmetry in proteins was apparent from the earliest days of X-ray structure studies, leading to the hypothesis that gene duplication and fusion was the underlying evolutionary mechanism. However, while such symmetry at the 3°- structure level is apparent, any 1°-structure symmetry when comparing repeating structural motifs is often largely absent. Such sequence analyses, along with theoretical considerations of folding frustration for exact repeating motifs, as well as natively-unstructured properties for peptides having reduced sequence complexity (as occurs with exact repeating motifs), resulted in a paradigm that proteins designed with exact symmetry are unlikely to be efficiently folding. Given that exact symmetry can substantially reduce the combinatorial explosion problem inherent to protein design, turning such problems from impossible to computationally tractable, the demonstration of efficient folding for purely symmetric de novo designed proteins took on considerable interest.

Structural and Dynamic Analysis on Disordered H4 Histone Tail by Modified AWSEM-MD

DNA compaction in eukaryotic cells is mediated by positively charged octamers comprised of histone proteins. The latter consists of well folded core segments, that come together to form a central cylinder, and flexible tails, protruding out from cylinder's rim. Despite being disordered, histone tails play an important role in bridging interactions between neighboring nucleosomes, regulating folding structure and dynamics of chromatin fibers. Histone tails, in turn, are mainly regulated via post-translational modifications, such as methylation and acetylation at various positions. Because of their flexibility and disordered nature, it has been difficult to investigate histone tails both computationally and experimentally. In particular, it is desirable to develop coarse-grained, yet accurate models of histone tails, such that subsequent nucleosomal and polynucleosmal simulations could be carried out within feasible times. To achieve this goal, we added new interactions to the associative memory, water mediated, structure and energy model (AWSEM-MD), which is typically used for folding of globular proteins or binding studies. We found that modified AWSEM-MD reproduces well the complex conformational ensemble of the H4 histone tail, obtained from atomistic simulations with explicit solvent. In particular, the cumulative and site-specific effects of various acetylation combinations are consistent with the all-atom results. Our proposed extension of AWSEM-MD may allow simulating intrinsically disordered proteins with high accuracy and computational efficiency.

Fast Closure of Long Loops at the Initiation of Folding of Globular Proteins Studied by FRET Based Methods

The rate of protein folding is determined by the rate of passage through the transition state, however major structural transition precede the TSE formation. We hypothesize that few non-local interactions are effective in the early phases of the folding transition prior to the cooperative transition. These interactions loosely stabilize few closed loops which form the folding noncontiguous nucleus, reduce the chain entropy and determine the course of the folding pathway (the “loop hypothesis”). We study the order of formation of secondary structure elements and loop closure during the early phases of the folding of E. coli adenylate kinase (AK) by combination of rapid mixing methods and time resolved FRET spectroscopy. We find that at the initiation of folding of the AK molecule the two N terminal closed loop structures in the CORE domain reach native end to end distance within the first 300 microseconds of the transition. Three representative CORE domain β-strands have non-native end to end distance during the first 15 ms and undergo slow change (3 sec) to native distance. Not all helical segments in the CORE domain fold at the slow rate as the beta strands. We conclude that non local interactions are essential factor at the early phases of the folding transition and that the folding of sub-domain elements is context dependent and should be studied in the whole molecule, in situ.

Some insights into the transition state ensemble of the folding of globular proteins

Some insights into the transition state ensemble of the folding of globular proteins

Predictive energy landscapes for folding α-helical transmembrane proteins.

We explore the hypothesis that the folding landscapes of membrane proteins are funneled once the proteins' topology within the membrane is established. We extend a protein folding model, the associative memory, water-mediated, structure, and energy model (AWSEM) by adding an implicit membrane potential and reoptimizing the force field to account for the differing nature of the interactions that stabilize proteins within lipid membranes, yielding a model that we call AWSEM-membrane. Once the protein topology is set in the membrane, hydrophobic attractions play a lesser role in finding the native structure, whereas polar-polar attractions are more important than for globular proteins. We examine both the quality of predictions made with AWSEM-membrane when accurate knowledge of the topology and secondary structure is available and the quality of predictions made without such knowledge, instead using bioinformatically inferred topology and secondary structure based on sequence alone. When no major errors are made by the bioinformatic methods used to assign the topology of the transmembrane helices, these two types of structure predictions yield roughly equivalent quality structures. Although the predictive energy landscape is transferable and not structure based, within the correct topological sector we find the landscape is indeed very funneled: Thermodynamic landscape analysis indicates that both the total potential energy and the contact energy decrease as native contacts are formed. Nevertheless the near symmetry of different helical packings with respect to native contact formation can result in multiple packings with nearly equal thermodynamic occupancy, especially at temperatures just below collapse.

Long-range effects of tag sequence on marginally stabilized structure in HIV-1 p24 capsid protein monitored using NMR

N-terminal domain of HIV-1 p24 capsid protein is a globular fold composed of seven helices and two β-strands with a flexible structure including the α4–5 loop and both N- and C-terminal ends. However, the protein shows a high tendency (48%) for an intrinsically disordered structure based on the PONDR VL-XT prediction from the primary sequence. To assess the possibility of marginally stabilized structure under physiological conditions, the N-terminal domain of p24 was destabilized by the addition of an artificial flexible tag to either N- or C-terminal ends, and it was analyzed using T1, T2, hetero-nuclear NOE, and amide-proton exchange experiments. When the C-terminal tag (12 residues) was attached, the regions of the α3–4 loop and helix 6 as well as the α4–5 loop attained the flexible structures. Furthermore, in the protein containing the N-terminal tag (27 residues), helix 4 in addition to the above-mentioned area including α3–4 and α4–5 loops as well as helix 6 exhibited highly disordered structures. Thus, the long-range effects of the existence of tag sequence was observed in the stepwise manner of the appearance of disordered structures (step 1: α4–5 loop, step 2: α3–4 loop and helix 6, and step 3: helix 4). Furthermore, the disordered regions in tagged proteins were consistent with the PONDR VL-XT disordered prediction. The dynamic structure located in the middle part (α3–4 loop to helix 6) of the protein shown in this study may be related to the assembly of the viral particle.

Quantitative thermodynamic model for globular protein folding

We present a statistical mechanics formalism for theoretical description of the process of protein folding ↔ unfolding transition in water environment. The formalism is based on the construction of the partition function of a protein obeying two-stage-like folding kinetics. Using the statistical mechanics model of solvation of hydrophobic hydrocarbons we obtain the partition function of infinitely diluted solution of proteins in water environment. The calculated dependencies of the protein heat capacities upon temperature are compared with the corresponding results of experimental measurements for staphylococcal nuclease and metmyoglobin.

Mechanism of Folding and Binding of an Intrinsically Disordered Protein As Revealed by ab Initio Simulations.

A complex of the phosphorylated kinase-inducible domain (pKID) with its interacting domain (KIX) is a model system for studies of mechanisms by which intrinsically unfolded proteins perform their functions. These mechanisms are not fully understood. Using an efficient coarse-grained model, ab initio simulations were performed of the coupled folding and binding of the pKID to the KIX. The simulations start from an unbound, randomly positioned and disordered pKID structure. During the simulations the pKID chain and its position remain completely unrestricted, while the KIX backbone is limited to near-native fluctuations. Ab initio simulations of such large-scale conformational transitions, unaffected by any knowledge about the bound pKID structure, remain inaccessible to classical simulations. Our simulations recover an ensemble of transient encounter complexes in good agreement with experimental results. We find that a key folding and binding step is linked to the formation of weak native interactions between a preformed nativelike fragment of a pKID helix and KIX surface. Once that nucleus forms, the pKID chain may condense from a largely disordered encounter ensemble to a natively bound and ordered conformation. The observed mechanism is reminiscent of a nucleation-condensation model, a common scenario for folding of globular proteins.