StructureDistiller: Structural relevance scoring identifies the most informative entries of a contact map

Protein folding and structure prediction are two sides of the same coin. Contact maps and the related techniques of constraint-based structure reconstruction can be considered as unifying aspects of both processes. We present the Structural Relevance (SR) score which quantifies the information content of individual contacts and residues in the context of the whole native structure. The physical process of protein folding is commonly characterized with spatial and temporal resolution: some residues are Early Folding while others are Highly Stable with respect to unfolding events. We employ the proposed SR score to demonstrate that folding initiation and structure stabilization are subprocesses realized by distinct sets of residues. The example of cytochrome c is used to demonstrate how StructureDistiller identifies the most important contacts needed for correct protein folding. This shows that entries of a contact map are not equally relevant for structural integrity. The proposed StructureDistiller algorithm identifies contacts with the highest information content; these entries convey unique constraints not captured by other contacts. Identification of the most informative contacts effectively doubles resilience toward contacts which are not observed in the native contact map. Furthermore, this knowledge increases reconstruction fidelity on sparse contact maps significantly by 0.4 Å.

Unbiased Atomistic Insight in the Mechanisms and Solvent Role for Globular Protein Dimer Dissociation

Association and dissociation of proteins are fundamental processes in nature. While this process is simple to understand conceptually, the details of the underlying mechanism and role of the solvent are poorly understood. Here we investigate the mechanism and solvent role for the dissociation of the hydrophilic β-lactoglobulin dimer by employing transition path sampling. Analysis of the sampled path ensembles indicates that dissociation (and association) occurs via a variety of mechanisms: 1) a direct aligned dissociation 2) a hopping and rebinding transition followed by unbinding 3) a sliding transition before unbinding. Reaction coordinate and transition state analysis predicts that, besides native contact and vicinity salt-bridge interactions, solvent degrees of freedom play an important role in the dissociation process. Analysis of the structure and dynamics of the solvent molecules reveals that the dry native interface induces enhanced populations of both disordered hydration water and hydration water with higher tetrahedrality, mainly nearby hydrophobic residues. Bridging waters, hydrogen bonded to both proteins, support contacts, and exhibit a faster decay and reorientation dynamics in the transition state than in the native state interface, which renders the proteins more mobile and assists in rebinding. While not exhaustive, our sampling of rare un-biased reactive molecular dynamics trajectories shows in full detail how proteins can dissociate via complex pathways including (multiple) rebinding events. The atomistic insight obtained assists in further understanding and control of the dynamics of protein-protein interaction including the role of solvent.PACS numbers:

Influences of heterogeneous native contact energy and many-body interactions on the prediction of protein folding mechanisms

Combining heterogenous native contact energies and many-body interactions could improve the prediction of Brønsted plots using a structure-based model.

Sequence, structure, and cooperativity in folding of elementary protein structural motifs

Residue-level unfolding of two helix-turn-helix proteins—one naturally occurring and one de novo designed—is reconstructed from multiple sets of site-specific 13C isotopically edited infrared (IR) and circular dichroism (CD) data using Ising-like statistical-mechanical models. Several model variants are parameterized to test the importance of sequence-specific interactions (approximated by Miyazawa–Jernigan statistical potentials), local structural flexibility (derived from the ensemble of NMR structures), interhelical hydrogen bonds, and native contacts separated by intervening disordered regions (through the Wako–Saitô–Muñoz–Eaton scheme, which disallows such configurations). The models are optimized by directly simulating experimental observables: CD ellipticity at 222 nm for model proteins and their fragments and 13C-amide I′ bands for multiple isotopologues of each protein. We find that data can be quantitatively reproduced by the model that allows two interacting segments flanking a disordered loop (double sequence approximation) and incorporates flexibility in the native contact maps, but neither sequence-specific interactions nor hydrogen bonds are required. The near-identical free energy profiles as a function of the global order parameter are consistent with expected similar folding kinetics for nearly identical structures. However, the predicted folding mechanism for the two motifs is different, reflecting the order of local stability. We introduce free energy profiles for “experimental” reaction coordinates—namely, the degree of local folding as sensed by site-specific 13C-edited IR, which highlight folding heterogeneity and contrast its overall, average description with the detailed, local picture.