Stereochemistry of Polypeptide Conformation in Coarse Grained Analysis

The conformations available to polypeptides are determined by the interatomic forces acting on the peptide units, whereby backbone torsion angles are restricted as described by the Ramachandran plot. Although typical proteins are composed predominantly from {\alpha}-helices and {\beta}-sheets, they nevertheless adopt diverse tertiary structure, each folded as dictated by its unique amino-acid sequence. Despite such uniqueness, however, the functioning of many proteins involves changes between quite different conformations. The study of large-scale conformational changes, particularly in large systems, is facilitated by a coarse-grained representation such as provided by virtually bonded C{\alpha} atoms. We have developed a virtual atom molecular mechanics (VAMM) force field to describe conformational dynamics in proteins and a VAMM-based algorithm for computing conformational transition pathways. Here we describe the stereochemical analysis of proteins in this coarse-grained representation, comparing the relevant plots in coarse-grained conformational space to the corresponding Ramachandran plots, having contoured each at levels determined statistically from residues in a large database. The distributions shown for an all-{\alpha} protein, two all-{\beta} proteins and one {\alpha}+{\beta} protein serve to relate the coarse-grained distributions to the familiar Ramachandran plot.Comment: 12 pages, 3 figures, Postprint of book chapter submitted to the Biomolecular Forms and Functions, M. Bansal and N. Srinivasan, Eds. copyright (2013) [copyright World Scientific Publishing Company

An effective all-atom potential for proteins

We describe and test an implicit solvent all-atom potential for simulations of protein folding and aggregation. The potential is developed through studies of structural and thermodynamic properties of 17 peptides with diverse secondary structure. Results obtained using the final form of the potential are presented for all these peptides. The same model, with unchanged parameters, is furthermore applied to a heterodimeric coiled-coil system, a mixed alpha/beta protein and a three-helix-bundle protein, with very good results. The computational efficiency of the potential makes it possible to investigate the free-energy landscape of these 49--67-residue systems with high statistical accuracy, using only modest computational resources by today's standards

Monte Carlo Update for Chain Molecules: Biased Gaussian Steps in Torsional Space

We develop a new elementary move for simulations of polymer chains in torsion angle space. The method is flexible and easy to implement. Tentative updates are drawn from a (conformation-dependent) Gaussian distribution that favors approximately local deformations of the chain. The degree of bias is controlled by a parameter b. The method is tested on a reduced model protein with 54 amino acids and the Ramachandran torsion angles as its only degrees of freedom, for different b. Without excessive fine tuning, we find that the effective step size can be increased by a factor of three compared to the unbiased b=0 case. The method may be useful for kinetic studies, too.Comment: 14 pages, 4 figure

Peptide redesign for inhibition of the complement system: Targeting age-related macular degeneration.

PurposeTo redesign a complement-inhibiting peptide with the potential to become a therapeutic for dry and wet age-related macular degeneration (AMD).MethodsWe present a new potent peptide (Peptide 2) of the compstatin family. The peptide is developed by rational design, based on a mechanistic binding hypothesis, and structural and physicochemical properties derived from molecular dynamics (MD) simulation. The inhibitory activity, efficacy, and solubility of Peptide 2 are evaluated using a hemolytic assay, a human RPE cell-based assay, and ultraviolet (UV) absorption properties, respectively, and compared to the respective properties of its parent peptide (Peptide 1).ResultsThe sequence of Peptide 2 contains an arginine-serine N-terminal extension (a characteristic of parent Peptide 1) and a novel 8-polyethylene glycol (PEG) block C-terminal extension. Peptide 2 has significantly improved aqueous solubility compared to Peptide 1 and comparable complement inhibitory activity. In addition, Peptide 2 is more efficacious in inhibiting complement activation in a cell-based model that mimics the pathobiology of dry AMD.ConclusionsWe have designed a new peptide analog of compstatin that combines N-terminal polar amino acid extensions and C-terminal PEGylation extensions. This peptide demonstrates significantly improved aqueous solubility and complement inhibitory efficacy, compared to the parent peptide. The new peptide overcomes the aggregation limitation for clinical translation of previous compstatin analogs and is a candidate to become a therapeutic for the treatment of AMD

Investigation of sequence features of hinge-bending regions in proteins with domain movements using kernel logistic regression

Background: Hinge-bending movements in proteins comprising two or more domains form a large class of functional movements. Hinge-bending regions demarcate protein domains and collectively control the domain movement. Consequently, the ability to recognise sequence features of hinge-bending regions and to be able to predict them from sequence alone would benefit various areas of protein research. For example, an understanding of how the sequence features of these regions relate to dynamic properties in multi-domain proteins would aid in the rational design of linkers in therapeutic fusion proteins. Results: The DynDom database of protein domain movements comprises sequences annotated to indicate whether the amino acid residue is located within a hinge-bending region or within an intradomain region. Using statistical methods and Kernel Logistic Regression (KLR) models, this data was used to determine sequence features that favour or disfavour hinge-bending regions. This is a difficult classification problem as the number of negative cases (intradomain residues) is much larger than the number of positive cases (hinge residues). The statistical methods and the KLR models both show that cysteine has the lowest propensity for hinge-bending regions and proline has the highest, even though it is the most rigid amino acid. As hinge-bending regions have been previously shown to occur frequently at the terminal regions of the secondary structures, the propensity for proline at these regions is likely due to its tendency to break secondary structures. The KLR models also indicate that isoleucine may act as a domain-capping residue. We have found that a quadratic KLR model outperforms a linear KLR model and that improvement in performance occurs up to very long window lengths (eighty residues) indicating long-range correlations. Conclusion: In contrast to the only other approach that focused solely on interdomain hinge-bending regions, the method provides a modest and statistically significant improvement over a random classifier. An explanation of the KLR results is that in the prediction of hinge-bending regions a long-range correlation is at play between a small number amino acids that either favour or disfavour hinge-bending regions. The resulting sequence-based prediction tool, HingeSeek, is available to run through a webserver at hingeseek.cmp.uea.ac.uk

Omics assisted N-terminal proteoform and protein expression profiling on methionine aminopeptidase 1 (MetAP1) deletion

Excision of the N-terminal initiator methionine (iMet) residue from nascent peptide chains is an essential and omnipresent protein modification carried out by methionine aminopeptidases (MetAPs) that accounts for a major source of N-terminal proteoform diversity. Although MetAP2 is known to be implicated in processes such as angiogenesis and proliferation in mammals, the physiological role of MetAP1 is much less clear. In this report we studied the omics-wide effects of human MetAP1 deletion and general MetAP inhibition. The levels of iMet retention are inversely correlated with cellular proliferation rates. Further, despite the increased MetAP2 expression on MetAP1 deletion, MetAP2 was unable to restore processing of Met-Ser-, Met-Pro-, and Met-Ala- starting N termini as inferred from the iMet retention profiles observed, indicating a higher activity of MetAP1 over these N termini. Proteome and transcriptome expression profiling point to differential expression of proteins implicated in lipid metabolism, cytoskeleton organization, cell proliferation and protein synthesis upon perturbation of MetAP activity

The Genetic Code as a Periodic Table: Algebraic Aspects

The systematics of indices of physico-chemical properties of codons and amino acids across the genetic code are examined. Using a simple numerical labelling scheme for nucleic acid bases, data can be fitted as low-order polynomials of the 6 coordinates in the 64-dimensional codon weight space. The work confirms and extends recent studies by Siemion of amino acid conformational parameters. The connections between the present work, and recent studies of the genetic code structure using dynamical symmetry algebras, are pointed out.Comment: 26 pages Latex, 10 figures (4 ps, 6 Tex). Refereed version, small changes to discussion (conclusion unaltered). Minor alterations to format of figures and tables. To appear in BioSystem

Actions of Agonists, Fipronil and Ivermectin on the Predominant In Vivo Splice and Edit Variant (RDLbd, I/V) of the Drosophila GABA Receptor Expressed in Xenopus laevis Oocytes

Ionotropic GABA receptors are the targets for several classes of insecticides. One of the most widely-studied insect GABA receptors is RDL (resistance to dieldrin), originally isolated from Drosophila melanogaster. RDL undergoes alternative splicing and RNA editing, which influence the potency of GABA. Most work has focussed on minority isoforms. Here, we report the first characterisation of the predominant native splice variant and RNA edit, combining functional characterisation with molecular modelling of the agonist-binding region. The relative order of agonist potency is GABA> muscimol> TACA> β-alanine. The I/V edit does not alter the potency of GABA compared to RDLbd. Docking calculations suggest that these agonists bind and activate RDLbdI/V through a similar binding mode. TACA and β-alanine are predicted to bind with lower affinity than GABA, potentially explaining their lower potency, whereas the lower potency of muscimol and isoguvacine cannot be explained structurally from the docking calculations. The A301S (resistance to dieldrin) mutation reduced the potency of antagonists picrotoxin, fipronil and pyrafluprole but the I/V edit had no measurable effect. Ivermectin suppressed responses to GABA of RDLbdI/V, RDLbd and RDLbdI/VA301S. The dieldrin resistant variant also showed reduced sensitivity to Ivermectin. This study of a highly abundant insect GABA receptor isoform will help the design of new insecticides

Specialized dynamical properties of promiscuous residues revealed by simulated conformational ensembles

The ability to interact with different partners is one of the most important features in proteins. Proteins that bind a large number of partners (hubs) have been often associated with intrinsic disorder. However, many examples exist of hubs with an ordered structure, and evidence of a general mechanism promoting promiscuity in ordered proteins is still elusive. An intriguing hypothesis is that promiscuous binding sites have specific dynamical properties, distinct from the rest of the interface and pre-existing in the protein isolated state. Here, we present the first comprehensive study of the intrinsic dynamics of promiscuous residues in a large protein data set. Different computational methods, from coarse-grained elastic models to geometry-based sampling methods and to full-atom Molecular Dynamics simulations, were used to generate conformational ensembles for the isolated proteins. The flexibility and dynamic correlations of interface residues with a different degree of binding promiscuity were calculated and compared considering side chain and backbone motions, the latter both on a local and on a global scale. The study revealed that (a) promiscuous residues tend to be more flexible than nonpromiscuous ones, (b) this additional flexibility has a higher degree of organization, and (c) evolutionary conservation and binding promiscuity have opposite effects on intrinsic dynamics. Findings on simulated ensembles were also validated on ensembles of experimental structures extracted from the Protein Data Bank (PDB). Additionally, the low occurrence of single nucleotide polymorphisms observed for promiscuous residues indicated a tendency to preserve binding diversity at these positions. A case study on two ubiquitin-like proteins exemplifies how binding promiscuity in evolutionary related proteins can be modulated by the fine-tuning of the interface dynamics. The interplay between promiscuity and flexibility highlighted here can inspire new directions in protein-protein interaction prediction and design methods. © 2013 American Chemical Society

Deletion of TSPO resulted in change of metabolomic profile in retinal pigment epithelial cells

Age-related macular degeneration is the main cause of vision loss in the aged population worldwide. Drusen, extracellular lesions formed underneath the retinal pigment epithelial (RPE) cells, are a clinical feature of AMD and associated with AMD progression. RPE cells support photoreceptor function by providing nutrition, phagocytosing outer segments and removing metabolic waste. Dysfunction and death of RPE cells are early features of AMD. The translocator protein, TSPO, plays an important role in RPE cholesterol efflux and loss of TSPO results in increased intracellular lipid accumulation and reactive oxygen species (ROS) production. This study aimed to investigate the impact of TSPO knockout on RPE cellular metabolism by identifying the metabolic differences between wildtype and knockout RPE cells, with or without treatment with oxidized low density lipoprotein (oxLDL). Using liquid chromatography mass spectrometry (LC/MS), we differentiated several metabolic pathways among wildtype and knockout cells. Lipids amongst other intracellular metabolites were the most influenced by loss of TSPO and/or oxLDL treatment. Glucose, amino acid and nucleotide metabolism was also affected. TSPO deletion led to up-regulation of fatty acids and glycerophospholipids, which in turn possibly affected the cell membrane fluidity and stability. Higher levels of glutathione disulphide (GSSG) were found in TSPO knockout RPE cells, suggesting TSPO regulates mitochondrial-mediated oxidative stress. These data provide biochemical insights into TSPO-associated function in RPE cells and may shed light on disease mechanisms in AMD