Plasticity within the obligatory folding nucleus of an immunoglobulin-like domain

Pathogenic variants in the fibronectin type III domain of leptin receptor: Molecular dynamics simulation and structural analysis

Several case reports have identified leptin receptor (LEPR) variants associated with severe obesity in humans. However, the structure of LEPR has only been partially understood until recently, and few studies have investigated the detrimental effects of these variants on the protein's three-dimensional structure. Notably, fibronectin type III (FnIII) domains play a crucial role in signal transduction. In this study, we examined the impact of 10 variants within the FnIII domains on LEPR structure using molecular dynamics (MD) simulations and structural analysis. Our 300 ns MD simulations revealed that the C604S variant, which disrupts a key disulfide bond, significantly increased the overall root-mean-square deviation (RMSD) of the FnIII-2 and FnIII-3 domains, indicating destabilization of the interdomain rigidity required for proper signaling. Variants such as P639L, N718S, and W646C also induced abnormal bending and rotational misalignment between the FnIII domains, contributing to interdomain destabilization. Structural analysis identified folding nuclei and demonstrated that L662S, W664R, H684P, and S723F destabilize the internal domain. Variants affecting interdomain resulted in lower-than-expected damage prediction scores by bioinformatics tools. This study is expected to contribute to the elucidation of the disease-causing mechanisms of missense variants in the leptin receptor.

From covalent transition states in chemistry to noncovalent in biology: from β- to Φ-value analysis of protein folding

Solving the mechanism of a chemical reaction requires determining the structures of all the ground states on the pathway and the elusive transition states linking them. 2024 is the centenary of Brønsted's landmark paper that introduced the β-value and structure-activity studies as the only experimental means to infer the structures of transition states. It involves making systematic small changes in the covalent structure of the reactants and analysing changes in activation and equilibrium-free energies. Protein engineering was introduced for an analogous procedure, Φ-value analysis, to analyse the noncovalent interactions in proteins central to biological chemistry. The methodology was developed first by analysing noncovalent interactions in transition states in enzyme catalysis. The mature procedure was then applied to study transition states in the pathway of protein folding - 'part (b) of the protein folding problem'. This review describes the development of Φ-value analysis of transition states and compares and contrasts the interpretation of β- and Φ-values and their limitations. Φ-analysis afforded the first description of transition states in protein folding at the level of individual residues. It revealed the nucleation-condensation folding mechanism of protein domains with the transition state as an expanded, distorted native structure, containing little fully formed secondary structure but many weak tertiary interactions. A spectrum of transition states with various degrees of structural polarisation was then uncovered that spanned from nucleation-condensation to the framework mechanism of fully formed secondary structure. Φ-analysis revealed how movement of the expanded transition state on an energy landscape accommodates the transition from framework to nucleation-condensation mechanisms with a malleability of structure as a unifying feature of folding mechanisms. Such movement follows the rubric of analysis of classical covalent chemical mechanisms that began with Brønsted. Φ-values are used to benchmark computer simulation, and Φ and simulation combine to describe folding pathways at atomic resolution.

Protein folding problem: enigma, paradox, solution

The ability of protein chains to spontaneously form their three-dimensional structures is a long-standing mystery in molecular biology. The most conceptual aspect of this mystery is how the protein chain can find its native, "working" spatial structure (which, for not too big protein chains, corresponds to the global free energy minimum) in a biologically reasonable time, without exhaustive enumeration of all possible conformations, which would take billions of years. This is the so-called "Levinthal's paradox." In this review, we discuss the key ideas and discoveries leading to the current understanding of protein folding kinetics, including folding landscapes and funnels, free energy barriers at the folding/unfolding pathways, and the solution of Levinthal's paradox. A special role here is played by the "all-or-none" phase transition occurring at protein folding and unfolding and by the point of thermodynamic (and kinetic) equilibrium between the "native" and the "unfolded" phases of the protein chain (where the theory obtains the simplest form). The modern theory provides an understanding of key features of protein folding and, in good agreement with experiments, it (i) outlines the chain length-dependent range of protein folding times, (ii) predicts the observed maximal size of "foldable" proteins and domains. Besides, it predicts the maximal size of proteins and domains that fold under solely thermodynamic (rather than kinetic) control. Complementarily, a theoretical analysis of the number of possible protein folding patterns, performed at the level of formation and assembly of secondary structures, correctly outlines the upper limit of protein folding times.

Genotype-Structure-Phenotype Correlations of Disease-Associated IGF1R Variants and Similarities to Those of INSR Variants

We previously reported genotype-phenotype correlations in 12 missense variants causing severe insulin resistance, located in the second and third fibronectin type III (FnIII) domains of the insulin receptor (INSR), containing the α-β cleavage and part of insulin-binding sites. This study aimed to identify genotype-phenotype correlations in FnIII domain variants of IGF1R, a structurally related homolog of INSR, which may be associated with growth retardation, using the recently reported crystal structures of IGF1R. A structural bioinformatics analysis of five previously reported disease-associated heterozygous missense variants and a likely benign variant in the FnIII domains of IGF1R predicted that the disease-associated variants would severely impair the hydrophobic core formation and stability of the FnIII domains or affect the α-β cleavage site, while the likely benign variant would not affect the folding of the domains. A functional analysis of these variants in CHO cells showed impaired receptor processing and autophosphorylation in cells expressing the disease-associated variants but not in those expressing the wild-type form or the likely benign variant. These results demonstrated genotype-phenotype correlations in the FnIII domain variants of IGF1R, which are presumably consistent with those of INSR and would help in the early diagnosis of patients with disease-associated IGF1R variants.

Nanobody stability engineering by employing the ΔTm shift; a comparison with apparent rate constants of heat-induced aggregation

The antigen-binding domains of camelid heavy-chain antibodies, also called nanobodies, gained strong attention because of their unique functional and biophysical properties. They gave rise to an entire spectrum of applications in biotechnology, research and medicine. Despite several reports about reversibly refolding nanobodies, protein aggregation plays a major role in nanobody thermoresistance, asking for strategies to engineer their refolding behavior. Here, we use measurements of nanobody aggregation kinetics to validate structural features in the nanobody fold that are suppressing heat-induced nanobody aggregation. Furthermore, the kinetic measurements yielded a detailed insight into the concept of the ΔTm shift, a metric for protein aggregation propensities obtained from differential scanning fluorimetry measurements. By relating the equilibrium measurements of the ΔTm shift to the kinetic measurements of heat-induced nanobody aggregation, a distinct relationship could be identified that allows a prediction of nanobody aggregation rates from a simple equilibrium measurement of ΔTm.

Investigating the Effect of Chain Connectivity on the Folding of a Beta-Sheet Protein On and Off the Ribosome

Determining the relationship between protein folding pathways on and off the ribosome remains an important area of investigation in biology. Studies on isolated domains have shown that alteration of the separation of residues in a polypeptide chain, while maintaining their spatial contacts, may affect protein stability and folding pathway. Due to the vectorial emergence of the polypeptide chain from the ribosome, chain connectivity may have an important influence upon cotranslational folding. Using MATH, an all β-sandwich domain, we investigate whether the connectivity of residues and secondary structure elements is a key determinant of when cotranslational folding can occur on the ribosome. From Φ-value analysis, we show that the most structured region of the transition state for folding in MATH includes the N and C terminal strands, which are located adjacent to each other in the structure. However, arrest peptide force-profile assays show that wild-type MATH is able to fold cotranslationally, while some C-terminal residues remain sequestered in the ribosome, even when destabilized by 2–3 kcal mol⁻¹. We show that, while this pattern of Φ-values is retained in two circular permutants in our studies of the isolated domains, one of these permutants can fold only when fully emerged from the ribosome. We propose that in the case of MATH, onset of cotranslational folding is determined by the ability to form a sufficiently stable folding nucleus involving both β-sheets, rather than by the location of the terminal strands in the ribosome tunnel.

•

Adjacent N and C terminal strands are most structured region in the transition state.

•

Two circular permutants retain the same folding pathway as wild-type MATH.

•

On the ribosome, early emergence of terminal strands does not promote earlier folding.

•

Formation of both β-sheets is energetically critical for folding on the ribosome.

•

Folding pathway minimizes formation of partly structured states prone to mis-folding.

Nonnative Energetic Frustrations in Protein Folding at Residual Level: A Simulation Study of Homologous Immunoglobulin-like β-Sandwich Proteins

Nonnative interactions cause energetic frustrations in protein folding and were found to dominate key events in folding intermediates. However, systematically characterizing energetic frustrations that are caused by nonnative intra-residue interactions at residual resolution is still lacking. Recently, we studied the folding of a set of homologous all-α proteins and found that nonnative-contact-based energetic frustrations are highly correlated to topology of the protein native-contact network. Here, we studied the folding of nine homologous immunoglobulin-like (Ig-like) β-sandwich proteins, and examined nonnative-contact-based energetic frustrations Gō-like model. Our calculations showed that nonnative-interaction-based energetic frustrations in β-sandwich proteins are much more complicated than those in all-α proteins, and they exhibit highly heterogeneous effects on the folding of secondary structures. Further, the nonnative interactions introduced distinct correlations in the folding of different folding-patches of β-sandwich proteins. Taken together, a strong interplay might exist between nonnative-interaction energetic frustrations and the protein native-contact networks, which ensures that β-sandwich domains adopt a common folding mechanism.

Structural Basis and Genotype-Phenotype Correlations of INSR Mutations Causing Severe Insulin Resistance

The insulin receptor (INSR) gene was analyzed in four patients with severe insulin resistance, revealing five novel mutations and a deletion that removed exon 2. A patient with Donohue syndrome (DS) had a novel p.V657F mutation in the second fibronectin type III domain (FnIII-2), which contains the α-β cleavage site and part of the insulin-binding site. The mutant INSR was expressed in Chinese hamster ovary cells, revealing that it reduced insulin proreceptor processing and impaired activation of downstream signaling cascades. Using online databases, we analyzed 82 INSR missense mutations and demonstrated that mutations causing DS were more frequently located in the FnIII domains than those causing the milder type A insulin resistance (P = 0.016). In silico structural analysis revealed that missense mutations predicted to severely impair hydrophobic core formation and stability of the FnIII domains all caused DS, whereas those predicted to produce localized destabilization and to not affect folding of the FnIII domains all caused the less severe Rabson-Mendenhall syndrome. These results suggest the importance of the FnIII domains, provide insight into the molecular mechanism of severe insulin resistance, will aid early diagnosis, and will provide potential novel targets for treating extreme insulin resistance.

There and back again: Two views on the protein folding puzzle

The ability of protein chains to spontaneously form their spatial structures is a long-standing puzzle in molecular biology. Experimentally measured folding times of single-domain globular proteins range from microseconds to hours: the difference (10-11 orders of magnitude) is the same as that between the life span of a mosquito and the age of the universe. This review describes physical theories of rates of overcoming the free-energy barrier separating the natively folded (N) and unfolded (U) states of protein chains in both directions: "U-to-N" and "N-to-U". In the theory of protein folding rates a special role is played by the point of thermodynamic (and kinetic) equilibrium between the native and unfolded state of the chain; here, the theory obtains the simplest form. Paradoxically, a theoretical estimate of the folding time is easier to get from consideration of protein unfolding (the "N-to-U" transition) rather than folding, because it is easier to outline a good unfolding pathway of any structure than a good folding pathway that leads to the stable fold, which is yet unknown to the folding protein chain. And since the rates of direct and reverse reactions are equal at the equilibrium point (as follows from the physical "detailed balance" principle), the estimated folding time can be derived from the estimated unfolding time. Theoretical analysis of the "N-to-U" transition outlines the range of protein folding rates in a good agreement with experiment. Theoretical analysis of folding (the "U-to-N" transition), performed at the level of formation and assembly of protein secondary structures, outlines the upper limit of protein folding times (i.e., of the time of search for the most stable fold). Both theories come to essentially the same results; this is not a surprise, because they describe overcoming one and the same free-energy barrier, although the way to the top of this barrier from the side of the unfolded state is very different from the way from the side of the native state; and both theories agree with experiment. In addition, they predict the maximal size of protein domains that fold under solely thermodynamic (rather than kinetic) control and explain the observed maximal size of the "foldable" protein domains.

Fold and flexibility: what can proteins' mechanical properties tell us about their folding nucleus?

The determination of a protein's folding nucleus, i.e. a set of native contacts playing an important role during its folding process, remains an elusive yet essential problem in biochemistry. In this work, we investigate the mechanical properties of 70 protein structures belonging to 14 protein families presenting various folds using coarse-grain Brownian dynamics simulations. The resulting rigidity profiles combined with multiple sequence alignments show that a limited set of rigid residues, which we call the consensus nucleus, occupy conserved positions along the protein sequence. These residues' side chains form a tight interaction network within the protein's core, thus making our consensus nuclei potential folding nuclei. A review of experimental and theoretical literature shows that most (above 80%) of these residues were indeed identified as folding nucleus member in earlier studies.

The outer-membrane export signal of Porphyromonas gingivalis type IX secretion system (T9SS) is a conserved C-terminal β-sandwich domain

In the recently characterized Type IX Secretion System (T9SS), the conserved C-terminal domain (CTD) in secreted proteins functions as an outer membrane translocation signal for export of virulence factors to the cell surface in the Gram-negative Bacteroidetes phylum. In the periodontal pathogen Porphyromonas gingivalis, the CTD is cleaved off by PorU sortase in a sequence-independent manner, and anionic lipopolysaccharide (A-LPS) is attached to many translocated proteins, thus anchoring them to the bacterial surface. Here, we solved the atomic structure of the CTD of gingipain B (RgpB) from P. gingivalis, alone and together with a preceding immunoglobulin-superfamily domain (IgSF). The CTD was found to possess a typical Ig-like fold encompassing seven antiparallel β-strands organized in two β-sheets, packed into a β-sandwich structure that can spontaneously dimerise through C-terminal strand swapping. Small angle X-ray scattering (SAXS) revealed no fixed orientation of the CTD with respect to the IgSF. By introducing insertion or substitution of residues within the inter-domain linker in the native protein, we were able to show that despite the region being unstructured, it nevertheless is resistant to general proteolysis. These data suggest structural motifs located in the two adjacent Ig-like domains dictate the processing of CTDs by the T9SS secretion pathway.

Boosting protein stability with the computational design of β-sheet surfaces

β-sheets often have one face packed against the core of the protein and the other facing solvent. Mutational studies have indicated that the solvent-facing residues can contribute significantly to protein stability, and that the preferred amino acid at each sequence position is dependent on the precise structure of the protein backbone and the identity of the neighboring amino acids. This suggests that the most advantageous methods for designing β-sheet surfaces will be approaches that take into account the multiple energetic factors at play including side chain rotamer preferences, van der Waals forces, electrostatics, and desolvation effects. Here, we show that the protein design software Rosetta, which models these energetic factors, can be used to dramatically increase protein stability by optimizing interactions on the surfaces of small β-sheet proteins. Two design variants of the β-sandwich protein from tenascin were made with 7 and 14 mutations respectively on its β-sheet surfaces. These changes raised the thermal midpoint for unfolding from 45°C to 64°C and 74°C. Additionally, we tested an empirical approach based on increasing the number of potential salt bridges on the surfaces of the β-sheets. This was not a robust strategy for increasing stability, as three of the four variants tested were unfolded.

Transient misfolding dominates multidomain protein folding

Neighbouring domains of multidomain proteins with homologous tandem repeats have divergent sequences, probably as a result of evolutionary pressure to avoid misfolding and aggregation, particularly at the high cellular protein concentrations. Here we combine microfluidic-mixing single-molecule kinetics, ensemble experiments and molecular simulations to investigate how misfolding between the immunoglobulin-like domains of titin is prevented. Surprisingly, we find that during refolding of tandem repeats, independent of sequence identity, more than half of all molecules transiently form a wide range of misfolded conformations. Simulations suggest that a large fraction of these misfolds resemble an intramolecular amyloid-like state reported in computational studies. However, for naturally occurring neighbours with low sequence identity, these transient misfolds disappear much more rapidly than for identical neighbours. We thus propose that evolutionary sequence divergence between domains is required to suppress the population of long-lived, potentially harmful misfolded states, whereas large populations of transient misfolded states appear to be tolerated.

Single molecule kinetics investigations and molecular simulations are useful tools in elucidating protein assembly mechanisms. Here, the authors use these to show that even naturally occurring tandem repeats undergo transient misfolding and that assembly is much more complex than we previously understood.

Structural and dynamic properties that govern the stability of an engineered fibronectin type III domain

Consensus protein design is a rapid and reliable technique for the improvement of protein stability, which relies on the use of homologous protein sequences. To enhance the stability of a fibronectin type III (FN3) domain, consensus design was employed using an alignment of 2123 sequences. The resulting FN3 domain, FN3con, has unprecedented stability, with a melting temperature >100°C, a ΔG_D−N of 15.5 kcal mol⁻¹ and a greatly reduced unfolding rate compared with wild-type. To determine the underlying molecular basis for stability, an X-ray crystal structure of FN3con was determined to 2.0 Å and compared with other FN3 domains of varying stabilities. The structure of FN3con reveals significantly increased salt bridge interactions that are cooperatively networked, and a highly optimized hydrophobic core. Molecular dynamics simulations of FN3con and comparison structures show the cooperative power of electrostatic and hydrophobic networks in improving FN3con stability. Taken together, our data reveal that FN3con stability does not result from a single mechanism, but rather the combination of several features and the removal of non-conserved, unfavorable interactions. The large number of sequences employed in this study has most likely enhanced the robustness of the consensus design, which is now possible due to the increased sequence availability in the post-genomic era. These studies increase our knowledge of the molecular mechanisms that govern stability and demonstrate the rising potential for enhancing stability via the consensus method.

Contribution to the prediction of the fold code: application to immunoglobulin and flavodoxin cases

Background

Folding nucleus of globular proteins formation starts by the mutual interaction of a group of hydrophobic amino acids whose close contacts allow subsequent formation and stability of the 3D structure. These early steps can be predicted by simulation of the folding process through a Monte Carlo (MC) coarse grain model in a discrete space. We previously defined MIRs (Most Interacting Residues), as the set of residues presenting a large number of non-covalent neighbour interactions during such simulation. MIRs are good candidates to define the minimal number of residues giving rise to a given fold instead of another one, although their proportion is rather high, typically [15-20]% of the sequences. Having in mind experiments with two sequences of very high levels of sequence identity (up to 90%) but different folds, we combined the MIR method, which takes sequence as single input, with the “fuzzy oil drop” (FOD) model that requires a 3D structure, in order to estimate the residues coding for the fold. FOD assumes that a globular protein follows an idealised 3D Gaussian distribution of hydrophobicity density, with the maximum in the centre and minima at the surface of the “drop”. If the actual local density of hydrophobicity around a given amino acid is as high as the ideal one, then this amino acid is assigned to the core of the globular protein, and it is assumed to follow the FOD model. Therefore one obtains a distribution of the amino acids of a protein according to their agreement or rejection with the FOD model.

Results

We compared and combined MIR and FOD methods to define the minimal nucleus, or keystone, of two populated folds: immunoglobulin-like (Ig) and flavodoxins (Flav). The combination of these two approaches defines some positions both predicted as a MIR and assigned as accordant with the FOD model. It is shown here that for these two folds, the intersection of the predicted sets of residues significantly differs from random selection. It reduces the number of selected residues by each individual method and allows a reasonable agreement with experimentally determined key residues coding for the particular fold. In addition, the intersection of the two methods significantly increases the specificity of the prediction, providing a robust set of residues that constitute the folding nucleus.

Amino acid distribution rules predict protein fold: protein grammar for beta-strand sandwich-like structures

We present an alternative approach to protein 3D folding prediction based on determination of rules that specify distribution of "favorable" residues, that are mainly responsible for a given fold formation, and "unfavorable" residues, that are incompatible with that fold, in polypeptide sequences. The process of determining favorable and unfavorable residues is iterative. The starting assumptions are based on the general principles of protein structure formation as well as structural features peculiar to a protein fold under investigation. The initial assumptions are tested one-by-one for a set of all known proteins with a given structure. The assumption is accepted as a "rule of amino acid distribution" for the protein fold if it holds true for all, or near all, structures. If the assumption is not accepted as a rule, it can be modified to better fit the data and then tested again in the next step of the iterative search algorithm, or rejected. We determined the set of amino acid distribution rules for a large group of beta sandwich-like proteins characterized by a specific arrangement of strands in two beta sheets. It was shown that this set of rules is highly sensitive (~90%) and very specific (~99%) for identifying sequences of proteins with specified beta sandwich fold structure. The advantage of the proposed approach is that it does not require that query proteins have a high degree of homology to proteins with known structure. So long as the query protein satisfies residue distribution rules, it can be confidently assigned to its respective protein fold. Another advantage of our approach is that it allows for a better understanding of which residues play an essential role in protein fold formation. It may, therefore, facilitate rational protein engineering design.

Mechanistic studies of the biogenesis and folding of outer membrane proteins in vitro and in vivo: what have we learned to date?

Research into the mechanisms by which proteins fold into their native structures has been on-going since the work of Anfinsen in the 1960s. Since that time, the folding mechanisms of small, water-soluble proteins have been well characterised. By contrast, progress in understanding the biogenesis and folding mechanisms of integral membrane proteins has lagged significantly because of the need to create a membrane mimetic environment for folding studies in vitro and the difficulties in finding suitable conditions in which reversible folding can be achieved. Improved knowledge of the factors that promote membrane protein folding and disfavour aggregation now allows studies of folding into lipid bilayers in vitro to be performed. Consequently, mechanistic details and structural information about membrane protein folding are now emerging at an ever increasing pace. Using the panoply of methods developed for studies of the folding of water-soluble proteins. This review summarises current knowledge of the mechanisms of outer membrane protein biogenesis and folding into lipid bilayers in vivo and in vitro and discusses the experimental techniques utilised to gain this information. The emerging knowledge is beginning to allow comparisons to be made between the folding of membrane proteins with current understanding of the mechanisms of folding of water-soluble proteins.

NMR solution structure of the terminal immunoglobulin-like domain from the leptospira host-interacting outer membrane protein, LigB

A number of surface proteins specific to pathogenic strains of Leptospira have been identified. The Lig protein family has shown promise as a marker in typing leptospiral isolates for pathogenesis and as an antigen in vaccines. We used NMR spectroscopy to solve the solution structure of the twelfth immunoglobulin-like (Ig-like) repeat domain from LigB (LigB-12). The fold is similar to that of other bacterial Ig-like domains and comprised mainly of β-strands that form a β-sandwich based on a Greek-key folding arrangement. Based on sequence analysis and conservation of structurally important residues, homology models for the other LigB Ig-like domains were generated. The set of LigB models illustrates the electrostatic differences between the domains as well as the possible interactions between neighboring domains. Understanding the structure of the extracellular portion of LigB and related proteins is important for developing diagnostic methods and new therapeutics directed toward leptospirosis.

The folding of a family of three-helix bundle proteins: spectrin R15 has a robust folding nucleus, unlike its homologous neighbours

Three homologous spectrin domains have remarkably different folding characteristics. We have previously shown that the slow-folding R16 and R17 spectrin domains can be altered to resemble the fast folding R15, in terms of speed of folding (and unfolding), landscape roughness and folding mechanism, simply by substituting five residues in the core. Here we show that, by contrast, R15 cannot be engineered to resemble R16 and R17. It is possible to engineer a slow-folding version of R15, but our analysis shows that this protein neither has a rougher energy landscape nor does change its folding mechanism. Quite remarkably, R15 appears to be a rare example of a protein with a folding nucleus that does not change in position or in size when its folding nucleus is disrupted. Thus, while two members of this protein family are remarkably plastic, the third has apparently a restricted folding landscape.

Antibody variable domain interface and framework sequence requirements for stability and function by high-throughput experiments

Protein structural stability and biological functionality are dictated by the formation of intradomain cores and interdomain interfaces, but the intricate sequence-structure-function interrelationships in the packing of protein cores and interfaces remain difficult to elucidate due to the intractability of enumerating all packing possibilities and assessing the consequences of all the variations. In this work, groups of β strand residues of model antibody variable domains were randomized with saturated mutagenesis and the functional variants were selected for high-throughput sequencing and high-throughput thermal stability measurements. The results show that the sequence preferences of the intradomain hydrophobic core residues are strikingly flexible among hydrophobic residues, implying that these residues are coupled indirectly with antigen binding through energetic stabilization of the protein structures. By contrast, the interdomain interface residues are directly coupled with antigen binding. The interdomain interface should be treated as an integral part of the antigen-binding site.

Loop-sequence features and stability determinants in antibody variable domains by high-throughput experiments

Protein loops are frequently considered as critical determinants in protein structure and function. Recent advances in high-throughput methods for DNA sequencing and thermal stability measurement have enabled effective exploration of sequence-structure-function relationships in local protein regions. Using these data-intensive technologies, we investigated the sequence-structure-function relationships of six complementarity-determining regions (CDRs) and ten non-CDR loops in the variable domains of a model vascular endothelial growth factor (VEGF)-binding single-chain antibody variable fragment (scFv) whose sequence had been optimized via a consensus-sequence approach. The results show that only a handful of residues involving long-range tertiary interactions distant from the antigen-binding site are strongly coupled with antigen binding. This implies that the loops are passive regions in protein folding; the essential sequences of these regions are dictated by conserved tertiary interactions and the consensus local loop-sequence features contribute little to protein stability and function.

Hydrophobic core in domains of immunoglobulin-like fold

This work analyzes proteins which contain an immunoglobulin fold, focusing on their hydrophobic core structure. The "fuzzy oil drop" model was used to measure the regularity of hydrophobicity distribution in globular domains belonging to proteins which exhibit the above-mentioned fold. Light-chain IgG domains are found to frequently contain regular hydrophobic cores, unlike the corresponding heavy-chain domains. Enzymes and DNA binding proteins present in the data-set are found to exhibit poor accordance with the hydrophobic core model.

The loop hypothesis: contribution of early formed specific non-local interactions to the determination of protein folding pathways

The extremely fast and efficient folding transition (in seconds) of globular proteins led to the search for some unifying principles embedded in the physics of the folding polypeptides. Most of the proposed mechanisms highlight the role of local interactions that stabilize secondary structure elements or a folding nucleus as the starting point of the folding pathways, i.e., a "bottom-up" mechanism. Non-local interactions were assumed either to stabilize the nucleus or lead to the later steps of coalescence of the secondary structure elements. An alternative mechanism was proposed, an "up-down" mechanism in which it was assumed that folding starts with the formation of very few non-local interactions which form closed long loops at the initiation of folding. The possible biological advantage of this mechanism, the "loop hypothesis", is that the hydrophobic collapse is associated with ordered compactization which reduces the chance for degradation and misfolding. In the present review the experiments, simulations and theoretical consideration that either directly or indirectly support this mechanism are summarized. It is argued that experiments monitoring the time-dependent development of the formation of specifically targeted early-formed sub-domain structural elements, either long loops or secondary structure elements, are necessary. This can be achieved by the time-resolved FRET-based "double kinetics" method in combination with mutational studies. Yet, attempts to improve the time resolution of the folding initiation should be extended down to the sub-microsecond time regime in order to design experiments that would resolve the classes of proteins which first fold by local or non-local interactions.

Interplay between drying and stability of a TIM barrel protein: a combined simulation-experimental study

Recent molecular dynamics simulations have suggested important roles for nanoscale dewetting in the stability, function, and folding dynamics of proteins. Using a synergistic simulation-experimental approach on the αTS TIM barrel protein, we validated this hypothesis by revealing the occurrence of drying inside hydrophobic amino acid clusters and its manifestation in experimental measures of protein stability and structure. Cavities created within three clusters of branched aliphatic amino acids [isoleucine, leucine, and valine (ILV) clusters] were found to experience strong water density fluctuations or intermittent dewetting transitions in simulations. Individually substituting 10 residues in the large ILV cluster at the N-terminus with less hydrophobic alanines showed a weakening or diminishing effect on dewetting that depended on the site of the mutation. Our simulations also demonstrated that replacement of buried leucines with isosteric, polar asparagines enhanced the wetting of the N- and C-terminal clusters. The experimental results on the stability, secondary structure, and compactness of the native and intermediate states for the asparagine variants are consistent with the preferential drying of the large N-terminal cluster in the intermediate. By contrast, the region encompassing the small C-terminal cluster experiences only partial drying in the intermediate, and its structure and stability are unaffected by the asparagine substitution. Surprisingly, the structural distortions required to accommodate the replacement of leucine by asparagine in the N-terminal cluster revealed the existence of alternative stable folds in the native basin. This combined simulation-experimental study demonstrates the critical role of drying within hydrophobic ILV clusters in the folding and stability of the αTS TIM barrel.

Take home lessons from studies of related proteins

The 'Fold Approach' involves a detailed analysis of the folding of several topologically, structurally and/or evolutionarily related proteins. Such studies can reveal determinants of the folding mechanism beyond the gross topology, and can dissect the residues required for folding from those required for stability or function. While this approach has not yet matured to the point where we can predict the native conformation of any polypeptide chain in silico, it has been able to highlight, amongst others, the specific residues that are responsible for nucleation, pathway malleability, kinetic intermediates, chain knotting, internal friction and Paracelsus switches. Some of the most interesting discoveries have resulted from the attempt to explain differences between homologues.

Folding mechanism of an extremely thermostable (βα)(8)-barrel enzyme: a high kinetic barrier protects the protein from denaturation

HisF, the cyclase subunit of imidazole glycerol phosphate synthase (ImGPS) from Thermotoga maritima, is an extremely thermostable (βα)(8)-barrel protein. We elucidated the unfolding and refolding mechanism of HisF. Its unfolding transition is reversible and adequately described by the two-state model, but 6 weeks is necessary to reach equilibrium (at 25 °C). During refolding, initially a burst-phase off-pathway intermediate is formed. The subsequent productive folding occurs in two kinetic phases with time constants of ~3 and ~20 s. They reflect a sequential process via an on-pathway intermediate, as revealed by stopped-flow double-mixing experiments. The final step leads to native HisF, which associates with the glutaminase subunit HisH to form the functional ImGPS complex. The conversion of the on-pathway intermediate to the native protein results in a 10(6)-fold increase of the time constant for unfolding from 89 ms to 35 h (at 4.0 M GdmCl) and thus establishes a high energy barrier to denaturation. We conclude that the extra stability of HisF is used for kinetic protection against unfolding. In its refolding mechanism, HisF resembles other (βα)(8)-barrel proteins.

Big domains are novel Ca²+-binding modules: evidences from big domains of Leptospira immunoglobulin-like (Lig) proteins

BACKGROUND

Many bacterial surface exposed proteins mediate the host-pathogen interaction more effectively in the presence of Ca²+. Leptospiral immunoglobulin-like (Lig) proteins, LigA and LigB, are surface exposed proteins containing Bacterial immunoglobulin like (Big) domains. The function of proteins which contain Big fold is not known. Based on the possible similarities of immunoglobulin and βγ-crystallin folds, we here explore the important question whether Ca²+ binds to a Big domains, which would provide a novel functional role of the proteins containing Big fold.

PRINCIPAL FINDINGS

We selected six individual Big domains for this study (three from the conserved part of LigA and LigB, denoted as Lig A3, Lig A4, and LigBCon5; two from the variable region of LigA, i.e., 9(th) (Lig A9) and 10(th) repeats (Lig A10); and one from the variable region of LigB, i.e., LigBCen2. We have also studied the conserved region covering the three and six repeats (LigBCon1-3 and LigCon). All these proteins bind the calcium-mimic dye Stains-all. All the selected four domains bind Ca²+ with dissociation constants of 2-4 µM. Lig A9 and Lig A10 domains fold well with moderate thermal stability, have β-sheet conformation and form homodimers. Fluorescence spectra of Big domains show a specific doublet (at 317 and 330 nm), probably due to Trp interaction with a Phe residue. Equilibrium unfolding of selected Big domains is similar and follows a two-state model, suggesting the similarity in their fold.

CONCLUSIONS

We demonstrate that the Lig are Ca²+-binding proteins, with Big domains harbouring the binding motif. We conclude that despite differences in sequence, a Big motif binds Ca²+. This work thus sets up a strong possibility for classifying the proteins containing Big domains as a novel family of Ca²+-binding proteins. Since Big domain is a part of many proteins in bacterial kingdom, we suggest a possible function these proteins via Ca²+ binding.

What lessons can be learned from studying the folding of homologous proteins?

The studies of the folding of structurally related proteins have proved to be a very important tool for investigating protein folding. Here we review some of the insights that have been gained from such studies. Our highlighted studies show just how such an investigation should be designed and emphasise the importance of the synergy between experiment and theory. We also stress the importance of choosing the right system carefully, exploiting the excellent structural and sequence databases at our disposal.

Topological frustration in beta alpha-repeat proteins: sequence diversity modulates the conserved folding mechanisms of alpha/beta/alpha sandwich proteins

The thermodynamic hypothesis of Anfinsen postulates that structures and stabilities of globular proteins are determined by their amino acid sequences. Chain topology, however, is known to influence the folding reaction, in that motifs with a preponderance of local interactions typically fold more rapidly than those with a larger fraction of nonlocal interactions. Together, the topology and sequence can modulate the energy landscape and influence the rate at which the protein folds to the native conformation. To explore the relationship of sequence and topology in the folding of beta alpha-repeat proteins, which are dominated by local interactions, we performed a combined experimental and simulation analysis on two members of the flavodoxin-like, alpha/beta/alpha sandwich fold. Spo0F and the N-terminal receiver domain of NtrC (NT-NtrC) have similar topologies but low sequence identity, enabling a test of the effects of sequence on folding. Experimental results demonstrated that both response-regulator proteins fold via parallel channels through highly structured submillisecond intermediates before accessing their cis prolyl peptide bond-containing native conformations. Global analysis of the experimental results preferentially places these intermediates off the productive folding pathway. Sequence-sensitive Gō-model simulations conclude that frustration in the folding in Spo0F, corresponding to the appearance of the off-pathway intermediate, reflects competition for intra-subdomain van der Waals contacts between its N- and C-terminal subdomains. The extent of transient, premature structure appears to correlate with the number of isoleucine, leucine, and valine (ILV) side chains that form a large sequence-local cluster involving the central beta-sheet and helices alpha2, alpha 3, and alpha 4. The failure to detect the off-pathway species in the simulations of NT-NtrC may reflect the reduced number of ILV side chains in its corresponding hydrophobic cluster. The location of the hydrophobic clusters in the structure may also be related to the differing functional properties of these response regulators. Comparison with the results of previous experimental and simulation analyses on the homologous CheY argues that prematurely folded unproductive intermediates are a common property of the beta alpha-repeat motif.

How antibodies fold

B cells use unconventional strategies for the production of a seemingly unlimited number of antibodies from a very limited amount of DNA. These methods dramatically increase the likelihood of producing proteins that cannot fold or assemble appropriately. B cells are therefore particularly dependent on 'quality control' mechanisms to oversee antibody production. Recent in vitro experiments demonstrate that Ig domains have evolved diverse folding strategies ranging from robust spontaneous folding to intrinsically disordered domains that require assembly with their partner domains to fold; in vivo experiments reveal that these different folding characteristics form the basis for cellular checkpoints in Ig transport. Taken together, these reports provide a detailed understanding of how B cells monitor and ensure the functional fidelity of Ig proteins.

The folding pathway of the antibody V(L) domain

Antibodies are modular proteins consisting of domains that exhibit a beta-sandwich structure, the so-called immunoglobulin fold. Despite structural similarity, differences in folding and stability exist between different domains. In particular, the variable domain of the light chain V(L) is unusual as it is associated with misfolding diseases, including the pathologic assembly of the protein into fibrillar structures. Here, we have analysed the folding pathway of a V(L) domain with a view to determine features that may influence the relationship between productive folding and fibril formation. The V(L) domain from MAK33 (murine monoclonal antibody of the subtype kappa/IgG1) has not previously been associated with fibrillisation but is shown here to be capable of forming fibrils. The folding pathway of this V(L) domain is complex, involving two intermediates in different pathways. An obligatory early molten globule-like intermediate with secondary structure but only loose tertiary interactions is inferred. The native state can then be formed directly from this intermediate in a phase that can be accelerated by the addition of prolyl isomerases. However, an alternative pathway involving a second, more native-like intermediate is also significantly populated. Thus, the protein can reach the native state via two distinct folding pathways. Comparisons to the folding pathways of other antibody domains reveal similarities in the folding pathways; however, in detail, the folding of the V(L) domain is striking, with two intermediates populated on different branches of the folding pathway, one of which could provide an entry point for molecules diverted into the amyloid pathway.

Prediction of the protein folding core: application to the immunoglobulin fold

We propose an algorithm that allows predicting residues important for the formation of the structure of globular proteins. It relies on a simulation that detects the amino acids presenting a maximum number of neighbours during the early steps of the folding process. They have been called MIR (Most Interacting Residues). Independently, description of the protein structures in fragments with closed ends shows the correlation between these extremities and the core of the globules. These fragments are of rather constant length, typically between 20 and 25 amino acids, and we have previously shown that their extremities are preferentially occupied by MIR. Introduction of rules derived from this fragment analysis of tertiary structures allows to smooth the distribution of MIR, for a better match between TEF ends and MIR. In order to assess this prediction of the folding core, a large family of structures has been used, with sequences as different as possible. A dataset of 56 immunoglobulin structures of various functions but common fold has been used in this study. This fold was chosen because it is one of the most populated with a large amount of data available on its nucleus. In the immunoglobulin domain, "functional and structural load is clearly separated: loops are responsible for binding and recognition while interactions between several residues of the buried core provide stability and fast folding"[1]. We then determined the positions susceptible of high importance for the folding process to occur and compared them to published data, either to High Throw Out Order (HTOO), Conservatism of Conservatism (CoC) or Phi value experiments. It results a reasonable agreement between the positions that we predict and experimental data. Besides, our prediction goes beyond the simple use of a null solvent accessibility of amino acids as a criterion to predict the core. We find the same quality of our prediction on the flavodoxin like superfamily.

The folding mechanism of BBL: Plasticity of transition-state structure observed within an ultrafast folding protein family

Studies on members of protein families with similar structures but divergent sequences provide insights into the effects of sequence composition on the mechanism of folding. Members of the peripheral subunit-binding domain (PSBD) family fold ultrafast and approach the smallest size for cooperatively folding proteins. Phi-Value analysis of the PSBDs E3BD and POB reveals folding via nucleation-condensation through structurally very similar, polarized transition states. Here, we present a Phi-value analysis of the family member BBL and found that it also folds by a nucleation-condensation mechanism. The mean Phi values of BBL, E3BD, and POB were near identical, indicating similar fractions of non-covalent interactions being formed in the transition state. Despite the overall conservation of folding mechanism in this protein family, however, the pattern of Phi values determined for BBL revealed a larger dispersion of the folding nucleus across the entire structure, and the transition state was less polarized. The observed plasticity of transition-state structure can be rationalized by the different helix-forming propensities of PSBD sequences. The very strong helix propensity in the first helix of BBL, relative to E3BD and POB, appears to recruit more structure formation in that helix in the transition state at the expense of weaker interactions in the second helix. Differences in sequence composition can modulate transition-state structure of even the smallest natural protein domains.

Topology is the principal determinant in the folding of a complex all-alpha Greek key death domain from human FADD

In order to elucidate the relative importance of secondary structure and topology in determining folding mechanism, we have carried out a phi-value analysis of the death domain (DD) from human FADD. FADD DD is a 100 amino acid domain consisting of six anti-parallel alpha helices arranged in a Greek key structure. We asked how does the folding of this domain compare with that of (a) other all-alpha-helical proteins and (b) other Greek key proteins? Is the folding pathway determined mainly by secondary structure or is topology the principal determinant? Our Φ-value analysis reveals a striking resemblance to the all-beta Greek key immunoglobulin-like domains. Both fold via diffuse transition states and, importantly, long-range interactions between the four central elements of secondary structure are established in the transition state. The elements of secondary structure that are less tightly associated with the central core are less well packed in both cases. Topology appears to be the dominant factor in determining the pathway of folding in all Greek key domains.

Crystal structure of the C-terminal domain of a flagellar hook-capping protein from Xanthomonas campestris

The crystal structure of the C-terminal domain of a hook-capping protein FlgD from the plant pathogen Xanthomonas campestris (Xc) has been determined to a resolution of ca 2.5 A using X-ray crystallography. The monomer of whole FlgD comprises 221 amino acids with a molecular mass of 22.7 kDa, but the flexible N-terminus is cleaved for up to 75 residues during crystallization. The final structure of the C-terminal domain reveals a novel hybrid comprising a tudor-like domain interdigitated with a fibronectin type III domain. The C-terminal domain of XcFlgD forms three types of dimers in the crystal. In agreement with this, analytical ultracentrifugation and gel filtration experiments reveal that they form a stable dimer in solution. From these results, we propose that the Xc flagellar hook cap protein FlgD comprises two individual domains, a flexible N-terminal domain that cannot be detected in the current study and a stable C-terminal domain that forms a stable dimer.

Folding of a LysM domain: entropy-enthalpy compensation in the transition state of an ideal two-state folder

Protein-engineering methods (Phi-values) were used to investigate the folding transition state of a lysin motif (LysM) domain from Escherichia coli membrane-bound lytic murein transglycosylase D. This domain consists of just 48 structured residues in a symmetrical betaalphaalphabeta arrangement and is the smallest alphabeta protein yet investigated using these methods. An extensive mutational analysis revealed a highly robust folding pathway with no detectable transition state plasticity, indicating that LysM is an example of an ideal two-state folder. The pattern of Phi-values denotes a highly polarised transition state, with significant formation of the helices but no structure within the beta-sheet. Remarkably, this transition state remains polarised after circularisation of the domain, and exhibits an identical Phi-value pattern; however, the interactions within the transition state are uniformly weaker in the circular variant. This observation is supported by results from an Eyring analysis of the folding rates of the two proteins. We propose that the folding pathway of LysM is dominated by enthalpic rather than entropic considerations, and suggest that the lower entropy cost of formation of the circular transition state is balanced, to some extent, by the lower enthalpy of contacts within this structure.