The influence of intrinsic folding mechanism of an unfolded protein on the coupled folding-binding process during target recognition

Intrinsically disordered proteins (IDPs) are extensively involved in dynamic signaling processes which require a high association rate and a high dissociation rate for rapid binding/unbinding events and at the same time a sufficient high affinity for specific recognition. Although the coupled folding-binding processes of IDPs have been extensively studied, it is still impossible to predict whether an unfolded protein is suitable for molecular signaling via coupled folding-binding. In this work, we studied the interplay between intrinsic folding mechanisms and coupled folding-binding process for unfolded proteins through molecular dynamics simulations. We first studied the folding process of three representative IDPs with different folded structures, that is, c-Myb, AF9, and E3 rRNase. We found the folding free energy landscapes of IDPs are downhill or show low barriers. To further study the influence of intrinsic folding mechanism on the binding process, we modulated the folding mechanism of barnase via circular permutation and simulated the coupled folding-binding process between unfolded barnase permutant and folded barstar. Although folding of barnase was coupled to target binding, the binding kinetics was significantly affected by the intrinsic folding free energy barrier, where reducing the folding free energy barrier enhances binding rate up to two orders of magnitude. This accelerating effect is different from previous results which reflect the effect of structure flexibility on binding kinetics. Our results suggest that coupling the folding of an unfolded protein with no/low folding free energy barrier with its target binding may provide a way to achieve high specificity and rapid binding/unbinding kinetics simultaneously.     

Conformational changes associated with L16P and T118M mutations in the membrane-embedded PMP22 protein,consequential in CMT-1A

Peripheral myelin protein 22 (PMP22) resides in the plasma membrane and is required for myelin formation in the peripheral nervous system. Excess PMP22 mutants accumulate in the endoplasmic reticulum (ER) resulting in the inherited neuropathies of Charcot–Marie–Tooth disease. However, there was no evidence of the structure of PMP22 or how mutations affect its folding. Therefore, in this study, we combined bioinformatics and homology modeling approaches to obtain three-dimensional native and mutated PMP22 models and its anchoring to a POPC membrane, submitted to .5-μs MD simulations, to determine how the L16P and T118M mutations affect the conformational behavior of PMP22. In addition, we investigated the ability of the native and mutated species to accumulate in the ER, via interaction with RER1, by combining protein–protein docking and MD simulations, taking the conformations that were most representative of the native and mutated PMP22 systems and RER1 conformations. Principal component analysis over MD simulations revealed that L16P and T118M mutations resulted in increased structural instability compared to the native form, which is consistent with previous experimental findings of increased structural fluctuations along a loop connecting transmembrane α-helix1 and α-helix2. Docking and MD simulations coupled with the MMGBSA approach allowed the identification that the binding interface for the PMP22-RER1 complex takes place through transmembrane α-helix1 and α-helix2, with higher effective binding free energy values between the mutated PMP22 systems and RER1 than for the native PMP22, mainly through van der Waals interactions.     

Ab initio simulation of a 57-residue protein in explicit solvent reproduces the native conformation in the lowest free-energy cluster

An enhanced conformational sampling method, multicanonical molecular dynamics (McMD), was applied to the ab intio folding of the 57-residue first repeat of human glutamyl- prolyl-tRNA synthetase (EPRS-R1) in explicit solvent. The simulation started from a fully extended structure of EPRS-R1 and did not utilize prior structural knowledge. A canonical ensemble, which is a conformational ensemble thermodynamically probable at an arbitrary temperature, was constructed by reweighting the sampled structures. Conformational clusters were obtained from the canonical ensemble at 300 K, and the largest cluster (i.e., the lowest free-energy cluster), which contained 34% of the structures in the ensemble, was characterized by the highest similarity to the NMR structure relative to all alternative clusters. This lowest free-energy cluster included native-like structures composed of two anti-parallel α-helices. The canonical ensemble at 300 K also showed that a short Gly-containing segment, which adopts an α-helix in the native structure, has a tendency to be structurally disordered. Atomic-level analyses demonstrated clearly that inter-residue hydrophobic interactions drive the helix formation of the Gly-containing segment, and that increasing the hydrophobic contacts accompanies exclusion of water molecules from the vicinity of this segment. This study has shown, for the first time, that the free-energy landscape of a structurally well-ordered protein of about 60 residues is obtainable with an all atom model in explicit water without prior structural knowledge.     

The surface of β-sheet proteins contains amphiphilic regions which may provide clues about protein folding

A major bottleneck in the field of biochemistry is our limited understanding of the processes by which a protein folds into its native conformation. Much of the work on this issue has focused on the conserved core of the folded protein. However, one might imagine that a ubiquitous motif for unaided folding or for the recognition of chaperones may involve regions on the surface of the native structure. We explore this possibility by an analysis of the spatial distribution of regions with amphiphilic α-helical potential on the surface of β-sheet proteins. All proteins, Including β-sheet proteins, contain regions with amphiphilic α-helical potential. That is, any α-helix formed by that region would be amphiphilic, having both hydrophobic and hydrophilic surfaces. In the three-dimensional structure of all β-sheet proteins analyzed, we have found a distinct pattern in the spatial distribution of sequences with amphiphilic α-helical potential. The amphiphilic regions occur in ring shaped clusters approximately 20 to 30 Å in diameter on the surface of the protein. In addition, these regions have a strong preference for positively charged amino acids and a lower preference for residues not favorable to α-helix formation. Although the purpose of these amphiphilic regions which are not associated with naturally occurring α-helix is unknown, they may play a critical role in highly conserved processes such as protein folding. © 1996 Wiley-Liss, Inc.     

Modulation of an IDP binding mechanism and rates by helix propensity and non-native interactions: association of HIF1α with CBP

Intrinsically disordered proteins that acquire their three dimensional structures only upon binding to their targets are very important in cellular signal regulation. While experimental studies have been made on the structures of both bound (structured) and unbound (disordered) states, less is known about the actual folding-binding transition. Coarse grained simulations using native-centric (i.e. Gō) potentials have been particularly useful in addressing this problem, given the large search space for IDP binding, but have well-known deficiencies in reproducing the unfolded state structure and dynamics. Here, we investigate the interaction of HIF1α with CBP using a hierarchy of coarse-grained models, in each case matching the binding affinity at 300 K to the experimental value. Starting from a pure Gō-like model based on the native structure of the complex we go on to consider a more realistic model of helix propensity in the HIF1α, and finally the effect of non-native interactions between binding partners. We find structural disorder (i.e."fuzziness") in the bound state of HIF1α in all models which is supported by the results of atomistic simulations. Correcting the over-stabilized helices in the unbound state gives rise to a more cooperative folding-binding transition (destabilizing partially bound intermediates). Adding non-native contacts lowers the free energy barrier for binding to an almost barrierless scenario, leading to higher binding/unbinding rates relative to the other models, in better agreement with the near diffusion-limited binding rates measured experimentally. Transition state structures for the three models are highly disordered, supporting a fly-casting mechanism for binding.     

Ultrafast Protein Folding in Membrane-Mimetic Environments

Proteins fold on timescales from hours to microseconds. In addition to protein size, sequence, and topology, the environment represents an equally important factor in determining folding speed. This is particularly relevant for proteins that require a lipid membrane or a membrane mimic to fold. However, only little is known about how properties of such a hydrophilic/hydrophobic interface modulate the folding landscape of membrane-interacting proteins. Here, we studied the influence of different membrane-mimetic micellar environments on the folding and unfolding kinetics of the helical-bundle protein Mistic. Devising a single-molecule fluorescence spectroscopy approach, we extracted folding and unfolding rates under equilibrium conditions and dissected the contributions from different detergent moieties to the free-energy landscape. While both polar and nonpolar moieties contribute to stability, they exert differential effects on the free-energy barrier: Hydrophobic burial stabilizes the folded state but not the transition state in reference to a purely aqueous environment; by contrast, zwitterionic headgroup moieties stabilize the folded state and, additionally, lower the free-energy barrier to accelerate the folding of Mistic to achieve ultrafast folding times down to 35 μs.     

Binding of Two Intrinsically Disordered Peptides to a Multi-Specific Protein: A Combined Monte Carlo and Molecular Dynamics Study

The unique ability of intrinsically disordered proteins (IDPs) to fold upon binding to partner molecules makes them functionally well-suited for cellular communication networks. For example, the folding-binding of different IDP sequences onto the same surface of an ordered protein provides a mechanism for signaling in a many-to-one manner. Here, we study the molecular details of this signaling mechanism by applying both Molecular Dynamics and Monte Carlo methods to S100B, a calcium-modulated homodimeric protein, and two of its IDP targets, p53 and TRTK-12. Despite adopting somewhat different conformations in complex with S100B and showing no apparent sequence similarity, the two IDP targets associate in virtually the same manner. As free chains, both target sequences remain flexible and sample their respective bound, natively -helical states to a small extent. Association occurs through an intermediate state in the periphery of the S100B binding pocket, stabilized by nonnative interactions which are either hydrophobic or electrostatic in nature. Our results highlight the importance of overall physical properties of IDP segments, such as net charge or presence of strongly hydrophobic amino acids, for molecular recognition via coupled folding-binding.     

Molecular recognition via coupled folding and binding in a TPR domain

The majority of known tetratricopeptide repeat (TPR) domains consist of three copies of the helix-turn-helix TPR motif, together with a seventh C-terminal helix. TPR domains function as protein-protein recognition modules in intracellular signalling. This function is exemplified by the TPR domain of protein phosphatase 5 (PP5), which binds to the C terminus of the chaperone protein Hsp90. Here, we report NMR and CD spectroscopic studies that reveal that this domain is largely unfolded at physiological temperatures, and that interaction with an MEEVD pentapeptide derived from Hsp90 stabilises a folded structure. This complex, coupled folding-binding mechanism is characterised further by its observed enthalpy change on binding (determined by isothermal titration calorimetry), which displays a markedly non-linear relationship with temperature. A nested Gibbs-Helmholtz model is used in a novel combined analysis of the CD and ITC data to determine separately the thermodynamic contributions of the intrinsic folding and binding events to the overall coupled process. The analysis shows that, despite the expected large entropic opposition to the folding process, a nearly equal favourable folding enthalpy means the net effect of coupled folding on the observed affinity is small across a broad range of temperature. We hypothesise that a coupled folding-binding mechanism is common in this class of domains.     

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.     

Structure of the leech protein saratin and characterization of its binding to collagen

The leech protein Saratin from Hirudo medicinalis prevents thrombocyte aggregation by interfering with the first binding step of the thrombocytes to collagen by binding to collagen. We solved the three-dimensional structure of the leech protein Saratin in solution and identified its collagen binding site by NMR titration experiments. The NMR structure of Saratin consists of one α-helix and a five-stranded β-sheet arranged in the topology ββαβββ. The C-terminal region, of about 20 amino acids in length, adopts no regular structure. NMR titration experiments with collagen peptides show that the collagen interaction of Saratin takes place in a kind of notch that is formed by the end of the α-helix and the β-sheet. NMR data-driven docking experiments to collagen model peptides were used to elucidate the putative binding mode of Saratin and collagen. Mainly, parts of the first and the end of the fifth β-strand, the loop connecting the α-helix and the third β-strand, and a short part of the loop connecting the fourth and fifth β-strand participate in binding.     

Distinct phases of free α-synuclein--a Monte Carlo study

The α-synuclein protein (αS), implicated in Parkinson's disease, shows conformational versatility. It aggregates into β-sheet-rich fibrils, occurs in helical membrane-bound forms, is disordered as a free monomer, and has recently been suggested to have a folded helical tetramer as its main physiological form. Here, we use implicit solvent all-atom Monte Carlo methods to explore the conformational ensemble sampled by the free αS monomer. We analyze secondary structure propensities, size, and topological properties and compare with existing experimental data. Our study suggests that free αS has two distinct phases. One phase has the expected disordered character. The other phase also shows large conformational variability. However, in this phase, the β-strand content is substantial, and the backbone fold shows statistical similarities with that in αS fibrils. Presence of this phase is consistent with data from low-temperature experiments. Conversion of disordered αS to this fibril-like form requires the crossing of a rather large apparent free-energy barrier.     

Modeling of Arylamide Helix Mimetics in the p53 Peptide Binding Site of hDM2 Suggests Parallel and Anti-Parallel Conformations Are Both Stable

The design of novel α-helix mimetic inhibitors of protein-protein interactions is of interest to pharmaceuticals and chemical genetics researchers as these inhibitors provide a chemical scaffold presenting side chains in the same geometry as an α-helix. This conformational arrangement allows the design of high affinity inhibitors mimicking known peptide sequences binding specific protein substrates. We show that GAFF and AutoDock potentials do not properly capture the conformational preferences of α-helix mimetics based on arylamide oligomers and identify alternate parameters matching solution NMR data and suitable for molecular dynamics simulation of arylamide compounds. Results from both docking and molecular dynamics simulations are consistent with the arylamides binding in the p53 peptide binding pocket. Simulations of arylamides in the p53 binding pocket of hDM2 are consistent with binding, exhibiting similar structural dynamics in the pocket as simulations of known hDM2 binders Nutlin-2 and a benzodiazepinedione compound. Arylamide conformations converge towards the same region of the binding pocket on the 20 ns time scale, and most, though not all dihedrals in the binding pocket are well sampled on this timescale. We show that there are two putative classes of binding modes for arylamide compounds supported equally by the modeling evidence. In the first, the arylamide compound lies parallel to the observed p53 helix. In the second class, not previously identified or proposed, the arylamide compound lies anti-parallel to the p53 helix.     

Characterization and comparison of two binding sites on obscurin for small ankyrin 1

Obscurin A, an ～720 kDa modular protein of striated muscles, binds to small ankyrin 1 (sAnk1, Ank 1.5), an integral protein of the sarcoplasmic reticulum, through two distinct carboxy-terminal sequences, Obsc(6316-6436) and Obsc(6236-6260). We hypothesized that these sequences differ in affinity but that they compete for the same binding site on sAnk1. We show that the sequence within Obsc(6316-6436) that binds to sAnk1 is limited to residues 6316-6345. Comparison of Obsc(6231-6260) to Obsc(6316-6345) reveals that Obsc(6316-6345) binds sAnk1 with an affinity (133 ± 43 nM) comparable to that of the Obsc(6316-6436) fusion protein, whereas Obsc(6231-6260) binds with lower affinity (384 ± 53 nM). Oligopeptides of each sequence compete for binding with both sites at half-maximal inhibitory concentrations consistent with the affinities measured directly. Five of six site-directed mutants of sAnk1 showed similar reductions in binding to each binding site on obscurin, suggesting that they dock to many of the same residues of sAnk1. Circular dichroism (CD) analysis of the synthetic oligopeptides revealed a 2-fold greater α-helical content in Obsc(6316-6346), ～35%, than Obsc(6231-6260,) ～17%. Using these data, structural prediction algorithms, and homology modeling, we predict that Obsc(6316-6345) contains a bent α-helix of 12 amino acids, flanked by short disordered regions, and that Obsc(6231-6260) has a short, N-terminal α-helix of 4-5 residues followed by a long disordered region. Our results are consistent with a model in which both sequences of obscurin differ significantly in structure but bind to the ankyrin-like repeat motifs of sAnk1 in a similar though not identical manner.     

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.     

A physical reference state unifies the structure-derived potential of mean force for protein folding and binding

Extracting knowledge-based statistical potential from known structures of proteins is proved to be a simple, effective method to obtain an approximate free-energy function. However, the different compositions of amino acid residues at the core, the surface, and the binding interface of proteins prohibited the establishment of a unified statistical potential for folding and binding despite the fact that the physical basis of the interaction (water-mediated interaction between amino acids) is the same. Recently, a physical state of ideal gas, rather than a statistically averaged state, has been used as the reference state for extracting the net interaction energy between amino acid residues of monomeric proteins. Here, we find that this monomer-based potential is more accurate than an existing all-atom knowledge-based potential trained with interfacial structures of dimers in distinguishing native complex structures from docking decoys (100% success rate vs. 52% in 21 dimer/trimer decoy sets). It is also more accurate than a recently developed semiphysical empirical free-energy functional enhanced by an orientation-dependent hydrogen-bonding potential in distinguishing native state from Rosetta docking decoys (94% success rate vs. 74% in 31 antibody-antigen and other complexes based on Z score). In addition, the monomer potential achieved a 93% success rate in distinguishing true dimeric interfaces from artificial crystal interfaces. More importantly, without additional parameters, the potential provides an accurate prediction of binding free energy of protein-peptide and protein-protein complexes (a correlation coefficient of 0.87 and a root-mean-square deviation of 1.76 kcal/mol with 69 experimental data points). This work marks a significant step toward a unified knowledge-based potential that quantitatively captures the common physical principle underlying folding and binding. A Web server for academic users, established for the prediction of binding free energy and the energy evaluation of the protein-protein complexes, may be found at http://theory.med.buffalo.edu.     

Mapping the Anopheles gambiae Odorant Binding Protein 1 (AgamOBP1) using modeling techniques,site directed mutagenesis,circular dichroism and ligand binding assays

The major malaria vector in Sub-Saharan Africa is the Anopheles gambiae mosquito. This species is a key target of malaria control measures. Mosquitoes find humans primarily through olfaction, yet the molecular mechanisms associated with host-seeking behavior remain largely unknown. To further understand the functionality of A. gambiae odorant binding protein 1 (AgamOBP1), we combined in silico protein structure modeling and site-directed mutagenesis to generate 16 AgamOBP1 protein analogues containing single point mutations of interest. Circular dichroism (CD) and ligand-binding assays provided data necessary to probe the effects of the point mutations on ligand binding and the overall structure of AgamOBP1. Far-UV CD spectra of mutated AgamOBP1 variants displayed both substantial decreases to ordered α-helix structure (up to22%) and increases to disordered α-helix structure(up to 15%) with only minimal changes in random coil (unordered) structure. In mutations Y54A, Y122A and W114Q, aromatic side chain removal from the binding site significantly reduced N-phenyl-1-naphthylamine binding. Several non-aromatic mutations (L15T, L19T, L58T, L58Y, M84Q, M84K, H111A, Y122A and L124T) elicited changes to protein conformation with subsequent effects on ligand binding. This study provides empirical evidence for the in silico predicted functions of specific amino acids in AgamOBP1 folding and ligand binding characteristics.     

Multifaceted recognition of vertebrate Rev1 by translesion polymerases ζ and κ

Translesion synthesis is a fundamental biological process that enables DNA replication across lesion sites to ensure timely duplication of genetic information at the cost of replication fidelity, and it is implicated in development of cancer drug resistance after chemotherapy. The eukaryotic Y-family polymerase Rev1 is an essential scaffolding protein in translesion synthesis. Its C-terminal domain (CTD), which interacts with translesion polymerase ζ through the Rev7 subunit and with polymerases κ, ι, and η in vertebrates through the Rev1-interacting region (RIR), is absolutely required for function. We report the first solution structures of the mouse Rev1 CTD and its complex with the Pol κ RIR, revealing an atypical four-helix bundle. Using yeast two-hybrid assays, we have identified a Rev7-binding surface centered at the α2-α3 loop and N-terminal half of α3 of the Rev1 CTD. Binding of the mouse Pol κ RIR to the Rev1 CTD induces folding of the disordered RIR peptide into a three-turn α-helix, with the helix stabilized by an N-terminal cap. RIR binding also induces folding of a disordered N-terminal loop of the Rev1 CTD into a β-hairpin that projects over the shallow α1-α2 surface and creates a deep hydrophobic cavity to interact with the essential FF residues juxtaposed on the same side of the RIR helix. Our combined structural and biochemical studies reveal two distinct surfaces of the Rev1 CTD that separately mediate the assembly of extension and insertion translesion polymerase complexes and provide a molecular framework for developing novel cancer therapeutics to inhibit translesion synthesis.     

Mapping the Anopheles gambiae Odorant Binding Protein 1 (AgamOBP1) using modeling techniques, site directed mutagenesis, circular dichroism and ligand binding assays

The major malaria vector in Sub-Saharan Africa is the Anopheles gambiae mosquito. This species is a key target of malaria control measures. Mosquitoes find humans primarily through olfaction, yet the molecular mechanisms associated with host-seeking behavior remain largely unknown. To further understand the functionality of A. gambiae odorant binding protein 1 (AgamOBP1), we combined in silico protein structure modeling and site-directed mutagenesis to generate 16 AgamOBP1 protein analogues containing single point mutations of interest. Circular dichroism (CD) and ligand-binding assays provided data necessary to probe the effects of the point mutations on ligand binding and the overall structure of AgamOBP1. Far-UV CD spectra of mutated AgamOBP1 variants displayed both substantial decreases to ordered α-helix structure (up to22%) and increases to disordered α-helix structure(up to 15%) with only minimal changes in random coil (unordered) structure. In mutations Y54A, Y122A and W114Q, aromatic side chain removal from the binding site significantly reduced N-phenyl-1-naphthylamine binding. Several non-aromatic mutations (L15T, L19T, L58T, L58Y, M84Q, M84K, H111A, Y122A and L124T) elicited changes to protein conformation with subsequent effects on ligand binding. This study provides empirical evidence for the in silico predicted functions of specific amino acids in AgamOBP1 folding and ligand binding characteristics.     

Partner-Mediated Polymorphism of an Intrinsically Disordered Protein

Intrinsically disordered proteins (IDPs) recognize their partners through molecular recognition elements (MoREs). The MoRE of the C-terminal intrinsically disordered domain of the measles virus nucleoprotein (N_TAIL) is partly pre-configured as an α-helix in the free form and undergoes α-helical folding upon binding to the X domain (XD) of the viral phosphoprotein. Beyond XD, N_TAIL also binds the major inducible heat shock protein 70 (hsp70). So far, no structural information is available for the N_TAIL/hsp70 complex. Using mutational studies combined with a protein complementation assay based on green fluorescent protein reconstitution, we have investigated both N_TAIL/XD and N_TAIL/hsp70 interactions. Although the same N_TAIL region binds the two partners, the binding mechanisms are different. Hsp70 binding is much more tolerant of MoRE substitutions than XD, and the majority of substitutions lead to an increased N_TAIL/hsp70 interaction strength. Furthermore, while an increased and a decreased α-helicity of the MoRE lead to enhanced and reduced interaction strength with XD, respectively, the impact on hsp70 binding is negligible, suggesting that the MoRE does not adopt an α-helical conformation once bound to hsp70. Here, by showing that the α-helical conformation sampled by the free form of the MoRE does not systematically commit it to adopt an α-helical conformation in the bound form, we provide an example of partner-mediated polymorphism of an IDP and of the relative insensitiveness of the bound structure to the pre-recognition state. The present results therefore contribute to shed light on the molecular mechanisms by which IDPs recognize different partners.