Finding of residues crucial for supersecondary structure formation

Stable Substructures in Proteins and How to Find Them Using Single-Molecule Force Spectroscopy

Three-dimensional structures of proteins are a source of fascination for scientists, due to the beauty of their sequence-encoded architectures and their highly diverse range of functions. These functions include acting as powerful catalysts, signal receptors, and versatile molecular motors as well as being building blocks for macroscopic structures, thus defining the shape of multicellular organisms. How protein structure is organized and assembled at the sub-nanometer scale is of great current interest. Specifically, the discovery of stable substructures and supersecondary structures has inspired research into their potential use in rationally engineered proteins with tailor-made properties. Here, we show how the search for stable substructures in large proteins can benefit from recent advances in single-molecule force spectroscopy using highly sensitive dual-beam optical tweezers. Our chapter provides a step-by-step description of the experimental workflow for (1) preparing proteins for mechanical interrogation, (2) interpreting the data, and (3) avoiding the most commonly occurring mistakes.

Contribution to the prediction of the fold code: application to immunoglobulin and flavodoxin cases

Background

Folding nucleus of globular proteins formation starts by the mutual interaction of a group of hydrophobic amino acids whose close contacts allow subsequent formation and stability of the 3D structure. These early steps can be predicted by simulation of the folding process through a Monte Carlo (MC) coarse grain model in a discrete space. We previously defined MIRs (Most Interacting Residues), as the set of residues presenting a large number of non-covalent neighbour interactions during such simulation. MIRs are good candidates to define the minimal number of residues giving rise to a given fold instead of another one, although their proportion is rather high, typically [15-20]% of the sequences. Having in mind experiments with two sequences of very high levels of sequence identity (up to 90%) but different folds, we combined the MIR method, which takes sequence as single input, with the “fuzzy oil drop” (FOD) model that requires a 3D structure, in order to estimate the residues coding for the fold. FOD assumes that a globular protein follows an idealised 3D Gaussian distribution of hydrophobicity density, with the maximum in the centre and minima at the surface of the “drop”. If the actual local density of hydrophobicity around a given amino acid is as high as the ideal one, then this amino acid is assigned to the core of the globular protein, and it is assumed to follow the FOD model. Therefore one obtains a distribution of the amino acids of a protein according to their agreement or rejection with the FOD model.

Results

We compared and combined MIR and FOD methods to define the minimal nucleus, or keystone, of two populated folds: immunoglobulin-like (Ig) and flavodoxins (Flav). The combination of these two approaches defines some positions both predicted as a MIR and assigned as accordant with the FOD model. It is shown here that for these two folds, the intersection of the predicted sets of residues significantly differs from random selection. It reduces the number of selected residues by each individual method and allows a reasonable agreement with experimentally determined key residues coding for the particular fold. In addition, the intersection of the two methods significantly increases the specificity of the prediction, providing a robust set of residues that constitute the folding nucleus.

Amino acid distribution rules predict protein fold: protein grammar for beta-strand sandwich-like structures

We present an alternative approach to protein 3D folding prediction based on determination of rules that specify distribution of "favorable" residues, that are mainly responsible for a given fold formation, and "unfavorable" residues, that are incompatible with that fold, in polypeptide sequences. The process of determining favorable and unfavorable residues is iterative. The starting assumptions are based on the general principles of protein structure formation as well as structural features peculiar to a protein fold under investigation. The initial assumptions are tested one-by-one for a set of all known proteins with a given structure. The assumption is accepted as a "rule of amino acid distribution" for the protein fold if it holds true for all, or near all, structures. If the assumption is not accepted as a rule, it can be modified to better fit the data and then tested again in the next step of the iterative search algorithm, or rejected. We determined the set of amino acid distribution rules for a large group of beta sandwich-like proteins characterized by a specific arrangement of strands in two beta sheets. It was shown that this set of rules is highly sensitive (~90%) and very specific (~99%) for identifying sequences of proteins with specified beta sandwich fold structure. The advantage of the proposed approach is that it does not require that query proteins have a high degree of homology to proteins with known structure. So long as the query protein satisfies residue distribution rules, it can be confidently assigned to its respective protein fold. Another advantage of our approach is that it allows for a better understanding of which residues play an essential role in protein fold formation. It may, therefore, facilitate rational protein engineering design.

Identification of an ideal-like fingerprint for a protein fold using overlapped conserved residues based approach

Design of an efficient fingerprint that detects homologous proteins at distant sequence identity has been a great challenge. This paper proposes a strategy to extract an ideal-like fingerprint with high specificity and sensitivity from a group of sequences related to a fold. The approach is devised based on the assumptions that the critical residues for a protein fold may be conserved in three aspects, i.e. sequence, structure, and intramolecular interaction, and embedded in secondary structures. We hypothesized that the residues satisfying such conditions simultaneously may work as an efficient fingerprint. This idea was tested on protein folds of various classes, such as beta-strand rich, alpha + beta proteins and alpha/beta proteins with discrete sequence similarities. The fingerprint for each fold was generated by selecting the overlapped conserved residues (OCR) from the conserved residues obtained using independent three alignment methods, i.e. multiple sequence alignment, structure-based alignment, and alignment based on the interstrand hydrogen-bonds. The OCR fingerprints showed more than 90% detection efficiency for all the folds tested and were identified to be almost the minimal fingerprints composed of only critical residues. This study is expected to provide an important conceptual improvement in the identification or design of ideal fingerprints for a protein fold.

Creating supersecondary structures with BuildBeta

BuildBeta is a feature of the ProteinShop software designed to thoroughly sample a protein conformational space given the protein's sequence of amino acids and secondary structure predictions. It targets proteins with beta sheets because they are particularly challenging to predict due to the complexity of sampling long-range strand pairings. Here we discuss some of the most difficult targets in the recent 9th Community Wide Experiment on the Critical Assessment of Techniques for Protein Structure Prediction (CASP9) and show how BuildBeta can leverage some of the most successful methods in the category "template-free modeling" by augmenting their sampling capabilities. We also discuss ongoing efforts to improve the quality of the supersecondary structures it generates.