Protein intrachain contact prediction with most interacting residues (MIR)

Abstract

The transition state ensemble during the folding process of globular proteins occurs when a sufficient number of intrachain contacts are formed, mainly, but not exclusively, due to hydrophobic interactions. These contacts are related to the folding nucleus, and they contribute to the stability of the native structure, although they may disappear after the energetic barrier of transition states has been passed. A number of structure and sequence analyses, as well as protein engineering studies, have shown that the signature of the folding nucleus is surprisingly present in the native three-dimensional structure, in the form of closed loops, and also in the early folding events. These findings support the idea that the residues of the folding nucleus become buried in the very first folding events, therefore helping the formation of closed loops that act as anchor structures, speed up the process, and overcome the Levinthal paradox. We present here a review of an algorithm intended to simulate in a discrete space the early steps of the folding process. It is based on a Monte Carlo simulation where perturbations, or moves, are randomly applied to residues within a sequence. In contrast with many technically similar approaches, this model does not intend to fold the protein but to calculate the number of non-covalent neighbors of each residue, during the early steps of the folding process. Amino acids along the sequence are categorized as most interacting residues (MIRs) or least interacting residues. The MIR method can be applied under a variety of circumstances. In the cases tested thus far, MIR has successfully identified the exact residue whose mutation causes a switch in conformation. This follows with the idea that MIR identifies residues that are important in the folding process. Most MIR positions correspond to hydrophobic residues; correspondingly, MIRs have zero or very low accessible surface area. Alongside the review of the MIR method, we present a new postprocessing method called smoothed MIR (SMIR), which refines the original MIR method by exploiting the knowledge of residue hydrophobicity. We review known results and present new ones, focusing on the ability of MIR to predict structural changes, secondary structure, and the improved precision with the SMIR method.

Contribution to the Prediction of the Fold Code: Application to Immunoglobulin and Flavodoxin Cases

“…In order to get rid of statistical fluctuations, one compares the NCN smoothed distribution of these two domains (available at http://sprouts.rpbs.univ-paris-diderot.fr/mir.html ). Peak in the NCN distribution of PDB code 2LHC at L45 for GA disappears when mutated to Tyrosine in GB while the global shape of the distribution is maintained [ 75 , 76 ]. Combination of MIR and FOD methods reinforces the hypothesis of stability: in GA, two positions are both registered as MIR and H+2, L45 and I49, while in GB, there are three, at L7, F30 and F52.…”