Significance: Quantitative understanding of how coupled binding and folding occurs for intrinsically disordered molecular recognition element (MoRE) critical for function is still challenging. By developing an integrated approach, we provided physical quantification of folding and binding energy landscapes of a MoRE at atomistic level. Our predictions are in remarkable agreements with the experiments, and lead to the recognition mechanism via conformational selection, followed by induced folding. We provided an explanation for the underlying connections among “downhill folding,” “molten globule,” and “intrinsic disorder.” We proposed a “kinetic divide-and-conquer” mechanism to understand the high specificity without high affinity in intrinsic disordered protein binding.

Many proteins contain multiple folded domains separated by flexible linkers, and the ability to describe the structure and conformational heterogeneity of such flexible systems pushes the limits of structural biology. Using the three-domain protein TIA-1 as an example, we here combine coarse-grained molecular dynamics simulations with previously measured small-angle scattering data to study the conformation of TIA-1 in solution. We show that while the coarse-grained potential (Martini) in itself leads to too compact conformations, increasing the strength of protein-water interactions results in ensembles that are in very good agreement with experiments. We show how these ensembles can be refined further using a Bayesian/Maximum Entropy approach, and examine the robustness to errors in the energy function. In particular we find that as long as the initial simulation is relatively good, reweighting against experiments is very robust. We also study the relative information in X-ray and neutron scattering experiments and find that refining against the SAXS experiments leads to improvement in the SANS data. Our results suggest a general strategy for studying the conformation of multi-domain proteins in solution that combines coarse-grained simulations with small-angle X-ray scattering data that are generally most easy to obtain. These results may in turn be used to design further small-angle neutron scattering experiments that exploit contrast variation through ¹H/²H isotope substitutions. Many proteins contain multiple folded domains separated by flexible linkers, and in order to understand how such multi-domain proteins function, we need to be able to describe how these domains are oriented in space. We have used the three-domain protein TIA-1 as an example to combine molecular simulations with biophysical experiments to describe the structural and dynamical properties of a multi-domain protein. We show that while standard simulations do not lead to good agreement with the experimental data, we can improve the agreement substantially by tuning a single parameter in the model that describes the interaction between protein and water. We can gain further information about the system by a more direct integration of the data, and we find that we can provide a detailed and robust description of the relative location of the different domains in TIA-1. The method is general and will be useful to study the relationship between structure, dynamics and function in multi-domain proteins in other systems.

The exchange of metabolites between the mitochondrial matrix and the cytosol depends on β-barrel channels in the outer membrane and α-helical carrier proteins in the inner membrane. The essential translocase of the inner membrane (TIM) chaperones escort these proteins through the intermembrane space, but the structural and mechanistic details remain elusive. We have used an integrated structural biology approach to reveal the functional principle of TIM chaperones. Multiple clamp-like binding sites hold the mitochondrial membrane proteins in a translocation-competent elongated form, thus mimicking characteristics of co-translational membrane insertion. The bound preprotein undergoes conformational dynamics within the chaperone binding clefts, pointing to a multitude of dynamic local binding events. Mutations in these binding sites cause cell death or growth defects associated with impairment of carrier and β-barrel protein biogenesis. Our work reveals how a single mitochondrial "transfer-chaperone" system is able to guide α-helical and β-barrel membrane proteins in a "nascent chain-like" conformation through a ribosome-free compartment.

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.