Mapping the Protein Fold Universe Using the CamTube Force Field in Molecular Dynamics Simulations

It has been recently shown that the coarse-graining of the structures of polypeptide chains as self-avoiding tubes can provide an effective representation of the conformational space of proteins. In order to fully exploit the opportunities offered by such a \u2018tube model\u2019 approach, we present here a strategy to combine it with molecular dynamics simulations. This strategy is based on the incorporation of the \u2018CamTube\u2019 force field into the Gromacs molecular dynamics package. By considering the case of a 60-residue polyvaline chain, we show that CamTube molecular dynamics simulations can comprehensively explore the conformational space of proteins. We obtain this result by a 20 \u3bcs metadynamics simulation of the polyvaline chain that recapitulates the currently known protein fold universe. We further show that, if residue-specific interaction potentials are added to the CamTube force field, it is possible to fold a protein into a topology close to that of its native state. These results illustrate how the CamTube force field can be used to explore efficiently the universe of protein folds with good accuracy and very limited computational cost

Structural characterization of the interaction of α-synuclein nascent chains with the ribosomal surface and trigger factor

The ribosome is increasingly becoming recognized as a key hub for integrating quality control processes associated with protein biosynthesis and cotranslational folding (CTF). The molecular mechanisms by which these processes take place, however, remain largely unknown, in particular in the case of intrinsically disordered proteins (IDPs). To address this question, we studied at a residue-specific level the structure and dynamics of ribosome-nascent chain complexes (RNCs) of α-synuclein (αSyn), an IDP associated with Parkinson’s disease (PD). Using solution-state nuclear magnetic resonance (NMR) spectroscopy and coarse-grained molecular dynamics (MD) simulations, we find that, although the nascent chain (NC) has a highly disordered conformation, its N-terminal region shows resonance broadening consistent with interactions involving specific regions of the ribosome surface. We also investigated the effects of the ribosome-associated molecular chaperone trigger factor (TF) on αSyn structure and dynamics using resonance broadening to define a footprint of the TF–RNC interactions. We have used these data to construct structural models that suggest specific ways by which emerging NCs can interact with the biosynthesis and quality control machinery

Toward an accurate prediction of inter-residue distances in proteins using 2D recursive neural networks

BACKGROUND: Protein inter-residue contact maps provide a translation and rotation invariant topological representation of a protein. They can be used as an intermediary step in protein structure predictions. However, the prediction of contact maps represents an unbalanced problem as far fewer examples of contacts than non-contacts exist in a protein structure. In this study we explore the possibility of completely eliminating the unbalanced nature of the contact map prediction problem by predicting real-value distances between residues. Predicting full inter-residue distance maps and applying them in protein structure predictions has been relatively unexplored in the past. RESULTS: We initially demonstrate that the use of native-like distance maps is able to reproduce 3D structures almost identical to the targets, giving an average RMSD of 0.5Å. In addition, the corrupted physical maps with an introduced random error of ±6Å are able to reconstruct the targets within an average RMSD of 2Å. After demonstrating the reconstruction potential of distance maps, we develop two classes of predictors using two-dimensional recursive neural networks: an ab initio predictor that relies only on the protein sequence and evolutionary information, and a template-based predictor in which additional structural homology information is provided. We find that the ab initio predictor is able to reproduce distances with an RMSD of 6Å, regardless of the evolutionary content provided. Furthermore, we show that the template-based predictor exploits both sequence and structure information even in cases of dubious homology and outperforms the best template hit with a clear margin of up to 3.7Å. Lastly, we demonstrate the ability of the two predictors to reconstruct the CASP9 targets shorter than 200 residues producing the results similar to the state of the machine learning art approach implemented in the Distill server. CONCLUSIONS: The methodology presented here, if complemented by more complex reconstruction protocols, can represent a possible path to improve machine learning algorithms for 3D protein structure prediction. Moreover, it can be used as an intermediary step in protein structure predictions either on its own or complemented by NMR restraints

NMR characterization of the conformational fluctuations of the human lymphocyte function-associated antigen-1 I-domain

Lymphocyte function-associated antigen-1 (LFA-1) is an integrin protein that transmits information across the plasma membrane through the so-called inside-out and outside-in signaling mechanisms. To investigate these mechanisms, we carried out an NMR analysis of the dynamics of the LFA-1 I-domain, which has enabled us to characterize the motions of this domain on a broad range of timescales. We studied first the internal motions on the nanosecond timescale by spin relaxation measurements and model-free analysis. We then extended this analysis to the millisecond timescale motions by measuring 15N-1H residual dipolar couplings of the backbone amide groups. We analyzed these results in the context of the three major conformational states of the I-domain using their corresponding X-ray crystallographic structures. Our results highlight the importance of the low-frequency motions of the LFA-1 I-domain in the inactive apo-state. We found in particular that α-helix 7 is in a position in the apo-closed state that cannot be fully described by any of the existing X-ray structures, as it appears to be in dynamic exchange between different conformations. This type of motion seems to represent an inherent property of the LFA-1 I-domain and might be relevant for controlling the access to the allosteric binding pocket, as well as for the downward displacement of α-helix 7 that is required for the activation of LFA-1