Coupled folding-binding in a hydrophobic/polar protein model: impact of synergistic folding and disordered flanks

Combined prediction and design reveals the target recognition mechanism of an intrinsically disordered protein interaction domain

An increasing number of protein interaction domains and their targets are being found to be intrinsically disordered proteins (IDPs). The corresponding target recognition mechanisms are mostly elusive because of challenges in performing detailed structural analysis of highly dynamic IDP-IDP complexes. Here, we show that by combining recently developed computational approaches with experiments, the structure of the complex between the intrinsically disordered C-terminal domain (CTD) of protein 4.1G and its target IDP region in NuMA can be dissected at high resolution. First, we carry out systematic mutational scanning using dihydrofolate reductase-based protein complementarity analysis to identify essential interaction regions and key residues. The results are found to be highly consistent with an α/β-type complex structure predicted by AlphaFold2 (AF2). We then design mutants based on the predicted structure using a deep learning protein sequence design method. The solved crystal structure of one mutant presents the same core structure as predicted by AF2. Further computational prediction and experimental assessment indicate that the well-defined core structure is conserved across complexes of 4.1G CTD with other potential targets. Thus, we reveal that an intrinsically disordered protein interaction domain uses an α/β-type structure module formed through synergistic folding to recognize broad IDP targets. Moreover, we show that computational prediction and experiment can be jointly applied to segregate true IDP regions from the core structural domains of IDP-IDP complexes and to uncover the structure-dependent mechanisms of some otherwise elusive IDP-IDP interactions.

Understanding protein diffusion on force-induced stretched DNA conformation

DNA morphology is subjected to environmental conditions and is closely coupled with its function. For example, DNA experiences stretching forces during several biological processes, including transcription and genome transactions, that significantly alter its conformation from that of B-DNA. Indeed, a well-defined 1.5 times extended conformation of dsDNA, known as Σ-DNA, has been reported in DNA complexes with proteins such as Rad51 and RecA. A striking feature in Σ-DNA is that the nucleobases are partitioned into triplets of three locally stacked bases separated by an empty rise gap of ∼ 5 Å. The functional role of such a DNA base triplet was hypothesized to be coupled with the ease of recognition of DNA bases by DNA-binding proteins (DBPs) and the physical origin of three letters (codon/anti-codon) in the genetic code. However, the underlying mechanism of base-triplet formation and the ease of DNA base-pair recognition by DBPs remain elusive. To investigate, here, we study the diffusion of a protein on a force-induced stretched DNA using coarse-grained molecular dynamics simulations. Upon pulling at the 3' end of DNA by constant forces, DNA exhibits a conformational transition from B-DNA to a ladder-like S-DNA conformation via Σ-DNA intermediate. The resulting stretched DNA conformations exhibit non-uniform base-pair clusters such as doublets, triplets, and quadruplets, of which triplets are energetically more stable than others. We find that protein favors the triplet formation compared to its unbound form while interacting non-specifically along DNA, and the relative population of it governs the ruggedness of the protein-DNA binding energy landscape and enhances the efficiency of DNA base recognition. Furthermore, we analyze the translocation mechanism of a DBP under different force regimes and underscore the significance of triplet formation in regulating the facilitated diffusion of protein on DNA. Our study, thus, provides a plausible framework for understanding the structure-function relationship between triplet formation and base recognition by a DBP and helps to understand gene regulation in complex regulatory processes.

Crowding-induced protein destabilization in the absence of soft attractions

It is generally assumed that volume exclusion by macromolecular crowders universally stabilizes the native states of proteins and destabilization suggests soft attractions between crowders and protein. Here we show that proteins can be destabilized even by crowders that are purely repulsive. With a coarse-grained sequence-based model, we study the folding thermodynamics of two sequences with different native folds, a helical hairpin and a β-barrel, in a range of crowder volume fractions, φ_c. We find that the native state, N, remains structurally unchanged under crowded conditions, while the size of the unfolded state, U, decreases monotonically with φ_c. Hence, for all φ_c>0, U is entropically disfavored relative to N. This entropy-centric view holds for the helical hairpin protein, which is stabilized under all crowded conditions as quantified by changes in either the folding midpoint temperature, T_m, or the free energy of folding. We find, however, that the β-barrel protein is destabilized under low-T, low-φ_c conditions. This destabilization can be understood from two characteristics of its folding: 1) a relatively compact U at T<T_m, such that U is only weakly disfavored entropically by the crowders; and 2) a transient, compact, and relatively low-energy nonnative state that has a maximum population of only a few percent at φ_c=0, but increasing monotonically with φ_c. Overall, protein destabilization driven by hard-core effects appears possible when a compaction of U leads to even a modest population of compact nonnative states that are energetically competitive with N.

Molecular dynamics simulations and biochemical characterization of Pf14-3-3 and PfCDPK1 interaction towards its role in growth of human malaria parasite

Scaffold proteins play pivotal role as modulators of cellular processes by operating as multipurpose conformation clamps. 14-3-3 proteins are gold-standard scaffold modules that recognize phosphoSer/Thr (pS/pT) containing conserved motifs, and confer conformational changes leading to modulation of functional parameters of their target proteins. Modulation in functional activity of kinases has been attributed to their interaction with 14-3-3 proteins. Herein, we have annotated and characterized PF3D7_0818200 as 14-3-3 isoform I in Plasmodium falciparum 3D7, and its interaction with one of the key kinases of the parasite, Calcium-Dependent Protein Kinase 1 (CDPK1) by performing various analytical biochemistry and biophysical assays. Molecular dynamics simulation studies indicated that CDPK1 polypeptide sequence (61KLGpS64) behaves as canonical Mode I-type (RXXpS/pT) consensus 14-3-3 binding motif, mediating the interaction. The 14-3-3I/CDPK1 interaction was validated in vitro with ELISA and SPR, which confirmed that the interaction is phosphorylation dependent, with binding affinity constant of 670 ± 3.6 nM. The interaction of 14-3-3I with CDPK1 was validated with well characterized optimal 14-3-3 recognition motifs: Mode I-type ARSHpSYPA and Mode II-type RLYHpSLPA, by simulation studies and ITC. This interaction was found to marginally enhance CDPK1 functional activity. Furthermore, interaction antagonizing peptidomimetics showed growth inhibitory impact on the parasite indicating crucial physiological role of 14-3-3/CDPK1 interaction. Overall, this study characterizes 14-3-3I as a scaffold protein in the malaria parasite and unveils CDPK1 as its previously unidentified target. This sets a precedent for the rational design of 14-3-3 based PPI inhibitors by utilizing 14-3-3 recognition motif peptides, as a potential antimalarial strategy.

Mechanism of Dynamic Binding of Replication Protein A to ssDNA

Replication protein A (RPA) serves as a hub protein inside eukaryotic cells, where it coordinates crucial DNA metabolic processes and activates the DNA-damage response system. A characteristic feature of its action is to associate with single-stranded DNA (ssDNA) intermediates before handing them over to downstream proteins. The length of ssDNA intermediates differs for different pathways. This means that RPA must have mechanisms for selective processing of ssDNA intermediates based on their length, the knowledge of which is fundamental to elucidate when and how DNA repair and replication processes are symphonized. By employing extensive molecular dynamics simulations, we investigated the mechanism of binding of RPA to ssDNA of different lengths. We show that the binding involves dynamic equilibrium with a stable intermediate, the population of which increases with the length of ssDNA. The vital underlying factors are decoded through collective variable principal component analysis. It suggests a differently orchestrated set of interactions that define the action of RPA based on the length of ssDNA intermediates. We further estimated the association kinetics that matches excellently well with previous experimental studies and probed the diffusion mechanism of RPA to ssDNA. RPA diffuses on short ssDNA through progressive "bulge" formation. With long ssDNA, we observed a conformational change in ssDNA coupled with its binding to RPA in a cooperative fashion. This unanticipated binding mechanism successfully explains how the "short-lived", long ssDNA intermediates are processed quickly in vivo. This study thus reveals the molecular basis of several recent experimental observations related to RPA binding to ssDNA and provides novel insights into the RPA functioning in DNA repair and replication.

Effects of Topology and Sequence in Protein Folding Linked via Conformational Fluctuations

Experiments have compared the folding of proteins with different amino acid sequences but the same basic structure, or fold. Results indicate that folding is robust to sequence variations for proteins with some nonlocal folds, such as all-β, whereas the folding of more local, all-α proteins typically exhibits a stronger sequence dependence. Here, we use a coarse-grained model to systematically study how variations in sequence perturb the folding energy landscapes of three model sequences with 3α, 4β + α, and β-barrel folds, respectively. These three proteins exhibit folding features in line with experiments, including expected rank order in the cooperativity of the folding transition and stability-dependent shifts in the location of the free-energy barrier to folding. Using a generalized-ensemble simulation approach, we determine the thermodynamics of around 2000 sequence variants representing all possible hydrophobic or polar single- and double-point mutations. From an analysis of the subset of stability-neutral mutations, we find that folding is perturbed in a topology-dependent manner, with the β-barrel protein being the most robust. Our analysis shows, in particular, that the magnitude of mutational perturbations of the transition state is controlled in part by the size or "width" of the underlying conformational ensemble. This result suggests that the mutational robustness of the folding of the β-barrel protein is underpinned by its conformationally restricted transition state ensemble, revealing a link between sequence and topological effects in protein folding.

Disparity in anomalous diffusion of proteins searching for their target DNA sites in a crowded medium is controlled by the size, shape and mobility of macromolecular crowders

Using extensive computer simulations, we analyzed the role of physical properties of molecular crowding agents such as size, shape and mobility in the target search dynamics of DNA binding proteins. Our main result is that the sub-diffusive dynamics of a protein inside a crowded medium strongly depends on the crowder properties and also on the protein's mode of diffusion. For instance, while scanning the DNA one-dimensionally, the protein dynamics does not vary with the change in crowder properties. Conversely, the diffusion exponent varies non-monotonically during 3D diffusion and is maximally affected when the crowders match the protein physically. The investigation shows that the effect stems from the ruggedness of the associated potential energy landscape, which is regulated by the protein-crowder and DNA-crowder interactions. Our findings have broad significance in understanding the target search dynamics of proteins on DNA in crowded cellular milieu and selecting appropriate crowding agents when designing in vitro experiments.

Modulation of the aggregation of an amyloidogenic sequence by flanking-disordered region in the intrinsically disordered antigen merozoite surface protein 2

The abundant Plasmodium falciparum merozoite surface protein MSP2, a potential malaria vaccine candidate, is an intrinsically disordered protein with some nascent secondary structure present in its conserved N-terminal region. This relatively ordered region has been implicated in both membrane interactions and amyloid-like aggregation of the protein, while the significance of the flanking-disordered region is unclear. In this study, we show that aggregation of the N-terminal conserved region of MSP2 is influenced in a length- and sequence-dependent fashion by the disordered central variable sequences. Intriguingly, MSP2 peptides containing the conserved region and the first five residues of the variable disordered regions aggregated more rapidly than a peptide corresponding to the conserved region alone. In contrast, MSP2 peptides extending 8 or 12 residues into the disordered region aggregated more slowly, consistent with the expected inhibitory effect of flanking-disordered sequences on the aggregation of amyloidogenic ordered sequences. Computational analyses indicated that the helical propensity of the ordered region of MSP2 was modulated by the adjacent disordered five residues in a sequence-dependent manner. Nuclear magnetic resonance and circular dichroism spectroscopic studies with synthetic peptides confirmed the computational predictions, emphasizing the correlation between aggregation propensity and conformation of the ordered region and the effects thereon of the adjacent disordered region. These results show that the effects of flanking-disordered sequences on a more ordered sequence may include enhancement of aggregation through modulation of the conformational properties of the more ordered sequence.

Intrinsic Disorder in a Well-Folded Globular Protein

The folded structure of the heterodimeric sweet protein monellin mimics single-chain proteins with topology β1-α1-β2-β3-β4-β5 (chain A: β3-β4-β5; chain B: β1-α1-β2). Furthermore, like naturally occurring single-chain proteins of a similar size, monellin folds cooperatively with no detectable intermediates. However, the two monellin chains, A and B, are marginally structured in isolation and fold only upon binding to each other. Thus, monellin presents a unique opportunity to understand the design of intrinsically disordered proteins that fold upon binding. Here, we study the folding of a single-chain variant of monellin (scMn) using simulations of an all heavy-atom structure-based model. These simulations can explain mechanistic details derived from scMn experiments performed using several different structural probes. scMn folds cooperatively in our structure-based simulations, as is also seen in experiments. We find that structure formation near the transition-state ensemble of scMn is not uniformly distributed but is localized to a hairpin-like structure which contains one strand from each chain (β2, β3). Thus, the sequence and the underlying energetics of heterodimeric monellin promote the early formation of the interchain interface (β2-β3). By studying computational scMn mutants whose "interchain" interactions are deleted, we infer that this energy distribution allows the two protein chains to remain largely disordered when this interface is not folded. From these results, we suggest that cutting the protein backbone of a globular protein between residues which lie within its folding nucleus may be one way to construct two disordered fragments which fold upon binding.

Probing the dynamics of disorder

Intrinsically disordered proteins (IDPs) play an important role in many diseases. IDPs are a large and important class of proteins; estimated to represent a significant fraction of many genomes. In contrast to protein-protein interactions between well-folded proteins, IDPs typically bind to targets using short consecutive stretches of amino acids. Structures of IDPs complexed with a target have shown great diversity in binding modes. However, how this binding diversity is achieved at the molecular level is not well understood. Unfortunately, the prediction and detailed characterization of IDPs experimentally is still a very challenging task; however molecular mechanics based molecular dynamics simulation are well suited for studying the dynamic behavior of IDPs. We look into the current state for force fields for simulating IDPs and an example of how these methods have been applied to the p53 protein. p53 is one of the most extensively studied IDPs, with multiple intrinsically disordered regulatory domains that mediate its interactions with many other proteins engaged in multiple biological pathways. We show how molecular dynamics simulations can be used to elucidate on the mechanisms involved in selection of the different binding partners.

Fuzzy regions in an intrinsically disordered protein impair protein-protein interactions

Despite the partial disorder-to-order transition that intrinsically disordered proteins often undergo upon binding to their partners, a considerable amount of residual disorder may be retained in the bound form, resulting in a fuzzy complex. Fuzzy regions flanking molecular recognition elements may enable partner fishing through non-specific, transient contacts, thereby facilitating binding, but may also disfavor binding through various mechanisms. So far, few computational or experimental studies have addressed the effect of fuzzy appendages on partner recognition by intrinsically disordered proteins. In order to shed light onto this issue, we used the interaction between the intrinsically disordered C-terminal domain of the measles virus (MeV) nucleoprotein (NTAIL ) and the X domain (XD) of the viral phosphoprotein as model system. After binding to XD, the N-terminal region of NTAIL remains conspicuously disordered, with α-helical folding taking place only within a short molecular recognition element. To study the effect of the N-terminal fuzzy region on NTAIL /XD binding, we generated N-terminal truncation variants of NTAIL , and assessed their binding abilities towards XD. The results revealed that binding increases with shortening of the N-terminal fuzzy region, with this also being observed with hsp70 (another MeV NTAIL binding partner), and for the homologous NTAIL /XD pairs from the Nipah and Hendra viruses. Finally, similar results were obtained when the MeV NTAIL fuzzy region was replaced with a highly dissimilar artificial disordered sequence, supporting a sequence-independent inhibitory effect of the fuzzy region.

Smooth functional transition along a mutational pathway with an abrupt protein fold switch

Recent protein design experiments have demonstrated that proteins can migrate between folds through the accumulation of substitution mutations without visiting disordered or nonfunctional points in sequence space. To explore the biophysical mechanism underlying such transitions we use a three-letter continuous protein model with seven atoms per amino acid to provide realistic sequence-structure and sequence-function mappings through explicit simulation of the folding and interaction of model sequences. We start from two 16-amino-acid sequences folding into an α-helix and a β-hairpin, respectively, each of which has a preferred binding partner with 35 amino acids. We identify a mutational pathway between the two folds, which features a sharp fold switch. By contrast, we find that the transition in function is smooth. Moreover, the switch in preferred binding partner does not coincide with the fold switch. Discovery of new folds in evolution might therefore be facilitated by following fitness slopes in sequence space underpinned by binding-induced conformational switching.

Knotted proteins: A tangled tale of Structural Biology

Knotted proteins have their native structures arranged in the form of an open knot. In the last ten years researchers have been making significant efforts to reveal their folding mechanism and understand which functional advantage(s) knots convey to their carriers. Molecular simulations have been playing a fundamental role in this endeavor, and early computational predictions about the knotting mechanism have just been confirmed in wet lab experiments. Here we review a collection of simulation results that allow outlining the current status of the field of knotted proteins, and discuss directions for future research.

Subunit Interactions within the Carbon-Phosphorus Lyase Complex from Escherichia coli

Phosphonates are a large class of organophosphorus compounds with a characteristic carbon-phosphorus bond. The genes responsible for phosphonate utilization in Gram-negative bacteria are arranged in an operon of 14 genes. The carbon-phosphorus lyase complex, encoded by the genes phnGHIJKLM, catalyzes the cleavage of the stable carbon-phosphorus bond of organophosphonates to the corresponding hydrocarbon and inorganic phosphate. Recently, complexes of this enzyme containing five subunits (PhnG-H-I-J-K), four subunits (PhnG-H-I-J), and two subunits (PhnG-I) were purified after expression in Escherichia coli ( Proc. Natl. Acad. Sci., U. S. A. 2011 , 108 , 11393 ). Here we demonstrated using mass spectrometry, ultracentrifugation, and chemical cross-linking experiments that these complexes are formed from a PhnG2I2 core that is further elaborated by the addition of two copies each of PhnH and PhnJ to generate PhnG2H2I2J2. This complex adds an additional subunit of PhnK to form PhnG2H2I2J2K. Chemical cross-linking of the five-component complex demonstrated that PhnJ physically interacts with both PhnG and PhnI. We were unable to demonstrate the interaction of PhnH or PhnK with any other subunits by chemical cross-linking. Hydrogen-deuterium exchange was utilized to probe for alterations in the dynamic properties of individual subunits within the various complexes. Significant regions of PhnG become less accessible to hydrogen/deuterium exchange from solvent within the PhnG2I2 complex compared with PhnG alone. Specific regions of PhnI exhibited significant differences in the H/D exchange rates in PhnG2I2 and PhnG2H2I2J2K.

Theoretical perspectives on nonnative interactions and intrinsic disorder in protein folding and binding

The diverse biological functions of intrinsically disordered proteins (IDPs) have markedly raised our appreciation of protein conformational versatility, whereas the existence of energetically favorable yet functional detrimental nonnative interactions underscores the physical limitations of evolutionary optimization. Here we survey recent advances in using biophysical modeling to gain insight into experimentally observed nonnative behaviors and IDP properties. Simulations of IDP interactions to date focus mostly on coupled folding-binding, which follows essentially the same organizing principle as the local-nonlocal coupling mechanism in cooperative folding of monomeric globular proteins. By contrast, more innovative theories of electrostatic and aromatic interactions are needed for the conceptually novel but less-explored 'fuzzy' complexes in which the functionally bound IDPs remain largely disordered.

Hybrid Monte Carlo with non-uniform step size

The Hybrid Monte Carlo method offers a rigorous and potentially efficient approach to the simulation of dense systems, by combining numerical integration of Newton's equations of motion with a Metropolis accept-or-reject step. The Metropolis step corrects for sampling errors caused by the discretization of the equations of motion. The integration is usually performed using a uniform step size. Here, we present simulations of the Lennard-Jones system showing that the use of smaller time steps in the tails of each integration trajectory can reduce errors in energy. The acceptance rate is 10-15 percentage points higher in these runs, compared to simulations with the same trajectory length and the same number of integration steps but a uniform step size. We observe similar effects for the harmonic oscillator and a coarse-grained peptide model, indicating generality of the approach.

Effects of desolvation barriers and sidechains on local–nonlocal coupling and chevron behaviors in coarse-grained models of protein folding

Coarse-grained protein chain models with desolvation barriers or sidechains lead to stronger local–nonlocal coupling and more linear chevron plots.

Polycation-π interactions are a driving force for molecular recognition by an intrinsically disordered oncoprotein family

Molecular recognition by intrinsically disordered proteins (IDPs) commonly involves specific localized contacts and target-induced disorder to order transitions. However, some IDPs remain disordered in the bound state, a phenomenon coined “fuzziness”, often characterized by IDP polyvalency, sequence-insensitivity and a dynamic ensemble of disordered bound-state conformations. Besides the above general features, specific biophysical models for fuzzy interactions are mostly lacking. The transcriptional activation domain of the Ewing's Sarcoma oncoprotein family (EAD) is an IDP that exhibits many features of fuzziness, with multiple EAD aromatic side chains driving molecular recognition. Considering the prevalent role of cation-π interactions at various protein-protein interfaces, we hypothesized that EAD-target binding involves polycation- π contacts between a disordered EAD and basic residues on the target. Herein we evaluated the polycation-π hypothesis via functional and theoretical interrogation of EAD variants. The experimental effects of a range of EAD sequence variations, including aromatic number, aromatic density and charge perturbations, all support the cation-π model. Moreover, the activity trends observed are well captured by a coarse-grained EAD chain model and a corresponding analytical model based on interaction between EAD aromatics and surface cations of a generic globular target. EAD-target binding, in the context of pathological Ewing's Sarcoma oncoproteins, is thus seen to be driven by a balance between EAD conformational entropy and favorable EAD-target cation-π contacts. Such a highly versatile mode of molecular recognition offers a general conceptual framework for promiscuous target recognition by polyvalent IDPs.

Understanding how proteins recognize each other is central to deciphering the inner workings of living things and for biomedical research. It has long been known that the sequence of a protein, which is a string of different amino acids, can dictate how a protein molecule folds into a well-defined shape required for biological tasks. Many folded proteins recognize and bind with each other by a tight geometric fit similar to that between a lock and its key. Recently, however, it has become clear that some proteins function as a flexible string, in constant motion, without forming a stable shape. Understanding how such “disordered” proteins work is challenging. To gain insight, we studied a disordered protein region that causes a large family of human cancers. Employing an innovative combination of experimental and theoretical techniques, we describe a new mode of protein interaction based on multiple simple contacts between one type of amino acid (aromatic) in the disordered protein and another type (positively charged) on the partner protein. Because this mechanism also underlies the ability of the disordered protein to cause cancer, further investigation of this unprecedented mode of protein-protein interaction may open up new avenues for cancer therapy.

Multiscaled exploration of coupled folding and binding of an intrinsically disordered molecular recognition element in measles virus nucleoprotein

Numerous relatively short regions within intrinsically disordered proteins (IDPs) serve as molecular recognition elements (MoREs). They fold into ordered structures upon binding to their partner molecules. Currently, there is still a lack of in-depth understanding of how coupled binding and folding occurs in MoREs. Here, we quantified the unbound ensembles of the α-MoRE within the intrinsically disordered C-terminal domain of the measles virus nucleoprotein. We developed a multiscaled approach by combining a physics-based and an atomic hybrid model to decipher the mechanism by which the α-MoRE interacts with the X domain of the measles virus phosphoprotein. Our multiscaled approach led to remarkable qualitative and quantitative agreements between the theoretical predictions and experimental results (e.g., chemical shifts). We found that the free α-MoRE rapidly interconverts between multiple discrete partially helical conformations and the unfolded state, in accordance with the experimental observations. We quantified the underlying global folding-binding landscape. This leads to a synergistic mechanism in which the recognition event proceeds via (minor) conformational selection, followed by (major) induced folding. We also provided evidence that the α-MoRE is a compact molten globule-like IDP and behaves as a downhill folder in the induced folding process. We further provided a theoretical explanation for the inherent connections between "downhill folding," "molten globule," and "intrinsic disorder" in IDP-related systems. Particularly, we proposed that binding and unbinding of IDPs proceed in a stepwise way through a "kinetic divide-and-conquer" strategy that confers them high specificity without high affinity.

Effects of knots on protein folding properties

This work explores the impact of knots, knot depth and motif of the threading terminus in protein folding properties (kinetics, thermodynamics and mechanism) via extensive Monte Carlo simulations of lattice models. A knotted backbone has no effect on protein thermodynamic stability but it may affect key aspects of folding kinetics. In this regard, we found clear evidence for a functional advantage of knots: knots enhance kinetic stability because a knotted protein unfolds at a distinctively slower rate than its unknotted counterpart. However, an increase in knot deepness does not necessarily lead to more effective changes in folding properties. In this regard, a terminus with a non-trivial conformation (e.g. hairpin) can have a more dramatic effect in enhancing kinetic stability than knot depth. Nevertheless, our results suggest that the probability of the denatured ensemble to keep knotted is higher for proteins with deeper knots, indicating that knot depth plays a role in determining the topology of the denatured state. Refolding simulations starting from denatured knotted conformations show that not every knot is able to nucleate folding and further indicate that the formation of the knotting loop is a key event in the folding of knotted trefoils. They also show that there are specific native contacts within the knotted core that are crucial to keep a native knotting loop in denatured conformations which otherwise have no detectable structure. The study of the knotting mechanism reveals that the threading of the knotting loop generally occurs towards late folding in conformations that exhibit a significant degree of structural consolidation.

Inherent relationships among different biophysical prediction methods for intrinsically disordered proteins

Intrinsically disordered proteins do not have stable secondary and/or tertiary structures but still function. More than 50 prediction methods have been developed and inherent relationships may be expected to exist among them. To investigate this, we conducted molecular simulations and algorithmic analyses on a minimal coarse-grained polypeptide model and discovered a common basis for the charge-hydropathy plot and packing-density algorithms that was verified by correlation analysis. The correlation analysis approach was applied to realistic datasets, which revealed correlations among some physical-chemical properties (charge-hydropathy plot, packing density, pairwise energy). The correlations indicated that these biophysical methods find a projected direction to discriminate ordered and disordered proteins. The optimized projection was determined and the ultimate accuracy limit of the existing algorithms is discussed.

Dissecting partner recognition by an intrinsically disordered protein using descriptive random mutagenesis

In view of getting insights into the molecular determinants of the binding efficiency of intrinsically disordered proteins (IDPs), we used random mutagenesis. As a proof of concept, we chose the interaction between the intrinsically disordered C-terminal domain of the measles virus nucleoprotein (NTAIL) and the X domain (XD) of the viral phosphoprotein and assessed how amino acid substitutions introduced at random within NTAIL affect partner recognition. In contrast with directed evolution approaches, we did not apply any selection and used the gene library approach not for production purposes but for achieving a better understanding of the NTAIL/XD interaction. For that reason, and to differentiate our approach from similar approaches that make use of systematic (i.e., targeted) mutagenesis, we propose to call it "descriptive random mutagenesis" (DRM). NTAIL variants generated by error-prone PCR were picked at random in the absence of selection pressure and were characterized in terms of sequence and binding abilities toward XD. DRM not only identified determinants of NTAIL/XD interaction that were in good agreement with previous work but also provided new insights. In particular, we discovered that the primary interaction site is poorly evolvable in terms of binding abilities toward XD. We also identified a critical NTAIL residue whose role in stabilizing the NTAIL/XD complex had previously escaped detection, and we identified NTAIL regulatory sites that dampen the interaction while being located outside the primary interaction site. Results show that DRM is a valuable approach to study binding abilities of IDPs.

How difficult is it to fold a knotted protein? In silico insights from surface-tethered folding experiments

We explore the effect of surface tethering on the folding process of a lattice protein that contains a trefoil knot in its native structure via Monte Carlo simulations. We show that the outcome of the tethering experiment depends critically on which terminus is used to link the protein to a chemically inert plane. In particular, if surface tethering occurs at the bead that is closer to the knotted core the folding rate becomes exceedingly slow and the protein is not able to find the native structure in all the attempted folding trajectories. Such low folding efficiency is also apparent from the analysis of the probability of knot formation, p(knot), as a function of nativeness. Indeed, p(knot) increases abruptly from ∼0 to ∼1 only when the protein has more than 80% of its native contacts formed, showing that a highly compact conformation must undergo substantial structural re-arrangement in order to get effectively knotted. When the protein is surface tethered by the bead that is placed more far away from the knotted core p(knot) is higher than in the other folding setups (including folding in the bulk), especially if conformations are highly native-like. These results show that the mobility of the terminus closest to the knotted core is critical for successful folding of trefoil proteins, which, in turn, highlights the importance of a knotting mechanism that is based on a threading movement of this terminus through a knotting loop. The results reported here predict that if this movement is blocked, knotting occurs via an alternative mechanism, the so-called spindle mechanism, which is prone to misfolding. Our simulations show that in the three considered folding setups the formation of the knot is typically a late event in the folding process. We discuss the implications of our findings for co-translational folding of knotted trefoils.

Modulation of folding kinetics of repeat proteins: interplay between intra- and interdomain interactions

Repeat proteins have unique elongated structures that, unlike globular proteins, are quite modular. Despite their simple one-dimensional structure, repeat proteins exhibit intricate folding behavior with a complexity similar to that of globular proteins. Therefore, repeat proteins allow one to quantify fundamental aspects of the biophysics of protein folding. One important feature of repeat proteins is the interfaces between the repeating units. In particular, the distribution of stabilities within and between the repeats was previously suggested to affect their folding characteristics. In this study, we explore how the interface affects folding kinetics and cooperativity by investigating two families of repeat proteins, namely, the Ankyrin and tetratricopeptide repeat proteins, which differ in the number of interfacial contacts that are formed between their units as well as in their folding behavior. By using simple topology-based models, we show that modulating the energetic strength of the interface relative to that of the repeat itself can drastically change the protein stability, folding rate, and cooperativity. By further dissecting the interfacial contacts into several subsets, we isolated the effects of each of these groups on folding kinetics. Our study highlights the importance of interface connectivity in determining the folding behavior.