Methods for the refinement of protein structure 3D models

The refinement of predicted 3D protein models is crucial in bringing them closer towards experimental accuracy for further computational studies. Refinement approaches can be divided into two main stages: The sampling and scoring stages. Sampling strategies, such as the popular Molecular Dynamics (MD)-based protocols, aim to generate improved 3D models. However, generating 3D models that are closer to the native structure than the initial model remains challenging, as structural deviations from the native basin can be encountered due to force-field inaccuracies. Therefore, different restraint strategies have been applied in order to avoid deviations away from the native structure. For example, the accurate prediction of local errors and/or contacts in the initial models can be used to guide restraints. MD-based protocols, using physics-based force fields and smart restraints, have made significant progress towards a more consistent refinement of 3D models. The scoring stage, including energy functions and Model Quality Assessment Programs (MQAPs) are also used to discriminate near-native conformations from non-native conformations. Nevertheless, there are often very small differences among generated 3D models in refinement pipelines, which makes model discrimination and selection problematic. For this reason, the identification of the most native-like conformations remains a major challenge

MERMAID: dedicated web server to prepare and run coarse-grained membrane protein dynamics

Atomistic molecular dynamics simulations of membrane proteins have been shown to be extremely useful for characterizing the molecular features underlying their function, but require high computational power, limiting the understanding of complex events in membrane proteins, e.g. ion channels gating, GPCRs activation. To overcome this issue, it has been shown that coarse-grained approaches, although requiring less computational power, are still capable of correctly describing molecular events underlying big conformational changes in biological systems. Here, we present the Martini coarse-grained membrane protein dynamics (MERMAID), a publicly available web interface that allows the user to prepare and run coarse-grained molecular dynamics (CGMD) simulations and to analyse the trajectories

Homotypic and Heterotypic Self-Assembly of Claudin Family of Tight Junction Proteins

Tight junctions, found in all epithelial cells, are selective barriers that restrict the diffusion of ions and molecules within the paracellular space. They are important for maintaining cell polarity as well as for compartmentalization and establishing homeostasis within the human body. Tight junctions are comprised of complex protein assemblies. Within this protein assemble are a family of transmembrane proteins known as claudins that play a crucial role in establishing the tight junction network. Claudins are also influential in controlling the tight junction permeability. When mutations or malfunctioning occurs within a claudin gene, tight junction function is impaired. Disruption of their function is associated with a variety of human conditions, such as brain disease, deafness, renal failure, and various cancers. A deeper understanding of claudins and of their contribution towards tight junction’s function will provide researchers additional insight as to how to eventually approach creating therapeutics to treat tight junction-related diseases. In this thesis, homotypic and heterotypic cis self-assembly of claudin-claudin interactions were studied for both classic and non-classic claudins (-2, -4, -11, -14, -16, -18, -19, and -23). Homology modeling was utilized to generate structures for each of the eight claudins studied, which were then equilibrated and refined in a DOPC (1,2-Dioleoyl-sn-glycero-3-phosphocholine) lipid bilayer system. Consequently, self-assemble simulations were carried out to study the cis interaction in either homotypic or heterotypic fashion. Each self-assembly system contained 72 monomers and was simulated for 4 µs. Results showed aggregation of claudin monomers into strand-like assemblies, which were then analyzed through a dimer distribution and orientation analysis. Four dimer types (dimers A, B, C, D) were identified and dimer populations were calculated in each of the claudin self-assembly systems. An additional new analysis method developed by a colleague and still in the refining phase is also introduced and discussed. Results for a single test system are discussed, but nonetheless provide an alternative means of analyzing dimers, representing them through various energy state profiles rather than of population density probabilities

Current state-of-the-art of the research conducted in mapping protein cavities – binding sites of bioactive compounds, peptides or other proteins

Ο σκοπός της διπλωματικής εργασίας είναι η διερεύνηση και αποτύπωση των ερευνητικών μελετών που αφορούν στον χαρακτηρισμό μιας πρωτεϊνικής κοιλότητας – κέντρου πρόσδεσης βιοδραστικών ενώσεων, πεπτιδίων ή άλλων πρωτεϊνών. Στην παρούσα εργασία χρησιμοποιήθηκε η μέθοδος της βιβλιογραφικής επισκόπησης. Παρουσιάζονται τα κυριότερα ευρήματα προηγούμενων ερευνών που σχετίζονται με τη διαδικασία σχεδιασμού φαρμάκων και τον εντοπισμό φαρμακοφόρων με βάση ένα σύνολο προσδετών. Στη συνέχεια συγκρίνονται διαδικασίες επεξεργασίας και ανάλυσης της πρωτεϊνικής κοιλότητας προγενέστερων ερευνών με τη προσέγγιση που προτάθηκε από τους Παπαθανασίου και Φωτόπουλου το 2015. Αναδεικνύονται βασικά πλεονεκτήματα της προσέγγισης αυτής, όπως η εφαρμογή του αλγορίθμου πολυδιάστατη k-means ομαδοποίηση (multidimensional k-means clustering). Η εύρεση βιβλιογραφίας βασίστηκε σε αναζήτηση επιστημονικών άρθρων σε ξενόγλωσσα επιστημονικά περιοδικά, σε κεφάλαια βιβλίων και σε διάφορα άρθρα σε ηλεκτρονικούς ιστότοπους σχετικά με τον σχεδιασμό φαρμάκων και τις κοιλότητες που απαντώνται στις πρωτεΐνες. Στην παρούσα εργασία παρουσιάζονται εν συντομία εργαλεία που εντοπίστηκαν χρησιμοποιώντας λέξεις κλειδιά όπως για παράδειγμα δυναμική πρωτεϊνικής κοιλότητας, καταλυτικό κέντρο ενός ενζύμου, πρόσδεση, πρωτεϊνική θήκη κλπ. Στη συνέχεια συγκροτήθηκε κατάλογος με τα εργαλεία βιοπληροφορικής ανάλυσης που βρέθηκαν και ακολούθησε εκτενής αναφορά επιλεκτικά σε κάποια από αυτά. Κριτήριο επιλογής αυτών των εργαλείων αποτέλεσε η ημερομηνία δημοσίευσής τους, οι αλγόριθμοι και η μεθοδολογία που χρησιμοποιούν. Τα εργαλεία αυτά κατηγοριοποιήθηκαν με βάση τις λέξεις κλειδιά που χρησιμοποιήθηκαν για την εξόρυξη των δεδομένων από την βιβλιογραφία. Τέλος πραγματοποιήθηκε συγκριτική μελέτη αυτών αναδεικνύοντας τα πλεονεκτήματα και εστιάζοντας στην περαιτέρω αξιοποίησή τους.The aim of this thesis was to report on the current state-of-the-art of the research conducted concerning mapping of protein cavities with a potential function role as binding sites of bioactive compounds, peptides or other proteins. A literature review was performed with emphasis on the relevant tools developed during the last decade. In addition, the main research findings regarding drug design and druggable targets based on binding sites are presented. Processes performed in protein cavity detection and analysis, of previous research articles, are compared with the approach described by Anaxagoras Fotopoulos and Athanasios Papathanasiou (2015). The results showed that a competitive advantage of their approach is the multidimensional k-means algorithm for clustering. For the bibliographic review the scientific knowledgebase has been used, which includes international articles and journals, book chapters, as well as online articles regarding drug design and protein cavity. Search keywords such as protein cavity dynamics, catalytic sites of enzymes, protein pocket etc. were used to identify bioinformatics tools with text mining. A catalogue of the most recently developed tools is presented followed by a brief description of selected tools. The selection criteria imposed for preparing the catalogue and the detailed description included the publication date, as well as the algorithms and the methods they use. The tools were then classified according to the search keywords. The findings of this research are discussed, and the algorithms and methods they use are compared, highlighting the advantages of protein cavity detection

Structural analysis of plasmid-mediated verotoxin gene producing non-o157 escherichia coli using molecular dynamic simulations

Plasmid-mediated Non-O157 VTEC is a pathogenic E. coli serotype, are responsible for many life-threatening diseases such as diarrhoea. The thermostable crystal structure has been widely sought after for industry and therapeutic applications through the structural analysis. The National Center for Biotechnology Information (NCBI) Genbank Database has been sourced to obtain plasmid-mediated Verotoxin genes producing non-O157 MN 696158 (Vtx1-1) and MN688720 (Vtx2) sequences. The tertiary structure of MN696158 (Vtx1-1) and MN688720 (Vtx2) was generated by operating MODELLER. The result exposed multiple templates during modelling processes have improved the local stereochemical quality of the produced models. The structural analysis also disclosed the similarities and differences between the models (Vtx1-1 and Vtx2). Furthermore, the thermal stability profile of plasmid-mediated non-O157 VTEC were studied. The molecular dynamics simulations of plasmid-mediated non-O157 VTEC structures (Vtx1-1 and Vtx2) illustrated the interactions between amino acids. On the other hand, Vtx1-1 and Vtx2 showed noticeable differences in their relative conformational flexibility and stability at elevated temperature. However, it is expected that the information of the thermal stable of plasmid-mediated non-O157 VTEC models can be used for potential vaccine candidate through protein engineering in future

Evaluation of model refinement in CASP13.

Performance in the model refinement category of the 13th round of Critical Assessment of Structure Prediction (CASP13) is assessed, showing that some groups consistently improve most starting models whereas the majority of participants continue to degrade the starting model on average. Using the ranking formula developed for CASP12, it is shown that only 7 of 32 groups perform better than a "naïve predictor" who just submits the starting model. Common features in their approaches include a dependence on physics-based force fields to judge alternative conformations and the use of molecular dynamics to relax models to local minima, usually with some restraints to prevent excessively large movements. In addition to the traditional CASP metrics that focus largely on the quality of the overall fold, alternative metrics are evaluated, including comparisons of the main-chain and side-chain torsion angles, and the utility of the models for solving crystal structures by the molecular replacement method. It is proposed that the introduction of these metrics, as well as consideration of the accuracy of coordinate error estimates, would improve the discrimination between good and very good models.Wellcome Trust Marie Sklowdowska-Curie grain for EU Horizon 202

Mass & secondary structure propensity of amino acids explain their mutability and evolutionary replacements

Why is an amino acid replacement in a protein accepted during evolution? The answer given by bioinformatics relies on the frequency of change of each amino acid by another one and the propensity of each to remain unchanged. We propose that these replacement rules are recoverable from the secondary structural trends of amino acids. A distance measure between high-resolution Ramachandran distributions reveals that structurally similar residues coincide with those found in substitution matrices such as BLOSUM: Asn Asp, Phe Tyr, Lys Arg, Gln Glu, Ile Val, Met → Leu; with Ala, Cys, His, Gly, Ser, Pro, and Thr, as structurally idiosyncratic residues. We also found a high average correlation (\overline{R} R = 0.85) between thirty amino acid mutability scales and the mutational inertia (I X ), which measures the energetic cost weighted by the number of observations at the most probable amino acid conformation. These results indicate that amino acid substitutions follow two optimally-efficient principles: (a) amino acids interchangeability privileges their secondary structural similarity, and (b) the amino acid mutability depends directly on its biosynthetic energy cost, and inversely with its frequency. These two principles are the underlying rules governing the observed amino acid substitutions. © 2017 The Author(s)

Common synaptic phenotypes arising from diverse mutations in the human NMDA receptor subunit GluN2A

Dominant mutations in the human gene GRIN2A, encoding NMDA receptor (NMDAR) subunit GluN2A, make a significant and growing contribution to the catalogue of published single-gene epilepsies. Understanding the disease mechanism in these epilepsy patients is complicated by the surprising diversity of effects that the mutations have on NMDARs. Here we have examined the cell-autonomous effect of five GluN2A mutations, 3 loss-of-function and 2 gain-of-function, on evoked NMDAR-mediated synaptic currents (NMDA-EPSCs) in CA1 pyramidal neurons in cultured hippocampal slices. Despite the mutants differing in their functional incorporation at synapses, prolonged NMDA-EPSC current decays (with only marginal changes in charge transfer) were a common effect for both gain- and loss-of-function mutants. Modelling NMDA-EPSCs with mutant properties in a CA1 neuron revealed that the effect of GRIN2A mutations can lead to abnormal temporal integration and spine calcium dynamics during trains of concerted synaptic activity. Investigations beyond establishing the molecular defects of GluN2A mutants are much needed to understand their impact on synaptic transmission

A domain based protein structural modelling platform applied in the analysis of alternative splicing

Functional families (FunFams) are a sub-classification of CATH protein domain superfamilies that cluster relatives likely to have very similar structures and functions. The functional purity of FunFams has been demonstrated by comparing against experimentally determined Enzyme Commission annotations and by checking whether known functional sites coincide with highly conserved residues in the multiple sequence alignments of FunFams. We hypothesised that clustering relatives into FunFams may help in protein structure modelling. In the first work chapter, we demonstrate the structural coherence of domains in FunFams. We then explore the usage of FunFams in protein monomer modelling. The FunFam based protocol produced higher percentages of good models compared to an HHsearch (the state-of-the-art HMM based sequence search tool) based protocol for both close and remote homologs. We developed a modelling pipeline that, utilises the FunFam protocol, and is able to model up to 70% of domain sequences from human and fly genomes. In the second work chapter, we explore the usage of FunFams in protein complex modelling. Our analysis demonstrated that domain-domain interfaces in FunFams tend to be conserved. The FunFam based complex modelling protocol produced significantly more good quality models when compared to a BLAST based protocol and slightly better than a HHsearch based protocol. In the final work chapter, we employ the FunFam based structural modelling tool to understand the implications of alternative splicing. We focused on isoforms derived from mutually exclusively exons (MXEs) for which there is more enriched in proteomics data. MXEs which could be mapped to structure show a significant tendency to be exposed to the solvent, are likely to exhibit a significant change in their physiochemical property and to lie close to a known/predicted functional sites. Our results suggest that MXE events may have a number of important roles in cells generally

Local Protein Structure Refinement via Molecular Dynamics Simulations with locPREFMD

A method for the local refinement of protein structures that targets improvements in local stereochemistry while preserving the overall fold is presented. The method uses force field-based minimization and sampling via molecular dynamics simulations with a modified force field to bring bonds, angles, and torsion angles into an acceptable range for high-resolution protein structures. The method is implemented in the locPREFMD web server and was tested on computational models submitted to CASP11. Using MolProbity scores as the main assessment criterion, the locPREFMD method significantly improves the stereochemical quality of given input models close to the quality expected for experimental structures while maintaining the Cα coordinates of the initial model