Analysis of Integrin αIIb Subunit Dynamics Reveals Long-Range Effects of Missense Mutations on Calf Domains.

Integrin αIIbß3, a glycoprotein complex expressed at the platelet surface, is involved in platelet aggregation and contributes to primary haemostasis. Several integrin αIIbß3 polymorphisms prevent the aggregation that causes haemorrhagic syndromes, such as Glanzmann thrombasthenia (GT). Access to 3D structure allows understanding the structural effects of polymorphisms related to GT. In a previous analysis using Molecular Dynamics (MD) simulations of αIIbCalf-1 domain structure, it was observed that GT associated with single amino acid variation affects distant loops, but not the mutated position. In this study, experiments are extended to Calf-1, Thigh, and Calf-2 domains. Two loops in Calf-2 are unstructured and therefore are modelled expertly using biophysical restraints. Surprisingly, MD revealed the presence of rigid zones in these loops. Detailed analysis with structural alphabet, the Proteins Blocks (PBs), allowed observing local changes in highly flexible regions. The variant P741R located at C-terminal of Calf-1 revealed that the Calf-2 presence did not affect the results obtained with isolated Calf-1 domain. Simulations for Calf-1 + Calf-2, and Thigh + Calf-1 variant systems are designed to comprehend the impact of five single amino acid variations in these domains. Distant conformational changes are observed, thus highlighting the potential role of allostery in the structural basis of GT.

Investigating the Product Profiles and Structural Relationships of New Levansucrases with Conventional and Non-Conventional Substrates.

The synthesis of complex oligosaccharides is desired for their potential as prebiotics, and their role in the pharmaceutical and food industry. Levansucrase (LS, EC 2.4.1.10), a fructosyl-transferase, can catalyze the synthesis of these compounds. LS acquires a fructosyl residue from a donor molecule and performs a non-Lenoir transfer to an acceptor molecule, via ß-(2→6)-glycosidic linkages. Genome mining was used to uncover new LS enzymes with increased transfructosylating activity and wider acceptor promiscuity, with an initial screening revealing five LS enzymes. The product profiles and activities of these enzymes were examined after their incubation with sucrose. Alternate acceptor molecules were also incubated with the enzymes to study their consumption. LSs from Gluconobacter oxydans and Novosphingobium aromaticivorans synthesized fructooligosaccharides (FOSs) with up to 13 units in length. Alignment of their amino acid sequences and substrate docking with homology models identified structural elements causing differences in their product spectra. Raffinose, over sucrose, was the preferred donor molecule for the LS from Vibrio natriegens, N. aromaticivorans, and Paraburkolderia graminis. The LSs examined were found to have wide acceptor promiscuity, utilizing monosaccharides, disaccharides, and two alcohols to a high degree.

Investigation of the impact of PTMs on the protein backbone conformation.

Post-translational modifications (PTMs) are known to play a critical role in the regulation of protein functions. Their impact on protein structures and their link to disorder regions have already been spotted in the past decade. Nonetheless, the high diversity of PTM types and the multiple schemes of protein modifications (multiple PTMs, of different types, at different time, etc.) make difficult the direct confrontation of PTM annotations and protein structure data. Therefore, we analyzed the impact of the residue modifications on the protein structures at the local level. Thanks to a dedicated structure database, namely PTM-SD, a large screen of PTMs have been done and analyzed at local protein conformation levels using the structural alphabet protein blocks (PBs). We investigated the relation between PTMs and the backbone conformation of modified residues, of their local environment, and at the level of the complete protein structure. The two main PTM types (N-glycosylation and phosphorylation) have been studied in non-redundant datasets, and then four different proteins were focused, covering three types of PTMs: N-glycosylation in renin endopeptidase and liver carboxylesterase, phosphorylation in cyclin-dependent kinase 2 (CDK2), and methylation in actin. We observed that PTMs could either stabilize or destabilize the backbone structure, at a local and global scale, and that these effects depend on the PTM types.

Shaking the ß-Bulges.

ß-bulges are irregularities inside the ß-sheets. They represent more than 3 percent of the protein residues, i.e., they are as frequent as 3.10 helices. In terms of evolution, ß-bulges are not more conserved than any other local protein conformations within homologous protein structures. In a first of its kind study, we have investigated the dynamical behaviour of ß-bulges using the largest known set of protein molecular dynamics simulations. We observed that more than 50 percent of the existing ß-bulges in protein crystal structures remained stable during dynamics while more than1/6th were not stable at all and disappeared entirely. Surprisingly, 1.1 percent of ß-bulges that appeared remained stable. ß-bulges have been categorized in different subtypes. The most common ß-bulges' types are the smallest insertion in ß-strands (namely AC and AG); they are found as stable as the whole ß-bulges dataset. Low occurring types (namely PC and AS), that have the largest insertions, are significantly more stable than expected. Thus, this pioneer study allowed to precisely quantify the stability of the ß-bulges, demonstrating their structural robustness, with few unexpected cases raising structural questions.

Discrete analysis of camelid variable domains: sequences, structures, and in-silico structure prediction.

Antigen binding by antibodies requires precise orientation of the complementarity- determining region (CDR) loops in the variable domain to establish the correct contact surface. Members of the family Camelidae have a modified form of immunoglobulin gamma (IgG) with only heavy chains, called Heavy Chain only Antibodies (HCAb). Antigen binding in HCAbs is mediated by only three CDR loops from the single variable domain (VHH) at the N-terminus of each heavy chain. This feature of the VHH, along with their other important features, e.g., easy expression, small size, thermo-stability and hydrophilicity, made them promising candidates for therapeutics and diagnostics. Thus, to design better VHH domains, it is important to thoroughly understand their sequence and structure characteristics and relationship. In this study, sequence characteristics of VHH domains have been analysed in depth, along with their structural features using innovative approaches, namely a structural alphabet. An elaborate summary of various studies proposing structural models of VHH domains showed diversity in the algorithms used. Finally, a case study to elucidate the differences in structural models from single and multiple templates is presented. In this case study, along with the above-mentioned aspects of VHH, an exciting view of various factors in structure prediction of VHH, like template framework selection, is also discussed.

iPBAvizu: a PyMOL plugin for an efficient 3D protein structure superimposition approach.

BACKGROUND: Protein 3D structure is the support of its function. Comparison of 3D protein structures provides insight on their evolution and their functional specificities and can be done efficiently via protein structure superimposition analysis. Multiple approaches have been developed to perform such task and are often based on structural superimposition deduced from sequence alignment, which does not take into account structural features. Our methodology is based on the use of a Structural Alphabet (SA), i.e. a library of 3D local protein prototypes able to approximate protein backbone. The interest of a SA is to translate into 1D sequences into the 3D structures. RESULTS: We used Protein blocks (PB), a widely used SA consisting of 16 prototypes, each representing a conformation of the pentapeptide skeleton defined in terms of dihedral angles. Proteins are described using PB from which we have previously developed a sequence alignment procedure based on dynamic programming with a dedicated PB Substitution Matrix. We improved the procedure with a specific two-step search: (i) very similar regions are selected using very high weights and aligned, and (ii) the alignment is completed (if possible) with less stringent parameters. Our approach, iPBA, has shown to perform better than other available tools in benchmark tests. To facilitate the usage of iPBA, we designed and implemented iPBAvizu, a plugin for PyMOL that allows users to run iPBA in an easy way and analyse protein superimpositions. CONCLUSIONS: iPBAvizu is an implementation of iPBA within the well-known and widely used PyMOL software. iPBAvizu enables to generate iPBA alignments, create and interactively explore structural superimposition, and assess the quality of the protein alignments.

In silico prediction of protein flexibility with local structure approach.

Flexibility is an intrinsic essential feature of protein structures, directly linked to their functions. To this day, most of the prediction methods use the crystallographic data (namely B-factors) as the only indicator of protein's inner flexibility and predicts them as rigid or flexible. PredyFlexy stands differently from other approaches as it relies on the definition of protein flexibility (i) not only taken from crystallographic data, but also (ii) from Root Mean Square Fluctuation (RMSFs) observed in Molecular Dynamics simulations. It also uses a specific representation of protein structures, named Long Structural Prototypes (LSPs). From Position-Specific Scoring Matrix, the 120 LSPs are predicted with a good accuracy and directly used to predict (i) the protein flexibility in three categories (flexible, intermediate and rigid), (ii) the normalized B-factors, (iii) the normalized RMSFs, and (iv) a confidence index. Prediction accuracy among these three classes is equivalent to the best two class prediction methods, while the normalized B-factors and normalized RMSFs have a good correlation with experimental and in silico values. Thus, PredyFlexy is a unique approach, which is of major utility for the scientific community. It support parallelization features and can be run on a local cluster using multiple cores.

Protein flexibility in the light of structural alphabets.

Protein structures are valuable tools to understand protein function. Nonetheless, proteins are often considered as rigid macromolecules while their structures exhibit specific flexibility, which is essential to complete their functions. Analyses of protein structures and dynamics are often performed with a simplified three-state description, i.e., the classical secondary structures. More precise and complete description of protein backbone conformation can be obtained using libraries of small protein fragments that are able to approximate every part of protein structures. These libraries, called structural alphabets (SAs), have been widely used in structure analysis field, from definition of ligand binding sites to superimposition of protein structures. SAs are also well suited to analyze the dynamics of protein structures. Here, we review innovative approaches that investigate protein flexibility based on SAs description. Coupled to various sources of experimental data (e.g., B-factor) and computational methodology (e.g., Molecular Dynamic simulation), SAs turn out to be powerful tools to analyze protein dynamics, e.g., to examine allosteric mechanisms in large set of structures in complexes, to identify order/disorder transition. SAs were also shown to be quite efficient to predict protein flexibility from amino-acid sequence. Finally, in this review, we exemplify the interest of SAs for studying flexibility with different cases of proteins implicated in pathologies and diseases.