Effects of core mutations on the folding of a beta-sheet protein: implications for backbone organization in the I-state

Accurate prediction of protein tertiary structural changes induced by single-site mutations with equivariant graph neural networks

Predicting the change of protein tertiary structure caused by singlesite mutations is important for studying protein structure, function, and interaction. Even though computational protein structure prediction methods such as AlphaFold can predict the overall tertiary structures of most proteins rather accurately, they are not sensitive enough to accurately predict the structural changes induced by single-site amino acid mutations on proteins. Specialized mutation prediction methods mostly focus on predicting the overall stability or function changes caused by mutations without attempting to predict the exact mutation-induced structural changes, limiting their use in protein mutation study. In this work, we develop the first deep learning method based on equivariant graph neural networks (EGNN) to directly predict the tertiary structural changes caused by single-site mutations and the tertiary structure of any protein mutant from the structure of its wild-type counterpart. The results show that it performs substantially better in predicting the tertiary structures of protein mutants than the widely used protein structure prediction method AlphaFold.

Repair of topoisomerase 1-induced DNA damage by tyrosyl-DNA phosphodiesterase 2 (TDP2) is dependent on its magnesium binding

Topoisomerases are enzymes that relax DNA supercoiling during replication and transcription. Camptothecin, a topoisomerase 1 (TOP1) inhibitor, and its analogs trap TOP1 at the 3'-end of DNA as a DNA-bound intermediate, resulting in DNA damage that can kill cells. Drugs with this mechanism of action are widely used to treat cancers. It has previously been shown that tyrosyl-DNA phosphodiesterase 1 (TDP1) repairs TOP1-induced DNA damage generated by camptothecin. In addition, tyrosyl-DNA phosphodiesterase 2 (TDP2) plays critical roles in repairing topoisomerase 2 (TOP2)-induced DNA damage at the 5'-end of DNA and in promoting the repair of TOP1-induced DNA damage in the absence of TDP1. However, the catalytic mechanism by which TDP2 processes TOP1-induced DNA damage has not been elucidated. In this study, we found that a similar catalytic mechanism underlies the repair of TOP1- and TOP2-induced DNA damage by TDP2, with Mg2+-TDP2 binding playing a role in both repair mechanisms. We show chain-terminating nucleoside analogs are incorporated into DNA at the 3'-end and abort DNA replication to kill cells. Furthermore, we found that Mg2+-TDP2 binding also contributes to the repair of incorporated chain-terminating nucleoside analogs. Overall, these findings reveal the role played by Mg2+-TDP2 binding in the repair of both 3'- and 5'-blocking DNA damage.

How a single mutation alters the protein structure: a simulation investigation on protein tyrosine phosphatase SHP2

Protein tyrosine phosphatase SHP2 is a key regulator modulating several signaling pathways. The oncogenic mutation E76K in SHP2 releases the enzyme from an autoinhibited, closed conformation into an active, open conformation. Here, we investigated the conformational dynamics of SHP2 and the effect of the E76K mutation on its conformational ensemble via extensive molecular dynamics (MD) and metadynamics (MetaD) simulations. Our simulations provide atomistic details on how the E76K mutated SHP2 prefers the open state and also reveal that the transition between the closed and the open states is highly collective. Several intermediate metastable states during the conformational transition between the closed and the open states were also investigated. Understanding how the single E76K mutation induces the conformational change in SHP2 could facilitate the further design of SHP2 inhibitors.

Molecular insights into the interaction of HPV-16 E6 variants against MAGI-1 PDZ1 domain

Oncogenic protein E6 from Human Papilloma Virus 16 (HPV-16) mediates the degradation of Membrane-associated guanylate kinase with inverted domain structure-1 (MAGI-1), throughout the interaction of its protein binding motif (PBM) with the Discs-large homologous regions 1 (PDZ1) domain of MAG1-1. Generic variation in the E6 gene that translates to changes in the protein’s amino acidic sequence modifies the interaction of E6 with the cellular protein MAGI-1. MAGI-1 is a scaffolding protein found at tight junctions of epithelial cells, where it interacts with a variety of proteins regulating signaling pathways. MAGI-1 is a multidomain protein containing two WW (rsp-domain-9), one guanylate kinase-like, and six PDZ domains. PDZ domains played an important role in the function of MAGI-1 and served as targets for several viral proteins including the HPV-16 E6. The aim of this work was to evaluate, with an in silico approach, employing molecular dynamics simulation and protein–protein docking, the interaction of the intragenic variants E-G350 (L83V), E-C188/G350 (E29Q/L83V), E-A176/G350 (D25N/L83V), E6-AAa (Q14H/H78Y/83V) y E6-AAc (Q14H/I27RH78Y/L83V) and E6-reference of HPV-16 with MAGI-1. We found that variants E-G350, E-C188/G350, E-A176/G350, AAa and AAc increase their affinity to our two models of MAGI-1 compared to E6-reference.

In silico characterization of mutations circulating in SARS-CoV-2 structural proteins

SARS-CoV-2 has recently emerged as a pandemic that has caused more than 2.4 million deaths worldwide. Since the onset of infections, several full-length sequences of viral genome have been made available which have been used to gain insights into viral dynamics. We utilised a meta-data driven comparative analysis tool for sequences (Meta-CATS) algorithm to identify mutations in 829 SARS-CoV-2 genomes from around the world. The algorithm predicted sixty-one mutations among SARS-CoV-2 genomes. We observed that most of the mutations were concentrated around three protein coding genes viz nsp3 (non-structural protein 3), RdRp (RNA-directed RNA polymerase) and Nucleocapsid (N) proteins of SARS-CoV-2. We used various computational tools including normal mode analysis (NMA), C-α discrete molecular dynamics (DMD) and all-atom molecular dynamic simulations (MD) to study the effect of mutations on functionality, stability and flexibility of SARS-CoV-2 structural proteins including envelope (E), N and spike (S) proteins. PredictSNP predictor suggested that four mutations (L37H in E, R203K and P344S in N and D614G in S) out of seven were predicted to be neutral whilst the remaining ones (P13L, S197L and G204R in N) were predicted to be deleterious in nature thereby impacting protein functionality. NMA, C-α DMD and all-atom MD suggested some mutations to have stabilizing roles (P13L, S197L and R203K in N protein) where remaining ones were predicted to destabilize mutant protein. In summary, we identified significant mutations in SARS-CoV-2 genomes as well as used computational approaches to further characterize the possible effect of highly significant mutations on SARS-CoV-2 structural proteins.

Communicated by Ramaswamy H. Sarma

Influence of Disease-Causing Mutations on Protein Structural Networks

The interactions between residues in a protein tertiary structure can be studied effectively using the approach of protein structure network (PSN). A PSN is a node-edge representation of the structure with nodes representing residues and interactions between residues represented by edges. In this study, we have employed weighted PSNs to understand the influence of disease-causing mutations on proteins of known 3D structures. We have used manually curated information on disease mutations from UniProtKB/Swiss-Prot and their corresponding protein structures of wildtype and disease variant from the protein data bank. The PSNs of the wildtype and disease-causing mutant are compared to analyse variation of global and local dissimilarity in the overall network and at specific sites. We study how a mutation at a given site can affect the structural network at a distant site which may be involved in the function of the protein. We have discussed specific examples of the disease cases where the protein structure undergoes limited structural divergence in their backbone but have large dissimilarity in their all atom networks and vice versa, wherein large conformational alterations are observed while retaining overall network. We analyse the effect of variation of network parameters that characterize alteration of function or stability.

ProThermDB: thermodynamic database for proteins and mutants revisited after 15 years

ProThermDB is an updated version of the thermodynamic database for proteins and mutants (ProTherm), which has ∼31 500 data on protein stability, an increase of 84% from the previous version. It contains several thermodynamic parameters such as melting temperature, free energy obtained with thermal and denaturant denaturation, enthalpy change and heat capacity change along with experimental methods and conditions, sequence, structure and literature information. Besides, the current version of the database includes about 120 000 thermodynamic data obtained for different organisms and cell lines, which are determined by recent high throughput proteomics techniques using whole-cell approaches. In addition, we provided a graphical interface for visualization of mutations at sequence and structure levels. ProThermDB is cross-linked with other relevant databases, PDB, UniProt, PubMed etc. It is freely available at https://web.iitm.ac.in/bioinfo2/prothermdb/index.html without any login requirements. It is implemented in Python, HTML and JavaScript, and supports the latest versions of major browsers, such as Firefox, Chrome and Safari.

Functional impact of cancer patient-associated Bcl-xL mutations

Bcl-xL, an antiapoptotic protein, is frequently overexpressed in cancer to promote survival of tumor cells. However, we have previously shown that Bcl-xL promotes migration, invasion, and metastasis independent of its antiapoptotic function in mitochondria. The pro-metastatic function of Bcl-xL may require its translocation into the nucleus. Besides overexpression, patient-associated mutations of Bcl-xL have been identified in large-scale cancer genomics projects. Understanding the functions of these mutations will guide the development of precision medicine. Here, we selected four patient-associated Bcl-xL mutations, R132W, N136K, R165W, and A201T, to investigate their impacts on antiapoptosis, migration, and nuclear translocation. We found that all four mutation proteins could be detected in both the nucleus and cytosol. Although all four mutations disrupted the antiapoptosis function, one of these mutants, N136K, significantly improved the ability to promote cell migration. These data suggest the importance of developing novel Bcl-xL inhibitors to ablate both antiapoptotic and pro-metastatic functions of Bcl-xL in cancer.

Conserved tryptophan mutation disrupts structure and function of immunoglobulin domain revealing unusual tyrosine fluorescence

Immunoglobulin (Ig) domains are the most prevalent protein domain structure and share a highly conserved folding pattern; however, this structural family of proteins is also the most diverse in terms of biological roles and tissue expression. Ig domains vary significantly in amino acid sequence but share a highly conserved tryptophan in the hydrophobic core of this beta-stranded protein. Palladin is an actin binding and bundling protein that has five Ig domains and plays an important role in normal cell adhesion and motility. Mutation of the core tryptophan in one Ig domain of palladin has been identified in a pancreatic cancer cell line, suggesting a crucial role for this sole tryptophan in palladin Ig domain structure, stability, and function. We found that actin binding and bundling was not completely abolished with removal of this tryptophan despite a partially unfolded structure and significantly reduced stability of the mutant Ig domain as shown by circular dichroism investigations. In addition, this mutant palladin domain displays a tryptophan-like fluorescence attributed to an anomalous tyrosine emission at 341 nm. Our results indicate that this emission originates from a tyrosinate that may be formed in the excited ground state by proton transfer to a nearby aspartic acid residue. Furthermore, this study emphasizes the importance of tryptophan in protein structural stability and illustrates how tyrosinate emission contributions may be overlooked during the interpretation of the fluorescence properties of proteins.

Mutation effect estimation on protein-protein interactions using deep contextualized representation learning

The functional impact of protein mutations is reflected on the alteration of conformation and thermodynamics of protein–protein interactions (PPIs). Quantifying the changes of two interacting proteins upon mutations is commonly carried out by computational approaches. Hence, extensive research efforts have been put to the extraction of energetic or structural features on proteins, followed by statistical learning methods to estimate the effects of mutations on PPI properties. Nonetheless, such features require extensive human labors and expert knowledge to obtain, and have limited abilities to reflect point mutations. We present an end-to-end deep learning framework, MuPIPR (Mutation Effects in Protein–protein Interaction PRediction Using Contextualized Representations), to estimate the effects of mutations on PPIs. MuPIPR incorporates a contextualized representation mechanism of amino acids to propagate the effects of a point mutation to surrounding amino acid representations, therefore amplifying the subtle change in a long protein sequence. On top of that, MuPIPR leverages a Siamese residual recurrent convolutional neural encoder to encode a wild-type protein pair and its mutation pair. Multi-layer perceptron regressors are applied to the protein pair representations to predict the quantifiable changes of PPI properties upon mutations. Experimental evaluations show that, with only sequence information, MuPIPR outperforms various state-of-the-art systems on estimating the changes of binding affinity for SKEMPI v1, and offers comparable performance on SKEMPI v2. Meanwhile, MuPIPR also demonstrates state-of-the-art performance on estimating the changes of buried surface areas. The software implementation is available at https://github.com/guangyu-zhou/MuPIPR.

Recent Advances on the Structure and Function of RNA Acetyltransferase Kre33/NAT10

Ribosome biogenesis is one of the most energy demanding processes in the cell. In eukaryotes, the main steps of this process occur in the nucleolus and include pre-ribosomal RNA (pre-rRNA) processing, post-transcriptional modifications, and assembly of many non-ribosomal factors and ribosomal proteins in order to form mature and functional ribosomes. In yeast and humans, the nucleolar RNA acetyltransferase Kre33/NAT10 participates in different maturation events, such as acetylation and processing of 18S rRNA, and assembly of the 40S ribosomal subunit. Here, we review the structural and functional features of Kre33/NAT10 RNA acetyltransferase, and we underscore the importance of this enzyme in ribosome biogenesis, as well as in acetylation of non-ribosomal targets. We also report on the role of human NAT10 in Hutchinson-Gilford progeria syndrome.

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.

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.

Predicting the functional impact of protein mutations: application to cancer genomics

As large-scale re-sequencing of genomes reveals many protein mutations, especially in human cancer tissues, prediction of their likely functional impact becomes important practical goal. Here, we introduce a new functional impact score (FIS) for amino acid residue changes using evolutionary conservation patterns. The information in these patterns is derived from aligned families and sub-families of sequence homologs within and between species using combinatorial entropy formalism. The score performs well on a large set of human protein mutations in separating disease-associated variants (∼19 200), assumed to be strongly functional, from common polymorphisms (∼35 600), assumed to be weakly functional (area under the receiver operating characteristic curve of ∼0.86). In cancer, using recurrence, multiplicity and annotation for ∼10 000 mutations in the COSMIC database, the method does well in assigning higher scores to more likely functional mutations ('drivers'). To guide experimental prioritization, we report a list of about 1000 top human cancer genes frequently mutated in one or more cancer types ranked by likely functional impact; and, an additional 1000 candidate cancer genes with rare but likely functional mutations. In addition, we estimate that at least 5% of cancer-relevant mutations involve switch of function, rather than simply loss or gain of function.

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.

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.

Crosstalk between the protein surface and hydrophobic core in a core-swapped fibronectin type III domain

Two homologous fibronectin type III (fnIII) domains, FNfn10 (the 10th fnIII domain of human fibronectin) and TNfn3 (the third fnIII domain of human tenascin), have essentially the same backbone structure, although they share only ∼ 24% sequence identity. While they share a similar folding mechanism with a common core of key residues in the folding transition state, they differ in many other physical properties. We use a chimeric protein, FNoTNc, to investigate the molecular basis for these differences. FNoTNc is a core-swapped protein, containing the “outside” (surface and loops) of FNfn10 and the hydrophobic core of TNfn3. Remarkably, FNoTNc retains the structure of the parent proteins despite the extent of redesign, allowing us to gain insight into which components of each parent protein are responsible for different aspects of its behaviour. Naively, one would expect properties that appear to depend principally on the core to be similar to TNfn3, for example, the response to mutations, folding kinetics and side-chain dynamics, while properties apparently determined by differences in the surface and loops, such as backbone dynamics, would be more like FNfn10. While this is broadly true, it is clear that there are also unexpected crosstalk effects between the core and the surface. For example, the anomalous response of FNfn10 to mutation is not solely a property of the core as we had previously suggested.

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

A number of β-sandwich immunoglobulin-like domains have been shown to fold using a set of structurally equivalent residues that form a folding nucleus deep within the core of the protein. Formation of this nucleus is sufficient to establish the complex Greek key topology of the native state. These nucleating residues are highly conserved within the immunoglobulin superfamily, but are less well conserved in the fibronectin type III (fnIII) superfamily, where the requirement is simply to have four interacting hydrophobic residues. However, there are rare examples where this nucleation pattern is absent. In this study, we have investigated the folding of a novel member of the fnIII superfamily whose nucleus appears to lack one of the four buried hydrophobic residues. We show that the folding mechanism is unaltered, but the folding nucleus has moved within the hydrophobic core.

Protein folding kinetics provides a context-independent assessment of beta-strand propensity in the Fyn SH3 domain

Structural database-derived propensities for amino acids to adopt particular local protein structures, such as alpha-helix and beta-strand, have long been recognized and effectively exploited for the prediction of protein secondary structure. However, the experimental verification of database-derived propensities using mutagenesis studies has been problematic, especially for beta-strand propensities, because local structural preferences are often confounded by non-local interactions arising from formation of the native tertiary structure. Thus, the overall thermodynamic stability of a protein is not always altered in a predictable manner by changes in local structural propensity at a single position. In this study, we have undertaken an investigation of the relationship between beta-strand propensity and protein folding kinetics. By characterizing the effects of a wide variety of amino acid substitutions at two different beta-strand positions in an SH3 domain, we have found that the observed changes in protein folding rates are very well correlated to beta-strand propensities for almost all of the substitutions examined. In contrast, there is little correlation between propensities and unfolding rates. These data indicate that beta-strand conformation is well formed in the structured portion of the SH3 domain transition state, and that local structure propensity strongly influences the stability of the transition state. Since the transition state is known to be packed more loosely than the native state and likely lacks many of the non-local stabilizing interactions seen in the native state, we suggest that folding kinetics studies may generally provide an effective means for the experimental validation of database-derived local structural propensities.

The folding and evolution of multidomain proteins

Analyses of genomes show that more than 70% of eukaryotic proteins are composed of multiple domains. However, most studies of protein folding focus on individual domains and do not consider how interactions between domains might affect folding. Here, we address this by analysing the three-dimensional structures of multidomain proteins that have been characterized experimentally and observe that where the interface is small and loosely packed, or unstructured, the folding of the domains is independent. Furthermore, recent studies indicate that multidomain proteins have evolved mechanisms to minimize the problems of interdomain misfolding.

Folding of Cu/Zn superoxide dismutase suggests structural hotspots for gain of neurotoxic function in ALS: parallels to precursors in amyloid disease

Amyotrophic lateral sclerosis (ALS) is a neurodegenerative disease linked to misfolding of the ubiquitous enzyme Cu/Zn superoxide dismutase (SOD). In contrast to other protein-misfolding disorders with similar neuropathogenesis, ALS is not always associated with the in vivo deposition of protein aggregates. Thus, under the assumption that all protein-misfolding disorders share at primary level a similar disease mechanism, ALS constitutes an interesting disease model for identifying the yet-mysterious precursor states from which the cytotoxic pathway emerges. In this study, we have mapped out the conformational repertoire of the apoSOD monomer through analysis of its folding behavior. The results allow us to target the regions of the SOD structure that are most susceptible to unfolding locally under physiological conditions, leading to the exposure of structurally promiscuous interfaces that are normally hidden in the protein's interior. The structure of this putative ALS precursor is strikingly similar to those implicated in amyloid disease.

Phi-analysis of the folding of the B domain of protein A using multiple optical probes

We examined the co-operativity of ultra-fast folding of a protein and whether the Phi-value analysis of its transition state depended on the location of the optical probe. We incorporated in turn a tryptophan residue into each of the three helices of the B domain of Protein A. Each Trp mutant of the three-helix bundle protein was used as a pseudo-wild-type parent for Phi-analysis in which the intrinsic Trp fluorescence probed the formation of each helix during the transition state. Apart from local effects in the immediate vicinity of the probe, the three separate sets of Phi-values were in excellent agreement, demonstrating the overall co-operativity of folding and the robustness of the Phi-analysis. The transition state of folding of Protein A contains the second helix being well formed with many stabilizing tertiary hydrophobic interactions. In contrast, the first and the third helices are more poorly structured in the transition state. The mechanism of folding thus involves the concurrent formation of secondary and tertiary interactions, and is towards the nucleation-condensation extreme in the nucleation-condensation-framework continuum of mechanism, with helix 2 being the nucleus. We provide an error analysis of Phi-values derived purely from the kinetics of two-state chevron plots.

Dependence of the size of the initially collapsed form during the refolding of barstar on denaturant concentration: evidence for a continuous transition

Two-site fluorescence resonance energy transfer (FRET) measurements have been made to determine how two intra-molecular distances contract in the sub-millisecond collapse reaction that occurs initially during the refolding of the small protein barstar. FRET measurements were made on two, single-Cys and single-Trp-containing mutant forms of barstar, Cys25 and Cys62, in each of which a thionitrobenzoate (TNB) adduct was attached to the cysteine thiol. In each protein, the core tryptophan, Trp53, acted as the FRET donor, and the TNB adduct, located either at C25 or at C62, acted as the FRET acceptor. The stabilities as well as observable folding kinetics of the Cys25 and Cys62 mutant proteins were found to be identical. The presence of the TNB adduct on the cysteine did not alter the stability or folding kinetics of either protein. Thus, the FRET-monitored changes in the two labeled mutant proteins, Cys25-TNB and Cys62-TNB, could be compared directly. Refolding was commenced from unfolded protein in 8M urea, and both the Trp53 to C25-TNB distance and the Trp53 to C62-TNB distance were found to contract upon dilution of urea. The extent of contraction of each distance, which was measured at a few milliseconds of refolding, was dependent continuously on the concentration of urea present during refolding, and was different for the two distances. For either FRET pair, the gradual contraction of distance with a decrease in the concentration of urea in which refolding occurs, was continuous with the contraction of the polypeptide chain that is seen with a decrease in the concentration of urea in the range in which the protein remains completely unfolded. It therefore appears that the products of the initial sub-millisecond refolding reaction of barstar are collapsed forms, whose dimensions do not change cooperatively in an all-or-none manner, but instead, change gradually with a change in concentration of urea. Thus, the sub-millisecond polypeptide chain collapse reaction of barstar upon denaturant dilution, appears to be a continuous structural transition.

Topology characterization of a benzodiazepine-binding beta-rich domain of the GABAA receptor alpha1 subunit

Structural investigation of GABAA receptors has been limited by difficulties imposed by its trans-membrane-complex nature. In the present study, the topology of a membrane-proximal beta-rich (MPB) domain in the C139-L269 segment of the receptor alpha1 subunit was probed by mapping the benzodiazepine (BZ)-binding and epitopic sites, as well as fluorescence resonance energy transfer (FRET) analysis. Ala-scanning and semiconservative substitutions within this segment revealed the contribution of the phenyl rings of Y160 and Y210, the hydroxy group of S186 and the positive charge on R187 to BZ-binding. FRET with the bound BZ ligand indicated the proximity of Y160, S186, R187, and S206 to the BZ-binding site. On the other hand, epitope-mapping using the monoclonal antibodies (mAbs) against the MPB domain established a clustering of T172, R173, E174, Q196, and T197. Based on the lack of FRET between Trp substitutionally placed at R173 or V198 and bound BZ, this epitope-mapped cluster is located on a separate end of the folded protein from the BZ-binding site. Mutations of the five conserved Cys and Trp residues in the MPB domain gave rise to synergistic and rescuing effects on protein secondary structures and unfolding stability that point to a CCWCW-pentad, reminiscent to the CWC-triad "pin" of immunoglobulin (Ig)-like domains, important for the structural maintenance. These findings, together with secondary structure and fold predictions suggest an anti-parallel beta-strand topology with resemblance to Ig-like fold, having the BZ-binding and the epitopic residues being clustered at two different ends of the fold.

Proteins QSAR with Markov average electrostatic potentials

Classic physicochemical and topological indices have been largely used in small molecules QSAR but less in proteins QSAR. In this study, a Markov model is used to calculate, for the first time, average electrostatic potentials xik for an indirect interaction between aminoacids placed at topologic distances k within a given protein backbone. The short-term average stochastic potential xi1 for 53 Arc repressor mutants was used to model the effect of Alanine scanning on thermal stability. The Arc repressor is a model protein of relevance for biochemical studies on bioorganics and medicinal chemistry. A linear discriminant analysis model developed correctly classified 43 out of 53, 81.1% of proteins according to their thermal stability. More specifically, the model classified 20/28, 71.4% of proteins with near wild-type stability and 23/25, 92.0% of proteins with reduced stability. Moreover, predictability in cross-validation procedures was of 81.0%. Expansion of the electrostatic potential in the series xi0, xi1, xi2, and xi3, justified the use of the abrupt truncation approach, being the overall accuracy >70.0% for xi0 but equal for xi1, xi2, and xi3. The xi1 model compared favorably with respect to others based on D-Fire potential, surface area, volume, partition coefficient, and molar refractivity, with less than 77.0% of accuracy [Ramos de Armas, R.; González-Díaz, H.; Molina, R.; Uriarte, E. Protein Struct. Func. Bioinf.2004, 56, 715]. The xi1 model also has more tractable interpretation than others based on Markovian negentropies and stochastic moments. Finally, the model is notably simpler than the two models based on quadratic and linear indices. Both models, reported by Marrero-Ponce et al., use four-to-five time more descriptors. Introduction of average stochastic potentials may be useful for QSAR applications; having xik amenable physical interpretation and being very effective.

The Structure of the Major Transition State for Folding of an FF Domain from Experiment and Simulation

We have analysed the transition state of folding of the four-helix FF domain from HYPA/FBP11 by high-resolution experiment and simulation as part of a continuing effort to understand the principles of folding and the refinement of predictive methods. The major transition state for folding was subjected to a Phi-value analysis utilising 50 mutants. The transition state contained a nucleus for folding centred around the end of helix 1 (H1) and the beginning of helix 2 (H2). Secondary structure in this region was fully formed (PhiF=0.9-1) and tertiary interactions were well developed. Interactions in the distal part of the native structure were weak (PhiF=0-0.2). The hydrophobic core and other parts of the protein displayed intermediate Phi-values, suggesting that interactions coalesce as the end of H1 and beginning of H2 are in the process of being formed. The distribution of Phi-values resembled that of barnase, which folds via an intermediate, rather than that of CI2 which folds by a concerted nucleation-condensation mechanism. The overall picture of the transition state structure identified in molecular dynamics simulations is in qualitative agreement, with the turn connecting H1 and H2 being formed, a loosened core, and H4 partially unfolded and detached from the core. There are some differences in the details and interpretation of specific Phi-values.

Predicting stability of Arc repressor mutants with protein stochastic moments

As more and more protein structures are determined and applied to drug manufacture, there is increasing interest in studying their stability. In this study, the stochastic moments ((SR)pi(k)) of 53 Arc repressor mutants were introduced as molecular descriptors modeling protein stability. The Linear Discriminant Analysis model developed correctly classified 43 out of 53, 81.13% of proteins according to their thermal stability. More specifically, the model classified 20/28 (71.4%) proteins with near wild-type stability and 23/25 (92%) proteins with reduced stability. Moreover, validation of the model was carried out by re-substitution procedures (81.0%). In addition, the stochastic moments based model compared favorably with respect to others based on physicochemical and geometric parameters such as D-Fire potential, surface area, volume, partition coefficient, and molar refractivity, which presented less than 77% of accuracy. This result illustrates the possibilities of the stochastic moments' method for the study of bioorganic and medicinal chemistry relevant proteins.

Inter-residue interactions in protein folding and stability

During the process of protein folding, the amino acid residues along the polypeptide chain interact with each other in a cooperative manner to form the stable native structure. The knowledge about inter-residue interactions in protein structures is very helpful to understand the mechanism of protein folding and stability. In this review, we introduce the classification of inter-residue interactions into short, medium and long range based on a simple geometric approach. The features of these interactions in different structural classes of globular and membrane proteins, and in various folds have been delineated. The development of contact potentials and the application of inter-residue contacts for predicting the structural class and secondary structures of globular proteins, solvent accessibility, fold recognition and ab initio tertiary structure prediction have been evaluated. Further, the relationship between inter-residue contacts and protein-folding rates has been highlighted. Moreover, the importance of inter-residue interactions in protein-folding kinetics and for understanding the stability of proteins has been discussed. In essence, the information gained from the studies on inter-residue interactions provides valuable insights for understanding protein folding and de novo protein design.

Localized nature of the transition-state structure in goat alpha-lactalbumin folding

To investigate whether the structure partially formed in the molten globule folding intermediate of goat alpha-lactalbumin is further organized in the transition state of folding, we constructed a number of mutant proteins and performed Phi-value analysis on them. For this purpose, we measured the equilibrium unfolding transitions and kinetic refolding and unfolding reactions of the mutants using equilibrium and stopped-flow kinetic circular dichroism techniques. The results show that the mutants with mutations located in the A-helix (V8A, L12A), the B-helix (V27A), the beta-domain (L52A, W60A), the C-helix (K93A, L96A), the C-D loop (Y103F), the D-helix (L105A, L110A), and the C-terminal 3(10)-helix (W118F), have low Phi-values, less than 0.2. On the other hand, D87N, which is located on the Ca(2+)-binding site, has a high Phi-value, 0.91, indicating that tight packing of the side-chain around Asp87 occurs in the transition state. One beta-domain mutant (I55V) and three C-helix mutants (I89V, V90A, and I95V) demonstrated intermediate Phi-values, between 0.4 and 0.7. These results indicate that the folding nucleus in the transition state of goat alpha-LA is not extensively distributed over the alpha-domain of the protein, but very localized in a region that contains the Ca(2+)-binding site and the interface between the C-helix and the beta-domain. This is apparently in contrast with the fact that the molten globule state of alpha-lactalbumin has a partially formed structure inside the alpha-domain. It is concluded that the specific docking of the alpha and beta-domains at a domain interface is necessary for this protein to organize its native structure from the molten globule intermediate.

Beyond the EX1 limit: probing the structure of high-energy states in protein unfolding

Hydrogen exchange kinetics in native solvent conditions have been used to explore the conformational fluctuations of an immunoglobulin domain (CD2.domain1). The global folding/unfolding kinetics of the protein are unaltered between pH 4.5 and pH 9.5, allowing us to use the pH-dependence of amide hydrogen/deuterium exchange to characterise conformational states with energies up to 7.2kcal/mol higher than the folded ground state. The study was intended to search for discreet unfolding intermediates in this region of the energy spectrum, their presence being revealed by the concerted exchange behaviour of subsets of amide groups that become accessible at a given free energy, i.e. the spectrum would contain discreet groupings. Protection factors for 58 amide groups were measured across the pH range and the hydrogen-exchange energy profile is described. More interestingly, exchange behaviour could be grouped into three categories; the first two unremarkable, the third unexpected. (1) In 33 cases, amide exchange was dominated by rapid fluctuation, i.e. the free energy difference between the ground state and the rapidly accessed open state is sufficiently low that the contribution from crossing the unfolding barrier is negligible. (2) In 18 cases exchange is dominated by the global folding transition barrier across the whole pH range measured. The relationship between hydroxyl ion concentration and observed exchange rate is hyperbolic, with the limiting rate being that for global unfolding; the so-called EX1 limit. For these, the free energy difference between the folded ground state and any rapidly-accessed open state is too great for the proton to be exchanged through such fluctuations, even at the highest pH employed in this study. (3) For the third group, comprising five cases, we observe a behaviour that has not been described. In this group, as in category 2, the rate of exchange reaches a plateau; the EX1 limit. However, as the intrinsic exchange rate (k(int)) is increased, this limit is breached and the rate begins to rise again. This unintuitive behaviour does not result from pH instability, rather it is a consequence of amide groups experiencing two processes; rapid fluctuation of structure and crossing the global barrier for unfolding. The boundary at which the EX1 limit is overcome is determined by the equilibrium distribution of the fluctuating open and closed states (K(O/C)) and the rate constant for unfolding (k(u)). This critical boundary is reached when k(int)K(O/C)=k(u). Given that, in a simple transition state formalism: k(u)=K(#)k' (where K(#) describes the equilibrium distribution between the transition and ground state and k' describes the rate of a barrierless rearrangement), it follows that if the pH is raised to a level where k(int)=k', then the entire free energy spectrum from ground state to transition state could be sampled.

Quantitative Observation of Backbone Disorder in Native Elastin

Elastin is a key protein in soft tissue function and pathology. Establishing a structural basis for understanding its reversible elasticity has proven to be difficult. Complementary to structure is the important aspect of flexibility and disorder in elastin. We have used solid-state NMR methods to examine polypeptide and hydrate ordering in both elastic (hydrated) and brittle (dry) elastin fibers and conclude (i) that tightly bound waters are absent in both dry and hydrated elastin and (ii) that the backbone in the hydrated protein is highly disordered with large amplitude motions. The hydrate was studied by (2)H and (17)O NMR, and the polypeptide by (13)C and (2)H NMR. Using a two-dimensional (13)C MAS method, an upper limit of S < 0.1 was determined for the backbone carbonyl group order parameter in hydrated elastin. For comparison, S approximately approximately 0.9 in most proteins. The former result is substantiated by two additional observations: the absence of the characteristic (2)H spectrum for stationary amides and "solution-like" (13)C magic angle spinning spectra at 75 degrees C, at which the material retains elasticity. Comparison of the observed shifts with accepted values for alpha-helices, beta-sheets, or random coils indicates a random coil structure at all carbons. These conclusions are discussed in the context of known thermodynamic properties of elastin and, more generally, protein folding. Because coacervation is an entropy-driven process, it is enhanced by the observed backbone disorder, which, we suggest, is the result of high proline content. This view is supported by recent studies of recombinant elastin polypeptides with systematic proline substitutions.

Mechanical unfolding of a titin Ig domain: structure of transition state revealed by combining atomic force microscopy, protein engineering and molecular dynamics simulations

Titin I27 shows a high resistance to unfolding when subject to external force. To investigate the molecular basis of this mechanical stability, protein engineering Phi-value analysis has been combined with atomic force microscopy to investigate the structure of the barrier to forced unfolding. The results indicate that the transition state for forced unfolding is significantly structured, since highly destabilising mutations in the core do not affect the force required to unfold the protein. As has been shown before, mechanical strength lies in the region of the A' and G-strands but, contrary to previous suggestions, the results indicate clearly that side-chain interactions play a significant role in maintaining mechanical stability. Since Phi-values calculated from molecular dynamics simulations are the same as those determined experimentally, we can, with confidence, use the molecular dynamics simulations to analyse the structure of the transition state in detail, and are able to show loss of interactions between the A' and G-strands with associated A-B and E-F loops in the transition state. The key event is not a simple case of loss of hydrogen bonding interactions between the A' and G-strands alone. Comparison with Phi-values from traditional folding studies shows differences between the force and "no-force" transition states but, nevertheless, the region important for kinetic stability is the same in both cases. This explains the correspondence between hierarchy of kinetic stability (measured in stopped-flow denaturant studies) and mechanical strength in these titin domains.

Stability scale and atomic solvation parameters extracted from 1023 mutation experiments

The stability scale of 20 amino acid residues is derived from a database of 1023 mutation experiments on 35 proteins. The resulting scale of hydrophobic residues has an excellent correlation with the octanol-to-water transfer free energy corrected with an additional Flory-Huggins molar-volume term (correlation coefficient r = 0.95, slope = 1.05, and a near zero intercept). Thus, hydrophobic contribution to folding stability is characterized remarkably well by transfer experiments. However, no corresponding correlation is found for hydrophilic residues. Both the hydrophilic portion and the entire scale, however, correlate strongly with average burial accessible surface (r = 0.76 and 0.97, respectively). Such a strong correlation leads to a near uniform value of the atomic solvation parameters for atoms C, S, O/N, O(-0.5), and N(+0.5,1). All are in the range of 12-28 cal x mol(-1) A(-2), close to the original estimate of hydrophobic contribution of 25-30 cal x mol(-1) A(-2) to folding stability. Without any adjustable parameters, the new stability scale and new atomic solvation parameters yielded an accurate prediction of protein-protein binding free energy for a separate database of 21 protein-protein complexes (r = 0.80 and slope = 1.06, and r = 0.83 and slope = 0.93, respectively).

Fast and slow intermediate accumulation and the initial barrier mechanism in protein folding

Do stable intermediates form very early in the protein folding process? New results and a quantity of literature that bear on this issue are examined here. Results available provide little support for early intermediate accumulation before an initial search-dependent nucleation barrier.

Key interactions in neocarzinostatin, a protein of the immunoglobulin fold family

Neocarzinostatin (NCS) is a seven-stranded beta-sandwich protein, the folding of which is similar to that of the variable domains of immunoglobulins (Ig). The investigation of the backbone dynamics of apo-NCS [Izadi-Pruneyre et al. (2001) Protein Sci., 10, 2228-2240] enabled us to identify the involvement of long side-chain residues in maintaining the rigidity of this beta-protein. In the perspective of using this protein for drug targeting, this raises the following question: do these residues also play a key role in the stabilization of the beta-sheet? To investigate this problem, various genetically engineered variants were constructed by mutating these residues to amino acids with shorter aliphatic side chains. These substitutions have no effects on the global fold. However, an important destabilization of the protein, higher than that expected for a simple 'large-to-small' substitution of buried hydrophobic residues, is observed for three mutants, V34A, V21A and V95A. Interestingly, the nature of the residues in these positions is highly conserved in the other Ig-like proteins. The absence of an evolutionary relationship between NCS and the other Ig-like proteins strongly suggests that this hydrophobic core is characteristic of the Ig-fold itself.

Thermodynamics and folding kinetics analysis of the SH3 domain form discrete molecular dynamics

We perform a detailed analysis of the thermodynamics and folding kinetics of the SH3 domain fold with discrete molecular dynamic simulations. We propose a protein model that reproduces some of the experimentally observed thermodynamic and folding kinetic properties of proteins. Specifically, we use our model to study the transition state ensemble of the SH3 fold family of proteins, a set of unstable conformations that fold to the protein native state with probability 1/2. We analyze the participation of each secondary structure element formed at the transition state ensemble. We also identify the folding nucleus of the SH3 fold and test extensively its importance for folding kinetics. We predict that a set of amino acid contacts between the RT-loop and the distal hairpin are the critical folding nucleus of the SH3 fold and propose a hypothesis that explains this result.

Stiffness of the distal loop restricts the structural heterogeneity of the transition state ensemble in SH3 domains

Protein engineering experiments and Phi(F)-value analysis of SH3 domains reveal that their transition state ensemble (TSE) is conformationally restricted, i.e. the fluctuations in the transition state (TS) structures are small. In the TS of src SH3 and alpha-spectrin SH3 the distal loop and the associated hairpin are fully structured, while the rest of the protein is relatively disordered. If native structure predominantly determines the folding mechanism, the findings for SH3 folds raise the question: What are the features of the native topology that determine the nature of the TSE? We propose that the presence of stiff loops in the native state that connect local structural elements (such as the distal hairpin in SH3 domains) conformationally restricts TSE. We validate this hypothesis using the simulations of a "control" system (16 residue beta-hairpin forming C-terminal fragment of the GBl protein) and its variants. In these fragments the role of bending rigidity in determining the nature of the TSE can be directly examined without complications arising from interactions with the rest of the protein. The TSE structures in the beta-hairpins are determined computationally using cluster analysis and limited Phi(F)-value analysis. Both techniques prove that the conformational heterogeneity decreases as the bending rigidity of the loop increases. To extend this finding to SH3 domains a measure of bending rigidity based on loop curvature, which utilizes native structures in the Protein Data Bank (PDB), is introduced. Using this measure we show that, with few exceptions, the ordering of stiffness of the distal, n-src, and RT loops in the 29 PDB structures of SH3 domains is conserved. Combining the simulation results for beta-hairpins and the analysis of PDB structures for SH3 domains, we propose that the stiff distal loop restricts the conformational fluctuations in the TSE. We also predict that constraining the distal loop to be preformed in the denatured ensemble should not alter the nature of TSE. On the other hand, if the amino and carboxy terminals are cross-linked to form a circular polypeptide chain, the pathways and TSs are altered. These contrasting scenarios are illustrated using simulations of cross-linked WT beta-hairpin fragments. Computations of bending rigidities for immunoglobulin-like domain proteins reveal no clear separation in the stiffness of their loops. In the beta-sandwich proteins, which have large fractions of non-local native contacts, the nature of the TSE cannot be apparently determined using purely local structural characteristics. Nevertheless, the measure of loop stiffness still provides qualitative predictions of the ordered regions in the TSE of Ig27 and TenFn3.

Residues participating in the protein folding nucleus do not exhibit preferential evolutionary conservation

To what extent does natural selection act to optimize the details of protein folding kinetics? In an effort to address this question, the relationship between an amino acid's evolutionary conservation and its role in protein folding kinetics has been investigated intensively. Despite this effort, no consensus has been reached regarding the degree to which residues involved in native-like transition state structure (the folding nucleus) are conserved. Here we report the results of an exhaustive, systematic study of sequence conservation among residues known to participate in the experimentally (Phi-value) defined folding nuclei of all of the appropriately characterized proteins reported to date. We observe no significant evidence that these residues exhibit any anomalous sequence conservation. We do observe, however, a significant bias in the existing kinetic data: the mean sequence conservation of the residues that have been the subject of kinetic characterization is greater than the mean sequence conservation of all residues in 13 of 14 proteins studied. This systematic experimental bias gives rise to the previous observation that the median conservation of residues reported to participate in the folding nucleus is greater than the median conservation of all of the residues in a protein. When this bias is corrected (by comparing, for example, the conservation of residues known to participate in the folding nucleus with that of other, kinetically characterized residues) the previously reported preferential conservation is effectively eliminated. In contrast to well-established theoretical expectations, both poorly and highly conserved residues are apparently equally likely to participate in the protein-folding nucleus.

Role of backbone solvation in determining thermodynamic beta propensities of the amino acids

There is a paradox concerning the beta propensities of the amino acids: the amino acids with the highest beta propensities such as valine and isoleucine have the highest tendency to desolvate the peptide backbone, which should result in a loss of stability. Nevertheless, backbone solvation, calculated as electrostatic solvation free energy (ESF), is highly correlated with mutant stability in the zinc-finger system studied by Kim and Berg [Kim, C. A. & Berg, J. M. (1993) Nature (London) 362, 267-270], and valine and isoleucine are among the most stabilizing amino acids. This inverse correlation between stability and ESF can be explained, because the mutant ESF differences in the unfolded protein are larger than in the native protein. Consequently, mutations such as Ala to Val destabilize the unfolded form more than the native protein. By comparing mutant Delta ESF values in isolated beta-strands versus beta-sheets, we conclude that amino acids with high beta propensities should exert their stabilizing effects at early stages in folding. This deduction agrees with the studies by Clarke and coworkers [Lorch, M., Mason, J. M., Clarke, A. R. & Parker, M. J. (1999) Biochemistry 38, 1377-1385, and Lorch, M., Mason, J. M., Sessions, R. B. & Clarke, A. R. (2000) Biochemistry 39, 3480-3485] of the thermodynamics of folding of the beta-sheet protein CD2.d1.

Understanding hierarchical protein evolution from first principles

We propose a model that explains the hierarchical organization of proteins in fold families. The model, which is based on the evolutionary selection of proteins by their native state stability, reproduces patterns of amino acids conserved across protein families. Due to its dynamic nature, the model sheds light on the evolutionary time-scales. By studying the relaxation of the correlation function between consecutive mutations at a given position in proteins, we observe separation of the evolutionary time-scales: at short time intervals families of proteins with similar sequences and structures are formed, while at long time intervals the families of structurally similar proteins that have low sequence similarity are formed. We discuss the evolutionary implications of our model. We provide a "profile" solution to our model and find agreement between predicted patterns of conserved amino acids and those actually observed in nature.

Protein folding theory: from lattice to all-atom models

This review focuses on recent advances in understanding protein folding kinetics in the context of nucleation theory. We present basic concepts such as nucleation, folding nucleus, and transition state ensemble and then discuss recent advances and challenges in theoretical understanding of several key aspects of protein folding kinetics. We cover recent topology-based approaches as well as evolutionary studies and molecular dynamics approaches to determine protein folding nucleus and analyze other aspects of folding kinetics. Finally, we briefly discuss successful all-atom Monte-Carlo simulations of protein folding and conclude with a brief outlook for the future.

Mapping the Folding Pathway of an Immunoglobulin Domain

BACKGROUND

Do proteins that have the same structure fold by the same pathway even when they are unrelated in sequence? To address this question, we are comparing the folding of a number of different immunoglobulin-like proteins. Here, we present a detailed protein engineering phi value analysis of the folding pathway of TI I27, an immunoglobulin domain from human cardiac titin.

RESULTS

TI I27 folds rapidly via a kinetic intermediate that is destabilized by most mutations. The transition state for folding is remarkably native-like in terms of solvent accessibility. We use phi value analysis to map this transition state and show that it is highly structured; only a few residues close to the N-terminal region of the protein remain completely unfolded. Interestingly, most mutations cause the transition state to become less native-like. This anti-Hammond behavior can be used as a novel means of obtaining additional structural information about the transition state.

CONCLUSIONS

The residues that are involved in nucleating the folding of TI I27 are structurally equivalent to the residues that form the folding nucleus in an evolutionary unrelated fibronectin type III protein. These residues form part of the common structural core of Ig-like domains. The data support the hypothesis that interactions essential for defining the structure of these beta sandwich proteins are also important in nucleation of folding.

Ab initio protein structure prediction using physicochemical potentials and a simplified off-lattice model

This study describes a computational method for ab inito protein structure prediction. Protein conformation has been modeled by using six optimized backbone torsion angles and fixed side chains approximating rotationally averaged real side chains. The approximations aim to keep complexity of the structure description to a minimum without seriously compromising the accuracy of the structural representation. An evolutionary Monte Carlo algorithm has been developed to search through this restricted conformational space to locate low-energy protein structures. A simple physicochemical force field has been developed to assess the energies of different conformations within this structural description. The corresponding residue interaction energies are based on hydrophobic, hydrophilic, steric, and hydrogen-bonding potentials. The search procedure has been used to locate native energy minima from primary sequence alone. The 3-D structures of polypeptides up to 38 residues with both beta and alpha secondary structural elements have been accurately predicted. The search procedure has been found to be highly efficient and follows an energetically and structurally plausible pathway to locate native populations. The simple force field described in the study has been compared with a more complex all-atom model and been found to be similarly effective in predicting the structures of proposed independent folding units. Proteins 2001;43:186-202.

Evolutionary conservation of the folding nucleus

Here, we present statistical analysis of conservation profiles in families of homologous sequences for nine proteins whose folding nucleus was determined by protein engineering methods. We show that in all but one protein (AcP) folding nucleus residues are significantly more conserved than the rest of the protein. Two aspects of our study are especially important: (i) grouping of amino acid residues into classes according to their physical-chemical properties and (ii) proper normalization of amino acid probabilities that reflects the fact that evolutionary pressure to conserve some amino acid types may itself affect concentration of various amino acid types in protein families. Neglect of any of those two factors may make physical and biological "signals" from conservation profiles disappear.

Using chimeric immunity proteins to explore the energy landscape for alpha-helical protein folding

To address the role of sequence in the folding of homologous proteins, the folding and unfolding kinetics of the all-helical bacterial immunity proteins Im2 and Im9 were characterised, together with six chimeric derivatives of these proteins. We show that both Im2 and Im9 fold rapidly (k(UN)(H(2)O)) approximately 2000 s(-1) at pH 7.0, 25 degrees C) in apparent two-state transitions, through rate-limiting transition states that are highly compact (beta(TS)0.93 and 0.96, respectively). Whilst the folding and unfolding properties of three of the chimeras (Im2 (1-44)(Im9), Im2 (1-64)(Im9 )and Im2 (25-44)(Im9)) are similar to their parental counterparts, in other chimeric proteins the introduced sequence variation results in altered kinetic behaviour. At low urea concentrations, Im2 (1-29)(Im9) and Im2 (56-64)(Im9) fold in two-state transitions via transition states that are significantly less compact (beta(TS) approximately 0.7) than those characterised for the other immunity proteins presented here. At higher urea concentrations, however, the rate-limiting transition state for these two chimeras switches or moves to a more compact species (beta(TS) approximately 0.9). Surprisingly, Im2 (30-64)(Im9) populates a highly collapsed species (beta(I)=0.87) in the dead-time (2.5 ms) of stopped flow measurements. These data indicate that whilst topology may place significant constraints on the folding process, specific inter-residue interactions, revealed here through multiple sequence changes, can modulate the ruggedness of the folding energy landscape.

Forced unfolding modulated by disulfide bonds in the Ig domains of a cell adhesion molecule

Cell adhesion molecules (CAMs) mediate cell attachment and stress transfer through extracellular domains. Here we forcibly unfold the Ig domains of a prototypical Ig superfamily CAM that contains intradomain disulfide bonds. The Ig domains of all such CAMs have conformations homologous to cadherin extracellular domains, titin Ig-type domains, and fibronectin type-III (FNIII) domains. Atomic force microscopy has been used to extend the five Ig domains of Mel-CAM (melanoma CAM)--a protein that is overexpressed in metastatic melanomas--under conditions where the disulfide bonds were either left intact or disrupted through reduction. Under physiological conditions where intradomain disulfide bonds are intact, partial unfolding was observed at forces far smaller than those reported previously for either titin's Ig-type domains or tenascin's FNIII domains. This partial unfolding under low force may be an important mechanism for imparting elasticity to cell-cell contacts, as well as a regulatory mechanism for adhesive interactions. Under reducing conditions, Mel-CAM's Ig domains were found to fully unfold through a partially folded state and at slightly higher forces. The results suggest that, in divergent evolution of all such domains, stabilization imparted by disulfide bonds relaxes requirements for strong, noncovalent, folded-state interactions.