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

Crowding by Poly(ethylene glycol) Destabilizes Chemotaxis Protein Y (CheY)

The prevailing understanding of various aspects of biochemical processes, including folding, stability, intermolecular interactions, and the binding of metals, substrates, and inhibitors, is derived from studies carried out under dilute and homogeneous conditions devoid of a crowding-related environment. The effect of crowding-induced modulation on the structure and stability of native and magnesium-dependent Chemotaxis Y (CheY), a bacterial signaling protein, was probed in the presence and absence of poly(ethylene glycol) (PEG). A combined analysis from circular dichroism, intrinsic and extrinsic fluorescence, and tryptophan fluorescence lifetime changes indicates that PEG perturbs the structure but leaves the thermal stability largely unchanged. Intriguingly, while the stability of the protein is enhanced in the presence of magnesium under dilute buffer conditions, PEG-induced crowding leads to reduced thermal stability in the presence of magnesium. Nuclear magnetic resonance (NMR) chemical shift perturbations and resonance broadening for a subset of residues indicate that PEG interacts specifically with a subset of hydrophilic and hydrophobic residues found predominantly in α helices, β strands, and in the vicinity of the metal-binding region. Thus, PEG prompted conformational perturbation, presumably provides a different situation for magnesium interaction, thereby perturbing the magnesium-prompted stability. In summary, our results highlight the dominance of enthalpic contributions between PEG and CheY via both hydrophilic and hydrophobic interactions, which can subtly affect the conformation, modulating the metal-protein interaction and stability, implying that in the context of cellular situation, structure, stability, and magnesium binding thermodynamics of CheY may be different from those measured in dilute solution.

The Folding Pathway of 6aJL2

One-third of the reported cases of light chain amyloidosis are related to the germ line λ6 family; remarkably, healthy individuals express this type of protein in just 2% of the peripheral blood and bone marrow B-cells. The appearance of the disease has been related to the inherent properties of this protein family. A recombinant representative model for λ6 proteins called 6aJL2 containing the amino acid sequence encoded by the 6a and JL2 germ line genes was previously designed and synthesized to study the properties of this family. Previous work on 6aJL2 suggested a simple two-state folding model at 25 °C; no intermediate could be identified either by kinetics or by fluorescence and circular dichroism equilibrium studies, although the presence of an intermediate that is populated at ∼2.4 M urea was suggested by size exclusion chromatography. In this study we employed classic equilibrium and kinetic experiments and analysis to elucidate the detailed folding mechanism of this protein. We identify species that are kinetically accessible and/or are populated at equilibrium. We describe the presence of intermediate and native-like species and propose a five-species folding mechanism at 25 °C at short incubation times, similar to and consistent with those observed in other proteins of this fold. The formation of intermediates in the mechanism of 6aJL2 is faster than that proposed for a V_κ light chain, which could be an important distinction in the amyloidogenic potential of both germ lines.

Computer design of microfluidic mixers for protein/RNA folding studies

Kinetic studies of biological macromolecules increasingly use microfluidic mixers to initiate and monitor reaction progress. A motivation for using microfluidic mixers is to reduce sample consumption and decrease mixing time to microseconds. Some applications, such as small-angle x-ray scattering, also require large (>10 micron) sampling areas to ensure high signal-to-noise ratios and to minimize parasitic scattering. Chaotic to marginally turbulent mixers are well suited for these applications because this class of mixers provides a good middle ground between existing laminar and turbulent mixers. In this study, we model various chaotic to marginally turbulent mixing concepts such as flow turning, flow splitting, and vortex generation using computational fluid dynamics for optimization of mixing efficiency and observation volume. Design iterations show flow turning to be the best candidate for chaotic/marginally turbulent mixing. A qualitative experimental test is performed on the finalized design with mixing of 10 M urea and water to validate the flow turning unsteady mixing concept as a viable option for RNA and protein folding studies. A comparison of direct numerical simulations (DNS) and turbulence models suggests that the applicability of turbulence models to these flow regimes may be limited.

Optimal allosteric stabilization sites using contact stabilization analysis

Proteins can be destabilized by a number of environmental factors such as temperature, pH, and mutation. The ability to subsequently restore function under these conditions by adding small molecule stabilizers, or by introducing disulfide bonds, would be a very powerful tool, but the physical principles that drive this stabilization are not well understood. The first problem lies is in choosing an appropriate binding site or disulfide bond location to best confer stability to the active site and restore function. Here, we present a general framework for predicting which allosteric binding sites correlate with stability in the active site. Using the Karanicolas-Brooks Gō-like model, we examine the dynamics of the enzyme β-glucuronidase using an Umbrella Sampling method to thoroughly sample the conformational landscape. Each intramolecular contact is assigned a score termed a "stabilization factor" that measures its correlation with structural changes in the active site. We have carried out this analysis for three different scaling strengths for the intramolecular contacts, and we examine how the calculated stabilization factors depend on the ensemble of destabilized conformations. We further examine a locally destabilized mutant of β-glucuronidase that has been characterized experimentally, and show that this brings about local changes in the stabilization factors. We find that the proximity to the active site is not sufficient to determine which contacts can confer active site stability. © 2016 Wiley Periodicals, Inc.

Folding of proteins with a flavodoxin-like architecture

The flavodoxin-like fold is a protein architecture that can be traced back to the universal ancestor of the three kingdoms of life. Many proteins share this α-β parallel topology and hence it is highly relevant to illuminate how they fold. Here, we review experiments and simulations concerning the folding of flavodoxins and CheY-like proteins, which share the flavodoxin-like fold. These polypeptides tend to temporarily misfold during unassisted folding to their functionally active forms. This susceptibility to frustration is caused by the more rapid formation of an α-helix compared to a β-sheet, particularly when a parallel β-sheet is involved. As a result, flavodoxin-like proteins form intermediates that are off-pathway to native protein and several of these species are molten globules (MGs). Experiments suggest that the off-pathway species are of helical nature and that flavodoxin-like proteins have a nonconserved transition state that determines the rate of productive folding. Folding of flavodoxin from Azotobacter vinelandii has been investigated extensively, enabling a schematic construction of its folding energy landscape. It is the only flavodoxin-like protein of which cotranslational folding has been probed. New insights that emphasize differences between in vivo and in vitro folding energy landscapes are emerging: the ribosome modulates MG formation in nascent apoflavodoxin and forces this polypeptide toward the native state.

Sequence-, structure-, and dynamics-based comparisons of structurally homologous CheY-like proteins

We recently introduced a physically based approach to sequence comparison, the property factor method (PFM). In the present work, we apply the PFM approach to the study of a challenging set of sequences-the bacterial chemotaxis protein CheY, the N-terminal receiver domain of the nitrogen regulation protein NT-NtrC, and the sporulation response regulator Spo0F. These are all response regulators involved in signal transduction. Despite functional similarity and structural homology, they exhibit low sequence identity. PFM sequence comparison demonstrates a statistically significant qualitative difference between the sequence of CheY and those of the other two proteins that is not found using conventional alignment methods. This difference is shown to be consonant with structural characteristics, using distance matrix comparisons. We also demonstrate that residues participating strongly in native contacts during unfolding are distributed differently in CheY than in the other two proteins. The PFM result is also in accord with dynamic simulation results of several types. Molecular dynamics simulations of all three proteins were carried out at several temperatures, and it is shown that the dynamics of CheY are predicted to differ from those of NT-NtrC and Spo0F. The predicted dynamic properties of the three proteins are in good agreement with experimentally determined B factors and with fluctuations predicted by the Gaussian network model. We pinpoint the differences between the PFM and traditional sequence comparisons and discuss the informatic basis for the ability of the PFM approach to detect physical differences between these sequences that are not apparent from traditional alignment-based comparison.

Using the folding landscapes of proteins to understand protein function

Proteins fold on a biologically-relevant timescale because of a funnel-shaped energy landscape. This landscape is sculpted through evolution by selecting amino-acid sequences that stabilize native interactions while suppressing stable non-native interactions that occur during folding. However, there is strong evolutionary selection for functional residues and these cannot be chosen to optimize folding. Their presence impacts the folding energy landscape in a variety of ways. Here, we survey the effects of functional residues on folding by providing several examples. We then review how such effects can be detected computationally and be used as assays for protein function. Overall, an understanding of how functional residues modulate folding should provide insights into the design of natural proteins and their homeostasis.

Evolutionary relationship of two ancient protein superfolds

Proteins are the molecular machines of the cell that fold into specific three-dimensional structures to fulfill their functions. To improve our understanding of how the structure and function of proteins arises, it is crucial to understand how evolution has generated the structural diversity we observe today. Classically, proteins that adopt different folds are considered to be nonhomologous. However, using state-of-the-art tools for homology detection, we found evidence of homology between proteins of two ancient and highly populated protein folds, the (βα)8-barrel and the flavodoxin-like fold. We detected a family of sequences that show intermediate features between both folds and determined what is to our knowledge the first representative crystal structure of one of its members, giving new insights into the evolutionary link of two of the earliest folds. Our findings contribute to an emergent vision where protein superfolds share common ancestry and encourage further approaches to complete the mapping of structure space onto sequence space.

Modulation of frustration in folding by sequence permutation

Folding of globular proteins can be envisioned as the contraction of a random coil unfolded state toward the native state on an energy surface rough with local minima trapping frustrated species. These substructures impede productive folding and can serve as nucleation sites for aggregation reactions. However, little is known about the relationship between frustration and its underlying sequence determinants. Chemotaxis response regulator Y (CheY), a 129-amino acid bacterial protein, has been shown previously to populate an off-pathway kinetic trap in the microsecond time range. The frustration has been ascribed to premature docking of the N- and C-terminal subdomains or, alternatively, to the formation of an unproductive local-in-sequence cluster of branched aliphatic side chains, isoleucine, leucine, and valine (ILV). The roles of the subdomains and ILV clusters in frustration were tested by altering the sequence connectivity using circular permutations. Surprisingly, the stability and buried surface area of the intermediate could be increased or decreased depending on the location of the termini. Comparison with the results of small-angle X-ray-scattering experiments and simulations points to the accelerated formation of a more compact, on-pathway species for the more stable intermediate. The effect of chain connectivity in modulating the structures and stabilities of the early kinetic traps in CheY is better understood in terms of the ILV cluster model. However, the subdomain model captures the requirement for an intact N-terminal domain to access the native conformation. Chain entropy and aliphatic-rich sequences play crucial roles in biasing the early events leading to frustration in the folding of CheY.

Sequence analysis on the information of folding initiation segments in ferredoxin-like fold proteins

BACKGROUND

While some studies have shown that the 3D protein structures are more conservative than their amino acid sequences, other experimental studies have shown that even if two proteins share the same topology, they may have different folding pathways. There are many studies investigating this issue with molecular dynamics or Go-like model simulations, however, one should be able to obtain the same information by analyzing the proteins' amino acid sequences, if the sequences contain all the information about the 3D structures. In this study, we use information about protein sequences to predict the location of their folding segments. We focus on proteins with a ferredoxin-like fold, which has a characteristic topology. Some of these proteins have different folding segments.

RESULTS

Despite the simplicity of our methods, we are able to correctly determine the experimentally identified folding segments by predicting the location of the compact regions considered to play an important role in structural formation. We also apply our sequence analyses to some homologues of each protein and confirm that there are highly conserved folding segments despite the homologues' sequence diversity. These homologues have similar folding segments even though the homology of two proteins' sequences is not so high.

CONCLUSION

Our analyses have proven useful for investigating the common or different folding features of the proteins studied.

Energetic frustrations in protein folding at residue resolution: a homologous simulation study of Im9 proteins

Energetic frustration is becoming an important topic for understanding the mechanisms of protein folding, which is a long-standing big biological problem usually investigated by the free energy landscape theory. Despite the significant advances in probing the effects of folding frustrations on the overall features of protein folding pathways and folding intermediates, detailed characterizations of folding frustrations at an atomic or residue level are still lacking. In addition, how and to what extent folding frustrations interact with protein topology in determining folding mechanisms remains unclear. In this paper, we tried to understand energetic frustrations in the context of protein topology structures or native-contact networks by comparing the energetic frustrations of five homologous Im9 alpha-helix proteins that share very similar topology structures but have a single hydrophilic-to-hydrophobic mutual mutation. The folding simulations were performed using a coarse-grained Gō-like model, while non-native hydrophobic interactions were introduced as energetic frustrations using a Lennard-Jones potential function. Energetic frustrations were then examined at residue level based on φ-value analyses of the transition state ensemble structures and mapped back to native-contact networks. Our calculations show that energetic frustrations have highly heterogeneous influences on the folding of the four helices of the examined structures depending on the local environment of the frustration centers. Also, the closer the introduced frustration is to the center of the native-contact network, the larger the changes in the protein folding. Our findings add a new dimension to the understanding of protein folding the topology determination in that energetic frustrations works closely with native-contact networks to affect the protein folding.

Balancing bond, nonbond, and gō-like terms in coarse grain simulations of conformational dynamics

Characterization of the protein conformational landscape remains a challenging problem, whether it concerns elucidating folding mechanisms, predicting native structures or modeling functional transitions. Coarse-grained molecular dynamics simulation methods enable exhaustive sampling of the energetic landscape at resolutions of biological interest. The general utility of structure-based models is reviewed along with their differing levels of approximation. Simple Gō models incorporate attractive native interactions and repulsive nonnative contacts, resulting in an ideal smooth landscape. Non-Gō coarse-grained models reduce the parameter set as needed but do not include bias to any desired native structure. While non-Gō models have achieved limited success in protein coarse-graining, they can be combined with native structured-based potentials to create a balanced and powerful force field. Recent applications of such Gō-like models have yielded insight into complex folding mechanisms and conformational transitions in large macromolecules. The accuracy and usefulness of reduced representations are also revealed to be a function of the mathematical treatment of the intrinsic bonded topology.

Assessment of local friction in protein folding dynamics using a helix cross-linker

Internal friction arising from local steric hindrance and/or the excluded volume effect plays an important role in controlling not only the dynamics of protein folding but also conformational transitions occurring within the native state potential well. However, experimental assessment of such local friction is difficult because it does not manifest itself as an independent experimental observable. Herein, we demonstrate, using the miniprotein trp-cage as a testbed, that it is possible to selectively increase the local mass density in a protein and hence the magnitude of local friction, thus making its effect directly measurable via folding kinetic studies. Specifically, we show that when a helix cross-linker, m-xylene, is placed near the most congested region of the trp-cage it leads to a significant decrease in both the folding rate (by a factor of 3.8) and unfolding rate (by a factor of 2.5 at 35 °C) but has little effect on protein stability. Thus, these results, in conjunction with those obtained with another cross-linked trp-cage and two uncross-linked variants, demonstrate the feasibility of using a nonperturbing cross-linker to help quantify the effect of internal friction. In addition, we estimate that a m-xylene cross-linker could lead to an increase in the roughness of the folding energy landscape by as much as 0.4-1.0k(B)T.

Sequence-dependent base-stacking stabilities guide tRNA folding energy landscapes

The folding of bacterial tRNAs with disparate sequences has been observed to proceed in distinct folding mechanisms despite their structural similarity. To explore the folding landscapes of tRNA, we performed ion concentration-dependent coarse-grained TIS model MD simulations of several E. coli tRNAs to compare their thermodynamic melting profiles to the classical absorbance spectra of Crothers and co-workers. To independently validate our findings, we also performed atomistic empirical force field MD simulations of tRNAs, and we compared the base-to-base distances from coarse-grained and atomistic MD simulations to empirical base-stacking free energies. We then projected the free energies to the secondary structural elements of tRNA, and we observe distinct, parallel folding mechanisms whose differences can be inferred on the basis of their sequence-dependent base-stacking stabilities. In some cases, a premature, nonproductive folding intermediate corresponding to the Ψ hairpin loop must backtrack to the unfolded state before proceeding to the folded state. This observation suggests a possible explanation for the fast and slow phases observed in tRNA folding kinetics.

Understanding the folding-function tradeoff in proteins

When an amino-acid sequence cannot be optimized for both folding and function, folding can get compromised in favor of function. To understand this tradeoff better, we devise a novel method for extracting the "function-less" folding-motif of a protein fold from a set of structurally similar but functionally diverse proteins. We then obtain the β-trefoil folding-motif, and study its folding using structure-based models and molecular dynamics simulations. CompariA protein sequence serves two purpson with the folding of wild-type β-trefoil proteins shows that function affects folding in two ways: In the slower folding interleukin-1β, binding sites make the fold more complex, increase contact order and slow folding. In the faster folding hisactophilin, residues which could have been part of the folding-motif are used for function. This reduces the density of native contacts in functional regions and increases folding rate. The folding-motif helps identify subtle structural deviations which perturb folding. These may then be used for functional annotation. Further, the folding-motif could potentially be used as a first step in the sequence design of function-less scaffold proteins. Desired function can then be engineered into these scaffolds.

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.

Native-state heterogeneity of β(2)-microglobulin as revealed by kinetic folding and real-time NMR experiments

The kinetic folding of β(2)-microglobulin from the acid-denatured state was investigated by interrupted-unfolding and interrupted-refolding experiments using stopped-flow double-jump techniques. In the interrupted unfolding, we first unfolded the protein by a pH jump from pH7.5 to pH2.0, and the kinetic refolding assay was carried out by the reverse pH jump by monitoring tryptophan fluorescence. Similarly, in the interrupted refolding, we first refolded the protein by a pH jump from pH2.0 to pH7.5 and used a guanidine hydrochloride (GdnHCl) concentration jump as well as the reverse pH jump as unfolding assays. Based on these experiments, the folding is represented by a parallel-pathway model, in which the molecule with the correct Pro32 cis isomer refolds rapidly with a rate constant of 5-6 s(-1), while the molecule with the Pro32 trans isomer refolds more slowly (pH7.5 and 25°C). At the last step of folding, the native-like trans conformer produced on the latter pathway isomerizes very slowly (0.001-0.002 s(-1)) into the native cis conformer. In the GdnHCl-induced unfolding assays in the interrupted refolding, the native-like trans conformer unfolded remarkably faster than the native cis conformer, and the direct GdnHCl-induced unfolding was also biphasic, indicating that the native-like trans conformer is populated at a significant level under the native condition. The one-dimensional NMR and the real-time NMR experiments of refolding further indicated that the population of the trans conformer increases up to 7-9% under a more physiological condition (pH7.5 and 37°C).

Multidomain protein solves the folding problem by multifunnel combined landscape: theoretical investigation of a Y-family DNA polymerase

Approximately three-fourths of eukaryotic proteins are composed of multiple independently folded domains. However, much of our understanding is based on single domain proteins or isolated domains whose studies directly lead to well-known energy landscape theory in which proteins fold by navigating through a funneled energy landscape toward native structure ensembles. The degrees of freedom for proteins with multiple domains are many orders of magnitude larger than that for single domain proteins. Now, the question arises: How do the multidomain proteins solve the "protein folding problem"? Here, we specifically address this issue by exploring the structure folding relationship of Sulfolobus solfataricus DNA polymerase IV (DPO4), a prototype Y-family DNA polymerase which contains a polymerase core consisting of a palm (P domain), a finger (F domain), and a thumb domain (T domain) in addition to a little finger domain (LF domain). The theoretical results are in good agreement with the experimental data and lead to several theoretical predictions. Finally, we propose that for rapid folding into well-defined conformations which carry out the biological functions, four-domain DPO4 employs a divide-and-conquer strategy, that is, combining multiple individual folding funnels into a single funnel (domains fold independently and then coalesce). In this way, the degrees of freedom for multidomain proteins are polynomial rather than exponential, and the conformational search process can be reduced effectively from a large to a smaller time scale.

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

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

Environment-dependent long-range structural distortion in a temperature-sensitive point mutant

Extensive environment-dependent rearrangement of the helix-turn-helix DNA recognition region and adjacent L-tryptophan binding pocket is reported in the crystal structure of dimeric E. coli trp aporepressor with point mutation Leu75Phe. In one of two subunits, the eight residues immediately C-terminal to the mutation are shifted forward in helical register by three positions, and the five following residues form an extrahelical loop accommodating the register shift. In contrast, the second subunit has wildtype-like conformation, as do both subunits in an isomorphous wildtype control structure. Treated together as an ensemble pair, the distorted and wildtype-like conformations of the mutant apoprotein agree more fully than either conformation alone with previously reported NOE measurements, and account more completely for its diverse biochemical and biophysical properties. The register-shifted segment Ile79-Ala80-Thr81-Ile82-Thr83 is helical in both conformations despite low helical propensity, suggesting an important structural role for the steric constraints imposed by β-branched residues in helical conformation.

Topological constraints: using RNA secondary structure to model 3D conformation, folding pathways, and dynamic adaptation

Accompanying recent advances in determining RNA secondary structure is the growing appreciation for the importance of relatively simple topological constraints, encoded at the secondary structure level, in defining the overall architecture, folding pathways, and dynamic adaptability of RNA. A new view is emerging in which tertiary interactions do not define RNA 3D structure, but rather, help select specific conformers from an already narrow, topologically pre-defined conformational distribution. Studies are providing fundamental insights into the nature of these topological constraints, how they are encoded by the RNA secondary structure, and how they interplay with other interactions, breathing new meaning to RNA secondary structure. New approaches have been developed that take advantage of topological constraints in determining RNA backbone conformation based on secondary structure, and a limited set of other, easily accessible constraints. Topological constraints are also providing a much-needed framework for rationalizing and describing RNA dynamics and structural adaptation. Finally, studies suggest that topological constraints may play important roles in steering RNA folding pathways. Here, we review recent advances in our understanding of topological constraints encoded by the RNA secondary structure.

Multi-scaled explorations of binding-induced folding of intrinsically disordered protein inhibitor IA3 to its target enzyme

Biomolecular function is realized by recognition, and increasing evidence shows that recognition is determined not only by structure but also by flexibility and dynamics. We explored a biomolecular recognition process that involves a major conformational change – protein folding. In particular, we explore the binding-induced folding of IA3, an intrinsically disordered protein that blocks the active site cleft of the yeast aspartic proteinase saccharopepsin (YPrA) by folding its own N-terminal residues into an amphipathic alpha helix. We developed a multi-scaled approach that explores the underlying mechanism by combining structure-based molecular dynamics simulations at the residue level with a stochastic path method at the atomic level. Both the free energy profile and the associated kinetic paths reveal a common scheme whereby IA3 binds to its target enzyme prior to folding itself into a helix. This theoretical result is consistent with recent time-resolved experiments. Furthermore, exploration of the detailed trajectories reveals the important roles of non-native interactions in the initial binding that occurs prior to IA3 folding. In contrast to the common view that non-native interactions contribute only to the roughness of landscapes and impede binding, the non-native interactions here facilitate binding by reducing significantly the entropic search space in the landscape. The information gained from multi-scaled simulations of the folding of this intrinsically disordered protein in the presence of its binding target may prove useful in the design of novel inhibitors of aspartic proteinases.

The intrinsically disordered peptide IA3 is the endogenous inhibitor for the enzyme named yeast aspartic proteinase saccharopepsin (YPrA). In the presence of YPrA, IA3 folds itself into an amphipathic helix that blocks the active site cleft of the enzyme. We developed a multi-scaled approach to explore the underlying mechanism of this binding-induced ordering transition. Our approach combines a structure-based molecular dynamics model at the residue level with a stochastic path method at the atomic level. Our simulations suggest that IA3 inhibits YPrA through an induced-fit mechanism where the enzyme (YPrA) induces conformational change of its inhibitor (IA3). This expands the definition of an induced-fit model from its original meaning that the binding of substrate (IA3) drives conformational change in the protein (YPrA). Our result is consistent with recent kinetic experiments and provides a microscopic explanation for the underlying mechanism. We also discuss the important roles of non-native interactions and backtracking. These results enrich our understanding of the enzyme-inhibition mechanism and may have value in the design of drugs.

The folding of single domain proteins--have we reached a consensus?

Rather than stressing the most recent advances in the field, this review highlights the fundamental topics where disagreement remains and where adequate experimental data are lacking. These topics include properties of the denatured state and the role of residual structure, the nature of the fundamental steps and barriers, the extent of pathway heterogeneity and non-native interactions, recent comparisons between theory and experiment, and finally, dynamical properties of the folding reaction.