1
|
Yu B, Bolik-Coulon N, Rangadurai AK, Kay LE, Iwahara J. Gadolinium-Based NMR Spin Relaxation Measurements of Near-Surface Electrostatic Potentials of Biomolecules. J Am Chem Soc 2024; 146:20788-20801. [PMID: 39028837 PMCID: PMC11295196 DOI: 10.1021/jacs.4c04433] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/31/2024] [Revised: 06/09/2024] [Accepted: 06/28/2024] [Indexed: 07/21/2024]
Abstract
NMR spectroscopy is an important tool for the measurement of the electrostatic properties of biomolecules. To this point, paramagnetic relaxation enhancements (PREs) of 1H nuclei arising from nitroxide cosolutes in biomolecular solutions have been used to measure effective near-surface electrostatic potentials (ϕENS) of proteins and nucleic acids. Here, we present a gadolinium (Gd)-based NMR method, exploiting Gd chelates with different net charges, for measuring ϕENS values and demonstrate its utility through applications to a number of biomolecular systems. The use of Gd-based cosolutes offers several advantages over nitroxides for ϕENS measurements. First, unlike nitroxide compounds, Gd chelates enable electrostatic potential measurements on oxidation-sensitive proteins that require reducing agents. Second, the large electron spin quantum number of Gd (7/2) results in notably larger PREs for Gd chelates when used at the same concentrations as nitroxide radicals. Thus, it is possible to measure ϕENS values exclusively from + and - charged compounds even for highly charged biomolecules, avoiding the use of neutral cosolutes that, as we further establish here, limits the accuracy of the measured electrostatic potentials. In addition, the smaller concentrations of cosolutes required minimize potential binding to sites on macromolecules. Fourth, the closer proximity of the paramagnetic center and charged groups within Gd chelates, in comparison to the corresponding nitroxide compounds, enables more accurate predictions of ϕENS potentials for cross-validation of the experimental results. The Gd-based method described here, thus, broadens the applicability of studies of biomolecular electrostatics using solution NMR spectroscopy.
Collapse
Affiliation(s)
- Binhan Yu
- Department
of Biochemistry & Molecular Biology, Sealy Center for Structural
Biology & Molecular Biophysics, University
of Texas Medical Branch, Galveston, Texas 77555-1068, United States
| | - Nicolas Bolik-Coulon
- Department
of Molecular Genetics, University of Toronto, Toronto, Ontario M5S 1A8, Canada
- Department
of Chemistry, University of Toronto, Toronto, Ontario M5S 3H6, Canada
- Department
of Biochemistry, University of Toronto, Toronto, Ontario M5S 3H6, Canada
| | - Atul K. Rangadurai
- Department
of Molecular Genetics, University of Toronto, Toronto, Ontario M5S 1A8, Canada
- Department
of Chemistry, University of Toronto, Toronto, Ontario M5S 3H6, Canada
- Department
of Biochemistry, University of Toronto, Toronto, Ontario M5S 3H6, Canada
- Program
in Molecular Medicine, Hospital for Sick
Children Research Institute, Toronto, Ontario M5G 0A4, Canada
| | - Lewis E. Kay
- Department
of Molecular Genetics, University of Toronto, Toronto, Ontario M5S 1A8, Canada
- Department
of Chemistry, University of Toronto, Toronto, Ontario M5S 3H6, Canada
- Department
of Biochemistry, University of Toronto, Toronto, Ontario M5S 3H6, Canada
- Program
in Molecular Medicine, Hospital for Sick
Children Research Institute, Toronto, Ontario M5G 0A4, Canada
| | - Junji Iwahara
- Department
of Biochemistry & Molecular Biology, Sealy Center for Structural
Biology & Molecular Biophysics, University
of Texas Medical Branch, Galveston, Texas 77555-1068, United States
| |
Collapse
|
2
|
Iwahara J, Pettitt BM, Yu B. Direct measurements of biomolecular electrostatics through experiments. Curr Opin Struct Biol 2023; 82:102680. [PMID: 37573815 PMCID: PMC10947535 DOI: 10.1016/j.sbi.2023.102680] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/28/2023] [Revised: 07/19/2023] [Accepted: 07/19/2023] [Indexed: 08/15/2023]
Abstract
Biomolecular electrostatics has been a subject of computational investigations based on 3D structures. This situation is changing because emerging experimental tools allow us to quantitatively investigate biomolecular electrostatics without any use of structure information. Now, electrostatic potentials around biomolecules can directly be measured for many residues simultaneously by nuclear magnetic resonance (NMR) spectroscopy. This NMR method can be used to study electrostatic aspects of various processes, including macromolecular association and liquid-liquid phase separation. Applications to structurally flexible biomolecules such as intrinsically disordered proteins are particularly useful. The new tools also facilitate examination of theoretical models and methods for biomolecular electrostatics.
Collapse
Affiliation(s)
- Junji Iwahara
- Department of Biochemistry & Molecular Biology, Sealy Center for Structural Biology & Molecular Biophysics, University of Texas Medical Branch, Galveston, TX 77555, USA.
| | - B Montgomery Pettitt
- Department of Biochemistry & Molecular Biology, Sealy Center for Structural Biology & Molecular Biophysics, University of Texas Medical Branch, Galveston, TX 77555, USA
| | - Binhan Yu
- Department of Biochemistry & Molecular Biology, Sealy Center for Structural Biology & Molecular Biophysics, University of Texas Medical Branch, Galveston, TX 77555, USA
| |
Collapse
|
3
|
Chen WH, Strych U, Bottazzi ME, Lin YP. Past, present, and future of Lyme disease vaccines: antigen engineering approaches and mechanistic insights. Expert Rev Vaccines 2022; 21:1405-1417. [PMID: 35836340 DOI: 10.1080/14760584.2022.2102484] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/04/2022]
Abstract
Introduction: Transmitted by ticks, Lyme disease is the most common vector-borne disease in the Northern hemisphere. Despite the geographical expansion of human Lyme disease cases, no effective preventive strategies are currently available. Developing an efficacious and safe vaccine is therefore urgently needed. Efforts have previously been taken to identify vaccine targets in the causative pathogen (Borrelia burgdorferi sensu lato) and arthropod vector (Ixodes spp.). However, progress was impeded due to a lack of consumer confidence caused by the myth of undesired off-target responses, low immune responses, a limited breadth of immune reactivity, as well as by the complexities of the vaccine process development.Area covered: In this review, we summarize the antigen engineering approaches that have been applied to overcome those challenges and the underlying mechanisms that can be exploited to improve both safety and efficacy of future Lyme disease vaccines.Expert opinion: Over the past two decades, several new genetically redesigned Lyme disease vaccine candidates have shown success in both preclinical and clinical settings and built a solid foundation for further development. These studies have greatly informed the protective mechanisms of reducing Lyme disease burdens and ending the endemic of this disease.
Collapse
Affiliation(s)
- Wen-Hsiang Chen
- Department of Pediatrics, National School of Tropical Medicine, Baylor College of Medicine, Houston, TX, USA.,Texas Children's Hospital Center for Vaccine Development, Houston, TX, USA
| | - Ulrich Strych
- Department of Pediatrics, National School of Tropical Medicine, Baylor College of Medicine, Houston, TX, USA.,Texas Children's Hospital Center for Vaccine Development, Houston, TX, USA
| | - Maria Elena Bottazzi
- Department of Pediatrics, National School of Tropical Medicine, Baylor College of Medicine, Houston, TX, USA.,Texas Children's Hospital Center for Vaccine Development, Houston, TX, USA.,Department of Biology, Baylor University, Waco, TX, United States
| | - Yi-Pin Lin
- Division of Infectious Diseases, Wadsworth Center, NYSDOH, Albany, NY, USA.,Department of Biomedical Sciences, SUNY Albany, Albany, NY, USA
| |
Collapse
|
4
|
Yu B, Pletka CC, Iwahara J. Protein Electrostatics Investigated through Paramagnetic NMR for Nonpolar Groups. J Phys Chem B 2022; 126:2196-2202. [PMID: 35266708 PMCID: PMC8973454 DOI: 10.1021/acs.jpcb.1c10930] [Citation(s) in RCA: 6] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Experimental validation of theoretical models for protein electrostatics remains rare. Recently, we have developed a paramagnetic NMR-based method for de novo determination of effective near-surface electrostatic potentials, which allows for straightforward examination of electrostatic models for biomolecules. In the current work, we expand this method and demonstrate that effective near-surface electrostatic potentials can readily be determined from 1H paramagnetic relaxation enhancement (PRE) data for protein CαH and CH3 groups. The experimental data were compared with those predicted from the Poisson-Boltzmann theory. The impact of structural dynamics on the effective near-surface electrostatic potentials was also assessed. The agreement between the experimental and theoretical data was particularly good for methyl 1H nuclei. Compared to the conventional pKa-based validation, our paramagnetic NMR-based approach can provide a far larger number of experimental data that can directly be used to examine the validity of theoretical electrostatic models for proteins.
Collapse
|
5
|
Liang S, Li Z, Zhan J, Zhou Y. De novo protein design by an energy function based on series expansion in distance and orientation dependence. Bioinformatics 2021; 38:86-93. [PMID: 34406339 DOI: 10.1093/bioinformatics/btab598] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/05/2021] [Revised: 08/11/2021] [Accepted: 08/16/2021] [Indexed: 02/03/2023] Open
Abstract
MOTIVATION Despite many successes, de novo protein design is not yet a solved problem as its success rate remains low. The low success rate is largely because we do not yet have an accurate energy function for describing the solvent-mediated interaction between amino acid residues in a protein chain. Previous studies showed that an energy function based on series expansions with its parameters optimized for side-chain and loop conformations can lead to one of the most accurate methods for side chain (OSCAR) and loop prediction (LEAP). Following the same strategy, we developed an energy function based on series expansions with the parameters optimized in four separate stages (recovering single-residue types without and with orientation dependence, selecting loop decoys and maintaining the composition of amino acids). We tested the energy function for de novo design by using Monte Carlo simulated annealing. RESULTS The method for protein design (OSCAR-Design) is found to be as accurate as OSCAR and LEAP for side-chain and loop prediction, respectively. In de novo design, it can recover native residue types ranging from 38% to 43% depending on test sets, conserve hydrophobic/hydrophilic residues at ∼75%, and yield the overall similarity in amino acid compositions at more than 90%. These performance measures are all statistically significantly better than several protein design programs compared. Moreover, the largest hydrophobic patch areas in designed proteins are near or smaller than those in native proteins. Thus, an energy function based on series expansion can be made useful for protein design. AVAILABILITY AND IMPLEMENTATION The Linux executable version is freely available for academic users at http://zhouyq-lab.szbl.ac.cn/resources/.
Collapse
Affiliation(s)
- Shide Liang
- Department of R & D, Bio-Thera Solutions, Guangzhou 510530, China
| | - Zhixiu Li
- Institute of Health and Biomedical Innovation, Queensland University of Technology at Translational Research Institute, Woolloongabba, QLD 3001, Australia
| | - Jian Zhan
- Institute for Glycomics and School of Information and Communication Technology, Griffith University, Gold Coast Campus, Southport, QLD 4222, Australia.,Institute for Systems and Physical Biology, Shenzhen Bay Laboratory, Shenzhen 518055, China
| | - Yaoqi Zhou
- Institute for Systems and Physical Biology, Shenzhen Bay Laboratory, Shenzhen 518055, China.,Peking University Shenzhen Graduate School, Shenzhen 518055, China
| |
Collapse
|
6
|
Abstract
Electrostatic potentials computed from three-dimensional structures of biomolecules by solving the Poisson-Boltzmann equation are widely used in molecular biophysics, structural biology, and medicinal chemistry. Despite the approximate nature of the Poisson-Boltzmann theory, validation of the computed electrostatic potentials around biological macromolecules is rare and methodologically limited. Here, we present a unique and powerful NMR method that allows for straightforward and extensive comparison with electrostatic models for biomolecules and their complexes. This method utilizes paramagnetic relaxation enhancement arising from analogous cationic and anionic cosolutes whose spatial distributions around biological macromolecules reflect electrostatic potentials. We demonstrate that this NMR method enables de novo determination of near-surface electrostatic potentials for individual protein residues without using any structural information. We applied the method to ubiquitin and the Antp homeodomain-DNA complex. The experimental data agreed well with predictions from the Poisson-Boltzmann theory. Thus, our experimental results clearly support the validity of the theory for these systems. However, our experimental study also illuminates certain weaknesses of the Poisson-Boltzmann theory. For example, we found that the theory predicts stronger dependence of near-surface electrostatic potentials on ionic strength than observed in the experiments. Our data also suggest that conformational flexibility or structural uncertainties may cause large errors in theoretical predictions of electrostatic potentials, particularly for highly charged systems. This NMR-based method permits extensive assessment of near-surface electrostatic potentials for various regions around biological macromolecules and thereby may facilitate improvement of the computational approaches for electrostatic potentials.
Collapse
|
7
|
Sacchi S, Rabattoni V, Miceli M, Pollegioni L. Yin and Yang in Post-Translational Modifications of Human D-Amino Acid Oxidase. Front Mol Biosci 2021; 8:684934. [PMID: 34041270 PMCID: PMC8141710 DOI: 10.3389/fmolb.2021.684934] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/24/2021] [Accepted: 04/23/2021] [Indexed: 11/30/2022] Open
Abstract
In the central nervous system, the flavoprotein D-amino acid oxidase is responsible for catabolizing D-serine, the main endogenous coagonist of N-methyl-D-aspartate receptor. Dysregulation of D-serine brain levels in humans has been associated with neurodegenerative and psychiatric disorders. This D-amino acid is synthesized by the enzyme serine racemase, starting from the corresponding L-enantiomer, and degraded by both serine racemase (via an elimination reaction) and the flavoenzyme D-amino acid oxidase. To shed light on the role of human D-amino acid oxidase (hDAAO) in D-serine metabolism, the structural/functional relationships of this enzyme have been investigated in depth and several strategies aimed at controlling the enzymatic activity have been identified. Here, we focused on the effect of post-translational modifications: by using a combination of structural analyses, biochemical methods, and cellular studies, we investigated whether hDAAO is subjected to nitrosylation, sulfhydration, and phosphorylation. hDAAO is S-nitrosylated and this negatively affects its activity. In contrast, the hydrogen sulfide donor NaHS seems to alter the enzyme conformation, stabilizing a species with higher affinity for the flavin adenine dinucleotide cofactor and thus positively affecting enzymatic activity. Moreover, hDAAO is phosphorylated in cerebellum; however, the protein kinase involved is still unknown. Taken together, these findings indicate that D-serine levels can be also modulated by post-translational modifications of hDAAO as also known for the D-serine synthetic enzyme serine racemase.
Collapse
Affiliation(s)
- Silvia Sacchi
- "The Protein Factory 2.0", Dipartimento di Biotecnologie e Scienze della Vita, Università degli Studi Dell'Insubria, Varese, Italy
| | - Valentina Rabattoni
- "The Protein Factory 2.0", Dipartimento di Biotecnologie e Scienze della Vita, Università degli Studi Dell'Insubria, Varese, Italy
| | - Matteo Miceli
- "The Protein Factory 2.0", Dipartimento di Biotecnologie e Scienze della Vita, Università degli Studi Dell'Insubria, Varese, Italy
| | - Loredano Pollegioni
- "The Protein Factory 2.0", Dipartimento di Biotecnologie e Scienze della Vita, Università degli Studi Dell'Insubria, Varese, Italy
| |
Collapse
|
8
|
Mignon D, Druart K, Michael E, Opuu V, Polydorides S, Villa F, Gaillard T, Panel N, Archontis G, Simonson T. Physics-Based Computational Protein Design: An Update. J Phys Chem A 2020; 124:10637-10648. [DOI: 10.1021/acs.jpca.0c07605] [Citation(s) in RCA: 9] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/06/2023]
Affiliation(s)
- David Mignon
- Laboratoire de Biologie Structurale de la Cellule (CNRS UMR7654), Ecole Polytechnique, 91128 Palaiseau, France
| | - Karen Druart
- Laboratoire de Biologie Structurale de la Cellule (CNRS UMR7654), Ecole Polytechnique, 91128 Palaiseau, France
| | - Eleni Michael
- Department of Physics, University of Cyprus, PO20537, CY1678 Nicosia, Cyprus
| | - Vaitea Opuu
- Laboratoire de Biologie Structurale de la Cellule (CNRS UMR7654), Ecole Polytechnique, 91128 Palaiseau, France
| | - Savvas Polydorides
- Department of Physics, University of Cyprus, PO20537, CY1678 Nicosia, Cyprus
| | - Francesco Villa
- Laboratoire de Biologie Structurale de la Cellule (CNRS UMR7654), Ecole Polytechnique, 91128 Palaiseau, France
| | - Thomas Gaillard
- Laboratoire de Biologie Structurale de la Cellule (CNRS UMR7654), Ecole Polytechnique, 91128 Palaiseau, France
| | - Nicolas Panel
- Laboratoire de Biologie Structurale de la Cellule (CNRS UMR7654), Ecole Polytechnique, 91128 Palaiseau, France
| | - Georgios Archontis
- Department of Physics, University of Cyprus, PO20537, CY1678 Nicosia, Cyprus
| | - Thomas Simonson
- Laboratoire de Biologie Structurale de la Cellule (CNRS UMR7654), Ecole Polytechnique, 91128 Palaiseau, France
| |
Collapse
|
9
|
Gütlein P, Blumberger J, Oberhofer H. An Iterative Fragment Scheme for the ACKS2 Electronic Polarization Model: Application to Molecular Dimers and Chains. J Chem Theory Comput 2020; 16:5723-5735. [PMID: 32701273 DOI: 10.1021/acs.jctc.0c00151] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/14/2023]
Abstract
The treatment of electrostatic interactions is a key ingredient in the force field-based simulation of condensed phase systems. Most approaches used fixed, site-specific point charges. Yet, it is now clear that many applications of force fields (FFs) demand more sophisticated treatments, prompting the implementation of charge equilibration methods in polarizable FFs to allow the redistribution of charge within the system. One approach allowing both, charge redistribution and site-specific polarization, while at the same time solving methodological shortcomings of earlier methods, is the first-principles-derived atom-condensed Kohn-Sham density functional theory method approximated to the second order (ACKS2). In this work, we present two fragment approaches to ACKS2, termed f-ACKS2 and a self-consistent version, scf-ACKS2, that treat condensed phase systems as a collection of electronically polarizable molecular fragments. The fragmentation approach to ACKS2 not only leads to a more transferable and less system-specific collection of electronic response parameters but also opens up the method to large condensed phase systems. We validate the accuracies of f-ACKS2 and scf-ACKS2 by comparing polarization energies and induced dipole moments for a number of charged hydrocarbon dimers against DFT reference calculations. Finally, we also apply both fragmented ACKS2 variants to calculate the polarization energy for electron-hole pair separation along a chain of anthracene molecules and find excellent agreement with reference DFT calculations.
Collapse
Affiliation(s)
- Patrick Gütlein
- Chair for Theoretical Chemistry and Catalysis Research Center, Technische Universität München, Lichtenbergstrasse 4, D-85747 Garching, Germany
| | - Jochen Blumberger
- Department of Physics and Astronomy, University College London, London WC1E 6BT, U.K.,Institute for Advanced Study, Technische Universität München, Lichtenbergstrasse 2 a, D-85748 Garching, Germany
| | - Harald Oberhofer
- Chair for Theoretical Chemistry and Catalysis Research Center, Technische Universität München, Lichtenbergstrasse 4, D-85747 Garching, Germany
| |
Collapse
|
10
|
Chemometric Models of Differential Amino Acids at the Na vα and Na vβ Interface of Mammalian Sodium Channel Isoforms. Molecules 2020; 25:molecules25153551. [PMID: 32756517 PMCID: PMC7435598 DOI: 10.3390/molecules25153551] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/01/2020] [Revised: 07/21/2020] [Accepted: 07/22/2020] [Indexed: 12/19/2022] Open
Abstract
(1) Background: voltage-gated sodium channels (Navs) are integral membrane proteins that allow the sodium ion flux into the excitable cells and initiate the action potential. They comprise an α (Navα) subunit that forms the channel pore and are coupled to one or more auxiliary β (Navβ) subunits that modulate the gating to a variable extent. (2) Methods: after performing homology in silico modeling for all nine isoforms (Nav1.1α to Nav1.9α), the Navα and Navβ protein-protein interaction (PPI) was analyzed chemometrically based on the primary and secondary structures as well as topological or spatial mapping. (3) Results: our findings reveal a unique isoform-specific correspondence between certain segments of the extracellular loops of the Navα subunits. Precisely, loop S5 in domain I forms part of the PPI and assists Navβ1 or Navβ3 on all nine mammalian isoforms. The implied molecular movements resemble macroscopic springs, all of which explains published voltage sensor effects on sodium channel fast inactivation in gating. (4) Conclusions: currently, the specific functions exerted by the Navβ1 or Navβ3 subunits on the modulation of Navα gating remain unknown. Our work determined functional interaction in the extracellular domains on theoretical grounds and we propose a schematic model of the gating mechanism of fast channel sodium current inactivation by educated guessing.
Collapse
|
11
|
Marashiyan M, Kalhor H, Ganji M, Rahimi H. Effects of tosyl-l-arginine methyl ester (TAME) on the APC/c subunits: An in silico investigation for inhibiting cell cycle. J Mol Graph Model 2020; 97:107563. [PMID: 32066079 DOI: 10.1016/j.jmgm.2020.107563] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/27/2019] [Revised: 01/11/2020] [Accepted: 02/01/2020] [Indexed: 11/28/2022]
Abstract
The anaphase-promoting complex/cyclosome (APC/c) is requisite for controlling mitosis, which is activated by Cdh1 and Cdc20 activators. Dysregulation of APC/c is observed in many cancers and is known as a targeted drug particularly in cancer drug resistance. It was shown that tosyl-l-arginine methyl ester (TAME), via mimicking isoleucine-arginine (IR) tail of co-activators, inhibits APC/c functions. However, structure details and interaction of TAME with APC/c are poorly defined. In the current study, a well-established set of computational methods was used to identify the best binding pocket in order to inhibit APC activity. Therefore, the interaction of IR tail and Cbox of co-activators, as well as TAME as an inhibitor, as an inhibitor, with APC3 and APC8 subunits of APC/c were analyzed, regarding structure, molecular docking, molecular dynamics, and free binding energy. The results indicated that TAME bound to APC3 with a higher binding affinity (∼-7.3 kcal/mol) than APC8 (∼-5.7 kcal/mol). Also, the binding free energy value obtained for the APC3-TAME was -22.25 ± 1.12 kcal/mol. According to binding free energies, van der Waals energy was the major favorable contributor to the ligand binding. These results offer that TAME had more affinity to interact with the APC3 subunit, at the IR binding pocket than the APC8 subunit at the Cbox binding pocket. In conclusion, IR binding pocket can serve as an appropriate potential target for TAME as an inhibitor of APC/c.
Collapse
Affiliation(s)
- Mahya Marashiyan
- Molecular Medicine Department, Biotechnology Research Center, Pasteur Institute of Iran, Tehran, Iran
| | - Hourieh Kalhor
- Cellular and Molecular Research Center,Qom University of Medical Sciences, Qom, Iran
| | - Maziar Ganji
- Department of Medical Genetics, School of Medicine, Shahid Beheshti University of Medical Sciences, Tehran, Iran
| | - Hamzeh Rahimi
- Molecular Medicine Department, Biotechnology Research Center, Pasteur Institute of Iran, Tehran, Iran.
| |
Collapse
|
12
|
Gaillard T, Simonson T. Full Protein Sequence Redesign with an MMGBSA Energy Function. J Chem Theory Comput 2017; 13:4932-4943. [DOI: 10.1021/acs.jctc.7b00202] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Thomas Gaillard
- Laboratoire de Biochimie
(CNRS UMR7654), Department of Biology, Ecole Polytechnique, 91128 Palaiseau, France
| | - Thomas Simonson
- Laboratoire de Biochimie
(CNRS UMR7654), Department of Biology, Ecole Polytechnique, 91128 Palaiseau, France
| |
Collapse
|
13
|
Knowledge-based entropies improve the identification of native protein structures. Proc Natl Acad Sci U S A 2017; 114:2928-2933. [PMID: 28265078 DOI: 10.1073/pnas.1613331114] [Citation(s) in RCA: 20] [Impact Index Per Article: 2.9] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022] Open
Abstract
Evaluating protein structures requires reliable free energies with good estimates of both potential energies and entropies. Although there are many demonstrated successes from using knowledge-based potential energies, computing entropies of proteins has lagged far behind. Here we take an entirely different approach and evaluate knowledge-based conformational entropies of proteins based on the observed frequencies of contact changes between amino acids in a set of 167 diverse proteins, each of which has two alternative structures. The results show that charged and polar interactions break more often than hydrophobic pairs. This pattern correlates strongly with the average solvent exposure of amino acids in globular proteins, as well as with polarity indices and the sizes of the amino acids. Knowledge-based entropies are derived by using the inverse Boltzmann relationship, in a manner analogous to the way that knowledge-based potentials have been extracted. Including these new knowledge-based entropies almost doubles the performance of knowledge-based potentials in selecting the native protein structures from decoy sets. Beyond the overall energy-entropy compensation, a similar compensation is seen for individual pairs of interacting amino acids. The entropies in this report have immediate applications for 3D structure prediction, protein model assessment, and protein engineering and design.
Collapse
|
14
|
Abstract
Computational protein design (CPD) has established itself as a leading field in basic and applied science with a strong coupling between the two. Proteins are computationally designed from the level of amino acids to the level of a functional protein complex. Design targets range from increased thermo- (or other) stability to specific requested reactions such as protein-protein binding, enzymatic reactions, or nanotechnology applications. The design scheme may encompass small regions of the proteins or the entire protein. In either case, the design may aim at the side-chains or at the full backbone conformation. Herein, the main framework for the process is outlined highlighting key elements in the CPD iterative cycle. These include the very definition of CPD, the diverse goals of CPD, components of the CPD protocol, methods for searching sequence and structure space, scoring functions, and augmenting the CPD with other optimization tools. Taken together, this chapter aims to introduce the framework of CPD.
Collapse
Affiliation(s)
- Ilan Samish
- Department of Plants and Environmental Sciences, Weizmann Institute of Science, Rehovot, Israel.
- Department of Biotechnology Engineering, Braude Academic College of Engineering, Karmiel, Israel.
- Amai Proteins Ltd., Ashdod, Israel.
| |
Collapse
|
15
|
Nix MA, Kaelin CB, Palomino R, Miller JL, Barsh GS, Millhauser GL. Electrostatic Similarity Analysis of Human β-Defensin Binding in the Melanocortin System. Biophys J 2016; 109:1946-58. [PMID: 26536271 DOI: 10.1016/j.bpj.2015.09.005] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/13/2015] [Revised: 08/27/2015] [Accepted: 09/03/2015] [Indexed: 12/13/2022] Open
Abstract
The β-defensins are a class of small cationic proteins that serve as components of numerous systems in vertebrate biology, including the immune and melanocortin systems. Human β-defensin 3 (HBD3), which is produced in the skin, has been found to bind to melanocortin receptors 1 and 4 through complementary electrostatics, a unique mechanism of ligand-receptor interaction. This finding indicates that electrostatics alone, and not specific amino acid contact points, could be sufficient for function in this ligand-receptor system, and further suggests that other small peptide ligands could interact with these receptors in a similar fashion. Here, we conducted molecular-similarity analyses and functional studies of additional members of the human β-defensin family, examining their potential as ligands of melanocortin-1 receptor, through selection based on their electrostatic similarity to HBD3. Using Poisson-Boltzmann electrostatic calculations and molecular-similarity analysis, we identified members of the human β-defensin family that are both similar and dissimilar to HBD3 in terms of electrostatic potential. Synthesis and functional testing of a subset of these β-defensins showed that peptides with an HBD3-like electrostatic character bound to melanocortin receptors with high affinity, whereas those that were anticorrelated to HBD3 showed no binding affinity. These findings expand on the central role of electrostatics in the control of this ligand-receptor system and further demonstrate the utility of employing molecular-similarity analysis. Additionally, we identified several new potential ligands of melanocortin-1 receptor, which may have implications for our understanding of the role defensins play in melanocortin physiology.
Collapse
Affiliation(s)
- Matthew A Nix
- Department of Chemistry and Biochemistry, University of California, Santa Cruz, Santa Cruz, California
| | - Christopher B Kaelin
- HudsonAlpha Institute for Biotechnology, Huntsville, Alabama; Department of Genetics, Stanford University, Stanford, California
| | - Rafael Palomino
- Department of Chemistry and Biochemistry, University of California, Santa Cruz, Santa Cruz, California
| | - Jillian L Miller
- Department of Chemistry and Biochemistry, University of California, Santa Cruz, Santa Cruz, California
| | - Gregory S Barsh
- HudsonAlpha Institute for Biotechnology, Huntsville, Alabama; Department of Genetics, Stanford University, Stanford, California.
| | - Glenn L Millhauser
- Department of Chemistry and Biochemistry, University of California, Santa Cruz, Santa Cruz, California.
| |
Collapse
|
16
|
Parmar AS, James JK, Grisham DR, Pike DH, Nanda V. Dissecting Electrostatic Contributions to Folding and Self-Assembly Using Designed Multicomponent Peptide Systems. J Am Chem Soc 2016; 138:4362-7. [PMID: 26966815 DOI: 10.1021/jacs.5b10304] [Citation(s) in RCA: 25] [Impact Index Per Article: 3.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
We investigate formation of nano- to microscale peptide fibers and sheets where assembly requires association of two distinct collagen mimetic peptides (CMPs). The multicomponent nature of these designs allows the decoupling of amino acid contributions to peptide folding versus higher-order assembly. While both arginine and lysine containing CMP sequences can favor triple-helix folding, only arginine promotes rapid supramolecular assembly in each of the three two-component systems examined. Unlike lysine, the polyvalent guanidyl group of arginine is capable of both intra- and intermolecular contacts, promoting assembly. This is consistent with the supramolecular diversity of CMP morphologies observed throughout the literature. It also connects CMP self-assembly with a broad range of biomolecular interaction phenomena, providing general principles for modeling and design.
Collapse
Affiliation(s)
- Avanish S Parmar
- Department of Physics, Indian Institute of Technology (Banaras Hindu University) , Varanasi 221005, Uttar Pradesh, India
| | - Jose K James
- Center for Advanced Biotechnology and Medicine, Department of Biochemistry and Molecular Biology, Robert Wood Johnson Medical School, Rutgers University , 679 Hoes Lane West, Piscataway, New Jersey 08854, United States
| | - Daniel R Grisham
- Center for Advanced Biotechnology and Medicine, Department of Biochemistry and Molecular Biology, Robert Wood Johnson Medical School, Rutgers University , 679 Hoes Lane West, Piscataway, New Jersey 08854, United States
| | - Douglas H Pike
- Center for Advanced Biotechnology and Medicine, Department of Biochemistry and Molecular Biology, Robert Wood Johnson Medical School, Rutgers University , 679 Hoes Lane West, Piscataway, New Jersey 08854, United States
| | - Vikas Nanda
- Center for Advanced Biotechnology and Medicine, Department of Biochemistry and Molecular Biology, Robert Wood Johnson Medical School, Rutgers University , 679 Hoes Lane West, Piscataway, New Jersey 08854, United States
| |
Collapse
|
17
|
Estrada J, Echenique P, Sancho J. Predicting stabilizing mutations in proteins using Poisson-Boltzmann based models: study of unfolded state ensemble models and development of a successful binary classifier based on residue interaction energies. Phys Chem Chem Phys 2015; 17:31044-54. [PMID: 26530878 DOI: 10.1039/c5cp04348d] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/02/2023]
Abstract
In many cases the stability of a protein has to be increased to permit its biotechnological use. Rational methods of protein stabilization based on optimizing electrostatic interactions have provided some fine successful predictions. However, the precise calculation of stabilization energies remains challenging, one reason being that the electrostatic effects on the unfolded state are often neglected. We have explored here the feasibility of incorporating Poisson-Boltzmann model electrostatic calculations performed on representations of the unfolded state as large ensembles of geometrically optimized conformations calculated using the ProtSA server. Using a data set of 80 electrostatic mutations experimentally tested in two-state proteins, the predictive performance of several such models has been compared to that of a simple one that considers an unfolded structure of non-interacting residues. The unfolded ensemble models, while showing correlation between the predicted stabilization values and the experimental ones, are worse than the simple model, suggesting that the ensembles do not capture well the energetics of the unfolded state. A more attainable goal is classifying potential mutations as either stabilizing or non-stabilizing, rather than accurately calculating their stabilization energies. To implement a fast classification method that can assist in selecting stabilizing mutations, we have used a much simpler electrostatic model based only on the native structure and have determined its precision using different stabilizing energy thresholds. The binary classifier developed finds 7 true stabilizing mutants out of every 10 proposed candidates and can be used as a robust tool to propose stabilizing mutations.
Collapse
Affiliation(s)
- Jorge Estrada
- Departamento de Bioquímica y Biología Molecular y Celular, Facultad de Ciencias, Universidad de Zaragoza, Pedro Cerbuna 12, 50009 Zaragoza, Spain. and Biocomputation and Complex Systems Physics Institute (BIFI), Joint Unit BIFI-IQFR (CSIC), Mariano Esquillor s/n, Edificio I+D, 50018, Zaragoza, Spain
| | - Pablo Echenique
- Biocomputation and Complex Systems Physics Institute (BIFI), Joint Unit BIFI-IQFR (CSIC), Mariano Esquillor s/n, Edificio I+D, 50018, Zaragoza, Spain and Instituto de Química Física "Rocasolano", CSIC, Serrano 119, 28006, Madrid, Spain
| | - Javier Sancho
- Departamento de Bioquímica y Biología Molecular y Celular, Facultad de Ciencias, Universidad de Zaragoza, Pedro Cerbuna 12, 50009 Zaragoza, Spain. and Biocomputation and Complex Systems Physics Institute (BIFI), Joint Unit BIFI-IQFR (CSIC), Mariano Esquillor s/n, Edificio I+D, 50018, Zaragoza, Spain
| |
Collapse
|
18
|
Druart K, Palmai Z, Omarjee E, Simonson T. Protein:Ligand binding free energies: A stringent test for computational protein design. J Comput Chem 2015; 37:404-15. [PMID: 26503829 DOI: 10.1002/jcc.24230] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/01/2015] [Revised: 10/01/2015] [Accepted: 10/02/2015] [Indexed: 01/29/2023]
Abstract
A computational protein design method is extended to allow Monte Carlo simulations where two ligands are titrated into a protein binding pocket, yielding binding free energy differences. These provide a stringent test of the physical model, including the energy surface and sidechain rotamer definition. As a test, we consider tyrosyl-tRNA synthetase (TyrRS), which has been extensively redesigned experimentally. We consider its specificity for its substrate l-tyrosine (l-Tyr), compared to the analogs d-Tyr, p-acetyl-, and p-azido-phenylalanine (ac-Phe, az-Phe). We simulate l- and d-Tyr binding to TyrRS and six mutants, and compare the structures and binding free energies to a more rigorous "MD/GBSA" procedure: molecular dynamics with explicit solvent for structures and a Generalized Born + Surface Area model for binding free energies. Next, we consider l-Tyr, ac- and az-Phe binding to six other TyrRS variants. The titration results are sensitive to the precise rotamer definition, which involves a short energy minimization for each sidechain pair to help relax bad contacts induced by the discrete rotamer set. However, when designed mutant structures are rescored with a standard GBSA energy model, results agree well with the more rigorous MD/GBSA. As a third test, we redesign three amino acid positions in the substrate coordination sphere, with either l-Tyr or d-Tyr as the ligand. For two, we obtain good agreement with experiment, recovering the wildtype residue when l-Tyr is the ligand and a d-Tyr specific mutant when d-Tyr is the ligand. For the third, we recover His with either ligand, instead of wildtype Gln.
Collapse
Affiliation(s)
- Karen Druart
- Laboratoire De Biochimie (UMR CNRS 7654), Department of Biology, Ecole Polytechnique, Palaiseau, France
| | - Zoltan Palmai
- Laboratoire De Biochimie (UMR CNRS 7654), Department of Biology, Ecole Polytechnique, Palaiseau, France
| | - Eyaz Omarjee
- Laboratoire De Biochimie (UMR CNRS 7654), Department of Biology, Ecole Polytechnique, Palaiseau, France
| | - Thomas Simonson
- Laboratoire De Biochimie (UMR CNRS 7654), Department of Biology, Ecole Polytechnique, Palaiseau, France
| |
Collapse
|
19
|
Computational tools for epitope vaccine design and evaluation. Curr Opin Virol 2015; 11:103-12. [PMID: 25837467 DOI: 10.1016/j.coviro.2015.03.013] [Citation(s) in RCA: 47] [Impact Index Per Article: 5.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/26/2014] [Revised: 03/13/2015] [Accepted: 03/16/2015] [Indexed: 12/15/2022]
Abstract
Rational approaches will be required to develop universal vaccines for viral pathogens such as human immunodeficiency virus, hepatitis C virus, and influenza, for which empirical approaches have failed. The main objective of a rational vaccine strategy is to design novel immunogens that are capable of inducing long-term protective immunity. In practice, this requires structure-based engineering of the target neutralizing epitopes and a quantitative readout of vaccine-induced immune responses. Therefore, computational tools that can facilitate these two areas have played increasingly important roles in rational vaccine design in recent years. Here we review the computational techniques developed for protein structure prediction and antibody repertoire analysis, and demonstrate how they can be applied to the design and evaluation of epitope vaccines.
Collapse
|
20
|
Koehl P, Poitevin F, Orland H, Delarue M. Modified Poisson–Boltzmann equations for characterizing biomolecular solvation. JOURNAL OF THEORETICAL & COMPUTATIONAL CHEMISTRY 2014. [DOI: 10.1142/s021963361440001x] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/18/2022]
Abstract
Methods for computing electrostatic interactions often account implicitly for the solvent, due to the much smaller number of degrees of freedom involved. In the Poisson–Boltzmann (PB) approach the electrostatic potential is obtained by solving the Poisson–Boltzmann equation (PBE), where the solvent region is modeled as a homogeneous medium with a high dielectric constant. PB however is not exempt of problems. It does not take into account for example the sizes of the ions in the atmosphere surrounding the solute, nor does it take into account the inhomogeneous dielectric response of water due to the presence of a highly charged surface. In this paper we review two major modifications of PB that circumvent these problems, namely the size-modified PB (SMPB) equation and the Dipolar Poisson–Boltzmann Langevin (DPBL) model. In SMPB, steric effects between ions are accounted for with a lattice gas model. In DPBL, the solvent region is no longer modeled as a homogeneous dielectric media but rather as an assembly of self-orienting interacting dipoles of variable density. This model results in a dielectric profile that transits smoothly from the solute to the solvent region as well as in a variable solvent density that depends on the charges of the solute. We show successful applications of the DPBL formalism to computing the solvation free energies of isolated ions in water. Further developments of more accurately modified PB models are discussed.
Collapse
Affiliation(s)
- Patrice Koehl
- Department of Computer Science and Genome Center, University of California, Davis, CA 95616, USA
| | - Frederic Poitevin
- Unité de Dynamique Structurale des Macromolécules, UMR 3528 du CNRS, Institut Pasteur, 75015 Paris, France
| | - Henri Orland
- Service de Physique Théorique, CEA-Saclay, 91191 Gif/Yvette Cedex, France
| | - Marc Delarue
- Unité de Dynamique Structurale des Macromolécules, UMR 3528 du CNRS, Institut Pasteur, 75015 Paris, France
| |
Collapse
|
21
|
An accurate binding interaction model in de novo computational protein design of interactions: If you build it, they will bind. J Struct Biol 2014; 185:136-46. [DOI: 10.1016/j.jsb.2013.03.012] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/22/2012] [Revised: 03/15/2013] [Accepted: 03/21/2013] [Indexed: 01/07/2023]
|
22
|
Schwans JP, Sunden F, Gonzalez A, Tsai Y, Herschlag D. Uncovering the determinants of a highly perturbed tyrosine pKa in the active site of ketosteroid isomerase. Biochemistry 2013; 52:7840-55. [PMID: 24151972 PMCID: PMC3890242 DOI: 10.1021/bi401083b] [Citation(s) in RCA: 35] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/03/2023]
Abstract
Within the idiosyncratic enzyme active-site environment, side chain and ligand pKa values can be profoundly perturbed relative to their values in aqueous solution. Whereas structural inspection of systems has often attributed perturbed pKa values to dominant contributions from placement near charged groups or within hydrophobic pockets, Tyr57 of a Pseudomonas putida ketosteroid isomerase (KSI) mutant, suggested to have a pKa perturbed by nearly 4 units to 6.3, is situated within a solvent-exposed active site devoid of cationic side chains, metal ions, or cofactors. Extensive comparisons among 45 variants with mutations in and around the KSI active site, along with protein semisynthesis, (13)C NMR spectroscopy, absorbance spectroscopy, and X-ray crystallography, was used to unravel the basis for this perturbed Tyr pKa. The results suggest that the origin of large energetic perturbations are more complex than suggested by visual inspection. For example, the introduction of positively charged residues near Tyr57 raises its pKa rather than lowers it; this effect, and part of the increase in the Tyr pKa from the introduction of nearby anionic groups, arises from accompanying active-site structural rearrangements. Other mutations with large effects also cause structural perturbations or appear to displace a structured water molecule that is part of a stabilizing hydrogen-bond network. Our results lead to a model in which three hydrogen bonds are donated to the stabilized ionized Tyr, with these hydrogen-bond donors, two Tyr side chains, and a water molecule positioned by other side chains and by a water-mediated hydrogen-bond network. These results support the notion that large energetic effects are often the consequence of multiple stabilizing interactions rather than a single dominant interaction. Most generally, this work provides a case study for how extensive and comprehensive comparisons via site-directed mutagenesis in a tight feedback loop with structural analysis can greatly facilitate our understanding of enzyme active-site energetics. The extensive data set provided may also be a valuable resource for those wishing to extensively test computational approaches for determining enzymatic pKa values and energetic effects.
Collapse
Affiliation(s)
- Jason P. Schwans
- Department of Biochemistry, Stanford University, Stanford, California 94305
| | - Fanny Sunden
- Department of Biochemistry, Stanford University, Stanford, California 94305
| | - Ana Gonzalez
- Stanford Synchrotron Radiation Lightsource, SLAC National Accelerator Laboratory, Menlo Park, California 94025
| | - Yingssu Tsai
- Department of Chemistry, Stanford University, Stanford, California 94305
- Stanford Synchrotron Radiation Lightsource, SLAC National Accelerator Laboratory, Menlo Park, California 94025
| | - Daniel Herschlag
- Department of Biochemistry, Stanford University, Stanford, California 94305
- Department of Chemistry, Stanford University, Stanford, California 94305
| |
Collapse
|
23
|
Simonson T. What Is the Dielectric Constant of a Protein When Its Backbone Is Fixed? J Chem Theory Comput 2013; 9:4603-8. [DOI: 10.1021/ct400398e] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Thomas Simonson
- Laboratoire de Biochimie
(CNRS UMR7654), Department of Biology, Ecole Polytechnique, 91128 Palaiseau, France
| |
Collapse
|
24
|
Muddana HS, Sapra NV, Fenley AT, Gilson MK. The electrostatic response of water to neutral polar solutes: implications for continuum solvent modeling. J Chem Phys 2013; 138:224504. [PMID: 23781802 PMCID: PMC3695974 DOI: 10.1063/1.4808376] [Citation(s) in RCA: 18] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/13/2013] [Accepted: 05/10/2013] [Indexed: 12/20/2022] Open
Abstract
Continuum solvation models are widely used to estimate the hydration free energies of small molecules and proteins, in applications ranging from drug design to protein engineering, and most such models are based on the approximation of a linear dielectric response by the solvent. We used explicit-water molecular dynamics simulations with the TIP3P water model to probe this linear response approximation in the case of neutral polar molecules, using miniature cucurbituril and cyclodextrin receptors and protein side-chain analogs as model systems. We observe supralinear electrostatic solvent responses, and this nonlinearity is found to result primarily from waters' being drawn closer and closer to the solutes with increased solute-solvent electrostatic interactions; i.e., from solute electrostriction. Dielectric saturation and changes in the water-water hydrogen bonding network, on the other hand, play little role. Thus, accounting for solute electrostriction may be a productive approach to improving the accuracy of continuum solvation models.
Collapse
Affiliation(s)
- Hari S Muddana
- Skaggs School of Pharmacy and Pharmaceutical Sciences, University of California San Diego, La Jolla, California 92093-0736, USA
| | | | | | | |
Collapse
|
25
|
Wijma HJ, Janssen DB. Computational design gains momentum in enzyme catalysis engineering. FEBS J 2013; 280:2948-60. [DOI: 10.1111/febs.12324] [Citation(s) in RCA: 53] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/20/2013] [Revised: 04/19/2013] [Accepted: 04/24/2013] [Indexed: 01/19/2023]
Affiliation(s)
- Hein J. Wijma
- Department of Biochemistry; Groningen Biomolecular Sciences and Biotechnology Institute; University of Groningen; The Netherlands
| | - Dick B. Janssen
- Department of Biochemistry; Groningen Biomolecular Sciences and Biotechnology Institute; University of Groningen; The Netherlands
| |
Collapse
|
26
|
Li Z, Yang Y, Zhan J, Dai L, Zhou Y. Energy functions in de novo protein design: current challenges and future prospects. Annu Rev Biophys 2013; 42:315-35. [PMID: 23451890 DOI: 10.1146/annurev-biophys-083012-130315] [Citation(s) in RCA: 65] [Impact Index Per Article: 5.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/26/2022]
Abstract
In the past decade, a concerted effort to successfully capture specific tertiary packing interactions produced specific three-dimensional structures for many de novo designed proteins that are validated by nuclear magnetic resonance and/or X-ray crystallographic techniques. However, the success rate of computational design remains low. In this review, we provide an overview of experimentally validated, de novo designed proteins and compare four available programs, RosettaDesign, EGAD, Liang-Grishin, and RosettaDesign-SR, by assessing designed sequences computationally. Computational assessment includes the recovery of native sequences, the calculation of sizes of hydrophobic patches and total solvent-accessible surface area, and the prediction of structural properties such as intrinsic disorder, secondary structures, and three-dimensional structures. This computational assessment, together with a recent community-wide experiment in assessing scoring functions for interface design, suggests that the next-generation protein-design scoring function will come from the right balance of complementary interaction terms. Such balance may be found when more negative experimental data become available as part of a training set.
Collapse
Affiliation(s)
- Zhixiu Li
- School of Informatics, Indiana University-Purdue University, Indianapolis, Indiana 46202, USA
| | | | | | | | | |
Collapse
|
27
|
Designing electrostatic interactions in biological systems via charge optimization or combinatorial approaches: insights and challenges with a continuum electrostatic framework. Theor Chem Acc 2012. [DOI: 10.1007/s00214-012-1252-5] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
|
28
|
Yu CM, Peng HP, Chen IC, Lee YC, Chen JB, Tsai KC, Chen CT, Chang JY, Yang EW, Hsu PC, Jian JW, Hsu HJ, Chang HJ, Hsu WL, Huang KF, Ma AC, Yang AS. Rationalization and design of the complementarity determining region sequences in an antibody-antigen recognition interface. PLoS One 2012; 7:e33340. [PMID: 22457753 PMCID: PMC3310866 DOI: 10.1371/journal.pone.0033340] [Citation(s) in RCA: 37] [Impact Index Per Article: 3.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/02/2011] [Accepted: 02/14/2012] [Indexed: 12/01/2022] Open
Abstract
Protein-protein interactions are critical determinants in biological systems. Engineered proteins binding to specific areas on protein surfaces could lead to therapeutics or diagnostics for treating diseases in humans. But designing epitope-specific protein-protein interactions with computational atomistic interaction free energy remains a difficult challenge. Here we show that, with the antibody-VEGF (vascular endothelial growth factor) interaction as a model system, the experimentally observed amino acid preferences in the antibody-antigen interface can be rationalized with 3-dimensional distributions of interacting atoms derived from the database of protein structures. Machine learning models established on the rationalization can be generalized to design amino acid preferences in antibody-antigen interfaces, for which the experimental validations are tractable with current high throughput synthetic antibody display technologies. Leave-one-out cross validation on the benchmark system yielded the accuracy, precision, recall (sensitivity) and specificity of the overall binary predictions to be 0.69, 0.45, 0.63, and 0.71 respectively, and the overall Matthews correlation coefficient of the 20 amino acid types in the 24 interface CDR positions was 0.312. The structure-based computational antibody design methodology was further tested with other antibodies binding to VEGF. The results indicate that the methodology could provide alternatives to the current antibody technologies based on animal immune systems in engineering therapeutic and diagnostic antibodies against predetermined antigen epitopes.
Collapse
Affiliation(s)
- Chung-Ming Yu
- Genomics Research Center, Academia Sinica, Taipei, Taiwan
| | - Hung-Pin Peng
- Genomics Research Center, Academia Sinica, Taipei, Taiwan
- Institute of Biomedical Informatics, National Yang-Ming University, Taipei, Taiwan
- Bioinformatics Program, Taiwan International Graduate Program, Institute of Information Science, Academia Sinica, Taipei, Taiwan
| | - Ing-Chien Chen
- Genomics Research Center, Academia Sinica, Taipei, Taiwan
| | - Yu-Ching Lee
- Genomics Research Center, Academia Sinica, Taipei, Taiwan
| | - Jun-Bo Chen
- Genomics Research Center, Academia Sinica, Taipei, Taiwan
- Department of Computer Science, National Tsing-Hua University, Hsinchu, Taiwan
| | | | - Ching-Tai Chen
- Genomics Research Center, Academia Sinica, Taipei, Taiwan
- Institute of Bioinformatics and Systems Biology, National Chiao-Tung University, Hsinchu, Taiwan
- Institute of Information Sciences, Academia Sinica, Taipei, Taiwan
- Bioinformatics Program, Taiwan International Graduate Program, Institute of Information Science, Academia Sinica, Taipei, Taiwan
| | - Jeng-Yih Chang
- Genomics Research Center, Academia Sinica, Taipei, Taiwan
| | - Ei-Wen Yang
- Genomics Research Center, Academia Sinica, Taipei, Taiwan
- Institute of Information Sciences, Academia Sinica, Taipei, Taiwan
| | - Po-Chiang Hsu
- Genomics Research Center, Academia Sinica, Taipei, Taiwan
| | - Jhih-Wei Jian
- Genomics Research Center, Academia Sinica, Taipei, Taiwan
- Institute of Biomedical Informatics, National Yang-Ming University, Taipei, Taiwan
- Bioinformatics Program, Taiwan International Graduate Program, Institute of Information Science, Academia Sinica, Taipei, Taiwan
| | - Hung-Ju Hsu
- Genomics Research Center, Academia Sinica, Taipei, Taiwan
- Graduate Institute of Life Sciences, National Defense University, Taipei, Taiwan
| | - Hung-Ju Chang
- Genomics Research Center, Academia Sinica, Taipei, Taiwan
- Institute of Biochemical Science, National Taiwan University, Taipei, Taiwan
- Chemical Biology and Molecular Biophysics Program, Taiwan International Graduate Program, Institute of Biological Chemistry, Academia Sinica, Taipei, Taiwan
| | - Wen-Lian Hsu
- Institute of Information Sciences, Academia Sinica, Taipei, Taiwan
| | - Kai-Fa Huang
- Institute of Biological Chemistry, Academia Sinica, Taipei, Taiwan
| | - Alex Che Ma
- Genomics Research Center, Academia Sinica, Taipei, Taiwan
| | - An-Suei Yang
- Genomics Research Center, Academia Sinica, Taipei, Taiwan
| |
Collapse
|
29
|
Verstraelen T, Pauwels E, De Proft F, Van Speybroeck V, Geerlings P, Waroquier M. Assessment of Atomic Charge Models for Gas-Phase Computations on Polypeptides. J Chem Theory Comput 2012; 8:661-76. [PMID: 26596614 DOI: 10.1021/ct200512e] [Citation(s) in RCA: 57] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/21/2023]
Abstract
The concept of the atomic charge is extensively used to model the electrostatic properties of proteins. Atomic charges are not only the basis for the electrostatic energy term in biomolecular force fields but are also derived from quantum mechanical computations on protein fragments to get more insight into their electronic structure. Unfortunately there are many atomic charge schemes which lead to significantly different results, and it is not trivial to determine which scheme is most suitable for biomolecular studies. Therefore, we present an extensive methodological benchmark using a selection of atomic charge schemes [Mulliken, natural, restrained electrostatic potential, Hirshfeld-I, electronegativity equalization method (EEM), and split-charge equilibration (SQE)] applied to two sets of penta-alanine conformers. Our analysis clearly shows that Hirshfeld-I charges offer the best compromise between transferability (robustness with respect to conformational changes) and the ability to reproduce electrostatic properties of the penta-alanine. The benchmark also considers two charge equilibration models (EEM and SQE), which both clearly fail to describe the locally charged moieties in the zwitterionic form of penta-alanine. This issue is analyzed in detail because charge equilibration models are computationally much more attractive than the Hirshfeld-I scheme. Based on the latter analysis, a straightforward extension of the SQE model is proposed, SQE+Q(0), that is suitable to describe biological systems bearing many locally charged functional groups.
Collapse
Affiliation(s)
- Toon Verstraelen
- Center For Molecular Modeling, Ghent University , Technologiepark 903, 9050 Zwijnaarde, Belgium (Member of the QCMM Ghent-Brussels Alliance Group)
| | - Ewald Pauwels
- Center For Molecular Modeling, Ghent University , Technologiepark 903, 9050 Zwijnaarde, Belgium (Member of the QCMM Ghent-Brussels Alliance Group)
| | - Frank De Proft
- Department of General Chemistry (ALGC), Free University of Brussels-VUB , Pleinlaan 2, 1050 Brussels, Belgium (Member of the QCMM Ghent-Brussels Alliance Group)
| | - Veronique Van Speybroeck
- Center For Molecular Modeling, Ghent University , Technologiepark 903, 9050 Zwijnaarde, Belgium (Member of the QCMM Ghent-Brussels Alliance Group)
| | - Paul Geerlings
- Department of General Chemistry (ALGC), Free University of Brussels-VUB , Pleinlaan 2, 1050 Brussels, Belgium (Member of the QCMM Ghent-Brussels Alliance Group)
| | - Michel Waroquier
- Center For Molecular Modeling, Ghent University , Technologiepark 903, 9050 Zwijnaarde, Belgium (Member of the QCMM Ghent-Brussels Alliance Group)
| |
Collapse
|
30
|
Koehl P, Orland H, Delarue M. Adapting Poisson-Boltzmann to the self-consistent mean field theory: application to protein side-chain modeling. J Chem Phys 2011; 135:055104. [PMID: 21823735 DOI: 10.1063/1.3621831] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We present an extension of the self-consistent mean field theory for protein side-chain modeling in which solvation effects are included based on the Poisson-Boltzmann (PB) theory. In this approach, the protein is represented with multiple copies of its side chains. Each copy is assigned a weight that is refined iteratively based on the mean field energy generated by the rest of the protein, until self-consistency is reached. At each cycle, the variational free energy of the multi-copy system is computed; this free energy includes the internal energy of the protein that accounts for vdW and electrostatics interactions and a solvation free energy term that is computed using the PB equation. The method converges in only a few cycles and takes only minutes of central processing unit time on a commodity personal computer. The predicted conformation of each residue is then set to be its copy with the highest weight after convergence. We have tested this method on a database of hundred highly refined NMR structures to circumvent the problems of crystal packing inherent to x-ray structures. The use of the PB-derived solvation free energy significantly improves prediction accuracy for surface side chains. For example, the prediction accuracies for χ(1) for surface cysteine, serine, and threonine residues improve from 68%, 35%, and 43% to 80%, 53%, and 57%, respectively. A comparison with other side-chain prediction algorithms demonstrates that our approach is consistently better in predicting the conformations of exposed side chains.
Collapse
Affiliation(s)
- Patrice Koehl
- Department of Biological Sciences, National University of Singapore, Singapore.
| | | | | |
Collapse
|
31
|
Bardhan JP, Knepley MG. Mathematical analysis of the boundary-integral based electrostatics estimation approximation for molecular solvation: exact results for spherical inclusions. J Chem Phys 2011; 135:124107. [PMID: 21974512 DOI: 10.1063/1.3641485] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We analyze the mathematically rigorous BIBEE (boundary-integral based electrostatics estimation) approximation of the mixed-dielectric continuum model of molecular electrostatics, using the analytically solvable case of a spherical solute containing an arbitrary charge distribution. Our analysis, which builds on Kirkwood's solution using spherical harmonics, clarifies important aspects of the approximation and its relationship to generalized Born models. First, our results suggest a new perspective for analyzing fast electrostatic models: the separation of variables between material properties (the dielectric constants) and geometry (the solute dielectric boundary and charge distribution). Second, we find that the eigenfunctions of the reaction-potential operator are exactly preserved in the BIBEE model for the sphere, which supports the use of this approximation for analyzing charge-charge interactions in molecular binding. Third, a comparison of BIBEE to the recent GBε theory suggests a modified BIBEE model capable of predicting electrostatic solvation free energies to within 4% of a full numerical Poisson calculation. This modified model leads to a projection-framework understanding of BIBEE and suggests opportunities for future improvements.
Collapse
Affiliation(s)
- Jaydeep P Bardhan
- Department of Molecular Biophysics and Physiology, Rush University Medical Center, Chicago, Illinois 60612, USA.
| | | |
Collapse
|
32
|
Leherte L, Vercauteren DP. Charge density distributions derived from smoothed electrostatic potential functions: design of protein reduced point charge models. J Comput Aided Mol Des 2011; 25:913-30. [DOI: 10.1007/s10822-011-9471-8] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/19/2011] [Accepted: 08/30/2011] [Indexed: 11/24/2022]
|
33
|
Abstract
Accurate computational methods for predicting electrostatic energies are of major importance for our understanding of protein energetics in general for computer-aided drug design as well as for the design of novel biocatalysts and protein therapeutics. Electrostatic energies are of particular importance in such applications as virtual screening, drug design and protein-protein docking due to the high charge density of protein ligands and small-molecule drugs, and the frequent protonation state changes observed when drugs bind to their protein targets. Therefore, the development of a reliable and fast algorithm for the evaluation of electrostatic free energies, as an important contributor to the overall protein energy function, has been the focus for many scientists over the past three decades. In this review we describe the current state-of-the-art in modeling electrostatic effects in proteins and protein-ligand complexes. We focus mainly on the merits and drawbacks of the continuum methodology, and speculate on future directions in refining algorithms for calculating electrostatic energies in proteins using experimental data.
Collapse
|
34
|
Electrostatics of folded and unfolded bovine β-lactoglobulin. Amino Acids 2011; 42:2019-30. [DOI: 10.1007/s00726-011-0933-z] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/31/2011] [Accepted: 04/30/2011] [Indexed: 11/26/2022]
|
35
|
Polydorides S, Amara N, Aubard C, Plateau P, Simonson T, Archontis G. Computational protein design with a generalized Born solvent model: application to Asparaginyl-tRNA synthetase. Proteins 2011; 79:3448-68. [PMID: 21563215 DOI: 10.1002/prot.23042] [Citation(s) in RCA: 19] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/29/2010] [Revised: 02/25/2011] [Accepted: 03/03/2011] [Indexed: 12/13/2022]
Abstract
Computational Protein Design (CPD) is a promising method for high throughput protein and ligand mutagenesis. Recently, we developed a CPD method that used a polar-hydrogen energy function for protein interactions and a Coulomb/Accessible Surface Area (CASA) model for solvent effects. We applied this method to engineer aspartyl-adenylate (AspAMP) specificity into Asparaginyl-tRNA synthetase (AsnRS), whose substrate is asparaginyl-adenylate (AsnAMP). Here, we implement a more accurate function, with an all-atom energy for protein interactions and a residue-pairwise generalized Born model for solvent effects. As a first test, we compute aminoacid affinities for several point mutants of Aspartyl-tRNA synthetase (AspRS) and Tyrosyl-tRNA synthetase and stability changes for three helical peptides and compare with experiment. As a second test, we readdress the problem of AsnRS aminoacid engineering. We compare three design criteria, which optimize the folding free-energy, the absolute AspAMP affinity, and the relative (AspAMP-AsnAMP) affinity. The sequences and conformations are improved with respect to our previous, polar-hydrogen/CASA study: For several designed complexes, the AspAMP carboxylate forms three interactions with a conserved arginine and a designed lysine, as in the active site of the AspRS:AspAMP complex. The conformations and interactions are well maintained in molecular dynamics simulations and the sequences have an inverted specificity, favoring AspAMP over AsnAMP. The method is not fully successful, since experimental measurements with the seven most promising sequences show that they do not catalyze at a detectable level the adenylation of Asp (or Asn) with ATP. This may be due to weak AspAMP binding and/or disruption of transition-state stabilization.
Collapse
|
36
|
Samish I, MacDermaid CM, Perez-Aguilar JM, Saven JG. Theoretical and Computational Protein Design. Annu Rev Phys Chem 2011; 62:129-49. [DOI: 10.1146/annurev-physchem-032210-103509] [Citation(s) in RCA: 119] [Impact Index Per Article: 9.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/09/2022]
Affiliation(s)
| | | | | | - Jeffery G. Saven
- Department of Chemistry, University of Pennsylvania, Philadelphia, Pennsylvania 19104;
| |
Collapse
|
37
|
Schmid FX. Lessons about Protein Stability from in vitro Selections. Chembiochem 2011; 12:1501-7. [DOI: 10.1002/cbic.201100018] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/11/2011] [Indexed: 11/07/2022]
|
38
|
Zarić BL, Jovanović VB, Stojanović SĐ. Non-covalent interactions across subunit interfaces in Sm proteins. J Theor Biol 2011; 271:18-26. [DOI: 10.1016/j.jtbi.2010.11.025] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/15/2010] [Revised: 10/11/2010] [Accepted: 11/18/2010] [Indexed: 11/29/2022]
|
39
|
Yang CM, Zhang J. Insights into intramolecular Trp and His side-chain orientation and stereospecific π interactions surrounding metal centers: an investigation using protein metal-site mimicry in solution. Chemistry 2011; 16:10854-65. [PMID: 20669189 DOI: 10.1002/chem.200903149] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/10/2022]
Abstract
Metal-binding scaffolds incorporating a Trp/His-paired epitope are instrumental in giving novel insights into the physicochemical basis of functional and mechanistic versatility conferred by the Trp-His interplay at a metal site. Herein, by coupling biometal site mimicry and (1)H and (13)C NMR spectroscopy experiments, modular constructs EDTA-(L-Trp, L-His) (EWH; EDTA=ethylenediamino tetraacetic acid) and DTPA-(L-Trp, L-His) (DWH; DTPA=diethylenetriamino pentaacetic acid) were employed to dissect the static and transient physicochemical properties of hydrophobic/hydrophilic aromatic interactive modes surrounding biometal centers. The binding feature and identities of the stoichiometric metal-bound complexes in solution were investigated by using (1)H and (13)C NMR spectroscopy, which facilitated a cross-validation of the carboxylate, amide oxygen, and tertiary amino groups as the primary ligands and indole as the secondary ligand, with the imidazole (Im) N3 nitrogen being weakly bound to metals such as Ca(2+) owing to a multivalency effect. Surrounding the metal centers, the stereospecific orientation of aromatic rings in the diastereoisomerism is interpreted with the Ca(2+)-EWH complex. With respect to perturbed Trp side-chain rotamer heterogeneity, drastically restricted Trp side-chain flexibility and thus a dynamically constrained rotamer interconversion due to π interactions is evident from the site-selective (13)C NMR spectroscopic signal broadening of the Trp indolyl C3 atom. Furthermore, effects of Trp side-chain fluctuation on indole/Im orientation were the subject of a 2D NMR spectroscopy study by using the Ca(2+)-bound state; a C-H2(indolyl)/C-H5(Im(+)) connectivity observed in the NOESY spectra captured direct evidence that the N-H1 of the Ca(2+)-Im(+) unit interacted with the pyrrole ring of the indole unit in Ca(2+)-bound EWH but not in DWH, which is assignable to a moderately static, anomalous, T-shaped, interplanar π(+)-π stacking alignment. Nevertheless, a comparative (13)C NMR spectroscopy study of the two homologous scaffolds revealed that the overall response of the indole unit arises predominantly from global attractions between the indole ring and the entire positively charged first coordination sphere. The study thus demonstrates the coordination-sphere/geometry dependence of the Trp/His side-chain interplay, and established that π interactions allow (13)C NMR spectroscopy to offer a new window for investigating Trp rotamer heterogeneity near metal-binding centers.
Collapse
Affiliation(s)
- Chi Ming Yang
- Neurochemistry & Biophysical Organic Chemistry, College of Chemistry, Nankai University, Tianjin 300071, China.
| | | |
Collapse
|
40
|
Sabirianov RF, Rubinstein A, Namavar F. Enhanced initial protein adsorption on engineered nanostructured cubic zirconia. Phys Chem Chem Phys 2011; 13:6597-609. [DOI: 10.1039/c0cp02389b] [Citation(s) in RCA: 30] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/02/2023]
|
41
|
Wallnoefer HG, Lingott T, Gutiérrez JM, Merfort I, Liedl KR. Backbone flexibility controls the activity and specificity of a protein-protein interface: specificity in snake venom metalloproteases. J Am Chem Soc 2010; 132:10330-7. [PMID: 20617834 DOI: 10.1021/ja909908y] [Citation(s) in RCA: 43] [Impact Index Per Article: 3.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022]
Abstract
Protein-protein interfaces have crucial functions in many biological processes. The large interaction areas of such interfaces show complex interaction motifs. Even more challenging is the understanding of (multi)specificity in protein-protein binding. Many proteins can bind several partners to mediate their function. A perfect paradigm to study such multispecific protein-protein interfaces are snake venom metalloproteases (SVMPs). Inherently, they bind to a variety of basement membrane proteins of capillaries, hydrolyze them, and induce profuse bleeding. However, despite having a high sequence homology, some SVMPs show a strong hemorrhagic activity, while others are (almost) inactive. We present computer simulations indicating that the activity to induce hemorrhage, and thus the capability to bind the potential reaction partners, is related to the backbone flexibility in a certain surface region. A subtle interplay between flexibility and rigidity of two loops seems to be the prerequisite for the proteins to carry out their damaging function. Presumably, a significant alteration in the backbone dynamics makes the difference between SVMPs that induce hemorrhage and the inactive ones.
Collapse
Affiliation(s)
- Hannes G Wallnoefer
- Institute of General, Inorganic and Theoretical Chemistry, Faculty of Chemistry and Pharmacy, University of Innsbruck, Innrain 52a, A-6020 Innsbruck, Austria
| | | | | | | | | |
Collapse
|
42
|
Dai L, Yang Y, Kim HR, Zhou Y. Improving computational protein design by using structure-derived sequence profile. Proteins 2010; 78:2338-48. [PMID: 20544969 DOI: 10.1002/prot.22746] [Citation(s) in RCA: 21] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Abstract
Designing a protein sequence that will fold into a predefined structure is of both practical and fundamental interest. Many successful, computational designs in the last decade resulted from improved understanding of hydrophobic and polar interactions between side chains of amino acid residues in stabilizing protein tertiary structures. However, the coupling between main-chain backbone structure and local sequence has yet to be fully addressed. Here, we attempt to account for such coupling by using a sequence profile derived from the sequences of five residue fragments in a fragment library that are structurally matched to the five-residue segments contained in a target structure. We further introduced a term to reduce low complexity regions of designed sequences. These two terms together with optimized reference states for amino-acid residues were implemented in the RosettaDesign program. The new method, called RosettaDesign-SR, makes a 12% increase (from 34 to 46%) in fraction of proteins whose designed sequences are more than 35% identical to wild-type sequences. Meanwhile, it reduces 8% (from 22% to 14%) to the number of designed sequences that are not homologous to any known protein sequences according to psi-blast. More importantly, the sequences designed by RosettaDesign-SR have 2-3% more polar residues at the surface and core regions of proteins and these surface and core polar residues have about 4% higher sequence identity to wild-type sequences than by RosettaDesign. Thus, the proteins designed by RosettaDesign-SR should be less likely to aggregate and more likely to have unique structures due to more specific polar interactions.
Collapse
Affiliation(s)
- Liang Dai
- School of Informatics, Indiana University Purdue University, Indianapolis, Indiana 46202, USA
| | | | | | | |
Collapse
|
43
|
Lassila JK. Conformational diversity and computational enzyme design. Curr Opin Chem Biol 2010; 14:676-82. [PMID: 20829099 PMCID: PMC2953567 DOI: 10.1016/j.cbpa.2010.08.010] [Citation(s) in RCA: 28] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/04/2010] [Revised: 08/06/2010] [Accepted: 08/06/2010] [Indexed: 11/22/2022]
Abstract
The application of computational protein design methods to the design of enzyme active sites offers potential routes to new catalysts and new reaction specificities. Computational design methods have typically treated the protein backbone as a rigid structure for the sake of computational tractability. However, this fixed-backbone approximation introduces its own special challenges for enzyme design and it contrasts with an emerging picture of natural enzymes as dynamic ensembles with multiple conformations and motions throughout a reaction cycle. This review considers the impact of conformational variation and dynamics on computational enzyme design and it highlights new approaches to addressing protein conformational diversity in enzyme design including recent advances in multi-state design, backbone flexibility, and computational library design.
Collapse
Affiliation(s)
- Jonathan K Lassila
- Department of Biochemistry, Stanford University, Stanford, CA 94305, USA.
| |
Collapse
|
44
|
Barakat NH, Barakat NH, Love JJ. Combined use of experimental and computational screens to characterize protein stability. Protein Eng Des Sel 2010; 23:799-807. [PMID: 20805093 DOI: 10.1093/protein/gzq052] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/12/2022] Open
Abstract
One of the primary goals of protein design is to engineer proteins with improved stability. Protein stability is a key issue for chemical, biotechnology and pharmaceutical industries. The development of robust proteins/enzymes with the ability to withstand the potentially harsh conditions of industrial operations is of high importance. A number of strategies are currently being employed to achieve this goal. Two particular approaches, (i) directed evolution and (ii) computational protein design, are quite powerful yet have only recently been combined or applied and analyzed in parallel. In directed evolution, libraries of variants are searched experimentally for clones possessing the desired properties. With computational methods, protein design algorithms are utilized to perform in silico screening for stable protein sequences. Here, we used gene libraries of an unstable variant of streptococcal protein G (Gbeta1) and an in vivo screening method to identify stabilized variants. Many variants with notably increased thermal stabilities were isolated and characterized. Concomitantly, computational techniques and protein design algorithms were used to perform in silico screening of the same destabilized variant of Gbeta1. The combined use, and critical analysis, of these methods promises to advance the field of protein design.
Collapse
Affiliation(s)
- Nora H Barakat
- Department of Chemistry and Biochemistry, San Diego State University, 5500 Campanile Dr, San Diego, CA 92182-1030, USA
| | | | | |
Collapse
|
45
|
Rubinstein A, Sabirianov RF, Mei WN, Namavar F, Khoynezhad A. Effect of the ordered interfacial water layer in protein complex formation: A nonlocal electrostatic approach. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2010; 82:021915. [PMID: 20866845 DOI: 10.1103/physreve.82.021915] [Citation(s) in RCA: 10] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 04/22/2010] [Indexed: 05/29/2023]
Abstract
Using a nonlocal electrostatic approach that incorporates the short-range structure of the contacting media, we evaluated the electrostatic contribution to the energy of the complex formation of two model proteins. In this study, we have demonstrated that the existence of an ordered interfacial water layer at the protein-solvent interface reduces the charging energy of the proteins in the aqueous solvent, and consequently increases the electrostatic contribution to the protein binding (change in free energy upon the complex formation of two proteins). This is in contrast with the finding of the continuum electrostatic model, which suggests that electrostatic interactions are not strong enough to compensate for the unfavorable desolvation effects.
Collapse
Affiliation(s)
- A Rubinstein
- Department of Biomedical Sciences and Surgery, Creighton University Medical Center, Omaha, Nebraska 68131, USA.
| | | | | | | | | |
Collapse
|
46
|
Aleksandrov A, Polydorides S, Archontis G, Simonson T. Predicting the Acid/Base Behavior of Proteins: A Constant-pH Monte Carlo Approach with Generalized Born Solvent. J Phys Chem B 2010; 114:10634-48. [DOI: 10.1021/jp104406x] [Citation(s) in RCA: 28] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Alexey Aleksandrov
- Laboratoire de Biochimie (CNRS UMR7654), Department of Biology, Ecole Polytechnique, 91128 Palaiseau, France, and Department of Physics, University of Cyprus, PO20537, CY1678, Nicosia, Cyprus
| | - Savvas Polydorides
- Laboratoire de Biochimie (CNRS UMR7654), Department of Biology, Ecole Polytechnique, 91128 Palaiseau, France, and Department of Physics, University of Cyprus, PO20537, CY1678, Nicosia, Cyprus
| | - Georgios Archontis
- Laboratoire de Biochimie (CNRS UMR7654), Department of Biology, Ecole Polytechnique, 91128 Palaiseau, France, and Department of Physics, University of Cyprus, PO20537, CY1678, Nicosia, Cyprus
| | - Thomas Simonson
- Laboratoire de Biochimie (CNRS UMR7654), Department of Biology, Ecole Polytechnique, 91128 Palaiseau, France, and Department of Physics, University of Cyprus, PO20537, CY1678, Nicosia, Cyprus
| |
Collapse
|
47
|
Han S, Yin S, Yi H, Mouhat S, Qiu S, Cao Z, Sabatier JM, Wu Y, Li W. Protein−Protein Recognition Control by Modulating Electrostatic Interactions. J Proteome Res 2010; 9:3118-25. [DOI: 10.1021/pr100027k] [Citation(s) in RCA: 31] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/17/2023]
Affiliation(s)
- Song Han
- State Key Laboratory of Virology, College of Life Sciences, Wuhan University, Wuhan, Hubei 430072, P. R. China, and ERT 62 ≪Ingénierie des peptides à visée thérapeutique≫, Université de la Méditerranée - Ambrilia Biopharma S.A., Faculté de Médecine Nord, Boulevard Pierre Dramard, 13916 Marseille Cédex 20, France
| | - Shijin Yin
- State Key Laboratory of Virology, College of Life Sciences, Wuhan University, Wuhan, Hubei 430072, P. R. China, and ERT 62 ≪Ingénierie des peptides à visée thérapeutique≫, Université de la Méditerranée - Ambrilia Biopharma S.A., Faculté de Médecine Nord, Boulevard Pierre Dramard, 13916 Marseille Cédex 20, France
| | - Hong Yi
- State Key Laboratory of Virology, College of Life Sciences, Wuhan University, Wuhan, Hubei 430072, P. R. China, and ERT 62 ≪Ingénierie des peptides à visée thérapeutique≫, Université de la Méditerranée - Ambrilia Biopharma S.A., Faculté de Médecine Nord, Boulevard Pierre Dramard, 13916 Marseille Cédex 20, France
| | - Stéphanie Mouhat
- State Key Laboratory of Virology, College of Life Sciences, Wuhan University, Wuhan, Hubei 430072, P. R. China, and ERT 62 ≪Ingénierie des peptides à visée thérapeutique≫, Université de la Méditerranée - Ambrilia Biopharma S.A., Faculté de Médecine Nord, Boulevard Pierre Dramard, 13916 Marseille Cédex 20, France
| | - Su Qiu
- State Key Laboratory of Virology, College of Life Sciences, Wuhan University, Wuhan, Hubei 430072, P. R. China, and ERT 62 ≪Ingénierie des peptides à visée thérapeutique≫, Université de la Méditerranée - Ambrilia Biopharma S.A., Faculté de Médecine Nord, Boulevard Pierre Dramard, 13916 Marseille Cédex 20, France
| | - Zhijian Cao
- State Key Laboratory of Virology, College of Life Sciences, Wuhan University, Wuhan, Hubei 430072, P. R. China, and ERT 62 ≪Ingénierie des peptides à visée thérapeutique≫, Université de la Méditerranée - Ambrilia Biopharma S.A., Faculté de Médecine Nord, Boulevard Pierre Dramard, 13916 Marseille Cédex 20, France
| | - Jean-Marc Sabatier
- State Key Laboratory of Virology, College of Life Sciences, Wuhan University, Wuhan, Hubei 430072, P. R. China, and ERT 62 ≪Ingénierie des peptides à visée thérapeutique≫, Université de la Méditerranée - Ambrilia Biopharma S.A., Faculté de Médecine Nord, Boulevard Pierre Dramard, 13916 Marseille Cédex 20, France
| | - Yingliang Wu
- State Key Laboratory of Virology, College of Life Sciences, Wuhan University, Wuhan, Hubei 430072, P. R. China, and ERT 62 ≪Ingénierie des peptides à visée thérapeutique≫, Université de la Méditerranée - Ambrilia Biopharma S.A., Faculté de Médecine Nord, Boulevard Pierre Dramard, 13916 Marseille Cédex 20, France
| | - Wenxin Li
- State Key Laboratory of Virology, College of Life Sciences, Wuhan University, Wuhan, Hubei 430072, P. R. China, and ERT 62 ≪Ingénierie des peptides à visée thérapeutique≫, Université de la Méditerranée - Ambrilia Biopharma S.A., Faculté de Médecine Nord, Boulevard Pierre Dramard, 13916 Marseille Cédex 20, France
| |
Collapse
|
48
|
Ishida T. Effects of Point Mutation on Enzymatic Activity: Correlation between Protein Electronic Structure and Motion in Chorismate Mutase Reaction. J Am Chem Soc 2010; 132:7104-18. [DOI: 10.1021/ja100744h] [Citation(s) in RCA: 27] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/28/2022]
Affiliation(s)
- Toyokazu Ishida
- Research Institute for Computational Sciences (RICS), National Institute of Advanced Industrial Science and Technology (AIST), Tsukuba Central 2, 1-1-1 Umezono, Tsukuba 305-8568, Japan
| |
Collapse
|
49
|
Abstract
Predictive methods for the computational design of proteins search for amino acid sequences adopting desired structures that perform specific functions. Typically, design of 'function' is formulated as engineering new and altered binding activities into proteins. Progress in the design of functional protein-protein interactions is directed toward engineering proteins to precisely control biological processes by specifically recognizing desired interaction partners while avoiding competitors. The field is aiming for strategies to harness recent advances in high-resolution computational modeling-particularly those exploiting protein conformational variability-to engineer new functions and incorporate many functional requirements simultaneously.
Collapse
Affiliation(s)
- Daniel J Mandell
- Graduate Program in Bioinformatics and Computational Biology, California Institute for Quantitative Biosciences, and Department of Bioengineering and Therapeutic Sciences, University of California San Francisco, USA
| | | |
Collapse
|
50
|
Tian P. Computational protein design, from single domain soluble proteins to membrane proteins. Chem Soc Rev 2010; 39:2071-82. [DOI: 10.1039/b810924a] [Citation(s) in RCA: 22] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
|