1
|
Lee J, Ryu B, Kim T, Kim KK. Cryo-EM structure of a 16.5-kDa small heat-shock protein from Methanocaldococcus jannaschii. Int J Biol Macromol 2024; 258:128763. [PMID: 38103675 DOI: 10.1016/j.ijbiomac.2023.128763] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/23/2023] [Revised: 12/07/2023] [Accepted: 12/11/2023] [Indexed: 12/19/2023]
Abstract
The small heat-shock protein (sHSP) from the archaea Methanocaldococcus jannaschii, MjsHSP16.5, functions as a broad substrate ATP-independent holding chaperone protecting misfolded proteins from aggregation under stress conditions. This protein is the first sHSP characterized by X-ray crystallography, thereby contributing significantly to our understanding of sHSPs. However, despite numerous studies assessing its functions and structures, the precise arrangement of the N-terminal domains (NTDs) within this sHSP cage remains elusive. Here we present the cryo-electron microscopy (cryo-EM) structure of MjsHSP16.5 at 2.49-Å resolution. The subunits of MjsHSP16.5 in the cryo-EM structure exhibit lesser compaction compared to their counterparts in the crystal structure. This structural feature holds particular significance in relation to the biophysical properties of MjsHSP16.5, suggesting a close resemblance to this sHSP native state. Additionally, our cryo-EM structure unveils the density of residues 24-33 within the NTD of MjsHSP16.5, a feature that typically remains invisible in the majority of its crystal structures. Notably, these residues show a propensity to adopt a β-strand conformation and engage in antiparallel interactions with strand β1, both intra- and inter-subunit modes. These structural insights are corroborated by structural predictions, disulfide bond cross-linking studies of Cys-substitution mutants, and protein disaggregation assays. A comprehensive understanding of the structural features of MjsHSP16.5 expectedly holds the potential to inspire a wide range of interdisciplinary applications, owing to the renowned versatility of this sHSP as a nanoscale protein platform.
Collapse
Affiliation(s)
- Joohyun Lee
- Department of Precision Medicine, Graduate School of Basic Medical Science (GSBMS), Institute for Antimicrobial Resistance Research and Therapeutics, Sungkyunkwan University School of Medicine, Suwon 16419, Republic of Korea
| | - Bumhan Ryu
- Research Solution Center, Institute for Basic Science (IBS), Daejeon 34126, Republic of Korea
| | - Truc Kim
- Department of Precision Medicine, Graduate School of Basic Medical Science (GSBMS), Institute for Antimicrobial Resistance Research and Therapeutics, Sungkyunkwan University School of Medicine, Suwon 16419, Republic of Korea.
| | - Kyeong Kyu Kim
- Department of Precision Medicine, Graduate School of Basic Medical Science (GSBMS), Institute for Antimicrobial Resistance Research and Therapeutics, Sungkyunkwan University School of Medicine, Suwon 16419, Republic of Korea.
| |
Collapse
|
2
|
Math BA, Waibl F, Lamp LM, Fernández‐Quintero ML, Liedl KR. Cross-linking disulfide bonds govern solution structures of diabodies. Proteins 2023; 91:1316-1328. [PMID: 37376973 PMCID: PMC10952579 DOI: 10.1002/prot.26509] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/12/2022] [Accepted: 04/19/2023] [Indexed: 06/29/2023]
Abstract
In the last years, antibodies have emerged as a promising new class of therapeutics, due to their combination of high specificity with long serum half-life and low risk of side-effects. Diabodies are a popular novel antibody format, consisting of two Fv domains connected with short linkers. Like IgG antibodies, they simultaneously bind two target proteins. However, they offer altered properties, given their smaller size and higher rigidity. In this study, we conducted the-to our knowledge-first molecular dynamics (MD) simulations of diabodies and find a surprisingly high conformational flexibility in the relative orientation of the two Fv domains. We observe rigidifying effects through the introduction of disulfide bonds in the Fv -Fv interface and characterize the effect of different disulfide bond locations on the conformation. Additionally, we compare VH -VL orientations and paratope dynamics between diabodies and an antigen binding fragment (Fab) of the same sequence. We find mostly consistent structures and dynamics, indicating similar antigen binding properties. The most significant differences can be found within the CDR-H2 loop dynamics. Of all CDR loops, the CDR-H2 is located closest to the artificial Fv -Fv interface. All examined diabodies show similar VH -VL orientations, Fv -Fv packing and CDR loop conformations. However, the variant with a P14C-K64C disulfide bond differs most from the Fab in our measures, including the CDR-H3 loop conformational ensemble. This suggests altered antigen binding properties and underlines the need for careful validation of the disulfide bond locations in diabodies.
Collapse
Affiliation(s)
- Barbara A. Math
- Institute of General, Inorganic and Theoretical Chemistry, and Center for Molecular Biosciences Innsbruck (CMBI)University of InnsbruckInnsbruckAustria
| | - Franz Waibl
- Institute of General, Inorganic and Theoretical Chemistry, and Center for Molecular Biosciences Innsbruck (CMBI)University of InnsbruckInnsbruckAustria
| | - Leonida M. Lamp
- Institute of General, Inorganic and Theoretical Chemistry, and Center for Molecular Biosciences Innsbruck (CMBI)University of InnsbruckInnsbruckAustria
| | - Monica L. Fernández‐Quintero
- Institute of General, Inorganic and Theoretical Chemistry, and Center for Molecular Biosciences Innsbruck (CMBI)University of InnsbruckInnsbruckAustria
| | - Klaus R. Liedl
- Institute of General, Inorganic and Theoretical Chemistry, and Center for Molecular Biosciences Innsbruck (CMBI)University of InnsbruckInnsbruckAustria
| |
Collapse
|
3
|
Abanades B, Wong WK, Boyles F, Georges G, Bujotzek A, Deane CM. ImmuneBuilder: Deep-Learning models for predicting the structures of immune proteins. Commun Biol 2023; 6:575. [PMID: 37248282 DOI: 10.1038/s42003-023-04927-7] [Citation(s) in RCA: 50] [Impact Index Per Article: 50.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/15/2022] [Accepted: 05/11/2023] [Indexed: 05/31/2023] Open
Abstract
Immune receptor proteins play a key role in the immune system and have shown great promise as biotherapeutics. The structure of these proteins is critical for understanding their antigen binding properties. Here, we present ImmuneBuilder, a set of deep learning models trained to accurately predict the structure of antibodies (ABodyBuilder2), nanobodies (NanoBodyBuilder2) and T-Cell receptors (TCRBuilder2). We show that ImmuneBuilder generates structures with state of the art accuracy while being far faster than AlphaFold2. For example, on a benchmark of 34 recently solved antibodies, ABodyBuilder2 predicts CDR-H3 loops with an RMSD of 2.81Å, a 0.09Å improvement over AlphaFold-Multimer, while being over a hundred times faster. Similar results are also achieved for nanobodies, (NanoBodyBuilder2 predicts CDR-H3 loops with an average RMSD of 2.89Å, a 0.55Å improvement over AlphaFold2) and TCRs. By predicting an ensemble of structures, ImmuneBuilder also gives an error estimate for every residue in its final prediction. ImmuneBuilder is made freely available, both to download ( https://github.com/oxpig/ImmuneBuilder ) and to use via our webserver ( http://opig.stats.ox.ac.uk/webapps/newsabdab/sabpred ). We also make available structural models for ~150 thousand non-redundant paired antibody sequences ( https://doi.org/10.5281/zenodo.7258553 ).
Collapse
Affiliation(s)
| | - Wing Ki Wong
- Large Molecule Research, Roche Pharma Research and Early Development, Roche Innovation Center Munich, Penzberg, Germany
| | - Fergus Boyles
- Department of Statistics, University of Oxford, Oxford, UK
| | - Guy Georges
- Large Molecule Research, Roche Pharma Research and Early Development, Roche Innovation Center Munich, Penzberg, Germany
| | - Alexander Bujotzek
- Large Molecule Research, Roche Pharma Research and Early Development, Roche Innovation Center Munich, Penzberg, Germany
| | | |
Collapse
|
4
|
Chenchiliyan M, Kübel J, Ooi SA, Salvadori G, Mennucci B, Westenhoff S, Maj M. Ground-state heterogeneity and vibrational energy redistribution in bacterial phytochrome observed with femtosecond 2D IR spectroscopy. J Chem Phys 2023; 158:085103. [PMID: 36859103 DOI: 10.1063/5.0135268] [Citation(s) in RCA: 4] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/25/2023] Open
Abstract
Phytochromes belong to a group of photoreceptor proteins containing a covalently bound biliverdin chromophore that inter-converts between two isomeric forms upon photoexcitation. The existence and stability of the photocycle products are largely determined by the protein sequence and the presence of conserved hydrogen-bonding interactions in the vicinity of the chromophore. The vibrational signatures of biliverdin, however, are often weak and obscured under more intense protein bands, limiting spectroscopic studies of its non-transient signals. In this study, we apply isotope-labeling techniques to isolate the vibrational bands from the protein-bound chromophore of the bacterial phytochrome from Deinococcus radiodurans. We elucidate the structure and ultrafast dynamics of the chromophore with 2D infra-red (IR) spectroscopy and molecular dynamics simulations. The carbonyl stretch vibrations of the pyrrole rings show the heterogeneous distribution of hydrogen-bonding structures, which exhibit distinct ultrafast relaxation dynamics. Moreover, we resolve a previously undetected 1678 cm-1 band that is strongly coupled to the A- and D-ring of biliverdin and demonstrate the presence of complex vibrational redistribution pathways between the biliverdin modes with relaxation-assisted measurements of 2D IR cross peaks. In summary, we expect 2D IR spectroscopy to be useful in explaining how point mutations in the protein sequence affect the hydrogen-bonding structure around the chromophore and consequently its ability to photoisomerize to the light-activated states.
Collapse
Affiliation(s)
- Manoop Chenchiliyan
- Department of Chemistry and Molecular Biology, University of Gothenburg, Box 462, 40530 Gothenburg, Sweden
| | - Joachim Kübel
- Department of Chemistry and Molecular Biology, University of Gothenburg, Box 462, 40530 Gothenburg, Sweden
| | - Saik Ann Ooi
- Department of Chemistry and Molecular Biology, University of Gothenburg, Box 462, 40530 Gothenburg, Sweden
| | - Giacomo Salvadori
- Department of Chemistry and Industrial Chemistry, University of Pisa, Via G. Moruzzi 13, 56126 Pisa, Italy
| | - Benedetta Mennucci
- Department of Chemistry and Industrial Chemistry, University of Pisa, Via G. Moruzzi 13, 56126 Pisa, Italy
| | - Sebastian Westenhoff
- Department of Chemistry and Molecular Biology, University of Gothenburg, Box 462, 40530 Gothenburg, Sweden
| | - Michał Maj
- Department of Chemistry-Ångström Laboratory, Uppsala University, SE-75120 Uppsala, Sweden
| |
Collapse
|
5
|
Řezáč J, Stewart JJP. How well do semiempirical QM methods describe the structure of proteins? J Chem Phys 2023; 158:044118. [PMID: 36725526 DOI: 10.1063/5.0135091] [Citation(s) in RCA: 6] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/12/2023] Open
Abstract
Semiempirical quantum-mechanical (QM) computational methods are an increasingly popular tool for the study of biomolecular systems. They were, however, developed and tested mostly on small model molecules. In this work, we explore one topic fundamental to these applications: the ability of the methods to describe the structure of proteins. In a set of 19 proteins for which a crystal structure with very high resolution is available, we analyze the properties of the protein geometries optimized using several semiempirical QM methods including PM6-D3H4, PM7, and GFN2-xTB. Some of the methods provide a very good description of the general structural features of the protein, yielding results better than or comparable to the AMBER ff03 force field. However, PM7 and PM6-D3H4 optimizations introduce artificial close contacts in the structure, which is partially remediated by reparameterization.
Collapse
Affiliation(s)
- J Řezáč
- Institute of Organic Chemistry and Biochemistry, Czech Academy of Sciences, 16000 Prague, Czech Republic
| | - J J P Stewart
- Stewart Computational Chemistry, 15210 Paddington Circle, Colorado Springs, Colorado 80921, USA
| |
Collapse
|
6
|
Shen Y, Parks JM, Smith JC. HLA Class I Supertype Classification Based on Structural Similarity. JOURNAL OF IMMUNOLOGY (BALTIMORE, MD. : 1950) 2023; 210:103-114. [PMID: 36453976 DOI: 10.4049/jimmunol.2200685] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 09/13/2022] [Accepted: 10/31/2022] [Indexed: 12/24/2022]
Abstract
HLA class I proteins, a critical component in adaptive immunity, bind and present intracellular Ags to CD8+ T cells. The extreme polymorphism of HLA genes and associated peptide binding specificities leads to challenges in various endeavors, including neoantigen vaccine development, disease association studies, and HLA typing. Supertype classification, defined by clustering functionally similar HLA alleles, has proven helpful in reducing the complexity of distinguishing alleles. However, determining supertypes via experiments is impractical, and current in silico classification methods exhibit limitations in stability and functional relevance. In this study, by incorporating three-dimensional structures we present a method for classifying HLA class I molecules with improved breadth, accuracy, stability, and flexibility. Critical for these advances is our finding that structural similarity highly correlates with peptide binding specificity. The new classification should be broadly useful in peptide-based vaccine development and HLA-disease association studies.
Collapse
Affiliation(s)
- Yue Shen
- UT-ORNL Graduate School of Genome Science and Technology, University of Tennessee, Knoxville, TN
| | - Jerry M Parks
- Biosciences Division, Oak Ridge National Laboratory, Oak Ridge, TN; and
| | - Jeremy C Smith
- UT-ORNL Graduate School of Genome Science and Technology, University of Tennessee, Knoxville, TN.,Biosciences Division, Oak Ridge National Laboratory, Oak Ridge, TN; and.,Department of Biochemistry & Cellular and Molecular Biology, University of Tennessee, Knoxville, TN
| |
Collapse
|
7
|
Liu N, Mikhailovskii O, Skrynnikov NR, Xue Y. Simulating diffraction photographs based on molecular dynamics trajectories of a protein crystal: a new option to examine structure-solving strategies in protein crystallography. IUCRJ 2023; 10:16-26. [PMID: 36598499 PMCID: PMC9812212 DOI: 10.1107/s2052252522011198] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 06/06/2022] [Accepted: 11/21/2022] [Indexed: 06/17/2023]
Abstract
A molecular dynamics (MD)-based pipeline has been designed and implemented to emulate the entire process of collecting diffraction photographs and calculating crystallographic structures of proteins from them. Using a structure of lysozyme solved in-house, a supercell comprising 125 (5 × 5 × 5) crystal unit cells containing a total of 1000 protein molecules and explicit interstitial solvent was constructed. For this system, two 300 ns MD trajectories at 298 and 250 K were recorded. A series of snapshots from these trajectories were then used to simulate a fully realistic set of diffraction photographs, which were further fed into the standard pipeline for structure determination. The resulting structures show very good agreement with the underlying MD model not only in terms of coordinates but also in terms of B factors; they are also consistent with the original experimental structure. The developed methodology should find a range of applications, such as optimizing refinement protocols to solve crystal structures and extracting dynamics information from diffraction data or diffuse scattering.
Collapse
Affiliation(s)
- Ning Liu
- School of Life Sciences, Tsinghua University, Beijing 100084, People’s Republic of China
| | - Oleg Mikhailovskii
- Laboratory of Biomolecular NMR, St Petersburg State University, St Petersburg, Russian Federation
- Department of Chemistry, Purdue University, West Lafayette, IN 47907, USA
| | - Nikolai R. Skrynnikov
- Laboratory of Biomolecular NMR, St Petersburg State University, St Petersburg, Russian Federation
- Department of Chemistry, Purdue University, West Lafayette, IN 47907, USA
| | - Yi Xue
- School of Life Sciences, Tsinghua University, Beijing 100084, People’s Republic of China
- Beijing Advanced Innovation Center for Structural Biology, Tsinghua University, Beijing 100084, People’s Republic of China
- Tsinghua University–Peking University Joint Center for Life Sciences, Tsinghua University, Beijing 100084, People’s Republic of China
| |
Collapse
|
8
|
Richardson BC, Shek R, Van Voorhis WC, French JB. Structure of Klebsiella pneumoniae adenosine monophosphate nucleosidase. PLoS One 2022; 17:e0275023. [PMID: 36264993 PMCID: PMC9584410 DOI: 10.1371/journal.pone.0275023] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/13/2022] [Accepted: 09/09/2022] [Indexed: 11/05/2022] Open
Abstract
Klebsiella pneumoniae is a bacterial pathogen that is increasingly responsible for hospital-acquired pneumonia and sepsis. Progressive development of antibiotic resistance has led to higher mortality rates and creates a need for novel treatments. Because of the essential role that nucleotides play in many bacterial processes, enzymes involved in purine and pyrimidine metabolism and transport are ideal targets for the development of novel antibiotics. Herein we describe the structure of K. pneumoniae adenosine monophosphate nucleosidase (KpAmn), a purine salvage enzyme unique to bacteria, as determined by cryoelectron microscopy. The data detail a well conserved fold with a hexameric overall structure and clear density for the putative active site residues. Comparison to the crystal structures of homologous prokaryotic proteins confirms the presence of many of the conserved structural features of this protein yet reveals differences in distal loops in the absence of crystal contacts. This first cryo-EM structure of an Amn enzyme provides a basis for future structure-guided drug development and extends the accuracy of structural characterization of this family of proteins beyond this clinically relevant organism.
Collapse
Affiliation(s)
- Brian C. Richardson
- The Hormel Institute, University of Minnesota, Austin, Minnesota, United States of America
| | - Roger Shek
- Division of Allergy and Infectious Diseases, Department of Medicine, Center for Emerging and Re-emerging Infectious Diseases, University of Washington School of Medicine, Seattle, WA, United States of America
| | - Wesley C. Van Voorhis
- Division of Allergy and Infectious Diseases, Department of Medicine, Center for Emerging and Re-emerging Infectious Diseases, University of Washington School of Medicine, Seattle, WA, United States of America
| | - Jarrod B. French
- The Hormel Institute, University of Minnesota, Austin, Minnesota, United States of America
| |
Collapse
|
9
|
Laurents DV. AlphaFold 2 and NMR Spectroscopy: Partners to Understand Protein Structure, Dynamics and Function. Front Mol Biosci 2022; 9:906437. [PMID: 35655760 PMCID: PMC9152297 DOI: 10.3389/fmolb.2022.906437] [Citation(s) in RCA: 25] [Impact Index Per Article: 12.5] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/28/2022] [Accepted: 04/25/2022] [Indexed: 11/29/2022] Open
Abstract
The artificial intelligence program AlphaFold 2 is revolutionizing the field of protein structure determination as it accurately predicts the 3D structure of two thirds of the human proteome. Its predictions can be used directly as structural models or indirectly as aids for experimental structure determination using X-ray crystallography, CryoEM or NMR spectroscopy. Nevertheless, AlphaFold 2 can neither afford insight into how proteins fold, nor can it determine protein stability or dynamics. Rare folds or minor alternative conformations are also not predicted by AlphaFold 2 and the program does not forecast the impact of post translational modifications, mutations or ligand binding. The remaining third of human proteome which is poorly predicted largely corresponds to intrinsically disordered regions of proteins. Key to regulation and signaling networks, these disordered regions often form biomolecular condensates or amyloids. Fortunately, the limitations of AlphaFold 2 are largely complemented by NMR spectroscopy. This experimental approach provides information on protein folding and dynamics as well as biomolecular condensates and amyloids and their modulation by experimental conditions, small molecules, post translational modifications, mutations, flanking sequence, interactions with other proteins, RNA and virus. Together, NMR spectroscopy and AlphaFold 2 can collaborate to advance our comprehension of proteins.
Collapse
|
10
|
Gaber A, Pavšič M. Modeling and Structure Determination of Homo-Oligomeric Proteins: An Overview of Challenges and Current Approaches. Int J Mol Sci 2021; 22:9081. [PMID: 34445785 PMCID: PMC8396596 DOI: 10.3390/ijms22169081] [Citation(s) in RCA: 9] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/30/2021] [Revised: 08/20/2021] [Accepted: 08/20/2021] [Indexed: 12/12/2022] Open
Abstract
Protein homo-oligomerization is a very common phenomenon, and approximately half of proteins form homo-oligomeric assemblies composed of identical subunits. The vast majority of such assemblies possess internal symmetry which can be either exploited to help or poses challenges during structure determination. Moreover, aspects of symmetry are critical in the modeling of protein homo-oligomers either by docking or by homology-based approaches. Here, we first provide a brief overview of the nature of protein homo-oligomerization. Next, we describe how the symmetry of homo-oligomers is addressed by crystallographic and non-crystallographic symmetry operations, and how biologically relevant intermolecular interactions can be deciphered from the ordered array of molecules within protein crystals. Additionally, we describe the most important aspects of protein homo-oligomerization in structure determination by NMR. Finally, we give an overview of approaches aimed at modeling homo-oligomers using computational methods that specifically address their internal symmetry and allow the incorporation of other experimental data as spatial restraints to achieve higher model reliability.
Collapse
|
11
|
Kochert M, Nocek BP, Habeeb Mohammad TS, Gild E, Lovato K, Heath TK, Holz RC, Olsen KW, Becker DP. Atomic-Resolution 1.3 Å Crystal Structure, Inhibition by Sulfate, and Molecular Dynamics of the Bacterial Enzyme DapE. Biochemistry 2021; 60:908-917. [PMID: 33721990 DOI: 10.1021/acs.biochem.0c00926] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
We report the atomic-resolution (1.3 Å) X-ray crystal structure of an open conformation of the dapE-encoded N-succinyl-l,l-diaminopimelic acid desuccinylase (DapE, EC 3.5.1.18) from Neisseria meningitidis. This structure [Protein Data Bank (PDB) entry 5UEJ] contains two bound sulfate ions in the active site that mimic the binding of the terminal carboxylates of the N-succinyl-l,l-diaminopimelic acid (l,l-SDAP) substrate. We demonstrated inhibition of DapE by sulfate (IC50 = 13.8 ± 2.8 mM). Comparison with other DapE structures in the PDB demonstrates the flexibility of the interdomain connections of this protein. This high-resolution structure was then utilized as the starting point for targeted molecular dynamics experiments revealing the conformational change from the open form to the closed form that occurs when DapE binds l,l-SDAP and cleaves the amide bond. These simulations demonstrated closure from the open to the closed conformation, the change in RMS throughout the closure, and the independence in the movement of the two DapE subunits. This conformational change occurred in two phases with the catalytic domains moving toward the dimerization domains first, followed by a rotation of catalytic domains relative to the dimerization domains. Although there were no targeting forces, the substrate moved closer to the active site and bound more tightly during the closure event.
Collapse
Affiliation(s)
- Matthew Kochert
- Department of Chemistry and Biochemistry, Loyola University Chicago, 1032 West Sheridan Road, Chicago, Illinois 60660, United States
| | - Boguslaw P Nocek
- The Center for Structural Genomics of Infectious Diseases, Computation Institute, The University of Chicago, Chicago, Illinois 60637, United States
| | - Thahani S Habeeb Mohammad
- Department of Chemistry and Biochemistry, Loyola University Chicago, 1032 West Sheridan Road, Chicago, Illinois 60660, United States
| | - Elliot Gild
- Department of Chemistry and Biochemistry, Loyola University Chicago, 1032 West Sheridan Road, Chicago, Illinois 60660, United States
| | - Kaitlyn Lovato
- Department of Chemistry and Biochemistry, Loyola University Chicago, 1032 West Sheridan Road, Chicago, Illinois 60660, United States
| | - Tahirah K Heath
- Department of Chemistry and Biochemistry, Loyola University Chicago, 1032 West Sheridan Road, Chicago, Illinois 60660, United States
| | - Richard C Holz
- Department of Chemistry, Colorado School of Mines, 1500 Illinois Street, Golden, Colorado 80401, United States
| | - Kenneth W Olsen
- Department of Chemistry and Biochemistry, Loyola University Chicago, 1032 West Sheridan Road, Chicago, Illinois 60660, United States
| | - Daniel P Becker
- Department of Chemistry and Biochemistry, Loyola University Chicago, 1032 West Sheridan Road, Chicago, Illinois 60660, United States
| |
Collapse
|
12
|
Reinknecht C, Riga A, Rivera J, Snyder DA. Patterns in Protein Flexibility: A Comparison of NMR "Ensembles", MD Trajectories, and Crystallographic B-Factors. Molecules 2021; 26:molecules26051484. [PMID: 33803249 PMCID: PMC7967184 DOI: 10.3390/molecules26051484] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/31/2020] [Revised: 02/18/2021] [Accepted: 02/28/2021] [Indexed: 11/16/2022] Open
Abstract
Proteins are molecular machines requiring flexibility to function. Crystallographic B-factors and Molecular Dynamics (MD) simulations both provide insights into protein flexibility on an atomic scale. Nuclear Magnetic Resonance (NMR) lacks a universally accepted analog of the B-factor. However, a lack of convergence in atomic coordinates in an NMR-based structure calculation also suggests atomic mobility. This paper describes a pattern in the coordinate uncertainties of backbone heavy atoms in NMR-derived structural “ensembles” first noted in the development of FindCore2 (previously called Expanded FindCore: DA Snyder, J Grullon, YJ Huang, R Tejero, GT Montelione, Proteins: Structure, Function, and Bioinformatics 82 (S2), 219–230) and demonstrates that this pattern exists in coordinate variances across MD trajectories but not in crystallographic B-factors. This either suggests that MD trajectories and NMR “ensembles” capture motional behavior of peptide bond units not captured by B-factors or indicates a deficiency common to force fields used in both NMR and MD calculations.
Collapse
|
13
|
Ferreira C, Pinto MF, Macedo-Ribeiro S, Pereira PJB, Rocha FA, Martins PM. Protein crystals as a key for deciphering macromolecular crowding effects on biological reactions. Phys Chem Chem Phys 2020; 22:16143-16149. [PMID: 32638771 DOI: 10.1039/d0cp02469d] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/30/2022]
Abstract
When placed in the same environment, biochemically unrelated macromolecules influence each other's biological function through macromolecular crowding (MC) effects. This has been illustrated in vitro by the effects of inert polymers on protein stability, protein structure, enzyme kinetics and protein aggregation kinetics. While a unified way to quantitatively characterize MC is still lacking, we show that the crystal solubility of lysozyme can be used to predict the influence of crowding agents on the catalytic efficiency of this enzyme. In order to capture general enthalpic effects, as well as hard entropic effects that are specific of large molecules, we tested sucrose and its cross-linked polymer Ficoll-70 as additives. Despite the different conditions of pH and ionic strength adopted, both the crystallization and the enzymatic assays point to an entropic contribution of approximately -1 kcal mol-1 caused by MC. Our results demonstrate that the thermodynamic activity of proteins is markedly increased by the reduction of accessible volume caused by the presence of macromolecular cosolutes. Unlike what is observed in protein folding studies, this MC effect cannot be reproduced using equivalent concentrations of monomeric crowding units. Applicable to any crystallizable protein, the thermodynamic interpretation of MC based on crystal solubility is expected to help in elucidating the full extent and importance of hard-type interactions in the crowded environment of the cell.
Collapse
Affiliation(s)
- Cecília Ferreira
- LEPABE - Laboratory for Process Engineering, Environment, Biotechnology and Energy, Faculty of Engineering, University of Porto, Rua Dr. Roberto Frias, 4200-465 Porto, Portugal
| | - Maria Filipa Pinto
- LEPABE - Laboratory for Process Engineering, Environment, Biotechnology and Energy, Faculty of Engineering, University of Porto, Rua Dr. Roberto Frias, 4200-465 Porto, Portugal and ICBAS - Instituto de Ciências Biomédicas Abel Salazar, Universidade do Porto, Porto, Portugal and IBMC - Instituto de Biologia Molecular e Celular, Universidade do Porto, Porto, Portugal. and Instituto de Investigação e Inovação em Saúde, Universidade do Porto, Porto, Portugal
| | - Sandra Macedo-Ribeiro
- IBMC - Instituto de Biologia Molecular e Celular, Universidade do Porto, Porto, Portugal. and Instituto de Investigação e Inovação em Saúde, Universidade do Porto, Porto, Portugal
| | - Pedro José Barbosa Pereira
- IBMC - Instituto de Biologia Molecular e Celular, Universidade do Porto, Porto, Portugal. and Instituto de Investigação e Inovação em Saúde, Universidade do Porto, Porto, Portugal
| | - Fernando Alberto Rocha
- LEPABE - Laboratory for Process Engineering, Environment, Biotechnology and Energy, Faculty of Engineering, University of Porto, Rua Dr. Roberto Frias, 4200-465 Porto, Portugal
| | - Pedro Miguel Martins
- ICBAS - Instituto de Ciências Biomédicas Abel Salazar, Universidade do Porto, Porto, Portugal and IBMC - Instituto de Biologia Molecular e Celular, Universidade do Porto, Porto, Portugal. and Instituto de Investigação e Inovação em Saúde, Universidade do Porto, Porto, Portugal
| |
Collapse
|
14
|
Schmitz S, Seibert J, Ostermeir K, Hansen A, Göller AH, Grimme S. Quantum Chemical Calculation of Molecular and Periodic Peptide and Protein Structures. J Phys Chem B 2020; 124:3636-3646. [PMID: 32275425 DOI: 10.1021/acs.jpcb.0c00549] [Citation(s) in RCA: 25] [Impact Index Per Article: 6.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/24/2022]
Abstract
Special-purpose classical force fields (FFs) provide good accuracy at very low computational cost, but their application is limited to systems for which potential energy functions are available. This excludes most metal-containing proteins or those containing cofactors. In contrast, the GFN2-xTB semiempirical quantum chemical method is parametrized for almost the entire periodic table. The accuracy of GFN2-xTB is assessed for protein structures with respect to experimental X-ray data. Furthermore, the results are compared with those of two special-purpose FFs, HF-3c, PM6-D3H4X, and PM7. The test sets include proteins without any prosthetic groups as well as metalloproteins. Crystal packing effects are examined for a set of smaller proteins to validate the molecular approach. For the proteins without prosthetic groups, the special purpose FF OPLS-2005 yields the smallest overall RMSD to the X-ray data but GFN2-xTB provides similarly good structures with even better bond-length distributions. For the metalloproteins with up to 5000 atoms, a good overall structural agreement is obtained with GFN2-xTB. The full geometry optimizations of protein structures with on average 1000 atoms in wall-times below 1 day establishes the GFN2-xTB method as a versatile tool for the computational treatment of various biomolecules with a good accuracy/computational cost ratio.
Collapse
Affiliation(s)
- Sarah Schmitz
- Mulliken Center for Theoretical Chemistry, Institut für Physikalische und Theoretische Chemie der Universität Bonn, Beringstraße 4, D-53115 Bonn, Germany
| | - Jakob Seibert
- Mulliken Center for Theoretical Chemistry, Institut für Physikalische und Theoretische Chemie der Universität Bonn, Beringstraße 4, D-53115 Bonn, Germany
| | - Katja Ostermeir
- Mulliken Center for Theoretical Chemistry, Institut für Physikalische und Theoretische Chemie der Universität Bonn, Beringstraße 4, D-53115 Bonn, Germany
| | - Andreas Hansen
- Mulliken Center for Theoretical Chemistry, Institut für Physikalische und Theoretische Chemie der Universität Bonn, Beringstraße 4, D-53115 Bonn, Germany
| | - Andreas H Göller
- Computational Molecular Design, Pharmaceuticals, R&D, Bayer AG, D-42096 Wuppertal, Germany
| | - Stefan Grimme
- Mulliken Center for Theoretical Chemistry, Institut für Physikalische und Theoretische Chemie der Universität Bonn, Beringstraße 4, D-53115 Bonn, Germany
| |
Collapse
|
15
|
Farhadi Z, Farhadi T, Hashemian SM. Virtual screening for potential inhibitors of β(1,3)-D-glucan synthase as drug candidates against fungal cell wall. J Drug Assess 2020; 9:52-59. [PMID: 32284908 PMCID: PMC7144292 DOI: 10.1080/21556660.2020.1734010] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/06/2019] [Accepted: 02/07/2020] [Indexed: 01/17/2023] Open
Abstract
Background To enhance the outcome in patients with invasive candidiasis, initiation of an efficient antifungal treatment in a suitable dosage is necessary. Echinocandins (e.g. caspofungin) inhibit the enzyme β(1,3)-D-glucan synthase of the fungal cell wall. Compared to azoles and other antifungal agents, echinocandins have lower adverse effects and toxicity in humans. Echinocandins are available in injectable (intravenous) form. Methods In this study, to identify the novel oral drug-like compounds that affect the fungal cell wall, downloaded oral drug-like compounds from the ZINC database were processed with a virtual screening procedure. The docking free energies were calculated and compared with the known inhibitor caspofungin. Four molecules were selected as the most potent ligands and subjected to hydrogen bonds analysis. Results Considering the hydrogen bond analysis, two compounds (ZINC71336662 and ZINC40910772) were predicted to better interact with the active site of β(1,3)-D-glucan synthase compared with caspofungin. Conclusion The introduced compound in this study may be valuable to analyze experimentally as a novel oral drug candidate targeting fungal cell walls.
Collapse
Affiliation(s)
- Zinat Farhadi
- Chronic Respiratory Diseases Research Center (CRDRC), National Research Institute of Tuberculosis and Lung Diseases (NRITLD), Shahid Beheshti University of Medical Sciences, Tehran, Iran.,Behavioral Disease Counseling Center, Marvdasht Health Center, Shiraz University of Medical Sciences, Shiraz, Iran.,Department of Microbiology, Shiraz Branch, Islamic Azad University, Shiraz, Iran
| | - Tayebeh Farhadi
- Chronic Respiratory Diseases Research Center (CRDRC), National Research Institute of Tuberculosis and Lung Diseases (NRITLD), Shahid Beheshti University of Medical Sciences, Tehran, Iran
| | - Seyed MohammadReza Hashemian
- Chronic Respiratory Diseases Research Center (CRDRC), National Research Institute of Tuberculosis and Lung Diseases (NRITLD), Shahid Beheshti University of Medical Sciences, Tehran, Iran.,Critical Care Department, Farhikhtegan Hospital, Tehran Medical Branch, Islamic Azad University, Tehran, Iran
| |
Collapse
|
16
|
Bala S, Shinya S, Srivastava A, Ishikawa M, Shimada A, Kobayashi N, Kojima C, Tama F, Miyashita O, Kohda D. Crystal contact-free conformation of an intrinsically flexible loop in protein crystal: Tim21 as the case study. Biochim Biophys Acta Gen Subj 2019; 1864:129418. [PMID: 31449839 DOI: 10.1016/j.bbagen.2019.129418] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/18/2019] [Revised: 08/02/2019] [Accepted: 08/22/2019] [Indexed: 02/06/2023]
Abstract
BACKGROUND In protein crystals, flexible loops are frequently deformed by crystal contacts, whereas in solution, the large motions result in the poor convergence of such flexible loops in NMR structure determinations. We need an experimental technique to characterize the structural and dynamic properties of intrinsically flexible loops of protein molecules. METHODS We designed an intended crystal contact-free space (CCFS) in protein crystals, and arranged the flexible loop of interest in the CCFS. The yeast Tim 21 protein was chosen as the model protein, because one of the loops (loop 2) is distorted by crystal contacts in the conventional crystal. RESULTS Yeast Tim21 was fused to the MBP protein by a rigid α-helical linker. The space created between the two proteins was used as the CCFS. The linker length provides adjustable freedom to arrange loop 2 in the CCFS. We re-determined the NMR structure of yeast Tim21, and conducted MD simulations for comparison. Multidimensional scaling was used to visualize the conformational similarity of loop 2. We found that the crystal contact-free conformation of loop 2 is located close to the center of the ensembles of the loop 2 conformations in the NMR and MD structures. CONCLUSIONS Loop 2 of yeast Tim21 in the CCFS adopts a representative, dominant conformation in solution. GENERAL SIGNIFICANCE No single powerful technique is available for the characterization of flexible structures in protein molecules. NMR analyses and MD simulations provide useful, but incomplete information. CCFS crystallography offers a third route to this goal.
Collapse
Affiliation(s)
- Siqin Bala
- Division of Structural Biology, Medical Institute of Bioregulation, Kyushu University, Maidashi 3-1-1, Higashi-ku, Fukuoka 812-8582, Japan
| | - Shoko Shinya
- Institute for Protein Research, Osaka University, Yamadaoka 3-2, Suita, Osaka 565-0871, Japan
| | - Arpita Srivastava
- Department of Physics, Graduate School of Science, Nagoya University, Furo-cho, Chikusa-ku, Nagoya, Aichi 464-8602, Japan
| | - Marie Ishikawa
- Division of Structural Biology, Medical Institute of Bioregulation, Kyushu University, Maidashi 3-1-1, Higashi-ku, Fukuoka 812-8582, Japan
| | - Atsushi Shimada
- Division of Structural Biology, Medical Institute of Bioregulation, Kyushu University, Maidashi 3-1-1, Higashi-ku, Fukuoka 812-8582, Japan
| | - Naohiro Kobayashi
- Institute for Protein Research, Osaka University, Yamadaoka 3-2, Suita, Osaka 565-0871, Japan
| | - Chojiro Kojima
- Institute for Protein Research, Osaka University, Yamadaoka 3-2, Suita, Osaka 565-0871, Japan; Graduate School of Engineering Science, Yokohama National University, 79-5 Tokiwadai, Hodogaya-ku, Yokohama 240-8501, Japan
| | - Florence Tama
- Department of Physics, Graduate School of Science, Nagoya University, Furo-cho, Chikusa-ku, Nagoya, Aichi 464-8602, Japan; Institute of Transformative Bio-Molecules (WPI-ITbM), Nagoya University, Furo-cho, Chikusa-ku, Nagoya, Aichi 464-8602, Japan; Center for Computational Science, RIKEN, 6-7-1 Minatojima-minami-machi, Chuo-ku, Kobe, Hyogo 650-0047, Japan
| | - Osamu Miyashita
- Center for Computational Science, RIKEN, 6-7-1 Minatojima-minami-machi, Chuo-ku, Kobe, Hyogo 650-0047, Japan
| | - Daisuke Kohda
- Division of Structural Biology, Medical Institute of Bioregulation, Kyushu University, Maidashi 3-1-1, Higashi-ku, Fukuoka 812-8582, Japan.
| |
Collapse
|
17
|
Holt ME, Salay LE, O’Brien E, Barton JK, Chazin WJ. Functional and structural similarity of human DNA primase [4Fe4S] cluster domain constructs. PLoS One 2018; 13:e0209345. [PMID: 30562384 PMCID: PMC6298731 DOI: 10.1371/journal.pone.0209345] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/25/2018] [Accepted: 11/14/2018] [Indexed: 01/15/2023] Open
Abstract
The regulatory subunit of human DNA primase has a C-terminal domain (p58C) that contains a [4Fe4S] cluster and binds DNA. Previous electrochemical analysis of a p58C construct revealed that its affinity for DNA is sensitive to the redox state of the [4Fe4S] cluster. Concerns about the validity of this conclusion have been raised, based in part on differences in X-ray crystal structures of the p58C272-464 construct used for that study and that of a N-terminally shifted p58C266-456 construct and consequently, an assumption that p58C272-464 has abnormal physical and functional properties. To address this controversy, a new p58C266-464 construct containing all residues was crystallized under the conditions previously used for crystallizing p58C272-464, and the solution structures of both constructs were assessed using circular dichroism and NMR spectroscopy. In the new crystal structure, p58C266-464 exhibits the same elements of secondary structure near the DNA binding site as observed in the crystal structure of p58C272-464. Moreover, in solution, circular dichroism and 15N,1H-heteronuclear single quantum coherence (HSQC) NMR spectra show there are no significant differences in the distribution of secondary structures or in the tertiary structure or the two constructs. To validate that the two constructs have the same functional properties, binding of a primed DNA template was measured using a fluorescence-based DNA binding assay, and the affinities for this substrate were the same (3.4 ± 0.5 μM and 2.7 ± 0.3 μM, respectively). The electrochemical properties of p58C266-464 were also measured and this p58C construct was able to engage in redox switching on DNA with the same efficiency as p58C272-464. Together, these results show that although p58C can be stabilized in different conformations in the crystalline state, in solution there is effectively no difference in the structure and functional properties of p58C constructs of different lengths.
Collapse
Affiliation(s)
- Marilyn E. Holt
- Departments of Biochemistry and Chemistry, and Center for Structural Biology, Vanderbilt University, Nashville, Tennessee, United States of America
| | - Lauren E. Salay
- Departments of Biochemistry and Chemistry, and Center for Structural Biology, Vanderbilt University, Nashville, Tennessee, United States of America
| | - Elizabeth O’Brien
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California, United States of America
| | - Jacqueline K. Barton
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California, United States of America
| | - Walter J. Chazin
- Departments of Biochemistry and Chemistry, and Center for Structural Biology, Vanderbilt University, Nashville, Tennessee, United States of America
- * E-mail:
| |
Collapse
|
18
|
Kumar A, Karthikeyan S. Crystal structure of the MSMEG_4306 gene product from Mycobacterium smegmatis. Acta Crystallogr F Struct Biol Commun 2018; 74:166-173. [PMID: 29497021 PMCID: PMC5947703 DOI: 10.1107/s2053230x18002236] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/04/2017] [Accepted: 02/06/2018] [Indexed: 12/31/2022] Open
Abstract
The MSMEG_4306 gene from Mycobacterium smegmatis encodes a protein of unknown function with 242 amino-acid residues that contains a conserved zinc-ribbon domain at its C-terminus. Here, the crystal structure of MSMEG_4306 determined by the single-wavelength anomalous dispersion method using just one zinc ion co-purified with the protein is reported. The crystal structure of MSMEG_4306 shows a coiled-coil helix domain in the N-terminal region and a zinc-ribbon domain in the C-terminal region. A structural similarity search against the Protein Data Bank using MSMEG_4306 as a query revealed two similar structures, namely CT398 from Chlamydia trachomatis and HP0958 from Helicobacter pylori, although they share only ∼15% sequence identity with MSMEG_4306. Based on comparative analysis, it is predicted that MSMEG_4306 may be involved in secretion systems, possibly by interacting with multiple proteins or nucleic acids.
Collapse
Affiliation(s)
- Adarsh Kumar
- CSIR – Institute of Microbial Technology, Council of Scientific and Industrial Research (CSIR), Sector 39A, Chandigarh 160 036, India
| | - Subramanian Karthikeyan
- CSIR – Institute of Microbial Technology, Council of Scientific and Industrial Research (CSIR), Sector 39A, Chandigarh 160 036, India
| |
Collapse
|
19
|
Condic-Jurkic K, Subramanian N, Mark AE, O’Mara ML. The reliability of molecular dynamics simulations of the multidrug transporter P-glycoprotein in a membrane environment. PLoS One 2018; 13:e0191882. [PMID: 29370310 PMCID: PMC5785007 DOI: 10.1371/journal.pone.0191882] [Citation(s) in RCA: 31] [Impact Index Per Article: 5.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/15/2017] [Accepted: 01/13/2018] [Indexed: 11/19/2022] Open
Abstract
Despite decades of research, the mechanism of action of the ABC multidrug transporter P-glycoprotein (P-gp) remains elusive. Due to experimental limitations, many researchers have turned to molecular dynamics simulation studies in order to investigate different aspects of P-gp function. However, such studies are challenging and caution is required when interpreting the results. P-gp is highly flexible and the time scale on which it can be simulated is limited. There is also uncertainty regarding the accuracy of the various crystal structures available, let alone the structure of the protein in a physiologically relevant environment. In this study, three alternative structural models of mouse P-gp (3G5U, 4KSB, 4M1M), all resolved to 3.8 Å, were used to initiate sets of simulations of P-gp in a membrane environment in order to determine: a) the sensitivity of the results to differences in the starting configuration; and b) the extent to which converged results could be expected on the times scales commonly simulated for this system. The simulations suggest that the arrangement of the nucleotide binding domains (NBDs) observed in the crystal structures is not stable in a membrane environment. In all simulations, the NBDs rapidly associated (within 10 ns) and changes within the transmembrane helices were observed. The secondary structure within the transmembrane domain was best preserved in the 4M1M model under the simulation conditions used. However, the extent to which replicate simulations diverged on a 100 to 200 ns timescale meant that it was not possible to draw definitive conclusions as to which structure overall was most stable, or to obtain converged and reliable results for any of the properties examined. The work brings into question the reliability of conclusions made in regard to the nature of specific interactions inferred from previous simulation studies on this system involving similar sampling times. It also highlights the need to demonstrate the statistical significance of any results obtained in simulations of large flexible proteins, especially where the initial structure is uncertain.
Collapse
Affiliation(s)
- Karmen Condic-Jurkic
- Research School of Chemistry, The Australian National University, Canberra, ACT 2601, Australia
- School of Chemistry and Molecular Biosciences, University of Queensland, Brisbane, ACT 4072, Australia
| | - Nandhitha Subramanian
- Research School of Chemistry, The Australian National University, Canberra, ACT 2601, Australia
- School of Chemistry and Molecular Biosciences, University of Queensland, Brisbane, ACT 4072, Australia
| | - Alan E. Mark
- School of Chemistry and Molecular Biosciences, University of Queensland, Brisbane, ACT 4072, Australia
| | - Megan L. O’Mara
- Research School of Chemistry, The Australian National University, Canberra, ACT 2601, Australia
| |
Collapse
|
20
|
Monzon AM, Zea DJ, Marino-Buslje C, Parisi G. Homology modeling in a dynamical world. Protein Sci 2017; 26:2195-2206. [PMID: 28815769 DOI: 10.1002/pro.3274] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/06/2017] [Revised: 08/09/2017] [Accepted: 08/09/2017] [Indexed: 12/31/2022]
Abstract
A key concept in template-based modeling (TBM) is the high correlation between sequence and structural divergence, with the practical consequence that homologous proteins that are similar at the sequence level will also be similar at the structural level. However, conformational diversity of the native state will reduce the correlation between structural and sequence divergence, because structural variation can appear without sequence diversity. In this work, we explore the impact that conformational diversity has on the relationship between structural and sequence divergence. We find that the extent of conformational diversity can be as high as the maximum structural divergence among families. Also, as expected, conformational diversity impairs the well-established correlation between sequence and structural divergence, which is nosier than previously suggested. However, we found that this noise can be resolved using a priori information coming from the structure-function relationship. We show that protein families with low conformational diversity show a well-correlated relationship between sequence and structural divergence, which is severely reduced in proteins with larger conformational diversity. This lack of correlation could impair TBM results in highly dynamical proteins. Finally, we also find that the presence of order/disorder can provide useful beforehand information for better TBM performance.
Collapse
Affiliation(s)
- Alexander Miguel Monzon
- Departamento de Ciencia y Tecnología, Universidad Nacional de Quilmes, CONICET, B1876BXD, Bernal, Argentina
| | - Diego Javier Zea
- Structural Bioinformatics Unit, Fundación Instituto Leloir, CONICET, C1405BWE Ciudad Autónoma de Buenos Aires, Argentina
| | - Cristina Marino-Buslje
- Structural Bioinformatics Unit, Fundación Instituto Leloir, CONICET, C1405BWE Ciudad Autónoma de Buenos Aires, Argentina
| | - Gustavo Parisi
- Departamento de Ciencia y Tecnología, Universidad Nacional de Quilmes, CONICET, B1876BXD, Bernal, Argentina
| |
Collapse
|
21
|
Proteome scale identification, classification and structural analysis of iron-binding proteins in bread wheat. J Inorg Biochem 2017; 170:63-74. [DOI: 10.1016/j.jinorgbio.2017.02.012] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/02/2016] [Revised: 01/23/2017] [Accepted: 02/10/2017] [Indexed: 12/26/2022]
|
22
|
Virtual Screening for Potential Inhibitors of CTX-M-15 Protein of Klebsiella pneumoniae. Interdiscip Sci 2017; 10:694-703. [PMID: 28374117 DOI: 10.1007/s12539-017-0222-y] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/06/2016] [Revised: 02/28/2017] [Accepted: 03/09/2017] [Indexed: 12/31/2022]
Abstract
The Gram-negative bacterium Klebsiella pneumoniae, responsible for a wide variety of nosocomial infections in immuno-deficient patients, involves the respiratory, urinary and gastrointestinal tract infections and septicemia. Extended spectrum β-lactamases (ESBL) belong to β-lactamases capable of conferring antibiotic resistance in Gram-negative bacteria. CTX-M-15, a prevalent ESBL reported from Enterobacteriaceae including K. pneumoniae, was selected as a potent anti-bacterial target. To identify the novel drug-like compounds, structure-based screening procedure was employed against downloaded drug-like compounds from ZINC database. An acronym for "ZINC" is not commercial. The docking free energy values were investigated and compared to the known inhibitor Avibactam. Six best novel drug-like compounds were selected and their hydrogen bindings with the receptor were determined. Based on the binding efficiency mode, three among these six identified most potential inhibitors, ZINC21811621, ZINC93091917 and ZINC19488569, were predicted as potential competitive inhibitors against CTX-M-15 compared to Avibactam. These three inhibitors may provide a framework for the experimental studies to develop anti-Klebsiella novel drug candidates targeting CTX-M-15.
Collapse
|
23
|
Quistgaard EM, Martinez Molledo M, Löw C. Structure determination of a major facilitator peptide transporter: Inward facing PepTSt from Streptococcus thermophilus crystallized in space group P3121. PLoS One 2017; 12:e0173126. [PMID: 28264013 PMCID: PMC5338821 DOI: 10.1371/journal.pone.0173126] [Citation(s) in RCA: 24] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/15/2016] [Accepted: 02/15/2017] [Indexed: 12/03/2022] Open
Abstract
Major facilitator superfamily (MFS) peptide transporters (typically referred to as PepT, POT or PTR transporters) mediate the uptake of di- and tripeptides, and so play an important dietary role in many organisms. In recent years, a better understanding of the molecular basis for this process has emerged, which is in large part due to a steep increase in structural information. Yet, the conformational transitions underlying the transport mechanism are still not fully understood, and additional data is therefore needed. Here we report in detail the detergent screening, crystallization, experimental MIRAS phasing, and refinement of the peptide transporter PepTSt from Streptococcus thermophilus. The space group is P3121, and the protein is crystallized in a monomeric inward facing form. The binding site is likely to be somewhat occluded, as the lobe encompassing transmembrane helices 10 and 11 is markedly bent towards the central pore of the protein, but the extent of this potential occlusion could not be determined due to disorder at the apex of the lobe. Based on structural comparisons with the seven previously determined P212121 and C2221 structures of inward facing PepTSt, the structural flexibility as well as the conformational changes mediating transition between the inward open and inward facing occluded states are discussed. In conclusion, this report improves our understanding of the structure and conformational cycle of PepTSt, and can furthermore serve as a case study, which may aid in supporting future structure determinations of additional MFS transporters or other integral membrane proteins.
Collapse
Affiliation(s)
- Esben M. Quistgaard
- Centre for Structural Systems Biology (CSSB), DESY and European Molecular Biology Laboratory Hamburg, Hamburg, Germany
- Department of Medical Biochemistry and Biophysics, Karolinska Institutet, Stockholm, Sweden
| | - Maria Martinez Molledo
- Centre for Structural Systems Biology (CSSB), DESY and European Molecular Biology Laboratory Hamburg, Hamburg, Germany
| | - Christian Löw
- Centre for Structural Systems Biology (CSSB), DESY and European Molecular Biology Laboratory Hamburg, Hamburg, Germany
- Department of Medical Biochemistry and Biophysics, Karolinska Institutet, Stockholm, Sweden
- * E-mail:
| |
Collapse
|
24
|
Structural variability of E. coli thioredoxin captured in the crystal structures of single-point mutants. Sci Rep 2017; 7:42343. [PMID: 28181556 PMCID: PMC5299410 DOI: 10.1038/srep42343] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/19/2016] [Accepted: 01/09/2017] [Indexed: 11/11/2022] Open
Abstract
Thioredoxin is a ubiquitous small protein that catalyzes redox reactions of protein thiols. Additionally, thioredoxin from E. coli (EcTRX) is a widely-used model for structure-function studies. In a previous paper, we characterized several single-point mutants of the C-terminal helix (CTH) that alter global stability of EcTRX. However, spectroscopic signatures and enzymatic activity for some of these mutants were found essentially unaffected. A comprehensive structural characterization at the atomic level of these near-invariant mutants can provide detailed information about structural variability of EcTRX. We address this point through the determination of the crystal structures of four point-mutants, whose mutations occurs within or near the CTH, namely L94A, E101G, N106A and L107A. These structures are mostly unaffected compared with the wild-type variant. Notably, the E101G mutant presents a large region with two alternative traces for the backbone of the same chain. It represents a significant shift in backbone positions. Enzymatic activity measurements and conformational dynamics studies monitored by NMR and molecular dynamic simulations show that E101G mutation results in a small effect in the structural features of the protein. We hypothesize that these alternative conformations represent samples of the native-state ensemble of EcTRX, specifically the magnitude and location of conformational heterogeneity.
Collapse
|
25
|
Ahlstrom LS, Vorontsov II, Shi J, Miyashita O. Effect of the Crystal Environment on Side-Chain Conformational Dynamics in Cyanovirin-N Investigated through Crystal and Solution Molecular Dynamics Simulations. PLoS One 2017; 12:e0170337. [PMID: 28107510 PMCID: PMC5249168 DOI: 10.1371/journal.pone.0170337] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/14/2016] [Accepted: 01/03/2017] [Indexed: 11/18/2022] Open
Abstract
Side chains in protein crystal structures are essential for understanding biochemical processes such as catalysis and molecular recognition. However, crystal packing could influence side-chain conformation and dynamics, thus complicating functional interpretations of available experimental structures. Here we investigate the effect of crystal packing on side-chain conformational dynamics with crystal and solution molecular dynamics simulations using Cyanovirin-N as a model system. Side-chain ensembles for solvent-exposed residues obtained from simulation largely reflect the conformations observed in the X-ray structure. This agreement is most striking for crystal-contacting residues during crystal simulation. Given the high level of correspondence between our simulations and the X-ray data, we compare side-chain ensembles in solution and crystal simulations. We observe large decreases in conformational entropy in the crystal for several long, polar and contacting residues on the protein surface. Such cases agree well with the average loss in conformational entropy per residue upon protein folding and are accompanied by a change in side-chain conformation. This finding supports the application of surface engineering to facilitate crystallization. Our simulation-based approach demonstrated here with Cyanovirin-N establishes a framework for quantitatively comparing side-chain ensembles in solution and in the crystal across a larger set of proteins to elucidate the effect of the crystal environment on protein conformations.
Collapse
Affiliation(s)
- Logan S. Ahlstrom
- Department of Molecular, Cellular, and Developmental Biology, University of Michigan, Ann Arbor, Michigan, United States of America
| | - Ivan I. Vorontsov
- Department of Chemistry and Biochemistry, University of Arizona, Tucson, Arizona, United States of America
| | - Jun Shi
- Department of Chemistry and Biochemistry, University of Arizona, Tucson, Arizona, United States of America
| | - Osamu Miyashita
- RIKEN Advanced Institute for Computational Science, Chuo-ku, Kobe, Hyogo, Japan
- * E-mail:
| |
Collapse
|
26
|
Cunha ES, Sfriso P, Rojas AL, Roversi P, Hospital A, Orozco M, Abrescia NG. Mechanism of Structural Tuning of the Hepatitis C Virus Human Cellular Receptor CD81 Large Extracellular Loop. Structure 2017; 25:53-65. [DOI: 10.1016/j.str.2016.11.003] [Citation(s) in RCA: 17] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/26/2016] [Revised: 09/21/2016] [Accepted: 10/27/2016] [Indexed: 12/12/2022]
|
27
|
He W, Liang Z, Teng M, Niu L. LibME-automatic extraction of 3D ligand-binding motifs for mechanistic analysis of protein-ligand recognition. FEBS Open Bio 2016; 6:1331-1340. [PMID: 28255540 PMCID: PMC5324770 DOI: 10.1002/2211-5463.12150] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/07/2016] [Revised: 10/26/2016] [Accepted: 10/27/2016] [Indexed: 11/23/2022] Open
Abstract
Identifying conserved binding motifs is an efficient way to study protein–ligand recognition. Most 3D binding motifs only contain information from the protein side, and so motifs that combine information from both protein and ligand sides are desired. Here, we propose an algorithm called LibME (Ligand‐binding Motif Extractor), which automatically extracts 3D binding motifs composed of the target ligand and surrounding conserved residues. We show that the motifs extracted by LibME for ATP and its analogs are highly similar to well‐known motifs reported by previous studies. The superiority of our method to handle flexible ligands was also demonstrated using isocitric acid as an example. Finally, we show that these motifs, together with their visual exhibition, permit better investigating and understanding of protein–ligand recognition process.
Collapse
Affiliation(s)
- Wei He
- Hefei National Laboratory for Physical Sciences at Microscale and School of Life Sciences University of Science and Technology of China Anhui China
| | - Zhi Liang
- Hefei National Laboratory for Physical Sciences at Microscale and School of Life Sciences University of Science and Technology of China Anhui China
| | - MaiKun Teng
- Hefei National Laboratory for Physical Sciences at Microscale and School of Life Sciences University of Science and Technology of China Anhui China
| | - LiWen Niu
- Hefei National Laboratory for Physical Sciences at Microscale and School of Life Sciences University of Science and Technology of China Anhui China
| |
Collapse
|
28
|
North RA, Watson AJA, Pearce FG, Muscroft-Taylor AC, Friemann R, Fairbanks AJ, Dobson RCJ. Structure and inhibition of N-acetylneuraminate lyase from methicillin-resistant Staphylococcus aureus. FEBS Lett 2016; 590:4414-4428. [PMID: 27943302 DOI: 10.1002/1873-3468.12462] [Citation(s) in RCA: 16] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/19/2016] [Revised: 10/04/2016] [Accepted: 10/09/2016] [Indexed: 01/07/2023]
Abstract
N-Acetylneuraminate lyase is the first committed enzyme in the degradation of sialic acid by bacterial pathogens. In this study, we analyzed the kinetic parameters of N-acetylneuraminate lyase from methicillin-resistant Staphylococcus aureus (MRSA). We determined that the enzyme has a relatively high KM of 3.2 mm, suggesting that flux through the catabolic pathway is likely to be controlled by this enzyme. Our data indicate that sialic acid alditol, a known inhibitor of N-acetylneuraminate lyase enzymes, is a stronger inhibitor of MRSA N-acetylneuraminate lyase than of Clostridium perfringens N-acetylneuraminate lyase. Our analysis of the crystal structure of ligand-free and 2R-sialic acid alditol-bound MRSA N-acetylneuraminate lyase suggests that subtle dynamic differences in solution and/or altered binding interactions within the active site may account for species-specific inhibition.
Collapse
Affiliation(s)
- Rachel A North
- Biomolecular Interaction Centre, School of Biological Sciences, University of Canterbury, Christchurch, New Zealand
| | - Andrew J A Watson
- Department of Chemistry, University of Canterbury, Christchurch, New Zealand
| | - F Grant Pearce
- Biomolecular Interaction Centre, School of Biological Sciences, University of Canterbury, Christchurch, New Zealand
| | - Andrew C Muscroft-Taylor
- Protein Science and Engineering, Callaghan Innovation, University of Canterbury, Christchurch, New Zealand
| | - Rosmarie Friemann
- Department of Chemistry and Molecular Biology, University of Gothenburg, Sweden.,Centre for Antibiotic Resistance Research (CARe), University of Gothenburg, Sweden.,Department of Structural Biology, School of Medicine, Stanford University, CA, USA
| | - Antony J Fairbanks
- Biomolecular Interaction Centre, School of Biological Sciences, University of Canterbury, Christchurch, New Zealand.,Department of Chemistry, University of Canterbury, Christchurch, New Zealand
| | - Renwick C J Dobson
- Biomolecular Interaction Centre, School of Biological Sciences, University of Canterbury, Christchurch, New Zealand.,Department of Biochemistry and Molecular Biology, Bio21 Molecular Science and Biotechnology Institute, University of Melbourne, Parkville, Australia
| |
Collapse
|
29
|
Forli S. Charting a Path to Success in Virtual Screening. Molecules 2015; 20:18732-58. [PMID: 26501243 PMCID: PMC4630810 DOI: 10.3390/molecules201018732] [Citation(s) in RCA: 52] [Impact Index Per Article: 5.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/13/2015] [Revised: 10/07/2015] [Accepted: 10/12/2015] [Indexed: 12/27/2022] Open
Abstract
Docking is commonly applied to drug design efforts, especially high-throughput virtual screenings of small molecules, to identify new compounds that bind to a given target. Despite great advances and successful applications in recent years, a number of issues remain unsolved. Most of the challenges and problems faced when running docking experiments are independent of the specific software used, and can be ascribed to either improper input preparation or to the simplified approaches applied to achieve high-throughput speed. Being aware of approximations and limitations of such methods is essential to prevent errors, deal with misleading results, and increase the success rate of virtual screening campaigns. In this review, best practices and most common issues of docking and virtual screening will be discussed, covering the journey from the design of the virtual experiment to the hit identification.
Collapse
Affiliation(s)
- Stefano Forli
- Molecular Graphics Laboratory, Department of Integrative Structural and Computational Biology, The Scripps Research Institute, 10550 North Torrey Pines Road, La Jolla, CA 92037, USA.
| |
Collapse
|
30
|
Messih MA, Lepore R, Tramontano A. LoopIng: a template-based tool for predicting the structure of protein loops. Bioinformatics 2015; 31:3767-72. [PMID: 26249814 PMCID: PMC4653384 DOI: 10.1093/bioinformatics/btv438] [Citation(s) in RCA: 16] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/07/2015] [Accepted: 07/21/2015] [Indexed: 12/31/2022] Open
Abstract
Motivation: Predicting the structure of protein loops is very challenging, mainly because they are not necessarily subject to strong evolutionary pressure. This implies that, unlike the rest of the protein, standard homology modeling techniques are not very effective in modeling their structure. However, loops are often involved in protein function, hence inferring their structure is important for predicting protein structure as well as function. Results: We describe a method, LoopIng, based on the Random Forest automated learning technique, which, given a target loop, selects a structural template for it from a database of loop candidates. Compared to the most recently available methods, LoopIng is able to achieve similar accuracy for short loops (4–10 residues) and significant enhancements for long loops (11–20 residues). The quality of the predictions is robust to errors that unavoidably affect the stem regions when these are modeled. The method returns a confidence score for the predicted template loops and has the advantage of being very fast (on average: 1 min/loop). Availability and implementation:www.biocomputing.it/looping Contact:anna.tramontano@uniroma1.it Supplementary information:Supplementary data are available at Bioinformatics online.
Collapse
Affiliation(s)
| | - Rosalba Lepore
- Department of Physics, Sapienza University, 00185 Rome, Italy and
| | - Anna Tramontano
- Department of Physics, Sapienza University, 00185 Rome, Italy and Istituto Pasteur-Fondazione Cenci Bolognetti, Viale Regina Elena 291, 00161 Rome, Italy
| |
Collapse
|
31
|
Martin BP, Brandon CJ, Stewart JJP, Braun‐Sand SB. Accuracy issues involved in modeling in vivo protein structures using PM7. Proteins 2015; 83:1427-35. [PMID: 25973843 PMCID: PMC4744657 DOI: 10.1002/prot.24826] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/12/2014] [Revised: 04/20/2015] [Accepted: 05/02/2015] [Indexed: 11/07/2022]
Abstract
Using the semiempirical method PM7, an attempt has been made to quantify the error in prediction of the in vivo structure of proteins relative to X-ray structures. Three important contributory factors are the experimental limitations of X-ray structures, the difference between the crystal and solution environments, and the errors due to PM7. The geometries of 19 proteins from the Protein Data Bank that had small R values, that is, high accuracy structures, were optimized and the resulting drop in heat of formation was calculated. Analysis of the changes showed that about 10% of this decrease in heat of formation was caused by faults in PM7, the balance being attributable to the X-ray structure and the difference between the crystal and solution environments. A previously unknown fault in PM7 was revealed during tests to validate the geometries generated using PM7. Clashscores generated by the Molprobity molecular mechanics structure validation program showed that PM7 was predicting unrealistically close contacts between nonbonding atoms in regions where the local geometry is dominated by very weak noncovalent interactions. The origin of this fault was traced to an underestimation of the core-core repulsion between atoms at distances smaller than the equilibrium distance.
Collapse
Affiliation(s)
- Benjamin P. Martin
- Department of Chemistry and BiochemistryUniversity of ColoradoBoulderColorado80918
| | | | - James J. P. Stewart
- Department of Chemistry and BiochemistryUniversity of ColoradoBoulderColorado80918
- Stewart Computational Chemistry15210 Paddington CircleColorado SpringsCO80921
| | - Sonja B. Braun‐Sand
- Department of Chemistry and BiochemistryUniversity of ColoradoBoulderColorado80918
- Stewart Computational Chemistry15210 Paddington CircleColorado SpringsCO80921
| |
Collapse
|
32
|
Ahmed MH, Koparde VN, Safo MK, Neel Scarsdale J, Kellogg GE. 3d interaction homology: The structurally known rotamers of tyrosine derive from a surprisingly limited set of information-rich hydropathic interaction environments described by maps. Proteins 2015; 83:1118-36. [DOI: 10.1002/prot.24813] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/05/2014] [Revised: 02/24/2015] [Accepted: 04/05/2015] [Indexed: 12/20/2022]
Affiliation(s)
- Mostafa H. Ahmed
- Institute for Structural Biology and Drug Discovery, Virginia Commonwealth University; Richmond VA 23298 USA
- Department of Medicinal Chemistry; Virginia Commonwealth University; Richmond VA 23298 USA
| | - Vishal N. Koparde
- Institute for Structural Biology and Drug Discovery, Virginia Commonwealth University; Richmond VA 23298 USA
- Department of Medicinal Chemistry; Virginia Commonwealth University; Richmond VA 23298 USA
| | - Martin K. Safo
- Institute for Structural Biology and Drug Discovery, Virginia Commonwealth University; Richmond VA 23298 USA
- Department of Medicinal Chemistry; Virginia Commonwealth University; Richmond VA 23298 USA
| | - J. Neel Scarsdale
- Institute for Structural Biology and Drug Discovery, Virginia Commonwealth University; Richmond VA 23298 USA
- Center for the Study of Biological Complexity; Virginia Commonwealth University; Richmond VA 23298 USA
| | - Glen E. Kellogg
- Institute for Structural Biology and Drug Discovery, Virginia Commonwealth University; Richmond VA 23298 USA
- Department of Medicinal Chemistry; Virginia Commonwealth University; Richmond VA 23298 USA
| |
Collapse
|
33
|
García-Fernández R, Perbandt M, Rehders D, Ziegelmüller P, Piganeau N, Hahn U, Betzel C, Chávez MDLÁ, Redecke L. Three-dimensional Structure of a Kunitz-type Inhibitor in Complex with an Elastase-like Enzyme. J Biol Chem 2015; 290:14154-65. [PMID: 25878249 DOI: 10.1074/jbc.m115.647586] [Citation(s) in RCA: 16] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/25/2015] [Indexed: 11/06/2022] Open
Abstract
Elastase-like enzymes are involved in important diseases such as acute pancreatitis, chronic inflammatory lung diseases, and cancer. Structural insights into their interaction with specific inhibitors will contribute to the development of novel anti-elastase compounds that resist rapid oxidation and proteolysis. Proteinaceous Kunitz-type inhibitors homologous to the bovine pancreatic trypsin inhibitor (BPTI) provide a suitable scaffold, but the structural aspects of their interaction with elastase-like enzymes have not been elucidated. Here, we increased the selectivity of ShPI-1, a versatile serine protease inhibitor from the sea anemone Stichodactyla helianthus with high biomedical and biotechnological potential, toward elastase-like enzymes by substitution of the P1 residue (Lys(13)) with leucine. The variant (rShPI-1/K13L) exhibits a novel anti-porcine pancreatic elastase (PPE) activity together with a significantly improved inhibition of human neuthrophil elastase and chymotrypsin. The crystal structure of the PPE·rShPI-1/K13L complex determined at 2.0 Å resolution provided the first details of the canonical interaction between a BPTI-Kunitz-type domain and elastase-like enzymes. In addition to the essential impact of the variant P1 residue for complex stability, the interface is improved by increased contributions of the primary and secondary binding loop as compared with similar trypsin and chymotrypsin complexes. A comparison of the interaction network with elastase complexes of canonical inhibitors from the chelonian in family supports a key role of the P3 site in ShPI-1 in directing its selectivity against pancreatic and neutrophil elastases. Our results provide the structural basis for site-specific mutagenesis to further improve the binding affinity and/or direct the selectivity of BPTI-Kunitz-type inhibitors toward elastase-like enzymes.
Collapse
Affiliation(s)
- Rossana García-Fernández
- From the Centro de Estudio de Proteínas, Facultad de Biología, Universidad de la Habana, 20146 Habana, Cuba
| | - Markus Perbandt
- the Institute of Biochemistry and Molecular Biology, Department of Chemistry, University of Hamburg, 22761 Hamburg, Germany, the Hamburg Centre for Ultrafast Imaging, 22761 Hamburg, Germany, and
| | - Dirk Rehders
- the Joint Laboratory for Structural Biology of Infection and Inflammation, Institute of Biochemistry and Molecular Biology, University of Hamburg, 20146 Hamburg, Germany, and Institute of Biochemistry, University of Lübeck, c/o Deutsches Elektronen Synchrotron (DESY), 22603 Hamburg, Germany
| | - Patrick Ziegelmüller
- the Institute of Biochemistry and Molecular Biology, Department of Chemistry, University of Hamburg, 22761 Hamburg, Germany
| | - Nicolas Piganeau
- the Institute of Biochemistry and Molecular Biology, Department of Chemistry, University of Hamburg, 22761 Hamburg, Germany
| | - Ulrich Hahn
- the Institute of Biochemistry and Molecular Biology, Department of Chemistry, University of Hamburg, 22761 Hamburg, Germany
| | - Christian Betzel
- the Institute of Biochemistry and Molecular Biology, Department of Chemistry, University of Hamburg, 22761 Hamburg, Germany
| | | | - Lars Redecke
- the Joint Laboratory for Structural Biology of Infection and Inflammation, Institute of Biochemistry and Molecular Biology, University of Hamburg, 20146 Hamburg, Germany, and Institute of Biochemistry, University of Lübeck, c/o Deutsches Elektronen Synchrotron (DESY), 22603 Hamburg, Germany
| |
Collapse
|
34
|
A critical evaluation of in silico methods for detection of membrane protein intrinsic disorder. Biophys J 2014; 106:1638-49. [PMID: 24739163 DOI: 10.1016/j.bpj.2014.02.025] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/29/2013] [Revised: 02/03/2014] [Accepted: 02/25/2014] [Indexed: 11/23/2022] Open
Abstract
Intrinsically disordered regions in proteins possess important biological roles including transcriptional regulation, molecular recognition, and provision of sites for posttranslational modification. In three-dimensional crystallization of both soluble and membrane proteins, identification and removal of disordered regions is often necessary for obtaining crystals possessing sufficient long-range order for structure determination. Disordered regions can be identified experimentally, with techniques such as limited proteolysis coupled with mass spectrometry, or computationally, by using disorder prediction programs, of which many are available. Although these programs use various methods to predict disorder from a protein's primary sequence, they all were developed using information derived from soluble protein structures. Therefore, their performance and accuracy when applied to integral membrane proteins remained an open question. We evaluated the performance of 13 disorder prediction programs on a dataset containing 343 membrane proteins, and upon subdatasets containing only α-helical or β-barrel proteins. These programs were ranked using multiple metrics, including metrics specifically created for membrane proteins. Analysis of these data shows a clear distinction between programs that accurately predict disordered regions in membrane proteins and programs which perform poorly, and allows for the robust integration of in silico disorder prediction into our PSI:Biology membrane protein structural genomics pipeline.
Collapse
|
35
|
Inhester T, Rarey M. Protein-ligand interaction databases: advanced tools to mine activity data and interactions on a structural level. WILEY INTERDISCIPLINARY REVIEWS-COMPUTATIONAL MOLECULAR SCIENCE 2014. [DOI: 10.1002/wcms.1192] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/29/2022]
Affiliation(s)
- Therese Inhester
- Center for Bioinformatics; University of Hamburg; Hamburg Germany
| | - Matthias Rarey
- Center for Bioinformatics; University of Hamburg; Hamburg Germany
| |
Collapse
|
36
|
Solis AD. Deriving high-resolution protein backbone structure propensities from all crystal data using the information maximization device. PLoS One 2014; 9:e94334. [PMID: 24896099 PMCID: PMC4045576 DOI: 10.1371/journal.pone.0094334] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/22/2014] [Accepted: 03/12/2014] [Indexed: 11/28/2022] Open
Abstract
The most informative probability distribution functions (PDFs) describing the Ramachandran phi-psi dihedral angle pair, a fundamental descriptor of backbone conformation of protein molecules, are derived from high-resolution X-ray crystal structures using an information-theoretic approach. The Information Maximization Device (IMD) is established, based on fundamental information-theoretic concepts, and then applied specifically to derive highly resolved phi-psi maps for all 20 single amino acid and all 8000 triplet sequences at an optimal resolution determined by the volume of current data. The paper shows that utilizing the latent information contained in all viable high-resolution crystal structures found in the Protein Data Bank (PDB), totaling more than 77,000 chains, permits the derivation of a large number of optimized sequence-dependent PDFs. This work demonstrates the effectiveness of the IMD and the superiority of the resulting PDFs by extensive fold recognition experiments and rigorous comparisons with previously published triplet PDFs. Because it automatically optimizes PDFs, IMD results in improved performance of knowledge-based potentials, which rely on such PDFs. Furthermore, it provides an easy computational recipe for empirically deriving other kinds of sequence-dependent structural PDFs with greater detail and precision. The high-resolution phi-psi maps derived in this work are available for download.
Collapse
Affiliation(s)
- Armando D. Solis
- Biological Sciences Department, New York City College of Technology, The City University of New York, Brooklyn, New York, United States of America
- * E-mail:
| |
Collapse
|
37
|
Frappier V, Najmanovich RJ. A coarse-grained elastic network atom contact model and its use in the simulation of protein dynamics and the prediction of the effect of mutations. PLoS Comput Biol 2014; 10:e1003569. [PMID: 24762569 PMCID: PMC3998880 DOI: 10.1371/journal.pcbi.1003569] [Citation(s) in RCA: 72] [Impact Index Per Article: 7.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/20/2013] [Accepted: 02/25/2014] [Indexed: 11/18/2022] Open
Abstract
Normal mode analysis (NMA) methods are widely used to study dynamic aspects of protein structures. Two critical components of NMA methods are coarse-graining in the level of simplification used to represent protein structures and the choice of potential energy functional form. There is a trade-off between speed and accuracy in different choices. In one extreme one finds accurate but slow molecular-dynamics based methods with all-atom representations and detailed atom potentials. On the other extreme, fast elastic network model (ENM) methods with Cα-only representations and simplified potentials that based on geometry alone, thus oblivious to protein sequence. Here we present ENCoM, an Elastic Network Contact Model that employs a potential energy function that includes a pairwise atom-type non-bonded interaction term and thus makes it possible to consider the effect of the specific nature of amino-acids on dynamics within the context of NMA. ENCoM is as fast as existing ENM methods and outperforms such methods in the generation of conformational ensembles. Here we introduce a new application for NMA methods with the use of ENCoM in the prediction of the effect of mutations on protein stability. While existing methods are based on machine learning or enthalpic considerations, the use of ENCoM, based on vibrational normal modes, is based on entropic considerations. This represents a novel area of application for NMA methods and a novel approach for the prediction of the effect of mutations. We compare ENCoM to a large number of methods in terms of accuracy and self-consistency. We show that the accuracy of ENCoM is comparable to that of the best existing methods. We show that existing methods are biased towards the prediction of destabilizing mutations and that ENCoM is less biased at predicting stabilizing mutations.
Collapse
Affiliation(s)
- Vincent Frappier
- Department of Biochemistry, Faculty of Medicine and Health Sciences, Université de Sherbrooke, Sherbrooke, Quebec, Canada
| | - Rafael J Najmanovich
- Department of Biochemistry, Faculty of Medicine and Health Sciences, Université de Sherbrooke, Sherbrooke, Quebec, Canada
| |
Collapse
|
38
|
Abstract
By focusing on essential features, while averaging over less important details, coarse-grained (CG) models provide significant computational and conceptual advantages with respect to more detailed models. Consequently, despite dramatic advances in computational methodologies and resources, CG models enjoy surging popularity and are becoming increasingly equal partners to atomically detailed models. This perspective surveys the rapidly developing landscape of CG models for biomolecular systems. In particular, this review seeks to provide a balanced, coherent, and unified presentation of several distinct approaches for developing CG models, including top-down, network-based, native-centric, knowledge-based, and bottom-up modeling strategies. The review summarizes their basic philosophies, theoretical foundations, typical applications, and recent developments. Additionally, the review identifies fundamental inter-relationships among the diverse approaches and discusses outstanding challenges in the field. When carefully applied and assessed, current CG models provide highly efficient means for investigating the biological consequences of basic physicochemical principles. Moreover, rigorous bottom-up approaches hold great promise for further improving the accuracy and scope of CG models for biomolecular systems.
Collapse
Affiliation(s)
- W G Noid
- Department of Chemistry, The Pennsylvania State University, University Park, Pennsylvania 16802, USA
| |
Collapse
|
39
|
Hernández-Alcántara G, Torres-Larios A, Enríquez-Flores S, García-Torres I, Castillo-Villanueva A, Méndez ST, de la Mora-de la Mora I, Gómez-Manzo S, Torres-Arroyo A, López-Velázquez G, Reyes-Vivas H, Oria-Hernández J. Structural and functional perturbation of Giardia lamblia triosephosphate isomerase by modification of a non-catalytic, non-conserved region. PLoS One 2013; 8:e69031. [PMID: 23894402 PMCID: PMC3718800 DOI: 10.1371/journal.pone.0069031] [Citation(s) in RCA: 20] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/15/2013] [Accepted: 06/04/2013] [Indexed: 11/18/2022] Open
Abstract
BACKGROUND We have previously proposed triosephosphate isomerase of Giardia lamblia (GlTIM) as a target for rational drug design against giardiasis, one of the most common parasitic infections in humans. Since the enzyme exists in the parasite and the host, selective inhibition is a major challenge because essential regions that could be considered molecular targets are highly conserved. Previous biochemical evidence showed that chemical modification of the non-conserved non-catalytic cysteine 222 (C222) inactivates specifically GlTIM. The inactivation correlates with the physicochemical properties of the modifying agent: addition of a non-polar, small chemical group at C222 reduces the enzyme activity by one half, whereas negatively charged, large chemical groups cause full inactivation. RESULTS In this work we used mutagenesis to extend our understanding of the functional and structural effects triggered by modification of C222. To this end, six GlTIM C222 mutants with side chains having diverse physicochemical characteristics were characterized. We found that the polarity, charge and volume of the side chain in the mutant amino acid differentially alter the activity, the affinity, the stability and the structure of the enzyme. The data show that mutagenesis of C222 mimics the effects of chemical modification. The crystallographic structure of C222D GlTIM shows the disruptive effects of introducing a negative charge at position 222: the mutation perturbs loop 7, a region of the enzyme whose interactions with the catalytic loop 6 are essential for TIM stability, ligand binding and catalysis. The amino acid sequence of TIM in phylogenetic diverse groups indicates that C222 and its surrounding residues are poorly conserved, supporting the proposal that this region is a good target for specific drug design. CONCLUSIONS The results demonstrate that it is possible to inhibit species-specifically a ubiquitous, structurally highly conserved enzyme by modification of a non-conserved, non-catalytic residue through long-range perturbation of essential regions.
Collapse
Affiliation(s)
- Gloria Hernández-Alcántara
- Laboratorio de Bioquímica-Genética, Instituto Nacional de Pediatría, Secretaría de Salud, Mexico City, Mexico
| | - Alfredo Torres-Larios
- Instituto de Fisiología Celular, Universidad Nacional Autónoma de México, Mexico City, Mexico
| | - Sergio Enríquez-Flores
- Laboratorio de Bioquímica-Genética, Instituto Nacional de Pediatría, Secretaría de Salud, Mexico City, Mexico
| | - Itzhel García-Torres
- Laboratorio de Bioquímica-Genética, Instituto Nacional de Pediatría, Secretaría de Salud, Mexico City, Mexico
| | - Adriana Castillo-Villanueva
- Laboratorio de Bioquímica-Genética, Instituto Nacional de Pediatría, Secretaría de Salud, Mexico City, Mexico
| | - Sara T. Méndez
- Laboratorio de Bioquímica-Genética, Instituto Nacional de Pediatría, Secretaría de Salud, Mexico City, Mexico
| | | | - Saúl Gómez-Manzo
- Laboratorio de Bioquímica-Genética, Instituto Nacional de Pediatría, Secretaría de Salud, Mexico City, Mexico
| | - Angélica Torres-Arroyo
- Laboratorio de Bioquímica-Genética, Instituto Nacional de Pediatría, Secretaría de Salud, Mexico City, Mexico
| | - Gabriel López-Velázquez
- Laboratorio de Bioquímica-Genética, Instituto Nacional de Pediatría, Secretaría de Salud, Mexico City, Mexico
| | - Horacio Reyes-Vivas
- Laboratorio de Bioquímica-Genética, Instituto Nacional de Pediatría, Secretaría de Salud, Mexico City, Mexico
- * E-mail: (JOH); (HRV)
| | - Jesús Oria-Hernández
- Laboratorio de Bioquímica-Genética, Instituto Nacional de Pediatría, Secretaría de Salud, Mexico City, Mexico
- * E-mail: (JOH); (HRV)
| |
Collapse
|
40
|
Ezhevskaya M, Bordignon E, Polyhach Y, Moens L, Dewilde S, Jeschke G, Van Doorslaer S. Distance determination between low-spin ferric haem and nitroxide spin label using DEER: the neuroglobin case. Mol Phys 2013. [DOI: 10.1080/00268976.2013.813592] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022]
|
41
|
Abaturov LV, Nosova NG. Structure of crambin in solution, crystal and in the trajectories of molecular dynamics simulations. Biophysics (Nagoya-shi) 2013. [DOI: 10.1134/s0006350913030020] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/23/2022] Open
|
42
|
Garcia BL, Summers BJ, Ramyar KX, Tzekou A, Lin Z, Ricklin D, Lambris JD, Laity JH, Geisbrecht BV. A structurally dynamic N-terminal helix is a key functional determinant in staphylococcal complement inhibitor (SCIN) proteins. J Biol Chem 2013; 288:2870-81. [PMID: 23233676 PMCID: PMC3554951 DOI: 10.1074/jbc.m112.426858] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/08/2023] Open
Abstract
Complement is a network of interacting circulatory and cell surface proteins that recognizes, marks, and facilitates clearance of microbial invaders. To evade complement attack, the pathogenic organism Staphylococcus aureus expresses a number of secreted proteins that interfere with activation and regulation of the complement cascade. Staphylococcal complement inhibitors (SCINs) are one important class of these immunomodulators and consist of three active members (SCIN-A/-B/-C). SCINs inhibit a critical enzymatic complex, the alternative pathway C3 convertase, by targeting a functional "hot spot" on the central opsonin of complement, C3b. Although N-terminal truncation mutants of SCINs retain complement inhibitory properties, they are significantly weaker binders of C3b. To provide a structural basis for this observation, we undertook a series of crystallographic and NMR dynamics studies on full-length SCINs. This work reveals that N-terminal SCIN domains are characterized by a conformationally dynamic helical motif. C3b binding and functional experiments further demonstrate that this sequence-divergent N-terminal region of SCINs is both functionally important and context-dependent. Finally, surface plasmon resonance data provide evidence for the formation of inhibitor·enzyme·substrate complexes ((SCIN·C3bBb)·C3). Similar to the (SCIN·C3bBb)(2) pseudodimeric complexes, ((SCIN·C3bBb)·C3) interferes with the interaction of complement receptors and C3b. This activity provides an additional mechanism by which SCIN couples convertase inhibition to direct blocking of phagocytosis. Together, these data suggest that tethering multi-host protein complexes by small modular bacterial inhibitors may be a global strategy of immune evasion used by S. aureus. The work presented here provides detailed structure-activity relationships and improves our understanding of how S. aureus circumvents human innate immunity.
Collapse
Affiliation(s)
- Brandon L Garcia
- Division of Cell Biology and Biophysics, School of Biological Sciences, University of Missouri, Kansas City, Missouri 64110, USA
| | | | | | | | | | | | | | | | | |
Collapse
|
43
|
Computational study and peptide inhibitors design for the CDK9 - cyclin T1 complex. J Mol Model 2013; 19:1711-25. [PMID: 23296566 DOI: 10.1007/s00894-012-1735-2] [Citation(s) in RCA: 10] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/17/2012] [Accepted: 12/13/2012] [Indexed: 01/08/2023]
Abstract
Cyclin dependent kinase 9 (CDK9) is a protein that belongs to the cyclin-dependent kinases family, and its main role is in the regulation of the cell transcription processes. Since the increased activity of CDK9 is connected with the development of pathological processes such as tumor growth and survival and HIV-1 replication, inhibition of the CDK9 could be of particular interest for treating such diseases. The activation of CDK9 is initiated by the formation of CDK9/cyclin T1 complex, therefore disruption of its formation could be a promising strategy for the design of CDK9 inhibitors. In order to assist in the design of potential inhibitors of CDK9/cyclin T1 complex formation, a computational study of the CDK9/cyclin T1 interface was conducted. Ten peptides were designed using the information from the analysis of the complex, hot spot residues and fragment based design. The designed peptides were docked to CDK9 structures obtained by molecular dynamics simulations of CDK9/cyclin T1 complex and the CDK9 alone and their binding affinities were evaluated using molecular mechanics Poisson Boltzman surface area (MM-PBSA) method and steered molecular dynamics (SMD). Designed peptide sequences LQTLGF and ESIILQ, both derived from the surface of cyclin T1, as well as the peptide sequence PRWPE, derived from fragment based design, showed the most favorable binding properties and were selected for our further studies.
Collapse
|
44
|
Ruvinsky AM, Kirys T, Tuzikov AV, Vakser IA. Structure fluctuations and conformational changes in protein binding. J Bioinform Comput Biol 2012; 10:1241002. [PMID: 22809338 DOI: 10.1142/s0219720012410028] [Citation(s) in RCA: 28] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022]
Abstract
Structure fluctuations and conformational changes accompany all biological processes involving macromolecules. The paper presents a classification of protein residues based on the normalized equilibrium fluctuations of the residue centers of mass in proteins and a statistical analysis of conformation changes in the side-chains upon binding. Normal mode analysis and an elastic network model were applied to a set of protein complexes to calculate the residue fluctuations and develop the residue classification. Comparison with a classification based on normalized B-factors suggests that the B-factors may underestimate protein flexibility in solvent. Our classification shows that protein loops and disordered fragments are enriched with highly fluctuating residues and depleted with weakly fluctuating residues. Strategies for engineering thermostable proteins are discussed. To calculate the dihedral angles distribution functions, the configuration space was divided into cells by a cubic grid. The effect of protein association on the distribution functions depends on the amino acid type and a grid step in the dihedral angles space. The changes in the dihedral angles increase from the near-backbone dihedral angle to the most distant one, for most residues. On average, one fifth of the interface residues change the rotamer state upon binding, whereas the rest of the interface residues undergo local readjustments within the same rotamer.
Collapse
Affiliation(s)
- Anatoly M Ruvinsky
- Center for Bioinformatics, University of Kansas, Lawrence, KS 66047, USA.
| | | | | | | |
Collapse
|
45
|
Dilanian RA, Streltsov VA, Quiney HM, Nugent KA. Continuous X-ray diffractive field in protein nanocrystallography. Acta Crystallogr A 2012; 69:108-18. [PMID: 23250067 DOI: 10.1107/s0108767312042535] [Citation(s) in RCA: 18] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/11/2012] [Accepted: 10/11/2012] [Indexed: 11/11/2022] Open
Abstract
The recent development of X-ray free-electron laser sources has created new opportunities for the structural analysis of protein nanocrystals. The extremely small sizes of the crystals, as well as imperfections of the crystal structure, result in an interference phenomenon in the diffraction pattern. With decreasing crystallite size the structural imperfections play a role in the formation of the diffraction pattern that is comparable in importance to the size effects and should be taken into account during the data analysis and structure reconstruction processes. There now exists a need to develop new methods of protein structure determination that do not depend on the availability of good-quality crystals and that can treat proteins under conditions close to the active form. This paper demonstrates an approach that is specifically tailored to nanocrystalline samples and offers a unique crystallographic solution.
Collapse
Affiliation(s)
- Ruben A Dilanian
- ARC Centre of Excellence for Coherent X-ray Science, School of Physics, The University of Melbourne, Victoria, Australia.
| | | | | | | |
Collapse
|
46
|
Proliferating cell nuclear antigen (PCNA) interactions in solution studied by NMR. PLoS One 2012; 7:e48390. [PMID: 23139781 PMCID: PMC3491057 DOI: 10.1371/journal.pone.0048390] [Citation(s) in RCA: 21] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/25/2012] [Accepted: 09/24/2012] [Indexed: 01/25/2023] Open
Abstract
PCNA is an essential factor for DNA replication and repair. It forms a ring shaped structure of 86 kDa by the symmetric association of three identical protomers. The ring encircles the DNA and acts as a docking platform for other proteins, most of them containing the PCNA Interaction Protein sequence (PIP-box). We have used NMR to characterize the interactions of PCNA with several other proteins and fragments in solution. The binding of the PIP-box peptide of the cell cycle inhibitor p21 to PCNA is consistent with the crystal structure of the complex. A shorter p21 peptide binds with reduced affinity but retains most of the molecular recognition determinants. However the binding of the corresponding peptide of the tumor suppressor ING1 is extremely weak, indicating that slight deviations from the consensus PIP-box sequence dramatically reduce the affinity for PCNA, in contrast with a proposed less stringent PIP-box sequence requirement. We could not detect any binding between PCNA and the MCL-1 or the CDK2 protein, reported to interact with PCNA in biochemical assays. This suggests that they do not bind directly to PCNA, or they do but very weakly, with additional unidentified factors stabilizing the interactions in the cell. Backbone dynamics measurements show three PCNA regions with high relative flexibility, including the interdomain connector loop (IDCL) and the C-terminus, both of them involved in the interaction with the PIP-box. Our work provides the basis for high resolution studies of direct ligand binding to PCNA in solution.
Collapse
|
47
|
Kirys T, Ruvinsky AM, Tuzikov AV, Vakser IA. Correlation analysis of the side-chains conformational distribution in bound and unbound proteins. BMC Bioinformatics 2012; 13:236. [PMID: 22984947 PMCID: PMC3479416 DOI: 10.1186/1471-2105-13-236] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/23/2012] [Accepted: 09/11/2012] [Indexed: 01/30/2023] Open
Abstract
BACKGROUND Protein interactions play a key role in life processes. Characterization of conformational properties of protein-protein interactions is important for understanding the mechanisms of protein association. The rapidly increasing amount of experimentally determined structures of proteins and protein-protein complexes provides foundation for research on protein interactions and complex formation. The knowledge of the conformations of the surface side chains is essential for modeling of protein complexes. The purpose of this study was to analyze and compare dihedral angle distribution functions of the side chains at the interface and non-interface areas in bound and unbound proteins. RESULTS To calculate the dihedral angle distribution functions, the configuration space was divided into grid cells. Statistical analysis showed that the similarity between bound and unbound interface and non-interface surface depends on the amino acid type and the grid resolution. The correlation coefficients between the distribution functions increased with the grid spacing increase for all amino acid types. The Manhattan distance showing the degree of dissimilarity between the distribution functions decreased accordingly. Short residues with one or two dihedral angles had higher correlations and smaller Manhattan distances than the longer residues. Met and Arg had the slowest growth of the correlation coefficient with the grid spacing increase. The correlations between the interface and non-interface distribution functions had a similar dependence on the grid resolution in both bound and unbound states. The interface and non-interface differences between bound and unbound distribution functions, caused by biological protein-protein interactions or crystal contacts, disappeared at the 70° grid spacing for interfaces and 30° for non-interface surface, which agrees with an average span of the side-chain rotamers. CONCLUSIONS The two-fold difference in the critical grid spacing indicates larger conformational changes upon binding at the interface than at the rest of the surface. At the same time, transitions between rotamers induced by interactions across the interface or the crystal packing are rare, with most side chains having local readjustments that do not change the rotameric state. The analysis is important for better understanding of protein interactions and development of flexible docking approaches.
Collapse
Affiliation(s)
- Tatsiana Kirys
- Center for Bioinformatics, The University of Kansas, Lawrence, KS 66047, USA
| | | | | | | |
Collapse
|
48
|
Xin F, Radivojac P. Post-translational modifications induce significant yet not extreme changes to protein structure. ACTA ACUST UNITED AC 2012; 28:2905-13. [PMID: 22947645 DOI: 10.1093/bioinformatics/bts541] [Citation(s) in RCA: 99] [Impact Index Per Article: 8.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022]
Abstract
MOTIVATION A number of studies of individual proteins have shown that post-translational modifications (PTMs) are associated with structural rearrangements of their target proteins. Although such studies provide critical insights into the mechanics behind the dynamic regulation of protein function, they usually feature examples with relatively large conformational changes. However, with the steady growth of Protein Data Bank (PDB) and available PTM sites, it is now possible to more systematically characterize the role of PTMs as conformational switches. In this study, we ask (1) what is the expected extent of structural change upon PTM, (2) how often are those changes in fact substantial, (3) whether the structural impact is spatially localized or global and (4) whether different PTMs have different signatures. RESULTS We exploit redundancy in PDB and, using root-mean-square deviation, study the conformational heterogeneity of groups of protein structures corresponding to identical sequences in their unmodified and modified forms. We primarily focus on the two most abundant PTMs in PDB, glycosylation and phosphorylation, but show that acetylation and methylation have similar tendencies. Our results provide evidence that PTMs induce conformational changes at both local and global level. However, the proportion of large changes is unexpectedly small; only 7% of glycosylated and 13% of phosphorylated proteins undergo global changes >2 Å. Further analysis suggests that phosphorylation stabilizes protein structure by reducing global conformational heterogeneity by 25%. Overall, these results suggest a subtle but common role of allostery in the mechanisms through which PTMs affect regulatory and signaling pathways. CONTACT predrag@indiana.edu. SUPPLEMENTARY INFORMATION Supplementary data are available at Bioinformatics online.
Collapse
Affiliation(s)
- Fuxiao Xin
- School of Informatics and Computing, Indiana University, Bloomington, IN, USA
| | | |
Collapse
|
49
|
Abstract
Structure-based drug design has become an essential tool for rapid lead discovery and optimization. As available structural information has increased, researchers have become increasingly aware of the importance of protein flexibility for accurate description of the native state. Typical protein-ligand docking efforts still rely on a single rigid receptor, which is an incomplete representation of potential binding conformations of the protein. These rigid docking efforts typically show the best performance rates between 50 and 75%, while fully flexible docking methods can enhance pose prediction up to 80-95%. This review examines the current toolbox for flexible protein-ligand docking and receptor surface mapping. Present limitations and possibilities for future development are discussed.
Collapse
Affiliation(s)
- Katrina W. Lexa
- Department of Medicinal Chemistry, University of Michigan, 428 Church Street, Ann Arbor, MI 48109-1065, USA
| | - Heather A. Carlson
- Department of Medicinal Chemistry, University of Michigan, 428 Church Street, Ann Arbor, MI 48109-1065, USA
| |
Collapse
|
50
|
Issack BB, Berjanskii M, Wishart DS, Stepanova M. Exploring the essential collective dynamics of interacting proteins: application to prion protein dimers. Proteins 2012; 80:1847-65. [PMID: 22488640 DOI: 10.1002/prot.24082] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/30/2011] [Revised: 03/07/2012] [Accepted: 03/18/2012] [Indexed: 11/11/2022]
Abstract
Essential collective dynamics is a promising and robust approach for analysing the slow motions of macromolecules from short molecular dynamics trajectories. In this study, an extension of the method to treat a collection of interacting protein molecules is presented. The extension is applied to investigate the effects of dimerization on the collective dynamics of ovine prion protein molecules in two different arrangements. Examination of the structural plasticity shows that aggregation has a restricting effect on the local mobility of the prion protein molecules in the interfacial regions. Domain motions of the two dimeric ovine prion protein conformations are distinctly different and can be related to interatomic correlations at the interface. Correlated motions are among the slow collective modes extensively analysed by considering both main-chain and side-chain atoms. Correlation maps reveal the existence of a vast network of dynamically correlated side groups, which extends beyond individual subunits via interfacial interconnections. The network is formed by a core of hydrophobic side chains surrounded by hydrophilic groups at the periphery. The relevance of these findings are discussed in the context of mutations associated with prion diseases. The binding free energy of the dimeric conformations is evaluated to probe their thermodynamic stability. The descriptions afforded by the analysis of the essential collective dynamics of the prion dimers are consistent with their binding free energies. The agreement validates the extension of the methodology and provides a means of interpreting the collective dynamics in terms of the thermodynamic stability of ovine prion proteins.
Collapse
Affiliation(s)
- Bilkiss B Issack
- National Institute for Nanotechnology, National research Council, Edmonton, AB, Canada
| | | | | | | |
Collapse
|